Genetic consequences of postglacial colonization by the endemic Yarkand hare (Lepus yarkandensis) of the arid Tarim Basin

Selective sweep analysis of the adaptability of the Yarkand hare (Lepus yarkandensis) to hot arid environments using SLAF-seq

The Yarkand hare (Lepus yarkandensis) inhabits arid desert areas and is endemic to China. It has evolved various adaptations to survive in hot arid environments, including stress responses, the ability to maintain water homeostasis and heat tolerance. Here, we performed a selective sweep analysis to identify the candidate genes for adaptation to hot arid environments in the Yarkand hare. A total of 397 237 single-nucleotide polymorphisms were obtained from 80 Yarkand hares, which inhabit hot arid environments, and 36 Tolai hares (Lepus tolai), which inhabit environments with a mild climate, via specific-locus amplified fragment sequencing. We identified several candidate genes that were associated with the heat stress response (HSPE1), oxidative stress response (SLC23A and GLRX2), immune response (IL1R1 and IRG1), central nervous system development (FGF13, THOC2, FMR1 and MECP2) and regulation of water homeostasis (CDK1) according to fixation index values and θπ ratios in the selective sweep analysis, and six of these genes (GLRX2, IRG1, FGF13, FMR1, MECP2 and CDK1) are newly discovered genes. To the best of our knowledge, this is the first study to identify candidate genes for adaptation to hot arid environments in the Yarkand hare. The results of this study enhance our understanding of the adaptation of the Yarkand hare to hot arid environments and will aid future studies aiming to functionally verify these candidate genes.

Population genetics analysis of Tolai hares (Lepus tolai) in Xinjiang, China using genome-wide SNPs from SLAF-seq and mitochondrial markers

The main topic of population genetics and evolutionary biology is the influence of the ecological environment, geographical isolation, and climatic factors on population structure and history. Here, we estimated the genetic diversity, genetic structure, and population history of two subspecies of Tolai hares (Lepus tolai Pallas, 1778), L. t. lehmanni inhabiting Northern and Northwest Xinjiang and L. t. centrasiaticus inhabiting Central and Eastern Xinjiang using SNP of specific-length amplified fragment sequencing (SLAF-seq) and four mitochondrial DNA (mtDNA). Our results showed a relatively high degree of genetic diversity for Tolai hares, and the diversity of L. t. lehmanni was slightly higher than that of L. t. centrasiaticus, likely due to the more favorable ecological environment, such as woodlands and plains. Phylogenetic analysis from SNP and mtDNA indicated a rough phylogeographical distribution pattern among Tolai hares. Strong differentiation was found between the two subspecies and the two geographical groups in L. t. centrasiaticus, possibly due to the geographical isolation of mountains, basins, and deserts. However, gene flow was also detected between the two subspecies, which might be attributed to the Tianshan Corridor and the strong migration ability of hares. Tolai hare population differentiation occurred at approximately 1.2377 MYA. Population history analysis based on SNP and mtDNA showed that the Tolai hare population has a complex history and L. t. lehmanni was less affected by the glacial event, possibly because its geographic location and terrain conditions weaken the drastic climate fluctuations. In conclusion, our results indicated that the joint effect of ecological environment, geographic events, and climatic factors might play important roles in the evolutionary process of L. t. lehmanni and L. t. centrasiaticus, thus resulting in differentiation, gene exchange, and different population history.

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.

Complete mitochondrial genome sequence of Lepus yarkandensis Günther, 1875 (Lagomorpha, Leporidae): characterization and phylogenetic analysis

