Most of our current understanding of the underlying mechanisms of aging comes from short-lived model species. It is, however, still largely unclear to what extent insights obtained from short-lived organisms can be transferred to long-lived species, such as humans (Parker et al., 2004; Keller and Jemielity, 2006). Comparative approaches, involving species with particularly long healthy lives and seeking the causative mechanisms that distinguish them from shorter-lived relatives, try to overcome this limitation (Austad, 2009). Many studies that involved organisms with particularly long lifespans, for example, queens in social hymenoptera, birds, bats, African mole-rats, and primates, have produced findings that were not always congruent with established aging theories (Keller and Jemielity, 2006; Austad, 2009; Salmon et al., 2009; Austad, 2011; Dammann, 2017; Bens et al., 2018). Species comparisons, however, also have their limitations. Many observed differences between species with differing lifespans are influenced by phylogenetic constraints, ecophysiological differences, or both, rather than being causal for the species-specific differences in aging and longevity.

Bimodal aging occurs naturally in the genus Fukomys from the rodent family Bathyergidae (African mole-rats). These animals live in families (often called colonies) of usually consisting of 9–16 individuals (Sichilima et al., 2008; Šklíba et al., 2012), although single families may occasionally grow considerably larger in some species (Jarvis and Bennett, 1993; Scharff et al., 2001). Regardless of group size, an established family typically consists of only one breeding pair (the founders of the family, often called king/queen) and their progeny from multiple litters (often called workers). Because of strict avoidance of incest (Burda, 1995), the progeny do not engage in sexual activity in the confines of their natal family, even after reaching sexual maturity. Hence, grown Fukomys families are characterized by a subdivision into breeders (the founder pair) and non-breeders (all other family members). Interestingly, breeders reach the age of 20 years or more in captivity, whereas non-breeders usually die before their tenth birth date (Figure 1A). This divergence of survival probabilities between breeders and non-breeders is found in all Fukomys species studied so far, irrespective of sex. Because no difference in diet or workload has been observed between breeders and non-breeders in captivity, status-specific changes of gene expression after the transition from non-breeder to breeder are considered the most likely explanation of the differing lifespans (Dammann and Burda, 2006; Dammann et al., 2011).

Figure 1 with 1 supplement see all

Download asset Open asset

In the wild, non-breeders must meet a member of another family by chance to ascend to breeder status; in captivity, the establishment of new breeder pairs is subject to human control. Allowing an animal to breed in captivity can be regarded as a simple experimental intervention that results in an extension of life expectancy of approximately 100%. This extension is far more than most experimental interventions in vertebrates can achieve, for example, by caloric restriction (e.g., Carmona and Michan, 2016) or diets containing resveratrol or rapamycin (e.g., Valenzano et al., 2006; Johnson et al., 2013). Furthermore, this relative lifespan extension starts from a non-breeder lifespan that is already more than twice as long as that of the mammalian model organisms most widely used in aging research, such as mice or rats.

Until now, relatively few studies have addressed the potential mechanisms behind this natural status dependency of aging in Fukomys sp. Contrary to the predictions of both the advanced glycation formation theory (Morimoto and Cuervo, 2009) and the oxidative stress theory of aging (Harman, 1956; Harman, 2001), markers of protein cross-linking and -oxidation were surprisingly higher in breeders of Ansell's mole-rats (Fukomys anselli) than in age-matched non-breeders (Dammann et al., 2012). On the other hand, in the Damaraland mole-rat (Fukomys damarensis), oxidative damage to proteins and lipids was significantly lower in breeding females than in their non-reproductive counterparts (Schmidt et al., 2014), a finding that is compatible with the oxidative stress theory of aging. In good agreement with their overall longevity irrespective of social status, Ansell's mole-rats produce less thyroxine (T4) and recruit smaller proportions of their total T4 resources into the active unbound form than do euthyroid mammals. Still, nonetheless the levels of unbound T4 (fT4) do not explain the intraspecific longevity differences between F. anselli breeders and non-breeders because the levels of this hormone did not differ between the two cohorts (Henning et al., 2014). Closely connected to the topic of this paper is the finding that non-breeding giant mole-rats (Fukomys mechowii) maintain fairly stable gene expression into relative old ages, quite in contrast to the shorter-lived Norway rat (Sahm et al., 2018a). It is, however, still unclear what happens on the gene expression level when an individual attains breeding status.

Interestingly, a recent study by Bens et al., 2018 found that the longest-lived rodent, the naked mole-rat Heterocephalus glaber, tended globally to show opposite changes in the transition from non-breeders to breeders compared to shorter-lived guinea pigs. Heterocephalus and Fukomys are similar in their mating and social behavior, but differences appear to exist regarding the effect of breeding on longevity: until very recently, naked mole-rat non-breeders have been reported to be as long-lived as breeders (Sherman and Jarvis, 2002; Buffenstein, 2008). In 2018, a lifespan advantage of breeders over non-breeders was reported in females (but not males), yet the divergence of the two groups appeared to be considerably smaller than in Fukomys (Ruby et al., 2018) and underlying data is being debated (Dammann et al., 2019). In summary, status-dependent aging is either absent in Heterocephalus or less pronounced than in Fukomys. Nevertheless, both Fukomys and naked mole-rats can be seen as counterexamples to the disposable soma theory of aging, which assumes that the aging process is linked to an evolutionary tradeoff in which the preservation of the organism is essentially sacrificed to reproductive success (Kirkwood, 1977).

To gain insights into the molecular basis of the bimodal aging pattern in Fukomys, it would be of course interesting to compare sufficiently large cohorts of old breeders and non-breeders both at similar chronological and biological ages. But as old animals, particularly of a long-lived species, are extremely valuable due to obvious logistical, timely, and financial reasons, such an approach was impossible for us. As an alternative, we decided to study young breeders and non-breeders that are much easier to obtain and hypothesized that the status change marks the beginning of a slowdown in the aging process. In this paper, we studied two Fukomys species (F. mechowii and Fukomys micklemi, Figure 1A) by comparing the gene expression profiles of breeders (n = 24) and age-matched non-breeders (n = 22) in 16 organs or their substructures (hereinafter referred to as tissues, Figure 1—figure supplement 1, Figure 1B). Our main aim was to identify genes and pathways whose transcript levels are linked to the status-dependent longevities and to relate these patterns to insights into aging research obtained in shorter-lived species.

We measured gene expression differences between breeders and non-breeders in two African mole-rat species, F. mechowii and F. micklemi. Altogether, we performed RNA-seq for 636 tissue samples covering 16 tissue types from both species, sexes, and reproductive states (breeders and non-breeders). Each of the four groups (male/female breeders/non-breeders) of each species consisted of 5–7 animals (see Supplementary file 1a–e for sample sizes, animal data, and pairing schemes). For each tissue, we conducted a multifactorial analysis of differentially expressed genes (DEGs): the analysis was based on the variables reproductive state, sex, and species. During this exercise, we focused on the differences between slower-aging breeders and faster-aging non-breeders. This approach increases our statistical power by giving us a fourfold increase of sample size in comparison to species- and sex-specific breeder vs. non-breeder analyses. At the same time, we can additionally reduce the number of false-positive DEGs by restricting the analysis to those breeding status-related genes that show the same direction in both sexes and species. We deliberately focused on those genes to concentrate our study on universal mechanisms that hold for both sexes and species. The number of samples collected and processed for this work required splitting them into different batches for sequencing. To analyze the degree this procedure potentially biased our results, we systematically investigated potential batch effects. Our results suggest that the sequencing strategy has little effect on the outcomes reported in the following (see Control analyses).

To globally quantify the transcriptomic differences between the reproductive states, we performed three analyses: clustering of the samples based on pairwise correlation, principal variant component analysis, and an overview of the number of DEGs between reproductive states in comparison to DEGs between species and sex. Clustering of the samples based on pairwise correlations showed a full separation of the two species at the highest cluster level (Figure 2—figure supplement 1). Below that level, an almost complete separation according to tissues was observed. Within the tissue clusters, the samples did not show a clear-cut separation between sex or breeder/non-breeder status. Accordingly, a principal variance component analysis showed that species, tissue, and the combination of both variables accounted for 98.4% of the total variance in the data set; individual differences explained 1.4% of the variance, and only 0.004% was explained by breeder/non-breeder status (Figure 2A). Regarding the numbers of DEGs, we found – unsurprisingly considering the aforementioned facts – by far the highest number of DEGs in the species comparison (Figure 2B). Although in almost every examined tissue the numbers of detected DEGs were also high between sexes, most tissues exhibited very few DEGs due to breeder/non-breeder status. Exceptions were liver, spleen, ovary, and, especially, tissues of the endocrine system (adrenal gland, pituitary gland, thyroid), in which the number of DEGs between breeders and non-breeders ranged from more than 60 to several thousand.

Figure 2 with 2 supplements see all

Download asset Open asset

Since the change from non-breeder to breeder status apparently marks the beginning of a slowdown in the aging process, we first wanted to find out whether and where there are intersections of reproductive status DEGs with those whose expression level is known to change during aging. Therefore, we determined overlaps by using the Digital Aging Atlas (DAA) – a database of genes that show aging-related changes in humans (Craig et al., 2015). Across species and sexes, significant overlaps (false discovery rate [FDR] < 0.05, Fisher’s exact test) with the DAA were found in three tissues: adrenal gland, ovary, and pituitary gland (FDR=0.005, each; Figure 3A). Among these three endocrine tissues, the DEGs of the ovaries overlapped significantly with those from adrenal (p=2.8*10⁻²⁷) and pituitary glands (p=0.005), but there was no significant overlap between the two glands (Figure 3A). Thus, together, we found indications for aging-relevant expression changes after the transition from non-breeders to breeders in three tissues of the endocrine system, which presumably affect separate aspects of aging in adrenal and pituitary glands.

Figure 3

Download asset Open asset

Moreover, we compared the DEGs with respect to the reproductive status that we identified in Fukomys with regard to their direction to transcript-level changes observed in similar experiments using naked mole-rats and guinea pigs (Bens et al., 2018). The direction of the status-dependent DEGs regulation in Fukomys, as found in this study, was significantly more often the same rather than opposite as in the naked mole-rat (females, 60%, p=5.2*10⁻⁵⁸; males, 62%, p=10⁻⁴⁴ for females, Figure 3B, Supplementary file 1f). In the guinea pig, on the contrary, the Fukomys reproductive status DEGs were significantly more often regulated in the opposite direction (females, 57%, p=9*10⁻²⁵; males, 59%, p=4.7*10⁻²³, Figure 3B, Supplementary file 1f). Thus, at the single-gene level, the expression changes linked to reproductive status may affect lifespan differently in long-lived African mole-rats than in shorter-lived guinea pigs.

Beyond the single-gene level, we aimed to identify metabolic pathways and biological functions whose gene expression significantly depends on reproductive status. For this, we used Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Kanehisa et al., 2017) and Molecular Signatures Database (MSigDB) hallmarks (Liberzon et al., 2015) as concise knowledge bases. As standard approach, we used an established method combining all p-values of genes in a given pathway in a threshold-free manner (Figure 4). The advantage of this approach is that it bundles the p-values from test results of individual gene expression differences at the level of pathways (see Materials and methods for details). In addition, we applied a second enrichment method that aggregates fold-changes instead (Figure 4—figure supplement 11; for a comparison of the results of both approaches, see Control analyses). Altogether, the gene expression of 55 KEGG pathways and 41 MSigDB hallmarks was significantly affected by reproductive status in at least one tissue (Figure 4—figure supplements 1 and 2). Because the individual interpretation of each of these pathways/hallmarks would go beyond the scope of this study, we focus here on those 14 pathways and 13 hallmarks that were significantly different between breeders and non-breeders (FDR < 0.1) in a global analysis across all tissues (Figure 4). Because many pathways are driven mainly by gene expression in subsets of tissues, we weighted gene-wise the differential expression signals from the various tissues by the respective expression levels in the tissues. For instance, the expression level of the growth hormone (GH) gene GH1 is known to be almost exclusively expressed in the pituitary glands. In our data set, the GH1 level of the pituitary gland accounted for 99.96% of the total GH1 across all tissues. Accordingly, in pathways that contain GH, our weighted cross-tissue differential expression signal for this gene is almost exclusively determined by the pituitary gland. On the contrary, a differential expression signal of this gene in another tissue with a very low fraction of the gene’s total expression would have almost no impact on the weighted cross-tissue level – even if that signal were very strong (see Materials and methods for details).

