Phylogeography of ground tit (Pseudopodoces humilis) based on mtDNA: evidence of past fragmentation on the Tibetan Plateau

Patterns of endemism of the Chinese Larentiinae (Lepidoptera: Geometridae)

Chinese Larentiinae is exceptionally rich in species since it spans the Palearctic realm and the Oriental realm, comprising 715 species belonging to 130 genera. Among them, nearly half species are endemic. This taxon is considered an ideal object for the study of endemic patterns. In this study, parsimony analysis of endemicity was used to delineate areas of endemisms (AOEs), and track compatibility analysis (TCA) was used to perform the TCA. Eight AOEs were supported. Except for five AOEs around the Qinghai-Tibet plateau (QTP), two AOEs from northern China and Taiwan islands were found. Twenty-two generalized tracks were obtained from the combination of the track analysis and the AOEs. The uplift of the QTP was assumed to greatly influence these endemic patterns, presented an appearance of highly centralized endemism and diversity in the eastern and southern edges of the QTP. And the AOE of Mt. Hengduan was considered as a major center of development through recent geological history, where strong evolutionary radiation might have taken place, resulting in the migration routes described by the generalized tracks from this center to various regions in various directions.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau

Extremely heterogeneous topography and complex paleoclimate history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

Cenozoic aridization in Central Eurasia shaped diversification of toad-headed agamas (Phrynocephalus; Agamidae, Reptilia)

We hypothesize the phylogenetic relationships of the agamid genus Phrynocephalus to assess how past environmental changes shaped the evolutionary and biogeographic history of these lizards and especially the impact of paleogeography and climatic factors. Phrynocephalus is one of the most diverse and taxonomically confusing lizard genera. As a key element of Palearctic deserts, it serves as a promising model for studies of historical biogeography and formation of arid habitats in Eurasia. We used 51 samples representing 33 of 40 recognized species of Phrynocephalus covering all major areas of the genus. Molecular data included four mtDNA (COI, ND2, ND4, Cytb; 2,703 bp) and four nuDNA protein-coding genes (RAG1, BDNF, AKAP9, NKTR; 4,188 bp). AU-tests were implemented to test for significant differences between mtDNA- and nuDNA-based topologies. A time-calibrated phylogeny was estimated using a Bayesian relaxed molecular clock with nine fossil calibrations. We reconstructed the ancestral area of origin, biogeographic scenarios, body size, and the evolution of habitat preference. Phylogenetic analyses of nuDNA genes recovered a well-resolved and supported topology. Analyses detected significant discordance with the less-supported mtDNA genealogy. The position of Phrynocephalus mystaceus conflicted greatly between the two datasets. MtDNA introgression due to ancient hybridization best explained this result. Monophyletic Phrynocephalus contained three main clades: (I) oviparous species from south-western and Middle Asia; (II) viviparous species of Qinghai–Tibetan Plateau (QTP); and (III) oviparous species of the Caspian Basin, Middle and Central Asia. Phrynocephalus originated in late Oligocene (26.9 Ma) and modern species diversified during the middle Miocene (14.8–13.5 Ma). The reconstruction of ancestral areas indicated that Phrynocephalus originated in Middle East–southern Middle Asia. Body size miniaturization likely occurred early in the history of Phrynocephalus. The common ancestor of Phrynocephalus probably preferred sandy substrates with the inclusion of clay or gravel. The time of Agaminae radiation and origin of Phrynocephalus in the late Oligocene significantly precedes the landbridge between Afro-Arabia and Eurasia in the Early Miocene. Diversification of Phrynocephalus coincides well with the mid-Miocene climatic transition when a rapid cooling of climate drove progressing aridification and the Paratethys salinity crisis. These factors likely triggered the spreading of desert habitats in Central Eurasia, which Phrynocephalus occupied. The origin of the viviparous Tibetan clade has been associated traditionally with uplifting of the QTP; however, further studies are needed to confirm this. Progressing late Miocene aridification, the decrease of the Paratethys Basin, orogenesis, and Plio–Pleistocene climate oscillations likely promoted further diversification within Phrynocephalus. We discuss Phrynocephalus taxonomy in scope of the new analyses.

The geography and timing of genetic divergence in the lizard Phrynocephalus theobaldi on the Qinghai-Tibetan plateau

