Genetic consequences of climatic oscillations in the Quaternary

The role of Pleistocene dispersal in shaping species richness of sky island wintergreens from the Himalaya-Hengduan Mountains

In addition to topography and climate, biogeographic dispersal has been considered to influence plant diversity in the Himalaya-Hengduan Mountains (HHM), yet, the mode and tempo of sky island dispersal and its influence on species richness has been little explored. Through phylogenetic analysis of Gaultheria ser. Trichophyllae, a sky island alpine clade within the HHM, we test the hypothesis that dispersal has affected current local species richness. We inferred the dynamics of biogeographic dispersal with correlation tests on direction, distance, occurrence time, and regional species richness. We found that G. ser. Trichophyllae originated at the end of the Miocene and mostly dispersed toward higher longitudes (eastward). In particular, shorter intra-regional eastward dispersals and longer inter-regional westward dispersals were most frequently observed. We detected a prevalence of eastward intra-region dispersals in both glacial periods and interglacials. These dispersals may have been facilitated by the reorganization of paleo-drainages and monsoon intensification through time. We suggest that the timing of dispersal corresponding to glacial periods and the prevalence of intra-region dispersal, rather than dispersal frequency, most influenced the pattern of species richness of G. ser. Trichophyllae. This study facilitates a more comprehensive understanding of biodiversity in the sky islands within the HHM.

High dispersal ability versus migratory traditions: Fine-scale population structure and post-glacial colonisation in bar-tailed godwits

In migratory animals, high mobility may reduce population structure through increased dispersal and enable adaptive responses to environmental change, whereas rigid migratory routines predict low dispersal, increased structure, and limited flexibility to respond to change. We explore the global population structure and phylogeographic history of the bar-tailed godwit, Limosa lapponica, a migratory shorebird known for making the longest non-stop flights of any landbird. Using nextRAD sequencing of 14,318 single-nucleotide polymorphisms and scenario-testing in an Approximate Bayesian Computation framework, we infer that bar-tailed godwits existed in two main lineages at the last glacial maximum, when much of their present-day breeding range persisted in a vast, unglaciated Siberian-Beringian refugium, followed by admixture of these lineages in the eastern Palearctic. Subsequently, population structure developed at both longitudinal extremes: in the east, a genetic cline exists across latitude in the Alaska breeding range of subspecies L. l. baueri; in the west, one lineage diversified into three extant subspecies L. l. lapponica, taymyrensis, and yamalensis, the former two of which migrate through previously glaciated western Europe. In the global range of this long-distance migrant, we found evidence of both (1) fidelity to rigid behavioural routines promoting fine-scale geographic population structure (in the east) and (2) flexibility to colonise recently available migratory flyways and non-breeding areas (in the west). Our results suggest that cultural traditions in highly mobile vertebrates can override the expected effects of high dispersal ability on population structure, and provide insights for the evolution and flexibility of some of the world's longest migrations.

Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

Runs of homozygosity reveal contrasting histories of inbreeding across global lineages of the edible porcini mushroom, Boletus edulis

Inbreeding, the mating of individuals that are related through common ancestry, is of central importance in evolutionary and conservation biology due to its impacts on individual fitness and population dynamics. However, while advanced genomic approaches have revolutionised the study of inbreeding in animals, genomic studies of inbreeding are rare in plants and lacking in fungi. We investigated global patterns of inbreeding in the prized edible porcini mushroom Boletus edulis using 225 whole genomes from seven lineages distributed across the northern hemisphere. Genomic inbreeding was quantified using runs of homozygosity (ROHs). We found appreciable variation both among and within lineages, with some individuals having over 20% of their genomes in ROHs. Much of this variation could be explained by a combination of elevation and latitude, and to a lesser extent by predicted habitat suitability during the last glacial maximum. In line with this, the majority of ROHs were short, reflecting ancient common ancestry dating back approximately 200-1700 generations ago, while longer ROHs indicative of recent common ancestry (less than approximately 50 generations ago) were infrequent. Our study reveals the inbreeding legacy of major climatic events in a widely distributed forest mutualist, aligning with prevailing theories and empirical studies of the impacts of historical glaciation events on the dominant forest tree species of the northern hemisphere.

Montane diversification as a mechanism of speciation in neotropical butterflies

The mountains in the Atlantic Forest domain are environments that harbor a high biodiversity, including species adapted to colder climates that were probably influenced by the climatic variations of the Pleistocene. To understand the phylogeographic pattern and assess the taxonomic boundaries between two sister montane species, a genomic study of the butterflies Actinote mantiqueira and A. alalia (Nymphalidae: Acraeini) was conducted. Analyses based on partial sequences of the mitochondrial gene COI (barcode region) failed to recover any phylogenetic or genetic structure discriminating the two species or sampling localities. However, single nucleotide polymorphisms gathered using Genotyping-by-Sequencing provided a strong isolation pattern in all analyses (genetic distance, phylogenetic hypothesis, clustering analyses, and F ST statistics) which is consistent with morphology, separating all individuals of A. alalia from all populations of A. mantiqueira. The three sampled mountain ranges where A. mantiqueira populations occur-Serra do Mar, Serra da Mantiqueira, and Poços de Caldas Plateau-were identified as three isolated clusters. Paleoclimate simulations indicate that both species' distributions changed according to climatic oscillations in the Pleistocene period, with the two species potentially occurring in areas of lower altitude during glacial periods when compared to the interglacial periods (as the present). Besides, a potential path between their distribution through the Serra do Mar Mountain range was inferred. Therefore, the Pleistocene climatic fluctuation had a significant impact on the speciation process between A. alalia and A. mantiqueira, which was brought on by isolation at different mountain summits during interglacial periods, as shown by the modeled historical distribution and the observed genetic structure.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

Pushed waves, trailing edges, and extreme events: Eco-evolutionary dynamics of a geographic range shift in the owl limpet, Lottia gigantea

As climatic variation re-shapes global biodiversity, understanding eco-evolutionary feedbacks during species range shifts is of increasing importance. Theory on range expansions distinguishes between two different forms: "pulled" and "pushed" waves. Pulled waves occur when the source of the expansion comes from low-density peripheral populations, while pushed waves occur when recruitment to the expanding edge is supplied by high-density populations closer to the species' core. How extreme events shape pushed/pulled wave expansion events, as well as trailing-edge declines/contractions, remains largely unexplored. We examined eco-evolutionary responses of a marine invertebrate (the owl limpet, Lottia gigantea) that increased in abundance during the 2014-2016 marine heatwaves near the poleward edge of its geographic range in the northeastern Pacific. We used whole-genome sequencing from 19 populations across >11 degrees of latitude to characterize genomic variation, gene flow, and demographic histories across the species' range. We estimated present-day dispersal potential and past climatic stability to identify how contemporary and historical seascape features shape genomic characteristics. Consistent with expectations of a pushed wave, we found little genomic differentiation between core and leading-edge populations, and higher genomic diversity at range edges. A large and well-mixed population in the northern edge of the species' range is likely a result of ocean current anomalies increasing larval settlement and high-dispersal potential across biogeographic boundaries. Trailing-edge populations have higher differentiation from core populations, possibly driven by local selection and limited gene flow, as well as high genomic diversity likely as a result of climatic stability during the Last Glacial Maximum. Our findings suggest that extreme events can drive poleward range expansions that carry the adaptive potential of core populations, while also cautioning that trailing-edge extirpations may threaten unique evolutionary variation. This work highlights the importance of understanding how both trailing and leading edges respond to global change and extreme events.

Shared phylogeographic patterns and environmental responses of co-distributed soybean pests: Insights from comparative phylogeographic studies of Riptortus pedestris and Riptortus linearis in the subtropics of East Asia