Figure 4 with 11 supplements see all

Download asset Open asset

We found strong indications for increase in the activity of certain anabolic functions in breeders: ribosomal protein expression (hsa03010 Ribosome, hsa03008 Ribosome biogenesis in eukaryotes, Figure 4A) was elevated in most tissues (and accordingly also in the weighted cross-tissue analysis). In MSigDB hallmarks, the strongest enrichment signal came from MYC targets (HALLMARK_MYC_TARGETS_V1), which can largely be considered a reflection of enhanced ribosomal protein expression and the fact that MYC is a basal transcription factor up-regulating genes involved in protein translation (Hofmann et al., 2015; Figure 4B). In functional correspondence, we observed an increase in the expression of mitochondrial respiratory chain components (hsa00190 Oxidative phosphorylation, HALLMARK_OXIDATIVE_PHOSPHORYLATION, Figure 4A, B). We also found strong indications for increased protein degradation (hsa03050 PROTEASOME, Figure 4A). This weighted cross-tissue signal was, in contrast to the situation regarding ribosomes, driven mainly by two tissue types: the adrenal gland and the gonads.

To examine whether the simultaneous up-regulation of the ribosome, proteasome, and oxidative phosphorylation is a coordinated regulation, we performed a weighted gene co-expression network analysis (WGCNA) (Langfelder and Horvath, 2008) from our gene count data and examined the connectivity between pairs of those KEGG pathways flagged in the weighted cross-tissue analysis. We found that the expression of ribosomal genes (hsa03010) was significantly linked to those of ribosome biogenesis (hsa03008, FDR = 4.59*10⁻³), oxidative phosphorylation (hsa00190, FDR = 4.05*10⁻⁴), and proteasome (hsa03050, FDR = 4.59*10⁻³), whereas no other examined pathway pair exhibited a significant connectivity (Figure 4—figure supplement 3). Interestingly, ribosome, proteasome, and oxidative phosphorylation pathways also shared other characteristics of their differential expression signals: subtle fold-changes, that is, up-regulation of 3–9% on average. Thus, statistically significant signals at the pathway level resulted from relatively small shifts in all genes of these pathways in a seemingly coordinated manner and across multiple tissues (Source data 1a). In addition, ribosome (hsa03010, in ovary), proteasome (hsa03050, in ovary and adrenal gland), and RNA-transport (hsa03013, in adrenal gland) are enriched in those Fukomys status-dependent DEGs that show, in a similar experimental setting (Bens et al., 2018), the same direction in both naked mole-rat sexes and the opposite direction in both guinea pig sexes (FDR < 0.05, Fisher’s exact test).

The myogenesis hallmark (Figure 4B) was also found to be up-regulated in breeders. Expectedly, this weighted cross-tissue result was driven mainly by differential expression signals from muscle tissue: muscle from all tissues exhibited the lowest p-value (Figure 4A), and 15 of 20 up-regulated genes that contributed most to the weighted cross-tissue differential myogenesis signal exhibited their highest expression in muscle. These genes were involved mainly in calcium transport or part of the fast-skeletal muscle-troponin complex (Source data 1a). A clear exception of this muscle-dominated expression is found in the gene that exhibited the highest relative contribution to the differential myogenesis signal, insulin-like growth factor 1 (IGF1). This gene was found to be expressed most strongly in ovary and liver and was strongly up-regulated in the breeders' ovaries and adrenal glands (Table 1). IGF1 codes for a well-known key regulator of anabolic effects such as cell proliferation, myogenesis, and protein synthesis (Schiaffino and Mammucari, 2011; Jung and Suh, 2014) and has a tight functional relation to GH (gene: GH1) another key anabolic regulator upstream of IGF1; together, these factors form the so-called GH/IGF1 axis (Cannata et al., 2010; Junnila et al., 2013; Bodart et al., 2015; Raisingani et al., 2017; Carotti et al., 2018; Lozier et al., 2018). Also, GH1 was strongly up-regulated in breeders in its known principal place of synthesis, the pituitary gland (Table 1).

With xenobiotic metabolism and TNF-α signaling, two defense hallmarks were also found to be up-regulated in breeders by the weighted cross-tissue analysis (HALLMARK_XENOBIOTIC_METABLISM and HALLMARK_TNFA_SIGNALING_VIA_NFKB, Figure 4B). The up-regulation of the reactive oxygen species (ROS) hallmark comprising genes coding for proteins that detoxify ROS (HALLMARK_REACTIVE_OXYGEN_SPECIES_PATHWAY, Figure 4B) falls into a similar category.

Another interesting aspect that was found to be significantly altered in breeders is steroid hormone biosynthesis (hsa00140 Steroid hormone biosynthesis, Figure 4A). In this case, both up- and down-regulated genes were in the pathway, and their absolute fold-changes were roughly balanced. Steroid hormones, on the one hand, comprise sex steroids – because these hormones are important players in sexual reproduction, such differences should be expected given the experimental setup. On the other hand, the class of steroid hormones – corticosteroids – has regulatory functions in metabolism, growth, and the cardiovascular system, as well as in the calibration of the immune system and response to stress (Liu et al., 2013). The most influential contributor to the differential pathway signal by far on the weighted cross-tissue level was CYP11A1, which codes for the (single) enzyme that converts cholesterol to pregnenolone. This is the first and rate-limiting step in steroid hormone synthesis (Miller and Auchus, 2011). Because CYP11A1 was found to be up-regulated in breeders in its main places of synthesis – the gonads and the adrenal gland (Table 1) – it can be assumed that the total output of steroid hormone biosynthesis in breeders is increased. The pattern of up- and down-regulation on the KEGG pathway, however, suggests that sex steroids especially are produced at a higher rate in breeders, whereas circulating levels of glucocorticoids – such as cortisol – should be lower than in non-breeders (Figure 4—figure supplement 4; see also our discussion on ACTH-R below).

Finally, several pathways flagged by the weighted cross-tissue analysis seem to be derivatives of the abovementioned differentially expressed pathways instead of representing altered functions on their own. For example, Huntington’s (hsa05016), Parkinson’s (hsa05012), and Alzheimer’s (hsa05010) diseases could, in principle, be interpreted as highly relevant for aging and lifespan. A closer inspection of these pathways reveals, however, that the genes of the mitochondrial respiratory chain – which is the core of the oxidative phosphorylation pathway – are in all three cases the main contributors to the respective differential expression signals (Figure 4—figure supplements 5–7, Source data 1a). Similarly, we see in the case of fat digestion (hsa04975 Fat digestion and absorption) that two of the three largest contributors to the differential expression signal of that pathway – ABCG8 and SCARB1 – are directly involved in the transport of cholesterol (Liu et al., 2008; Wang et al., 2015). Therefore, it seems likely that this signal is an expression of the altered steroid hormone biosynthesis rather than indicating altered fat digestion.

Interestingly, three genes, which we had already mentioned as potential regulators during the pathway analysis, also appeared among the 10 most clearly altered genes on the weighted cross-tissue level: GH1, IGF1, and CYP11A1 (Table 1). The top 2 among these 10 are sulfotransferase family 2A member 1 (SULT2A1) and melanocortin two receptor (MC2R). SULT2A1 is the main catalyzer of the sulfonation of the steroid hormone dehydroepiandrosterone (DHEA) to its non-active form DHEA-S (Hammer et al., 2005). DHEA has repeatedly been proposed to be an ‘anti-aging hormone’ because its levels are negatively associated with chronological aging (Baulieu, 1996; Celec and Stárka, 2003; Rutkowski et al., 2014). We found that SULT2A1 is strongly down-regulated in breeders’ liver, which is also the main location of its enzymatic action. The second candidate, MC2R, encodes the adrenocorticotropin hormone (ACTH) receptor, which is the main inducer of glucocorticoid synthesis and a crucial component of the hypothalamic–pituitary–adrenal (HPA) axis (Walker et al., 2015). In humans and many other mammals, prolonged glucocorticoid excess leads to Cushing’s syndrome. Affected individuals exhibit muscle weakness, immune suppression, impairment of the GH/IGF1 axis, higher risk of diabetes, cardiovascular disease (hypertension), osteoporosis, decreased fertility, depression, and weight gain (Chabre, 2014; Ferraù and Korbonits, 2015). The large overlap of these symptoms with those of aging could explain to some extent that Cushing’s syndrome patients exhibit considerably higher mortality rates (Etxabe and Vazquez, 1994). We hypothesized that the increased expression of the ACTH receptor in Fukomys non-breeders can cause similar expression patterns and consequences (Figure 5).

Figure 5

Download asset Open asset

We subjected this hypothesis (Figure 5) to a first test by matching the global fold-changes of the breeder/non-breeder comparison with those of a gene expression comparison of human controls/Cushing patients (Hochberg et al., 2015) and found a significant correlation (R = 0.11, p=6.6*10⁻³⁶). Since the target tissue of Hochberg et al., subcutaneous adipose tissue was unfortunately not included in our study, we used our cross-tissue data for comparison. If we use instead our data from skin as adjacent tissue, a similar correlation results (R = 0.11, p=1.1*10⁻⁴⁰).

We tested the hypothesis further by checking five of its key predictions. Altered MC2R expression (Figure 6A) coincides with higher cortisol levels in hair samples from non-breeding F. mechowii than in those from breeders of the same species (Begall et al., 2021). Furthermore, glucocorticoids such as cortisol exert their effect by binding to the glucocorticoid receptor that, in turn, acts as a transcription factor for many genes (Gjerstad et al., 2018). We tested whether the expression of targets of the glucocorticoid receptor (NR3C1) was significantly altered throughout our data using two gene lists (Phuc Le et al., 2005): about 300 direct target genes of the receptor that were identified by chromatin immunoprecipitation (i) and about 1300 genes that were found to be differentially expressed depending on the presence or absence of exogenous glucocorticoid (ii). Both gene lists were found to be significantly affected by differential expression at the weighted cross-tissue level (i, p=0.001; ii, p<10⁻⁹) as well as in five (i) and eight (ii) single tissues (Supplementary file 1g, h). In line with our hypothesis, we observed that the weight gain in non-breeders was, on average, twice as strong compared to the weight gain in breeders during the experiment (p=7.49*10⁻³, type II ANOVA, Figure 6B). In addition, we found a subtle but significant influence of reproductive status on the density of the vertebrae: the vertebrae of breeders were slightly denser than those of age-matched non-breeders (p=0.03 for vertebra T12 only, and p=0.01 across all examined vertebrae L1, L2, and T12; ANOVA, Figure 6C).

Figure 6 with 1 supplement see all

Download asset Open asset

The vast lifespan differences between Fukomys breeders and non-breeders are, according to our RNA-seq data, associated with only subtle global pattern shifts in transcript levels. Concerning the tested explanatory variables (Fukomys species, sex, breeding status), we found by far the highest number of DEGs at the level of the species comparison. Although the number of DEGs between the sexes was comparably high in almost each examined tissue, only very few DEGs were found comparing breeders and non-breeders. One exception is the ovary, whose high number of DEGs corresponds well with the disparity in reproductive activity. Other exceptions are liver, spleen, and, especially, the tissues of the endocrine system (adrenal gland, pituitary gland, thyroid), in which the number of DEGs between breeders and non-breeders ranged from more than 100 to more than 2500.

