Lineage diversification and historical demography of a montane bird Garrulax elliotii--implications for the Pleistocene evolutionary history of the eastern Himalayas

Genetic insights for enhancing conservation strategies in captive and wild Asian elephants through improved non-invasive DNA-based individual identification

Asian elephant is a key umbrella species that plays a crucial role in maintaining biodiversity and ecological balance. As an iconic symbol of Thailand, it also significantly contributes to the nation's tourism industry. However, human activities pose serious threats to their long-term survival and population health. To tackle these challenges and develop effective conservation strategies, extensive genetic reference data were collected to enhance both captive and wild elephant conservation, improve non-invasive DNA-based individual identification, and assess genetic diversity using 18 microsatellite markers. High genetic diversity was observed across all populations; however, high levels of inbreeding were evident in NEI, EKS, BCEP, and wild elephant populations, except for the MEP population, which recorded low inbreeding levels. Significant variation in the gene pool estimates was observed across different populations, with three maternal haplogroups (α, β1, and a tentative β3) identified. A reduced panel of six microsatellite markers proved highly efficient for individual identification. Additionally, non-invasive DNA samples were tested using 18 microsatellite loci for individual identification. Using only 7 out of the 18 microsatellite loci tested, individuals were successfully identified, demonstrating enough discriminatory power to distinguish between individuals. Among these, four loci (LaT08, LaT13, FH19, and FH67) were both effective and efficient for reliable individual identification in fecal samples. These findings offer valuable insights for optimizing conservation efforts, including the design of tailored strategies to protect Asian elephants in Thailand and ensure the long-term viability of their populations.

Late Pleistocene Altitudinal Segregation and Demography Define Future Climate Change Distribution of the Peromyscus mexicanus Species Group: Conservation Implications

Mountains harbor a significant number of the World's biodiversity, both on tropical and temperate regions. Notably, one crucial gap in conservation is the consideration of historical and contemporary patterns influencing differential distribution in small mammal mountain species and how climate change will affect their distribution and survival. The mice Peromyscus mexicanus species group is distributed across mountains in Guatemala-Chiapas and Central America, which experienced significant effects of glacial and interglacial cycles. We determined phylogeographic and demographic patterns of lowlands and highlands mountain lineages, revealing that the radiation of modern P. mexicanus lineages occurred during the Pleistocene (ca. 2.6 mya) along Nuclear Central America. In concert with climatic cycles and the distribution of habitats, lowland and highland lineages showed recent population size increase and decrease, respectively. We also estimated the current and future distribution ranges for six lineages, finding marked area size increase for two lineages for which vegetation type and distribution would facilitate migrating towards higher elevations. Contrastingly, three lineages showed range size decrease; their ecological requirements make them highly susceptible to future habitat loss. Our findings are clear evidence of the negative impacts of future climate change, while our ability to manage and conserve these vulnerable ecosystems and mountain species is contingent on our understanding of the implications of climate change on the distribution, ecology, and genetics of wildlife populations.

Multiple concordant cytonuclear divergences and potential hybrid speciation within a species complex in Asia

Various environmental factors impact the distribution, population structure, demography and evolutionary trajectory of a bird species, leading to genetic and morphological divergences between populations across its distribution. The Paradoxornis webbianus species complex is found throughout much of East Asia, where its geographically distinct populations exhibit dramatic morphological variation. This has resulted in a hotly debated taxonomy. This study intended to identify genetic divergence patterns and their underlying contributing factors for this species complex. We collected 243 birds, whose data was combined with those available in GenBank to perform phylogeographic analyses using one mitochondrial and six nuclear loci. Six mitochondrial clades were observed in the species complex, while individual-based Bayesian clustering using nuclear markers showed multiple congruent breaks. Overall, the six molecular lineages could be recognized as independent species under the lineage species concept in view of genetic divergence, clade-specific morphological changes and distribution: P. webbianus, P. w. bulomachus, P. alphonsianus, P. a. ganluoensis, P. brunneus brunneus and P. b. ricketti. The estimated divergence times range from 0.46 to 3.36 million years ago, suggesting it was likely impacted by paleoclimatic changes. Interestingly, P. alphonsianus carries two divergent mitochondrial lineages shared with P. webbianus and P. a. ganluoensis, respectively, and analyses based on nuclear loci found a similar pattern. We discussed the various hypotheses for this pattern and argued that P. alphonsianus was likely the result of hybridization between P. webbianus and P. a. ganluoensis. Further data on genome, transcriptome and breeding ecology are needed to address the hypothesis of hybrid speciation and its underlying mechanisms.

The combination of genomic offset and niche modelling provides insights into climate change-driven vulnerability