The Qinghai-Tibetan Plateau (QTP) represents one of the earth’s most significant physical features and there is increasing interest in the historical generation of biodiversity within this region. We hypothesized that there should be clear geographically coherent genetic structuring within one of the world’s highest altitude lizards, Phrynocephalus theobaldi, due to considerable historical population fragmentation in this environment. This was tested using a major mitochondrial DNA (mtDNA) survey and sequencing of two nuclear markers (AME and RAG-1) from P. theobaldi, from across the southern QTP. A Bayesian method (BPEC) was used to detect four geographically structured mtDNA clusters. A Bayesian phylogenetic tree, together with associated dating analyses, supported four corresponding evolutionary lineages with a timing of 3.74–7.03 Ma for the most basal P. theobaldi split and Pliocene splits of 2.97–5.79 Ma and 2.40–5.39 Ma in the two daughter lineages. Himalayan uplift and changes in the Jilong basin may have contributed to these divergences, but uplift of the Gangdese mountains is rejected due to its timing. The nuclear markers appeared to be sorted between the four mtDNA groups, and species delimitation analyses supported the four phylogeographical groups as candidate species. The study contributes to our understanding of biodiversity on the QTP.

Evolution of body morphology and beak shape revealed by a morphometric analysis of 14 Paridae species

BACKGROUND

Morphological characters of birds reflect their adaptive evolution and ecological requirements and are also relevant to phylogenetic relationships within a group of related species. The tits (Paridae) are known to be outwardly homogeneous in shape, with one aberrant member, the Ground Tit (Pseudopodoces humilis), which is quite different from its relatives in both body morphology and beak shape. We combined traditional measurements and geometric morphometrics to quantify the variation in body morphology and beak shape of 14 Paridae species distributed in China. Based on these results, we sought to assess the contribution of phylogeny, altitude and species interactions to the evolution of morphological traits.

RESULTS

The basic features for discriminating among the 14 species studied here were overall body size, the ratio of body and tail length to culmen and tarsus length, and beak shape (long/slender/pointy vs. short/robust/blunt). These dimensions clearly separate Ps. humilis and Melanochlora sultanea from the other species in shape space. Body length and PC3 of beak shape (round outline vs. straight outline) show significant phylogenetic signals. Across 14 species, altitude is related to tarsus, culmen length and PC1 of beak shape. Within Parus major, altitude is related to body weight, body length, culmen length and PC1 of body morphology. Morphological distances and geographic distances among species are positively correlated.

CONCLUSIONS

The body morphology of Paridae species shows extensive evolutionary changes, while their beak has mainly evolved along the long/slender/pointy vs. short/robust/blunt dimension. Only body length and beak curvature show a phylogenetic signal. Altitude correlates with multiple traits both across and within species, suggesting that altitude is an important factor in promoting morphological divergence. The deviant appearance of Ps. humilis corresponds to its foraging and feeding adaptations to high-altitude steppe habitats. Our results also show a higher level of morphological divergence with greater difference in distribution ranges among the Paridae species involved in this study.

Species history and divergence times of viviparous and oviparous Chinese toad-headed sand lizards (Phrynocephalus) on the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) is an important biogeographical area and has recently become a focus for biodiversity studies. Phyrnocephalus lizards form a widespread Eurasian group with oviparous and viviparous reproductive modes, but two previous mtDNA studies of species from around the QTP have provided different phylogenetic hypotheses. We analysed three loci (mtDNA, RAG-1, AME) from all recognised Chinese Phrynocephalus species to reconstruct the speciation history of the group and to estimate species divergence times. The effects of mtDNA partitioning strategy on phylogenetic inference were examined. Bayes factor comparisons of marginal likelihoods (mLs) estimated using stepping-stone sampling revealed that partitioning strategy had a major impact on mL. Nevertheless, it had a negligible effect on the inferred tree topology. The impact of hard-bound uniform or equivalent soft-bound gamma speciation time calibration priors as well as the use of a fixed topology (as opposed to integration over all possible species histories) on divergence time estimation were also assessed, and found to have little impact on posterior estimates. All three gene trees and the species tree supported the hypothesis that the Chinese species form oviparous and viviparous sister clades. This was in agreement with an early mtDNA study but differed from a subsequent reanalysis of the mtDNA data. Inclusion of mtDNA from more widely distributed Phrynocephalus (from previous studies) indicates that the oviparous P. interscapularis from Central Asia lies outside the clade of Chinese viviparous and oviparous species, but that other Asian oviparous species lie within the Chinese oviparous clade. The median of the posterior on the divergence time of Chinese oviparous and viviparous species was 9.7 Ma ago (95% interval: 7.2-13.0 Ma ago), which coincides with major uplifting of the QTP and indicates that viviparity evolved when this clade became restricted to regions of high elevation. We also found that cladogenesis within the viviparous clade began around 5 Ma ago whereas those in the oviparous clade began around 8.6 Ma ago. We establish more robust estimates of divergence times and relationships within this important group and so provide improved insights into the origins of Phrynocephalus diversity across the QTP.

Effects of the Qinghai-Tibetan Plateau uplift and environmental changes on phylogeographic structure of the Daurian Partridge (Perdix dauuricae) in China

