Restricted Gene Flow Research Articles

Genetic analysis reveals spatial structure in an expanding introduced rusa deer population

Context Rusa deer (Cervus timorensis), originally introduced in the 1860s, are still spreading in eastern Australia. The expanding peri-urban rusa deer population in the Illawarra region of New South Wales, Australia is having undesirable impacts on human and ecological communities, but the spatial structure of this population has not been investigated. Genetic information on invasive species is potentially useful in identifying management units to mitigate undesirable impacts. Aims The aim of this study was to investigate population structure, characterise dispersal, and determine if natural and human-made landscape features affected gene flow in rusa deer invading the Illawarra region of New South Wales. Methods We used reduced representation sequencing (DArT-Seq) to analyse single nucleotide polymorphisms distributed throughout the genomic DNA of rusa deer culled during a management program. We used admixture and Principal Component Analyses to investigate population structure with respect to natural and human-made landscape features, and we investigated whether our genetic data supported the presence of sex-biased dispersal. Key results Genetic diversity was highest in the north, near the original introduction site. A railway line demarcated restricted gene flow. Surprisingly, the Illawarra escarpment, a prominent landscape feature, did not restrict gene flow. There was no evidence of sex-biased dispersal and seven individuals were identified as genetic outliers. Conclusions The genetic structure of the Illawarra rusa deer population is consistent with individuals spreading south from their introduction site in Royal National Park. The population is not panmictic, and a landscape feature associated with urbanisation was associated with increased spatial genetic structure. Outliers could indicate hybridisation or secondary incursion events. Implications Rusa deer can be expected to continue invading southwards in the Illawarra region, but landscape features associated with urbanisation might reduce dispersal across the landscape. The genetic structuring of the population identified three potential management units on which to prioritise ground shooting operations.

Phylogeography and diversification of Oriental weaverbirds (Ploceus spp.): A gradual increase of eurytopy

Weaverbirds are a speciose group of colorful passerines inhabiting the Old World Tropics. Nevertheless, the Oriental weaverbirds (Ploceus spp.), widespread across southern Asia, are much less diverse and restricted to a few ecological niches compared to their African counterpart. To investigate their phylogeography, we retrieved 101 samples of Baya Weaver (P. philippinus), Streaked Weaver (P. manyar), Black-Throated Weaver (P. benghalensis) and Asian Golden Weaver (P. hypoxanthus) along with GenBank sequences of Finn's Weaver (P. megarhynchus). We reconstructed the first molecular phylogeny based on a dataset consisting of both mitochondrial and nuclear genes, dating the most recent common ancestor of Oriental Ploceus to ∼11 mya. Subsequent speciation appears to have been a combination of divergence within the Indian subcontinent and dispersal across a barrier situated between the Indian subcontinent and the Indochinese region, which provided habitats with a varying degree of isolations and ultimately promoted divergences in allopatry. Two descendants of the earliest nodes, P. megarhynchus and P. hypoxanthus, are both rare and local, often found near large river systems, which perhaps reflects niche conservatism and a lack of adaptive potential. The three smaller species are all widespread, common and less habitat specific. The most recent divergence, between western and eastern P. philippinus populations, is supported by both phylogenetic and morphological evidence, pointing toward limited gene flow between them. However, a zone of intergradation may exist in Myanmar and Brahmaputra flood plains, thus preventing a recommendation for species level recognition without further study.

Phylogeography of the Atlantic Blue Crab Callinectes sapidus (Brachyura: Portunidae) in the Americas versus the Mediterranean Sea: Determining Origins and Genetic Connectivity of a Large-Scale Invasion.