Global warming is increasingly exacerbating biodiversity loss. Populations locally adapted to spatially heterogeneous environments may respond differentially to climate change, but this intraspecific variation has only recently been considered when modelling vulnerability under climate change. Here, we incorporate intraspecific variation in genomic offset and ecological niche modelling to estimate climate change-driven vulnerability in two bird species in the Sino-Himalayan Mountains. We found that the cold-tolerant populations show higher genomic offset but risk less challenge for niche suitability decline under future climate than the warm-tolerant populations. Based on a genome-niche index estimated by combining genomic offset and niche suitability change, we identified the populations with the least genome-niche interruption as potential donors for evolutionary rescue, i.e., the populations tolerant to climate change. We evaluated potential rescue routes via a landscape genetic analysis. Overall, we demonstrate that the integration of genomic offset, niche suitability modelling, and landscape connectivity can improve climate change-driven vulnerability assessments and facilitate effective conservation management.

A Phylogeographical Analysis of the Beetle Pest Species Callosobruchus chinensis (Linnaeus, 1758) in China

Callosobruchus chinensis (Coleoptera Bruchidae), is a pest of different varieties of legumes. In this paper, a phylogeographical analysis of C. chinensis was conducted to provide knowledge for the prevention and control of C. chinensis. A total of 224 concatenated mitochondrial sequences were obtained from 273 individuals. Suitable habitat shifts were predicted by the distribution modelling (SDM). Phylogeny, genetic structure and population demographic history were analyzed using multiple software. Finally, the least-cost path (LCP) method was used to identify possible dispersal corridors and genetic connectivity. The SDM results suggested that the distribution of C. chinensis experienced expansion and contraction with changing climate. Spatial distribution of mtDNA haplotypes showed there was partial continuity among different geographical populations of C. chinensis, except for the Hohhot (Inner Mongolia) population. Bayesian skyline plots showed that the population had a recent expansion during 0.0125 Ma and 0.025 Ma. The expansion and divergent events were traced back to Quaternary glaciations. The LCP method confirmed that there were no clear dispersal routes. Our findings indicated that climatic cycles of the Pleistocene glaciations, unsuitable climate and geographic isolation played important roles in the genetic differentiation of C. chinensis. Human activities weaken the genetic differentiation between populations. With the change in climate, the suitable areas of C. chinensis will disperse greatly in the future.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

Ring distribution patterns-diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin

Studying the genetic differentiation in a unique geographical area contributes to understanding the process of speciation. Here, we explore the spatial genetic structure and underlying formation mechanism of two congeneric small mammal species (Apodemus draco and A. chevrieri), which are mainly distributed in the mountains surrounding the lowland Sichuan Basin, southwest China. We applied a set of comparative phylogeographical analyses to determine their genetic diversification patterns, combining mitochondrial (Cytb and COI) and nuclear (microsatellite loci) markers, with dense sampling throughout the range (411 A. draco from 21 sites and 191 A. chevrieri from 22 sites). Moreover, we performed three complementary statistical methods to investigate the correlation between genotype and geographical and environmental components, and predicted the potential suitable distributional range under the present and historical climate conditions. Our results suggest that both species have experienced allopatric differentiation and admixture in historical periods, resulting in a ring-shape diversification, under the barrier effect of the Sichuan Basin. We infer that the tectonic events of the Qinghai-Tibetan Plateau and climatic oscillations during the Quaternary played an important role on the genetic divergence of the two species by providing environmental heterogeneity and geographical variation. Our study reveals a case of two sympatric small mammals following a ring-shaped diversification pattern and provides insight into the process of differentiation.

An Orchid in Retrograde: Climate-Driven Range Shift Patterns of Ophrys helenae in Greece

Climate change is regarded as one of the most important threats to plants. Already species around the globe are showing considerable latitudinal and altitudinal shifts. Helen's bee orchid (Ophrys helenae), a Balkan endemic with a distribution center in northwestern Greece, is reported to be expanding east and southwards. Since this southeastern movement goes against the usual expectations, we investigated via Species Distribution Modelling, whether this pattern is consistent with projections based on the species' response to climate change. We predicted the species' future distribution based on three different climate models in two climate scenarios. We also explored the species' potential distribution during the Last Interglacial and the Last Glacial Maximum. O. helenae is projected to shift mainly southeast and experience considerable area changes. The species is expected to become extinct in the core of its current distribution, but to establish a strong presence in the mid- and high-altitude areas of the Central Peloponnese, a region that could have provided shelter in previous climatic extremes.

Isolated alpine habitats reveal disparate ecological drivers of taxonomic and functional beta-diversity of small mammal assemblages