Comparative phylogeographic studies of closely related species sharing co-distribution areas can elucidate the role of shared historical factors and environmental changes in shaping their phylogeographic pattern. The bean bugs, Riptortus pedestris and Riptortus linearis, which both inhabit subtropical regions in East Asia, are recognized as highly destructive soybean pests. Many previous studies have investigated the biological characteristics, pheromones, chemicals and control mechanisms of these two pests, but few studies have explored their phylogeographic patterns and underlying factors. In this study, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) dataset to investigate phylogeographic patterns and construct ecological niche models (ENM) for both Riptortus species. Our findings revealed similar niche occupancies and population genetic structures between the two species, with each comprising two phylogeographic lineages (i.e., the mainland China and the Indochina Peninsula clades) that diverged approximately 0.1 and 0.3 million years ago, respectively. This divergence likely resulted from the combined effects of temperatures variation and geographical barriers in the mountainous regions of Southwest China. Further demographic history and ENM analyses suggested that both pests underwent rapid expansion prior to the Last Glacial Maximum (LGM). Furthermore, ENM predicts a northward shift of both pests into new soybean-producing regions due to global warming. Our study indicated that co-distribution soybean pests with overlapping ecological niches and similar life histories in subtropical regions of East Asia exhibit congruent phylogeographic and demographic patterns in response to shared historical biogeographic drivers.

Genomic data revealed inbreeding despite a geographically connected stable effective population size since the Holocene in the protected Formosan Long-Arm Scarab beetle, Cheirotonus formosanus

Biodiversity conservation is a top priority in the face of global environmental change, and the practical restoration of biodiversity has emerged as a key objective. Nevertheless, the question of how to effectively contribute to biodiversity restoration and identify suitable systems for such efforts continues to present major challenges. By using genome-wide SNP data, our study revealed that populations from different mountain ranges of the Formosan Long-Arm Scarab beetle, a flagship species that receives strict protection, exhibited a single genetic cluster with no subdivision. Additionally, our result implied an association between the demographic history and historical fluctuations in climate and environmental conditions. Furthermore, we showed that, despite a stable and moderately sized effective population over recent history, all the individuals we studied exhibited signs of genetic inbreeding. We argued that the current practice of protecting the species as one evolutionarily significant unit remains the best conservation plan and that recent habitat change may have led to the pattern of significant inbreeding. We closed by emphasizing the importance of conservation genetic studies in guiding policy decisions and highlighting the potential of genomic data for identifying ideal empirical systems for genetic rescue, or assisted gene flow studies.

Environmental and Biogeographic Drivers behind Alpine Plant Thermal Tolerance and Genetic Variation

In alpine ecosystems, elevation broadly functions as a steep thermal gradient, with plant communities exposed to regular fluctuations in hot and cold temperatures. These conditions lead to selective filtering, potentially contributing to species-level variation in thermal tolerance and population-level genetic divergence. Few studies have explored the breadth of alpine plant thermal tolerances across a thermal gradient or the underlying genetic variation thereof. We measured photosystem heat (Tcrit-hot) and cold (Tcrit-cold) thresholds of ten Australian alpine species across elevation gradients and characterised their neutral genetic variation. To reveal the biogeographical drivers of present-day genetic signatures, we also reconstructed temporal changes in habitat suitability across potential distributional ranges. We found intraspecific variation in thermal thresholds, but this was not associated with elevation, nor underpinned by genetic differentiation on a local scale. Instead, regional population differentiation and considerable homozygosity within populations may, in part, be driven by distributional contractions, long-term persistence, and migrations following habitat suitability. Our habitat suitability models suggest that cool-climate-distributed alpine plants may be threatened by a warming climate. Yet, the observed wide thermal tolerances did not reflect this vulnerability. Conservation efforts should seek to understand variations in species-level thermal tolerance across alpine microclimates.

The phylogenetic relationship and demographic history of rhesus macaques (Macaca mulatta) in subtropical and temperate regions, China

Pleistocene climatic oscillations exerted significant influences on the genetic structure and demography of rhesus macaque (Macaca mulatta) in eastern China. However, the evolutionary history of rhesus macaques in subtropical and temperate China remained unclear and/or controversial. Herein, we analyzed the autosomes, mitochondrial genomes, and Y-chromosomes from 84 individuals of Chinese rhesus macaque. The results revealed that (1) all individuals were clustered into pan-west and pan-east genetic groups, which exhibited Shaanxi Province as the northernmost region of western dispersal route of rhesus macaques in China; (2) in subtropical and temperate China, rhesus macaques were divided into four lineages (TH, DB, HS, and QL), and their divergence times corresponded to the Penultimate Glaciation (300-130 kya) and Last Glaciation (70-10 kya), respectively; (3) the individuals from Mt. Taihangshan (TH) are closely related to individuals from Mt. Dabashan (DB) in the autosomal tree, rather than individuals from Mt. Huangshan (HS) as indicated by the mitogenome tree, which supports the hypothesis that the ancestral rhesus macaques radiated into Mt. Taihangshan from Mt. Huangshan via Mt. Dabashan; and (4) the demographic scenario of the four lineages showed the ancestral rhesus macaques bottleneck and expansion corresponding to the suitable habitat reduction and expansion, which confirmed they had experienced northward recolonization and southward retreat events from Mt. Huangshan area via Northern China Plain to Northernmost China along with Pleistocene glacial cycles. This study provides a new insight into understanding how Pleistocene glaciation has influenced faunal diversity in subtropical and temperate China, especially for those exhibiting differential patterns of sex dispersal.

Phylogeography and demographic history of macaques, fascicularis species group, in East Asia: Inferred from multiple genomic markers

Climate changes at larger scales have influenced dispersal and range shifts of many taxa in East Asia. The fascicularis species group of macaques is composed of four species and is widely distributed in Southeast and East Asia. However, its phylogeography and demographic histories are currently poorly understood. Herein, we assembled autosomal, mitogenome, and Y-chromosome data for 106 individuals, and combined them with 174 mtDNA dloop haplotypes of this species group, with particular focus on the demographic histories and dispersal routes of Macaca fuscata, M. cyclopis, and M. mulatta. The results showed: (1) three monophyletic clades for M. fuscata, M. cyclopis, and M. mulatta based on the multiple genomics analyses; (2) the disparate demographic trajectories of the three species after their split ∼1.0 Ma revealed that M. cyclopis and M. fuscata were derived from an ancestral M. mulatta population; (3) the speciation time of M. cyclopis was later than that of M. fuscata, and their divergence time occurred at the beginning of "Ryukyu Coral Sea Stage" (1.0-0.2 Ma) when the East China Sea land bridge was completely submerged by the sea level rose; and (4) the three parallel rivers (Nujiang, Lancangjiang, and Jinshajiang) of Southwestern China divided M. mulatta into Indian and Chinese genetic populations ∼200 kya. These results shed light on understanding not only the evolutionary history of the fascicularis species group but also the formation mechanism of faunal diversity in East Asia during the Pleistocene.

Parallel patterns of genetic diversity and structure in circumboreal species of the Sphagnum capillifolium complex

PREMISE

Shared geographical patterns of population genetic variation among related species is a powerful means to identify the historical events that drive diversification. The Sphagnum capillifolium complex is a group of closely related peat mosses within the Sphagnum subgenus Acutifolia and contains several circumboreal species whose ranges encompass both glaciated and unglaciated regions across the northern hemisphere. In this paper, we (1) inferred the phylogeny of subg. Acutifolia and (2) investigated patterns of population structure and genetic diversity among five circumboreal species within the S. capillifolium complex.

METHODS

We generated RAD sequencing data from most species of the subg. Acutifolia and samples from across the distribution ranges of circumboreal species within the S. capillifolium complex.

RESULTS

We resolved at least 14 phylogenetic clusters within the S. capillifolium complex. Five circumboreal species show some common patterns: One population system comprises plants in eastern North America and Europe, and another comprises plants in the Pacific Northwest or around the Beringian and Arctic regions. Alaska appears to be a hotspot for genetic admixture, genetic diversity, and sometimes endemic subclades.

CONCLUSIONS

