Comparison of phylogeographic structure and population history of two Phrynocephalus species in the Tarim Basin and adjacent areas

All roads lead to Rome: the plasticity of gut microbiome drives the extensive adaptation of the Yarkand toad-headed agama (Phrynocephalus axillaris) to different altitudes

The gut microbiome was involved in a variety of physiological processes and played a key role in host environmental adaptation. However, the mechanisms of their response to altitudinal environmental changes remain unclear. In this study, we used 16S rRNA sequencing and LC-MS metabolomics to investigate the changes in the gut microbiome and metabolism of the Yarkand toad-headed agama (Phrynocephalus axillaris) at different altitudes (-80 m to 2000 m). The results demonstrated that Firmicutes, Bacteroidetes, and Proteobacteria were the dominant phylum, Lachnospiraceae and Oscillospiraceae were the most abundant family, and the low-altitude populations had higher richness than high-altitude populations; Akkermansiaceae appeared to be enriched in high-altitude populations and the relative abundance tended to increase with altitude. The gut microbiome of three populations of P. axillaris at different altitudes was clustered into two different enterotypes, low-altitude populations and high-altitude populations shared an enterotype dominated by Akkermansia, Kineothrix, Phocaeicola; intermediate-altitude populations had an enterotype dominated by Mesorhizobium, Bradyrhizobium. Metabolites involved in amino acid and lipid metabolism differed significantly at different altitudes. The above results suggest that gut microbiome plasticity drives the extensive adaptation of P. axillaris to multi-stress caused by different altitudes. With global warming, recognizing the adaptive capacity of wide-ranging species to altitude can help plan future conservation strategies.

Revealing the coexistence of differentiation and communication in an endemic hare, Lepus yarkandensis (Mammalia, Leporidae) using specific-length amplified fragment sequencing

BACKGROUND

The Yarkand hare (Lepus yarkandensis Günther, 1875) is endemic to oasis and desert areas around the Tarim Basin in the Xinjiang Uyghur Autonomous Region of northwest China; however, genome-wide information for this species remains limited. Moreover, the genetic variation, genetic structure, and phylogenetic relationships of Yarkand hare from the plateau mountain regions have not been reported. Thus, we used specific-length amplified fragment sequencing (SLAF-seq) technology to evaluate the genetic diversity of 76 Yarkand hares from seven geographic populations in the northern and southwestern parts of the Tarim Basin to investigate single-nucleotide polymorphism (SNP) marker-based population differentiation and evolutionary processes. Selective sweep analysis was conducted to identify genetic differences between populations.

RESULTS

Using SLAF-seq, a total of 1,835,504 SNPs were initially obtained, of which 308,942 high-confidence SNPs were selected for further analysis. Yarkand hares exhibited a relatively high degree of genetic diversity at the SNP level. Based on pairwise FST estimates, the north and southwest groups showed a moderate level of genetic differentiation. Phylogenetic tree and population structure analyses demonstrated evident systematic phylogeographical structure patterns consistent with the geographical distribution of the hares. Hierarchical analysis of molecular variation further indicated that genetic variation was mainly observed within populations. Low to moderate genetic differentiation also occurred among populations despite a common genomic background, likely due to geographical barriers, genetic drift, and differential selection pressure of distinct environments. Nevertheless, the observed lineage-mixing pattern, as indicated by the evolutionary tree, principal component analysis, population structure, and TreeMix analyses, suggests a certain degree of gene flow between the north and southwest groups. This may be related to the migration of hares to high-altitude water sources southwest of the basin during glacial climatic oscillations, as well as river re-diffusion and oasis restoration in the basin following the glacial period. We also identified candidate genes, and their associated gene ontology terms and pathways, related to the adaptation of Yarkand hares to different environmental habitats.

CONCLUSIONS

The identified genome-wide SNPs, genetic diversity, and population structure of Yarkand hares expand our understanding of the genetic background of this endemic species and provide valuable insights into its environmental adaptation, allowing for further exploration of the underlying mechanisms.

Environmental and geological changes in the Tarim Basin promoted the phylogeographic formation of Phrynocephalus forsythii (Squamata: Agamidae)

Environmental factors can promote genetic divergence among populations. The Tarim Basin has experienced a series of environmental and geological changes since the late Pliocene. Phrynocephalus forsythii (Agamidae) has a continuous and circular distribution around the Tarim Basin, and the evolutionary history of this species remains unclear. In the present study, mitochondrial DNA and single nucleotide polymorphism sequences were obtained from 195 P. forsythii samples in the Tarim Basin to examine the phylogeographic structure and evolutionary history of this species. All populations of P. forsythii formed three distinct clusters: the Minfeng, low-elevation, and high-elevation groups. The Minfeng group was the first to separate from all other groups at 4.26 Ma after the original desert environment emerged in the Tarim Basin. The uplift of the Kunlun Mountains created a new high-altitude environment along the southwestern edge of the basin and promoted the divergence of low-elevation and high-elevation groups at 3.67 Ma. Subsequently, the gradual development of deserts and rivers in the Tarim Basin accelerated the dispersal of P. forsythii, eventually forming the current distribution pattern. Our results indicate that ecological separation caused by environmental and geological changes in the Tarim Basin was the main cause for the genetic divergence of P. forsythii.

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.