The American blue crab Callinectes sapidus is a particularly successful invader in estuarine ecosystems worldwide. Despite increasing awareness of its potential harm, the invasion history and underlying genetic diversity of this species within the Mediterranean Sea remain unknown. This study constitutes the first large-scale approach to study phylogeographic patterns of C. sapidus in Europe, facilitated by the first comparison of all currently available COI sequence data. For this investigation, 71 individuals of C. sapidus were newly analyzed and the entire COI gene was sequenced and used for a comparative phylogeographic analyses. For the first time, two separately used adjacent regions of this gene were combined in a single dataset. This allowed emphasizing the prevalence of three geographically defined lineages within the native range: (1) eastern North America, including the Gulf of Mexico, (2) the Caribbean, and (3) Brazil. New data from the Mediterranean reveal that non-native populations of C. sapidus are characterized by a conspicuously low genetic diversity (except for Turkey, where stocking took place), and that there is surprisingly low connectivity among established populations. The occurrence of strong genetic bottlenecks suggests few founder individuals. This confirms that, even under a scenario of restricted large-scale gene flow, a very limited number of invasive individuals is sufficient for a massive impact.

A pollen expressed PME gene at Tcb1 locus confers maize unilateral cross-incompatibility.

A pollen expressed PME gene at Tcb1 locus confers maize unilateral cross-incompatibility.

Dermacentor reticulatus – a tick on its way from glacial refugia to a panmictic Eurasian population

The ornate dog tick (Dermacentor reticulatus) shows a recently expanding geographic distribution. Knowledge on its intraspecific variability, population structure, rate of genetic diversity and divergence, including its evolution and geographic distribution, is crucial to understand its dispersal capacity. All such information would help to evaluate the potential risk of future spread of associated pathogens of medical and veterinary concern. A set of 865 D. reticulatus ticks was collected from 65 localities across 21 countries, from Portugal in the west to Kazakhstan and southern Russia in the east. Cluster analyses of 16 microsatellite loci were combined with nuclear (ITS2, 18S) and mitochondrial (12S, 16S, COI) sequence data to uncover the ticks’ population structures and geographical patterns. Approximate Bayesian computation was applied to model evolutionary relationships among the found clusters. Low variability and a weak phylogenetic signal showing an east–west cline were detected both for mitochondrial and nuclear sequence markers. Microsatellite analyses revealed three genetic clusters, where the eastern and western cluster gradient was supplemented by a third, northern cluster. Alternative scenarios could explain such a tripartite population structure by independent formation of clusters in separate refugia, limited gene flow connected with isolation by distance causing a “bipolar pattern”, and the northern cluster deriving from admixture between the eastern and western populations. The best supported demographic scenario of this tick species indicates that the northern cluster derived from admixture between the eastern and western populations 441 (median) to 224 (mode) generations ago, suggesting a possible link with the end of the Little Ice Age in Europe.

Population structure of mud flounder Paralichthys orbignyanus from the south-western Atlantic Ocean.

The mud flounder Paralichthys orbignyanus (Pleuronectiformes, Paralichthyidae) inhabits shallow waters of low salinities and mud bottoms in the temperate marine coastal regions of the Bonaerensean Ecoregion of the Argentinean Biogeographic Province in the south-western Atlantic Ocean. Specimens of P. orbignyanus were collected from Lagoa dos Patos (LDP) (southern Brazil), Mar Chiquita (MCH) and Marisol (MAR) both located in Buenos Aires (Argentina), and San Antonio Oeste (SAO) in the San Matías Gulf, Rio Negro (Argentina). A fragment of the mitochondrial DNA of the Control Region and seven microsatellite loci were characterized. In the Control Region, P. orbignyanus showed high variability, low nucleotide diversity, mild population expansion and a coalescence time of 35,000 years before the present. Flounders provided evidence of a genetic structure between the sampling sites LDP, MCH, MAR vs. SAO. On the other hand, P. orbignyanus displayed a lower to moderate contemporary genetic structure among all samples except between LDP and MCH. With no evidence of isolation by distance, this analysis supports a model of limited gene flow that is likely to be associated with a consistent larvae retention in all sampling sites. In addition, the present connectivity is ascribed to a lower migration process from SAO in the San Matías Gulf congruent with the prevailing littoral drift.

How Many Abalone Species Live in the Mediterranean Sea?

