Search

ABSTRACT Giant clams (Tridacna and Hippopus) are large marine bivalves occupying tropical and subtropical reefs in the Indo‐Pacific. Giant clam populations have declined in many areas of the Indo‐Pacific and continue to be threatened by harvesting and environmental change. The small giant clam (Tridacna maxima) occurs throughout the Indo‐Pacific and has been subject to several phylogeographic studies across its range. However, given its broad range, there are several areas where the genetic diversity and connectivity of T. maxima populations has not been characterised. Here, we analyse the mitochondrial marker cytochrome oxidase 1 (CO1) to examine the genetic diversity and connectivity of T. maxima in two regions: Australia's Coral Sea Marine Park and the Cook Islands. Samples were collected from 13 reefs within the Coral Sea Marine Park and ten islands within the Cook Islands archipelago. Tridacna maxima across the sampled region of the Coral Sea did not display any population structure, whereas significant population structure was detected for T. maxima within the Cook Islands. For the Cook Islands, most pairwise comparisons involving an island in the northern group (Manihiki) were significant, as were comparisons for Palmerston (a more centrally located island) and the southern islands, Rarotonga and Mangaia. Both regions displayed high haplotype diversities (> 0.90), indicating that they are important repositories of genetic diversity. Additional CO1 data from throughout T. maxima's distribution showed that the Coral Sea clams belonged to the clade occurring in the South‐Western Pacific Ocean, whilst those from the Cook Islands belonged to a unique clade found in the Central Pacific Ocean. This clade extended from Fiji in the west to French Polynesia in the east and the atolls of Palmyra and Tarawa (Kiribati) in the north. Our assessment of genetic diversity and population structure in these regions will assist with management decisions for the species.

Abstract Genetic management is a critical component of threatened species conservation. Understanding spatial patterns of genetic diversity is essential for evaluating the resilience of fragmented populations to accelerating anthropogenic threats. Nowhere is this more relevant than on the Australian continent, which is experiencing an ongoing loss of biodiversity that exceeds any other developed nation. Using a proprietary genome complexity reduction‐based method (DArTSeq), we generated a data set of 3239 high quality Single Nucleotide Polymorphisms (SNPs) to investigate spatial patterns and indices of genetic diversity in the koala (Phascolarctos cinereus), a highly specialised folivorous marsupial that is experiencing rapid and widespread population declines across much of its former range. Our findings demonstrate that current management divisions across the state of New South Wales (NSW) do not fully represent the distribution of genetic diversity among extant koala populations, and that care must be taken to ensure that translocation paradigms based on these frameworks do not inadvertently restrict gene flow between populations and regions that were historically interconnected. We also recommend that koala populations should be prioritised for conservation action based on the scale and severity of the threatening processes that they are currently faced with, rather than placing too much emphasis on their perceived value (e.g., as reservoirs of potentially adaptive alleles), as our data indicate that existing genetic variation in koalas is primarily partitioned among individual animals. As such, the extirpation of koalas from any part of their range represents a potentially critical reduction of genetic diversity for this iconic Australian species.

Abstract Barley is one of the founder crops of the Neolithic transition in West Asia. While recent advances in genomics have provided a rather detailed picture of barley domestication, there are contradictory views on how the domestication process affected genetic diversity. We set out to revisit this question by integrating public DNA sequencing data from ancient barley and wide collections of extant wild and domesticated accessions. Using two previously overlooked approaches – analyses of chloroplast genomes and genome‐wide proportions of private variants – we found that the barley cultivated six millennia ago was genetically unique and more diverse when compared to extant landraces and cultivars. Moreover, the chloroplast genomes revealed a link between the ancient barley, an obscure wild genotype from north‐eastern Libya, and a distinct population of barley cultivated in Ethiopia/Eritrea. Based on these results, we hypothesize past existence of a wider North African population that included both wild and cultivated types and suffered from genetic erosion in the past six millennia, likely due to a rapid desertification that ended the Holocene African humid period. Besides providing clues about the origin of Ethiopian landraces, the hypothesis explains the post‐domestication loss of diversity observed in barley. Analyses of additional samples will be necessary to resolve the history of African barley and its contribution to the extant cultivated gene pool.

Abstract Global biodiversity loss continues unabated, and in Australia, the rate of recent mammal extinctions is among the worst in the world. Meanwhile, the diversity among and within many endemic mammal species remains undescribed. This information is crucial to delineate species boundaries and thus inform decision‐making for conservation. Sminthopsis virginiae (the red‐cheeked dunnart) is a small, dasyurid marsupial found in four disjunct populations around the northern coast of Australia and New Guinea. There are three currently recognized subspecies, each occupying a distinct geographic location. Sminthopsis v. virginiae occurs in Queensland, S. v. rufigenis is distributed across New Guinea and the Aru Islands, and S. v. nitela has populations in the Top End of the Northern Territory and the Kimberley region of Western Australia. Previous molecular work has suggested the current subspecies definitions are not aligned with DNA sequence data, though the sampling was limited. We undertook a comprehensive genetic and morphological review of S. virginiae to clarify relationships within the species. This included mitochondrial (CR, 12S, and cytb) and nuclear (omega‐globin, IRBP, and bfib7) loci, and morphometric analysis of skulls and whole wet‐preserved specimens held in museums. Maximum Likelihood and Bayesian phylogenetic analyses resolved samples into two distinct clades, demarcated by the Gulf of Carpentaria in Australia's north. Sminthopsis. v. nitela was consistently separated from S. v. virginiae and S. v. rufigenis, based on the overall body and skull size and craniodental features, while S. v. virginiae and S. v. rufigenis were more difficult to distinguish from each other. Thus, we redescribed S. virginiae, recognizing two species, S. nitela (raised from subspecies) and S. virginiae (now comprising the subspecies S. v. virginiae and S. v. rufigenis). This study highlights the importance of recognizing cryptic mammal fauna to help address the gap in our knowledge about diagnosing diversity during a time of conservation crisis.

Abstract Astydamia latifolia is the only species of the genus Astydamia, which forms an early‐diverging lineage of Apiaceae, subfamily Apioideae. This species is subendemic to the Canary Islands and one of the most representative species of the coastal environments of this archipelago. Astydamia displays diplochory, that is, diaspores with two long‐distance dispersal (LDD) syndromes. In particular, A. latifolia has both anemochorous and thalassochorous fruit traits (corky and winged mericarps). Although we expect this species to have a high dispersal capacity, there is no genetic study addressing it. The divergence time of this species from its sister taxon is also unknown. In this study, we aimed (i) to unveil the phylogenetic relationships and divergence times of A. latifolia; (ii) to reconstruct its phylogeographic structure across the Canary Islands; and (iii) to estimate the number of inter‐island colonization events. To these ends, we first sequenced the internal transcribed spacer (ITS) region for A. latifolia, reconstructed the phylogenetic relationships of Astydamia and Apiaceae relatives and estimated divergence times. Then, two plastid DNA regions (psaI‐aacD and psbK‐trnS) were sequenced for 158 individuals (from 36 Canarian population and one NW African population) to reconstruct a haplotype network. The recently developed method Phylogeographic Analysis of Island Colonization Events (PAICE) was used to estimate the number of inter‐island colonization events. Results show that A. latifolia is a phylogenetically isolated lineage that diverged from the most closely related genus (Molopospermum) in the Eocene–Miocene. It displays a low plastid DNA diversity (only four haplotypes detected), which is accompanied by a high degree of haplotype sharing between islands and highly linear rarefaction curves of colonization obtained in PAICE. These findings suggest a high colonization ability for this species, most likely related to the availability of two LDD syndromes.