Our results support the hypothesis that populations of five circumboreal species within the S. capillifolium complex survived in multiple refugia during the last glacial maximum. Long-distance dispersal out of refugia, population bottlenecks, and possible adaptations to conditions unique to each refugium could have contributed to current geographic patterns. These results indicate the important role of historical events in shaping the complex population structure of plants with broad distribution ranges.

Phylogeographic Pattern of the Assassin Bug Sycanus bifidus Inferred from Mitochondrial Genomes and Nuclear Genes

The assassin bug Sycanus bifidus has a wide distribution across southern China. This study explored its distribution and evolution by analyzing mitochondrial and nuclear ribosomal RNA genes, revealing how Pleistocene climate and geological changes shaped its phylogeography. We identified two main clades, A and B, that diverged in the Middle Pleistocene. Hainan Island's populations form a unique group within Clade A, suggesting that the Qiongzhou Strait served as a dispersal corridor during glaciation. Rising sea levels likely separated the Hainan population afterward. Ecological niche modeling showed that both populations have been viable since the last interglacial period, with demographic analyses indicating possible expansions during the Middle and Late Pleistocene, driven by favorable climates. This study highlights the significant effects of Pleistocene sea-level and climatic changes on the distribution and evolution of S. bifidus in China.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

Giraffe lineages are shaped by major ancient admixture events

Strong genetic structure has prompted discussion regarding giraffe taxonomy,1,2,3 including a suggestion to split the giraffe into four species: Northern (Giraffa c. camelopardalis), Reticulated (G. c. reticulata), Masai (G. c. tippelskirchi), and Southern giraffes (G. c. giraffa).4,5,6 However, their evolutionary history is not yet fully resolved, as previous studies used a simple bifurcating model and did not explore the presence or extent of gene flow between lineages. We therefore inferred a model that incorporates various evolutionary processes to assess the drivers of contemporary giraffe diversity. We analyzed whole-genome sequencing data from 90 wild giraffes from 29 localities across their current distribution. The most basal divergence was dated to 280 kya. Genetic differentiation, FST, among major lineages ranged between 0.28 and 0.62, and we found significant levels of ancient gene flow between them. In particular, several analyses suggested that the Reticulated lineage evolved through admixture, with almost equal contribution from the Northern lineage and an ancestral lineage related to Masai and Southern giraffes. These new results highlight a scenario of strong differentiation despite gene flow, providing further context for the interpretation of giraffe diversity and the process of speciation in general. They also illustrate that conservation measures need to target various lineages and sublineages and that separate management strategies are needed to conserve giraffe diversity effectively. Given local extinctions and recent dramatic declines in many giraffe populations, this improved understanding of giraffe evolutionary history is relevant for conservation interventions, including reintroductions and reinforcements of existing populations.

Range-Wide Phylogeography and Ecological Niche Modeling Provide Insights into the Evolutionary History of the Mongolian Racerunner (Eremias argus) in Northeast Asia

The Mongolian racerunner, Eremias argus, is a small lizard endemic to Northeast Asia that can serve as an excellent model for investigating how geography and past climate change have jointly influenced the evolution of biodiversity in this region. To elucidate the processes underlying its diversification and demography, we reconstructed the range-wide phylogeographic pattern and evolutionary trajectory, using phylogenetic, population genetic, landscape genetic, Bayesian phylogeographic reconstruction and ecological niche modeling approaches. Phylogenetic analyses of the mtDNA cyt b gene revealed eight lineages that were unbounded by geographic region. The genetic structure of E. argus was mainly determined by geographic distance. Divergence dating indicated that E. argus and E. brenchleyi diverged during the Mid-Pliocene Warm Period. E. argus was estimated to have coalesced at~0.4351 Ma (Marine Isotope Stage 19). Bayesian phylogeographic diffusion analysis revealed out-of-Inner Mongolia and rapid colonization events from the end of the Last Interglacial to the Last Glacial Maximum, which is consistent with the expanded suitable range of the Last Glacial Maximum. Pre-Last Glacial Maximum growth of population is presented for most lineages of E. argus. The Glacial Maximum contraction model and the previous multiple glacial refugia hypotheses are rejected. This may be due to an increase in the amount of climatically favorable habitats in Northeast Asia. Furthermore, E. argus barbouri most likely represents an invalid taxon. The present study is the first to report a range-wide phylogeography of reptiles over such a large region in Northeast Asia. Our results make a significant contribution towards understanding the biogeography of the entire Northeast Asia.

Response of Extremely Small Populations to Climate Change-A Case of Trachycarpus nanus in Yunnan, China

Climate change affects the geographical distribution of plant species. Rare Trachycarpus nanus with a narrow distribution range, high medicinal value and extremely small population is facing increasing extinction risks under global climate change. In this study, 96 recorded occurrences and 23 environmental factors are used to predict the potential suitable area of T. nanus based on the optimized MaxEnt (3.4.4) model and ArcGIS (10.7) software. The results show that when the parameters are FC = LQ and RM = 1, the MaxEnt model is optimal and AUC = 0.946. The distribution patterns were predicted in the past, present, and four future phases, i.e., 2021-2040 (2030), 2041-2060 (2050), 2061-2080 (2070), and 2081-2100 (2090). The main factors are the annual precipitation (bio12), mean temperature of the coldest quarter (bio11), temperature seasonality (bio4), precipitation of the wettest quarter (bio16), and isothermality (bio3). The potential distribution of T. nanus is primarily concentrated in central Chuxiong, encompassing a total potential suitable area of 5.65 × 104 km2. In historical periods, the total habitat area is smaller than that in the present. In the future, the potential suitable area is generally increased. The centroid analysis shows that T. nanus will move to a high-altitude area and to the southeast. But its dispersal capacity may not keep up with the climate change rate. Therefore, additional protection sites for this species should be appropriately established and the habitat connectivity should be enhanced.

Genomic analysis reveals deep population divergence in the water snake Trimerodytes percarinatus (Serpentes, Natricidae)

Although several phylogeographic studies of Asian snakes have been conducted, most have focused on pitvipers, with non-venomous snakes, such as colubrids or natricids, remaining poorly studied. The Chinese keelback water snake (Trimerodytes percarinatus Boulenger) is a widespread, semiaquatic, non-venomous species occurring in China and southeastern Asia. Based on mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, we explored the population genetic structure, genetic diversity, and evolutionary history of this species. MtDNA-based phylogenetic analysis showed that T. percarinatus was composed of five highly supported and geographically structured lineages. SNP-based phylogenetic analysis, principal component analysis, and population structure analysis consistently revealed four distinct, geographically non-overlapping lineages, which was different from the mtDNA-based analysis in topology. Estimation of divergence dates and ancestral area of origin suggest that T. percarinatus originated ~12.68 million years ago (95% highest posterior density: 10.36-15.96 Mya) in a region covering southwestern China and Vietnam. Intraspecific divergence may have been triggered by the Qinghai-Xizang Plateau uplift. Population demographics and ecological niche modeling indicated that the effective population size fluctuated during 0.5 Mya and 0.002 Mya. Based on the data collected here, we also comment on the intraspecific taxonomy of T. percarinatus and question the validity of the subspecies T. p. suriki.

Where east meets west: Phylogeography of the high Arctic North American brant goose

