Differential responses of Lasiopodomys mandarinus and Lasiopodomys brandtii to chronic hypoxia: a cross-species brain transcriptome analysis

Physiological indices and liver gene expression related to glucose supply in Brandt's vole (Lasiopodomys brandtii) exhibit species- and oxygen concentration-specific responses to hypoxia

Brandt's vole (Lasiopodomys brandtii) is a species with hypoxia tolerance, and glucose serves as the primary energy substrate under hypoxia. However, the glucose supply in Brandt's voles under hypoxia has not been studied. This study aimed to investigate characteristics in physiological indices and liver gene expression associated with glucose supply in Brandt's voles under hypoxia. Serum glucose of Brandt's voles remained stable under 10% O2, increased under 7.5% O2, and decreased under 5% O2. Serum lactate increased under 10% O2, decreased under 7.5% O2, increased at 6 h and decreased at 12 h under 5% O2. Liver glycogen increased under 10% O2, remained constant under 7.5% O2, and reduced under 5% O2. Pepck and G6pase expression associated with gluconeogenesis decreased under 10% O2, while Pepck expression decreased and G6pase expression increased under 7.5% and 5% O2. Regarding genes related to glycogen metabolism, Gys expression decreased at all oxygen concentrations, Phk expression increased under 5% O2, and Gp expression increased under 7.5% and 5% O2. The alterations in glucose, lactate, liver glycogen, and gene expression related to glycogenolysis in Kunming mice (Mus musculus, control species) are similar to discovery of Brandt's voles under 7.5% O2, but gene expression involved in gluconeogenesis and glycogen synthesis increased. The findings suggest that Brandt's voles are more tolerant to hypoxia than Kunming mice, and their physiological indices and liver gene expression related to glucose supply exhibit species- and oxygen concentration-specific responses to hypoxia. This research offers novel insights for studying hypoxia tolerance of Brandt's voles.

Weak gene-gene interaction facilitates the evolution of gene expression plasticity

BACKGROUND

Individual organisms may exhibit phenotypic plasticity when they acclimate to different conditions. Such plastic responses may facilitate or constrain the adaptation of their descendant populations to new environments, complicating their evolutionary trajectories beyond the genetic blueprint. Intriguingly, phenotypic plasticity itself can evolve in terms of its direction and magnitude during adaptation. However, we know little about what determines the evolution of phenotypic plasticity, including gene expression plasticity. Recent laboratory-based studies suggest dominance of reversing gene expression plasticity-plastic responses that move the levels of gene expression away from the new optima. Nevertheless, evidence from natural populations is still limited.

RESULTS

Here, we studied gene expression plasticity and its evolution in the montane and lowland populations of an elevationally widespread songbird-the Rufous-capped Babbler (Cyanoderma ruficeps)-with reciprocal transplant experiments and transcriptomic analyses; we set common gardens at altitudes close to these populations' native ranges. We confirmed the prevalence of reversing plasticity in genes associated with altitudinal adaptation. Interestingly, we found a positive relationship between magnitude and degree of evolution in gene expression plasticity, which was pertinent to not only adaptation-associated genes but also the whole transcriptomes from multiple tissues. Furthermore, we revealed that genes with weaker expressional interactions with other genes tended to exhibit stronger plasticity and higher degree of plasticity evolution, which explains the positive magnitude-evolution relationship.

CONCLUSIONS

Our experimental evidence demonstrates that species may initiate their adaptation to new habitats with genes exhibiting strong expression plasticity. We also highlight the role of expression interdependence among genes in regulating the magnitude and evolution of expression plasticity. This study illuminates how the evolution of phenotypic plasticity in gene expression facilitates the adaptation of species to challenging environments in nature.

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.

Antioxidant response to severe hypoxia in Brandt's vole Lasiopodomys brandtii

