Impact of Orogeny and Environmental Change on Genetic Divergence and Demographic History of Dipus sagitta (Dipodoidea, Dipodinae) since the Pliocene in Inland East Asia

Whole-genome sequencing reveals adaptations of hairy-footed jerboas (Dipus, Dipodidae) to diverse desert environments

BACKGROUND

Environmental conditions vary among deserts across the world, spanning from hyper-arid to high-elevation deserts. However, prior genomic studies on desert adaptation have focused on desert and non-desert comparisons overlooking the complexity of conditions within deserts. Focusing on the adaptation mechanisms to diverse desert environments will advance our understanding of how species adapt to extreme desert environments. The hairy-footed jerboas are well adapted to diverse desert environments, inhabiting high-altitude arid regions, hyper-arid deserts, and semi-deserts, but the genetic basis of their adaptation to different deserts remains unknown.

RESULTS

Here, we sequenced the whole genome of 83 hairy-footed jerboas from distinct desert zones in China to assess how they responded under contrasting conditions. Population genomics analyses reveal the existence of three species in hairy-footed jerboas distributed in China: Dipus deasyi, Dipus sagitta, and Dipus sowerbyi. Analyses of selection between high-altitude desert (elevation ≥ 3000m) and low-altitude desert (< 500m) populations identified two strongly selected genes, ATR and HIF1AN, associated with intense UV radiation and hypoxia in high-altitude environments. A number of candidate genes involved in energy and water homeostasis were detected in the comparative genomic analyses of hyper-arid desert (average annual precipitation < 70mm) and arid desert (< 200mm) populations versus semi-desert (> 360mm) populations. Hyper-arid desert animals also exhibited stronger adaptive selection in energy homeostasis, suggesting water and resource scarcity may be the main drivers of desert adaptation in hairy-footed jerboas.

CONCLUSIONS

Our study challenges the view of deserts as homogeneous environments and shows that distinct genomic adaptations can be found among desert animals depending on their habitats.

Cranial differences in three-toed jerboas (Dipodinae, Dipodidae, Rodentia) according to recent taxonomic revisions

Recent phylogenetic studies amended the taxonomy of three-toed jerboas (subfamily Dipodinae), including raising subspecies to full species. Here, we use geometric morphometrics to compare scaled-shape differences in dipodine crania while considering their revised taxonomy. We sampled Dipus deasyi, D. sagitta halli, D. s. sowerbyi, Jaculus blanfordi blanfordi, J. hirtipes, J. jaculus, J. loftusi, J. orientalis gerboa, J. o. mauritanicus, and Stylodipus andrewsi. Crania were not sexually dimorphic. Common allometry explained some of the shape variation, for example, reduced braincases in larger specimens. Most operational taxonomic unit pairs differed in both size and shape. Dipus and Stylodipus clustered together based on their cranial shape. Jaculus differed from the aforementioned genera by its larger tympanic bulla, broader braincase, larger infraorbital foramen, along with reduced molars and rostra. Jaculus orientalis differed from other Jaculus by its broader face versus reduced cranial vault. Jaculus blanfordi (subgenus Haltomys) resembles members of the subgenus Jaculus more than its consubgener (J. orientalis). Jaculus loftusi, previously considered a synonym of J. jaculus, clearly differed from the latter by its shorter rostrum, smaller infraorbital foramen, and more caudolaterally expanded tympanic bulla. Jaculus hirtipes, another recent synonym of J. jaculus, resembled J. blanfordi more in scaled cranial shape than it did J. jaculus. Dipus sagitta halli and D. s. sowerbyi were indistinguishable, but they clearly differed from D. deasyi (recently raised to full species) with the latter having a larger molar row, more inflated tympanic bulla, and shorter, slenderer rostrum. Ecological explanations for detected cranial shape differences are considered, including diet and habitat (particularly substrate).

Geography and past climate changes have shaped the evolution of a widespread lizard in arid Central Asia