Lepusyarkandensis is a national second-class protected animal endemic to China and distributed only in the hot and arid Tarim Basin in Xinjiang. We sequenced and described the complete mitogenome of L.yarkandensis to analyze its characteristics and phylogeny. The species’ DNA is a 17,047 bp circular molecule that includes 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. The overall base composition was as follows: A, 31.50%; T, 29.40%; G, 13.30% and C, 25.80%, with a high A+T bias of 60.9%. In the PCGs, ND6 had deviation ranges for AT skew (–0.303) and GC skew (0.636). The Ka/Ks values of ND1 (1.067) and ND6 (1.352) genes were >1, indicating positive selection, which might play an important role in the adaptation of L.yarkandensis to arid and hot environments. The conserved sequence block, the central conserved domain, and the extended termination-associated sequences of the control region and their features were identified and described. The phylogenetic tree based on the complete mitogenome showed that L.yarkandensis was closely related to the sympatric Lepustibetanuspamirensis. These novel datasets of L.yarkandensis can supply basic data for phylogenetic studies of Lepus spp., apart from providing essential and important resource for further genetic research and the protection of this species.

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.

Whole-genome sequencing of Tarim red deer (Cervus elaphus yarkandensis) reveals demographic history and adaptations to an arid-desert environment

Background

The initiation of desert conditions in the Tarim Basin in China since the late Miocene has led to the significant genetic structuring of local organisms. Tarim Red Deer (Cervus elaphus yarkandensis, TRD) have adapted to the harsh environmental conditions in this basin, including high solar radiation and temperature, aridity, and poor nutritional conditions. However, the underlying genetic basis of this adaptation is poorly understood.

Results

We sequenced the whole genomes of 13 TRD individuals, conducted comparative genomic analyses, and estimated demographic fluctuation. The ∂a∂i model estimated that the TRD and Tule elk (Cervus canadensis nannodes) populations diverged approximately 0.98 Mya. Analyses revealed a substantial influence of the Earth’s climate on the effective population size of TRD, associated with glacial advances and retreat, and human activities likely underlie a recent serious decline in population. A marked bottleneck may have profoundly affected the genetic diversity of TRD populations. We detected a set of candidate genes, pathways, and GO categories related to oxidative stress, water reabsorption, immune regulation, energy metabolism, eye protection, heat stress, respiratory system adaptation, prevention of high blood pressure, and DNA damage and repair that may directly or indirectly be involved in the adaptation of TRD to an arid-desert environment.

Conclusions

Our analyses highlight the role of historical global climates in the population dynamics of TRD. In light of ongoing global warming and the increasing incidence of droughts, our study offers insights into the genomic adaptations of animals, especially TRD, to extreme arid-desert environments and provides a valuable resource for future research on conservation design and biological adaptations to environmental change.

Higher expression levels of aquaporin (AQP)1 and AQP5 in the lungs of arid-desert living Lepus yarkandensis

Water transport across epithelial cells that line the airways and alveoli is a crucial component of lung physiology. Aquaporins (AQPs) facilitate water transport across the air space-capillary barrier in the distal lung. However, the roles of lung AQPs in desert animal adaptation to dry airstream environments are still unclear. A hare (Lepus yarkandensis) only lives in the Tarim Basin, and its living environment is an arid climate with rare precipitation. We studied cellular localization and expression levels of AQP1, AQP3, AQP4 and AQP5 in L. yarkandensis lungs by immunohistochemistry, quantitative real-time polymerase chain reaction and Western blot. The lung of rabbits (Oryctolagus cuniculus) that inhabit in mesic environment was similarly studied. Obtained results in two species of animals were compared to investigate whether AQPs in the lung altered expression in the animal living in arid region. AQP1 was localized to the endothelial cells in capillaries and venules surrounding terminal bronchioles and alveoli. AQP5 was localized to the ciliated columnar cells in terminal bronchioles and the alveolar type I cells in the alveolus. Quantitative real-time PCR analysis showed higher AQP1 and AQP5 mRNA levels in L. yarkandensis compared to O. cuniculus. Similar results were obtained by Western blot. These results revealed that the higher expression levels of AQP1 and AQP5 played a significant role in water transport in the lungs of arid-desert living L. yarkandensis and might accelerate water transport from capillary compartments to the airspace.

Higher Expression Levels of Aquaporin Family of Proteins in the Kidneys of Arid-Desert Living Lepus yarkandensis

