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

Impacts of Quaternary glaciation, geological history and geography on animal species history in continental East Asia: A phylogeographic review

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

Phylogeography of Triplophysa stenura (Nemacheilidae): Responded to the Mid-Pleistocene Climate Transition in the Qinghai-Tibetan Plateau

Pleistocene climatic oscillations have greatly influenced the evolutionary histories and distribution patterns of most extant species. In this study, the phylogeographic patterns and evolutionary history of Triplophysa stenura were investigated. A total of 359 individuals from 19 locations covering the species' distribution range were collected, and two mitochondrial genes (COI and Cyt b) and the recombination activating protein 1 gene (Rag1) were analyzed. Two deeply divergent lineages, corresponding geographically to the northern and southern Tanggula Mountain, were observed, and shallow divergences were found within the southern and central Qinghai-Tibetan plateau (QTP). The estimated divergence time between the northern and southern Tanggula lineages was around 1.07 Mya. Within the southern Tanggula Mountain lineages, the Nu and Subansiri rivers populations were split about 0.74 Mya, and the southwestern and central QTP populations were divided with the southeastern QTP ones approximately 0.67 Mya. The divergence times of the lineages were matched with two major glaciations in QTP (the Xixiabangma Glaciation, 1.17-0.80 Mya and Kunlun Glaciation, 0.72-0.62 Mya). Together with demographic history analyses, our results highlighted that major glacial cycles during the mid-Pleistocene played a major role in sculpting the distribution pattern of T. stenura, and led to the gene homogenization crossing the drainage populations in the southwestern and central QTP.

Evolutionary history of the oriental fire-bellied toad (Bombina orientalis) in Northeast China

The evolutionary history of a species is generally affected by the combination of geological events and climate fluctuations. By analyzing the population features, genetic structure and the effective population historical dynamics of existing species, the population evolutionary history can be reestablished. In recent years, geological evidence shows that the Yilan-Yitong fault zone located in Northeast Asia experienced strong and frequent geological changes in the late Quaternary period. Species population history has been shaped by the combination of the complex climatic conditions of the Quaternary and Pleistocene glacial interglacial cycles and palaeogeological events in Northeast Asia and it has become a research focus for evolutionary biology researchers. In this study, mitochondrial and microsatellite molecular markers were used to reveal the population features, genetic structure, and the effective population historical dynamics of the Oriental fire-bellied toad (Bombina orientalis). The results showed that the strong seismic activity of the Yilan-Yitong fault zone in the late Quaternary period was the main reason for the population differentiation of Oriental fire-bellied toad in northeast China. The Quaternary Pleistocene glacial interglacial cycles led to the significant bottleneck effect of the western population located in the Maoer mountain area. As a result, the western population has low genetic diversity. Recent gene flow between eastern and western populations and historical evidence of population expansion proved that the dispersal behavior of the western populations was the main cause of the low genetic diversity and mitochondrial and nuclear discordance. Human economic activity may be the mainly driving factor. These evidences showed that the comprehensive influence of geology, climate, human activities and other factors should be considered in the process of exploring the evolutionary history of species.

Whole-Genome Diversification Analysis of the Hornbeam Species Reveals Speciation and Adaptation Among Closely Related Species

Speciation is the key evolutionary process for generating biological diversity and has a central place in evolutionary and ecological research. How species diverge and adapt to different habitats is one of the most exciting areas in speciation studies. Here, we sequenced 55 individuals from three closely related species in the genus Carpinus: Carpinus tibetana, Carpinus monbeigiana, and Carpinus mollicoma to understand the strength and direction of gene flow and selection during the speciation process. We found low genetic diversity in C. tibetana, which reflects its extremely small effective population size. The speciation analysis between C. monbeigiana and C. mollicoma revealed that both species diverged ∼1.2 Mya with bidirectional gene flow. A total of 291 highly diverged genes, 223 copy number variants genes, and 269 positive selected genes were recovered from the two species. Genes associated with the diverged and positively selected regions were mainly involved in thermoregulation, plant development, and response to stress, which included adaptations to their habitats. We also found a great population decline and a low genetic divergence of C. tibetana, which suggests that this species is extremely vulnerable. We believe that the current diversification and adaption study and the important genomic resource sequenced herein will facilitate the speciation studies and serve as an important methodological reference for future research.

Historical biogeography of Acer L. (Sapindaceae): genetic evidence for Out-of-Asia hypothesis with multiple dispersals to North America and Europe