The complex orogenic history and structure of Central Asia, coupled with Pleistocene glacial cycles have generated its stepwise aridification. Such events would have significantly influenced the evolution of many mid-latitude species in arid Central Asia (ACA). In this study, we employed two mitochondrial genes (CO1 and ND2) and genome-wide SNPs, coupled with ecological niche modeling, to investigate the lineage diversification and historical demography within a widespread lizard Phrynocepahlus helioscopus, and their associations with geography and past climate change. We obtained themtDNA dataset for 300 individuals from 96 localities within the known range of the lizard, among which 51 individuals from 27 localities were selected for generating the SNP dataset via genotyping-by-sequencing approach. Phylogenetic analyses of the concatenated mtDNA dataset revealed eight geographically correlated lineages that diverged by 4.21-10.41% for the CO1 gene, which were estimated to have coalesced ∼4.47 million years ago. However, we observed mito-nuclear discordance pattern regarding the population of Clade V (P. helioscopus sergeevi) from Tajikistan. Ancestral area estimations suggested that P. helioscopus originated from the Fergana Valley and then dispersed into the adjacent areas in ACA along with a history of multiple allopatric divergence processes, suggesting that Fergana may have been the cradle of diversification of P. helioscopus. The intensification of aridification across Central Asia during the Late Pliocene may have facilitated the rapid radiation of this arid-adapted lizard throughout this vast territory. Subsequently, the geological events (e.g., uplift of the Hissar-Alay, transgressions of the Caspian Sea) and geographic barriers (e.g., Amu Darya River, Zerarshan River) during the Pleistocene triggered the progressive diversification of P. helioscopus. Interestingly, Clade VIII (P. helioscopus varius) experienced rapid population growth coupled with range expansion while Clade IV (P. helioscopus cameranoi) underwent drastic population expansion associated with range contraction during the Last Glacial Maximum. In Clade IV, but not in Clade VIII, environmental turnover contributes more to mitochondrial genetic distinctiveness than geographic distance. Overall, the SNP dataset demonstrates that geographic distance plays a greater role than environmental distance. Both the mtDNA dataset and the SNP dataset suggest local-scale genetic differentiation in Clade IV and Clade VIII, revealing potential geographic barriers in the Ili River Valley and the Junggar Basin, respectively. Twenty-seven outlier SNPs associated with environmental factors (precipitation and temperature) were identified, which supports the signature of local adaptation to the arid desert environment. Finally, our finding suggests taxonomic implications, such as support for full species status for P. saidalievi (Clade II) and P. meridionalis (Clade I). Future analyses based on further evidence and increased taxon and geographic sampling should be carried out to corroborate our findings.

Phylogeny of the Dipus sagitta Species Complex by Nuclear Gene Sequences

The northern three-toed jerboa Dipus sagitta had long been considered to be a single polytypic species. High genetic diversity of D. sagitta was earlier revealed on the basis of several mitochondrial and nuclear genes, and several separate species were hypothesized to occur within the taxon. However, the relationships between phylogenetic lineages have not been established because of the small sample size of nuclear genes. In the present work, a far larger set of nuclear DNA loci was used, and thus, a higher resolution of the phylogenetic tree was achieved for ten D. sagitta forms. The structure revealed for the species mainly confirmed the topology and relationships of the mtDNA lineages. Yet the mitochondrial and nuclear phylogenies were not completely consistent. Some of the D. sagitta genetic lineages were therefore assumed to be a product of reticular evolutionary processes. The taxon was concluded to be the diverse species complex D. sagitta sensu lato, in which long-diverged lineages are not always reproductively isolated.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

Phylogeny, taxonomic reassessment and ‘ecomorph’ relationship of the Orientallactaga sibirica complex (Rodentia: Dipodidae: Allactaginae)

The ecological gradient–morphological variation (‘ecomorph’) relationship has long interested ecologists and evolutionary biologists, but it is applied far less frequently than genetic differentiation in cryptic species detection and species identification. With integrative methods, we revise taxonomic uncertainties in the Orientallactaga sibirica complex (OSC), with 298 sequence specimens and 469 voucher specimens from 138 localities covering nearly the entire distribution of the OSC. Phylogenetic relationships are assessed by Bayesian inference and maximum likelihood using two mitochondrial and nine nuclear genes. We use species-delimitation approaches to divide and validate the ‘candidate species’. We evaluate correlations between ecological divergence and phylogenetic splits, and visualize geographical patterns of morphological variation. The OSC is divided into four phylogenetic groups, the Ognevi, Altay, Bogda and Sibirica groups, and the OSC exhibits a significant ecomorph relationship and ecological divergence pattern. Morphological variations not only follow the general regularity under a large gradient of ecological factors, but are also closely related to the local environment/habitat. We suggest considering the comprehensive ecomorph relationship to identify species. Molecular analyses reveal that the OSC more easily forms deeply divergent lineages in the foothills and this differentiation depth may be related to mountain system size.