Changes in the gene expression of the endocrine system are well known to play an important role both in sexual maturation and in aging and the development of aging-associated diseases, for example, diabetes and cardiovascular diseases (Tatar et al., 2003; Chahal and Drake, 2007; Jones and Boelaert, 2015). This finding fits well with the observation of substantial changes in the endocrine system after the transition from non-breeders to breeders in the related naked mole-rat – one of the key results of a recent study in that species (Bens et al., 2018). The dominance of differential expression in endocrine tissue is also plausible insofar as these tissues exert a strong control function for other tissues via hormone release.

Steroid hormone biosynthesis exhibits a bipartite pattern in breeders, with up-regulated sex steroid genes and a simultaneous down-regulation of corticosteroid synthesis genes. The former could be expected as a consequence of sexual activity in breeders. For aging, it is, however, interesting that SULT2A1, a gene that codes for the specialized sulfotransferase converting the sex steroid DHEA to its non-active form DHEA-S, was found to be heavily down-regulated in breeders. DHEA is the most abundant steroid hormone; it serves as a precursor for sex steroid biosynthesis but also has various metabolic functions on its own (Allolio and Arlt, 2002; Webb et al., 2006). DHEA levels decrease continuously during the human aging process to an extent that favors it as an aging biomarker (Maggio, 2007; Traish et al., 2011). As a result, an aggressive advertising of DHEA as ‘anti-aging hormone’ in the form of dietary supplements could be observed in recent years (Baulieu, 1996; Celec and Stárka, 2003; Webb et al., 2006; Rutkowski et al., 2014). Despite conflicting experimental data from various animal studies and clinical trials, a positive effect of DHEA on human health is frequently considered to be likely in the literature. Large-scale and long-term studies, however, are still pending (Traish et al., 2011; Rutkowski et al., 2014; Samaras et al., 2015). Interestingly, DHEA is not only described as an aging marker but also as a marker for clinically relevant glucocorticoid excess (Burkhardt et al., 2013).

Regarding corticosteroid synthesis, we hypothesize that the regulation of adrenal gland MC2R, coding for the ACTH receptor, causes a considerable proportion of the overall observed expression patterns and of the lifespan extension in breeders (Figure 5). As a critical component of the HPA axis, ACTH is a stress hormone that is produced by the pituitary gland and transported by the blood to the adrenal cortex, where it binds to the ACTH receptor (Fridmanis et al., 2017). Subsequently, the ACTH receptor induces the synthesis of glucocorticoids (Walker et al., 2015), for example, cortisol, which in turn cause immunosuppressive and various metabolic effects throughout the organism (Becker, 2013). In many mammals (such as humans, dogs, and guinea pigs) glucocorticoid excess disorder, Cushing’s syndrome, is caused by overproduction of the hormone ACTH. Our results for Fukomys suggest that the increased levels of the ACTH receptor may lead to symptoms and expression patterns that could resemble some of molecular and phenotypic aspects of this pathological condition (Figure 5).

We looked for evidence of a regulation upstream of the ACTH receptor. Given the known positive feedback loop between the ACTH receptor and its own ligand, decreased ACTH synthesis in breeders would have been an obvious explanation (Imai et al., 2001). We found that the expression of the ACTH polypeptide precursor gene (POMC) is unchanged. However, also other post-transcriptional or post-translational mechanisms such as cleaving may influence the ACTH levels in breeders and non-breeders. Of those genes known to be involved in the regulation of MC2R (Beuschlein et al., 2001; Lin et al., 2007), we found that only PRKAR1B was differentially regulated in the adrenal gland. However, this gene codes for only one of seven subunits of the involved protein kinase A. Alternative explanations could be that epigenetic modifications, other still unidentified regulators of the transcript level, or both are responsible for the differential expression of MC2R.

While the opportunities to pursue the hypothesis upstream of the ACTH receptor were limited with our data set, many of the predicted downstream effects (Figure 5) could be verified. One of the clearest symptoms of Cushing's syndrome known in many animal species is an increase in body fat and a consequent significant weight gain (Ferraù and Korbonits, 2015). Our recording of body weights at the beginning and end of the experiment showed that non-breeders gained twice as much weight on average as breeders (p=7.49*10⁻³, type II ANOVA, Figure 6B). However, it should be noted that our experimental design does not allow to distinguish between increased body fat content and other possible reasons (i.e., overall growth, muscle gain) for these differences in weight gain. Another known complication associated with glucocorticoid excess is osteoporosis (Lodish et al., 2010). A computed tomography measurement of three vertebrae (L1, L2, and T12), in accordance with this, showed significantly lower bone densities of non-breeders (p=0.01, ANOVA, Figure 6C). The absolute differences, however, are probably too small to assume that the non-breeders are already physically impaired by that. It seems possible that a longer experiment duration than the 1–2 years we conducted would have aggravated this difference, but this remains speculative. Directly from our gene expression analysis we know that the GH/IGF1 axis of the non-breeders is impaired, which is another well-known symptom of chronic glucocorticoid excess (Ferraù and Korbonits, 2015).

In addition, we have systematically compared our observed and expected – based on the Cushing’s syndrome hypothesis – expression patterns. A key element in this context is the glucocorticoid receptor (NR3C1), which, when activated by glucocorticoids such as cortisol, is considered to mediate their effects, including those of Cushing's syndrome in the case of chronic activation (Gjerstad et al., 2018). According to our hypothesis, the known, experimentally determined NR3C1 target genes (Phuc Le et al., 2005) should therefore show a clear pattern of expression differences between breeders and non-breeders. This was the case for both the 1300 indirect (p<10⁻⁹) and the 300 direct target genes (p=0.001). Furthermore, according to our hypothesis, the expression changes in the comparison of control individuals against patients with Cushing's syndrome should correspond to those of the comparison of breeders vs. non-breeders. A global comparison of fold-changes from human subcutaneous adipose tissue (Hochberg et al., 2015) against our cross-tissue data showed a significant correlation (R = 0.11, p=6.6*10⁻³⁶). A very similar result was obtained by comparing the human data with our data from skin as the tissue closest to subcutaneous adipose tissue that was examined in our study (R = 0.11, p=1.1*10⁻⁴⁰). Despite the clear significance of these correlations, the reported values of the correlation coefficients could be considered to be only modest. It should, however, be taken into account that two evolutionarily widely separated species and not exactly matching tissues were compared.

For further validation, we applied the combination of limma (Smyth, 2004) and voom (Law et al., 2014) as an alternative DEG detection method to DESeq2 to account for potential dispersion problems in our data. All single DEGs mentioned in the paper were confirmed by this validation. This specifically includes all genes supporting our HPA hypothesis, for example, SULT2A1 (liver), MC2R (adrenal gland), CYP11A1 (adrenal gland and ovary), GH (pituitary gland), and IGF1 (ovary). In a similar way, we checked the main findings obtained with our standard bioinformatic p-value-based enrichment approach (Fridley et al., 2010; Evangelou et al., 2012; Poole et al., 2016), using a second method that aggregates fold-changes instead. In particular, this step confirmed our initial finding of an enrichment of genes involved in the steroid hormone biosynthesis and ribosome metabolism as a central part of up-regulated anabolism in breeders, consistent with our HPA hypothesis (Figure 4—figure supplement 11) (see Control analyses).

In view of the different evidence presented, we consider it justified to assume that the transcriptional, metabolic and physiological condition of Fukomys non-breeders resembles that of chronic glucocorticoid excess described from other species. Since this condition shortens life in other species (Etxabe and Vazquez, 1994) and several of its symptoms also occur during aging, we hypothesize that it may shorten the lifespans in Fukomys non-breeders as well. The ultimate proof and the quantification of its impact on Fukomys lifespan, however, will be difficult in such a long-lived species, which cannot yet be genetically manipulated. Nevertheless, a recent study confirmed one of the most important predictions of the hypothesis: hair cortisol levels of non-breeders are on average about twice as high as those of breeders. However, a subgroup of non-breeders, defined by social environmental factor (one or both parents had died), exhibited cortisol levels nearly as low as in breeders (Begall et al., 2021). The drop of cortisol levels in orphaned non-breeders is accompanied by a substantial drop of mortality rates in this subgroup approaching the values of breeders (Figure 6—figure supplement 1). This means that across status groups (breeders, orphaned non-breeders, and non-breeders from intact families) mortality rates are clearly correlated with long-term cortisol levels. Taken together, all evidence suggests that the presence of breeders causes distress in non-breeders, resulting in constantly activating the HPA axis and contributing to reduced life expectancy.

In the circulatory system, represented by blood and spleen, down-regulation of coagulation factors was observed in breeders. Because coagulation factors are known to be up-regulated during aging in humans mice, rats, and even fish, this can be interpreted as a sign of a more juvenile breeders’ transcriptome (Ochi et al., 2016; Benayoun et al., 2019). Furthermore, the activity of coagulation factors is associated with a higher risk of coronary heart disease (Lowe and Rumley, 2014). However, we found no obvious histopathological lesions in the hearts or other organs (spleen, kidney, liver, lung) of non-breeders. Therefore, if down-regulation of coagulation attenuates the aging process in Fukomys, it seems to exert its effect only, if at all, at a latter age not investigated here.

Two defense mechanisms were also found to be up-regulated in breeders: xenobiotic metabolism and TNF-α signaling. Increased TNF-α signaling often leads to the induction of apoptosis (Annibaldi and Meier, 2018). In line with this, we found that apoptosis and P53 signaling were also up-regulated in breeders. Apoptosis is considered an important anticancer mechanism (Baig et al., 2016; Pistritto et al., 2016). We hypothesize this to compensate cancerogenic effects of the anabolic alterations described above (especially the up-regulation of the GH/IGF1 axis). In line with this, our more than 30 years’ breeding history with several Fukomys species in Germany and the Czech Republic suggests that breeders are as ‘cancer-proof’ as non-breeders despite their much longer lifespan (own unpublished data and R. Šumbera, personal communication).

Many anabolic pathways are up-regulated in breeders across tissues: protein biosynthesis, myogenesis, and the GH/IGF1 axis. In line with this finding and with the fact that protein synthesis consumes 30–40% of a cell’s ATP budget (Hands et al., 2009), we observed increased expression of mitochondrial respiratory chain components, implying an increase in the capacity for cellular ATP production. On the other hand, protein degradation and clearance by the proteasome are also up-regulated in breeders. The fact that the expression of proteasomal genes is significantly linked to the genes of ribosome biogenesis and oxidative phosphorylation indicates that those processes influence, or even trigger, each other and hence are regulated in a coordinated manner in Fukomys mole-rats.

An important positive regulator of ribosome biogenesis and protein synthesis is mTOR (Johnson et al., 2013). Based on this, one could expect an up-regulation of the corresponding gene in our scenario. At the same time, however, the inhibition of mTOR is one of the best-documented life-prolonging interventions from invertebrates to mammals (Table 2). Given these conflicting premises, we find mTOR in almost all Fukomys tissues as not significantly altered; the exception is the adrenal gland, where mTOR is significantly down-regulated in breeders. It is obvious that, during evolution of Fukomys mole-rats, both the extension of breeders’ lifespan and an increase in their anabolic processes provided fitness benefits. Consequently, in the underlying organismic framework mTOR may have acquired expression patterns and functions, different from those in organisms studied before. Thus, in the future it may be worth to study Fukomys mTOR biology in more detail, particularly in respect to its potential role in naturally evolved ways towards lifespan extension and healthy aging.