Morphological traits in Haliotidae may be highly variable and not consistently diagnostic for species identification, highlighting the need for an integrative approach to the taxonomy of the family, including genetic data. Four species of the genus Haliotis are currently reported for the Mediterranean Sea and the neighboring Atlantic Ocean: Haliotis tuberculata, the common European abalone with the widest Atlanto-Mediterranean range; Haliotis mykonosensis, from the Aegean, the Tyrrhenian, and the Adriatic; Haliotis stomatiaeformis, from Malta, Lampedusa, and southeastern Sicily; and the Lessepsian Haliotis pustulata, only known on the basis of few samples from the Levant. However, their taxonomic status still relies only on shell morphology. Here, sequences of two fragments of the mitochondrial molecular marker COI were obtained from 84 abalone specimens collected in the Mediterranean Sea and the neighboring Atlantic and analyzed in order to provide for the first time a genetic framework for species delimitation. This study’s results prove that H. mykonosensis is genetically identical to H. tuberculata, whereas H. stomatiaeformis is a distinct species, endemic to a restricted area of the southern Mediterranean Sea. Finally, Haliotis tuberculata coccinea from Macaronesia may deserve its status as a subspecies of H. tuberculata, with genetic signature of a limited gene flow found in specimens of the nominal subspecies (H. t. tuberculata) in both the Atlantic and the Mediterranean Sea.

Ice Age effects on genetic divergence of the American crocodile (Crocodylus acutus) in Panama: reconstructing limits of gene flow and environmental ranges: a reply to O'Dea et al.

O'Dea et al. (2022) (Pleistocene sea level changes and crocodile population histories on the isthmus of panama: a comment on Avila-Cervantes et al. (2020). Evolution, 76(11), 2778-2783. https://doi.org/10.1111/evo.14610) question our hypothesis that sea-level changes during the past glaciation played a role in restricting gene flow between Pacific and Caribbean Crocodylus acutus in Panama. They argue that an error in sea-level high-stand reconstruction during the last interglacial period (118-130 ka) does not support our hypothesis. Although they are correct in our high-stand reconstruction error, overlooked the point in that we were presenting a model of restricted gene flow across the Panamanian Isthmus during low sea levels. We review the assumptions of gene demographic methods, emphasizing that we were focusing on times of genetic divergence. We expand here why gene flow between these coastal populations could have been restricted during the last glacial maximum (19-26.5 ka) and the 50,000 years preceding it when sea levels were lower than today. O'Dea et al. suggest local climates may have played larger roles than sea levels. We demonstrate that paleoclimate estimates for the past 3.3 Ma in Panama are within the bounds of extant C. acutus. The importance of Ice Age Sea-level dynamics on Neotropical species was likely profound and should be incorporated into evolutionary studies of these taxa.

Flower morphology variation in five species of Penstemon (Plantaginaceae) displaying Hymenoptera pollination syndrome

Background: Geographic distance promotes phenotypic variation by facilitating environmental distance, limiting gene flow, and exposing plants to different pollen vectors. Therefore, understanding how plant morphology changes across a geographic range improves our understanding of the drivers of morphological diversification both on a macro- and micro-evolutionary scale.
 Questions: 1) How do geographic location and abiotic factors affect flower morphology between populations? 2) Is there a geographic pattern of flower morphology variation? and 3) How does yearly variation in temperature and precipitation affect flower morphology within populations?
 Studied species: Penstemon albidus, P. fruticosus, P. glandulosus, P. speciosus, and P. whippleanus
 Study site and dates: The continental USA, summers of 2017 and 2018
 Methods: Fifty-seven populations and 496 individuals were selected at random to measure ten floral traits. Bioclimatic variables were extracted from the WorldClim database and NOAA. Linear models, partial least squares regression, Mantel tests and canonical correlation analysis were used to analyze the data.
 Results: Geographic variables alone explained a significant portion of the variation in flower morphology in two species, while in others, flower morphology did not vary despite large geographic distances. Penstemon albidus and P. whippleanus flowers increase in size from south-north, while P. glandulosus and P. speciosus exhibited an east-west increasing trend. Additionally, mean annual precipitation was the most important variable influencing P. glandulosusflower morphology.
 Conclusions: Geographic distance facilitates isolation-by-distance and isolation-by-phenology as well as abiotic differences between populations; however, other factors such as pollinators might be keeping populations morphologically homogeneous despite large geographic distance.

