AbstractIn the subterranean rodent (Nanno)spalax galili, evolutionary adaptation to hypoxia is correlated with longevity and tumor resistance. Adapted gene-regulatory networks of Spalax might pinpoint strategies to maintain health in humans. Comparing liver, kidney and spleen transcriptome data from Spalax and rat at hypoxia and normoxia, we identified differentially expressed gene pathways common to multiple organs in both species. Body-wide interspecies differences affected processes like cell death, antioxidant defense, DNA repair, energy metabolism, immune response and angiogenesis, which may play a crucial role in Spalax’s adaptation to environmental hypoxia. In all organs, transcription of genes for genome stability maintenance and DNA repair was elevated in Spalax versus rat, accompanied by lower expression of aerobic energy metabolism and proinflammatory genes. These transcriptomic changes might account for the extraordinary lifespan of Spalax and its cancer resistance. The identified gene networks present candidates for further investigating the molecular basis underlying the complex Spalax phenotype.
IntroductionUnderstanding of principles of ageing and an extension of lifespan is a key interest of biomedical research1,2. Organisms with an extended lifespan might increase our insight into biological mechanisms of (healthy) ageing3. Rodents showing unusual longevity have consequently become valuable models for studying molecular mechanisms of geroprotection4. One of those emerging model taxa is Nannospalax galili (also called Spalax or blind mole rat, BMR), a member of the superspecies Spalax ehrenbergi from the muroid family of Spalacidae, a group of fossorial rodents that diverged from the phylogenetic lineages leading to mouse and rat approximately 47 million years ago5,6. While laboratory rats have an average lifespan of 3 years7, Spalax reaches up to 21 years without displaying clear signs of ageing or age-related disorders8. Moreover, in 50 years of Spalax research, there has never been a reported case of spontaneous tumor formation among thousands of observed individuals. Even exposition to potent carcinogens failed—with very rare exceptions—to induce tumor growth in Spalax, whereas 100% of exposed mice and rats in the same study developed cancer8. Interestingly, the blind mole rat shares these important traits with another fossorial rodent, the African naked mole rat (NMR) Heterocephalus glaber9,10, which is adapted to moderate, chronic hypoxia in its populated, eusocial underground colonies and which reaches lifespans of up to 30 years and develops cancer only in rare cases9. Since Spalax and H. glaber belong to different clades of the rodent phylogeny, they probably evolved these traits convergently during speciation in adaptation to hypoxia5,11. Hence, these fossorial rodents, which efficiently cope with the multiple stresses of underground life, are now studied to reveal the adaptive processes involved in tissue oxygenation, tumor formation and longevity9,12,13.Spalax individuals inhabit self-dug burrows to avoid predators and spend most of their lifes solitarily14. Due to flooding of its burrow system during seasonal rainfalls, N. galili overcomes the most severe fluctuations in oxygen and carbon dioxide among all S. ehrenbergi subspecies. O2 concentrations of only 7% were measured in flooded Spalax burrows15. Under laboratory conditions, N. galili survived 14 h at 3% O2, whereas rats died after 2-4 h16. This efficient hypoxia tolerance is realized in Spalax by a combination of physiological adaptations. As a shared feature of many subterranean rodents, Spalax displays a very low basal metabolic rate even under normoxic conditions, which limits oxygen consumption14. In addition, Spalax possesses a significantly smaller total skeletal muscle mass than comparatively sized rodents17 Oxygen distribution, on the other hand, is increased by elevated capillary and mitochondrial density18. In contrast to S. ehrenbergi subspecies living in arid, well-oxygenated and permeable soil, N. galili shows elevated concentrations of red blood cells and hemoglobin19. Together, these features are thought to guarantee unimpaired body functions even under challenging atmospheric oxygen conditions. However, it is suspected that many more adaptations are found on the molecular level which are hitherto incompletely understood.It has been hypothesized that the three combined phenotypes - hypoxia tolerance, longevity, tumor resistance - could be mechanistically linked at the genetic level20,21,22,23,24. Comparative genomics and molecular analyses have provided evidence that positively selected sequence changes in coding genes, but also changes in gene regulation have contributed to adaptation in the subterranean rodents5,13,25,26,27. In Spalax, for example, amino acid replacements in the crucial tumor suppressor p53 render the protein unable to activate some of its apoptosis-regulating target genes28, while over-activating targets mediating cell-cycle arrest. This suggests an adaptive strategy to prevent hypoxia-induced cell death, while allowing DNA repair to take place. Pseudogenization associated with visual perception and DNA repair is another evolutive feature of Spalax observed at the genomic level29. Other gene-focused studies revealed substantially increased mRNA and protein expression of key players involved in a variety of functions in Spalax organs in comparison to tissues from the non-hypoxia-adapted rat, e.g., elevated Hif1a30 and Epo levels16, overexpression of respiratory globin proteins31, ROS defense enzymes32, and tissue vascularization genes33. Large-scale transcriptome analyses of selected Spalax organs uncovered additional gene-regulatory adaptations, e.g., an attenuated transcriptional response of energy metabolism genes under normoxic conditions in brain26 and in response to hypoxia in liver tissue13. Both, Spalax liver and brain transcriptomes (in comparison to rat) revealed constitutively increased expression of important DNA damage repair genes, such as the base excision repair-associated glycosylase gene Endonuclease VIII-like 2 (Neil2), the Xeroderma pigmentosum-associated (Xpa) nucleotide excision repair gene and the progeria-associated helicase Werner (Wrn)13,26. In addition, replicative stress response genes interacting with Wrn, such as the Ataxia telangiectasia and Rad3 related gene (Atr), the Serine/Threonine Kinase Ataxia Telangiectasia Mutated gene (Atm) and members of the Fanconi-Anemia pathway showed increased mRNA expression in Spalax livers, compared to rat13. Altogether, these genetic and gene regulatory changes may ensure genome stability in Spalax when severe hypoxia and subsequent reoxygenation challenge tissues and cells, e.g., via production of toxic ROS26.The results above support the idea that gene expression networks in Spalax, which e.g. increase genome stability may explain decreased ageing and tumorigenesis as “side effects” of hypoxia adaptation. We hypothesize that such crucial adaptive features are found body-wide, i.e. in the transcriptomes of many vital tissues. To investigate this, we generated and analyzed RNA-Seq data from additional Spalax and rat organs (kidney, spleen) obtained from specimen subjected to hypoxia (6% O2) and normoxia as controls. Including previous RNA-seq data, we conducted meta-transcriptomic comparisons across Spalax tissues to elucidate tissue-overarching versus -specific gene regulatory responses and discuss our results with respect to molecular mechanisms of ageing and cancer.ResultsTranscriptome data metricsTo unveil the potential adaptive mechanisms that allow the long-lived rodent N. galili to tolerate severe hypoxia on the gene regulatory level, we subjected n = 3 individuals of both species to either 6% O2 or to normoxia (room air) for six hours and dissected liver, kidney, and spleen tissues. RNA was isolated to generate transcriptome-wide mRNA expression profiles by Illumina sequencing. Between 10 and 54 million reads per dataset were mapped to their respective reference genomes (Supplementary Dataset S1). For validation, expected organ-specific gene expression patterns as described by Yu et al. (2014) were clearly observable within our data (Supplementary Material Fig. 1). We detected expression in at least one of the biological replicates for 16359 genes in Spalax and 17741 genes in rat in all analyzed tissues under both oxygen conditions (Supplementary Dataset S2). When comparing the gene expression for each tissue via principal component analyses (PCA), we observed that replicate samples (n = 3) of the same species subjected to the same oxygen conditions clustered together as expected (Supplementary Material Fig. 2). In all tissue comparisons, the interspecies variance between Spalax and rat was larger than between the hypoxia- and normoxia-treated samples of the same species, implying that hypoxia induction had less overall impact on the transcriptional profile than species divergence.Quantitative transcriptome differences in the hypoxia response of Spalax and ratTo study the magnitude of the gene regulatory response after hypoxic stress, we identified differentially expressed genes (DEGs) between normoxia and hypoxia for all three tissues in both species (padj. < 0.05). In all tissues, the hypoxia response was notably stronger in rat than in Spalax (Fig. 1A). After hypoxia treatment, about 1.5 times more genes were dysregulated in rat liver, 1.45 times more in rat kidney and 1.6 times more in rat spleen, as compared to Spalax. Figure 1B–D shows tissue-wise scatterplots of the hypoxia-induced genes for rat and Spalax. In particular in spleen, many more dysregulated genes (red dots) were observed in rat (3326) compared to Spalax (2055). Together, these data indicate a stronger response to hypoxic stress of rat tissues than in the hypoxia-adapted Spalax.Fig. 1: Quantitative transcriptome differences in hypoxia-exposed Spalax and rat.A Number of Spalax and rat DEGs after hypoxic stress exposure in liver, kidney, and spleen. Red = upregulated, blue = downregulated genes. B–D Scatterplot visualization of DEGs following hypoxic treatment in Spalax and rat liver, kidney and spleen. In each graph, log2 TPM fold-changes are plotted. Differentially expressed genes (padj <0.05) are marked in red.Full size imageIn addition to the observed quantitative differences, the sets of genes induced by hypoxia differed qualitatively between the two species. From those genes differentially expressed under hypoxia in Spalax, only 46% (liver), 33% (kidney), and 37% (spleen) were also regulated in the corresponding rat organ (Supplementary Material Fig. 3). This suggests that both organisms transcribe a large set of species-specific hypoxia-responding genes. Ca. 32% (1677) and 39% (2897) of hypoxia-induced DEGs were common to at least two organs in Spalax and rat, respectively (Fig. 2). These genes might represent a systemic transcriptional response towards hypoxic stress acting across different tissues.Fig. 2: Venn diagrams illustrate the number of hypoxia-regulated genes compared between the three organs studied (padj < 0.05).Genes are counted separately for Spalax (A) and rat (B).Full size imageTo interpret the biological function of the hypoxia-responsive genes, pathway enrichment analyses were performed on the DEGs of all three organs in Spalax and rat (Supplementary Dataset S3). Identified pathways describing biological processes were summarized for comparison. In liver tissue of both species, enriched hypoxia-activated pathways were associated with traits involved in ageing, including cell death (e.g. “Senescence Pathway”), proinflammatory processes (e.g. “Acute Phase Response Signaling”), hypoxia stress response (e.g. “HIF1α Signaling”), extracellular matrix (ECM) and angiogenesis (e.g. “Tumor Microenvironment Pathway”) (Fig. 3). In rat, but not in Spalax, enrichment in metabolism-associated terms like “Gluconeogenesis I” and “Glycolysis I” was detected in the hypoxia-activated genes.Fig. 3: Selected pathways induced by hypoxia in the liver of Spalax and rat (padj < 0.05).Pathway enrichment was predicted based on the genes differentially expressed at hypoxia and normoxia (HxNx). Pathways are grouped into larger categories (left margin) for improved overview (orange = predicted activation (z-score > 1).Full size imageIn kidney, in addition to the regulation of general cellular signaling processes, hypoxic activation of the following pathways was observed in both species (Fig. 4): “Wound healing Signaling Pathway”, angiogenesis and ECM formation (”VEGF-Signaling”, “Tumor Microenvironment Pathway”), hypoxia- and oxidative stress-associated processes (“HIF1α Signaling”, “NRF2-mediated Oxidative Stress Response”), proinflammatory pathways (“Acute Phase Response Signaling”) and metabolism-associated terms (e.g. “Type II Diabetes Mellitus Signaling”). Particularly in Spalax, the activation of many ageing- and cell-death-associated terms was detected (e.g., “MYC Mediated Apoptosis Signaling).Fig. 4: Selected pathways induced by hypoxia in the kidney of Spalax and rat (padj < 0.05).Pathway enrichment was predicted based on the genes differentially expressed at hypoxia and normoxia (HxNx). Pathways are grouped into larger categories (left margin) for improved overview (orange = predicted activation (z-score > 1).Full size imageIn spleen tissue of both species (Fig. 5) hypoxia-responsive DEGs were enriched in terms related to the activation of signaling pathways that are associated with pathologies and damage (“Wound Healing Signaling Pathway”, “Osteoarthritis Pathway”, “IL6-Signaling”). In notable contrast to the other organs, hypoxia response pathways (“HIF1α Signaling”,” Nitric Oxide Signaling in the Cardiovascular System”) were specifically activated in rat, but not in Spalax.Fig. 5: Selected pathways induced by hypoxia in the spleen of Spalax and rat (padj < 0.05).Pathway enrichment was predicted based on the genes differentially expressed at hypoxia and normoxia (HxNx). Pathways are grouped into larger categories (left margin) for improved overview (orange = predicted activation (z-score > 1).Full size imageThe results indicated that many processes important for a hypoxia response are activated in both, the adapted rodent Spalax and the non-adapted rat. To determine whether this apparent functional similarity in the hypoxia response is orchestrated by the regulation of the same genes in Spalax and rat, we compared the regulated genes from hypoxia-responsive pathways like “HIF1α Signaling” or “NRF2-mediated Oxidative Stress Response” between Spalax and rat in all analyzed organs (Supplementary Dataset S4). In liver and kidney, both species regulated around half of the n = 139 analyzed hypoxia-associated genes, while in spleen only n = 50 genes were hypoxia-responsive in the rat and an even lower number (n = 19) in Spalax. Only n = 13 genes showed a hypoxia-induced regulation in both species and all three organs (irrespective of their direction of regulation) and might thus fulfill key roles in the stress response (Supplementary Dataset S5). In an organ-wise comparison, n = 28 genes in liver, n = 21 in kidney and no genes in spleen were regulated in the same direction in both species. On the other hand, n = 39 genes in liver, n = 40 in kidney and n = 42 in spleen were regulated in only one of the two species: e.g., Hypoxia-inducible factor-2alpha (Epas1) was only upregulated in rat, but not Spalax liver. Among this set were also genes that showed an opposing direction of regulation in Spalax and rat like the MAP kinase-interacting serine/threonine-protein kinase 1 (Mknk1), which was hypoxia-upregulated in rat liver, but downregulated in Spalax.In summary, the transcriptome data revealed that around 1.5-fold more genes are regulated in rat compared to Spalax in response to hypoxia. The majority of regulated genes differed between Spalax and rat, yet this resulted in the activation of similar biological pathways in liver and kidney. In case of the hypoxia and stress-response pathways, many genes showed species-specific response patterns. In spleen, the hypoxia-response was attenuated in both species, but especially in Spalax. Our results thus indicate an adapted gene regulatory hypoxia-response in Spalax, which involves less (and different) genes.Differences in gene expression under normoxic conditions between Spalax and ratDifferences between Spalax and rat in gene regulation may not only manifest themselves during an acute response to hypoxic stress, but also “constitutively” at normoxic conditions. In all three organs, more than half of the genes were differentially expressed between Spalax and rat at normoxia, the number being evenly distributed between higher- and lower-expressed genes (Fig. 6). Similar to the patterns observed in the transcriptional hypoxia response, most genes (54% of all DEGs, n = 8215 of 15196) were either up- or downregulated between the two species in at least two organs. Around 26% (n = 2000) and 23% (n = 1721) of the DEGs were upregulated and downregulated, respectively, in all three Spalax vs. rat organs, indicating body-wide adaptations on the gene regulatory level in the two species (Fig. 6).Fig. 6: Venn diagram showing the number of differentially expressed genes observed in kidney, liver, and spleen tissue in the interspecies comparison between Spalax and rat kept at normoxic conditions.A Genes showing higher expression in Spalax than in rat. B Genes showing lower in Spalax than in rat. Criterion for differential expression was a value of padj < 0.05. The large number of genes, which are differentially expressed in all three tissues, indicates across-organ adaptations on the gene regulatory level in Spalax compared to rat.Full size imageInterspecies DEGs at normoxia were analyzed for pathway enrichment to interpret their biological function (Supplementary Dataset S3). We observed activation of DNA repair-associated pathways (e.g., “Base excision repair”) in multiple Spalax organs, compared to rat (Fig. 7). In parallel, processes associated with proliferation (“mTOR signaling”), proinflammatory pathways (“Neuroinflammation Signaling Pathway”), coagulation (“Thrombin Signaling”) and extracellular matrix components (“Dermatan Sulfate Degradation”) were predicted as inactivated. Enrichment results in Spalax liver tissue specifically showed activation of metabolic pathways (“Glycogen Degradation II“, “Glycogen Degradation III“), but inactivation of hypoxia-associated pathways (“HIF1α Signaling”). In kidney, “HIF1α Signaling” was predicted as activated, while carbohydrate metabolism-associated pathways (“Glycolysis”,” Gluconeogenesis”, “Glycogen Degradation”) were predicted as inactivated. In kidney, we specifically observed enrichment of genes associated with lipid metabolism (“Triacylglycerol Biosynthesis” and “Fatty Acid β-oxidation I”). Notably, in spleen, we did not observe activity changes for DNA-repair- or metabolism-associated processes.Fig. 7: Selected pathways differentially expressed between Spalax and rat (padj < 0.05) in normoxic liver, kidney and spleen.Pathways are grouped into larger categories (left margin) for improved overview (orange = predicted activation in Spalax (z-score > 1), blue = predicted inactivation in Spalax (z-score < −1)).Full size imageIn conclusion, the enrichment analysis of pathways differentially expressed between Spalax and rat at normoxic conditions pinpointed several processes that showed the same activation status across Spalax organs. Detailed inspection of such pathways, as exemplified here by mTOR (Fig. 8), a key metabolic regulator implicated in ageing, in fact revealed that the same genes were differentially expressed in more than one Spalax organ with a similar direction and strength of regulation, suggesting that those pathways are regulated body-wide in Spalax.Fig. 8: Genes contributing to mTor signaling and their respective fold-changes in differential expression between Spalax and rat at normoxia.Orange = elevated expression in Spalax compared to rat (Log2FC > 0), blue = lower expression in Spalax (Log2FC < 0), white = differential expression not significant (padj. > 0.05).Full size imageValidation of the RNA-Seq analysis via qRT-PCRBased on the results of the functional enrichment analyses and selection of important candidate genes associated with ageing, cancer and hypoxia adaptation, we studied the hypoxia inducibility and interspecies expression differences of a selected gene set representing different functional categories (A2m, Atr, Cisd2, Fgf21, Wrn, Xpa, Rcan1, Gpnmb, Fen1, Hmox1, Vegfa, Pnkp) by qRT-PCR (Supplementary Material Tab. 1). We could confirm the direction of all relevant gene expressional changes in the respective organs.DiscussionSubterranean rodents like the blind mole rat (Nanno)Spalax spec. and the naked mole rat H. glaber exhibit a fascinating combination of phenotypes, linking their remarkable longevity and tumor resistance to an adaptation to environmental hypoxia9,12,13. Previous studies, addressing single genes and complete genome sequences revealed interesting mutations in candidate coding sequences, but also suggested that differences in gene regulation may underlie adaptations in the subterranean species13,27,34,35,36. Here, we generated and compared RNA-Seq datasets from vital organs of the hypoxia-tolerant blind mole rat S. galili and – for comparison – from the hypoxia-sensitive rat. We hypothesized that, in particular, organ-overarching patterns of gene expression might yield insights into the molecular pathways that are responsible for the complex phenotype of Spalax. Note as a limitation that, in this study, female Spalax individuals were compared to male rats. To identify and remove potential sex biases in gene expression, we analyzed publicly available RNA-Seq datasets from male and female rats (for details see Methods section). Pathways containing more than 10% sex-biased DEGs were then excluded from our subsequent analyses and data interpretation. However, the possibility of residual sex bias in our results cannot be entirely ruled out.Applying a threshold of 2-fold, the number of differentially expressed genes responding to hypoxia in liver, kidney and spleen of Spalax and rat, respectively, ranged between ca. 800 and 2100, in line with previous transcriptome studies of hypoxia-exposed rodents13,37,38. At normoxia, 6426 (liver), 7173 (kidney) and 6205 (spleen) genes were significantly differently expressed at least 2-fold (log2fc > 1; padj < 0.05) between the two species (Supplementary Dataset S6). This is mirrored by a liver transcriptome comparison of the hypoxia-adapted naked mole rat and its non-adapted hypoxia-sensitive relative, the guinea pig, reporting differential expression of > 4000 genes39. Considering differences in the experimental setup and the closer phylogenetic distance between naked mole rat and guinea pig (39.5 mya) compared to Spalax and rat (47.4 mya)40, our results on the interspecies expression differences are thus in the expected range. Measurements of RNA levels can serve as a reliable proxy for estimating gene expression and can account for approximately two thirds of the observed protein abundance in mammalian cells41. In the naked mole rat, a similar correlation between the RNA and the protein level in an interspecies comparison has been previously reported39. However, the same study showed that for some key adaptive processes such as oxidative phosphorylation, differential expression only manifested itself at the protein level. In the future, it will thus be necessary to integrate our transcriptomics approach with protein data to more comprehensively understand adaptation of gene regulation in Spalax.As a general trend in all three organs studied, a significantly higher number of genes responded to hypoxic stress in the hypoxia-sensitive rat as compared to the hypoxia-tolerant Spalax (comp. Figure 1). At the same time, the majority of hypoxia-responsive genes differed between the two rodents (Supplementary Material Fig. 3). Still, in many cases, pathways associated with similar biological functions such as response to hypoxia, proliferation control, angiogenesis or immune response were activated in both species (Fig. 9). Detailed investigation revealed that often different genes were regulated in the two species within a respective biological pathway, while in other cases the same genes – but showing opposing directions of regulation - were involved in Spalax and rat. These results indicate that, while classical hypoxia response pathways are activated in both rodents upon stress exposure, the adapted Spalax achieves a biological response by utilizing less genes. The involvement of partially different genes in Spalax compared to rat, as exemplified by the mTor signaling pathway, indicates a potentially more sophisticated gene regulatory response in the blind mole rat, which might conserve energy by utilizing less resources for gene expression.Fig. 9: Common and species-specific patterns of transcriptional gene regulation.A Often, similar biological pathways are significantly enriched and activated (plus symbol) in the organs of both, Spalax and rat after hypoxic stress. B A complex picture emerges when differentially expressed genes within pathways are inspected: expression of different genes in Spalax and rat may lead to activation of the same pathway, as exemplified by Epas1 and Mdm2 (red = upregulation). In contrast, the same gene (here: Mknk1) can be regulated in opposite directions in response to hypoxia in the two species (blue = downregulation).Full size imageDifferential gene expression under normoxic conditions is believed to be another important adaptive mechanism in Spalax. Elevated expression levels of important master regulators of hypoxic stress16 and antioxidant defense13,32,42 have been observed in the mole rat, suggesting that such higher expression levels of key genes and pathways, already present at normoxia, might enable Spalax to react substantially faster towards an acute onset of severe hypoxic stress, which the rodents meet in their underground habitat during heavy seasonal rainfalls13,15 This can clearly be a useful adaptive strategy for biological processes like antioxidative defense or DNA-repair, which guarantee genomic integrity. For other biological processes such as inflammation or energy metabolism, a constitutive activation could be wasteful in energy expenditure and would not be in line with the energy-conserving phenotype of many subterranean species including Spalax14. It is therefore not surprising that we predicted an inhibition of biological pathways associated with inflammation, energy metabolism or angiogenesis in Spalax compared to rat (Figs. 7, 10). These gene-regulatory patterns were in many cases observed in more than one organ, often utilizing the same pathway components (e.g. in mTor signaling; Fig. 8). This indeed confirms the hypothesis that body-wide adaptive gene-regulatory mechanisms are acting in Spalax on the transcriptional level.Fig. 10: Predicted status of pathways associated with the Spalax phenotype at normoxic, unstressed conditions.Activation is designated by a red plus sign, inactivation by a blue minus sign. Note that most pathways are affected in at least two or even three organs, indicating organ-overarching adaptations.Full size imageA more detailed look at the gene regulatory pathways involved in the differences between Spalax and rat revealed a number of observations that are of potential biological interest, in particular for the evolution of the complex long-lived phenotype. In both species and all three organs, we observed the enrichment of terms associated with the control of cell death and proliferation following hypoxic stress exposure. This is an expected response in sensitive species, since hypoxia often leads to the elimination of potentially damaged cells to prevent carcinogenesis43,44,45. Hypoxia-adapted taxa, in contrast, might want to avoid unfavorable cell loss when encountering environmental stress repeatedly. In fact, a tight control of apoptosis has evolved in Spalax via mutations in the key tumor suppressor p53, some of which mimic mutations found in P53 of human tumors28 or are convergently present in distantly related hypoxia-tolerant species46. These changes reduce Spalax p53’s ability to induce apoptosis-regulating target genes like Apaf1 and Noxa28. In addition, Mdm2, a negative regulator of p53, was significantly higher expressed in Spalax than in rat in hypoxic muscle, brain, and heart tissue47. In our study, we observed a similar overexpression of Mdm2 in Spalax liver and kidney, implying a body-wide involvement of this gene in controlling p53-mediated processes in the blind mole rat. While Spalax p53 appears to have a reduced stimulatory effect on pro-apoptotic genes, in vitro reporter gene assays revealed that it over-activates genes associated with cell-cycle arrest like Cdkn1a and Pten, suggesting an adaptive strategy that promotes temporary cell cycle arrest over apoptosis to prevent excess hypoxia-induced cell loss28. Interestingly, we detected a lower transcription of Pten and Cdkn1a in all three Spalax organs, possibly down-modulating p53-induced activation under normoxic conditions. Thus, to improve insight into complex pathways, mutations in master regulators should be interpreted together with mRNA levels of their respective target genes. An additional factor that plays an important role in the cellular response to hypoxia is the calcium influx, which in turn, activates apoptosis and affects the transcription landscape. A study on Spalax cultured cells revealed a significantly reduced calcium influx during hypoxia48, implying a crucial role of maintaining homeostasis, rather than ‘overreacting’ under stress.As expected, we observed an activation of classical hypoxia response terms like HIF1α signaling in liver and kidney in both species. Moreover, this pathway was predicted as already activated under normoxic conditions in Spalax kidney. Adaptations in gene regulation in the kidney have previously been proposed to play a role in Spalax hypoxia resistance, as the kidney is the major site of erythropoietin production. In this context, an elevated HIF1a and EPO production was measured in Spalax kidney16, which likely contributes to the mole-rat’s elevated erythrocyte numbers19. Under hypoxia, non-adapted rodent species additionally show stress-erythropoiesis in the spleen38, and the spleen serves as a reservoir of red blood cells in many diving species49. Interestingly, Spalax spleen showed a lower transcriptional response to hypoxia and fewer interspecies differences in pathway activation under normoxic conditions, as compared to the other organs. This hints at an obsolescence of splenic stress erythropoiesis in Spalax due to its well-adapted erythropoietic physiology.An important cause of cell and DNA damage and thus cancer- and ageing-related pathologies is oxidative stress. While basal levels of oxidative stress can be found in essentially all aerobic organisms, they are substantially increased in subterranean, diving, and high-altitude species50. Spalax faces acute cyclical changes in oxygen levels in its habitat when performing energy-consuming digging activities in confined tunnels at one moment and entering better-ventilated burrow areas at the next. Moreover, when its underground burrow systems are flooded by heavy rainfall, Spalax encounters oxidative stress in cycles of acute severe hypoxia followed by re-oxygenation, which produces harmful ROS51. Hypoxia-induced upregulation of anti-oxidant defense genes like Hmox1 can, therefore be observed in many species in response to hypoxia52. Our RNA-Seq data confirm earlier qPCR results32 showing elevated mRNA level of Hmox1 in Spalax and rat in response to hypoxia, indicating that despite the many species-specific differences in hypoxia gene regulation, also common, conserved transcriptomic responses exist.ROS that are not sufficiently buffered by the radical scavenging system represent a severe threat to genomic integrity, producing various types of DNA damage like nucleotide modifications or crosslinks53,54,55. Furthermore, severe hypoxia induces S-phase arrest and dNTP depletion, leading to additional replicative stress56,57,58. Consistent with this, we observed an enrichment of DNA damage response-associated processes among the genes that were differentially expressed already at normoxia between Spalax and rat. Of note, “Base excision repair” and “Role of BRCA1 in DNA damage response” were predicted as activated in liver and kidney of the mole rat. Previous analyses of RNA-Seq data from Spalax liver13 and brain26 had indicated an elevated expression in the mole rat of key genes from the Fanconi Anemia pathway, which orchestrates different DNA repair activities59. In the present study, we could extend these results to kidney and spleen (Fig. 11), which indicates an organ-overarching, body-wide adaptation of DNA repair capacity in Spalax. Looking further for candidate genes of special biological relevance, important DNA repair-associated genes like the master activator Ataxia telangiectasia and Rad3 related (Atr), the breast cancer type 2 susceptibility protein (Brca2) and the Werner syndrome ATP-dependent helicase (Wrn) were found overexpressed in Spalax liver, brain and spleen compared to rat (present study; as well as in brain tissue26, again suggesting body-wide elevated activities of their gene products. Recently, DNA damage repair assays performed on skin fibroblasts isolated from S. carmeli, another subspecies of the S. ehrenbergi taxon, indeed revealed an elevated resistance against H2O2, the topoisomerase inhibitor etoposide and UV-C radiation-induced damage compared to rat fibroblasts60, suggesting an accelerated repair of those lesions in Spalax. Since DNA damage is causal to both cancer development61) and ageing62, we hypothesize that an improved capacity of repairing DNA lesions may substantially contribute to Spalax’s anti-cancer and longevity phenotype.Fig. 11: Differential gene expression heatmap for the Fanconi anemia pathway in liver, kidney and spleen of Spalax (Nga) and rat (Rno).Respective organs are represented from left to right (liver>kidney>spleen) as squares inside the gene name boxes. Genes with higher expression in Spalax than in rat (marked orange) are overrepresented compared to genes with a lower expression in Spalax vs. rat. Differential expression is represented as the Log2 FC of expression levels between Spalax and rat at normoxia (Log2 FC (Nga Nx/Rno Nx)).Full size imageAn important physiological adaptation to its hypoxic habitat is Spalax’s low basal metabolic rate, which reduces its oxygen consumption and, thereby, oxidative stress17. Such hypometabolism, which is thought to promote longevity and prevent cancer63, was in fact previously inferred from transcriptome analyses of Spalax brain, where lower mRNA levels of genes associated with energy metabolism were found compared to rat13,26. In support of this, we detected differential expression of genes associated with cellular energy metabolism in liver, kidney and spleen between Spalax and rat at normoxia. For example, key genes contributing to oxidative respiration like cytochrome c oxidase showed lower normoxic mRNA expression for most subunits in all analyzed Spalax organs (Fig. 12).Fig. 12: Differential gene expression heatmap of genes contributing to oxidative phosphorylation in liver, kidney and spleen of Spalax (Nga) and rat (Rno).Respective organs are represented from left to right (liver>kidney>spleen) inside the gene name boxes. Genes with lower expression in Spalax than in rat (blue) are overrepresented compared to genes with a higher expression (orange). Differential expression is represented in form of the Log2 FC of expression levels between Spalax and rat at normoxia (Log2 FC (Nga Nx/Rno Nx)).Full size imageImportantly, the transcriptome data predicted an inactivation of mTor signaling in all three Spalax organs compared to rat. Due to its involvement in ageing-related processes like nutrient sensing, proteostasis and mitochondrial dysfunction, the mTOR cascade is considered a major player in life span regulation64,65,66 and mTOR inhibitors like rapamycin are researched for their potential use in lifespan extension67. Interestingly, recent studies in the hypoxia-adapted naked mole rat and red-eared slider turtles68 suggested an activation of mTOR in response to hypoxia, which may aid metabolic reprogramming in favor of anaerobic metabolism. Why this is not observed in Spalax remains unclear. However, specific amino-acid replacements in mTORs functional domains were observed in Spalax 13, suggesting potential adaptations on the functional protein level in the blind mole rat. A constitutive inactivation of the mTOR pathway, as inferred from the transcriptome level in Spalax, could mimic the effect of anti-cancer drugs69,66 and thus contribute to the mole rat’s cancer-resistance.Hypoxia also plays a key role in the regulation of immunity and inflammation since it is a hallmark of inflamed, infected or damaged tissue70. In this context, inflammation can induce the activity of hypoxia-response pathways, and hypoxia may modulate inflammatory signaling71. It is therefore not surprising that we detected an enrichment of genes associated with the regulation of inflammatory processes among the hypoxia-responsive genes in all organs of both species. The same was observed in brain tissue in a previous study72. Sterile inflammation, which is induced by the secretome of senescent cells, consisting of inflammatory cytokines, chemokines and growth factors, knowingly plays a role in age-related disorders and cancer73,74. Stress assays in Spalax fibroblasts suggest that the positive feedback loop of IL1α-NF‐κB, which is an upstream regulatory machinery of the inflammatory response, appears to be impaired75. Accordingly, we detected in all analyzed Spalax organs a decreased expression of proinflammatory factors, some of which are also part of the secretome of senescent cells (e.g. Il1a, Il6r and Cxcl1). We conclude that a lower expression of proinflammatory genes under normoxic conditions in healthy organs may produce an “inflammation-free” phenotype in Spalax, possibly contributing to the species’ longevity.As an important caveat to our (and other) transcriptomic interspecies comparisons, it might be argued that the observed differences between taxa simply reflect phylogenetic distance, but not necessarily adaptive processes. A way out of this analytical dilemma can be the inclusion of additional species with and without the studied phenotypes, which may reveal instances of phenotypic convergence, and thus indirectly point at “adaptation in action”. In rare cases, convergent adaptation is realized by ‘strict convergency’, e.g., by mutations affecting the same gene in two species adapted to the same selective pressure. One example is the presence of functionally similar amino acid replacements in the proton-gated nociceptor sodium channel Nav1.7 (Scn9a) of Spalax and several other species, which are adapted to hypercapnic environments5. More often, convergency is realized on a functional level by sometimes utilizing different genes (but from the same pathways) to achieve similar phenotypes. In the context of hypoxia, adapted species that also possess a longevity phenotype like naked mole rats or whales76 consistently show differential expression in genes associated with processes like oxidative stress response and ROS defense, metabolism, and DNA-repair compared to non-adapted, related taxa,34,39,77,78. Since we observed gene-regulatory changes influencing the very same biological processes in Spalax, we infer that these differences in gene expression are, to a high degree, indeed adaptive. There are even examples for strict convergence in overexpression of the same genes in Spalax and naked mole rat, as exemplified by the antiprotease alpha2-macroglobulin A2M13,79 and DNA-repair genes like ATR34. Studying and comparing functional and direct convergencies in adapted species like Spalax and naked-mole rat can thus provide valuable insight into the evolution of complex phenotypes and help to pinpoint different ways for modulating physiological processes like hypoxia response, tumor inhibition and geroprotection.We conclude that the complex cancer- and ageing-resistant phenotype of Spalax is likely influenced by transcriptomic changes, which act overarchingly across vital organs and tissues. These transcriptomic changes involve genes from fundamental biological processes like the control of cell death, ROS defense, DNA repair, energy metabolism, immune response, and angiogenesis. Many of these processes link the transcriptional hypoxia response to pathways determining the defense against cancer, the maintenance of genomic integrity, and the emergence of longevity. Our data thus provide additional evidence for the hypothesis that these biomedically most interesting phenotypes have evolved as consequences of Spalax’s adaptation to environmental hypoxia. Future studies on cis- and trans- regulatory mechanisms orchestrating adaptation on the gene regulatory level in Spalax will provide further insights in the evolution of the blind mole rat’s remarkably complex phenotype.Materials and MethodsSample preparationSix N. galili specimens were captured from basalt, heavy soil in fields of northern Israel and kept in the animal house of the Institute of Evolution, Haifa University in accordance with the regulations of the Israel Nature and Park Authority, Science and Conservation Unit. To investigate the impact of O2 deprivation on the transcriptomes, three Spalax (approximately 4 years old females, 150–250 g) and three rat individuals (Sprague Dawley, 4 month old males, 350–400 g) were exposed to 6% O2 for six hours in 70x70x50 cm chambers divided into separate cells at a gas flow rate of 3.5 l/min. No animals died during hypoxia exposure. Afterwards, they were sacrificed immediately by Ketaset CIII (Fort Dodge, USA) injection. From each individual, the kidney and spleen were dissected. Liver samples from the same animals were analyzed prior to this study13. Animal experiments were approved by the University of Haifa Ethics Committee (Permit #193/10). RNA was extracted from the kidney and spleen with the RNeasy lipid tissue kit (Qiagen, Hilden, Germany), including an RNase-Free DNase I treatment. RNA integrity was validated on a 2100 Bioanalyzer with the Agilent RNA 6000 Nano kit (Agilent Technologies, Santa Clara, USA). Sequencing libraries were generated with the Illumina TruSeq RNA library Prep kit v2 (StarSEQ, Mainz, Germany) and sequenced as 100 bp paired-end reads on an Illumina HiSeq2500 (Biology Department, Johannes Gutenberg University Mainz, Germany). Sequence data were uploaded for public availability to ENA (Acc. No. PRJEB64658).RNA-Seq analysesRNA-Seq data were quality-trimmed and filtered for sequencing adapter sequences with Trimmomatic v.0.3.680 using a sliding window approach with a window size of 4 nt and a quality score cutoff of 20. The first 13 nt of each read were cropped, and reads shorter than 20 bp discarded. The processed reads were mapped with the STAR aligner v.2.5.381 to the annotated reference genomes of Rattus norvegicus (6.0.92; https://ftp.ensembl.org/pub/release92/fasta/rattus_norvegicus/dna/Rattus_norvegicus.Rnor_6.0.dna.toplevel.fa.gz) and N. galili (1.0.92; https://ftp.ensembl.org/pub/release92/fasta/nannospalax_galili/dna/Nannospalax_galili.S.galili_v1.0.dna.toplevel.fa.gz). Count data files and index files were generated using an sjdbOverhang of 87. The mismatch maximum was set to 0.04 per base, and a multimap filter was applied with a maximum of 20 matches. Gene IDs of Spalax and rat were matched by identifying orthologue Ensemble IDs with BLASTP (BLOSUM45 matrix; cutoff length 100 amino acids; minimum identity 30%). Only genes that matched the top BLAST hits as judged by the E-value for both search directions, i.e., Spalax query searches in rat databases and vice versa, were used.Differential gene expression analysis between samples was conducted with DeSeq2 v.1.30.182. Statistical significance was determined by applying the default Wald-test of DESeq2. A padj. cutoff of < 0.05 was used to define differential gene expression. For interspecies comparison, counts were normalized with regard to the exon sum length using the normMatrix function as part of the DESeq2 package. Differentially expressed genes were mapped onto relevant KEGG pathways for visualization using the pathview R package83. Principal component analyses were generated with R, using the function plotPCA. Resulting gene lists were compared between species and between different O2-conditions using the R VennDiagram package84. We validated organ specificity by projecting our gene lists on known, organ-specific gene expression patterns85. RNA-Seq data were analyzed with the use of QIAGEN Ingenuity Pathway Analysis (QIAGEN Inc., https://digitalinsights.qiagen.com/IPA86); A padj. cutoff of p < 0.05 was used to define pathway enrichment and a z-score cutoff of < −1 or > 1 for inactivation and activation, respectively.A previously published liver RNA-Seq dataset 13 was re-evaluated using the same analytical tools and parameters to adequately compare the liver transcriptome to kidney and spleen. Potential biases in liver gene expression (sex, species, strain, age) were already accounted for by Schmidt et al. (2017) via comparison with public rat, mouse, and human RNA-Seq data. Based on these comparisons, a substantial effect of such biases on interspecies- and hypoxia-related gene expression changes was considered improbable. In the present work, we repeated and extended this analysis, in particular to identify the potential influence of sex-biased gene expression. Public RNA-Seq data of kidney, liver, and spleen tissue of n = 4 male and n = 4 female rats85 were processed in the same way as in our interspecies and hypoxia-normoxia comparisons in order to identify sex-biased genes. A gene was deemed sex-biased if it met the following criteria: padj < 0.05; |Log2FC|> 1. We then inspected the enriched pathways, which emerged in our interspecies and hypoxia-related comparisons (see Results), for the presence of differentially expressed genes (DEGs) with sex-bias. In the liver, 8 out of the 122 significantly enriched pathways contained more than 10% (and up to 17%) sex-biased DEGs. In kidney and spleen, only 5 out of 111 and 0 of 70 significantly enriched pathways met this criterion of sex-biased expression. Conversely, between 93% and 100% of the significantly enriched pathways in each tissue contained less than 10% of sex-biased DEGs (Supplementary Dataset S7). The few pathways that contained more than 10% sex-biased DEGs were not further discussed, and we therefore consider that our main conclusions are not confounded by marked sex-biased gene expression.Validation of RNA-Seq data by qRT-PCRThe differential expression of a subset of candidate genes was validated by quantitative real-time reverse-transcriptase PCR (qRT-PCR). 600 ng of quality-checked RNA of all samples was used to synthesize cDNA with the Superscript III RT kit (Thermo Fisher, Waltham, USA). To generate standard-curves for absolute RNA copy number quantification, plasmids containing corresponding amplicons were generated: PCRs were conducted in a peqSTAR 96 Thermocycler (Peqlab, Erlangen, Germany) with the TrueStart Taq DNA Polymerase kit (Thermo Fisher, Waltham, USA) and purified with the Wizard SV Gel and PCR clean up System (Promega, Madison, USA). PCR products were cloned into pGem-T Easy vectors (Promega), transformed in E. coli DH10B cells, plated out and purified with the GeneJET Plasmid Miniprep Kit (Thermo Fisher). Plasmid inserts were verified by Sanger sequencing (StarSEQ, Mainz, Germany). qRT-PCR was conducted with the GoTaq qPCR Master Mix (Promega) at a total volume of 10 µl in an ABI 7500 Fast Real Time PCR cycler (Applied Biosystems). Serial 10-fold dilutions of plasmids were applied for absolute quantification of samples, measured in triplicates.
Data availability
Sequencing data have been deposited in the European Nucleotide Archive (ENA) under project accession number PRJEB64658.
ReferencesToniolo, L., Sirago, G. & Giacomello, E. Experimental models for ageing research. Histol. Histopathol. 38, 597–605 (2023).PubMed
Google Scholar
Campisi, J. et al. From discoveries in ageing research to therapeutics for healthy ageing. Nature 571, 183–192 (2019).PubMed
PubMed Central
Google Scholar
Stenvinkel, P. & Shiels, P. G. Long-lived animals with negligible senescence: clues for ageing research. Biochem Soc. Trans. 47, 1157–1164 (2019).PubMed
Google Scholar
Gorbunova, V., Bozzella, M. J. & Seluanov, A. Rodents for comparative aging studies: from mice to beavers. Age 30, 111–119 (2008).PubMed
PubMed Central
Google Scholar
Fang, X. et al. Genome-wide adaptive complexes to underground stresses in blind mole rats Spalax. Nat. Commun. 5, 3966 (2014).PubMed
Google Scholar
Steppan, S., Adkins, R., Anderson, J. & Thorne, J. Phylogeny and divergence-date estimates of rapid radiations in muroid rodents based on multiple nuclear genes. Syst. Biol. 53, 533–553 (2004).PubMed
Google Scholar
Sengupta, P. The laboratory rat: relating its age with human’s. Int. J. Prev. Med. vol. 4 www.ijpm.ir (2013).Manov, I. et al. Pronounced cancer resistance in a subterranean rodent, the blind mole-rat, Spalax: in vivo and in vitro evidence. BMC Biol. 11, 91 (2013).PubMed
PubMed Central
Google Scholar
Buffenstein, R. et al. The naked truth: a comprehensive clarification and classification of current ‘myths’ in naked mole-rat biology. Biol. Rev. Camb. Philos. Soc. 97, 115–140 (2022).PubMed
Google Scholar
Shepard, A. & Kissil, J. L. The use of non-traditional models in the study of cancer resistance-the case of the naked mole rat. Oncogene 39, 5083–5097 (2020).PubMed
Google Scholar
Steppan, S. J., Adkins, R. M. & Anderson, J. Phylogeny and divergence-date estimates of rapid radiations in muroid rodents based on multiple nuclear genes. Syst. Biol. 53, 533–553 (2004).PubMed
Google Scholar
Lagunas-Rangel, F. A. Cancer-free aging: Insights from Spalax ehrenbergi superspecies. Ageing Res Rev. 47, 18–23 (2018).PubMed
Google Scholar
Schmidt, H. et al. Hypoxia tolerance, longevity and cancer-resistance in the mole rat Spalax - a liver transcriptomics approach. Sci. Rep. 7, 14348 (2017).PubMed
PubMed Central
Google Scholar
Lacey, E., Patton, J. & Cameron, G. Life Underground: The Biology Of Subterranean Rodents. Australian Mammalogy vol. 23 (The University of Chicago Press, Chicago, 2001).Shams, I., Avivi, A. & Nevo, E. Oxygen and carbon dioxide fluctuations in burrows of subterranean blind mole rats indicate tolerance to hypoxic-hypercapnic stresses. Comp. Biochem Physiol. A Mol. Integr. Physiol. 142, 376–382 (2005).PubMed
Google Scholar
Shams, I., Avivi, A. & Nevo, E. Hypoxic stress tolerance of the blind subterranean mole rat: expression of erythropoietin and hypoxia-inducible factor 1 alpha. Proc. Natl Acad. Sci. USA 101, 9698–9703 (2004).PubMed
PubMed Central
Google Scholar
Widmer, H. R., Hoppeler, H., Nevo, E., Taylor, C. R. & Weibel, E. R. Working underground: respiratory adaptations in the blind mole rat. Proc. Natl Acad. Sci. USA 94, 2062–2067 (1997).PubMed
PubMed Central
Google Scholar
Avivi, A., Resnick, M. B., Nevo, E., Joel, A. & Levy, A. P. Adaptive hypoxic tolerance in the subterranean mole rat Spalax ehrenbergi: The role of vascular endothelial growth factor. FEBS Lett. 452, 133–140 (1999).PubMed
Google Scholar
Arieli, R., Heth, G., Nevo, E. & Hoch, D. Hematocrit and hemoglobin concentration in four chromosomal species and some isolated populations of actively speciating subterranean mole rats in Israel. Experientia 42, 441–443 (1986).PubMed
Google Scholar
Mehta, R. et al. Proteasomal regulation of the hypoxic response modulates aging in C. elegans. Science 324, 1196–1198 (2009).PubMed
PubMed Central
Google Scholar
Boretto, J. M., Cabezas-Cartes, F. & Ibargüengoytía, N. R. Slow life histories in lizards living in the highlands of the Andes Mountains. J. Comp. Physiol. B 188, 491–503 (2018).PubMed
Google Scholar
Rascón, B. & Harrison, J. F. Lifespan and oxidative stress show a non-linear response to atmospheric oxygen in Drosophila. J. Exp. Biol. 213, 3441–3448 (2010).PubMed
Google Scholar
Snell, T. W. & Johnston, R. K. Glycerol extends lifespan of Brachionus manjavacas (Rotifera) and protects against stressors. Exp. Gerontol. 57, 47–56 (2014).PubMed
Google Scholar
Strzyz, P. Hypoxia prevents premature ageing in mice. Nat. Rev. Mol. Cell Biol. 24, 581 (2023).
Google Scholar
Gladyshev, V. N., Zhang, G. & Wang, J. The naked mole rat genome: understanding aging through genome analysis. Aging 3, 1124 (2011).PubMed
PubMed Central
Google Scholar
Malik, A. et al. Genome maintenance and bioenergetics of the long-lived hypoxia-tolerant and cancer-resistant blind mole rat, Spalax: a cross-species analysis of brain transcriptome. Sci. Rep. 6, 38624 (2016).PubMed
PubMed Central
Google Scholar
Davies, K. T. J., Bennett, N. C., Faulkes, C. G. & Rossiter, S. J. Limited evidence for parallel molecular adaptations associated with the subterranean niche in mammals: a comparative study of three superorders. Mol. Biol. Evol. 35, 2544–2559 (2018).PubMed
PubMed Central
Google Scholar
Ashur-Fabian, O. et al. Evolution of p53 in hypoxia-stressed Spalax mimics human tumor mutation. Proc. Natl Acad. Sci. USA 101, 12236–12241 (2004).PubMed
PubMed Central
Google Scholar
Zheng, Z., Hua, R., Xu, G., Yang, H. & Shi, P. Gene losses may contribute to subterranean adaptations in naked mole-rat and blind mole-rat. BMC Biol. 20, 44 (2022).PubMed
PubMed Central
Google Scholar
Shams, I., Nevo, E. & Avivi, A. Ontogenetic expression of erythropoietin and hypoxia-inducible factor-1 alpha genes in subterranean blind mole rats. FASEB J. 19, 307–309 (2005).PubMed
Google Scholar
Avivi, A. et al. Neuroglobin, cytoglobin, and myoglobin contribute to hypoxia adaptation of the subterranean mole rat Spalax. Proc. Natl Acad. Sci. 107, 21570–21575 (2010).PubMed
PubMed Central
Google Scholar
Schülke, S. et al. Living with stress: regulation of antioxidant defense genes in the subterranean, hypoxia-tolerant mole rat, Spalax. Gene 500, 199–206 (2012).PubMed
Google Scholar
Avivi, A. et al. Increased blood vessel density provides the mole rat physiological tolerance to its hypoxic subterranean habitat. FASEB J. 19, 1314–1316 (2005).PubMed
Google Scholar
MacRae, S. L. et al. DNA repair in species with extreme lifespan differences. Aging 7, 1171–1182 (2015).PubMed
PubMed Central
Google Scholar
Kurz, S. et al. The anti-tumorigenic activity of A2M-A lesson from the naked mole-rat. PLoS One 12, (2017).Fang, X. et al. Adaptations to a Subterranean Environment and Longevity Revealed by the Analysis of Mole Rat Genomes. Cell Rep. 8, 1354–1364 (2014).PubMed
PubMed Central
Google Scholar
Dong, Q. et al. Transcriptome analysis of the response provided by Lasiopodomys mandarinus to severe hypoxia includes enhancing DNA repair and damage prevention. Front. Zool. 17, 9 (2020).PubMed
PubMed Central
Google Scholar
Wang, H. et al. Exposure to hypoxia causes stress erythropoiesis and downregulates immune response genes in spleen of mice. BMC Genomics 22, 413 (2021).PubMed
PubMed Central
Google Scholar
Heinze, I. et al. Species comparison of liver proteomes reveals links to naked mole-rat longevity and human aging. BMC Biol. 16, 82 (2018).PubMed
PubMed Central
Google Scholar
Lewis, K. N. et al. Unraveling the message: insights into comparative genomics of the naked mole-rat. Mamm. Genome 27, 259–278 (2016).PubMed
PubMed Central
Google Scholar
Vogel, C. et al. Sequence signatures and mRNA concentration can explain two-thirds of protein abundance variation in a human cell line. Mol. Syst. Biol. 6, 400 (2010).PubMed
PubMed Central
Google Scholar
Schmidt, H., Hangmann, J., Shams, I., Avivi, A. & Hankeln, T. Molecular evolution of antioxidant and hypoxia response in long-lived, cancer-resistant blind mole rats: The Nrf2-Keap1 pathway. Gene 577, 293–298 (2016).PubMed
Google Scholar
Lipton, P. Ischemic Cell Death in Brain Neurons. Physiol. Rev. 79, 1431–1568 (1999).PubMed
Google Scholar
Yuan, J. Neuroprotective strategies targeting apoptotic and necrotic cell death for stroke. Apoptosis 14, 469–477 (2009).PubMed
PubMed Central
Google Scholar
Coimbra-Costa, D., Alva, N., Duran, M., Carbonell, T. & Rama, R. Oxidative stress and apoptosis after acute respiratory hypoxia and reoxygenation in rat brain. Redox Biol. 12, 216–225 (2017).PubMed
PubMed Central
Google Scholar
Domankevich, V. et al. Adaptive patterns in the p53 protein sequence of the hypoxia- and cancer-tolerant blind mole rat Spalax. BMC Evol. Biol. 16, 177 (2016).PubMed
PubMed Central
Google Scholar
Band, M., Ashur-Fabian, O. & Avivi, A. The expression of p53-target genes in the hypoxia-tolerant subterranean mole-rat is hypoxia-dependent and similar to expression patterns in solid tumors. Cell Cycle 9, 3367–3372 (2010).
Google Scholar
Salah-Hussiesy, S., Mamchur, A., Odeh, A., Domankevich, V. & Shams, I. Reduced calcium influx in the hypoxia-tolerant Spalax: The role of the erythropoietin receptor. Cell Calcium 74, 123–130 (2018).PubMed
Google Scholar
Cabanac, A., Folkow, L. P. & Blix, A. S. Volume capacity and contraction control of the seal spleen. J. Appl Physiol. (1985) 82, 1989–1994 (1997).PubMed
Google Scholar
Eaton, L. & Pamenter, M. E. What to do with low O2: Redox adaptations in vertebrates native to hypoxic environments. Comp. Biochem Physiol. A Mol. Integr. Physiol. 271, 111259 (2022).PubMed
Google Scholar
Guzy, R. D. & Schumacker, P. T. Oxygen sensing by mitochondria at complex III: the paradox of increased reactive oxygen species during hypoxia. Exp. Physiol. 91, 807–819 (2006).PubMed
Google Scholar
Panchenko, M. V., Farber, H. W. & Korn, J. H. Induction of heme oxygenase-1 by hypoxia and free radicals in human dermal fibroblasts. Am. J. Physiol. Cell Physiol. 278, C92–C101 (2000).PubMed
Google Scholar
Jena, N. R. DNA damage by reactive species: Mechanisms, mutation and repair. J. Biosci. 37, 503–517 (2012).PubMed
Google Scholar
Renaudin, X. Reactive oxygen species and DNA damage response in cancer. Int. Rev. Cell Mol. Biol. 364, 139–161 (2021).PubMed
Google Scholar
Habibi, P. et al. Effect of heat stress on DNA damage: a systematic literature review. Int. J. Biometeorol. 66, 2147–2158 (2022).PubMed
Google Scholar
Pires, I. M. et al. Effects of acute versus chronic hypoxia on DNA damage responses and genomic instability. Cancer Res. 70, 925–935 (2010).PubMed
PubMed Central
Google Scholar
Gelot, C., Magdalou, I. & Lopez, B. S. Replication stress in Mammalian cells and its consequences for mitosis. Genes 6, 267–298 (2015).PubMed
PubMed Central
Google Scholar
Ramachandran, S. et al. Hypoxia-induced SETX links replication stress with the unfolded protein response. Nat. Commun. 12, 3686(2021).Peake, J. D. & Noguchi, E. Fanconi anemia: current insights regarding epidemiology, cancer, and DNA repair. Hum. Genet 141, 1811–1836 (2022).PubMed
Google Scholar
Domankevich, V., Eddini, H., Odeh, A. & Shams, I. Resistance to DNA damage and enhanced DNA repair capacity in the hypoxia-tolerant blind mole rat Spalax carmeli. J. Exp. Biol. 221, jeb174540 (2018).Torgovnick, A. & Schumacher, B. DNA repair mechanisms in cancer development and therapy. Front. Genet. 6, 157 (2015).PubMed
PubMed Central
Google Scholar
Schumacher, B., Pothof, J., Vijg, J. & Hoeijmakers, J. H. J. The central role of DNA damage in the ageing process. Nature 592, 695–703 (2021).PubMed
PubMed Central
Google Scholar
Garasto, S. et al. Thyroid hormones in extreme longevity. Mech. Ageing Dev. 165, 98–106 (2017).PubMed
Google Scholar
Papadopoli, D. et al. Mtor as a central regulator of lifespan and aging. F1000Res 8, F1000 Faculty Rev-998 (2019).Weichhart, T. mTOR as regulator of lifespan, aging, and cellular senescence: a mini-review. Gerontology 64, 127–134 (2018).PubMed
Google Scholar
Bjedov, I. & Rallis, C. The target of Rapamycin signalling pathway in ageing and lifespan regulation. Genes 11, 1–20 (2020).
Google Scholar
Blagosklonny, M. V. Rapamycin for longevity: opinion article. Aging 11, 8048–8067 (2019).PubMed
PubMed Central
Google Scholar
Wu, C. W. & Storey, K. B. mTOR signaling in metabolic stress adaptation. Biomolecules 11, 681 (2021).PubMed
PubMed Central
Google Scholar
Hua, H. et al. Targeting mTOR for cancer therapy. J. Hematol. Oncol. 12, 71 (2019).PubMed
PubMed Central
Google Scholar
Villafuerte, F. C. et al. Article 676782 (2021) Hypoxia and Inflammation: insights from high-altitude physiology. www.frontiersin.orgFront. Physiol. 1, 676782 (2021).Wanderer, A. A. Hypoxia and Inflammation. N. Engl. J. Med. 364, 1976–1977 (2011).PubMed
Google Scholar
Malik, A. et al. Transcriptome analysis of the spalax hypoxia survival response includes suppression of apoptosis and tight control of angiogenesis. BMC Genomics 13, 615 (2012).PubMed
PubMed Central
Google Scholar
Franceschi, C., Garagnani, P., Parini, P., Giuliani, C. & Santoro, A. Inflammaging: a new immune-metabolic viewpoint for age-related diseases. Nat. Rev. Endocrinol. 14, 576–590 (2018).PubMed
Google Scholar
Teissier, T., Boulanger, E. & Cox, L. S. Interconnections between inflammageing and immunosenescence during ageing. Cells 11, 359 (2022).PubMed
PubMed Central
Google Scholar
Odeh, A., Dronina, M., Domankevich, V., Shams, I. & Manov, I. Downregulation of the inflammatory network in senescent fibroblasts and aging tissues of the long-lived and cancer-resistant subterranean wild rodent, Spalax. Aging Cell 19, e13045 (2020).PubMed
Google Scholar
Toren, D. et al. Gray whale transcriptome reveals longevity adaptations associated with DNA repair and ubiquitination. Aging Cell 19, e13158 (2020).PubMed
PubMed Central
Google Scholar
Yu, C. et al. RNA sequencing reveals differential expression of mitochondrial and oxidation reduction genes in the long-lived naked mole-rat when compared to mice. PLoS One 6, e26729 (2011).PubMed
PubMed Central
Google Scholar
Lewis, K. N., Andziak, B., Yang, T. & Buffenstein, R. The naked mole-rat response to oxidative stress: just deal with it. Antioxid. Redox Signal 19, 1388–1399 (2013).PubMed
PubMed Central
Google Scholar
Thieme, R. et al. Analysis of Alpha-2 macroglobulin from the long-lived and cancer-resistant naked mole-rat and human plasma. PLoS One 10, e0130470 (2015).PubMed
PubMed Central
Google Scholar
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).PubMed
PubMed Central
Google Scholar
Dobin, A. & Gingeras, T. R. Optimizing RNA-Seq mapping with STAR. in. Methods Mol. Biol. 1415, 245–262 (2016).PubMed
Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).PubMed
PubMed Central
Google Scholar
Luo, W. & Brouwer, C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics 29, 1830–1831 (2013).PubMed
PubMed Central
Google Scholar
Chen, H. & Boutros, P. C. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinforma. 12, 35 (2011).
Google Scholar
Yu, Y. et al. A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages. Nat. Commun. 5, 3230 (2014).PubMed
Google Scholar
Krämer, A., Green, J., Pollard, J. & Tugendreich, S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics 30, 523–530 (2014).PubMed
Google Scholar
Download referencesAcknowledgementsWe thank Alexandra Hassemer and Simon Steines for contributing to the qPCR analyses. This project was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) GRK2526/1, project no. 407023052 to T.H.FundingOpen Access funding enabled and organized by Projekt DEAL.Author informationAuthor notesSören Lukas HellmannPresent address: Nucleic Acids Core Facility, Faculty of Biology, Johannes Gutenberg-University, Mainz, GermanyAnne BickerPresent address: Department of Medicine I, Johannes Gutenberg University Medical Center, Mainz, GermanyThese authors jointly supervised this work: Aaron Avivi, Imad Shams, Anne Bicker, Thomas Hankeln.Authors and AffiliationsMolecular Genetics & Genome Analysis, Institute of Organismic and Molecular Evolution, Faculty of Biology, Johannes Gutenberg-University, Mainz, GermanyGesa Poetzsch, Luca Jelacic, Leon Dammer, Sören Lukas Hellmann, Michelle Balling, Anne Bicker & Thomas HankelnComputational Biology and Data Mining Group, Institute of Organismic and Molecular Evolution, Johannes Gutenberg University, Mainz, GermanyMiguel Andrade-NavarroInstitute of Evolution, University of Haifa, Mount Carmel, Haifa, IsraelAaron Avivi & Imad ShamsDepartment of Evolutionary and Environmental Biology, Faculty of Natural Sciences, University of Haifa, Mount Carmel, Haifa, IsraelImad ShamsAuthorsGesa PoetzschView author publicationsYou can also search for this author in
PubMed Google ScholarLuca JelacicView author publicationsYou can also search for this author in
PubMed Google ScholarLeon DammerView author publicationsYou can also search for this author in
PubMed Google ScholarSören Lukas HellmannView author publicationsYou can also search for this author in
PubMed Google ScholarMichelle BallingView author publicationsYou can also search for this author in
PubMed Google ScholarMiguel Andrade-NavarroView author publicationsYou can also search for this author in
PubMed Google ScholarAaron AviviView author publicationsYou can also search for this author in
PubMed Google ScholarImad ShamsView author publicationsYou can also search for this author in
PubMed Google ScholarAnne BickerView author publicationsYou can also search for this author in
PubMed Google ScholarThomas HankelnView author publicationsYou can also search for this author in
PubMed Google ScholarContributionsT.H., I.S., A.B., and A.A. designed and supervised the study. G.P, L.J., and I.S. performed experiments. G.P., L.J., L.D., A.B, L.H., M.B., M. A.-N., I.S., A.A., and T.H. analyzed data. G.P., L.J., A.B., and T.H. drafted the manuscript. All authors have read and approved the final version of the manuscript.Corresponding authorCorrespondence to
Thomas Hankeln.Ethics declarations
Competing interests
The authors declare no competing interests.
Additional informationPublisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary informationSupplementary materialData Set 1Data Set 2Data Set 3Data Set 4Data Set 5Data Set 6Data Set 7Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissionsAbout this articleCite this articlePoetzsch, G., Jelacic, L., Dammer, L. et al. Adaptation of the Spalax galili transcriptome to hypoxia may underlie the complex phenotype featuring longevity and cancer resistance.
npj Aging 11, 16 (2025). https://doi.org/10.1038/s41514-025-00206-3Download citationReceived: 17 July 2024Accepted: 19 February 2025Published: 06 March 2025DOI: https://doi.org/10.1038/s41514-025-00206-3Share this articleAnyone you share the following link with will be able to read this content:Get shareable linkSorry, a shareable link is not currently available for this article.Copy to clipboard
Provided by the Springer Nature SharedIt content-sharing initiative