Biogeography is the study of where, when, and how modern species evolved and diversified. Acer L. (maple) is one of the most diverse and widespread genera in the Northern Hemisphere. It comprises 124–156 species in the world, approximately 80% species of Acer are native in Asia. The current diversity center of Acer is not congruent with the distribution of the oldest fossils of the genus. Therefore, we herein used 84 species and subspecies to reconstruct the phylogeny and investigate the biogeographic history of Acer using nuclear ITS and three cpDNA fragments (psbA-trnH spacer, rpl16 intron, and trnL-trnF spacer) with maximum likelihood, maximum parsimony, and Bayesian inference methods. The analyses showed that the current diversity center and the origin center of Acer is Asia. Additionally, the North American and Euro-Mediterranean species originated from multiple sources from Asia via the North Atlantic Land Bridge and the Bering Land Bridge, and intercontinental migration has mainly occurred since the Miocene. This study not only provides a novel insight of the origin and dispersal routes of Acer but also exemplifies how past climatic changes affect the diversification-rates of Northern Hemisphere forest trees.

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.

Yangtze River, an insignificant genetic boundary in tufted deer (Elaphodus cephalophus): the evidence from a first population genetics study

Great rivers were generally looked at as the geographical barrier to gene flow for many taxonomic groups. The Yangtze River is the third largest river in the world, and flows across South China and into the East China Sea. Up until now, few studies have been carried out to evaluate its effect as a geographical barrier. In this study, we attempted to determine the barrier effect of the Yangtze River on the tufted deer (Elaphodus cephalophus) using the molecular ecology approach. Using mitochondrial DNA control region (CR) sequences and 13 nuclear microsatellite loci, we explored the genetic structure and gene flow in two adjacent tufted deer populations (Dabashan and Wulingshan populations), which are separated by the Yangtze River. Results indicated that there are high genetic diversity levels in the two populations, but no distinguishable haplotype group or potential genetic cluster was detected which corresponded to specific geographical population. At the same time, high gene flow was observed between Wulingshan and Dabashan populations. The tufted deer populations experienced population decrease from 0.3 to 0.09 Ma BP, then followed by a distinct population increase. A strong signal of recent population decline (T = 4,396 years) was detected in the Wulingshan population by a Markov-Switching Vector Autoregressions(MSVAR) process population demography analysis. The results indicated that the Yangtze River may not act as an effective barrier to gene flow in the tufted deer. Finally, we surmised that the population demography of the tufted deer was likely affected by Pleistocene climate fluctuations and ancient human activities.

Genetics, morphology and ecology reveal a cryptic pika lineage in the Sikkim Himalaya

Asian pika species are morphologically ∼similar and have overlapping ranges. This leads to uncertainty and species misidentification in the field. Phylogenetic analyses of such misidentified samples leads to taxonomic ambiguity. The ecology of many pika species remains understudied, particularly in the Himalaya, where sympatric species could be separated by elevation and/or substrate. We sampled, measured, and acquired genetic data from pikas in the Sikkim Himalaya. Our analyses revealed a cryptic lineage, Ochotona sikimaria, previously reported as a subspecies of O. thibetana. The results support the elevation of this lineage to the species level, as it is genetically divergent from O. thibetana, as well as sister species, O. cansus (endemic to central China) and O. curzoniae (endemic to the Tibetan plateau). The Sikkim lineage diverged from its sister species' about 1.7-0.8myrago, coincident with uplift events in the Himalaya. Our results add to the recent spate of cryptic diversity identified from the eastern Himalaya and highlight the need for further study within the Ochotonidae.

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri, Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.

The Phylogeography and Population Demography of the Yunnan Caecilian (Ichthyophis bannanicus): Massive Rivers as Barriers to Gene Flow

Ichthyophis bannanicus is the only caecilian species in China. In this study, the phylogeography and population demography of I. bannanicus were explored, based on the mitochondrial DNA genes (cyt b and ND2) and 15 polymorphic microsatellite loci. Altogether 158 individuals were collected from five populations in Yunnan province, Guangxi province, Guangdong province, and Northern Vietnam. Phylogeographical and population structure analysis identified either two groups (Xishuangbanna, Northern Vietnam-Yulin-Yangchun-Deqing) or three groups (Xishuangbanna, Northern Vietnam-Yulin-Yangchun, and Deqing), indicating that the Red River and Pearl River systems may have acted as gene-flow barriers for I. bannanicus. Historical population expansion that happened 15–17 Ka ago was detected for mtDNA data and was possibly triggered by warmer weather after the Last Glacial Maximum. However, the Bayesian simulations of population history based on microsatellite data pinpointed population decline in all populations since 19,123 to 1,029 years ago, demonstrating a significant influence of anthropogenic habitat alteration on I. bannanicus.

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.