Abstract The Yunnan‐Guizhou Plateau (YGP) is characterized by the distinctive isolated habitat of the limestone Karst Islands and features the Wumeng Mountains, which divide the YGP into the two Plateaus of Yunnan and Guizhou. This study aimed to assess the effects of geographic isolation and past climate fluctuation on the distribution of flora in the YGP. To achieve this, we carried out the phylogeographical pattern and genetic structure based on chloroplast and nuclear ribosomal DNA sequence in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modeling for Morella nana, an important wild plant resource and endemic to the YGP once considered a vulnerable species. The results suggest that the genetic and chlorotype network structures of M. nana are divided into at least two groups: cpDNA chlorotype H2 (or dominant nrDNA haplotypes h1 and h2), distributed primarily to the east of the Wumeng Mountains, and cpDNA chlorotypes H1 and H3–H10 (or dominant nrDNA haplotype h2 and h3), distributed to the west of the Wumeng Mountains. A deep genetic split was noted within the two groups to reach 25 steps, especially for the cpDNA fragment variation. This east–west divergence reveals the existence of a natural geographical isolation boundary in the form of the Wumeng Mountains, and supports the existence of at least two glacial refuges during the Quaternary glacial period, along with two genetic diversity center, and at least two large geographic protection units for the important species of M. nana. This study indicates that the phylogeographical pattern of M. nana can be attributed to geographic/environmental isolation caused by the Wumeng Mountains and climate fluctuation during the last glacial maximum, and proposes an effective strategy to protecting this important plant resource.

Abstract Ancient lakes are hotspots of species diversity, posing challenges and opportunities for exploration of the dynamics of endemic diversification. Lake Baikal in Siberia, the oldest lake in the world, hosts a particularly rich crustacean fauna, including the largest known species flock of harpacticoid copepods with some 70 species. Here, we focused on exploring the diversity and evolution within a single nominal species, Harpacticella inopinata Sars, 1908, using molecular markers (mitochondrial COI, nuclear ITS1 and 28S rRNA) and a set of qualitative and quantitative morphological traits. Five major mitochondrial lineages were recognized, with model‐corrected COI distances of 0.20–0.37. A concordant pattern was seen in the nuclear data set, and qualitative morphological traits also distinguish a part of the lineages. All this suggests the presence of several hitherto unrecognized cryptic taxa within the baikalian H. inopinata, with long independent histories. The abundances, distributions and inferred demographic histories were different among taxa. Two taxa, H. inopinata CE and H. inopinata CW, were widespread on the eastern and western coasts, respectively, and were largely allopatric. Patterns in mitochondrial variation, that is, shallow star‐like haplotype networks, suggest these taxa have spread through the lake relatively recently. Three other taxa, H. inopinata RE, RW and RW2, instead were rare and had more localized distributions on either coast, but showed deeper intraspecies genealogies, suggesting older regional presence. The rare taxa were often found in sympatry with the others and occasionally introgressed by mtDNA from the common ones. The mitochondrial divergence between and within the H. inopinata lineages is still unexpectedly deep, suggesting an unusually high molecular rate. The recognition of true systematic diversity in the evaluation and management of ecosystems is important in hotspots, as it is everywhere else, while the translation of the diversity into a formal taxonomy remains a challenge.

Abstract The climate fluctuations during the Late Pleistocene significantly influenced the phylogeographic structure and historical dynamics of marine fishes in the marginal seas of the western Pacific Ocean. The puffer fish, Lagocephalus spadiceus, holds substantial nutritional and economic value in the South China Sea. To investigate the demographic history and population structure of the L. spadiceus, the mitochondrial DNA COI and Cyt b gene datasets from 300 individuals across eight populations in the South China Sea were sequenced. Our findings revealed high haplotype diversity (0.874 ± 0.013) and low nucleotide diversity (0.00075 ± 0.00058). The phylogenetic tree and haplotype networks revealed no significant genetic differentiation along the northern coast of South China Sea. Neutrality tests, mismatch distribution analyses, and Bayesian skyline plots suggested that L. spadiceus underwent population expansion during the Late Pleistocene. Both ocean currents and climate change significantly influenced the geographical distribution and genetic population structure of L. spadiceus.

Abstract Erythroneurini is the largest tribe of the microleafhopper subfamily Typhlocybinae. Most prior research on this tribe has focused on traditional classification, phylogeny, and control of agricultural pests, and the phylogeography of the group remains poorly understood. In this study, the mitochondrial genomes of 10 erythroneurine species were sequenced, and sequences of four genes were obtained for 12 geographical populations of Seriana bacilla. The new sequence data were combined with previously available mitochondrial DNA sequence data and analyzed using Bayesian and Maximum‐Likelihood‐based phylogenetic methods to elucidate relationships among genera and species and estimate divergence times. Seriana was shown to be derived from within Empoascanara. Phylogeographic and population genetic analysis of the endemic Chinese species Seriana bacilla suggest that the species diverged about 54.85 Mya (95% HPD: 20.76–66.23 million years) in the Paleogene period and that population divergence occurred within the last 14 million years. Ancestral area reconstruction indicates that Seriana bacilla may have originated in the central region of Guizhou, and geographical barriers are the main factors affecting gene flow among populations. Ecological niche modeling using the MaxEnt model suggests that the distribution of the species was more restricted in the past but is likely to expand in the future years 2050 and 2070.

Abstract Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub‐patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high‐quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub‐populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility – crucial for the parasite's survival in diverse environments. Understanding these population‐level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Abstract Anthropogenic pressures on nature have been causing population declines for centuries. Intensified persecution of apex predators, like the golden eagle, resulted in population bottlenecks during the 19th and 20th centuries. To study population genetics and demographic history of the golden eagle throughout its distribution, we collected museum samples from previously underrepresented regions, such as Russia and Central Asia. We used 12 microsatellite loci and a fragment of the mitochondrial DNA control region to re‐evaluate phylogeography of Eurasian golden eagles and study the impacts of the population bottleneck. Our results revealed a north–south genetic gradient, expressed by the difference between Mediterranean and Holarctic lineages, as well as genetically distinct Northern Europe and Central Asia and Caucasus regions. Furthermore, Northern Europe exhibited the lowest, whereas Central Asia and Caucasus had the highest genetic diversity. Although golden eagles maintained relatively high genetic diversity, we detected genetic signatures of the recent bottleneck, including reduced genetic diversity and a decline in the effective female population size around the year 1975. Our study improves the knowledge of the genetic composition of Eurasian golden eagles and highlights the importance of understanding their historical population dynamics in the face of ongoing and future conservation efforts.