The results of differential expression of anabolic components such as the GH/IGF1 axis are surprising. They fall within a debate in aging research that has been highly controversial over time: based on the well-known fact that the expression and secretion of GH and IGF1 decline with age in humans and other mammals (Bartke, 2019), Rudman et al., 1990 administered synthetic GH to elderly subjects in 1990, thereby reversing a number of aging-associated effects such as expansion of adipose mass. This led to GH being celebrated as an anti-aging drug (Junnila et al., 2013), including dubious commercial offers. Today's aging research, on the contrary, strongly assumes that the enhanced activity of the GH/IGF1 axis accelerates aging and that its suppression could extend lifespan even in humans (López-Otín et al., 2013; Longo et al., 2015; Pitt and Kaeberlein, 2015). In addition to several studies of synthetic GH in humans yielding less convincing results than those of Rudman et al., the main reasons for this turn are the results of studies on short-lived model organisms. From worms to mice, the impairment of the GH/IGF1 axis by genetic intervention consistently led to longer lifespans (Table 2), for example, the up-regulation of Klotho – an IGF1 inhibitor – extended the mouse lifespan by as much as 30% (Kurosu et al., 2005). As with the impairment of the GH/IGF1 axis, reducing of protein synthesis by decreasing the expression of MYC, a basal transcription factor, extended the mouse lifespan by as much as 20% (Hofmann et al., 2015), whereas the impairment of the respiratory chain by rotenone resulted in prolongation of the killifish lifespan by 15% (Baumgart et al., 2016; Table 2).

Therefore, it is astonishing that massive up-regulation of these anabolic key components accompanies a lifespan extension of approximately 100% in long-lived mammals and potentially even contributes to it. Several points could help to resolve this apparent contradiction: first, the up-regulation of anabolic pathways and key genes is at least partially accompanied by the regulation of other mechanisms that could plausibly compensate for deleterious effects. For example, it is, widely assumed that the negative impact of enhanced protein synthesis on lifespan is to a large extent caused by the accumulation of damaged or misfolded proteins, which is also known to contribute to aging-associated neurodegenerative diseases (Hipkiss, 2007; Saez and Vilchez, 2014; Carmona and Michan, 2016). Up-regulation of the proteasome, as we observed in breeders in a weighted cross-tissue approach and especially in endocrine tissues, is known to counteract these effects by clearing damaged proteins, leading to lifespan extension in worm and fly (Saez and Vilchez, 2014). Enhanced proteasome activity has also been linked with higher longevity of the naked mole-rat compared to the laboratory mouse (Pérez et al., 2009; Rodriguez et al., 2012) and in comparative approaches across several mammalian lineages (Pride et al., 2015). We hypothesize that the simultaneously high anabolic synthesis and catabolic degradation of proteins will lead to a higher protein turnover rate in breeders and, accompanied with that, a reduced accumulation of damaged and misfolded proteins. Similarly, it seems plausible that the up-regulation of the mitochondrial respiratory chain (oxidative phosphorylation) in breeders is compensated for by simultaneous up-regulation of the reactive oxygen hallmark: the mitochondrial respiratory chain is the main source of cellular ROS that can damage DNA, proteins, and other cellular components (Balaban et al., 2005; Starkov, 2008); the reactive oxygen hallmark consists by definition of genes that are known to be up-regulated in response to ROS treatments. Unsurprisingly, at least 25% of these genes code for typical antioxidant enzymes such as thioredoxin, superoxide dismutase, peroxiredoxin, or catalase that can detoxify ROS. Furthermore, the known cancer-promoting effects of an enhanced GH/IGF1 axis (Junnila et al., 2013) could, to some extent, be compensated for by up-regulation of apoptosis and p53 signaling because these are the major mechanisms of cancer suppression (Bieging et al., 2014). More generally, potential lifespan-extending effects of moderate up-regulation of both the GH/IGF1 axis and ROS production can also be viewed in the light of the hormesis hypothesis (Ristow, 2014), which postulates that mild stressors can induce overall beneficial adaptive stress responses. In line with these arguments, we found higher resting metabolic rates in breeders compared to non-breeders in F. anselli (Schielke et al., 2017), a species closely related to F. micklemi.

A second point that could help to resolve this apparent contradiction concerns the time of intervention. The transition from non-breeder to breeder takes place in adulthood when by far the largest portion of the growth process has already been completed. In contrast, genetic interventions aimed at prolonging the lifespan by inhibiting the GH/IGF1 axis (Table 2) often affect the organisms throughout their entire lifespan, including infancy and youth. However, there are some exceptions: for example, late treatments with rapamycin prolong the life of mice almost as much as those that are started early (Johnson et al., 2013). Therefore, it is still under debate whether the up-regulation of translation and anabolic processes by the GH/IGF1 axis independently enhances growth and aging or enhances aging only secondarily as a consequence of accelerated growth (Bartke, 2017). Our results are an argument for the latter.

A third point is the question of the transferability of knowledge obtained in one species to other species. Most insights into current aging research originate from very short-lived model organisms (Table 2). It is clear, on the other hand, that the observed effects of lifespan-prolonging interventions listed in Table 2 are by far the smallest in the model organisms with the relatively longest lifespans: mice and rats. For example, knockout or knockdown of IGF1R in worms that normally live only a few days can result in a lifespan extensions of more than 500% (Houthoofd et al., 2004). For mice that have a maximum lifespan of 4 years (Tacutu et al., 2013), comparable interventions only lead to an extension of the lifespan by about 25% (Holzenberger et al., 2003).It seems reasonable to hypothesize that the effect size of such interventions, regardless of their direction, would diminish further in mole-rats or humans that can live more than 20 or 100 years, respectively, especially if the magnitude of gene expression change is small as in the comparison of breeders and non-breeders. Therefore, it may also be possible that those gene expression patterns caused by our lifespan-extending intervention in Fukomys mole-rats highlight particularly the differences in aging mechanisms between short-lived and long-lived species.

As a fourth perspective, one could interpret the down-regulation of the GH/IGF1 axis in non-breeders as a by-product of the apparent up-regulation of the HPA axis in non-breeders, which may well be adaptive in itself. In the wild, non-breeding Fukomys mole-rats can maximize their inclusive fitness by either supporting their kin in their natal family or by founding a new family elsewhere. It has therefore been suggested that the shorter lives of non-breeders could be adaptive on the ultimate level if longevity were traded against some other fitness traits, such as competitiveness, to defend colonies against intruders or to enhance their chances for successful dispersal (Dammann and Burda, 2006). A constitutively more activated HPA stress axis is expected to offer advantages for both family defense and dispersal but it carries the risk of status-specific aging symptoms, such as muscle weakness, lower GH/IGF1 axis activity, and lower bone density in the long run (Ferraù and Korbonits, 2015). This effect may become even more pronounced under laboratory conditions where grown non-breeders cannot decide to disperse even if they would like to.

However, even the down-regulation of the GH/IGF-axis may be adaptive in itself for non-breeders if it has the potential to protect them from further damage. Note that for today's conventional view that stronger activation of the GH/IGF1 axis accelerates aging (Table 2) it is generally challenging that decreasing activity is well documented to correlate with chronological age in a wide range of mammals, including mice, rats, dogs, and humans (Bartke, 2019). Also, decreasing activity correlates, under pathological conditions such as Cushing’s syndrome, with many symptoms akin to aging (Figure 5). It has been suggested that one solution to this apparent contradiction may be that the GH/IGF1 axis is adaptively down-regulated in aging organisms as a reaction to already accumulated aging-related symptoms so that additional damage can be avoided (Berryman et al., 2008; Milman et al., 2016).

Finally, recent findings of positively selected genes in African mole-rats (family Bathyergidae, containing also Fukomys and Heterocephalus) could partly explain some of the surprising results. It is striking that translation and oxidative phosphorylation were among the strongest differentially expressed molecular processes concerning the breeding status. Earlier, these processes were also reported to be the most affected by positive selection in the phylogeny of African mole-rats (Sahm et al., 2018b). Furthermore, IGF1 was one of 13 genes that were found to be under positive selection in the last common ancestor of the mole-rats. This may indicate that the corresponding mechanisms were evolutionarily adapted to be less detrimental and make their up-regulation more compatible with a long lifespan. Since the mere fact of positive selection does not permit to draw conclusions about the direction of the mechanistic effect, this hypothesis, however, needs to remain speculative.

Read datasets generated during the current study are available in the European Nucleotide Archive, study ID: PRJEB29798.

1. 1. Sahm A
   2. Matthias P
   3. Szafranski K
   4. Dammann P

   (2018) European Nucleotide Archive

   ID PRJEB29798. Transcriptome signatures of fast vs. slow aging in Fukomys mole-rat breeders vs. non-breeders.

   https://www.ebi.ac.uk/ena/browser/view/PRJEB29798

1. 1. Allolio B
   2. Arlt W

   (2002) DHEA treatment: myth or reality?

   Trends in Endocrinology & Metabolism 13:288–294.

   https://doi.org/10.1016/S1043-2760(02)00617-3

   • PubMed
   • Google Scholar
2. Software

   1. Andrews S

   (2020) FastQC a Quality Control Tool for High Throughput Sequence Data

   FastQC a Quality Control Tool for High Throughput Sequence Data.

   https://www.bioinformatics.babraham.ac.uk/projects/fastqc/
3. 1. Annibaldi A
   2. Meier P

   (2018) Checkpoints in TNF-Induced cell death: implications in inflammation and Cancer

   Trends in Molecular Medicine 24:49–65.

   https://doi.org/10.1016/j.molmed.2017.11.002

   • PubMed
   • Google Scholar
4. 1. Austad SN

   (2009) Comparative biology of aging

   The Journals of Gerontology Series A: Biological Sciences and Medical Sciences 64:199–201.

   https://doi.org/10.1093/gerona/gln060

   • PubMed
   • Google Scholar
5. 1. Austad SN

   (2011) Candidate bird species for use in aging research

   ILAR Journal 52:89–96.

   https://doi.org/10.1093/ilar.52.1.89

   • PubMed
   • Google Scholar
6. 1. Baig S
   2. Seevasant I
   3. Mohamad J
   4. Mukheem A
   5. Huri HZ
   6. Kamarul T

   (2016) Potential of apoptotic pathway-targeted Cancer therapeutic research: where do we stand?

   Cell Death & Disease 7:e2058.

   https://doi.org/10.1038/cddis.2015.275

   • PubMed
   • Google Scholar
7. 1. Balaban RS
   2. Nemoto S
   3. Finkel T

   (2005) Mitochondria, oxidants, and aging

   Cell 120:483–495.

   https://doi.org/10.1016/j.cell.2005.02.001

   • PubMed
   • Google Scholar
8. 1. Bartke A

   (2017) Somatic growth, aging, and longevity

   Npj Aging and Mechanisms of Disease 3:14.

   https://doi.org/10.1038/s41514-017-0014-y

   • Google Scholar
9. 1. Bartke A

   (2019) Growth hormone and aging: updated review

   The World Journal of Men's Health 37:19–30.

   https://doi.org/10.5534/wjmh.180018

   • PubMed
   • Google Scholar
10. 1. Baulieu EE

    (1996) Dehydroepiandrosterone (DHEA): a fountain of youth?

    The Journal of Clinical Endocrinology & Metabolism 81:3147–3151.

    https://doi.org/10.1210/jcem.81.9.8784058

    • PubMed
    • Google Scholar
11. 1. Baumgart M
    2. Priebe S
    3. Groth M
    4. Hartmann N
    5. Menzel U
    6. Pandolfini L
    7. Koch P
    8. Felder M
    9. Ristow M
    10. Englert C
    11. Guthke R
    12. Platzer M
    13. Cellerino A

    (2016) Longitudinal RNA-Seq analysis of vertebrate aging identifies mitochondrial complex I as a Small-Molecule-Sensitive modifier of lifespan

    Cell Systems 2:122–132.

    https://doi.org/10.1016/j.cels.2016.01.014

    • PubMed
    • Google Scholar
12. 1. Becker DE

    (2013) Basic and clinical pharmacology of glucocorticosteroids

    Anesthesia Progress 60:25–32.

    https://doi.org/10.2344/0003-3006-60.1.25

    • PubMed
    • Google Scholar
13. 1. Begall S
    2. Nappe R
    3. Ll. H
    4. Burda H
    5. Sahm A
    6. Szafranski K
    7. Dammann P
    8. Henning Y

    (2021) Life expectancy, family constellation, and stress in giant mole-rats (Fukomys mechowii)

    Philosophical Transactions of the Royal Society B 376:20200207.

    https://doi.org/10.1098/rstb.2020.0207

    • Google Scholar