Genetic variation in Arctic species is often influenced by vicariance during the Pleistocene, as ice sheets fragmented the landscape and displaced populations to low- and high-latitude refugia. The formation of secondary contact or suture zones during periods of ice sheet retraction has important consequences on genetic diversity by facilitating genetic connectivity between formerly isolated populations. Brant geese (Branta bernicla) are a maritime migratory waterfowl (Anseriformes) species that almost exclusively uses coastal habitats. Within North America, brant geese are characterized by two phenotypically distinct subspecies that utilize disjunct breeding and wintering areas in the northern Pacific and Atlantic. In the Western High Arctic of Canada, brant geese consist of individuals with an intermediate phenotype that are rarely observed nesting outside this region. We examined the genetic structure of brant geese populations from each subspecies and areas consisting of intermediate phenotypes using mitochondrial DNA (mtDNA) control region sequence data and microsatellite loci. We found a strong east-west partition in both marker types consistent with refugial populations. Within subspecies, structure was also observed at mtDNA while microsatellite data suggested the presence of only two distinct genetic clusters. The Western High Arctic (WHA) appears to be a secondary contact zone for both Atlantic and Pacific lineages as mtDNA and nuclear genotypes were assigned to both subspecies, and admixed individuals were observed in this region. The mtDNA sequence data outside WHA suggests no or very restricted intermixing between Atlantic and Pacific wintering populations which is consistent with published banding and telemetry data. Our study indicates that, although brant geese in the WHA are not a genetically distinct lineage, this region may act as a reservoir of genetic diversity and may be an area of high conservation value given the potential of low reproductive output in this species.

C4 photosynthesis provided an immediate demographic advantage to populations of the grass Alloteropsis semialata

C4 photosynthesis is a key innovation in land plant evolution, but its immediate effects on population demography are unclear. We explore the early impact of the C4 trait on the trajectories of C4 and non-C4 populations of the grass Alloteropsis semialata. We combine niche models projected into paleoclimate layers for the last 5 million years with demographic models based on genomic data. The initial split between C4 and non-C4 populations was followed by a larger expansion of the ancestral C4 population, and further diversification led to the unparalleled expansion of descendant C4 populations. Overall, C4 populations spread over three continents and achieved the highest population growth, in agreement with a broader climatic niche that rendered a large potential range over time. The C4 populations that remained in the region of origin, however, experienced lower population growth, rather consistent with local geographic constraints. Moreover, the posterior transfer of some C4-related characters to non-C4 counterparts might have facilitated the recent expansion of non-C4 populations in the region of origin. Altogether, our findings support that C4 photosynthesis provided an immediate demographic advantage to A. semialata populations, but its effect might be masked by geographic contingencies.

Ecological characteristics explain neutral genetic variation of three coastal sparrow species

Eco-phylogeographic approaches to comparative population genetic analyses allow for the inclusion of intrinsic influences as drivers of intraspecific genetic structure. This insight into microevolutionary processes, including changes within a species or lineage, provides better mechanistic understanding of species-specific interactions and enables predictions of evolutionary responses to environmental change. In this study, we used single nucleotide polymorphisms (SNPs) identified from reduced representation sequencing to compare neutral population structure, isolation by distance (IBD), genetic diversity and effective population size (Ne) across three closely related and co-distributed saltmarsh sparrow species differing along a specialization gradient-Nelson's (Ammospiza nelsoni subvirgata), saltmarsh (A. caudacuta) and seaside sparrows (A. maritima maritima). Using an eco-phylogeographic lens within a conservation management context, we tested predictions about species' degree of evolutionary history and ecological specialization to tidal marshes, habitat, current distribution and population status on population genetic metrics. Population structure differed among the species consistent with their current distribution and habitat factors, rather than degree of ecological specialization: seaside sparrows were panmictic, saltmarsh sparrows showed hierarchical structure and Nelson's sparrows were differentiated into multiple, genetically distinct populations. Neutral population genetic theory and demographic/evolutionary history predicted patterns of genetic diversity and Ne rather than degree of ecological specialization. Patterns of population variation and evolutionary distinctiveness (Shapely metric) suggest different conservation measures for long-term persistence and evolutionary potential in each species. Our findings contribute to a broader understanding of the complex factors influencing genetic variation, beyond specialist-generalist status and support the role of an eco-phylogeographic approach in population and conservation genetics.

Recent genetic, phenetic and ecological divergence across the Mesoamerican highlands: a study case with Diglossa baritula (Aves: Thraupidae)

The topographical, geological, climatic and biodiversity complexity of Mesoamerica has made it a primary research focus. The Mesoamerican highlands is a region with particularly high species richness and within-species variation. The Cinnamon-bellied Flowerpiercer, Diglossa baritula (Wagler, 1832), is a species endemic to the Mesoamerican highlands, with three allopatric subspecies currently recognized. To characterize divergence within this species, we integrated genomics, morphology, coloration and ecological niche modeling approaches, obtained from sampling individuals across the entire geographic distribution of the species. Our results revealed a clear genomic divergence between the populations to the east versus the west of the Isthmus of Tehuantepec. In contrast to the genomic results, morphology and coloration analyses showed intermediate levels of differentiation, indicating that population groups within D. baritula have probably been under similar selective pressures. Our morphology results indicated that the only sexually dimorphic morphological variable is the wing chord, with males having a longer wing chord than females. Finally, ecological data indicated that there are differences in ecological niche within D. baritula. Our data suggest that D. baritula could contain two or more incipient species at the intermediate phase of the speciation continuum. These results highlight the importance of the geographical barrier of the Isthmus of Tehuantepec and Pleistocene climatic events in driving isolation and population divergence in D. baritula. The present investigation illustrates the speciation potential of the D. baritula complex and the capacity of Mesoamerican highlands to create cryptic biodiversity and endemism.

Responses of Three Pedicularis Species to Geological and Climatic Changes in the Qinling Mountains and Adjacent Areas in East Asia

The Qinling Mountains in East Asia serve as the geographical boundary between the north and south of China and are also indicative of climatic differences, resulting in rich ecological and species diversity. However, few studies have focused on the responses of plants to geological and climatic changes in the Qinling Mountains and adjacent regions. Therefore, we investigated the evolutionary origins and phylogenetic relationships of three Pedicularis species in there to provide molecular evidence for the origin and evolution of plant species. Ecological niche modeling was used to predict the geographic distributions of three Pedicularis species during the last interglacial period, the last glacial maximum period, and current and future periods, respectively. Furthermore, the distribution patterns of climate fluctuations and the niche dynamics framework were used to assess the equivalence or difference of niches among three Pedicularis species. The results revealed that the divergence of three Pedicularis species took place in the Miocene and Holocene periods, which was significantly associated with the large-scale uplifts of the Qinling Mountains and adjacent regions. In addition, the geographic distributions of three Pedicularis species have undergone a northward migration from the past to the future. The most important environmental variables affecting the geographic distributions of species were the mean diurnal range and annual mean temperature range. The niche divergence analysis suggested that the three Pedicularis species have similar ecological niches. Among them, P. giraldiana showed the highest niche breadth, covering nearly all of the climatic niche spaces of P. dissecta and P. bicolor. In summary, this study provides novel insights into the divergence and origins of three Pedicularis species and their responses to climate and geological changes in the Qinling Mountains and adjacent regions. The findings have also provided new perspectives for the conservation and management of Pedicularis species.

Ensemble species distribution modeling and multilocus phylogeography provide insight into the spatial genetic patterns and distribution dynamics of a keystone forest species, Quercus glauca

BACKGROUND

Forests are essential for maintaining species diversity, stabilizing local and global climate, and providing ecosystem services. Exploring the impact of paleogeographic events and climate change on the genetic structure and distribution dynamics of forest keystone species could help predict responses to future climate change. In this study, we combined an ensemble species distribution model (eSDM) and multilocus phylogeography to investigate the spatial genetic patterns and distribution change of Quercus glauca Thunb, a keystone of East Asian subtropical evergreen broad-leaved forest.

RESULTS

A total of 781 samples were collected from 77 populations, largely covering the natural distribution of Q. glauca. The eSDM showed that the suitable habitat experienced a significant expansion after the last glacial maximum (LGM) but will recede in the future under a general climate warming scenario. The distribution centroid will migrate toward the northeast as the climate warms. Using nuclear SSR data, two distinct lineages split between east and west were detected. Within-group genetic differentiation was higher in the West than in the East. Based on the identified 58 haplotypes, no clear phylogeographic structure was found. Populations in the Nanling Mountains, Wuyi Mountains, and the southwest region were found to have high genetic diversity.