Continental Refugium in the Mongolian Plateau during Quaternary Glacial Oscillations: Phylogeography and Niche Modelling of the Endemic Desert Hamster, Phodopus roborovskii

The Mongolian Plateau (MP), which is situated in the interior of Asia and possesses a typical continental climate, experienced harsh climatic conditions during the Quaternary glacial fluctuations. Although these events likely had huge impacts on the local animal populations, the current effects have hardly been explored. To investigate whether the MP supported a refugium along an oceanic-continental gradient (ROCG), and whether this refugium was glacial or interglacial, we investigated the demographic and phylogeographic history of an endemic mammal species, the desert hamster Phodopus roborovskii. We reconstructed the demographic variation, the phylogeographic diffusion, and modelled the potential habitat during historical periods. The genetic diversity in the MP was the highest among all the localities, and the MP was a suitable habitat throughout the modelled historical periods. A phylogeographic diffusion analysis emphasized the importance of the MP as the centre of origin, preservation and spread for P. roborovskii. The homogeneous landscape provided the opportunity for a wide gene flow, which resulted in low resolution of the phylogenetic relationships. Moreover, P. roborovskii was favoured by the interglacial condition, with both its demographical and geographical ranges expanded within the interglacial periods. The range variation from the Last Glacial Maximum to the current condition reflects a distinct longitudinal shift, while both ranges largely contracted from that of the Last Interglacial. Our results support that the MP served as a refugium and spread centre for P. roborovskii during the Quaternary climate fluctuations. The interglacial expansion and the longitudinal shifts highlighted the important effects of precipitations on the distribution range of species adapted to arid and semi-arid during glacial oscillations.

Retreating or standing: responses of forest species and steppe species to climate change in arid Eastern Central Asia

Background

The temperature in arid Eastern Central Asia is projected to increase in the future, accompanied by increased variability of precipitation. To investigate the impacts of climate change on plant species in this area, we selected two widespread species as candidates, Clematis sibirica and C. songorica, from montane coniferous forest and arid steppe habitats respectively.

Methodology/Principal Findings

We employed a combined approach of molecular phylogeography and species distribution modelling (SDM) to predict the future responses of these two species to climate change, utilizing evidence of responses from the past. Genetic data for C. sibirica shows a significant phylogeographical signal (N_ST > F_ST, P<0.05) and demographic contraction during the glacial-interglacial cycles in the Pleistocene. This forest species would likely experience range reduction, though without genetic loss, in the face of future climate change. In contrast, SDMs predict that C. songorica, a steppe species, should maintain a consistently stable potential distribution under the Last Glacial Maximum (LGM) and the future climatic conditions referring to its existing potential distribution. Molecular results indicate that the presence of significant phylogeographical signal in this steppe species is rejected and this species contains a high level of genetic differentiation among populations in cpDNA, likely benefiting from stable habitats over a lengthy time period.

Conclusions/Significance

Evidence from the molecular phylogeography of these two species, the forest species is more sensitive to past climate changes than the steppe species. SDMs predict that the forest species will face the challenge of potential range contraction in the future more than the steppe species. This provides a perspective on ecological management in arid Eastern Central Asia, indicating that increased attention should be paid to montane forest species, due to their high sensitivity to disturbance.

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.

New insights into the intricate taxonomy and phylogeny of the Sylvia curruca complex

We use the mitochondrial cytochrome b from 213 individuals and the three nuclear introns BRM 15, myoglobin 2 and ODC 6-7 from a smaller subsample to evaluate the taxonomy of the Lesser Whitethroat Sylvia curruca (Aves, Passeriformes, Sylviidae) complex, which has long been controversial. We sequenced type material of the taxa althaea, blythi, margelanica and minula, and used topotypical material of caucasica, chuancheica, curruca and telengitica. The nuclear introns fail to resolve the complex, but cytochrome b recovers six major clades, revealing genetically identifiable populations corresponding to previously named taxa, and we propose that the names althaea, blythi, curruca, halimodendri, margelanica and minula, respectively, should be used for these. The margelanica clade is suggested to have a more extensive distribution than previously known, including both the taxon telengitica and a population in eastern Mongolia. The taxon minula is found to have a more restricted range than generally believed, only breeding in China. According to the mitochondrial gene tree, there is a basal dichotomy, with the taxa althaea, blythi, halimodendri and margelanica being part of one clade, well separated from a clade containing curruca and minula. Dating analysis suggests that a basal divergence separating curruca and minula from the other four taxa occurred between 4.2 and 7.2 mya; these two then diverged between 2.3 and 4.4 mya. The splits between the althaea, blythi, halimodendri and margelanica lineages is inferred to have occurred later, approximately between 1.0 and 2.5 mya (all 95% HPD). The nucleotide data suggest significant departure from demographic equilibrium in blythi (clade 1a), halimodendri (clade 2a) and minula, whereas tendencies are weaker for other clades. We propose that the names althaea, blythi, curruca, halimodendri, margelanica and minula should be used for the major clades. However, whether these are treated as subspecies or species is largely a matter of species definition and is not resolved by our data.