14. 1. Benayoun BA
    2. Pollina EA
    3. Singh PP
    4. Mahmoudi S
    5. Harel I
    6. Casey KM
    7. Dulken BW
    8. Kundaje A
    9. Brunet A

    (2019) Remodeling of epigenome and transcriptome landscapes with aging in mice reveals widespread induction of inflammatory responses

    Genome Research 29:697–709.

    https://doi.org/10.1101/gr.240093.118

    • PubMed
    • Google Scholar
15. 1. Benjamini Y
    2. Hochberg Y

    (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing

    Journal of the Royal Statistical Society: Series B 57:289–300.

    https://doi.org/10.1111/j.2517-6161.1995.tb02031.x

    • Google Scholar
16. 1. Bens M
    2. Szafranski K
    3. Holtze S
    4. Sahm A
    5. Groth M
    6. Kestler HA
    7. Hildebrandt TB
    8. Platzer M

    (2018) Naked mole-rat transcriptome signatures of socially suppressed sexual maturation and links of reproduction to aging

    BMC Biology 16:77.

    https://doi.org/10.1186/s12915-018-0546-z

    • PubMed
    • Google Scholar
17. 1. Berryman DE
    2. Christiansen JS
    3. Johannsson G
    4. Thorner MO
    5. Kopchick JJ

    (2008) Role of the GH/IGF-1 Axis in lifespan and healthspan: lessons from animal models

    Growth Hormone & IGF Research 18:455–471.

    https://doi.org/10.1016/j.ghir.2008.05.005

    • PubMed
    • Google Scholar
18. 1. Beuschlein F
    2. Fassnacht M
    3. Klink A
    4. Allolio B
    5. Reincke M

    (2001) ACTH-receptor expression, regulation and role in adrenocortial tumor formation

    European Journal of Endocrinology 144:199–206.

    https://doi.org/10.1530/eje.0.1440199

    • PubMed
    • Google Scholar
19. 1. Bieging KT
    2. Mello SS
    3. Attardi LD

    (2014) Unravelling mechanisms of p53-mediated tumour suppression

    Nature Reviews Cancer 14:359–370.

    https://doi.org/10.1038/nrc3711

    • PubMed
    • Google Scholar
20. 1. Bodart G
    2. Goffinet L
    3. Morrhaye G
    4. Farhat K
    5. de Saint-Hubert M
    6. Debacq-Chainiaux F
    7. Swine C
    8. Geenen V
    9. Martens HJ

    (2015) Somatotrope GHRH/GH/IGF-1 Axis at the crossroads between immunosenescence and frailty

    Annals of the New York Academy of Sciences 1351:61–67.

    https://doi.org/10.1111/nyas.12857

    • PubMed
    • Google Scholar
21. 1. Buffenstein R

    (2008) Negligible senescence in the longest living rodent, the naked mole-rat: insights from a successfully aging species

    Journal of Comparative Physiology B 178:439–445.

    https://doi.org/10.1007/s00360-007-0237-5

    • PubMed
    • Google Scholar
22. 1. Burda H

    (1995) Individual recognition and incest avoidance in eusocial common mole-rats rather than reproductive suppression by parents

    Experientia 51:411–413.

    https://doi.org/10.1007/BF01928906

    • PubMed
    • Google Scholar
23. 1. Burkhardt T
    2. Schmidt NO
    3. Vettorazzi E
    4. Aberle J
    5. Mengel M
    6. Flitsch J

    (2013) DHEA(S)--a novel marker in Cushing's disease

    Acta Neurochirurgica 155:479–484.

    https://doi.org/10.1007/s00701-012-1596-6

    • PubMed
    • Google Scholar
24. 1. Burnett C
    2. Valentini S
    3. Cabreiro F
    4. Goss M
    5. Somogyvári M
    6. Piper MD
    7. Hoddinott M
    8. Sutphin GL
    9. Leko V
    10. McElwee JJ
    11. Vazquez-Manrique RP
    12. Orfila AM
    13. Ackerman D
    14. Au C
    15. Vinti G
    16. Riesen M
    17. Howard K
    18. Neri C
    19. Bedalov A
    20. Kaeberlein M
    21. Soti C
    22. Partridge L
    23. Gems D

    (2011) Absence of effects of Sir2 overexpression on lifespan in C. elegans and Drosophila

    Nature 477:482–485.

    https://doi.org/10.1038/nature10296

    • PubMed
    • Google Scholar
25. Software

    1. Bushel P

    (2013) Pvca: Principal Variance Component Analysis (PVCA)

    Pvca: Principal Variance Component Analysis (PVCA).

    https://www.niehs.nih.gov/research/resources/software/biostatistics/pvca/index.cfm
26. 1. Cannata D
    2. Vijayakumar A
    3. Fierz Y
    4. LeRoith D

    (2010) The GH/IGF-1 Axis in growth and development: new insights derived from animal models

    Advances in Pediatrics 57:331–351.

    https://doi.org/10.1016/j.yapd.2010.09.003

    • PubMed
    • Google Scholar
27. 1. Carmona JJ
    2. Michan S

    (2016)

    Biology of healthy aging and longevity

    Revista De Investigacion Clinica; Organo Del Hospital De Enfermedades De La Nutricion 68:7–16.

    • PubMed
    • Google Scholar
28. 1. Carmona-Gutierrez D
    2. Hughes AL
    3. Madeo F
    4. Ruckenstuhl C

    (2016) The crucial impact of lysosomes in aging and longevity

    Ageing Research Reviews 32:2–12.

    https://doi.org/10.1016/j.arr.2016.04.009

    • PubMed
    • Google Scholar
29. 1. Carotti S
    2. Guarino MPL
    3. Valentini F
    4. Porzio S
    5. Vespasiani-Gentilucci U
    6. Perrone G
    7. Zingariello M
    8. Gallo P
    9. Cicala M
    10. Picardi A
    11. Morini S

    (2018) Impairment of GH/IGF-1 Axis in the liver of patients with HCV-Related chronic hepatitis

    Hormone and Metabolic Research 50:145–151.

    https://doi.org/10.1055/s-0043-118911

    • PubMed
    • Google Scholar
30. 1. Celec P
    2. Stárka L

    (2003)

    Dehydroepiandrosterone - is the fountain of youth drying out?

    Physiological Research 52:397–407.

    • PubMed
    • Google Scholar
31. 1. Chabre O

    (2014) [Cushing syndrome: physiopathology, etiology and principles of therapy]

    Presse Medicale 43:376–392.

    https://doi.org/10.1016/j.lpm.2014.02.001

    • PubMed
    • Google Scholar
32. 1. Chahal HS
    2. Drake WM

    (2007) The endocrine system and ageing

    The Journal of Pathology 211:173–180.

    https://doi.org/10.1002/path.2110

    • PubMed
    • Google Scholar
33. 1. Ching T
    2. Huang S
    3. Garmire LX

    (2014) Power analysis and sample size estimation for RNA-Seq differential expression

    RNA 20:1684–1696.

    https://doi.org/10.1261/rna.046011.114

    • PubMed
    • Google Scholar
34. 1. Copeland JM
    2. Cho J
    3. Lo T
    4. Hur JH
    5. Bahadorani S
    6. Arabyan T
    7. Rabie J
    8. Soh J
    9. Walker DW

    (2009) Extension of Drosophila life span by RNAi of the mitochondrial respiratory chain

    Current Biology 19:1591–1598.

    https://doi.org/10.1016/j.cub.2009.08.016

    • PubMed
    • Google Scholar
35. 1. Craig T
    2. Smelick C
    3. Tacutu R
    4. Wuttke D
    5. Wood SH
    6. Stanley H
    7. Janssens G
    8. Savitskaya E
    9. Moskalev A
    10. Arking R
    11. de Magalhães JP

    (2015) The digital ageing atlas: integrating the diversity of age-related changes into a unified resource

    Nucleic Acids Research 43:D873–D878.

    https://doi.org/10.1093/nar/gku843

    • PubMed
    • Google Scholar
36. 1. Dammann P
    2. Šumbera R
    3. Massmann C
    4. Scherag A
    5. Burda H

    (2011) Extended longevity of reproductives appears to be common in Fukomys mole-rats (Rodentia, Bathyergidae)

    PLOS ONE 6:e18757.

    https://doi.org/10.1371/journal.pone.0018757

    • PubMed
    • Google Scholar
37. 1. Dammann P
    2. Sell DR
    3. Begall S
    4. Strauch C
    5. Monnier VM

    (2012) Advanced glycation end-products as markers of aging and longevity in the long-lived Ansell's mole-rat (Fukomys anselli)

    The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 67:573–583.

    https://doi.org/10.1093/gerona/glr208

    • PubMed
    • Google Scholar
38. 1. Dammann P

    (2017) Slow aging in mammals-Lessons from african mole-rats and bats

    Seminars in Cell & Developmental Biology 70:154–163.

    https://doi.org/10.1016/j.semcdb.2017.07.006

    • PubMed
    • Google Scholar
39. 1. Dammann P
    2. Scherag A
    3. Zak N
    4. Szafranski K
    5. Holtze S
    6. Begall S
    7. Burda H
    8. Kestler HA
    9. Hildebrandt T
    10. Platzer M

    (2019) Comment on 'Naked mole-rat mortality rates defy Gompertzian laws by not increasing with age'

    eLife 8:e45415.

    https://doi.org/10.7554/eLife.45415

    • PubMed
    • Google Scholar
40. 1. Dammann P
    2. Burda H

    (2006) Sexual activity and reproduction delay ageing in a mammal

    Current Biology 16:R117–R118.

    https://doi.org/10.1016/j.cub.2006.02.012

    • PubMed
    • Google Scholar
41. 1. Defays D

    (1977) An efficient algorithm for a complete link method

    The Computer Journal 20:364–366.

    https://doi.org/10.1093/comjnl/20.4.364

    • Google Scholar
42. 1. Dillin A
    2. Hsu AL
    3. Arantes-Oliveira N
    4. Lehrer-Graiwer J
    5. Hsin H
    6. Fraser AG
    7. Kamath RS
    8. Ahringer J
    9. Kenyon C

    (2002) Rates of behavior and aging specified by mitochondrial function during development

    Science 298:2398–2401.

    https://doi.org/10.1126/science.1077780

    • PubMed
    • Google Scholar
43. 1. Edrey YH
    2. Salmon AB

    (2014) Revisiting an age-old question regarding oxidative stress

    Free Radical Biology and Medicine 71:368–378.

    https://doi.org/10.1016/j.freeradbiomed.2014.03.038

    • Google Scholar
44. 1. Elibol B
    2. Kilic U

    (2018) High levels of SIRT1 expression as a protective mechanism against Disease-Related conditions

    Frontiers in Endocrinology 9:614.

    https://doi.org/10.3389/fendo.2018.00614

    • Google Scholar
45. 1. Etxabe J
    2. Vazquez JA

    (1994) Morbidity and mortality in Cushing's disease: an epidemiological approach

    Clinical Endocrinology 40:479–484.

    https://doi.org/10.1111/j.1365-2265.1994.tb02486.x

    • Google Scholar
46. 1. Evangelou M
    2. Rendon A
    3. Ouwehand WH
    4. Wernisch L
    5. Dudbridge F

    (2012) Comparison of methods for competitive tests of pathway analysis

    PLOS ONE 7:e41018.

    https://doi.org/10.1371/journal.pone.0041018

    • PubMed
    • Google Scholar
47. 1. Ezkurdia I
    2. Rodriguez JM
    3. Carrillo-de Santa Pau E
    4. Vázquez J
    5. Valencia A
    6. Tress ML

    (2015) Most highly expressed protein-coding genes have a single dominant isoform

    Journal of Proteome Research 14:1880–1887.

    https://doi.org/10.1021/pr501286b

    • PubMed
    • Google Scholar
48. 1. Feng Z
    2. Lin M
    3. Wu R

    (2011) The regulation of aging and longevity: a new and complex role of p53

    Genes & Cancer 2:443–452.

    https://doi.org/10.1177/1947601911410223

    • PubMed
    • Google Scholar
49. 1. Ferraù F
    2. Korbonits M

    (2015) Metabolic comorbidities in Cushing's syndrome

    European Journal of Endocrinology 173:M133–M157.

    https://doi.org/10.1530/EJE-15-0354

    • PubMed
    • Google Scholar
50. 1. Fridley BL
    2. Jenkins GD
    3. Biernacka JM

    (2010) Self-contained gene-set analysis of expression data: an evaluation of existing and novel methods

    PLOS ONE 5:e12693.

    https://doi.org/10.1371/journal.pone.0012693

    • PubMed
    • Google Scholar
51. 1. Fridmanis D
    2. Roga A
    3. Klovins J

    (2017) ACTH receptor (MC2R) Specificity: what do we know about underlying molecular mechanisms?

    Frontiers in Endocrinology 8:13.

    https://doi.org/10.3389/fendo.2017.00013

    • Google Scholar
52. 1. Garcia Montero A
    2. Burda H
    3. Begall S

    (2015) Chemical restraint of african mole-rats (Fukomys sp.) with a combination of ketamine and xylazine

    Veterinary Anaesthesia and Analgesia 42:187–191.

    https://doi.org/10.1111/vaa.12180

    • PubMed
    • Google Scholar
53. 1. Giannakou ME
    2. Goss M
    3. Jünger MA
    4. Hafen E
    5. Leevers SJ
    6. Partridge L

    (2004) Long-lived Drosophila with overexpressed dFOXO in adult fat body

    Science 305:361.

    https://doi.org/10.1126/science.1098219

    • PubMed
    • Google Scholar
54. 1. Gjerstad JK
    2. Lightman SL
    3. Spiga F

    (2018) Role of glucocorticoid negative feedback in the regulation of HPA Axis pulsatility

    Stress 21:403–416.

    https://doi.org/10.1080/10253890.2018.1470238

    • PubMed
    • Google Scholar
55. 1. Hammer F
    2. Subtil S
    3. Lux P
    4. Maser-Gluth C
    5. Stewart PM
    6. Allolio B
    7. Arlt W

    (2005) No evidence for hepatic conversion of dehydroepiandrosterone (DHEA) sulfate to DHEA: in vivo and in vitro studies

    The Journal of Clinical Endocrinology & Metabolism 90:3600–3605.

    https://doi.org/10.1210/jc.2004-2386

    • PubMed
    • Google Scholar
56. 1. Hands SL
    2. Proud CG
    3. Wyttenbach A

    (2009) mTOR's role in ageing: protein synthesis or autophagy?

    Aging 1:586–597.

    https://doi.org/10.18632/aging.100070

    • PubMed
    • Google Scholar
57. 1. Harman D

    (1956) Aging: a theory based on free radical and radiation chemistry

    Journal of Gerontology 11:298–300.

    https://doi.org/10.1093/geronj/11.3.298

    • PubMed
    • Google Scholar
58. 1. Harman D

    (2001) Aging: overview

    Annals of the New York Academy of Sciences 928:1–21.

    https://doi.org/10.1111/j.1749-6632.2001.tb05631.x

    • PubMed
    • Google Scholar
59. 1. Heard NA
    2. Rubin-Delanchy P

    (2018) Choosing between methods of combining $p$-values

    Biometrika 105:239–246.

    https://doi.org/10.1093/biomet/asx076

    • Google Scholar
60. 1. Henning Y
    2. Vole C
    3. Begall S
    4. Bens M
    5. Broecker-Preuss M
    6. Sahm A
    7. Szafranski K
    8. Burda H
    9. Dammann P

    (2014) Unusual ratio between free thyroxine and free triiodothyronine in a long-lived mole-rat species with bimodal ageing

    PLOS ONE 9:e113698.

    https://doi.org/10.1371/journal.pone.0113698

    • PubMed
    • Google Scholar
61. 1. Hipkiss AR

    (2007) On why decreasing protein synthesis can increase lifespan

    Mechanisms of Ageing and Development 128:412–414.

    https://doi.org/10.1016/j.mad.2007.03.002

    • PubMed
    • Google Scholar
62. 1. Hochberg I
    2. Harvey I
    3. Tran QT
    4. Stephenson EJ
    5. Barkan AL
    6. Saltiel AR
    7. Chandler WF
    8. Bridges D

    (2015) Gene expression changes in subcutaneous adipose tissue due to Cushing's disease

    Journal of Molecular Endocrinology 55:81–94.

    https://doi.org/10.1530/JME-15-0119

    • PubMed
    • Google Scholar
63. 1. Hofmann JW
    2. Zhao X
    3. De Cecco M
    4. Peterson AL
    5. Pagliaroli L
    6. Manivannan J
    7. Hubbard GB
    8. Ikeno Y
    9. Zhang Y
    10. Feng B
    11. Li X
    12. Serre T
    13. Qi W
    14. Van Remmen H
    15. Miller RA
    16. Bath KG
    17. de Cabo R
    18. Xu H
    19. Neretti N
    20. Sedivy JM

    (2015) Reduced expression of MYC increases longevity and enhances healthspan

    Cell 160:477–488.

    https://doi.org/10.1016/j.cell.2014.12.016

    • PubMed
    • Google Scholar
64. 1. Holzenberger M
    2. Dupont J
    3. Ducos B
    4. Leneuve P
    5. Géloën A
    6. Even PC
    7. Cervera P
    8. Le Bouc Y

    (2003) IGF-1 receptor regulates lifespan and resistance to oxidative stress in mice

    Nature 421:182–187.

    https://doi.org/10.1038/nature01298

    • PubMed
    • Google Scholar
65. 1. Houthoofd K
    2. Braeckman BP
    3. Johnson TE
    4. Vanfleteren JR

    (2004) Extending life-span in C. elegans

    Science 305:1238–1239.

    https://doi.org/10.1126/science.305.5688.1238c

    • PubMed
    • Google Scholar
66. 1. Imai T
    2. Sarkar D
    3. Shibata A
    4. Funahashi H
    5. Morita-Matsuyama T
    6. Kikumori T
    7. Ohmori S
    8. Seo H

    (2001) Expression of adrenocorticotropin receptor gene in adrenocortical adenomas from patients with cushing syndrome: possible contribution for the autonomous production of cortisol

    Annals of Surgery 234:85–91.

    https://doi.org/10.1097/00000658-200107000-00013

    • PubMed
    • Google Scholar
67. 1. Jarvis JUM
    2. Bennett NC

    (1993) Eusociality has evolved independently in two genera of bathyergid mole-rats — but occurs in no other subterranean mammal

    Behavioral Ecology and Sociobiology 33:253–260.

    https://doi.org/10.1007/BF02027122

    • Google Scholar
68. 1. Johnson SC
    2. Rabinovitch PS
    3. Kaeberlein M

    (2013) mTOR is a key modulator of ageing and age-related disease

    Nature 493:338–345.

    https://doi.org/10.1038/nature11861

    • PubMed
    • Google Scholar
69. 1. Jones CM
    2. Boelaert K

    (2015) The endocrinology of ageing: a Mini-Review

    Gerontology 61:291–300.

    https://doi.org/10.1159/000367692

    • PubMed
    • Google Scholar
70. 1. Jung HJ
    2. Suh Y

    (2014) Regulation of IGF -1 signaling by microRNAs

    Frontiers in Genetics 5:472.

    https://doi.org/10.3389/fgene.2014.00472

    • PubMed
    • Google Scholar
71. 1. Junnila RK
    2. List EO
    3. Berryman DE
    4. Murrey JW
    5. Kopchick JJ

    (2013) The GH/IGF-1 Axis in ageing and longevity

    Nature Reviews Endocrinology 9:366–376.

    https://doi.org/10.1038/nrendo.2013.67

    • PubMed
    • Google Scholar
72. 1. Kanehisa M
    2. Furumichi M
    3. Tanabe M
    4. Sato Y
    5. Morishima K

    (2017) KEGG: new perspectives on genomes, pathways, diseases and drugs

    Nucleic Acids Research 45:D353–D361.

    https://doi.org/10.1093/nar/gkw1092

    • PubMed
    • Google Scholar
73. 1. Keller L
    2. Jemielity S

    (2006) Social insects as a model to study the molecular basis of ageing

    Experimental Gerontology 41:553–556.

    https://doi.org/10.1016/j.exger.2006.04.002

    • PubMed
    • Google Scholar
74. 1. Kirkwood TBL

    (1977) Evolution of ageing

    Nature 270:301–304.

    https://doi.org/10.1038/270301a0

    • Google Scholar
75. 1. Kooistra HS
    2. Galac S

    (2012) Recent advances in the diagnosis of Cushing's syndrome in dogs

    Topics in Companion Animal Medicine 27:21–24.

    https://doi.org/10.1053/j.tcam.2012.06.001

    • PubMed
    • Google Scholar
76. 1. Kurosu H
    2. Yamamoto M
    3. Clark JD
    4. Pastor JV
    5. Nandi A
    6. Gurnani P
    7. McGuinness OP
    8. Chikuda H
    9. Yamaguchi M
    10. Kawaguchi H
    11. Shimomura I
    12. Takayama Y
    13. Herz J
    14. Kahn CR
    15. Rosenblatt KP
    16. Kuro-o M

    (2005) Suppression of aging in mice by the hormone klotho

    Science 309:1829–1833.

    https://doi.org/10.1126/science.1112766

    • PubMed
    • Google Scholar
77. 1. Kwon Y
    2. Kim J
    3. Lee CY
    4. Kim H

    (2015) Expression of SIRT1 and SIRT3 varies according to age in mice

    Anatomy & Cell Biology 48:54–61.

    https://doi.org/10.5115/acb.2015.48.1.54

    • PubMed
    • Google Scholar
78. 1. Langfelder P
    2. Horvath S

    (2008) WGCNA: an R package for weighted correlation network analysis

    BMC Bioinformatics 9:559.

    https://doi.org/10.1186/1471-2105-9-559

    • PubMed
    • Google Scholar
79. 1. Law CW
    2. Chen Y
    3. Shi W
    4. Smyth GK

    (2014) Voom: precision weights unlock linear model analysis tools for RNA-seq read counts

    Genome Biology 15:R29.

    https://doi.org/10.1186/gb-2014-15-2-r29

    • PubMed
    • Google Scholar
80. 1. Li H
    2. Durbin R

    (2009) Fast and accurate short read alignment with Burrows-Wheeler transform

    Bioinformatics 25:1754–1760.

    https://doi.org/10.1093/bioinformatics/btp324

    • PubMed
    • Google Scholar
81. 1. Liberzon A
    2. Birger C
    3. Thorvaldsdóttir H
    4. Ghandi M
    5. Mesirov JP
    6. Tamayo P

    (2015) The molecular signatures database (MSigDB) hallmark gene set collection

    Cell Systems 1:417–425.

    https://doi.org/10.1016/j.cels.2015.12.004

    • PubMed
    • Google Scholar
82. 1. Lin L
    2. Hindmarsh PC
    3. Metherell LA
    4. Alzyoud M
    5. Al-Ali M
    6. Brain CE
    7. Clark AJL
    8. Dattani MT
    9. Achermann JC

    (2007) Severe loss-of-function mutations in the adrenocorticotropin receptor (ACTHR, MC2R) can be found in patients diagnosed with salt-losing adrenal hypoplasia

    Clinical Endocrinology 66:205–210.

    https://doi.org/10.1111/j.1365-2265.2006.02709.x

    • Google Scholar
83. 1. Liu Y
    2. Ordovas JM
    3. Gao G
    4. Province M
    5. Straka RJ
    6. Tsai MY
    7. Lai CQ
    8. Zhang K
    9. Borecki I
    10. Hixson JE
    11. Allison DB
    12. Arnett DK

    (2008) The SCARB1 gene is associated with lipid response to dietary and pharmacological interventions

    Journal of Human Genetics 53:709–717.

    https://doi.org/10.1007/s10038-008-0302-2

    • PubMed
    • Google Scholar
84. 1. Liu D
    2. Ahmet A
    3. Ward L
    4. Krishnamoorthy P
    5. Mandelcorn ED
    6. Leigh R
    7. Brown JP
    8. Cohen A
    9. Kim H

    (2013) A practical guide to the monitoring and management of the complications of systemic corticosteroid therapy

    Allergy, Asthma & Clinical Immunology 9:30.

    https://doi.org/10.1186/1710-1492-9-30

    • PubMed
    • Google Scholar
85. 1. Lodish MB
    2. Hsiao HP
    3. Serbis A
    4. Sinaii N
    5. Rothenbuhler A
    6. Keil MF
    7. Boikos SA
    8. Reynolds JC
    9. Stratakis CA

    (2010) Effects of cushing disease on bone mineral density in a pediatric population

    The Journal of Pediatrics 156:1001–1005.

    https://doi.org/10.1016/j.jpeds.2009.12.027

    • PubMed
    • Google Scholar
86. 1. Longo VD
    2. Antebi A
    3. Bartke A
    4. Barzilai N
    5. Brown-Borg HM
    6. Caruso C
    7. Curiel TJ
    8. de Cabo R
    9. Franceschi C
    10. Gems D
    11. Ingram DK
    12. Johnson TE
    13. Kennedy BK
    14. Kenyon C
    15. Klein S
    16. Kopchick JJ
    17. Lepperdinger G
    18. Madeo F
    19. Mirisola MG
    20. Mitchell JR
    21. Passarino G
    22. Rudolph KL
    23. Sedivy JM
    24. Shadel GS
    25. Sinclair DA
    26. Spindler SR
    27. Suh Y
    28. Vijg J
    29. Vinciguerra M
    30. Fontana L

    (2015) Interventions to slow aging in humans: are we ready?

    Aging Cell 14:497–510.

    https://doi.org/10.1111/acel.12338

    • PubMed
    • Google Scholar
87. 1. López-Otín C
    2. Blasco MA
    3. Partridge L
    4. Serrano M
    5. Kroemer G

    (2013) The hallmarks of aging

    Cell 153:1194–1217.

    https://doi.org/10.1016/j.cell.2013.05.039

    • PubMed
    • Google Scholar
88. 1. Love MI
    2. Huber W
    3. Anders S

    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

    Genome Biology 15:550.

    https://doi.org/10.1186/s13059-014-0550-8

    • PubMed
    • Google Scholar
89. 1. Lowe G
    2. Rumley A

    (2014) The relevance of coagulation in cardiovascular disease: what do the biomarkers tell Us?

    Thrombosis and Haemostasis 112:860–867.

    https://doi.org/10.1160/TH14-03-0199

    • PubMed
    • Google Scholar
90. 1. Lozier NR
    2. Kopchick JJ
    3. de Lacalle S

    (2018) Relative contributions of myostatin and the GH/IGF-1 Axis in body composition and muscle strength

    Frontiers in Physiology 9:1418.

    https://doi.org/10.3389/fphys.2018.01418

    • PubMed
    • Google Scholar
91. 1. Luo W
    2. Brouwer C

    (2013) Pathview: an R/Bioconductor package for pathway-based data integration and visualization

    Bioinformatics 29:1830–1831.

    https://doi.org/10.1093/bioinformatics/btt285

    • PubMed
    • Google Scholar
92. 1. Maggio M

    (2007) Relationship between low levels of anabolic hormones and 6-Year mortality in older men the aging in the chianti area (InCHIANTI) Study

    Archives of Internal Medicine 167:2249–2254.

    https://doi.org/10.1001/archinte.167.20.2249

    • Google Scholar
93. 1. Manimaran S
    2. Selby HM
    3. Okrah K
    4. Ruberman C
    5. Leek JT
    6. Quackenbush J
    7. Haibe-Kains B
    8. Bravo HC
    9. Johnson WE

    (2016) BatchQC: interactive software for evaluating sample and batch effects in genomic data

    Bioinformatics 32:3836–3838.

    https://doi.org/10.1093/bioinformatics/btw538

    • PubMed
    • Google Scholar
94. 1. McCue PM

    (2002) Equine Cushing's disease

    Veterinary Clinics of North America: Equine Practice 18:533–543.

    https://doi.org/10.1016/S0749-0739(02)00038-X

    • PubMed
    • Google Scholar
95. 1. Meijer JC
    2. Lubberink AA
    3. Gruys E

    (1978)

    Cushing's syndrome due to adrenocortical adenoma in a cat

    Tijdschrift Voor Diergeneeskunde 103:1048–1051.

    • PubMed
    • Google Scholar
96. 1. Miller WL
    2. Auchus RJ

    (2011) The molecular biology, biochemistry, and physiology of human steroidogenesis and its disorders

    Endocrine Reviews 32:81–151.

    https://doi.org/10.1210/er.2010-0013

    • PubMed
    • Google Scholar
97. 1. Milman S
    2. Huffman DM
    3. Barzilai N

    (2016) The somatotropic Axis in human aging: framework for the current state of knowledge and future research

    Cell Metabolism 23:980–989.

    https://doi.org/10.1016/j.cmet.2016.05.014

    • PubMed
    • Google Scholar
98. 1. Møller AP

    (2006) Sociality, age at first reproduction and senescence: comparative analyses of birds

    Journal of Evolutionary Biology 19:682–689.

    https://doi.org/10.1111/j.1420-9101.2005.01065.x

    • PubMed
    • Google Scholar
99. 1. Morimoto RI
    2. Cuervo AM

    (2009) Protein homeostasis and aging: taking care of proteins from the cradle to the grave

    The Journals of Gerontology Series A: Biological Sciences and Medical Sciences 64:167–170.

    https://doi.org/10.1093/gerona/gln071

    • PubMed
    • Google Scholar
100. 1. Ochi A
     2. Adachi T
     3. Inokuchi K
     4. Ogawa K
     5. Nakamura Y
     6. Chiba Y
     7. Kawasaki S
     8. Onishi Y
     9. Onuma Y
     10. Munetsugu Y
     11. Ito H
     12. Onuki T
     13. Minoura Y
     14. Watanabe N
     15. Kawamura M
     16. Asano T
     17. Kobayashi Y

     (2016) Effects of aging on the coagulation fibrinolytic system in outpatients of the cardiovascular department

     Circulation Journal 80:2133–2140.

     https://doi.org/10.1253/circj.CJ-16-0530

     • PubMed
     • Google Scholar
101. 1. Parker JD
     2. Parker KM
     3. Sohal BH
     4. Sohal RS
     5. Keller L

     (2004) Decreased expression of Cu-Zn superoxide dismutase 1 in ants with extreme lifespan

     PNAS 101:3486–3489.

     https://doi.org/10.1073/pnas.0400222101

     • PubMed
     • Google Scholar
102. 1. Pérez VI
     2. Buffenstein R
     3. Masamsetti V
     4. Leonard S
     5. Salmon AB
     6. Mele J
     7. Andziak B
     8. Yang T
     9. Edrey Y
     10. Friguet B
     11. Ward W
     12. Richardson A
     13. Chaudhuri A

     (2009) Protein stability and resistance to oxidative stress are determinants of longevity in the longest-living rodent, the naked mole-rat

     PNAS 106:3059–3064.

     https://doi.org/10.1073/pnas.0809620106

     • PubMed
     • Google Scholar
103. 1. Phuc Le P
     2. Friedman JR
     3. Schug J
     4. Brestelli JE
     5. Parker JB
     6. Bochkis IM
     7. Kaestner KH

     (2005) Glucocorticoid receptor-dependent gene regulatory networks

     PLOS Genetics 1:e16.

     https://doi.org/10.1371/journal.pgen.0010016

     • PubMed
     • Google Scholar
104. 1. Pistritto G
     2. Trisciuoglio D
     3. Ceci C
     4. Garufi A
     5. D'Orazi G

     (2016) Apoptosis as anticancer mechanism: function and dysfunction of its modulators and targeted therapeutic strategies

     Aging 8:603–619.

     https://doi.org/10.18632/aging.100934

     • PubMed
     • Google Scholar
105. 1. Pitt JN
     2. Kaeberlein M

     (2015) Why is aging conserved and what can we do about it?

     PLOS Biology 13:e1002131.

     https://doi.org/10.1371/journal.pbio.1002131

     • PubMed
     • Google Scholar
106. 1. Poole W
     2. Gibbs DL
     3. Shmulevich I
     4. Bernard B
     5. Knijnenburg TA

     (2016) Combining dependent P-values with an empirical adaptation of Brown's method

     Bioinformatics 32:i430–i436.

     https://doi.org/10.1093/bioinformatics/btw438

     • PubMed
     • Google Scholar
107. 1. Pride H
     2. Yu Z
     3. Sunchu B
     4. Mochnick J
     5. Coles A
     6. Zhang Y
     7. Buffenstein R
     8. Hornsby PJ
     9. Austad SN
     10. Pérez VI

     (2015) Long-lived species have improved proteostasis compared to phylogenetically-related shorter-lived species

     Biochemical and Biophysical Research Communications 457:669–675.

     https://doi.org/10.1016/j.bbrc.2015.01.046

     • PubMed
     • Google Scholar
108. 1. Raisingani M
     2. Preneet B
     3. Kohn B
     4. Yakar S

     (2017) Skeletal growth and bone mineral acquisition in type 1 diabetic children; abnormalities of the GH/IGF-1 Axis

     Growth Hormone & IGF Research 34:13–21.

     https://doi.org/10.1016/j.ghir.2017.04.003

     • PubMed
     • Google Scholar
109. 1. Ristow M

     (2014) Unraveling the truth about antioxidants: mitohormesis explains ROS-induced health benefits

     Nature Medicine 20:709–711.

     https://doi.org/10.1038/nm.3624

     • PubMed
     • Google Scholar
110. 1. Rodriguez KA
     2. Edrey YH
     3. Osmulski P
     4. Gaczynska M
     5. Buffenstein R

     (2012) Altered composition of liver proteasome assemblies contributes to enhanced proteasome activity in the exceptionally long-lived naked mole-rat

     PLOS ONE 7:e35890.

     https://doi.org/10.1371/journal.pone.0035890

     • PubMed
     • Google Scholar
111. 1. Ruby JG
     2. Smith M
     3. Buffenstein R

     (2018) Naked Mole-Rat mortality rates defy gompertzian laws by not increasing with age

     eLife 7:e31157.

     https://doi.org/10.7554/eLife.31157

     • PubMed
     • Google Scholar
112. 1. Rudman D
     2. Feller AG
     3. Nagraj HS
     4. Gergans GA
     5. Lalitha PY
     6. Goldberg AF
     7. Schlenker RA
     8. Cohn L
     9. Rudman IW
     10. Mattson DE

     (1990) Effects of human growth hormone in men over 60 years old

     New England Journal of Medicine 323:1–6.

     https://doi.org/10.1056/NEJM199007053230101

     • Google Scholar
113. 1. Rutkowski K
     2. Sowa P
     3. Rutkowska-Talipska J
     4. Kuryliszyn-Moskal A
     5. Rutkowski R

     (2014) Dehydroepiandrosterone (DHEA): hypes and hopes

     Drugs 74:1195–1207.

     https://doi.org/10.1007/s40265-014-0259-8

     • PubMed
     • Google Scholar
114. 1. Saez I
     2. Vilchez D

     (2014) The mechanistic links between proteasome activity, aging and Age-related diseases

     Current Genomics 15:38–51.

     https://doi.org/10.2174/138920291501140306113344

     • PubMed
     • Google Scholar
115. 1. Sahm A
     2. Bens M
     3. Henning Y
     4. Vole C
     5. Groth M
     6. Schwab M
     7. Hoffmann S
     8. Platzer M
     9. Szafranski K
     10. Dammann P

     (2018a) Higher gene expression stability during aging in long-lived giant mole-rats than in short-lived rats

     Aging 10:3938–3956.

     https://doi.org/10.18632/aging.101683

     • PubMed
     • Google Scholar
116. 1. Sahm A
     2. Bens M
     3. Szafranski K
     4. Holtze S
     5. Groth M
     6. Görlach M
     7. Calkhoven C
     8. Müller C
     9. Schwab M
     10. Kraus J
     11. Kestler HA
     12. Cellerino A
     13. Burda H
     14. Hildebrandt T
     15. Dammann P
     16. Platzer M

     (2018b) Long-lived rodents reveal signatures of positive selection in genes associated with lifespan

     PLOS Genetics 14:e1007272.

     https://doi.org/10.1371/journal.pgen.1007272

     • PubMed
     • Google Scholar
117. 1. Salmon AB
     2. Leonard S
     3. Masamsetti V
     4. Pierce A
     5. Podlutsky AJ
     6. Podlutskaya N
     7. Richardson A
     8. Austad SN
     9. Chaudhuri AR

     (2009) The long lifespan of two bat species is correlated with resistance to protein oxidation and enhanced protein homeostasis

     The FASEB Journal 23:2317–2326.

     https://doi.org/10.1096/fj.08-122523

     • PubMed
     • Google Scholar
118. 1. Samaras N
     2. Samaras D
     3. Forster A
     4. Frangos E
     5. Philippe J

     (2015)

     Age related dehydroepiandrosterone decrease: clinical significance and therapeutic interest]

     Revue Medicale Suisse 11:321–324.

     • PubMed
     • Google Scholar