The interpretation of patterns of biodiversity requires the disentanglement of geographical and environmental variables. Disjunct alpine communities are geographically isolated from one another but experience similar environmental impacts. Isolated homogenous habitats may promote speciation but constrain functional trait variation. In this study, we examined the hypothesis that dispersal limitation promotes taxonomic divergence, whereas habitat similarity in alpine mountains leads to functional convergence. We performed standardized field investigation to sample non-volant small mammals from 18 prominent alpine sites in the Three Parallel Rivers area. We estimated indices quantifying taxonomic and functional alpha- and beta-diversity, as well as beta-diversity components. We then assessed the respective importance of geographical and environmental predictors in explaining taxonomic and functional compositions. No evidence was found to show that species were more functionally similar than expected in local assemblages. However, the taxonomic turnover components were higher than functional ones (0.471±0.230 vs. 0.243±0.215), with nestedness components showing the opposite pattern (0.063±0.054 vs. 0.269±0.225). This indicated that differences in taxonomic compositions between sites occurred from replacement of functionally similar species. Geographical barriers were the key factor influencing both taxonomic total dissimilarity and turnover components, whereas functional beta-diversity was primarily explained by climatic factors such as minimum temperature of the coldest month. Our findings provide empirical evidence that taxonomic and functional diversity patterns can be independently driven by different ecological processes. Our results point to the importance of clarifying different components of beta-diversity to understand the underlying mechanisms of community assembly. These results also shed light on the assembly rules and ecological processes of terrestrial mammal communities in extreme environments.

The genetic differentiation of a cricket (Velarifictorus micado) with two modes of life cycle in East Asia after the middle Pleistocene and the invasion origin of the United States of America

The cricket Velarifictorus micado is widely distributed in East Asia and colonized the United States of America (the USA) in 1959. It has two life cycles: egg and nymph diapause. We aimed to investigate the biogeographic boundary between them and determine when and why V. micado diverged. Mitochondrial fragments including COI and CytB were used for haplotype network, demographic analysis, and divergence time estimation in individuals of East Asia. We selected several samples from the USA to find out the colonization origin. The haplotype network indicated there were three lineages based on COI, NE lineage (the egg diapause and mainly distributed in the northern regions), SE lineage (the egg diapause and mainly distributed in the southern regions), and SN lineage (the nymph diapause and mainly distributed in the southern regions). The molecular chronograms indicated that the first divergence of V. micado into two main lineages, NE and southern lineages (SE and SN), was essentially bounded by the Yangtze River. It occurred around ~0.79 Ma (95% HPD: 1.13-0.46 Ma) in the Middle Pleistocene Transition. This was followed by the divergence of the southern lineage into two sublineages, SE and SN lineage, occurred around ~0.50 Ma (95% HPD: 0.71-0.25 Ma), corresponding to the time of development of glaciers in various parts of the Qinghai-Tibet Plateau (QTP) (0.73-0.46 Ma). SE lineage might originate from southwestern China based on the comparison between the haplotype network based on COI and CytB. Our study suggested that divergences of lineages have twice co-occurred with tendency of cooling climatic in Asia after the Mid-Pleistocene, and the life-history strategy may play an important role in lineage diversification. Additionally, our results indicated that the USA populations were revealed at least twice separate Asian invasions. These both belonged to the egg diapause, which might provide a new perspective for invasion control.

Last Glacial Maximum led to community-wide population expansion in a montane songbird radiation in highland Papua New Guinea

BACKGROUND

Quaternary climate fluctuations are an engine of biotic diversification. Global cooling cycles, such as the Last Glacial Maximum (LGM), are known to have fragmented the ranges of higher-latitude fauna and flora into smaller refugia, dramatically reducing species ranges. However, relatively less is known about the effects of cooling cycles on tropical biota.

RESULTS

We analyzed thousands of genome-wide DNA markers across an assemblage of three closely related understorey-inhabiting scrubwrens (Sericornis and Aethomyias; Aves) from montane forest along an elevational gradient on Mt. Wilhelm, the highest mountain of Papua New Guinea. Despite species-specific differences in elevational preference, we found limited differentiation within each scrubwren species, but detected a strong genomic signature of simultaneous population expansions at 27-29 ka, coinciding with the onset of the LGM.

CONCLUSION

The remarkable synchronous timing of population expansions of all three species demonstrates the importance of global cooling cycles in expanding highland habitat. Global cooling cycles have likely had strongly different impacts on tropical montane areas versus boreal and temperate latitudes, leading to population expansions in the former and serious fragmentation in the latter.

Shared response to changes in drainage basin: Phylogeography of the Yunnan small narrow-mouthed frog, Glyphoglossus yunnanensis (Anura: Microhylidae)

AIM

With the late Cenozoic uplift of the Qinghai-Tibetan Plateau (QTP), drainage of the southeastern edge of the QTP changed significantly. However, the impact of this dramatic change on the geographical distribution and genetic diversity of endemic organisms is still poorly understood. Here, we examined the geographical patterns of genetic variation in the Yunnan small narrow-mouthed frog, Glyphoglossus yunnanensis (Microhylidae), and two alternative hypotheses were tested: That is, the geographical distribution of genetic variation was determined by either the contemporary drainage basin or historical drainage basins.