As the Goose Flies: Migration Routes and Timing Influence Patterns of Genetic Diversity in a Circumpolar Migratory Herbivore

Migration schedules and the timing of other annual events (e.g., pair formation and molt) can affect the distribution of genetic diversity as much as where these events occur. The greater white-fronted goose (Anser albifrons) is a circumpolar goose species, exhibiting temporal and spatial variation of events among populations during the annual cycle. Previous range-wide genetic assessments of the nuclear genome based on eight microsatellite loci suggest a single, largely panmictic population despite up to five subspecies currently recognized based on phenotypic differences. We used double digest restriction-site associated DNA (ddRAD-seq) and mitochondrial DNA (mtDNA) sequence data to re-evaluate estimates of spatial genomic structure and to characterize how past and present processes have shaped the patterns of genetic diversity and connectivity across the Arctic and subarctic. We uncovered previously undetected inter-population differentiation with genetic clusters corresponding to sampling locales associated with current management groups. We further observed subtle genetic clustering within each management unit that can be at least partially explained by the timing and directionality of migration events along with other behaviors during the annual cycle. The Tule Goose (A. a. elgasi) and Greenland subspecies (A. a. flavirostris) showed the highest level of divergence among all sampling locales investigated. The recovery of previously undetected broad and fine-scale spatial structure suggests that the strong cultural transmission of migratory behavior restricts gene flow across portions of the species’ range. Our data further highlight the importance of re-evaluating previous assessments conducted based on a small number of highly variable genetic markers in phenotypically diverse species.

Fine-scale genetic diversity and genet dynamics of the fairy ring fungus Floccularia luteovirens on the Qinghai–Tibet plateau

The fungus Floccularia luteovirens is mainly distributed in the alpine meadows of the Qinghai–Tibet plateau. Its fruit bodies tend to form fairy rings with a visible stimulating zone. Our previous studies have investigated the large-scale genetic structure among wild populations of F. luteovirens, but the mechanisms underlying the current genotype distribution pattern remain unknown. The balance between sexual and asexual reproduction affects the establishment and structure of populations. Measuring genet size and density is an effective approach to investigating the reproduction strategies of this species. In the current study, 234 fruit bodies and 79 soil samples were collected from three sampling sites over 3y, revealing that F. luteovirens exhibits relatively large genets. Very few new genets were detected over the 3y, illustrating that this species relies more on vegetative growth and can persist for long periods underground as mycelia. Moreover, the underground genet data showed a close relationship with the above-ground genet data. Our study found limited fine-scale gene flow, contrary to our previous large-scale genetic study of F. luteovirens, the present study found limited fine-scale gene flow of this species. The commercial trade of F. luteovirens fruit bodies by humans probably enhances the gene flow over QTP.

Molecular phylogenetics and systematics of two enteric helminth parasites (Baylisascaris laevis and Diandrya vancouverensis) in the Vancouver Island marmot (Marmota vancouverensis).