119. 1. Satoh A
     2. Brace CS
     3. Rensing N
     4. Cliften P
     5. Wozniak DF
     6. Herzog ED
     7. Yamada KA
     8. Imai S

     (2013) Sirt1 extends life span and delays aging in mice through the regulation of Nk2 homeobox 1 in the DMH and LH

     Cell Metabolism 18:416–430.

     https://doi.org/10.1016/j.cmet.2013.07.013

     • PubMed
     • Google Scholar
120. 1. Scharff A
     2. Locker-Grütjen O
     3. Kawalika M
     4. Burda H

     (2001) NATURAL HISTORY OF THE GIANT MOLE-RAT, CRYPTOMYS MECHOWI (RODENTIA: BATHYERGIDAE), FROM ZAMBIA

     Journal of Mammalogy 82:1003–1015.

     https://doi.org/10.1644/1545-1542(2001)082<1003:NHOTGM>2.0.CO;2

     • Google Scholar
121. 1. Schiaffino S
     2. Mammucari C

     (2011) Regulation of skeletal muscle growth by the IGF1-Akt/PKB pathway: insights from genetic models

     Skeletal Muscle 1:4.

     https://doi.org/10.1186/2044-5040-1-4

     • PubMed
     • Google Scholar
122. 1. Schielke CKM
     2. Burda H
     3. Henning Y
     4. Okrouhlík J
     5. Begall S

     (2017) Higher resting metabolic rate in long-lived breeding Ansell's mole-rats (Fukomys anselli)

     Frontiers in Zoology 14:45.

     https://doi.org/10.1186/s12983-017-0229-6

     • PubMed
     • Google Scholar