LOCATION

The Mountains of southwest China.

MATERIALS AND METHODS

Analyses were based on 417 specimens collected from across the distribution of the species. We reconstructed the genealogy (Bayesian and maximum parsimony methods) and assessed demographic history based on DNA sequencing data from mitochondrial and nuclear markers. We also mapped the genetic diversity and estimated the divergence times by a relaxed clock model.

RESULTS

The species has maintained a relatively stable population size without recent population expansion. Four major maternal lineages were identified with good support, one representing a possible cryptic species and the other three showing further subdivision. The distribution of these deeply differentiated lineages/sublineages corresponded well to geographical regions. The secondary contact zones and phylogeographic breaks in distinct lineages of G. yunnanensis were almost concordant with those of Nanorana yunnanensis.

MAIN CONCLUSIONS

Lineage division conformed to the hypothesis of drainage system evolution, that is, the phylogeographic pattern of G. yunnanensis was shaped by historical drainage patterns. Concordance in phylogeographic patterns may suggest a shared response to common hydrogeological history and also might indicate that there was more contribution of the drainage history than ecological or life-history traits in structuring genetic variation between these two disparate codistributed taxa G. yunnanensis and N. yunnanensis.

Sky islands as foci for divergence of fig trees and their pollinators in southwest China

The dynamics of populations and their divergence over time have shaped current levels of biodiversity and in the case of the "sky islands" of mountainous southwest (SW) China have resulted in an area of exceptional botanical diversity. Ficus tikoua is a prostrate fig tree subendemic to the area that displays unique intraspecific diversity, producing figs typical of different pollination modes in different parts of its range. By combining climate models, genetic variation in populations of the tree's obligate fig wasp pollinators and distributions of the different plant phenotypes, we examined how this unusual situation may have developed. We identified three genetically distinct groups of a single Ceratosolen pollinator species that have largely parapatric distributions. The complex topography of the region contributed to genetic divergence among the pollinators by facilitating geographical isolation and providing refugia. Migration along elevations in response to climate oscillations further enhanced genetic differentiation of the three pollinator groups. Their distributions loosely correspond to the distributions of the functionally significant morphological differences in the male figs of their host plants, but postglacial expansion of one group has not been matched by spread of its associated plant phenotype, possibly due to a major river barrier. The results highlight how interplay between the complex topography of the "sky island" complex and climate change has shaped intraspecies differentiation and relationships between the plant and its pollinator. Similar processes may explain the exceptional botanical diversity of SW China.

Evaluating Population Genetic Structure and Demographic History of Quercus spinosa (Fagaceae) Based on Specific Length Amplified Fragment Sequencing

Effectively identifying the genetic structure and related factors of a species can facilitate understanding the evolutionary history of the species. Phylogeographic patterns and genetic data are essential in investigating the species historical processes and diversification that response to environmental, climatic and geological influences. In this study, Specific Length Amplified Fragment Sequencing (SLAF-seq) data and ecological niche models (ENMs) are combined to identify the genetic structure and demographic modeling of Quercus spinosa, and evaluate the impacts of historical range shifts, climatic variation, and landscape factors on this species. The population topology and genetic divergence of the Cenozoic were inferred by a site frequency spectrum based composite-likelihood approach which is a novel strategy for maximizing the utility of linked SLAF markers. The overall genetic structure using model-based and model-free clustering methods was consistently identified as two geographically distinct genetic clusters. A deep divergence between two natural lineages (i.e., a western Himalaya-Hengduan Mountains lineage and an eastern Qin-ling Mountains lineage) was observed. The demographic modeling and Niche reconstruction indicated that the two groups were diverged in the late Miocene and then presented as two distinct genetic lineages. With the Quaternary glacial climate fluctuation, two groups had continuous asymmetrical secondary contact and gene exchange in the Sichuan Basin during the last glacial maximum. Besides, a significant relationship between genetic distance and geography in all individuals was identified by the Mantel test. Overall, this study 1) contributes to a better understanding of the role played by Quaternary climatic fluctuation in the present-day distributions of Q. spinosa; 2) provides a comprehensive view of the genome-wide variation of sclerophyllous forests in ecological adaptive evolution; 3) indicates that dispersal limitation and ecological divergence contribute to the genome-wide differentiation of Q. spinosa, which supports a hypothesis that complex geography and climatic changes strongly influence the evolutionary origin and history of the species.

Molecular phylogeographic analyses and species delimitations reveal that Leopoldamys edwardsi (Rodentia: Muridae) is a species complex

