Adaptive Transcriptome Profiling of Subterranean Zokor, Myospalax baileyi, to High- Altitude Stresses in Tibet

Genomic structural variation is associated with hypoxia adaptation in high-altitude zokors

Zokors, an Asiatic group of subterranean rodents, originated in lowlands and colonized high-elevational zones following the uplift of the Qinghai-Tibet plateau about 3.6 million years ago. Zokors live at high elevation in subterranean burrows and experience hypobaric hypoxia, including both hypoxia (low oxygen concentration) and hypercapnia (elevated partial pressure of CO2). Here we report a genomic analysis of six zokor species (genus Eospalax) with different elevational ranges to identify structural variants (deletions and inversions) that may have contributed to high-elevation adaptation. Based on an assembly of a chromosome-level genome of the high-elevation species, Eospalax baileyi, we identified 18 large inversions that distinguished this species from congeners native to lower elevations. Small-scale structural variants in the introns of EGLN1, HIF1A, HSF1 and SFTPD of E. baileyi were associated with the upregulated expression of those genes. A rearrangement on chromosome 1 was associated with altered chromatin accessibility, leading to modified gene expression profiles of key genes involved in the physiological response to hypoxia. Multigene families that underwent copy-number expansions in E. baileyi were enriched for autophagy, HIF1 signalling and immune response. E. baileyi show a significantly larger lung mass than those of other Eospalax species. These findings highlight the key role of structural variants underlying hypoxia adaptation of high-elevation species in Eospalax.

The hypoxia adaptation of small mammals to plateau and underground burrow conditions

Oxygen is one of the important substances for the survival of most life systems on the earth, and plateau and underground burrow systems are two typical hypoxic environments. Small mammals living in hypoxic environments have evolved different adaptation strategies, which include increased oxygen delivery, metabolic regulation of physiological responses and other physiological responses that change tissue oxygen utilization. Multi-omics predictions have also shown that these animals have evolved different adaptations to extreme environments. In particular, vascular endothelial growth factor (VEGF) and erythropoietin (EPO), which have specific functions in the control of O2 delivery, have evolved adaptively in small mammals in hypoxic environments. Naked mole-rats and blind mole-rats are typical hypoxic model animals as they have some resistance to cancer. This review primarily summarizes the main living environment of hypoxia tolerant small mammals, as well as the changes of phenotype, physiochemical characteristics and gene expression mode of their long-term living in hypoxia environment.

Phenotypic and genomic adaptations to the extremely high elevation in plateau zokor (Myospalax baileyi)

The evolutionary outcomes of high elevation adaptation have been extensively described. However, whether widely distributed high elevation endemic animals adopt uniform mechanisms during adaptation to different elevational environments remains unknown, especially with respect to extreme high elevation environments. To explore this, we analysed the phenotypic and genomic data of seven populations of plateau zokor (Myospalax baileyi) along elevations ranging from 2,700 to 4,300 m. Based on whole-genome sequencing data and demographic reconstruction of the evolutionary history, we show that two populations of plateau zokor living at elevations exceeding 3,700 m diverged from other populations nearly 10,000 years ago. Further, phenotypic comparisons reveal stress-dependent adaptation, as two populations living at elevations exceeding 3,700 m have elevated ratios of heart mass to body mass relative to other populations, and the highest population (4,300 m) displays alterations in erythrocytes. Correspondingly, genomic analysis of selective sweeps indicates that positive selection might contribute to the observed phenotypic alterations in these two extremely high elevation populations, with the adaptive cardiovascular phenotypes of both populations possibly evolving under the functional constrains of their common ancestral population. Taken together, phenotypic and genomic evidence demonstrates that heterogeneous stressors impact adaptations to extreme elevations and reveals stress-dependent and genetically constrained adaptation to hypoxia, collectively providing new insights into the high elevation adaptation.

Transcriptome analysis of the liver of Eospalax fontanierii under hypoxia

Hypoxia can induce cell damage, inflammation, carcinogenesis, and inhibit liver regeneration in non-adapted species. Because of their excellent hypoxia adaptation features, subterranean rodents have been widely studied to clarify the mechanism of hypoxia adaptation. Eospalax fontanierii, which is a subterranean rodent found in China, can survive for more than 10 h under 4% O₂ without observable injury, while Sprague-Dawley rats can survive for less than 6 h under the same conditions. To explore the potential mechanism of hypoxia responses in E. fontanierii, we performed RNA-seq analysis of the liver in E. fontanierii exposed to different oxygen levels (6.5% 6h, 10.5% 44h, and 21%). Based on the bioinformatics analysis, 39,439 unigenes were assembled, and 56.78% unigenes were annotated using public databases (Nr, GO, Swiss-Prot, KEGG, and Pfam). In total, 725 differentially expressed genes (DEGs) were identified in the response to hypoxia; six with important functions were validated by qPCR. Those DEGs were mainly involved in processes related to lipid metabolism, steroid catabolism, glycolysis/gluconeogenesis, and the AMPK and PPAR signaling pathway. By analyzing the expression patterns of important genes related to energy associated metabolism under hypoxia, we found that fatty acid oxidation and gluconeogenesis were increased, while protein synthesis and fatty acid synthesis were decreased. Furthermore, the upregulated expression of specific genes with anti-apoptosis or anti-oxidation functions under hypoxia may contribute to the mechanism by which E. fontanierii tolerates hypoxia. Our results provide an understanding of the response to hypoxia in E. fontanierii, and have potential value for biomedical studies.