Data from eight microsatellite loci were used to assess the evolutionary demographic processes of Daurian Partridges from 285 individuals distributed among 23 populations throughout much of the species' distribution range in China. Phylogenetic analysis using microsatellite DNA data indicated that the species were geographically structured and were split into two distinctive phylogroups A and B, splitting partridges geographically from high plateaus vs. lowlands in northern China with high bootstrap support in phylogenetic trees. Each of cluster A and cluster B was comprised with two subgroups respectively. Cluster A was further subdivided into the Loess Plateau group and the Qinghai-Tibetan group. Two subgroups of cluster B corresponded to populations of Liupan Mountains and northern China, respectively. This subdivision was confirmed by Bayesian clustering and population assignment analyses of microsatellite genotypes which indicated strong isolation of these groups/populations. Our results suggested that the vicariance patterns of genetic structures of Daurian Partridge may have resulted from: (1) the uplift of the Qinghai-Tibetan Plateau and the Loess Plateau followed by the increased aridification and desertification in northern China in Middle Pleistocene, (2) the monsoon and the existence of Qinling Mountains and Liupan Mountains, and (3) the glacial cycles from the late of Middle Pleistocene to early of Late Pleistocene. In addition, BOTTLENECK analysis indicated that Daurian Partridge had been experienced recent evolutionary bottlenecks.

Phylogeography of the South China field mouse (Apodemus draco) on the southeastern Tibetan Plateau reveals high genetic diversity and glacial refugia

The southeastern margin of the Tibetan Plateau (SEMTP) is a particularly interesting region due to its topographic complexity and unique geologic history, but phylogeographic studies that focus on this region are rare. In this study, we investigated the phylogeography of the South China field mouse, Apodemus draco, in order to assess the role of geologic and climatic events on the Tibetan Plateau in shaping its genetic structure. We sequenced mitochondrial cytochrome b (cyt b) sequences in 103 individuals from 47 sampling sites. In addition, 23 cyt b sequences were collected from GenBank for analyses. Phylogenetic, demographic and landscape genetic methods were conducted. Seventy-six cyt b haplotypes were found and the genetic diversity was extremely high (π = 0.0368; h = 0.989). Five major evolutionary clades, based on geographic locations, were identified. Demographic analyses implied subclade 1A and subclade 1B experienced population expansions at about 0.052-0.013 Mya and 0.014-0.004 Mya, respectively. The divergence time analysis showed that the split between clade 1 and clade 2 occurred 0.26 Mya, which fell into the extensive glacial period (EGP, 0.5-0.17 Mya). The divergence times of other main clades (2.20-0.55 Mya) were congruent with the periods of the Qingzang Movement (3.6-1.7 Mya) and the Kun-Huang Movement (1.2-0.6 Mya), which were known as the most intense uplift events in the Tibetan Plateau. Our study supported the hypothesis that the SEMTP was a large late Pleistocene refugium, and further inferred that the Gongga Mountain Region and Hongya County were glacial refugia for A. draco in clade 1. We hypothesize that the evolutionary history of A. draco in the SEMTP primarily occurred in two stages. First, an initial divergence would have been shaped by uplift events of the Tibetan Plateau. Then, major glaciations in the Pleistocene added complexity to its demographic history and genetic structure.

Profound climatic effects on two East Asian black-throated tits (Ave: Aegithalidae), revealed by ecological niche models and phylogeographic analysis

Although a number of studies have assessed the effects of geological and climatic changes on species distributions in East Asian, we still have limited knowledge of how these changes have impacted avian species in south-western and southern China. Here, we aim to study paleo-climatic effects on an East Asian bird, two subspecies of black-throated tit (A. c. talifuensis–concinnus) with the combined analysis of phylogeography and Ecological Niche Models (ENMs). We sequenced three mitochondrial DNA markers from 32 populations (203 individuals) and used phylogenetic inferences to reconstruct the intra-specific relationships among haplotypes. Population genetic analyses were undertaken to gain insight into the demographic history of these populations. We used ENMs to predict the distribution of target species during three periods; last inter-glacial (LIG), last glacial maximum (LGM) and present. We found three highly supported, monophyletic MtDNA lineages and different historical demography among lineages in A. c. talifuensis–concinnus. These lineages formed a narrowly circumscribed intra-specific contact zone. The estimated times of lineage divergences were about 2.4 Ma and 0.32 Ma respectively. ENMs predictions were similar between present and LGM but substantially reduced during LIG. ENMs reconstructions and molecular dating suggest that Pleistocene climate changes had triggered and shaped the genetic structure of black-throated tit. Interestingly, in contrast to profound impacts of other glacial cycles, ENMs and phylogeographic analysis suggest that LGM had limited effect on these two subspecies. ENMs also suggest that Pleistocene climatic oscillations enabled the formation of the contact zone and thus support the refuge theory.

Freshwater biogeography and limnological evolution of the Tibetan Plateau--insights from a plateau-wide distributed gastropod taxon (Radix spp.)