Island biogeography can promote rapid diversification and speciation via geographic isolation and novel selection pressures. These same factors can threaten the persistence of island endemics by limiting gene flow and suitable habitat. Host-parasite interactions on islands introduce another dimension of complexity as both species must simultaneously adapt to exogenous and endogenous factors. One example of host-parasite island biogeography is the critically endangered Vancouver Island (VI) marmot (Marmota vancouverensis) which is endemic to VI, Canada, and hosts two enteric helminth parasites: Baylisascaris laevis, an ascarid nematode common in tribe Marmotini, and Diandrya vancouverensis, an anoplocephalid cestode endemic to the VI marmot. Here, we aligned novel sequences from B. laevis (six genes) and D. vancouverensis (two genes) with congeneric sequences from GenBank. Phylogenies reconstructed using Bayesian and maximum parsimony approaches consistently placed B. laevis in a morphoclade, and D. vancouverensis in a monophyletic clade sister to D. composita. Mean pairwise sequence divergence between D. vancouverensis and D. composita (9.06±1.94%) surpassed commonly accepted thresholds for species delimitation, whereas divergence between VI and mainland populations of B. laevis (1.12±0.78%) was comparable to (or sometimes greater than) pairwise divergence values between other Baylisascaris species. Disparity in the genetic divergence of each parasite may reflect differences in their life cycle, host specificity, virulence, and the chronological extent of their isolation. Detailed descriptions of the population genetic structure and effects of both parasites on their shared host are crucial next steps in understanding the history of B. laevis and D. vancouverensis on VI and informing conservation efforts for the VI marmot and its enteric helminth parasites.

Commentary: using prey naïveté to inform ecological management.

Commentary: using prey naïveté to inform ecological management.

Allopatric Lineage Divergence of the East Asian Endemic Herb Conandron ramondioides Inferred from Low-Copy Nuclear and Plastid Markers

The evolutionary histories of ornamental plants have been receiving only limited attention. We examined the origin and divergence processes of an East Asian endemic ornamental plant, Conandron ramondioides. C. ramondioides is an understory herb occurring in primary forests, which has been grouped into two varieties. We reconstructed the evolutionary and population demography history of C. ramondioides to infer its divergence process. Nuclear and chloroplast DNA sequences were obtained from 21 Conandron populations on both sides of the East China Sea (ECS) to explore its genetic diversity, structure, and population differentiation. Interestingly, the reconstructed phylogeny indicated that the populations should be classified into three clades corresponding to geographical regions: the Japan (Honshu+Shikoku) clade, the Taiwan-Iriomote clade, and the Southeast China clade. Lineage divergence between the Japan clade and the Taiwan-Iriomote and Southeast China clades occured 1.14 MYA (95% HPD: 0.82-3.86), followed by divergence between the Taiwan-Iriomote and Southeast China clades approximately 0.75 MYA (95% HPD: 0.45-1.3). Furthermore, corolla traits (floral lobe length to tube length ratios) correlated with geographical distributions. Moreover, restricted gene flow was detected among clades. Lastly, the lack of potential dispersal routes across an exposed ECS seafloor during the last glacial maximum suggests that migration among the Conandron clades was unlikely. In summary, the extant Conandron exhibits a disjunct distribution pattern as a result of vicariance rather than long-distance dispersal. We propose that allopatric divergence has occurred in C. ramondioides since the Pleistocene. Our findings highlight the critical influence of species' biological characteristics on shaping lineage diversification of East Asian relic herb species during climate oscillations since the Quaternary.

Interpopulation variation in inbreeding is primarily driven by tolerance of mating with relatives in a spermcasting invertebrate.

The degree to which individuals inbreed is a fundamental aspect of population biology shaped by both passive and active processes. Yet, the relative influences of random and non-random mating on the overall magnitude of inbreeding are not well characterized for many taxa. We quantified variation in inbreeding among qualitatively accessible and isolated populations of a sessile marine invertebrate (the colonial ascidian Lissoclinum verrilli) in which hermaphroditic colonies cast sperm into the water column for subsequent uptake and internal fertilization. We compared estimates of inbreeding to simulations predicting random mating within sites to evaluate if levels of inbreeding were (1) less than expected because of active attempts to limit inbreeding, (2) as predicted by genetic subdivision and passive inbreeding tolerance, or (3) greater than simulations due to active attempts to promote inbreeding via self-fertilization or a preference for related mates. We found evidence of restricted gene flow and significant differences in the genetic diversity of L. verrilli colonies among sites, indicating that on average colonies were weakly related in accessible locations, but their levels of relatedness matched that of first cousins or half-siblings on isolated substrates. Irrespective of population size, progeny arrays revealed variation in the magnitude of inbreeding across sites that tracked with the mean relatedness of conspecifics. Biparental reproduction was confirmed in most offspring (86%) and estimates of total inbreeding largely overlapped with simulations of random mating, suggesting that interpopulation variation in mother-offspring resemblance was primarily due to genetic subdivision and passive tolerance of related mates. Our results highlight the influence of demographic isolation on the genetic composition of populations, and support theory predicting that tolerance of biparental inbreeding, even when mates are closely related, may be favoured under a broad set of ecological and evolutionary conditions.