Gut Microbiome Changes in Captive Plateau Zokors (Eospalax baileyi)

Wild-caught animals must cope with drastic lifestyle and dietary changes after being induced to captivity. How the gut microbiome structure of these animals will change in response receives increasing attention. The plateau zokor (Eospalax baileyi), a typic subterranean rodent endemic to the Qinghai-Tibet plateau, spends almost the whole life underground and is well adapted to the environmental pressures of both plateau and underground. However, how the gut microbiome of the plateau zokor will change in response to captivity has not been reported to date. This study compared the microbial community structure and functions of 22 plateau zokors before (the WS group) and after being kept in captivity for 15 days (the LS group, fed on carrots) using the 16S rRNA gene via high-throughput sequencing technology. The results showed that the LS group retained 973 of the 977 operational taxonomic units (OTUs) in the WS group, and no new OTUs were found in the LS group. The dominant bacterial phyla were Bacteroides and Firmicutes in both groups. In alpha diversity analysis, the Shannon, Sobs, and ACE indexes of the LS group were significantly lower than those of the WS group. A remarkable difference (P < 0.01) between groups was also detected in beta diversity analysis. The UPGMA clustering, NMDS, PCoA, and Anosim results all showed that the intergroup difference was significantly greater than the intragroup difference. And compared with the WS group, the intragroup difference of the gut microbiota in the LS group was much larger, which failed to support the assumption that similar diets should drive convergence of gut microbial communities. PICRUSt revealed that although some functional categories displayed significant differences between groups, the relative abundances of these categories were very close in both groups. Based on all the results, we conclude that as plateau zokors enter captivity for a short time, although the relative abundances of different gut microbiota categories shifted significantly, they can maintain almost all the OTUs and the functions of the gut microbiota in the wild. So, the use of wild-caught plateau zokors in gut microbial studies is acceptable if the time in captivity is short.

Cross-Species Insights Into Genomic Adaptations to Hypoxia

Over millions of years, vertebrate species populated vast environments spanning the globe. Among the most challenging habitats encountered were those with limited availability of oxygen, yet many animal and human populations inhabit and perform life cycle functions and/or daily activities in varying degrees of hypoxia today. Of particular interest are species that inhabit high-altitude niches, which experience chronic hypobaric hypoxia throughout their lives. Physiological and molecular aspects of adaptation to hypoxia have long been the focus of high-altitude populations and, within the past decade, genomic information has become increasingly accessible. These data provide an opportunity to search for common genetic signatures of selection across uniquely informative populations and thereby augment our understanding of the mechanisms underlying adaptations to hypoxia. In this review, we synthesize the available genomic findings across hypoxia-tolerant species to provide a comprehensive view of putatively hypoxia-adaptive genes and pathways. In many cases, adaptive signatures across species converge on the same genetic pathways or on genes themselves [i.e., the hypoxia inducible factor (HIF) pathway). However, specific variants thought to underlie function are distinct between species and populations, and, in most cases, the precise functional role of these genomic differences remains unknown. Efforts to standardize these findings and explore relationships between genotype and phenotype will provide important clues into the evolutionary and mechanistic bases of physiological adaptations to environmental hypoxia.

Plant diversity is closely related to the density of zokor mounds in three alpine rangelands on the Tibetan Plateau

Background

Plateau zokor (Myospalax baileyi) is a subterranean rodent endemic to the Tibetan Plateau. This species has been generally viewed as a pest in China due to the competition for food with livestock and also causing soil erosion. As a result, plateau zokor has been the target of widespread poisoning or trapping campaigns designed to control or eliminate it since 1970s. But there is little research on the effect of plateau zokor on plant diversity in alpine rangelands. Therefore, objectively evaluating the positive effects of the plateau zokors disturbance on their living environment and plant communities is of great significance to understand the function of plateau zokor in alpine ecosystem.

Methods

Here, we selected three rangelands (alpine meadow, alpine steppe and alpine shrub meadow) in which plateau zokors are typically distributed on the Tibetan Plateau, and five zokor mound density gradients were selected in each rangeland type to study the effects of the mounds on soil moisture and temperature related to plant species diversity.

Results

The results showed that, with the mound density increasing, the soil temperature decreased significantly in all three rangeland types, and the soil moisture significantly increased in all three rangeland types. In the alpine meadow, both the plant diversity and cumulative species richness increased significantly with increasing mound density. The increase in broad-leaved forbs is the main reason for the increase of plant diversity in the alpine meadow disturbed by zokor mounds. In the alpine steppe, the plant diversity decreased significantly with increasing mound density, while the cumulative species richness initially decreased and then increased. In the alpine shrub meadow, the plant diversity first increased and then decreased with increasing mound density as did the cumulative species richness. In conclusion, plateau zokor mounds dominated the distribution of soil moisture and temperature and significantly affected plant diversity in these three rangelands on Tibetan Plateau; the results further deepen our understanding toward a co-evolved process.