Divergent selection along elevational gradients promotes genetic and phenotypic disparities among small mammal populations

Species distributed along mountain slopes, facing contrasting habitats in short geographic scale, are of particular interest to test how ecologically based divergent selection promotes phenotypic and genetic disparities as well as to assess isolation-by-environment mechanisms. Here, we conduct the first broad comparative study of phenotypic variation along elevational gradients, integrating a large array of ecological predictors and disentangling population genetic driver processes. The skull form of nine ecologically distinct species distributed over a large altitudinal range (100-4200 m) was compared to assess whether phenotypic divergence is a common phenomenon in small mammals and whether it shows parallel patterns. We also investigated the relative contribution of biotic (competition and predation) and abiotic parameters on phenotypic divergence via mixed models. Finally, we assessed the population genetic structure of a rodent species (Niviventer confucianus) via analysis of molecular variance and FST along three mountain slopes and tested the isolation-by-environment hypothesis using Mantel test and redundancy analysis. We found a consistent phenotypic divergence and marked genetic structure along elevational gradients; however, the species showed mixed patterns of size and skull shape trends across mountain zones. Individuals living at lower altitudes differed greatly in both phenotype and genotype from those living at high elevations, while middle-elevation individuals showed more intermediate forms. The ecological parameters associated with phenotypic divergence along elevation gradients are partly related to species' ecological and evolutionary constraints. Fossorial and solitary animals are mainly affected by climatic factors, while terrestrial and more gregarious species are influenced by biotic and abiotic parameters. A novel finding of our study is that predator richness emerged as an important factor associated with the intraspecific diversification of the mammalian skull along elevational gradients, a previously overlooked parameter. Population genetic structure was mainly driven by environmental heterogeneity along mountain slopes, with no or a week spatial effect, fitting the isolation-by-environment scenario. Our study highlights the strong and multifaceted effects of heterogeneous steep habitats and ecologically based divergent selective forces in small mammal populations.

Disjunct distribution and distinct intraspecific diversification of Eothenomys melanogaster in South China

BACKGROUND

South China encompasses complex and diverse landforms, giving rise to high biological diversity and endemism from the Hengduan Mountains to Taiwan Island. Many species are widely distributed across South China with similar disjunct distribution patterns. To explore the causes of these disjunct distribution patterns and their genetic consequences, we investigated the endemic species Père David's Chinese Vole (Eothenomys melanogaster) by integrating geological and ecological factors. We analysed the genetic structure and divergence time of E. melanogaster based on fast-evolving mitochondrial and nuclear markers using Bayesian trees and coalescent species tree approaches. Historical scenarios of distribution range and demography were reconstructed based on spatial interpolations of genetic diversity and distance, extended Bayesian skyline plots, phylogeographic diffusion analysis, and ecological niche modelling (ENM) during different periods. We also assessed the relationships between geographical distance/ecological vicariance and genetic distance (isolation by distance, IBD; isolation by environment, IBE).

RESULTS

The genetic analysis revealed three deeply divergent clades-Southeast, Southwest and Central clades, centred on the Wuyi Mountains, the Yunnan-Guizhou Plateau (YGP) and the mountains around the Sichuan Basin, respectively-that have mostly developed since the Pleistocene. IBD played an important role in early divergence, and geological events (sedimentation of plains and linking of palaeo-rivers) and IBE further reinforced genetic differentiation. ENM shows the importance of suitable habitats and elevations.

CONCLUSIONS

Our results suggest that the primary cause of the disjunct distribution in E. melanogaster is the high dependence on middle-high-altitude habitat in the current period. Mountains in the occurence range have served as "sky islands" for E. melanogaster and hindered gene flow. Pleistocene climatic cycles facilitated genetic admixture in cold periods and genetic diversification in warm periods for inland clades. During cold periods, these cycles led to multiple colonization events between the mainland and Taiwan and erased genetic differentiation.