Leopoldamys edwardsi is a species with wide distribution ranges in southern China but is not discussed in studies on geographic variation and species differentiation. We used 2 mitochondrial (Cytb, CO1) and 3 nuclear (GHR, IRBP and RAG1) genes to clarify species phylogeography and geographical differentiation. Maximum likelihood (ML) and Bayesian phylogenetic inference (BI) trees consistently indicated that L. edwardsi is a species complex containing 3 main lineages with high Kimura-2-parameter (K2P) divergences (i.e. lineages L_N , L_S and L_HN ) found in the northern and southern China and Hainan Island, respectively. The 3 species delimitation methods, automated barcoding gap discovery, Bayesian poisson tree process analysis and Bayesian phylogenetics and phylogeography, consistently supported the existence of cryptic species. Divergence times among the main lineages were inferred to be during the Pleistocene, with L_HN /L_S split at 1.33 Ma and L_N /(L_HN +L_S ) at 2.61 Ma; the diversifications of L. edwardsi complex might be caused by the rapid uplifts of Tibetan Plateau, paleoclimate change and complex topography. The divergence between L_HN and L_S was probably related to the separation of Hainan Island from the mainland via the formation of the Qiongzhou Strait. Lineages L_N and (L_S +L_HN ) likely diverged due to the Wuyi-Nanling mountain range forming a dispersal barrier. Our results suggested that L. edwardsi complex contains at least 3 distinct species: L_HN represents L. hainanensis, endemic to Hainan Island and previously considered as a subspecies L. e. hainanensis; L_S represents a cryptic species distributed throughout the southern Chinese continent; and L_N represents the nominotypical species L. edwardsi.

Evolutionary history of a beautiful damselfly, Matrona basilaris, revealed by phylogeographic analyses: the first study of an odonate species in mainland China

Matrona basilaris Selys, 1853 is a damselfly distributed mainly in mainland China. A total of 423 individuals from 48 populations covering almost the entire range were sampled to explore the genetic diversity, phylogeographic structure, and demographic dynamics of the species using sequences of three mitochondrial genes (COI, COII, and ND1) and a nuclear (ITS1 + 5.8 S + ITS2) gene. Phylogenetic tree, median-joining network, and BAPS analyses indicated a four-group division of the entire population, and the divergence event was estimated to have occurred in the middle Pleistocene. The diverse terrain of mainland China as well as past climatic oscillations were assumed to have shaped the current phylogeographic pattern of M. basilaris. Multiple lines of evidence supported population expansion in Group 1 and Group 2 but not in Group 3 or Group 4. The expansion times corresponded to the transition phase from the LIG (∼0.14-0.12 Mya) to the LGM (∼0.021-0.018 Mya). The pre-LGM expansion model reflected a different pattern affecting the historical dynamics of the population of East Asian species caused by Pleistocene climatic changes. Interestingly, Group 2 exhibited a disjunctive distribution pattern. The possible reasons were introgression caused by female-biased dispersal or human phoresy during construction of the Forbidden City during the Ming Dynasty of China.

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon-primed intron-crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate-driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Phylogeographic structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China

BACKGROUND

Species with broad distributions frequently divide into multiple genetic forms and may therefore be viewed as "cryptic species". Here, we used the mitochondrial cytochrome b (Cytb) and 12 nuclear DNA loci to investigate phylogeographic structures of the sharpbelly (Hemiculter leucisculus) in rivers in southern China and explored how the geological and climatic factors have shaped the genetic diversity and evolutionary history of this species.

RESULTS

Our mitochondrial phylogenetic analysis identified three major lineages (lineages A, B, and C). Lineages B and C showed a relatively narrower geographic distribution, whereas lineage A was widely distributed in numerous drainages. Divergence dates suggested that H. leucisculus populations diverged between 1.61-2.38 Ma. Bayesian species delimitation analysis using 12 nuclear DNA loci indicated the three lineages probably represented three valid taxa. Isolation-with-migration (IM) analysis found substantial gene flow has occurred among the three lineages. Demographic analyses showed that lineages B and C have experienced rapid demographic expansion at 0.03 Ma and 0.08 Ma, respectively.

CONCLUSIONS

Hemiculter leucisculus populations in drainages in southern China comprise three mtDNA lineages, and each of which may represent a separate species. Intense uplift of the Qinghai-Tibetan Plateau, evolution of Asian monsoons, changes in paleo-drainages, and poor dispersal ability may have driven the divergence of the three putative species. However, gene flow occurs among the three lineages. Climatic fluctuations have a prominent impact on the populations from the lineages B and C, but exerted little influence on the lineage A.

Phylogeography of the Rice Spittle Bug (Callitettix versicolor) Implies Two Long-Term Mountain Barriers in South China