The antioxidant defense system is essential for animals to cope with homeostasis disruption and overcome oxidative stress caused by adverse environmental conditions such as hypoxia. However, our understanding of how this system works in subterranean rodents remains limited. In this study, Brandt's vole Lasiopodomys brandtii was exposed to normoxia (21% O₂ ) or hypoxia (mild or severe hypoxia: 10% or 5% O₂ ) for 6 h. Changes in key enzymes of the classic enzymatic antioxidant system at both mRNA and enzyme activity levels, and tissue antioxidant levels of the low-molecular-weight antioxidant system were determined in brain, liver, and kidney. Transcript levels of the upstream regulator NF-E2-related factor 2 (Nrf2) were also measured. We found that the mRNA expression of Nrf2 and its downstream antioxidant enzyme genes in L. brandtii were relatively conserved in response to hypoxia in most tissues and genes tested, except in the liver. Hepatic Nrf2, Cu/Zn SOD, GPx1, and GPx3 levels were significantly upregulated in response to mild hypoxia, whereas Mn SOD level decreased significantly in severe hypoxia. Unmatched with changes at the RNA level, constitutively high and relatively stable antioxidant enzyme activities were maintained throughout. For the low-molecular-weight antioxidant system, an abrupt increase of cerebral ascorbic acid (AA) levels in hypoxia indicated a tissue-specific antioxidant response. Although hypoxia did not cause significant oxidative damage in most tissues tested, the significant decrease in antioxidant enzyme activities (GPX and GR) and increase in lipid peroxidation in the kidney suggest that prolonged hypoxia may pose a critical threat to this species.

Diverse energy metabolism patterns in females in Neodon fuscus, Lasiopodomys brandtii, and Mus musculus revealed by comparative transcriptomics under hypoxic conditions

The effects of global warming and anthropogenic disturbance force animals to migrate from lower to higher elevations to find suitable new habitats. As such migrations increase hypoxic stress on the animals, it is important to understand how plateau- and plain-dwelling animals respond to low-oxygen environments. We used comparative transcriptomics to explore the response of Neodon fuscus, Lasiopodomys brandtii, and Mus musculus skeletal muscle tissues to hypoxic conditions. Results indicate that these species have adopted different oxygen transport and energy metabolism strategies for dealing with a hypoxic environment. N. fuscus promotes oxygen transport by increasing hemoglobin synthesis and reduces the risk of thrombosis through cooperative regulation of genes, including Fga, Fgb, Alb, and Ttr; genes such as Acs16, Gpat4, and Ndufb7 are involved in regulating lipid synthesis, fatty acid β-oxidation, hemoglobin synthesis, and electron-linked transmission, thereby maintaining a normal energy supply in hypoxic conditions. In contrast, the oxygen-carrying capacity and angiogenesis of red blood cells in L. brandtii are promoted by genes in the CYP and COL families; this species maintains its bodily energy supply by enhancing the pentose phosphate pathway and mitochondrial fatty acid synthesis pathway. However, under hypoxia, M. musculus cannot effectively transport additional oxygen; thus, its cell cycle, proliferation, and migration are somewhat affected. Given its lack of hypoxic tolerance experience, M. musculus also shows significantly reduced oxidative phosphorylation levels under hypoxic conditions. Our results suggest that the glucose capacity of M. musculus skeletal muscle does not provide sufficient energy during hypoxia; thus, we hypothesize that it supplements its bodily energy by synthesizing ketone bodies. For the first time, we describe the energy metabolism pathways of N. fuscus and L. brandtii skeletal muscle tissues under hypoxic conditions. Our findings, therefore, improve our understanding of how vertebrates thrive in high altitude and plain habitats when faced with hypoxic conditions.

Phylogeny and evolution of Lasiopodomys in subfamily Arvivolinae based on mitochondrial genomics