Limited genetic differentiation of Mycetomoellerius mikromelanos in Parque National Soberanía, Panama: Implications for queen dispersal

AbstractThe coevolutionary relationship between fungus‐growing ants (Formicidae: Attini: Attina) and their symbionts has been well studied in the Panamanian rain forests. To further understand the ecological context of these evolutionary relationships, we have examined the population‐genetic structure of the fungus‐growing ant species Mycetomoellerius mikromelanos Cardenas, Schultz, Adams 2021 in the Panama Canal Zone. We specifically investigated the presence of population structure, the significance of geographic features (i.e., creeks) limiting gene flow, and relatedness between ant colonies. To accomplish this, we genotyped 85 ant colonies from nine creeks across an approximately 30 km transect in Parque National Soberanía, Panama, using double‐digest restriction‐site‐associated DNA sequencing. We did not find distinct population structure using two genetic clustering methods; however, we did detect an effect of isolation by distance. Furthermore, related colonies were frequently detected on the same creek or neighboring creeks, and some at further geographic distances. Collectively, these findings demonstrate that new colonies tend to establish on natal creeks and occasionally on distant creeks following long‐distance dispersal events. We discuss how population‐genetic patterns reveal the natural history of M. mikromelanos in Parque National Soberanía and how these results fit into the context of fungus‐growing ant mutualisms.Abstract in Spanish is available with online material.

Fragmentation by major dams and implications for the future viability of platypus populations

The evolutionarily unique platypus (Ornithorhynchus anatinus) has experienced major declines and extinctions from a range of historical and recent interacting human-mediated threats. Although spending most of their time in the water, platypuses can move over land. Nevertheless, uncertainties remain whether dams are barriers to movement, thus limiting gene flow and dispersal, essential to evolution and ecology. Here we examined disruption of gene flow between platypus groups below and above five major dams, matched to four adjacent rivers without major dams. Genetic differentiation (FST) across dams was 4- to 20-fold higher than along similar stretches of adjacent undammed rivers; FST across dams was similar to differentiation between adjacent river systems. This indicates that major dams represent major barriers for platypus movements. Furthermore, FST between groups was correlated with the year in which the dam was built, increasing by 0.011 every generation, reflecting the effects of these barriers on platypus genetics. This study provides evidence of gene flow restriction, which jeopardises the long-term viability of platypus populations when groups are fragmented by major dams. Mitigation strategies, such as building of by-pass structures and translocation between upstream and downstream of the dam, should be considered in conservation and management planning.

Different Interspecies Demographic Histories within the Same Locality: A Case Study of Sea Cucumbers, Cuttlefish and Clams in Greek Waters

Coalescent methods in population genetics aim to detect biodiversity patterns, evolutionary mechanisms, and signatures of historical changes in effective population sizes with respect to the species fidelity. Restriction site-associated DNA sequencing (RADseq) was used to evaluate the population dynamics of invertebrate species within the same localities. New sequencing technologies, such as the ones employed by population genetics, could be used to improve the management and sustainability of marine and aquaculture resources. Sea cucumbers (Holothuria tubolosa) showed genetic differentiation patterns favoring limited gene flow between studied areas. Similar results for clams (Venus verrucosa) suggest local adaptation and low-dispersal abilities for sessile organisms. On the contrary, cuttlefish (Sepia officinalis) exhibited a panmictic pattern, resulting in a single genetic stock in the area. The larvae settlement duration may be responsible for such interspecies variations. Interspecies demographic modeling revealed different environmental pressures of historical events’ signatures with respect to the three invertebrates. Sea cucumbers favor a post-glacial bottleneck event followed by a more recent recovery, whereas cuttlefish favor an expansion before the late glacial maximum. Lastly, clams showed a constant effective population size in the area. The results of historical demographic changes in natural populations provide opportunities for critical evaluation and management in terms of the conservation of the species in the area.