South China is a region of remarkable topographic complexity. However, the impact of climate fluctuations in the Pleistocene on the local fauna and especially insects has not been extensively studied. We integrated mitochondrial DNA (mtDNA) and microsatellite data of the rice spittle bug, Callitettix versicolor, to determine the genetic structure, potential biogeographic barriers, and historical demography of this species. The mtDNA data revealed two distinct lineages (Western and Eastern) congruent with the geographically separated western and eastern sub-regions of the Hengduan Mountains. The Eastern lineage was subdivided into two sub-lineages, E1 and E2, congruent with the geographically separated northern and southern sub-regions of the Dabie Mountains. E2 was further subdivided into two sub-groups, E2-1 and E2-2, with a hybrid zone (Guizhou and Hunan Provinces) in which their areas were contiguous. The genetic structures constructed using mtDNA were corroborated by four clusters (G1-G4) of microsatellite data. The populations of each cluster were nearly consistent with a sub-lineage of the mtDNA gene tree (G1-G4 corresponded to the Western, E1, E2-2 and E2-1 lineages, respectively). The divergence time estimated between the Western and Eastern lineages was 1.17 (0.50-2.37) to 0.89 (0.39-1.78) Mya, indicating that the lineages diversified on both geographic and temporal scales. The historical demography of the Eastern lineage showed continuous population growth after the Last Interglacial (LIG) and a stable population during the Last Glacial Maximum (LGM) period. However, the Western lineage remained largely unchanged during the LIG and LGM periods. This suggests that the historical demography of C. versicolor is probably related not only to the paleoclimate of South China, but also to the geological restriction and specific habitat preferences of species.

Investigating the Genetic Diversity, Population Differentiation and Population Dynamics of Cycas segmentifida (Cycadaceae) Endemic to Southwest China by Multiple Molecular Markers

Climate change, species dispersal ability and habitat fragmentation are major factors influencing species distribution and genetic diversity, especially for the range-restricted and threatened taxa. Here, using four sequences of chloroplast DNAs (cpDNAs), three nuclear genes (nDNAs) and 12 nuclear microsatellites (SSRs), we investigated the genetic diversity, genetic structure, divergence time and population dynamics of Cycas segmentifida D. Y. Wang and C. Y. Deng, a threatened cycad species endemic to Southwest China. High levels of genetic diversity and genetic differentiation were revealed in C. segmentifida. Haplotypes of networks showed two evolutionary units in C. segmentifida, with the exception of the nuclear gene GTP network. Meanwhile, the UPGMA tree, structure and PCoA analyses suggested that 14 populations of C. segmentifida were divided into two clades. There was significant effect of isolation by distance (IBD) in this species. However, this species did not display a significant phylogeographic structure. The divergence time estimation suggested that its haplotypes diverged during the Middle Pleistocene. Additionally, the population dynamics inferred from different DNA sequences analyses were discordant. Bottleneck analysis showed that populations of C. segmentifida did not experience any recent bottleneck effect, but rather pointed to a contraction of its effective population size over time. Furthermore, our results suggested that the population BM which held an intact population structure and occupied undisturbed habitat was at the Hardy-Weinberg equilibrium, implying that this population is a free-mating system. These genetic features provide important information for the sustainable management of C. segmentifida.

Middle-Upper Pleistocene climate changes shaped the divergence and demography of Cycas guizhouensis (Cycadaceae): Evidence from DNA sequences and microsatellite markers

Climatic oscillations in the Pleistocene have had profound effects on the demography and genetic diversity of many extant species. Cycas guizhouensis Lan & R.F. Zou is an endemic and endangered species in Southwest China that is primarily distributed along the valleys of the Nanpan River. In this study, we used four chloroplast DNAs (cpDNA), three nuclear genes (nDNA) and 13 microsatellite (SSR) loci to investigate the genetic structure, divergence time and demographic history of 11 populations of C. guizhouensis. High genetic diversity and high levels of genetic differentiation among the populations were observed. Two evolutionary units were revealed based on network and Structure analysis. The divergence time estimations suggested that haplotypes of C. guizhouensis were diverged during the Middle-Upper Pleistocene. Additionally, the demographic histories deduced from different DNA sequences were discordant, but overall indicated that C. guizhouensis had experienced a recent population expansion during the post-glacial period. Microsatellite data revealed that there was a contraction in effective population size in the past. These genetic features allow conservation measures to be taken to ensure the protection of this endangered species from extinction.

Phylogeography of endemic Xantus' hummingbird (Hylocharis xantusii) shows a different history of vicariance in the Baja California Peninsula