The species of Lasiopodomys Lataste 1887 with their related genera remains undetermined owing to inconsistent morphological characteristics and molecular phylogeny. To investigate the phylogenetic relationship and speciation among species of the genus Lasiopodomys, we sequenced and annotated the whole mitochondrial genomes of three individual species, namely Lasiopodomys brandtii Radde 1861, L. mandarinus Milne-Edwards 1871, and Neodon (Lasiopodomys) fuscus Büchner 1889. The nucleotide sequences of the circular mitogenomes were identical for each individual species of L. brandtii, L. mandarinus, and N. fuscus. Each species contained 13 protein-coding genes (PCGs), 22 transfer RNAs, and 2 ribosomal RNAs, with mitochondrial genome lengths of 16,557 bp, 16,562 bp, and 16,324 bp, respectively. The mitogenomes and PCGs showed positive AT skew and negative GC skew. Mitogenomic phylogenetic analyses suggested that L. brandtii, L. mandarinus, and L. gregalis Pallas 1779 belong to the genus Lasiopodomys, whereas N. fuscus belongs to the genus Neodon grouped with N. irene. Lasiopodomys showed the closest relationship with Microtus fortis Büchner 1889 and M. kikuchii Kuroda 1920, which are considered as the paraphyletic species of genera Microtus. T_MRCA and niche model analysis revealed that Lasiopodomys may have first appeared during the early Pleistocene epoch. Further, L. gregalis separated from others over 1.53 million years ago (Ma) and then diverged into L. brandtii and L. mandarinus 0.76 Ma. The relative contribution of climatic fluctuations to speciation and selection in this group requires further research.

Regulation of HIF-1α and p53 in stress responses in the subterranean rodents Lasiopodomys mandarinus and Lasiopodomys brandtii (Rodentia: Cricetidae)

The response mechanism and interaction patterns of HIF-1α and p53 in animals in an hypoxic environment are crucial for their hypoxic tolerance and adaptation. Many studies have shown that underground rodents have better hypoxic adaptation characteristics. However, the mechanism by which HIF-1α and p53 in underground rodents respond to hypoxic environments compared with in ground rodents remains unclear. Further, whether a synergy between HIF-1α and p53 enables animals tolerate extremely hypoxic environments is unclear. We studied HIF-1α and p53 expression in the brain tissue and cell apoptosis in the hippocampal CA1 region during 6 hours of acute hypoxia (5% oxygen) in Lasiopodomys mandarinus (Milne-Edwards, 1871) and Lasiopodomys brandtii (Radde, 1861), two closely related small rodents with different life characteristics (underground and aboveground, respectively), using a comparative biology method to determine the mechanisms underlying their adaptation to this environment. Our results indicate that HIF-1α and p53 expression is more rapid in L. mandarinus than in L. brandtii under acute hypoxic environments, resulting in a significant synergistic effect in L. mandarinus. Correlation analysis revealed that HIF-1α expression and the apoptotic index of the hippocampal CA1 regions of the brain tissues of L. mandarinus and L. brandtii, both under hypoxia, were significantly negatively and positively correlated, respectively. Long-term existence in underground burrow systems could enable better adaptation to hypoxia in L. mandarinus than in L. brandtii. We speculate that L. mandarinus can quickly eliminate resulting damage via the synergistic effect of p53 and HIF-1α in response to acute hypoxic environments, helping the organism quickly return to a normal state after the stress.

Adaptive evolution of β-globin gene in subterranean in South America octodontid rodents

The convergent evolution of subterranean rodents is an excellent model to study how natural selection operates and the genetic bases of these adaptations, but the study on the different taxa has been very uneven and still insufficient. In the octodontoid caviomorph rodent superfamily there are two independent lineages where they have recently evolved into totally underground lifestyles: the genera Ctenomys (tuco-tucos) and Spalacopus (coruro). The underground habitat is characterized by an hypoxic and hypercapnic atmosphere, thus gas exchange is one of the most important challenges for these animals. The invasion of the underground niche could have modified the selective regimes of proteins involved in the respiration and transport of O₂ of these rodents, positively selecting mutations of higher affinity for O₂. Here we examine the sequence variation in the beta globin gene in these two lineages, within a robust phylogenetic context. Using different approaches (classical and Bayesian maximum likelihood (PAML/Datamonkey) and alternatives methods (TreeSAAP)) we found at least three sites with evidence of positive selection in underground lineages, especially the basal branch that leads to the Octodontidae family and the branch that leads to the coruro, suggesting some adaptive changes to the underground life. We also found a convergence with another underground rodent, which cannot be identified by the above methods.

Genome-wide adaptive evolution to underground stresses in subterranean mammals: Hypoxia adaption, immunity promotion, and sensory specialization