CONCLUSIONS

A significant negative correlation between habitat stability and heterozygosity might be explained by the mixing of different lineages in the expansion region after LGM and/or hybridization between Q. glauca and closely related species. The Nanling Mountains may be important for organisms as a dispersal corridor in the west-east direction and as a refugium during the glacial period. This study provided new insights into spatial genetic patterns and distribution dynamics of Q. glauca.

A renewed glance at the Palearctic golden eagle: Genetic variation in space and time

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.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Cryptic divergences and repeated hybridizations within the endangered "living fossil" dove tree (Davidia involucrata) revealed by whole genome resequencing

The identification and understanding of cryptic intraspecific evolutionary units (lineages) are crucial for planning effective conservation strategies aimed at preserving genetic diversity in endangered species. However, the factors driving the evolution and maintenance of these intraspecific lineages in most endangered species remain poorly understood. In this study, we conducted resequencing of 77 individuals from 22 natural populations of Davidia involucrata, a "living fossil" dove tree endemic to central and southwest China. Our analysis revealed the presence of three distinct local lineages within this endangered species, which emerged approximately 3.09 and 0.32 million years ago. These divergence events align well with the geographic and climatic oscillations that occurred across the distributional range. Additionally, we observed frequent hybridization events between the three lineages, resulting in the formation of hybrid populations in their adjacent as well as disjunct regions. These hybridizations likely arose from climate-driven population expansion and/or long-distance gene flow. Furthermore, we identified numerous environment-correlated gene variants across the total and many other genes that exhibited signals of positive evolution during the maintenance of two major local lineages. Our findings shed light on the highly dynamic evolution underlying the remarkably similar phenotype of this endangered species. Importantly, these results not only provide guidance for the development of conservation plans but also enhance our understanding of evolutionary past for this and other endangered species with similar histories.

Demographic patterns of two related desert shrubs with overlapping distributions in response to past climate changes

Numerous studies have revealed that past geological events and climatic fluctuations had profoundly affected the genetic structure and demographic patterns of species. However, related species with overlapping ranges may have responded to such environmental changes in different ways. In this study, we compared the genetic structure and population dynamics of two typical desert shrubs with overlapping distributions in northern China, Nitraria tangutorum and Nitraria sphaerocarpa, based on chloroplast DNA (cpDNA) variations and species distribution models. We sequenced two cpDNA fragments (trnH-trnA and atpH-atpI) in 633 individuals sampled from 52 natural populations. Twenty-four chlorotypes, including eight rare chlorotypes, were identified, and a single dominant haplotype (H4) widely occurred in the entire geographical ranges of the two species. There were also a few distinctive chlorotypes fixed in different geographical regions. Population structure analyses suggested that the two species had significantly different levels of total genetic diversity and interpopulation differentiation, which was highly likely correlated with the special habitat preferences of the two species. A clear phylogeographic structure was identified to exist among populations of N. sphaerocarpa, but not exist for N. tangutorum. The neutral tests, together with the distribution of pairwise differences revealed that N. tangutorum experienced a sudden demographic expansion, and its expansion approximately occurred between 21 and 7 Kya before present, while a rapid range expansion was not identified for N. sphaerocarpa. The ecological niche modeling (ENM) analysis indicated that the potential ranges of two species apparently fluctuated during the past and present periods, with obvious contraction in the Last Glacial Maximum (LGM) and recolonization in the present, respectively, comparing to the Last Interglacial (LIG). These findings suggest that the two species extensively occurred in the Northwest of China before the Quaternary, and the current populations of them originated from a few separated glacial refugia following their habitat fragmentation in the Quarternary. Our results provide new insights on the impact of past geological and climatic fluctuations on the population dynamics of desert plants in northwestern China, and further enforce the hypothesis that there were several independent glacial refugia for these species during the Quaternary glaciations.

Regional environmental differences significantly affect the genetic structure and genetic differentiation of Carpinus tientaiensis Cheng, an endemic and extremely endangered species from China

Differences in topography and environment greatly affect the genetic structure and genetic differentiation of species, and endemic or endangered species with limited geographic ranges seem to be more sensitive to changes in climate and other environmental factors. The complex topography of eastern China is likely to affect genetic differentiation of plants there. Carpinus tientaiensis Cheng is a native and endangered plants from China, and exploring its genetic diversity has profound significance for protection and the collection of germplasm resources. Based on AFLP markers, this study found that C. tientaiensis has low genetic diversity, which mainly came from within populations, while Shangshantou and Tiantai Mountain populations have relatively high genetic diversity. The Nei genetic distance was closely related to geographical distance, and temperature and precipitation notablely affected the genetic variation and genetic differentiation of C. tientaiensis. Based on cpDNA, this study indicated that C. tientaiensis exhibits a moderate level of genetic diversity, and which mainly came from among populations, while Tiantai Mountain population have the highest genetic diversity. It demonstrated that there was genetic differentiation between populations, which can be divided into two independent geographical groups, but there was no significant phylogeographic structure between them. The MaxEnt model showed that climate change significantly affects its distribution, and the suitable distribution areas in Zhejiang were primarily divided into two regions, eastern Zhejiang and southern Zhejiang, and there was niche differentiation in its suitable distribution areas. Therefore, this study speculated that the climate and the terrain of mountains and hills in East China jointly shape the genetic structure of C. tientaiensis, which gived rise to an obvious north-south differentiation trend of these species, and the populations located in the hilly areas of eastern Zhejiang and the mountainous areas of southern Zhejiang formed two genetic branches respectively.

At a crossroads: Genetic lineages and dispersal routes of Morimusasper (Sulzer, 1776) s.l. (Coleoptera, Cerambycidae) in Bulgaria

The present study fills a knowledge gap in the distribution and genetic variation of Morimus populations in the Balkans, by studiyng the representatives of the genus in Bulgaria - M.asperfunereus Mulsant, 1862, M.verecundusbulgaricus Danilevsky, 2016 and M.orientalis Reitter, 1894. Additional information is provided for Albania and northern Greece. The mitochondrial cytochrome C oxidase subunit I (COI) marker and the nuclear internal transcribed spacer 2 (ITS2) were used for the genetic analyses. Three of the previously-defined mitochondrial lineages (Lb/HgA, L2 and L3) were detected in Bulgaria, as well as a new lineage (Str) from the Strandzha Mountains (south-eastern Bulgaria). A total of 24 distinct haplotypes, 20 of them in Bulgaria, were found. Bulgarian populations of Morimus demonstrated relatively high nucleotide diversity. The L3 COI lineage was confirmed as the most diverse and frequent in the Balkans. The L3 lineage is dominant in most of Bulgaria, but was not identified in the easternmost parts near the Black Sea coast, where the L2 and Str lineages were found. New data highlighted two dispersal routes of the L2 mitochondrial lineage on the Balkan Peninsula: 1) northwards along the Black Sea coast and 2) westwards, across the Balkans where only disjunct populations remain. North-western Bulgaria seems to be the eastern limit of the basal lineage Lb/HgA distribution. Our results show high levels of genetic exchange between most of the mitochondrially defined lineages, yet some of the easternmost populations probably remained isolated for comparatively longer periods.

Co-occurrence of dual lineages within Simulium (Gomphostilbia) atratum De Meijere in the Indonesian Archipelago along Wallace's Line

Mitochondrial cytochrome c oxidase subunit I (COI) sequences were utilized to infer the population genetic structure of Simulium (Gomphostilbia) atratum De Meijere, an endemic simulid species to Indonesia. Both median-joining haplotype network and maximum-likelihood tree revealed two genetic lineages (A and B) within the species, with an overlap distribution in Lombok, which is situated along Wallace's line. Genetic differentiation and gene flow with varying frequencies (FST = 0.02-0.967; Nm = 0.01-10.58) were observed between populations of S. (G.) atratum, of which population pairs of different lineages showed high genetic differentiation. Notably, the high genetic distance of up to 5.92 % observed within S. (G.) atratum in Lombok was attributed to the existence of two genetically distinct lineages. The co-occurrence of distinct lineages in Lombok indicated that Wallace's line did not act as faunistic border for S. (G.) atratum in the present study. Moreover, both lineages also exhibited unimodal distributions and negative values of neutrality tests, suggesting a pattern of population expansion. The expansion and divergence time estimation suggested that the two lineages of S. (G.) atratum diverged and expanded during the Pleistocene era in Indonesia.