Lepus yarkandensis specifically lives in arid climate with rare precipitation of Tarim Basin in western China. Aquaporins (AQPs) are a family of channel proteins that facilitate water transportation across cell membranes. Kidney AQPs play vital roles in renal tubule water permeability and maintenance of body water homeostasis. This study aimed to investigate whether kidney AQPs exhibit higher expression in arid-desert living animals. Immunohistochemistry results revealed localization of AQP1 to the capillary endothelial cells in glomerulus and epithelial cells in proximal tubule and descending thin limbs, AQP2 to the apical plasma membrane of principal cells in the cortical collecting duct (CCD), outer medullary collecting duct (OMCD), and IMCD cells in the initial inner medullary collecting duct (IMCD1) and middle IMCD (IMCD2), and AQP3 and AQP4 to the basolateral plasma membrane of principal cells and IMCD cells in CCD, OMCD, IMCD1, and IMCD2 in L. yarkandensis kidneys. Quantitative real-time PCR analysis showed higher mRNA levels of AQP1, AQP2, AQP3, and AQP4 in L. yarkandensis kidneys compared with Oryctolagus cuniculus. Similar results were obtained by western blotting. Our results suggested that higher expression levels of AQP1, AQP2, AQP3, and AQP4 in L. yarkandensis kidneys favored for drawing more water from the tubular fluid.

Phylogeography and ecological niche modeling unravel the evolutionary history of the Yarkand hare, Lepus yarkandensis (Mammalia: Leporidae), through the Quaternary

Background

The Taklimakan Desert in China is characterized by unique geological and historical dynamics and endemic flora and fauna, but the influence of historical climate oscillations on the evolutionary history of endemic animals is poorly understood. Lepus yarkandensis is an oases-dependent Near Threatened species that lives in fragmented oasis habitats in the Taklimakan Desert, China. We investigated the geological and climatic impacts on its geographical differentiation, demographic history and influence of Pleistocene glacial-interglacial cycles on the evolutionary history of L. yarkandensis. Further, studied the impact of climatic oscillation based modification on phylogeography, distribution and diversification pattern of Yarkand hare by using Cytb (1140 bp), MGF (592 bp) and SPTBN1 (619 bp) markers. Ecological niche modeling (ENM) revealed the evolutionary history of this species in response to climate change during the Quaternary. Paleodistribution modeling was used to identify putative refugia and estimate their historical distributions.

Results

Both historical demographic analyses and climatic niche modeling revealed strong effects of glacial climate changes, suggesting recurrent range contractions and expansions. The EBSP results indicated clear population expansion of L. yarkandensis since the Pleistocene. In the “early Pleistocene”, the demographic expansion continued from 0.83 MYA to the last glacial period. The ENM analysis supported a wide distribution of Lepus yarkandensis at high altitudes during the last interglacial (LIG) period. During the last glacial maximum (LGM), the suitable climate was reduced and restricted to the western part of the Taklimakan Desert.

Conclusions

Inland aridification, oasis evolution and river flow played major roles in the population differentiation and demographic history of Yarkand hares. Historically, the large, continuous oases in the Taklimakan Desert contained a viable and unique population of L. yarkandensis. The fragmented desert environment might have caused low gene flow between individuals or groups, thus leading to predominant genetic differentiation. The Pleistocene climatic cycles triggered the diversification and expansion of this species during cold and warm periods, respectively, leading to multiple colonization events within the Taklimakan Desert. These events might be due to the expansion of the Taklimakan Desert during the Middle Pleistocene. Yarkand hare previously occupied vast areas at low and intermediate altitudes in Xinjiang, Gansu, Shanxi, Henan and Shaanxi Provinces in China. The past aridification, climate change-induced oasis modifications, changes in river volumes and flow directions, and human activities all affected the population demography and phylogeography of the Yarkand hare.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1426-z) contains supplementary material, which is available to authorized users.

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.

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.

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci

BACKGROUND

Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion.

RESULTS

Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species.

CONCLUSION

This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.