123. 1. Schmidt CM
     2. Blount JD
     3. Bennett NC

     (2014) Reproduction is associated with a tissue-dependent reduction of oxidative stress in eusocial female damaraland mole-rats (Fukomys damarensis)

     PLOS ONE 9:e103286.

     https://doi.org/10.1371/journal.pone.0103286

     • PubMed
     • Google Scholar
124. 1. Sherman PW
     2. Jarvis JUM

     (2002) Extraordinary life spans of naked mole-rats (Heterocephalus glaber)

     Journal of Zoology 258:307–311.

     https://doi.org/10.1017/S0952836902001437

     • Google Scholar
125. 1. Sichilima AM
     2. Faulkes CG
     3. Bennett NC

     (2008) Field evidence for aseasonality of reproduction and colony size in the Afrotropical giant mole-rat Fukomys mechowii (Rodentia: Bathyergidae)

     African Zoology 43:144–149.

     https://doi.org/10.1080/15627020.2008.11657231

     • Google Scholar
126. 1. Šklíba J
     2. Mazoch V
     3. Patzenhauerová H
     4. Hrouzková E
     5. Lövy M
     6. Kott O
     7. Šumbera R

     (2012) A maze-lover's dream: Burrow architecture, natural history and habitat characteristics of Ansell's mole-rat (Fukomys anselli)

     Mammalian Biology 77:420–427.

     https://doi.org/10.1016/j.mambio.2012.06.004

     • Google Scholar