Climate change threatens unique genetic diversity within the Balkan biodiversity hotspot – The case of the endangered stone crayfish

The stone crayfish Austropotamobius torrentium is globally endangered due to anthropogenic pressure on its habitats, climate change and the invasive crayfish species, particularly the signal crayfish Pacifastacus leniusculus . Aiming to guide A. torrentium conservation within the hotspot of its genetic diversity in the western Balkan Peninsula, we combined population genetics and species distribution models (SDMs) to reveal the impact of climate change and the invasive P. leniusculus on A. torrentium diversity. Population genetic analyses based on newly developed microsatellites revealed moderate within-population genetic diversity and high differentiation among populations, reflecting isolated populations with limited gene flow. Alongside strong genetic structuring, we discovered high level of inbreeding coefficient indicating homozygote excess within the majority of populations. The SDMs results predicted substantial reductions of suitable habitats for A. torrentium by 2070; 80 % of its currently suitable habitat is predicted to be lost under high‐warming climate change scenario. Obtained results indicated that 44 % of populations with high and/or unique genetic diversity, including three highly divergent and geographically restricted evolutionary lineages, are located in the areas predicted to become unsuitable in the future, highlighting their vulnerability to extinction. Further, SDMs revealed considerable decrease of future habitat suitability for P. leniusculus , suggesting that climate change represents greater threat to A. torrentium . Our study highlights the importance of conserving remnant A. torrentium populations, among other, through assisted migration and population mixing that could help populations overcome the risks of inbreeding and maladaptation, and thus enabling A. torrentium to withstand the ongoing climate change. • Native freshwater crayfish decline is a global conservation concern. • The first study of Austropotamobius torrentium genetic diversity in its hotspot (Balkan Peninsula). • Climate change poses a higher risk to A. torrentium than invasive crayfish species. • Transboundary conservation actions needed for species long-term preservation. • . Conservation strategies should include assisted migration and population reinforcement.

Comparative Phylogeography in the Taiwan–Luzon Volcanic Belt Indicates Fast Diversification History of Pachyrhynchus Weevils (Coleoptera: Curculionidae)

Abstract Geological isolation with limited gene flow can shape phenotypic diversification among different populations or species. From Luzon northward to Taiwan lies a series of isolated volcanic islands (the Taiwan–Luzon volcanic belt) where three Pachyrhynchus species [P. sarcitis (Behrens, 1887), P. nobilis (Heller, 1912) and P. semperi (Heller, 1912)] (Coleoptera: Curculionidae: Entiminae: Pachyrhynchini) are sympatrically distributed. With intraspecific color variation across different islands, these three weevil species provide a unique opportunity to study genetic divergence associated with geographic isolation and/or phenotypic differences. In this study, genome-wide single nucleotide polymorphism (SNP) data derived from ddRAD-seq was used to reconstruct the speciation history and to delimit species within each of the three species. Our results indicate all three species have ancient diversification histories in southern islands of the Taiwan–Luzon volcanic belt instead of more recent diversification histories in northern islands. Ancestral effective population size estimations also support a scenario of a series of founder colonization events from southern to northern islands in the volcanic belt. Most island populations are monophyletic, and early-divergent populations with deep genetic structures are supported statistically as distinct species. However, color variations could occur between populations with very recent diversification, indicating a fast-evolving rate in the change of coloration, possibly due to müllerian mimicry or founder effects across this region.