Studies of phylogeographic patterns provide insight into the processes driving lineage divergence in a particular region. To identify the processes that caused phylogeographic breaks, it is necessary to use historical information and a set of appropriate molecular data to explain current patterns. To understand the influence of geological or ecological processes on the phylogeography of the only species of hummingbird endemic to the Baja California Peninsula, Hylocharis xantusii, mitochondrial DNA sequences of three concatenated genes (Cyt-b, COI and ND2; 2297bp in total) in 100 individuals were analyzed. The spatial analyses of genetic variation showed phylogeographic structure consisting of a north, central and south regions. According to estimated divergence times, two vicariant events are supported, permanent separation of the peninsula and formation of the Gulf of California at 5mya and temporary isolation of the southern region at the Isthmus of La Paz at 3mya. The temporal frame of genetic differentiation of intraspecific haplotypes indicates that 90% of haplotypes diverged within the last 500,000years, with a population expansion 80,000years ago. Only four haplotypes diverged ∼2.2 my and occurred in the south (Hxan_36, 38 and 45), and north (Hxan_45 and 56) regions; only haplotype 45 is shared between south and north populations. These regions also have the most recent haplotypes from 12,500 to 16,200years ago, and together with high levels of genetic diversity, we suggest two refuge areas, the Northern and Southern regions. Our results indicate that the phylogeographic pattern first results from vicariance processes, then is followed by historical and recent climate fluctuations that influenced conditions on the peninsula, and it is also related to oases distribution. This study presents the first investigation of phylogeography of the peninsular' endemic Xantus' hummingbird.

Red River barrier and Pleistocene climatic fluctuations shaped the genetic structure of Microhyla fissipes complex (Anura: Microhylidae) in southern China and Indochina

South China and Indochina host striking species diversity and endemism. Complex tectonic and climatic evolutions appear to be the main drivers of the biogeographic patterns. In this study, based on the geologic history of this region, we test 2 hypotheses using the evolutionary history of Microhyla fissipes species complex. Using DNA sequence data from both mitochondrial and nuclear genes, we first test the hypothesis that the Red River is a barrier to gene flow and dispersal. Second, we test the hypothesis that Pleistocene climatic cycling affected the genetic structure and population history of these frogs. We detect 2 major genetic splits that associate with the Red River. Time estimation suggests that late Miocene tectonic movement associated with the Red River drove their diversification. Species distribution modeling (SDM) resolves significant ecological differences between sides of the Red River. Thus, ecological divergence also probably promoted and maintained the diversification. Genogeography, historical demography, and SDM associate patterns in southern China with climate changes of the last glacial maximum (LGM), but not Indochina. Differences in geography and climate between the 2 areas best explain the discovery. Responses to the Pleistocene glacial–interglacial cycling vary among species and regions.

Genetic Structure and Evolutionary History of Three Alpine Sclerophyllous Oaks in East Himalaya-Hengduan Mountains and Adjacent Regions

The East Himalaya-Hengduan Mountains (EH-HM) region has a high biodiversity and harbors numerous endemic alpine plants. This is probably the result of combined orographic and climate oscillations occurring since late Tertiary. Here, we determined the genetic structure and evolutionary history of alpine oak species (including Quercus spinosa, Quercus aquifolioides, and Quercus rehderiana) using both cytoplasmic-nuclear markers and ecological niche models (ENMs), and elucidated the impacts of climate oscillations and environmental heterogeneity on their population demography. Our results indicate there were mixed genetic structure and asymmetric contemporary gene flow within them. The ENMs revealed a similar demographic history for the three species expanded their ranges from the last interglacial (LIG) to the last glacial maximum (LGM), which was consistent with effective population sizes changes. Effects of genetic drift and fragmentation of habitats were responsible for the high differentiation and the lack of phylogeographic structure. Our results support that geological and climatic factors since Miocene triggered the differentiation, evolutionary origin and range shifts of the three oak species in the studied area and also emphasize that a multidisciplinary approach combining molecular markers, ENMs and population genetics can yield deep insights into diversification and evolutionary dynamics of species.

Genetic Structure and Evolutionary History of Three Alpine Sclerophyllous Oaks in East Himalaya-Hengduan Mountains and Adjacent Regions

The East Himalaya-Hengduan Mountains (EH-HM) region has a high biodiversity and harbors numerous endemic alpine plants. This is probably the result of combined orographic and climate oscillations occurring since late Tertiary. Here, we determined the genetic structure and evolutionary history of alpine oak species (including Quercus spinosa, Quercus aquifolioides, and Quercus rehderiana) using both cytoplasmic-nuclear markers and ecological niche models (ENMs), and elucidated the impacts of climate oscillations and environmental heterogeneity on their population demography. Our results indicate there were mixed genetic structure and asymmetric contemporary gene flow within them. The ENMs revealed a similar demographic history for the three species expanded their ranges from the last interglacial (LIG) to the last glacial maximum (LGM), which was consistent with effective population sizes changes. Effects of genetic drift and fragmentation of habitats were responsible for the high differentiation and the lack of phylogeographic structure. Our results support that geological and climatic factors since Miocene triggered the differentiation, evolutionary origin and range shifts of the three oak species in the studied area and also emphasize that a multidisciplinary approach combining molecular markers, ENMs and population genetics can yield deep insights into diversification and evolutionary dynamics of species.

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.

The potential drivers in forming avian biodiversity hotspots in the East Himalaya Mountains of Southwest China