Background

The Tibetan Plateau is not only the highest and largest plateau on earth; it is also home to numerous freshwater lakes potentially harbouring endemic faunal elements. As it remains largely unknown whether these lakes have continuously existed during the Last Glacial Maximum (LGM), questions arise as to whether taxa have been able to exist on the plateau since before the latest Pleistocene, from where and how often the plateau was colonized, and by which mechanisms organisms conquered remote high altitude lentic freshwater systems. In this study, species of the plateau-wide distributed freshwater gastropod genus Radix are used to answer these biogeographical questions.

Methodology/Principal Findings

Based on a broad spatial sampling of Radix spp. on the Tibetan Plateau, and phylogenetic analyses of mtDNA sequence data, three probably endemic and one widespread major Radix clade could be identified on the plateau. Two of the endemic clades show a remarkably high genetic diversity, indicating a relatively great phylogenetic age. Phylogeographical analyses of individuals belonging to the most widely distributed clade indicate that intra-plateau distribution cannot be explained by drainage-related dispersal alone.

Conclusions/Significance

Our study reveals that Radix spp. persisted throughout the LGM on the Tibetan Plateau. Therefore, we assume the continuous existence of suitable water bodies during that time. The extant Radix diversity on the plateau might have been caused by multiple colonization events combined with a relatively long intra-plateau evolution. At least one colonization event has a Palaearctic origin. In contrast to freshwater fishes, passive dispersal, probably by water birds, might be an important mechanism for conquering remote areas on the plateau. Patterns found in Radix spp. are shared with some terrestrial plateau taxa, indicating that Radix may be a suitable model taxon for inferring general patterns of biotic origin, dispersal and survival on the Tibetan Plateau.

Phylogeography of the endemic Gymnocypris chilianensis (Cyprinidae): sequential westward colonization followed by allopatric evolution in response to cyclical Pleistocene glaciations on the Tibetan Plateau

The schizothoracine Gymnocypris chilianensis is restricted to the Shiyang, Ruoshui and Shule Rivers, listed from east to west, along the northeast edge of the Tibetan Plateau. This distribution provides a valuable system to test hypotheses about postglacial colonization. We used mitochondrial DNA sequence data (a control region and the cytochrome b gene; 1894 bp) to assess the phylogeographic structure of this species based on 278 specimens sampled from throughout the species' entire geographical range. We found three lineages corresponding geographically to the three rivers, suggesting three independent glacial differentiation centers within the northeast edge of the Tibetan Plateau. The phylogenetic analysis suggested that the Shiyang River population forms a lineage that separated from the other populations of G. chilianensis at the basal phylogenetic split within this species. The molecular data further demonstrated a clear pattern of decreasing genetic diversity from the eastern Shiyang River towards the central Ruoshui River and western Shule River lineages, a pattern consistent with sequential western colonization. We therefore propose a phylogeographic scenario for G. chilianensis of a gradual westerly expansion from the Shiyang River population along the northeast edge of the Tibetan Plateau, with subsequent allopatric evolution at approximately 0.37 and 0.05 million years ago (Ma), through at least two glacial maxima. Together with the genetic evidence reported in other species, our findings suggest that this common biogeographic pattern emphasizes the importance of the northeastern edge region of the Tibetan Plateau as a hotspot of genetic diversity for some taxa.

Genetic homogeneity among colonies of the white-nest swiftlet (Aerodramus fuciphagus) in Thailand

The white-nest swiftlet, Aerodramus fuciphagus, originally lived in large colonies in natural caves, but now it also occurs in man-made buildings. We investigated the patterns of genetic differentiation in two mitochondrial DNA genes (cyt-b and ND2) and eight microsatellite loci among and within colonies of A. fuciphagus from across recently established man-made colonies in Thailand. Ten white-nest swiftlet colonies were sampled along the coast of the Gulf of Thailand and the Andaman Sea in Thailand during 2003-2006. The genetic diversity of mtDNA was very low, and few significant PhiST values were found between pairs of colonies. Analyses of haplotype relationships did not show genetic structure across the sampled distribution. The level of genetic diversity for microsatellite loci was high, but FST values were not significant. However, due to small sample sizes for some colonies that could limit conclusions on genetic differentiation from PhiST and FST, we also analyzed the microsatellite data using STRUCTURE and found that number of subpopulations of white-nest swiftlets in sampled colonies was one. The lack of genetic differentiation among swiftlet house colonies could be a result of high gene flow between colonies and large population sizes. Our results suggest that A. fuciphagus living in recently established man-made colonies in Thailand should be considered members of a single panmictic population. Future work will be necessary to determine whether this panmixia is stable or a temporary result of the recent explosive expansion of the number of colonies, and comparisons to natural colonies may provide an understanding of mechanisms producing the lack of genetic structure in swiftlet house colonies.