Abstract Camponotus japonicus (Hymenoptera: Formicidae) is an omnivorous social insect which builds sizable colonies in sparse woodlands or cropland and spreads across multiple climatic zones in the Chinese mainland. This study aims to reveal the role of climate changes and geological events in driving the genetic structure of social insect populations by investigating the phylogenetics and historical demography of C. japonicus in the Chinese mainland. Phylogenetic analyses were conducted based on the mitochondria DNA dataset using MrBayes and IQ‐TREE. We constructed a haplotype network, calculated analyses of molecular variance, estimated the divergence time, and reconstructed the maximum clade credibility tree. Mismatch distribution and Bayesian skyline plots were used to infer historical population fluctuations. Additionally, ecological niche modeling was employed to predict the potential distribution of the species during the present, mid‐holocene, and last glacial maximum periods in the Chinese mainland. The phylogenetic tree and median‐joining network analyses support the presence of four distinct lineages in C. japonicus. These lineages exhibit significant genetic differentiation and limited gene flow. The divergence among the four lineages began in the early Pleistocene, approximately 1.41 million years ago (Ma). Subsequently, the central lineage diverged from both the northern and southern lineages around 1.16 Ma, while the northern and southern lineages diverged from each other at approximately 1.07 Ma. Population expansion was observed in the southern, central, and northern lineages prior to the last glacial maximum, while the Yunnan‐Sichuan lineage experienced a slight increase in population size in more recent times. The predicted distribution of the species corresponds well with the actual distribution. Furthermore, the current suitable habitat areas in northern Xinjiang, southern Tibet, and the southeast coastal regions have significantly decreased compared to the last glacial maximum and the mid‐holocene periods. Our results suggest that climate oscillations and geological events play an important role in driving genetic patterns and differentiation of C. japonicus. Mountain barriers isolate populations from each other, hinder the flow of genes, and effectively prevent the spread of this species. But at the same time, it also formed refugia at low altitudes areas such as Qinling‐Bashan Mountains and Yanshan‐Taihang Mountains and provide suitable habitats during glaciation. This study provides a good model for understanding how complex climate changes and geological events affect population genetic differentiation of social insects in the Chinese mainland.

Abstract Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome‐wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid‐Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub‐central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post‐glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

Ecological restoration often requires translocating plant material from distant sites. Importing suitable plant material is important for successful establishment and persistence. Yet, published guidelines for seed transfer are available for very few species. Accurately predicting how transferred plants will perform requires multiyear and multi-environment field trials and comprehensive follow-up work, and is therefore infeasible given the number of species used in restoration programs. Alternative methods to predict the outcomes of seed transfer are valuable for species without published guidelines. In this study, we analyzed the genetic structure of an important shrub used in ecological restoration in the Southern Rocky Mountains called alder-leaf mountain mahogany (Cercocarpus montanus). We sequenced DNA from 1,440 plants in 48 populations across a broad geographic range. We found that genetic heterogeneity among populations reflected the complex climate and topography across which the species is distributed. We identified temperature and precipitation variables that were useful predictors of genetic differentiation and can be used to generate seed transfer recommendations. These results will be valuable for defining management and restoration practices for mountain mahogany.