Genetic consequences of Last Glacial-Holocene changes in snowfall regime in Arnica mallotopus populations: A plant confined to heavy-snow areas of Japan

PREMISE

Snow is an important environmental factor affecting plant distribution. Past changes in snowfall regimes may have controlled the demographies of snow-dependent plants. However, our knowledge of changes in the distribution and demographies of such plants is limited because of the lack of fossil records.

METHODS

Population genetic and landscape genetic analyses were used to investigate the response of population dynamics of Arnica mallotopus (Asteraceae)-a plant confined to heavy-snow areas of Japan-to changes in snowfall regimes from the Last Glacial Period to the Holocene.

RESULTS

The population genetic analysis suggested that the four geographic lineages diverged during the Last Glacial Period. The interaction between reduced snowfall and lower temperatures during this period likely triggered population isolation in separate refugia. Subpopulation differentiation in the northern group was lower than in the southern group. Our ecological niche model predicted that the current distribution was patchy in the southern region; that is, the populations were isolated by topologically flat and climatically unsuitable lowlands. The landscape genetic analysis suggested that areas with little snowfall acted as barriers to the Holocene expansion of species distribution and continued limiting gene flow between local populations.

CONCLUSIONS

These findings indicate that postglacial population responses vary among regions and are controlled by environmental and geographic factors. Thus, changes in snowfall regime played a major role in shaping the distribution and genetic structure of the snow-dependent plant.

The influence of climate oscillations and geological events on population differentiation of Camponotus japonicus in the Chinese mainland

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.

Mapping endemic freshwater fish richness to identify high-priority areas for conservation: An ecoregion approach

Freshwater ecosystems are experiencing accelerating global biodiversity loss. Thus, knowing where these unique ecosystems' species richness reaches a peak can facilitate their conservation planning. By hosting more than 290 freshwater fishes, Iran is a major freshwater fish hotspot in the Middle East. Considering the accelerating rate of biodiversity loss, there is an urgent need to identify species-rich areas and understand the mechanisms driving biodiversity distribution. In this study, we gathered distribution records of all endemic freshwater fishes of Iran (85 species) to develop their richness map and determine the most critical drivers of their richness patterns from an ecoregion approach. We performed a generalized linear model (GLM) with quasi-Poisson distribution to identify contemporary and historical determinants of endemic freshwater fish richness. We also quantified endemic fish similarity among the 15 freshwater ecoregions of Iran. Results showed that endemic freshwater fish richness is highest in the Zagros Mountains while a moderate level of richness was observed between Zagros and Alborz Mountains. High, moderate, and low richness of endemic freshwater fish match with Upper Tigris & Euphrates, Namak, and Kavir & Lut Deserts ecoregions respectively. Kura - South Caspian Drainages and Caspian Highlands were the most similar ecoregions and Orumiyeh was the most unique ecoregion according to endemic fish presence. Precipitation and precipitation change velocity since the Last Glacial Maximum were the most important predictors of endemic freshwater fish richness. Areas identified to have the highest species richness have high priority for the conservation of freshwater fish in Iran, therefore, should be considered in future protected areas development.

Evolutionary history of a cold-adapted limnephilid caddisfly: Effects of climate change and topography on genetic structure

The distribution of organisms is influenced by complex factors such as the phylogenetic evolutionary histories of species, the physiological and ecological characteristics of organisms, climate, and geographical and geohistorical features. In this study, we focused on a caddisfly, Asynarchus sachalinensis (Trichoptera: Limnephilidae), which has adapted to cold habitats. From phylogeographic analyses based on the mitochondrial DNA (mtDNA) cytochrome c oxidase subunit I (COI) and 16S rRNA regions and the nuclear DNA (nDNA) 18S rRNA, 28S rRNA, carbamoyl-phosphate synthetase (CAD), elongation factor-1 alpha (EF1-α), and RNA polymerase II (POLII) regions, two distinct genetic clades were detected. Clade I was shown to be widely distributed from Sakhalin to Honshu, whereas Clade II was only distributed within Honshu. The distributions of these clades overlapped in Honshu. The habitats were located at relatively lower altitudes for Clade I and higher altitudes for Clade II. The divergence time of these clades was estimated to be during the Pleistocene, indicating that repeated climatic changes facilitated distributional shifts. Haplotype network and demographic analyses based on the mtDNA COI region showed contrasting genetic structures in the two clades. It was indicated that the population sizes of Clade I had expanded rapidly in a recent period, whereas Clade II had maintained stable population sizes. The habitats of Clade II were typically isolated and scattered at high altitudes, resulting in restricted migration and dispersal because of their discontinuous "Sky Island" habitats. The habitats of Clade I were located at relatively low altitudes, and it was assumed that the populations were continuous, which resulted in a higher frequency of migration and dispersal between populations. Thus, differences in the spatial scale of the adapted habitats of each clade may have resulted in different patterns of population connectivity and fragmentation associated with repeated climatic changes during the Pleistocene. Our study provided new insight into the distributional patterns of cold-adapted aquatic insects in the Japanese Archipelago. Furthermore, the distributional shifts predicted by ecological niche modeling under future climatic change conditions were different for each clade. Therefore, different principles are required in the assessment of each clade to predict temporal changes in their distributions.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Differential distribution shifts in two subregions of East Asian subtropical evergreen broadleaved forests-a case of Magnoliaceae

Aim

East Asian subtropical evergreen broad-leaved forests (EBLFs) are composed of western and eastern subregions with different topographical and environmental conditions. The distribution shifts over time of plants in the two subregions are predicted to be different, but the difference has seldom been investigated.

Methods

Potential distributions of 53 Magnoliaceae species (22 in the western and 31 in the eastern subregion) during the last glacial maximum (LGM), present, and the 2070s were predicted using MaxEnt based on 58 environmental variables. The changes in the distribution range size and centroid over time were analyzed. Species-level potential habitats were overlaid to uncover species diversity distribution, and the distributions over time were overlaid to discover long-term refugia.

Results

At present, the potential distributions are significantly larger than those shown by the occurrence points. During the LGM, 20/22 species in the western subregion experienced increases in range size through downwards and southward migrations, while decreases in range size in the eastern subregion (27/31 species) were accompanied by northward and eastward migrations. In the future, range size declines and northward shifts will both be found; northwestward shifts will exist in most (20/22 species) species in the western subregion, while both northwest- and northeastward shifts will occur in the eastern subregion. The diversity hotspots experienced a slight southward shift in the past and upwards to the mountain region in the future in the western subregion; in the eastern subregion, shrinks occurred in eastern China in the past and shrinks were shown in all regions in the future. Long-term refugia-preserving diversity was found in the mountains across the entire EBLFs region.

Main conclusions

Significant differences in distribution shifts from past to present and similar distribution shifts from present to future are revealed in the two subregions. Species diversity in both subregions experienced no significant shifts from past to future, and Magnoliaceae plants could be preserved in mountainous regions throughout the EBLFs.

Genetic Structure and Population History of the Zaisan Toad-Headed Agama (Phrynocephalus melanurus) Inferred from Mitochondrial DNA