127. 1. Smyth GK

     (2004) Linear models and empirical bayes methods for assessing differential expression in microarray experiments

     Statistical Applications in Genetics and Molecular Biology 3:1–25.

     https://doi.org/10.2202/1544-6115.1027

     • Google Scholar
128. 1. Starkov AA

     (2008) The role of mitochondria in reactive oxygen species metabolism and signaling

     Annals of the New York Academy of Sciences 1147:37–52.

     https://doi.org/10.1196/annals.1427.015

     • PubMed
     • Google Scholar
129. 1. Tacutu R
     2. Craig T
     3. Budovsky A
     4. Wuttke D
     5. Lehmann G
     6. Taranukha D
     7. Costa J
     8. Fraifeld VE
     9. de Magalhães JP

     (2013) Human ageing genomic resources: integrated databases and tools for the biology and genetics of ageing

     Nucleic Acids Research 41:D1027–D1033.

     https://doi.org/10.1093/nar/gks1155

     • PubMed
     • Google Scholar
130. 1. Tatar M
     2. Bartke A
     3. Antebi A

     (2003) The endocrine regulation of aging by insulin-like signals

     Science 299:1346–1351.

     https://doi.org/10.1126/science.1081447

     • PubMed
     • Google Scholar
131. 1. Traish AM
     2. Kang HP
     3. Saad F
     4. Guay AT

     (2011) Dehydroepiandrosterone (DHEA)--a precursor steroid or an active hormone in human physiology

     The Journal of Sexual Medicine 8:2960–2982.

     https://doi.org/10.1111/j.1743-6109.2011.02523.x

     • PubMed
     • Google Scholar
132. 1. Valenzano DR
     2. Terzibasi E
     3. Genade T
     4. Cattaneo A
     5. Domenici L
     6. Cellerino A

     (2006) Resveratrol prolongs lifespan and retards the onset of age-related markers in a short-lived vertebrate

     Current Biology 16:296–300.

     https://doi.org/10.1016/j.cub.2005.12.038

     • PubMed
     • Google Scholar
133. 1. van Heemst D

     (2010)

     Insulin, IGF-1 and longevity

     Aging and Disease 1:147–157.

     • PubMed
     • Google Scholar
134. 1. Väremo L
     2. Nielsen J
     3. Nookaew I

     (2013) Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods

     Nucleic Acids Research 41:4378–4391.

     https://doi.org/10.1093/nar/gkt111

     • PubMed
     • Google Scholar
135. 1. Walker JJ
     2. Spiga F
     3. Gupta R
     4. Zhao Z
     5. Lightman SL
     6. Terry JR

     (2015) Rapid intra-adrenal feedback regulation of glucocorticoid synthesis

     Journal of the Royal Society Interface 12:20140875.

     https://doi.org/10.1098/rsif.2014.0875

     • Google Scholar
136. 1. Wang J
     2. Mitsche MA
     3. Lütjohann D
     4. Cohen JC
     5. Xie XS
     6. Hobbs HH

     (2015) Relative roles of ABCG5/ABCG8 in liver and intestine

     Journal of Lipid Research 56:319–330.

     https://doi.org/10.1194/jlr.M054544

     • PubMed
     • Google Scholar
137. 1. Webb SJ
     2. Geoghegan TE
     3. Prough RA
     4. Michael Miller KK

     (2006) The biological actions of dehydroepiandrosterone involves multiple receptors

     Drug Metabolism Reviews 38:89–116.

     https://doi.org/10.1080/03602530600569877

     • PubMed
     • Google Scholar
138. 1. Zeugswetter F
     2. Fenske M
     3. Hassan J
     4. Künzel F

     (2007) Cushing's syndrome in a guinea pig

     Veterinary Record 160:878–880.

     https://doi.org/10.1136/vr.160.25.878

     • Google Scholar

Author details

1. Arne Sahm

   Computational Biology Group, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

   Contribution

   Resources, Data curation, Software, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

   For correspondence

   arne.sahm@leibniz-fli.de

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7330-1790

2. Matthias Platzer

   Computational Biology Group, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

   Contribution

   Conceptualization, Resources, Supervision, Funding acquisition, Visualization, Writing - original draft, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-0596-8582

3. Philipp Koch

   Core Facility Life Science Computing, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

   Contribution

   Visualization, Writing - original draft, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-2825-7943

4. Yoshiyuki Henning

   Institute of Physiology, University Hospital, University of Duisburg-Essen, Essen, Germany

   Contribution

   Resources, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-0166-2204

5. Martin Bens

   Core Facility Sequencing, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

   Contribution

   Resources, Formal analysis, Visualization, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-7940-6668

6. Marco Groth

   Core Facility Sequencing, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

   Contribution

   Resources, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-9199-8990

7. Hynek Burda

   1. Department of General Zoology, Faculty of Biology, University of Duisburg-Essen, Essen, Germany
   2. Department of Game Management and Wildlife Biology, Faculty of Forestry and Wood Sciences, Czech University of Life Sciences, Prague, Czech Republic

   Contribution

   Resources, Writing - review and editing

   Competing interests

   No competing interests declared

8. Sabine Begall

   Department of General Zoology, Faculty of Biology, University of Duisburg-Essen, Essen, Germany

   Contribution

   Conceptualization, Resources, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-9907-6387

9. Saskia Ting

   Institute of Pathology and Neuropathology, University Hospital, University of Duisburg-Essen, Essen, Germany

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

10. Moritz Goetz

    Institute of Pathology and Neuropathology, University Hospital, University of Duisburg-Essen, Essen, Germany

    Contribution

    Investigation, Writing - review and editing

    Competing interests

    No competing interests declared

11. Paul Van Daele

    Department of Zoology, University of South Bohemia, České Budějovice, Czech Republic

    Contribution

    Resources, Writing - review and editing

    Competing interests

    No competing interests declared

12. Magdalena Staniszewska

    Department of Nuclear Medicine, University Hospital, University of Duisburg-Essen, Essen, Germany

    Contribution

    Investigation, Writing - review and editing

    Competing interests

    No competing interests declared

13. Jasmin Mona Klose

    Department of Nuclear Medicine, University Hospital, University of Duisburg-Essen, Essen, Germany

    Contribution

    Investigation, Writing - review and editing

    Competing interests

    No competing interests declared

14. Pedro Fragoso Costa

    Department of Nuclear Medicine, University Hospital, University of Duisburg-Essen, Essen, Germany

    Contribution

    Investigation, Writing - review and editing

    Competing interests

    No competing interests declared

15. Steve Hoffmann

    Computational Biology Group, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

    Contribution

    Formal analysis, Supervision, Investigation, Methodology, Writing - original draft, Writing - review and editing

    Contributed equally with

    Karol Szafranski and Philip Dammann

    Competing interests

    No competing interests declared

16. Karol Szafranski

    Core Facility Life Science Computing, Leibniz Institute on Aging – Fritz Lipmann Institute, Jena, Germany

    Contribution

    Conceptualization, Resources, Formal analysis, Funding acquisition, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

    Contributed equally with

    Steve Hoffmann and Philip Dammann

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0001-6391-1766

17. Philip Dammann

    1. Department of General Zoology, Faculty of Biology, University of Duisburg-Essen, Essen, Germany
    2. Central Animal Laboratory, University Hospital, University of Duisburg-Essen, Essen, Germany

    Contribution

    Conceptualization, Resources, Funding acquisition, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing

    Contributed equally with

    Steve Hoffmann and Karol Szafranski

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0002-0624-3965

• 2,781

  views

• 283

  downloads

• 22

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.