Abstract Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi‐locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo‐observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non‐fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline—high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (FS and ZnS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Abstract The biogeographic history of the Chihuahuan Desert is complex, driven by numerous physiographic events and climatic changes. This dynamic history would have influenced the flora and fauna of the region including the desert pocket gopher, Geomys arenarius, a subterranean rodent endemic to the northern Chihuahuan Desert. G. arenarius is restricted to sandy soils and are considered to have a disjunct distribution. Two subspecies are recognized: G. a. arenarius and G. a. brevirostris. We used multilocus nuclear (amplified fragment length polymorphisms) and mitochondrial DNA (ND2) sequence data to uncover patterns of genetic diversity within and among populations of G. arenarius. We evaluated correspondence of genetic patterns to traditionally accepted subspecies boundaries, mapped the distribution of potentially suitable soils to identify barriers or corridors to dispersal and to guide future survey efforts, provided evidence that could be used to recognize distinct population segments, and quantified genetic diversity within populations. Both datasets were largely concordant and demonstrated hierarchical patterns of genetic divergence. The greatest divergence was consistent with the two recognized subspecies. Mapping of potentially habitable soils revealed likely barriers to dispersal contributing to the allopatric pattern of geographic distribution and areas, which may be occupied by G. arenarius but not yet documented. Because G. arenarius is restricted to soils with high sand content, and these habitable soils are disjunct within the region occupied by this species, historical factors that impacted soil deposition and deflation likely contributed to the observed patterns of genetic divergence. Genetic diversity was higher within populations of the southern subspecies (G. a. arenarius) compared to G. a. brevirostris. This may be due to a greater availability of continuous suitable soils within the range of G. a. arenarius or higher density due to greater food availability (currently or historically)—both of which could allow for a higher effective population size.

Abstract Studies on host–parasite systems that have experienced distributional shifts, range fragmentation, and population declines in the past can provide information regarding how parasite community richness and genetic diversity will change as a result of anthropogenic environmental changes in the future. Here, we studied how sequential postglacial colonization, shifts in habitat, and reduced host population sizes have influenced species richness and genetic diversity of Corynosoma (Acanthocephala: Polymorphidae) parasites in northern European marine, brackish, and freshwater seal populations. We collected Corynosoma population samples from Arctic, Baltic, Ladoga, and Saimaa ringed seal subspecies and Baltic gray seals, and then applied COI barcoding and triple‐enzyme restriction‐site associated DNA (3RAD) sequencing to delimit species, clarify their distributions and community structures, and elucidate patterns of intraspecific gene flow and genetic diversity. Our results showed that Corynosoma species diversity reflected host colonization histories and population sizes, with four species being present in the Arctic, three in the Baltic Sea, two in Lake Ladoga, and only one in Lake Saimaa. We found statistically significant population‐genetic differentiation within all three Corynosoma species that occur in more than one seal (sub)species. Genetic diversity tended to be high in Corynosoma populations originating from Arctic ringed seals and low in the landlocked populations. Our results indicate that acanthocephalan communities in landlocked seal populations are impoverished with respect to both species and intraspecific genetic diversity. Interestingly, the loss of genetic diversity within Corynosoma species seems to have been less drastic than in their seal hosts, possibly due to their large local effective population sizes resulting from high infection intensities and effective intra‐host population mixing. Our study highlights the utility of genomic methods in investigations of community composition and genetic diversity of understudied parasites.

Abstract Species with wide‐range distributions usually display high genetic variation. This variation can be partly explained by historical lineages that were temporally isolated from each other and are back into secondary reproductive contact, and partly by local adaptations. The smooth newt (Lissotriton vulgaris) is one of the most widely distributed amphibians species across Eurasia and forms a species complex with a partially overlapping distribution and morphology. In the present study, we explored the population genomic structure of smooth newt lineages in the Carpathian Basin (CB) relying on single‐nucleotide polymorphisms. Our dataset included new and previously published data to study the secondary contact zone between lineages in the CB and also tested for the barrier effect of rivers to gene flow between these lineages. We confirmed the presence of the South L. v. vulgaris Lineage distributed in Transdanubia and we provided new distribution records of L. v. ampelensis inhabiting the eastern territories of the CB. High genetic diversity of smooth newts was observed, especially in the North Hungarian Mountains and at the interfluves of the main rivers in the South with four distinct lineages of L. v. vulgaris and one lineage of L. v. ampelensis showing a low level of admixture with the spatially closest L. v. vulgaris lineage. Moreover, admixture detected at the interfluve of the main rivers (i.e. Danube and Tisza) suggested a secondary contact zone in the area. Finally, we found that the river Danube has a very weak effect on population divergence, while the river Tisza is a geographical barrier limiting gene flow between smooth newt lineages. As the range boundaries of L. v. vulgaris and L. v. ampelensis in the CB coincide with the river Tisza, our study underpins the influence of rivers in lineage diversification.

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three-spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self-recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.

Abstract The Great American Biotic Interchange (GABI) was a key biogeographic event in the history of the Americas. The rising of the Panamanian land bridge ended the isolation of South America and ushered in a period of dispersal, mass extinction, and new community assemblages, which sparked competition, adaptation, and speciation. Diversification across many bird groups, and the elevational zonation of others, ties back to events triggered by the GABI. But the exact timing of these events is still being revealed, with recent studies suggesting a much earlier time window for faunal exchange, perhaps as early as 20 million years ago (Mya). Using a time‐calibrated phylogenetic tree, we show that the jay genus Cyanolyca is emblematic of bird dispersal trends, with an early, pre‐land bridge dispersal from Mesoamerica to South America 6.3–7.3 Mya, followed by a back‐colonization of C. cucullata to Mesoamerica 2.3–4.8 Mya, likely after the land bridge was complete. As Cyanolyca species came into contact in Mesoamerica, they avoided competition due to a prior shift to lower elevation in the ancestor of C. cucullata. This shift allowed C. cucullata to integrate itself into the Mesoamerican highland avifauna, which our time‐calibrated phylogeny suggests was already populated by higher‐elevation, congeneric dwarf‐jays (C. argentigula, C. pumilo, C. mirabilis, and C. nanus). The outcome of these events and fortuitous elevational zonation was that C. cucullata could continue colonizing new highland areas farther north during the Pleistocene. Resultingly, four C. cucullata lineages became isolated in allopatric, highland regions from Panama to Mexico, diverging in genetics, morphology, plumage, and vocalizations. At least two of these lineages are best described as species (C. mitrata and C. cucullata). Continued study will further document the influence of the GABI and help clarify how dispersal and vicariance shaped modern‐day species assemblages in the Americas.

Abstract The Sunda pangolin (Manis javanica) is the most widely distributed Asian pangolin species, occurring across much of Southeast Asia and in southern China. It is classified as Critically Endangered and is one of the most trafficked mammals in the world, which not only negatively impacts wild Sunda pangolin populations but also poses a potential disease risk to other species, including humans and livestock. Here, we aimed to investigate the species' phylogeography across its distribution to improve our understanding of the species' evolutionary history, elucidate any taxonomic uncertainties and enhance the species' conservation genetic management and potential wildlife forensics applications. We sequenced mtDNA genomes from 23 wild Sunda pangolins of known provenance originating from Malaysia to fill sampling gaps in previous studies, particularly in Borneo. To conduct phylogenetic and population genetic analyses of Sunda pangolins across their range, we integrated these newly generated mitochondrial genomes with previously generated mtDNA and nuclear DNA data sets (RAD‐seq SNP data). We identified an evolutionarily distinct mtDNA lineage in north Borneo, estimated to be ~1.6 million years divergent from lineages in west/south Borneo and the mainland, comparable to the divergence time from the Palawan pangolin. There appeared to be mitonuclear discordance, with no apparent genetic structure across Borneo based on analysis of nuclear SNPs. These findings are consistent with the ‘out of Borneo hypothesis’, whereby Sunda pangolins diversified in Borneo before subsequently migrating throughout Sundaland, and/or a secondary contact scenario between mainland and Borneo. We have elucidated possible taxonomic issues in the Sunda/Palawan pangolin complex and highlight the critical need for additional georeferenced samples to accurately apportion its range‐wide genetic variation into appropriate taxonomic and conservation units. Additionally, these data have improved forensic identification testing involving these species and permit the implementation of geographic provenance testing in some scenarios.

Abstract Pairs of species that exhibit broadly overlapping distributions, and multiple geographically isolated contact zones, provide opportunities to investigate the mechanisms of reproductive isolation. Such naturally replicated systems have demonstrated that hybridization rates can vary substantially among populations, raising important questions about the genetic basis of reproductive isolation. The topminnows, Fundulus notatus and F. olivaceus, are reciprocally monophyletic, and co‐occur in drainages throughout much of the central and southern United States. Hybridization rates vary substantially among populations in isolated drainage systems. We employed genome‐wide sampling to investigate geographic variation in hybridization, and to assess the possible importance of chromosome fusions to reproductive isolation among nine separate contact zones. The species differ by chromosomal rearrangements resulting from Robertsonian (Rb) fusions, so we hypothesized that Rb fusion chromosomes would serve as reproductive barriers, exhibiting steeper genomic clines than the rest of the genome. We observed variation in hybridization dynamics among drainages that ranged from nearly random mating to complete absence of hybridization. Contrary to predictions, our use of genomic cline analyses on mapped species‐diagnostic SNP markers did not indicate consistent patterns of variable introgression across linkage groups, or an association between Rb fusions and genomic clines that would be indicative of reproductive isolation. We did observe a relationship between hybridization rates and population phylogeography, with the lowest rates of hybridization tending to be found in populations inferred to have had the longest histories of drainage sympatry. Our results, combined with previous studies of contact zones between the species, support population history as an important factor in explaining variation in hybridization rates.

Abstract Madagascar exhibits exceptionally high levels of biodiversity and endemism. Models to explain the diversification and distribution of species in Madagascar stress the importance of historical variability in climate conditions which may have led to the formation of geographic barriers by changing water and habitat availability. The relative importance of these models for the diversification of the various forest‐adapted taxa of Madagascar has yet to be understood. Here, we reconstructed the phylogeographic history of Gerp's mouse lemur (Microcebus gerpi) to identify relevant mechanisms and drivers of diversification in Madagascar's humid rainforests. We used restriction site associated DNA (RAD) markers and applied population genomic and coalescent‐based techniques to estimate genetic diversity, population structure, gene flow and divergence times among M. gerpi populations and its two sister species M. jollyae and M. marohita. Genomic results were complemented with ecological niche models to better understand the relative barrier function of rivers and altitude. We show that M. gerpi diversified during the late Pleistocene. The inferred ecological niche, patterns of gene flow and genetic differentiation in M. gerpi suggest that the potential for rivers to act as biogeographic barriers depended on both size and elevation of headwaters. Populations on opposite sides of the largest river in the area with headwaters that extend far into the highlands show particularly high genetic differentiation, whereas rivers with lower elevation headwaters have weaker barrier functions, indicated by higher migration rates and admixture. We conclude that M. gerpi likely diversified through repeated cycles of dispersal punctuated by isolation to refugia as a result of paleoclimatic fluctuations during the Pleistocene. We argue that this diversification scenario serves as a model of diversification for other rainforest taxa that are similarly limited by geographic factors. In addition, we highlight conservation implications for this critically endangered species, which faces extreme habitat loss and fragmentation.

Abstract The genus Symplocarpus in basal Araceae includes both thermogenic and non/slightly thermogenic species that prefer cold environments. If floral thermogenesis of Symplocarpus contributes to cold adaptation, it would be expected that thermogenic species have a larger habitat than non/slightly thermogenic species during an ice age, leading to increased genetic diversity in the current population. To address this question, potential distribution in past environment predicted by ecological niche modeling (ENM), genetic diversity, and population structure of chloroplast and genome‐wide single nucleotide polymorphisms were compared between thermogenic Symplocarpus renifolius and non/slightly thermogenic Symplocarpus nipponicus. ENM revealed that the distribution of S. nipponicus decreased, whereas that of S. renifolius expanded in the Last Glacial Maximum. Phylogeographic analyses have shown that the population structures of the two species were genetically segmented and that the genetic diversity of S. renifolius was higher than that of S. nipponicus. The phylogenetic relationship between chloroplast and nuclear DNA is topologically different in the two species, which may be due to the asymmetric gene flow ubiquitously observed in plants. The results of this study imply that floral thermogenesis of Symplocarpus contributes to expanding the distribution during an ice age, resulting in increased genetic diversity due to cold adaptation.

Abstract For many arctic species, the spatial (re‐)colonization patterns after the last Pleistocene glaciation have been described. However, the temporal aspects of their colonization are largely missing. Did one route prevail early, while another was more important later? The high Arctic archipelago Svalbard represents a good model system to address timeframe of postglacial plant colonization. Svalbard was almost fully glaciated during last glacial maximum and (re‐)colonization of vascular plants began in early Holocene. Early Holocene climatic optimum (HCO) supported an expanded establishment of a partly thermophilic vegetation. Today, we find remnants of this vegetation in sheltered regions referred to as “Arctic biodiversity hotspots”. The oldest record of postglacial plant colonization to Svalbard is found in Ringhorndalen‐Flatøyrdalen. Even though thermophilic species could establish also later in Holocene, only HCO was favorable for vast colonization, and only hotspots offered stable conditions for thermophilic populations throughout Holocene. Thus, these relic populations may reflect colonization patterns of HCO. We investigate whether the colonization direction of thermophilic plants (Arnica angustifolia, Campanula uniflora, Pinguicula alpina, Tofieldia pusilla, and Vaccinium uliginosum ssp. microphyllum) in Ringhorndalen‐Flatøyrdalen was uniform and different from later colonization events in other localities and non‐thermophilic plants (Arenaria humifusa, Bistorta vivipara, Juncus biglumis, Oxyria digyna, and Silene acaulis). We analyzed plastid haplotypes of the 10 taxa from Ringhorndalen‐Flatøyrdalen, from later‐colonized localities in Svalbard, and from putative source regions outside Svalbard. Only rare and thermophilic taxa Campanula uniflora and Vaccinium uliginosum ssp. microphyllum provided results suggesting at least two colonization events from different source regions. Tofieldia pusilla and all the non‐thermophilic plants showed no clear phylogeographically differentiation within Svalbard. Two of the thermophilic species showed no sequence variation. Based on the results, a uniform colonization direction to Svalbard in early Holocene is not probable; several source areas and dispersal directions were contemporarily involved.

The genus Stasimopus is endemic to South Africa but has never undergone a phylogeographic review. This study aims to unravel the phylogeographic patterns and history of the many Stasimopus species which occur in the greater Karoo region. A fossil-calibrated phylogeny was produced based on three gene regions (CO1, 16S and EF-1ɣ) for Stasimopus (Cor-k-lid trapdoor spiders) specimens collected in the Karoo region, to infer dates of origin and diversification. Demographic analyses were performed on species with sufficient sample sizes (> 4). Haplotype networks were constructed for each gene region and plotted on a map to infer phylogeographic patterns. Lastly, Mantel tests were performed to test for isolation by distance. It was found that 15 species occur in the Karoo and that the genus radiation in the area is in the early Palaeocene. Most diversification occurred between the late Eocene and the Miocene, coinciding with significant changes in climate. Several species show signals of demographic expansions. Isolation by distance was detected, but only with a slight correlation. It is apparent that aridification has played a vital role in the diversification of the genus in the Karoo region. This is a shared biogeographic influence between the mygalomorph fauna of the Karoo and arid region of western Australia. Stasimopus has radiated from the late Eocene and through the Miocene resulting in 15 extant species in the region. The Tankwa Karoo has been identified as a possible Pleistocene glacial cycle refugia for the species S. leipoldti . Many of the species in the Karoo are short-range endemics, making them of high conservation concern. This study provided vital information as the Karoo is undergoing further desertification due to factors such as climate change, which may affect the future of short-range endemic spiders., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Sky islands provide insights on how glacial-interglacial cycles have shaped species distribution and help for predicting species' responses to climate warming. The alpine subnival belt of southwest China, especially in the Hengduan Mountains and adjacent areas, is sky island-like. Among them, the Yunnan-Kweichow Plateau harbors several isolated mountains with well-developed alpine subnival vegetation, sharing a similar species composition with the Hengduan Mountains. However, the relationship between the sky islands of the Hengduan Mountains and the Yunnan-Kweichow Plateau remains insufficiently explored. Pleurospermum foetens (Apiaceae) is a species endemic to the alpine screes of the Yunnan-Kweichow Plateau and the Hengduan Mountains. We used DNA sequence data from 59 individuals across 9 populations, combined with ecological niche modeling, to investigate the evolution history and future distribution of P. foetens within this sky island region. The results indicate the following: (1) P. foetens exhibits a significant phylogeographic structure and can be classified into three nrDNA clades and two cpDNA clades, respectively, (2) a nuclear-plastid discordance observed in P. foetens and its relatives based on phylogenetic analysis. P. foetens is monophyletic in the nrDNA phylogeny, while two major clades (HDM and YGP) are present in the cpDNA phylogeny, each forming a clade with other congeneric species. (3) Ecological niche modeling of P. foetens indicated that the species had the most extensive suitable habitat during the last glacial maximum (LGM). However, anticipated climate warming in the coming decades is expected to reduce the suitable range of P. foetens , posing a significant threat to isolated marginal populations (e.g., Shizi Mountain) with restricted alpine scree habitats. In conclusion, our study highlights the substantial effect of sky island and glacial-interglacial cycles on the population divergence of P. foetens . Conservation efforts for marginal populations of alpine plants in the Yunnan-Kweichow Plateau require increased attention and prioritization., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Giant clams ( Tridacna and Hippopus ) are large marine bivalves occupying tropical and subtropical reefs in the Indo-Pacific. Giant clam populations have declined in many areas of the Indo-Pacific and continue to be threatened by harvesting and environmental change. The small giant clam ( Tridacna maxima ) occurs throughout the Indo-Pacific and has been subject to several phylogeographic studies across its range. However, given its broad range, there are several areas where the genetic diversity and connectivity of T. maxima populations has not been characterised. Here, we analyse the mitochondrial marker cytochrome oxidase 1 (CO1) to examine the genetic diversity and connectivity of T. maxima in two regions: Australia's Coral Sea Marine Park and the Cook Islands. Samples were collected from 13 reefs within the Coral Sea Marine Park and ten islands within the Cook Islands archipelago. Tridacna maxima across the sampled region of the Coral Sea did not display any population structure, whereas significant population structure was detected for T. maxima within the Cook Islands. For the Cook Islands, most pairwise comparisons involving an island in the northern group (Manihiki) were significant, as were comparisons for Palmerston (a more centrally located island) and the southern islands, Rarotonga and Mangaia. Both regions displayed high haplotype diversities (> 0.90), indicating that they are important repositories of genetic diversity. Additional CO1 data from throughout T. maxima 's distribution showed that the Coral Sea clams belonged to the clade occurring in the South-Western Pacific Ocean, whilst those from the Cook Islands belonged to a unique clade found in the Central Pacific Ocean. This clade extended from Fiji in the west to French Polynesia in the east and the atolls of Palmyra and Tarawa (Kiribati) in the north. Our assessment of genetic diversity and population structure in these regions will assist with management decisions for the species., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Abstract Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life‐history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping‐stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Abstract The high mountain ranges of Western Europe had a profound effect on the biotic recolonization of Europe from glacial refugia. The Alps present a particularly interesting case because they form an absolute barrier to dispersal for most taxa, obstructing recolonization from multiple refugia in northern Italy. Here, we investigate the effect of the European Alps on the phylogeographic history of the European common frog Rana temporaria. Based on partial cytochrome b and COXI sequences from Switzerland, we find two mitochondrial lineages roughly north and south of the Alpine ridge, with contact zones between them in eastern and western Switzerland. The northern haplogroup falls within the previously identified Western European haplogroup, while the southern haplogroup is unique to Switzerland. We find that the lineages diverged ~110 kya, at approximately the onset of the last glacial glaciation; this indicates that they are from different glacial refugia. Phylogenetic analyses suggest that the northern and southern haplogroups colonized Switzerland via trans‐ and circum‐Alpine routes from at least two separate refugia in northern Italy. Our results illustrate how a complex recolonization history of the central European Alps can arise from the semi‐permeable barrier created by high mountains.

Abstract We assessed population structure and the spatio‐temporal pattern of diversification in the Glossy Antshrike Sakesphorus luctuosus (Aves, Thamnophilidae) to understand the processes shaping the evolutionary history of Amazonian floodplains and address unresolved taxonomic controversies surrounding its species limits. By targeting ultraconserved elements (UCEs) from 32 specimens of S. luctuosus, we identified independent lineages and estimated their differentiation, divergence times, and migration rates. We also estimated current and past demographic histories for each recovered lineage. We found evidence confirming that S. luctuosus consists of a single species, comprising at least four populations, with some highly admixed individuals and overall similar levels of migration between populations. We confirmed the differentiation of the Araguaia River basin population (S. l. araguayae) and gathered circumstantial evidence indicating that the taxon S. hagmanni may represent a highly introgressed population between three distinct phylogroups of S. luctuosus. Divergences between populations occurred during the last 1.2 mya. Signs of population expansions were detected for populations attributed to subspecies S. l. luctuosus, but not for the S. l. araguayae population. Our results support that S. luctuosus has had a complex population history, resulting from a high dependence on southeastern “clear water” seasonally flooded habitats and their availability through time. Spatial and demographic expansions toward the western “white water” flooded forests might be related to recent changes in connectivity and availability of these habitats. Our study reinforces the view that isolation due to absence of suitable habitat has been an important driver of population differentiation within Amazonian flooded forests, but also that differences between várzeas (“white water” floodplains, mostly in southwestern Amazonia) and igapós (“clear water” floodplains, especially located in the east) should be further explored as drivers of micro‐evolution for terrestrial species.

Abstract Although least chipmunks (Neotamias minimus) are a widely distributed North American species of least concern, the southernmost population, N. m. atristriatus (Peñasco least chipmunk), is imperiled and a candidate for federal listing as a subspecies. We conducted a phylogeographic analysis across the N. minimus range to assess genomic differentiation and distinctiveness of the N. m. atristriatus population. Additionally, we leveraged the historical component of sampling to conduct a temporal analysis of N. minimus genetic diversity and also considered climate change effects on range persistence probability by projecting a species distribution model into the IPCC5 RCP 2.6 and 8.5 scenarios. We identified three geographically structured groups (West, North, and South) that were supported by both mitochondrial and nuclear data. N. m. atristriatus grouped within a unique South subclade but were not reciprocally monophyletic from N. m. operarius, and nuclear genome analyses did not separate N. m. atristriatus, N. m. caryi, and N. m. operarius. Thus, while least chipmunks in the Southwest represent an evolutionary significant unit, subspecies distinctions were not supported and listing of the Peñasco population as a Distinct Population Segment of N. m. operarius may be warranted. Our results also support consideration of populations with North and West mitogenomes as two additional evolutionary significant units. We found that N. minimus genetic diversity declined by ~87% over the last century, and our models predicted substantial future habitat contraction, including the loss of the full contemporary ranges of N. m. atristriatus, N. m. arizonensis, and N. m. chuskaensis.

Abstract Geological events, landscape features, and climate fluctuations have shaped the distribution of genetic diversity and evolutionary history in freshwater fish, but little attention has been paid to that around the Gulf of Tonkin; therefore, we investigated the phylogeographic structure of the dwarf snakehead (Channa gachua) on Hainan Island and mainland China, as well as two populations in Vietnam. We attempted to elucidate the origins of freshwater fish in South Hainan by incorporating genetic data from DNA markers on both the mitochondrial cytochrome b gene (cyt b) and the nuclear recombination‐activating gene 1 (RAG‐1). Mitochondrial phylogenetic analysis identified two major lineages (lineages A and B), which may represent separate species. Divergence data suggested that C. gachua populations diverged between 0.516 and 2.376 myr. The divergence of the two cryptic species is congruent with sea‐level rise, which subsequently isolated Hainan from the mainland. During the Pleistocene glaciations, the entire region of the Gulf of Tonkin and the Qiongzhou Strait became part of the coastal plain of the Asian continent, which might have resulted in the current distribution patterns and dispersal routes of C. gachua populations. The formation of three sublineages in lineage A indicated that the Gulf of Tonkin was a geographical barrier between Hainan Island and mainland China but not between Vietnam and Hainan Island. The results of this study may help to elucidate the origins of freshwater fish in South Hainan and the phylogeographic structure of C. gachua.

Abstract The Northwest Pacific marginal waters comprising the South China Sea, East China Sea, Yellow Sea, and the Sea of Japan have unique geomorphic features. The Japanese flounder Paralichthys olivaceus, which is endemic to the Northwest Pacific, has high nutritional, economic, and ecological value. To allow the examination of the demographic history and population structure of the most common P. olivaceus species range over the five marginal seas (East China Sea, Yellow Sea, Bohai Sea, Northwest Pacific Ocean, and the Sea of Japan), the mitochondrial DNA control region of 91 individuals from six populations in China was sequenced. These sequences were combined with 233 sequences from four populations distributed in the Northwestern Pacific Ocean for analysis. Higher levels of nucleotide diversity (0.032 ± 0.016) and haplotype diversity (0.996 ± 0.001) were observed. The peripheral Fuqing population in the East China Sea had the relatively lowest genetic diversity and highest differentiation. Furthermore, when the results of the isolation by distance test, spatial analysis of molecular variation and geographic barrier analysis are also considered, there is a clear need to prioritize resource conservation and enhancement measures in this area. The phylogenetic trees, structure assignment test, and haplotypes network revealed no significant differences in the genealogical structure among ten populations. Mismatch distribution analysis, Bayesian skyline plots, and neutrality tests suggested that P. olivaceus experienced population expansion during the Pleistocene. Ocean currents and climate change play important roles in shaping the geographical distribution and genetic population structure of P. olivaceus.

Abstract Both intrinsic and extrinsic forces work together to shape connectivity and genetic variation in populations across the landscape. Here we explored how geography, breeding system traits, and environmental factors influence the population genetic patterns of Triodanis perfoliata, a widespread mix‐mating annual plant in the contiguous US. By integrating population genomic data with spatial analyses and modeling the relationship between a breeding system and genetic diversity, we illustrate the complex ways in which these forces shape genetic variation. Specifically, we used 4705 single nucleotide polymorphisms to assess genetic diversity, structure, and evolutionary history among 18 populations. Populations with more obligately selfing flowers harbored less genetic diversity (π: R2 = .63, p = .01, n = 9 populations), and we found significant population structuring (FST = 0.48). Both geographic isolation and environmental factors played significant roles in predicting the observed genetic diversity: we found that corridors of suitable environments appear to facilitate gene flow between populations, and that environmental resistance is correlated with increased genetic distance between populations. Last, we integrated our genetic results with species distribution modeling to assess likely patterns of connectivity among our study populations. Our landscape and evolutionary genetic results suggest that T. perfoliata experienced a complex demographic and evolutionary history, particularly in the center of its distribution. As such, there is no singular mechanism driving this species' evolution. Together, our analyses support the hypothesis that the breeding system, geography, and environmental variables shape the patterns of diversity and connectivity of T. perfoliata in the US.

Abstract This research investigates the extent and causal mechanisms of genetic population divergence in a poorly flighted passerine, the North Island Rifleman or Titipounamu (Acanthisitta chloris granti). While this species has a historically widespread distribution, anthropogenic forest clearance has resulted in a highly fragmented current distribution. We conducted analyses of mitochondrial DNA (COI and Control Region) and 12 nuclear DNA microsatellites to test for population divergence and estimate times of divergence. DiyAbc and BioGeoBears were then used to assess likely past dispersal scenarios based on both mtDNA and nDNA. The results reveal several significantly divergent lineages across the North Island of New Zealand and indicate that some populations have been isolated for extensive periods of time (0.7–4.9 mya). Modeling indicated a dynamic history of population connectivity, with a drastic restriction in gene flow between three geographic regions, followed by a more recent re‐establishment of connectivity. Our analyses indicate the dynamic influence of key geological and climatological events on the distribution of genetic diversity in this species, including support for the genetic impact of old biogeographic boundaries such as the Taupo Line and Cockayne's Line, rather than recent anthropogenic habitat fragmentation. These findings present a rare example of an avian species with a genetic history more like that of flightless taxa and so provide new general insights into vicariant processes affecting populations of passerines with limited dispersal.

Abstract The climatic and geological changes that occurred during the Quaternary, particularly the fluctuations during the glacial and interglacial periods of the Pleistocene, shaped the population demography and geographic distribution of many species. These processes have been studied in several groups of organisms in the Northern Hemisphere, but their influence on the evolution of Neotropical montane species and ecosystems remains unclear. This study contributes to the understanding of the effect of climatic fluctuations during the late Pleistocene on the evolution of Andean mountain forests. First, we describe the nuclear and plastidic DNA patterns of genetic diversity, structure, historical demography, and landscape connectivity of Quercus humboldtii, which is a typical species in northern Andean montane forests. Then, these patterns were compared with the palynological and evolutionary hypotheses postulated for montane forests of the Colombian Andes under climatic fluctuation scenarios during the Quaternary. Our results indicated that populations of Q. humboldtii have high genetic diversity and a lack of genetic structure and that they have experienced a historical increase in connectivity from the last glacial maximum (LGM) to the present. Furthermore, our results showed a dramatic reduction in the effective population size followed by an expansion before the LGM, which is consistent with the results found by palynological studies, suggesting a change in dominance in Andean forests that may be related to ecological factors rather than climate change.

Abstract Aim Evolutionary history of natural populations can be confounded by human intervention such as the case of decorator worm species Diopatra (Onuphidae), which have a history of being transported through anthropogenic activities. Because they build tubes and act as ecosystem engineers, they can have a large impact on the overall ecosystem in which they occur. One conspicuous member, Diopatra biscayensis, which was only described in 2012, has a fragmented distribution that includes the Bay of Biscay and the Normanno‐Breton Gulf in the English Channel. This study explores the origin of these worms in the Normanno‐Breton region, which has been debated to either be the result of a historic range contraction from a relic continuous population or a more recent introduction. Location Northeastern Atlantic, the Bay of Biscay, and the Normanno‐Breton Gulf. Methods We utilized a RAD‐tag‐based SNP approach to create a reduced genomic data set to recover fine‐scale population structure and infer which hypothesis best describes the D. biscayensis biogeographic distribution. The reduced genomic data set was used to calculate standard genetic diversities and genetic differentiation statistics, and utilized various clustering analyses, including PCAs, DAPC, and admixture. Results Clustering analyses were consistent with D. biscayensis as a single population spanning the Bay of Biscay to the Normanno‐Breton Gulf in the English Channel, although unexpected genetic substructure was recovered from Arcachon Bay, in the middle of its geographic range. Consistent with a hypothesized introduction, the isolated Sainte‐Anne locality in the Normanno‐Breton Gulf was recovered to be a subset of the diversity found in the rest of the Bay of Biscay. Main conclusions These results are congruent with previous simulations that did not support connectivity from the Bay of Biscay to the Normanno‐Breton Gulf by natural dispersal. These genomic findings, with support from previous climatic studies, further support the hypothesis that D. biscayensis phylogeographic connectivity is the result of introductions, likely through the regions’ rich shellfish aquaculture, and not of a historically held range contraction.

Abstract In human‐dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244‐bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos‐Dinara populations. Total effective population size (Ne) was estimated to be 199 individuals, and total combined population size (NC) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe's low‐latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long‐term persistence of the species in the Balkans.

Life-history traits, especially the mode and duration of larval development, are expected to strongly influence the population connectivity and phylogeography of marine species. Comparative analysis of sympatric, closely related species with differing life histories provides the opportunity to specifically investigate these mechanisms of evolution but have been equivocal in this regard. Here, we sample two sympatric sea stars across the same geographic range in temperate waters of Australia. Using a combination of mitochondrial DNA sequences, nuclear DNA sequences, and microsatellite genotypes, we show that the benthic-developing sea star, Parvulastra exigua, has lower levels of within- and among-population genetic diversity, more inferred genetic clusters, and higher levels of hierarchical and pairwise population structure than Meridiastra calcar, a species with planktonic development. While both species have populations that have diverged since the middle of the second glacial period of the Pleistocene, most P. exigua populations have origins after the last glacial maxima (LGM), whereas most M. calcar populations diverged long before the LGM. Our results indicate that phylogenetic patterns of these two species are consistent with predicted dispersal abilities; the benthic-developing P. exigua shows a pattern of extirpation during the LGM with subsequent recolonization, whereas the planktonic-developing M. calcar shows a pattern of persistence and isolation during the LGM with subsequent post-Pleistocene introgression.

Abstract Shared ancestral polymorphism and introgression are two main causes of chloroplast DNA (cpDNA) haplotype sharing among closely related angiosperms. In this study, we explored the roles of these two processes in shaping the phylogeographic patterns of East Asian Cerris oaks by examining the geographic distributions of randomly and locally distributed shared haplotypes, which coincide with the expectations of shared ancestry and introgression, respectively. We sequenced 1340 bp of non‐coding cpDNA from Quercus acutissima (n = 418) and Q. chenii (n = 183) and compiled previously published sequence data of Q. variabilis (n = 439). The phylogenetic relationships among haplotypes were examined using a median‐joining network. The geographic patterns of interspecifically shared haplotypes were assessed to test whether nearby populations have a higher degree of interspecific cpDNA sharing than distant ones. We identified a total of 27 haplotypes that were grouped into three non‐species‐specific lineages with overlapping distributions. Ancestral haplotypes were extensively shared and randomly distributed across populations of the three species. Some young haplotypes were locally shared in mountainous areas that may have been shared refugia. The local exchange of cpDNA resulted in an excess of similar haplotypes between nearby populations. Our study demonstrated that the haplotype sharing pattern among East Asian Cerris oaks reflected the imprints of both shared ancestral polymorphism and introgression. This pattern was also associated with the relatively stable climates and complex landscapes in East Asia, which not only allowed the long‐term persistence of ancestral lineages but also connected the survived populations across refugia.

Abstract Rhinogobius candidianus is a freshwater goby distributed in north, northwest, west, and south Taiwan, but this species has been introduced to east Taiwan and became dominant. To investigate its native population genetic diversity and structure and evaluate the sources and diversity of translocated populations, the mitochondrial DNA control region and cytochrome b gene (1981 bp) from 220 specimens were analyzed. These results indicated that (1) the east populations originated from two sources in west Taiwan; (2) translocated populations exist in east Taiwan and south Taiwan; (3) many populations have likely been moved secondarily by human intervention; (4) the effective size of the populations had declined greatly; (5) within the native populations, the ancestral populations colonized Taiwan during the land bridge phase in the Pleistocene through north Taiwan; (6) the landform changes in Taiwan shaped the population structure; and (7) the landforms of the coastline during glaciation also shaped the native range. The low‐level genetic diversity, high population differentiation, and population decline greatly suggest the need for resource management and conservation interventions. Four clades (α–δ) should be managed as four distinct evolutionarily significant units, while the translocated populations should be managed as separate management units. Moreover, the translocated populations in east Taiwan should be evaluated and monitored carefully.

Abstract Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish‐water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double‐digest restriction‐site‐associated DNA sequencing (ddRAD‐Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present‐day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype–phenotype‐environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.

Abstract Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five‐prime region of cytochrome c oxidase subunit I (COI‐5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI‐5P lineages. The distribution of COI‐5P haplotypes may not demonstrate significant phylogeographic structure (NST > GST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI‐5P haplotypes. For both COI‐5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica's three major mitochondrial COI‐5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015–0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species’ distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.

Abstract Many Northeast (NE) Pacific fishes and invertebrates survived Pleistocene glaciations in northern refugia, but the extent that kelps survived in northern areas is uncertain. Here, we test the hypothesis that populations of sugar kelp (Saccharina latissima) persisted in the Gulf of Alaska during ice‐age maxima when the western margin of the Cordilleran ice sheet covered coastal areas around the NE Pacific Ocean. We estimated genetic diversities within and phylogeographical relationships among 14 populations along 2,800 km in the NE Pacific and Bering Sea with partial sequences of mitochondrial DNA 5′‐cytochrome oxidase subunit I (COI, bp = 624, n = 543), chloroplast DNA ribulose‐1,5‐bisphosphate carboxylase large subunit‐3′ (rbcL, bp = 735, n = 514), and 11 microsatellite loci. Concatenated sequences of rbcL and COI showed moderate levels of within‐population genetic diversity (mean h = 0.200) but substantial differences among populations (ΦST = 0.834, p

Abstract Temporal genetic studies of low‐dispersing organisms are rare. Marine invertebrates lacking a planktonic larval stage are expected to have lower dispersal, low gene flow, and a higher potential for local adaptation than organisms with planktonic dispersal. Leptasterias is a genus of brooding sea stars containing several cryptic species complexes. Population genetic methods were used to resolve patterns of fine‐scale population structure in central California Leptasterias species using three loci from nuclear and mitochondrial genomes. Historic samples (collected between 1897 and 1998) were compared to contemporary samples (collected between 2008 and 2014) to delineate changes in species distributions in space and time. Phylogenetic analysis of contemporary samples confirmed the presence of a bay‐localized clade and revealed the presence of an additional bay‐localized and previously undescribed clade of Leptasterias. Analysis of contemporary and historic samples indicates two clades are experiencing a constriction in their southern range limit and suggests a decrease in clade‐specific abundance at sites at which they were once prevalent. Historic sampling revealed a dramatically different distribution of diversity along the California coastline compared to contemporary sampling and illustrates the importance of temporal genetic sampling in phylogeographic studies. These samples were collected prior to significant impacts of Sea Star Wasting Disease (SSWD) and represent an in‐depth analysis of genetic structure over 117 years prior to the SSWD‐associated mass die‐off of Leptasterias.

Abstract Leiolepis ocellata is a lizard species distributing in topographically diverse habitats in northern Thailand. To explore its evolutionary history, 113 samples of L. ocellata were collected from 11 localities covering its distributional range in northern Thailand, and sequenced for mtDNA fragments (Cyt b and ND2). Pairwise comparisons across sampling localities yielded significant genetic differentiation (FST and Jost's D) but no clear pattern of isolation by distance could be demonstrated based on the Mantel test. Phylogenetic and network analyses highlighted six haplogroups. Their divergence times were estimated to occur during the Pleistocene, much more recent than major orogenic events affecting northern Thailand. Instead, the results suggested that lineage divergences, of particularly eastern and western haplogroups of the region, coincided with the major rivers in the region (Yom river and Ping river, respectively), indicating vicariance in response to riverine barriers. Furthermore, ecological niche modeling suggested an expansion of suitable habitats of L. ocellata, when LGM‐liked conditions. This expansion potentially facilitated their dispersal among adjacent localities leading to lineage diversification and genetic admixture, after the riverine divergence.

Abstract We conducted a comprehensive analysis of the phylogenetic, phylogeographic, and demographic relationships of Caspian cobra (Naja oxiana; Eichwald, 1831) populations based on a concatenated dataset of two mtDNA genes (cyt b and ND4) across the species' range in Iran, Afghanistan, and Turkmenistan, along with other members of Asian cobras (i.e., subgenus Naja Laurenti, 1768). Our results robustly supported that the Asiatic Naja are monophyletic, as previously suggested by other studies. Furthermore, N. kaouthia and N. sagittifera were recovered as sister taxa to each other, and in turn sister clades to N. oxiana. Our results also highlighted the existence of a single major evolutionary lineage for populations of N. oxiana in the Trans‐Caspian region, suggesting a rapid expansion of this cobra from eastern to western Asia, coupled with a rapid range expansion from east of Iran toward the northeast. However, across the Iranian range of N. oxiana, subdivision of populations was not supported, and thus, a single evolutionary significant unit is proposed for inclusion in future conservation plans in this region.

Abstract Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (FST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.

Abstract Ecological restoration often requires translocating plant material from distant sites. Importing suitable plant material is important for successful establishment and persistence. Yet, published guidelines for seed transfer are available for very few species. Accurately predicting how transferred plants will perform requires multiyear and multi‐environment field trials and comprehensive follow‐up work, and is therefore infeasible given the number of species used in restoration programs. Alternative methods to predict the outcomes of seed transfer are valuable for species without published guidelines. In this study, we analyzed the genetic structure of an important shrub used in ecological restoration in the Southern Rocky Mountains called alder‐leaf mountain mahogany (Cercocarpus montanus). We sequenced DNA from 1,440 plants in 48 populations across a broad geographic range. We found that genetic heterogeneity among populations reflected the complex climate and topography across which the species is distributed. We identified temperature and precipitation variables that were useful predictors of genetic differentiation and can be used to generate seed transfer recommendations. These results will be valuable for defining management and restoration practices for mountain mahogany.