Little has been published to describe or interpret Asian biodiversity hotspots, including those in the East Himalayan Mountains of Southwest China (HMSC), thus making necessary a review of the current knowledge. The Pliocene and Pleistocene geological and glacial histories of the Asian continent differ from those of Europe and North America, suggesting different mechanisms of speciation and extinction, and, thus, different responses to climate changes during the Quaternary glaciations. This short review summarizes potential drivers in shaping and maintaining high species richness and endemism of birds in the HMSC. The geographical location at the junction of different biogeographical realms, the wide range of habitats and climates along the extensive elevational range, the complex topography and the distinct geological history of this region have probably contributed to the evolution of an exceptionally species-rich and endemic-rich, specialized montane avian fauna. The Mountain systems in the HMSC may have provided refugia where species survived during the glacial periods and barriers for preventing species dispersal after the glacial periods. More studies are required to further test this refugia hypothesis by comparing more cold-tolerent and warm-tolerent species.

Molecular phylogenetics and phylogeographic structure of Sorex bedfordiae based on mitochondrial and nuclear DNA sequences

The southeastern margin of the Tibetan Plateau is characterized by complex topography and a discontinuous landscape, creating a sky island situation. However, the way topography shapes genetic structures and demographic histories of endemic species has not been well studied. We examined the phylogeographic pattern and demographic histories of Sorex bedfordiae, a dispersal-limited small mammal, using three nuclear genes [1977bp] and two mitochondrial genes [1794bp] with comprehensive molecular approaches. We recovered five well-supported clades whose distributions are along mountain ridges and roughly subdivided by large rivers. Demographic expansions in the middle Pleistocene were strongly supported by both nuclear and mitochondrial genes. Our results support the hypothesis that sky island topography and river systems strongly affect the genetic structure of non-aquatic terrestrial species. We further clarify that S. bedfordiae and S. cylindricauda are valid sibling species, whereas S. excelsus is most likely a geographic subspecies of S. bedfordiae.

Molecular data and ecological niche modelling reveal the Pleistocene history of a semi-aquatic bug (Microvelia douglasi douglasi) in East Asia

This study investigated the Pleistocene history of a semi-aquatic bug, Microvelia douglasi douglasi Scott, 1874 (Hemiptera: Veliidae) in East Asia. We used M. douglasi douglasi as a model species to explore the effects of historical climatic fluctuations on montane semi-aquatic invertebrate species. Two hypotheses were developed using ecological niche models (ENMs). First, we hypothesized that M. douglasi douglasi persisted in suitable habitats in southern Guizhou, southern Yunnan, Hainan, Taiwan and southeast China during the LIG. After that, the populations expanded (Hypothesis 1). As the spatial prediction in the LGM was significantly larger than in the LIG, we then hypothesized that the population expanded during the LIG to LGM transition (Hypothesis 2). We tested these hypotheses using mitochondrial data (COI+COII) and nuclear data (ITS1+5.8S+ITS2). Young lineages, relatively deep splits, lineage differentiation among mountain ranges in central, south and southwest China and high genetic diversities were observed in these suitable habitats. Evidence of mismatch distributions and neutrality tests indicate that a population expansion occurred in the late Pleistocene. The Bayesian skyline plot (BSP) revealed an unusual population expansion that likely happened during the cooling transition between LIG and LGM. The results of genetic data were mostly consistent with the spatial predictions from ENM, a finding that can profoundly improve phylogeographic research. The ecological requirements of M. douglasi douglasi, together with the geographical heterogeneity and climatic fluctuations of Pleistocene in East Asia, could have shaped this unusual demographic history. Our study contributes to our knowledge of semi-aquatic bug/invertebrate responses to Pleistocene climatic fluctuations in East Asia.

Effect of Pleistocene climatic oscillations on the phylogeography and demography of red knobby newt (Tylototriton shanjing) from southwestern China

Factors that determine the genetic structure of species in southwestern China remain largely unknown. In this study, phylogeography and demography of Tylototriton shanjing was investigated from a mitochondrial perspective to address the role of the Quaternary ice ages in shaping phylogeographic history and genetic diversity of Yunnan. A total of 146 individuals from 19 populations across the entire range of the species were collected. We detected four maternal phylogenetic lineages corresponding to four population groups, and found that major glaciation events during the Pleistocene have triggered the intra-specific divergence. Coalescent simulations indicated that the populations retreated to different refugia located in southern Yunnan, northwestern Yunnan, the border region of western Yunnan with Myanmar, and middle-western Yunnan, respectively, during previous glacial periods in the Pleistocene, and these four refugia were not retained during the Last Glacial Maximum. Population expansions occurred during the last inter-glaciation, during which ice core and pollen data indicated that the temperature and precipitation gradually increased, and declines of population sizes started after the beginning of the Last Glacial Maximum when the climate became cooler and dryer. The paleo-drainage system had no contribution to the current genetic structure and the rivers were not dispersal barriers for this salamander.