Life underground has provided remarkable examples of adaptive evolution in subterranean mammals; however, genome-wide adaptive evolution to underground stresses still needs further research. There are approximately 250 species of subterranean mammals across three suborders and six families. These species not only inhabit hypoxic and dark burrows but also exhibit evolved adaptation to hypoxia, cancer resistance, and specialized sensory systems, making them an excellent model of evolution. The adaptive evolution of subterranean mammals has attracted great attention and needs further study. In the present study, phylogenetic analysis of 5,853 single-copy orthologous gene families of five subterranean mammals (Nannospalax galili, Heterocephalus glaber, Fukomys damarensis, Condylura cristata, and Chrysochloris asiatica) showed that they formed fou distinct clusters. This result is consistent with the traditional systematics of these species. Furthermore, comparison of the high-quality genomes of these five subterranean mammalian species led to the identification of the genomic signatures of adaptive evolution. Our results show that the five subterranean mammalian did not share positively selected genes but had similar functional enrichment categories, including hypoxia tolerance, immunity promotion, and sensory specialization, which adapted to the environment of underground stresses. Moreover, variations in soil hardness, climate, and lifestyles have resulted in different molecular mechanisms of adaptation to the hypoxic environment and different degrees of visual degradation. These results provide insights into the genome-wide adaptive evolution to underground stresses in subterranean mammals, with special focus on the characteristics of hypoxia adaption, immunity promotion, and sensory specialization response to the life underground.

Transcriptome analysis of the response provided by Lasiopodomys mandarinus to severe hypoxia includes enhancing DNA repair and damage prevention

Background

Severe hypoxia induces a series of stress responses in mammals; however, subterranean rodents have evolved several adaptation mechanisms of energy metabolisms and O₂ utilization for hypoxia. Mammalian brains show extreme aerobic metabolism. Following hypoxia exposure, mammals usually experience irreversible brain damage and can even develop serious diseases, such as hypoxic ischemic encephalopathy and brain edema. To investigate mechanisms underlying the responses of subterranean rodents to severe hypoxia, we performed a cross-species brain transcriptomic analysis using RNA sequencing and identified differentially expressed genes (DEGs) between the subterranean rodent Lasiopodomys mandarinus and its closely related aboveground species L. brandtii under severe hypoxia (5.0% O₂, 6 h) and normoxia (20.9% O₂, 6 h).

Results

We obtained 361 million clean reads, including 69,611 unigenes in L. mandarinus and 69,360 in L. brandtii. We identified 359 and 515 DEGs by comparing the hypoxic and normoxia groups of L. mandarinus and L. brandtii, respectively. Gene Ontology (GO) analysis showed that upregulated DEGs in both species displayed similar terms in response to severe hypoxia; the main difference is that GO terms of L. brandtii were enriched in the immune system. However, in the downregulated DEGs, GO terms of L. mandarinus were enriched in cell proliferation and protein transport and those of L. brandtii were enriched in nuclease and hydrolase activities, particularly in terms of developmental functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that upregulated DEGs in L. mandarinus were associated with DNA repair and damage prevention as well as angiogenesis and metastasis inhibition, whereas downregulated DEGs were associated with neuronal synaptic transmission and tumor-associated metabolic pathways. In L. brandtii, upregulated KEGG pathways were enriched in the immune, endocrine, and cardiovascular systems and particularly in cancer-related pathways, whereas downregulated DEGs were associated with environmental information processing and misregulation in cancers.

Conclusions

L. mandarinus has evolved hypoxia adaptation by enhancing DNA repair, damage prevention, and augmenting sensing, whereas L. brandtii showed a higher risk of tumorigenesis and promoted innate immunity toward severe hypoxia. These results reveal the hypoxic mechanisms of L. mandarinus to severe hypoxia, which may provide research clues for hypoxic diseases.

The complete mitochondrial genome of Lasiopodomys brandtii

In this study, the complete mitochondrial (mt) genome sequence of Lasiopodomys brandtii was determined using Illumina NovaSeq platform. The assembled genome was 16,557 bp in length and included 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. The total nucleotide composition frequencies present clearly the A-T skew (59.5%), which mostly in D-loop and PCGs regions. Whole mt genome phylogenetic analysis revealed a closely related among Lasiopodomys, Proedromys, and Microtus with high support. It would provide further evolutionary research for the subfamily Arvicolinae.