The agamid lizard Phrynocephalus melanurus is restricted to Northwest China (Dzungar Basin) and the adjacent Eastern Kazakhstan (Zaisan and Alakol basins). To elucidate the phylogeography of P. melanurus, we obtained the mitochondrial DNA COI segments of 175 sampled lizards from 44 localities across the whole distribution. Phylogenetic analyses revealed two main Clades comprising five geographically structured lineages (I, IIa, IIb1, IIb2, and IIb3) that fit an isolation-by-distance (IBD) model. The divergence from the most recent common ancestor was dated to ~1.87 million years ago (Ma). Demographic analyses demonstrated lineage-specific response to past climate change: stable population for Clade I, Subclade IIb1; past population expansion for IIb3 since 0.18 Ma, respectively. Bayesian phylogeographic diffusion analyses detected initial spreading at the Saur Mount vicinity, approximately 1.8 Ma. Historical species distribution model (SDM) projected expansion of the suitable habitat in the last interglacial and shift and contraction in the last glacial maximum and Holocene epochs. The SDM predicted a drastic reduction in suitable area throughout the range as a response to future climate change. Our findings suggest that the evolution of P. melanurus followed a parapatric divergence with subsequent dispersal and adaptation to cold and dry environments during the Quaternary. Overall, this work improves our understanding of the lineage diversification and population dynamics of P. melanurus, providing further insights into the evolutionary processes that occurred in Northwest China and adjacent Eastern Kazakhstan.

Conservation genetics and potential geographic distribution modeling of Corybas taliensis, a small 'sky Island' orchid species in China

BACKGROUND

Corybas taliensis is an endemic species of sky islands in China. Its habitat is fragile and unstable, and it is likely that the species is threatened. However, it is difficult to determine the conservation priority or unit without knowing the genetic background and the overall distribution of this species. In this study, we used double digest restriction-site associated DNA-sequencing (ddRAD-seq) to investigate the conservation genomics of C. taliensis. At the same time, we modeled the extent of suitable habitat for C. taliensis in present and future (2030 and 2090) habitat using the maximum-entropy (MaxEnt) model.

RESULTS

The results suggested that the related C. fanjingshanensis belongs to C. taliensis and should not be considered a separate species. All the sampling locations were divided into three genetic groups: the Sichuan & Guizhou population (SG population), the Hengduan Mountains population (HD population) and Himalayan population (HM population), and we found that there was complex gene flow between the sampling locations of HD population. MT was distinct genetically from the other sampling locations due to the unique environment in Motuo. The genetic diversity (π, He) of C. taliensis was relatively high, but its contemporary effective population size (Ne) was small. C. taliensis might be currently affected by inbreeding depression, although its large population density may be able to reduce the effect of this. The predicted areas of suitable habitat currently found in higher mountains will not change significantly in the future, and these suitable habitats are predicted to spread to other higher mountains under future climate change. However, suitable habitat in relatively low altitude areas may disappear in the future. This suggests that C. taliensis will be caught in a 'summit trap' in low altitude areas, however, in contrast, the high altitude of the Himalaya and the Hengduan Mountains are predicted to act as 'biological refuges' for C. taliensis in the future.

CONCLUSIONS

These results not only provide a new understanding of the genetic background and potential resource distribution of C. taliensis, but also lay the foundation for its conservation and management.

Patterns of speciation in a parapatric pair of Saturnia moths as revealed by target capture

The focus of this study has been to understand the evolutionary relationships and taxonomy of a widely distributed parapatric species pair of wild silk moths in Europe: Saturnia pavonia and Saturnia pavoniella (Lepidoptera: Saturniidae). To address species delimitation in these parapatric taxa, target enrichment and mtDNA sequencing was employed alongside phylogenetic, admixture, introgression, and species delimitation analyses. The dataset included individuals from both species close to and farther away from the contact zone as well as two hybrids generated in the lab. Nuclear markers strongly supported both S. pavonia and S. pavoniella as two distinct species, with hybrids forming a sister group to S. pavoniella. However, the Maximum Likelihood (ML) tree generated from mtDNA sequencing data presented a different picture, showing both taxa to be phylogenetically intermixed. This inconsistency is likely attributable to mitonuclear discordance, which can arise from biological factors (e.g., introgressive hybridization and/or incomplete lineage sorting). Our analyses indicate that past introgressions have taken place, but that there is no evidence to suggest an ongoing admixture between the two species, demonstrating that the taxa have reached full postzygotic reproductive isolation and hence represent two distinct biological species. Finally, we discuss our results from an evolutionary point of view taking into consideration the past climatic oscillations that have likely shaped the present dynamics between the two species. Overall, our study demonstrates the effectiveness of the target enrichment approach in resolving shallow phylogenetic relationships under complex evolutionary circumstances and that this approach is useful in establishing robust and well-informed taxonomic delimitations involving parapatric taxa.

Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub

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.

Phylogeography, hybrid zones and contemporary species boundaries in the south-eastern Australian smooth frogs (Anura: Myobatrachidae: Geocrinia)

Paleo-climatic fluctuations have driven episodic changes in species distributions, providing opportunities for populations to diverge in isolation and hybridise following secondary contact. Studies of phylogeographic diversity and patterns of gene flow across hybrid zones can provide insight into contemporary species boundaries and help to inform taxonomic and conservation inferences. Here we explore geographic diversity within the acoustically divergent yet morphologically conserved south-eastern Australian smooth frog complex and assess gene flow across a narrow hybrid zone using mitochondrial nucleotide sequences and nuclear genome-wide single nucleotide polymorphisms. Our analyses reveal the presence of an evolutionarily distinct taxon restricted to the Otway Plains and Ranges, Victoria, which forms a narrow (9-30 km wide), spatiotemporally stable (>50 years) hybrid zone with Geocrinia laevis, which we describe herein as a new species.

Conservation Priorities and Demographic History of Saussurea involucrata in the Tianshan Mountains and Altai Mountains

Rare and vulnerable endemic plants represent different evolutionary units that occur at different times, and protecting these species is a key issue in biological protection. Understanding the impact of the history of endangered plant populations on their genetic diversity helps to reveal evolutionary history and is crucial for guiding conservation efforts. Saussurea involucrata, a perennial alpine species mainly distributed in the Tianshan Mountains, is famous for its medicinal value but has become endangered due to over-exploitation. In the present study, we employed both nuclear and chloroplast DNA sequences to investigate the genetic distribution pattern and evolutionary history of S. involucrata. A total of 270 individuals covering nine S. involucrata populations were sampled for the amplification and sequencing of nrDNA Internal Transcribed Spacer (ITS) and chloroplast trnL-trnF, matK and ndhF-rpl32 sequences. Via calculation, we identified 7 nuclear and 12 plastid haplotypes. Among the nine populations, GL and BA were characterized by high haplotype diversity, whereas BG revealed the lowest haplotype diversity. Molecular dating estimations suggest that divergence among S. involucrata populations occurred around 0.75 Ma, coinciding with the uplift of Tianshan Mountains. Our results reveal that both isolation-by-distance (IBD) and isolation-by-resistance (IBR) have promoted genetic differentiation among populations of S. involucrata. The results from the ecological niche modeling analyses show a more suitable habitat for S. involucrata in the past than at present, indicating a historical distribution contraction of the species. This study provides new insight into understanding the genetic differentiation of S. involucrata, as well as the theoretical basis for conserving this species.

Investigation of phylogenetic relationships within Saxifraga diversifolia complex (Saxifragaceae) based on restriction-site associated DNA sequence markers

Subsect. Hirculoideae Engl. & Irmsch., belonging to Saxifraga sect. Ciliatae Haw., has high species richness. It can be divided into S. diversifolia, S. pseudohirculus, and S. sinomontana complexes based on morphological characteristics. The species with prominent leaf veins on the posterior leaf edge were placed in the S. diversifolia complex, which is mainly distributed on the eastern and southern margins of the Qinghai-Tibetan Plateau. In this study, 53 samples, representing 15 of the 33 described species in the S. diversifolia complex, were sequenced using the Restriction-site Associated DNA Sequence (RAD-seq) technique. A total of 111,938 high-quality SNP loci were screened to investigate the phylogenetic relationships within the S. diversifolia complex. The result of the neighbor-joining (NJ) tree shows that the S. diversifolia complex is a paraphyletic group. Despite of some inconsistencies as revealed by genetic structural analysis, clustering results of representative species reconstructed by both NJ and principal component analysis analyses support previous biogeographic and morphological evidences. In addition, long-distance gene flow events for 11 taxa were detected in the S. diversifolia complex, respectively from S. implicans 1 to S. implicans 2, S. diversifolia and S. maxionggouensis, and from S. maxionggouensis to S. nigroglandulifera. These findings may improve our comprehension of the phylogeny, classification, and evolution of the S. diversifolia complex.

Dispersal capacities of pollen, seeds and spores: insights from comparative analyses of spatial genetic structures in bryophytes and spermatophytes

Introduction

The dramatic fluctuations of climate conditions since the late Tertiary era have resulted in major species range shifts. These movements were conditioned by geographic barriers and species dispersal capacities. In land plants, gene flow occurs through the movement of male gametes (sperm cells, pollen grains), which carry nDNA, and diaspores (spores, seeds), which carry both cpDNA and nDNA, making them an ideal model to compare the imprints of past climate change on the spatial genetic structures of different genomic compartments. Based on a meta-analysis of cpDNA and nDNA sequence data in western Europe, we test the hypotheses that nDNA genetic structures are similar in bryophytes and spermatophytes due to the similar size of spores and pollen grains, whereas genetic structures derived from the analysis of cpDNA are significantly stronger in spermatophytes than in bryophytes due to the substantially larger size of seeds as compared to spores.

Methods

Sequence data at 1-4 loci were retrieved for 11 bryophyte and 17 spermatophyte species across their entire European range. Genetic structures between and within southern and northern populations were analyzed through F and N statistics and Mantel tests.

Results and discussion

Gst and Nst between southern and northern Europe derived from cpDNA were significantly higher, and the proportion of significant tests was higher in spermatophytes than in bryophytes. This suggests that in the latter, migrations across mountain ranges were sufficient to maintain a homogenous allelic structure across Europe, evidencing the minor role played by mountain ranges in bryophyte migrations. With nDNA, patterns of genetic structure did not significantly differ between bryophytes and spermatophytes, in line with the hypothesis that spores and pollen grains exhibit similar dispersal capacities due to their size similarity. Stronger levels of genetic differentiation between southern and northern Europe, and within southern Europe, in spermatophytes than in bryophytes, caused by higher long-distance dispersal capacities of spores as compared to seeds, may account for the strikingly higher levels of endemism in spermatophytes than in bryophytes in the Mediterranean biodiversity hotspot.

Phylogeography of Drosophila buzzatii (Diptera, Drosophilidae): responses of the species to Quaternary climates in tropical and subtropical South America

Drosophila buzzatii (Diptera: Drosophilidae) is a fly that breeds exclusively on decaying tissues of cacti species widely distributed in tropical and subtropical areas of South America. This distribution includes biomes in distinct climatic regimes (e.g., seasonal rain forest, semi-arid scrubs, savannas, and grasslands), which at first glance could might give the false impression that the species is not sensitive to either climate or vegetation physiognomies. However, detection of historical demographic events within D. buzzatii reveal the interplay between climate and the population structure of the species as the Late Quaternary climate changes occurred. To understand this process, we performed a phylogeographic analysis based on sequences of the mitochondrial gene COI for 128 individuals from 43 localities. Our analyses combined coalescent methods, population genetics, and paleodistributions estimation methods. Our study reveals that the COI haplotype diversity is geographically structured, with a decreasing cline from north to south. The results suggest an ancient range expansion, dated from 610k to 550k years before present, in the northernmost region of the species distribution, the Caatinga vegetation. More recently, an intense gene flow and a population expansion were detected in the central and south portions of its distribution. The demographic events detected date back to the glacial periods of the Quaternary.

Molecular phylogeography of Hipposideros pratti in China

Hipposideros pratti has low genetic diversity and was divided into two clades, the central-western clade and the eastern clade. We did not detect a clear east-to-west dispersal route along the longitudinal direction, and we found that the eastern clade spread outward from one population to another, while the central-western clade spread gradually. The glacial-interglacial period of the Quaternary influenced the migration and dispersal of H. pratti. H. pratti did not experience a significant population increase in the past, and the average population trajectory was decreasing. Given the convenient ecosystem services provided by bats, the preservation of bat populations is particularly critical. Nonetheless, we have discovered that the majority of H. pratti's distributional sites were not discovered in this study. Based on our results, it is important to apply in situ conservation measures for effective protection as soon as possible.

Distribution patterns of Quercus ilex from the last interglacial period to the future by ecological niche modeling

The plants' geographic distribution is affected by natural or human-induced climate change. Numerous studies at both the global and regional levels currently focus on the potential changes in plant distribution areas. Ecological niche modeling can help predict the likely distribution of species according to environmental variables under different climate scenarios. In this study, we predicted the potential geographic distributions of Quercus ilex L. (holm oak), a keystone species of the Mediterranean ecosystem, for the Last Interglacial period (LIG: ~130 Ka), the Last Glacial Maximum (LGM: ~22 Ka), mid-Holocene (MH: ~6 Ka), and future climate scenarios (Representative Concentration Pathway (RCP) 4.5 and 8.5 scenarios) for 2050-2070 obtained from CCSM4 and MIROC-ESM global climate scenarios respectively. The models were produced with algorithms from the R-package "biomod2" and assessed by AUC of the receiver operating characteristic plot and true skill statistics. Aside from BIOCLIM (SRE), all model algorithms performed similarly and produced projections that are supported by good evaluation scores, although random forest (RF) slightly outperformed all the others. Additionally, distribution maps generated for the past period were validated through a comparison with pollen data acquired from the Neotoma Pollen Database. The results revealed that southern areas of the Mediterranean Basin, particularly coastal regions, served as long-term refugia for Q. ilex, which was supported by fossil pollen data. Furthermore, the models suggest long-term refugia role for Anatolia and we argue that Anatolia may have served as a founding population for the species. Future climate scenarios indicated that Q. ilex distribution varied by region, with some areas experiencing range contractions and others range expands. This study provides significant insights into the vulnerability of the Q. ilex to future climate change in the Mediterranean ecosystem and highlights the crucial role of Anatolia in the species' historical distribution.

Can demographic histories explain long-term isolation and recent pulses of asymmetric gene flow between highly divergent grey fox lineages?

Secondary contact zones between deeply divergent, yet interfertile, lineages provide windows into the speciation process. North American grey foxes (Urocyon cinereoargenteus) are divided into western and eastern lineages that diverged approximately 1 million years ago. These ancient lineages currently hybridize in a relatively narrow zone of contact in the southern Great Plains, a pattern more commonly observed in smaller-bodied taxa, which suggests relatively recent contact after a long period of allopatry. Based on local ancestry inference with whole-genome sequencing (n = 43), we identified two distinct Holocene pulses of admixture. The older pulse (500-3500 YBP) reflected unidirectional gene flow from east to west, whereas the more recent pulse (70-200 YBP) of admixture was bi-directional. Augmented with genotyping-by-sequencing data from 216 additional foxes, demographic analyses indicated that the eastern lineage declined precipitously after divergence, remaining small throughout most of the late Pleistocene, and expanding only during the Holocene. Genetic diversity in the eastern lineage was highest in the southeast and lowest near the contact zone, consistent with a westward expansion. Concordantly, distribution modelling indicated that during their isolation, the most suitable habitat occurred far east of today's contact zone or west of the Great Plains. Thus, long-term isolation was likely caused by the small, distant location of the eastern refugium, with recent contact reflecting a large increase in suitable habitat and corresponding demographic expansion from the eastern refugium. Ultimately, long-term isolation in grey foxes may reflect their specialized bio-climatic niche. This system presents an opportunity for future investigation of potential pre- and post-zygotic isolating mechanisms.