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

Acceptance summary:

This study analyzed 636 RNA-seq samples across 15 tissues, which is a massive effort and a rich resource. Changes in the hypothalamic-pituitary-adrenal stress axis is found to be associated with the extended life expectancy of reproductive vs. non-reproductive mole-rats, which is consistent with existing knowledge, whereas other associations are contrary to existing knowledge, such as the up-regulation of the IGF1/GH axis and several other anabolic processes in the long lived species. The authors should tune down their claims of causality such as "play a key role", which was not demonstrated without genetic perturbations, by only concluding associations, which are essentially what were described.

Decision letter after peer review:

Thank you for submitting your article "The HPA stress axis shapes aging rates in long-lived, social mole-rats" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Jessica Tyler as the Senior Editor. The reviewers have opted to remain anonymous.

While we cannot invite revision of the work in its current form, we acknowledge the potential and uniqueness of the data resource you have generated. If in the future you could improve the computational analyses to arrive at more definitive conclusions and include at least limited experimental validation of hypotheses you are advancing, we would be interested to receive such a thoroughly revised submission. We would treat it as a new submission, but we would try to recruit the same reviewers.

Summary:

Differential gene expression between breeding and non-breeding individuals of the two long-lived Fukomys species investigated here is a very interesting topic, and the scale of collected data is impressive. Unfortunately, however, both the data analyses and experimental support for the conclusions are rather superficial and descriptive. Given the size and rarity of the dataset, this could be an important resource for the aging research community, if data quality can be clearly demonstrated by robust and extensive computational analyses and experimental validation.

Without demonstrating causality, the authors should refrain from making assertive mechanistic claims, such as "the HPA stress axis shapes aging rates".

Please see the detailed reviews for specific concerns.

Reviewer #1:

This is a primarily descriptive study reporting gene expression differences between breeding and non-breeding individuals of two long-lived Fukomys species. These mole rat species are interesting for this type of analysis because it has been previously shown that breeders live ~2-fold longer than non-breeders. Thus, it may be possible to learn about mechanisms that determine longevity by studying differences between breeders and non-breeders. Although the manuscript is primarily descriptive, there is value in observational data sets such as this, which can be hypothesis generating and can spur future, more mechanistic studies. I found the presentation to be a bit of a "data dump" with many different comparisons of differentially expressed genes, GO terms, etc., but without a lot in terms true biological insight. As it was largely limited to gene expression, it is also a bit one-dimensional. This work will likely be of interest and value to scientists who study the comparative biology of aging and systems biology of aging, but in its current form seems somewhat incremental and probably will not capture the interest of a broad scientific audience or even the broader gerontology community.

My largest concerns with the content of this manuscript are related to interpretation of the data. While this type of comparative gene expression analysis can generate hypotheses, it cannot strongly support or refute causal relationships without additional experimentation. The authors repeatedly appear to interpret correlative observations with causation and make claims that overreach the data. The title itself is a good example of this where the authors claim "the HPA stress axis shapes aging rates". There is no direct evidence to support this claim or others similar in nature throughout the manuscript.

Another area of overinterpretation is with respect to the relevance of these findings as a test of the validity of prior aging theories or mechanisms. It seems very likely that this is a somewhat unique evolutionary case where upon the switch from non-breeder to breeder there is a dramatic rewiring of physiology at many levels (transcriptional, translational, post-translational, metabolome, epigenome, etc.). Simply because this switch does not match patterns seen in longevity interventions such as caloric restriction, reduced GH/IGF-1 signaling, or mTOR inhibition does not refute or call into question literature describing potential mechanisms for how those interventions act. It likely simply reflects that this is a different path to achieve longevity.

Specific comments:

– The title is problematic as described above.

– The first sentence in the Abstract is problematic as it implies causality between sexual activity/reproduction and longevity. Sexual activity/reproduction are associated with life expectancy. It is interesting to consider whether these two things could perhaps be uncoupled in this animal, as has been shown for reduced fecundity and longevity in invertebrate models.

– The phrase "oppose crucial findings" in the Abstract is also problematic. First, the word "crucial" does not seem to makes sense in this context. What makes them crucial? Second, it is intuitive and expected, indeed perhaps required, that reproduction would be associated with anabolic processes, and this study does not show causality between the observed changes in these processes and longevity, so it does not actually "oppose" the prior findings in genetic models with reduced IGF-1/GH signaling where lifespan is extended.

– I have a problem with the premise that it is possible to "confirm or falsify" results from cross-species or intra-species studies through the type of approach taken here as implied in the text. Confirming or falsifying results implies something about the quality of the prior data itself. I assume what the authors mean is confirming or falsifying the underlying hypotheses or assumptions. Even that is questionable, however, since the mechanisms by which longevity are determined could simply be different within species versus across species versus cases like this where you have dramatically different life expectancies in breeders versus non-breeders.

– I thought the discussion of the GH/IGF-1 results was fairly balanced, but I would encourage the authors to consider more deeply the within species versus across species observations in the literature. The evidence for reduced GH/IGF-1 increasing lifespan comes from within species studies and appears to hold true from worms to dogs and likely in humans as well, although the correlation between body size and lifespan is a bit more complicated in humans for obvious reasons. Within species, smaller individuals who have reduced GH/IGF-1 tend to live longer – that's a correlation. In worms, flies, and mice a reduction in growth signaling has been shown to be sufficient – and very likely causal – for enhanced longevity through genetic and pharmacological studies. The comment that these interventions are all performed during development is not exactly true – in mice rapamycin at least works in adulthood and even when only given transiently or intermittently in adulthood. Across species, larger species tend to live longer. Perhaps the mechanisms going from non-breeders to breeders more resemble the evolutionary longevity strategies that have been taken at the species level rather than the mechanism that appear to determine longevity within species. Personally, I don't see this as a contradiction or a controversy within the field.

– I find the "short-lived" versus "long-lived" species argument to be overly speculative and arbitrary. Interventions such as CR, mTOR inhibition, reduced IGF-1, etc. extend lifespan at least from worms to mice, which is a >50-fold difference in lifespan. Mice to people is ~30-fold. Compared to mice, worms are shorter-lived (by a fold difference metric) than mice are compared to humans.

– The authors state that breeders and non-breeders have "massively diverging aging rates", but I think it is important to keep in mind that differences in lifespan do not necessarily imply differences in aging rate. Especially going from the long-lived state to the short-lived state. Perhaps there is good evidence that functional and molecular declines and diseases/pathologies of aging are accelerated in the non-breeders and delayed in the breeders, which would support this assertion, but this is unclear from the manuscript.

– The phrase "unilaterally described as harmful" in the conclusion is simply not true. GH/IGF-1 signaling limits lifespan in worms, flies, and mice but none of the papers cited unilaterally claim that it is harmful. In fact, some of those same papers note that high GH/IGF-1 signaling often confers a selective advantage in terms of faster maturation and reproduction. So, at the species level, this is beneficial. Even at the individual level, high GH/IGF-1 may be associated with better outcomes in the wild where predation is a factor. High GF/IGF-1 in these species is detrimental (only?) in the context of aging/longevity in a relatively safe laboratory environment at the individual level.

Reviewer #2:

Sahm et al. have analyzed gene expression of 22 samples across 15 tissues of Fukomys mole rat, and with these resources, they try to explain why breeders and non-breeders have different lifespan in those species. While I was intrigued by the idea and design of the research, I found the bioinformatic approaches, as well as the overall result of the study, fell far short of being acceptable for publication. Below I highlight key aspects of the approach (and related conclusions, or lack thereof) that I believe represent serious issues. One positive note, again, the questions and data of this paper that I believe are highly interesting and important if the author can pursue it in a more focused and solid way.

1) First, in my opinion, the discussion of the paper is not well synthesized, and the content does not help the reader for the aim of study and of what is discussed in the Title (e.g., HPA) and Abstract. While the entire discussion fail to build logics between HPA stress and aging in mole rats, it says story of intervention and hypothesis testing instead, with the last of results immediately jumps into comparisons positive selection genes and DEG (It is unclear why it is relevant even one gene is found under positive selection).

2) Given hundreds of samples have been sequenced in the paper, there is no extensive examination of batch effect, which could ultimately, in its present form, put all the results (e.g., DEG analysis, pathways enrichment, multifactors analysis) and conclusions at high risk. Particularly, the PCA analysis has indicated the species, tissue, and the combination of both variables accounted for 98.4 % of the total variance in the data set, it becomes more important to know how the author have organized the sequencing strategy.

3).Unsophisticated use of enrichment analysis on pathway likely lead to misleading conclusions – a large portion of the analyses presented use an unreasonable approach to identify genes with expression shift across tissues and species to make conclusions that I believe are largely misleading about genes important for longevity. The authors identify pathways that have an overall (gene-wide) changes in gene expression as p value looks like an potential indication of expression shift. While this approach MAY be lucky enough to catch a few of these “true positives", the VAST majority of what it will identify will be pathways with global accumulated changes (which is what the point would be), but rather small pathways or pathways fully of with wired p values or undetectable expression (and thus perhaps some of the least important genes). This overly coarse approach seems to me unreasonable and publishable for this purpose. Then, the extension of these approaches to pathways further muddies the waters of any discussions to the extent that I believe it is largely nonsense that overwhelms the results. The approaches for detecting affected global pathways have been the subject of a very large body of literature, and there are well developed hierarchical/empirical models for testing these hypotheses that the authors have completely and inexcusably ignored. And, of course, it is strong encouraged the author could formulate an new model/index to re-evaluate pathways analysis as the data and experimental design is unique to other study.

4) Last but not the least, small KEGG set, i.e., KEGG that with very few gene could be also confusing such whole categories analysis. The author should check this potential bias using random sampling “pseudo-KEGG set” with same gene number and/or identify a threshold to filer the each true KEGG categories that considered.

Reviewer #3:

The manuscript by Sahm et al. describes the transcriptomic comparison of breeders and non-breeders in two species of mole rats from genus Fukomys. The remarkable aspect of Fukomys mole rats is that the breeders live significantly longer than workers. The authors produced new breeder couples by pairing animals from different family groups and then compared their transcriptomes to non-breeders of the same age from the original colonies.

There were very few differences identified between breeders and non-breeders. Most transcriptomic changes were confined to gonads and endocrine glands. This is somewhat unsurprising and these organs become active for breeding. The pathways that became activated were related to ribosome biogenesis, protein translation, and MYC signaling, which all reflect physiological activation of these tissues in breeders.

Interestingly, some activation of these signaling pathways was also observed in non-gonadal tissues, which contradicts the common dogma that downregulation of these pathways is associated with lifespan extension. This is a remarkable result that suggests that short-lived model organisms may not correctly reflect signaling effects required for longevity in long-lived species. The authors also point out to higher glucocorticoid levels in workers and speculate that may be showing Cushing's like syndrome.

Overall, this is an important study of high interest.

Items to address:

1) It is not clear how long were the animals maintained in a breeder status prior to analysis.

2) Figure 3A, the rational for the comparison of changes in young breeder versus non breeder animals to changes that occur with age over time is unclear.

3) Figure 4, were the changes driven by gonads mainly? Hoe many were non-gonadal? For example, in Figure 2, muscle showed 0 DEGs by status, but in pathway analysis we see upregulation of pathways. Please explain.

4) Table 2: TOR is downregulated in breeders while ribosome processed are upregulated, please explain.

5) The finding of Cushing's syndrome needs more support. Please consider comparing Cushing's related transcriptome changes as a whole to the changes observed in mole rats, otherwise this conclusion may need to be toned down.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for submitting your article "Slower aging due to sexual activity in mole-rats is associated with transcriptional changes in the HPA stress axis" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Jessica Tyler as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Vera Gorbunova (Reviewer #3).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). Specifically, we are asking editors to accept without delay manuscripts, like yours, that they judge can stand as eLife papers without additional data, even if they feel that they would make the manuscript stronger. Thus the revisions requested below only address clarity, data analysis and presentation.

We are glad to see that you smoothed your HPA stress hypothesis according to our suggestions, but further revisions are still necessary to address our previous concerns about data quality and robustness of results have and remove the overinterpretation from the title.

Specifically you need to:

1) Control the batch effect using established packages, such as "BatchQC" or "BEclear";

2) Detect DEG using limma with voom to account for large variances and compare the results with DESeq2's output (see PMID: 23497356);

3) Redo the pathway enrichment with changes of expression values or other static indicators (see PMID: 23625889 and PMID: 32515539 for examples) but not p-value; to remove the bias and validate the HPA hypothesis.

4) Reconsider the title, which is still an over-interpretation of the data and therefore needs to be fixed as follows. The words "slower aging" imply that the same biological aging process is occurring in the breeders and non-breeders, which is not actually clear from the data. Indeed, it seems quite plausible that the non-breeders die prematurely due to too much stress (as you speculate). Whether the rate of aging is different or not remains unclear. Replacing "slower aging" with "Increased longevity" or something like that would address this issue.

Reviewer #1:

This is a primarily descriptive study reporting gene expression differences between breeding and non-breeding individuals of two long-lived Fukomys species. These mole rat species are interesting for this type of analysis because it has been previously shown that breeders live ~2-fold longer than non-breeders. Thus, it may be possible to learn about mechanisms that determine longevity by studying differences between breeders and non-breeders. Although the manuscript is primarily descriptive, there is value in observational data sets such as this, which can be hypothesis generating and can spur future, more mechanistic studies. I found the presentation to be a bit of a "data dump" with many different comparisons of differentially expressed genes, GO terms, etc., but without a lot in terms true biological insight. As it was largely limited to gene expression, it is also a bit one-dimensional. This work will likely be of interest and value to scientists who study the comparative biology of aging and systems biology of aging, but in its current form seems somewhat incremental and probably will not capture the interest of a broad scientific audience or even the broader gerontology community.

My largest concerns with the content of this manuscript are related to interpretation of the data. While this type of comparative gene expression analysis can generate hypotheses, it cannot strongly support or refute causal relationships without additional experimentation. The authors repeatedly appear to interpret correlative observations with causation and make claims that overreach the data. The title itself is a good example of this where the authors claim "the HPA stress axis shapes aging rates". There is no direct evidence to support this claim or others similar in nature throughout the manuscript.

Another area of overinterpretation is with respect to the relevance of these findings as a test of the validity of prior aging theories or mechanisms. It seems very likely that this is a somewhat unique evolutionary case where upon the switch from non-breeder to breeder there is a dramatic rewiring of physiology at many levels (transcriptional, translational, post-translational, metabolome, epigenome, etc.). Simply because this switch does not match patterns seen in longevity interventions such as caloric restriction, reduced GH/IGF-1 signaling, or mTOR inhibition does not refute or call into question literature describing potential mechanisms for how those interventions act. It likely simply reflects that this is a different path to achieve longevity.

We thank you for your outright comment that we have made too strong statements in light of the evidence presented. Your review was the first we received for this manuscript, and considerably helped us to realize the pitfalls of over-interpretation. We implemented your constructive criticism in two ways. First, we have, in accordance with your suggestions, toned down or deleted statements and generally made the correlative character of the facts presented more transparent. In addition, we are now discussing the previously collected and new, complementary evidence for the HPA stress hypothesis in more detail and with more differentiation than before.

Specific comments:

– The title is problematic as described above.

We agree and use a new title that reflects the correlative character of the underlying evidence: Slower aging due to sexual activity in mole-rats is associated with transcriptional changes in the HPA stress axis

– The first sentence in the Abstract is problematic as it implies causality between sexual activity/reproduction and longevity. Sexual activity/reproduction are associated with life expectancy. It is interesting to consider whether these two things could perhaps be uncoupled in this animal, as has been shown for reduced fecundity and longevity in invertebrate models.

Yes, we modified the Abstract according to your suggestions. Specifically, we made the following replacements:

“Sexual activity and/or reproduction are associated with a doubling of life expectancy in the long-lived rodent genus Fukomys.”

“This analysis suggests that changes in the regulation of the hypothalamic-pituitary-adrenal stress axis play a key role regarding the extended life expectancy of reproductive vs. non-reproductive mole-rats. This is substantiated by a corpus of independent evidence.”

– The phrase "oppose crucial findings" in the Abstract is also problematic. First, the word "crucial" does not seem to makes sense in this context. What makes them crucial? Second, it is intuitive and expected, indeed perhaps required, that reproduction would be associated with anabolic processes, and this study does not show causality between the observed changes in these processes and longevity, so it does not actually "oppose" the prior findings in genetic models with reduced IGF-1/GH signaling where lifespan is extended.

We agree and modified the Abstract according to your suggestions:

“On the other hand, several of our results are not consistent with knowledge about aging of short-lived model organisms.”

Moreover, we added to the Introduction:

“Nevertheless, both Fukomys and naked mole-rats can be seen as counter-examples 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).”

– I have a problem with the premise that it is possible to "confirm or falsify" results from cross-species or intra-species studies through the type of approach taken here as implied in the text. Confirming or falsifying results implies something about the quality of the prior data itself. I assume what the authors mean is confirming or falsifying the underlying hypotheses or assumptions. Even that is questionable, however, since the mechanisms by which longevity are determined could simply be different within species versus across species versus cases like this where you have dramatically different life expectancies in breeders versus non-breeders.

Yes, we agree and removed the respective lines from the manuscript.

– I thought the discussion of the GH/IGF-1 results was fairly balanced, but I would encourage the authors to consider more deeply the within species versus across species observations in the literature. The evidence for reduced GH/IGF-1 increasing lifespan comes from within species studies and appears to hold true from worms to dogs and likely in humans as well, although the correlation between body size and lifespan is a bit more complicated in humans for obvious reasons. Within species, smaller individuals who have reduced GH/IGF-1 tend to live longer – that's a correlation. In worms, flies, and mice a reduction in growth signaling has been shown to be sufficient – and very likely causal – for enhanced longevity through genetic and pharmacological studies.

That is, if we see it correctly, so far completely in line with the manuscript, which describes that today’s literature based on evidence in different species predominantly assumes that increased IGF1/GH signaling shortens life and decreased signaling prolongs it (Table2). If we have misunderstood you on this point and in your opinion important arguments or literature evidence is still missing, we would be happy to add it.

The comment that these interventions are all performed during development is not exactly true – in mice rapamycin at least works in adulthood and even when only given transiently or intermittently in adulthood.

Yes, we have supplemented and qualified the statement as follows:

“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…”

Across species, larger species tend to live longer. Perhaps the mechanisms going from non-breeders to breeders more resemble the evolutionary longevity strategies that have been taken at the species level rather than the mechanism that appear to determine longevity within species. Personally, I don't see this as a contradiction or a controversy within the field.

We understand your concerns, however, here we are making the case that there is no need to postulate two general strategies/mechanisms for longevity determination at the inter- and intraspecies levels. For example, in both intra- and interspecies comparisons, the growth rate is clearly anticorrelated with lifespan (https://pubmed.ncbi.nlm.nih.gov/11868658/, https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3574304/, https://pubmed.ncbi.nlm.nih.gov/12954479/). In respect to the well accepted positive correlation of body weight and longevity this means, that large long-lived species grow relative longer than small short-lived ones. Noteworthy, in mammals maximum life span is even stronger correlated to time to maturity than either to adult weight or growth rate (see, e.g., Figure 2 of https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4406664/). In respect to the manuscript, however, we believe that such a speculation and discussion would extend too far from the topics of our work.

– I find the "short-lived" versus "long-lived" species argument to be overly speculative and arbitrary. Interventions such as CR, mTOR inhibition, reduced IGF-1, etc. extend lifespan at least from worms to mice, which is a >50-fold difference in lifespan. Mice to people is ~30-fold. Compared to mice, worms are shorter-lived (by a fold difference metric) than mice are compared to humans.

We would like to point out that the concept of comparing short- versus long-lived has been applied successfully (see e.g. Austad SN, Comparative biology of aging, PMID: 19223603 and Gorbunova V, Comparative genetics of longevity and cancer, PMID: 24981598). Particular in species with a similar biology like mammals, where a strong correlation exists between body-weight and longevity, it is quite rational and powerful. Based on this correlation, short- and long-lived outliers can be objectively identified. Classification of Fukomys as a long-lived outlier is a basic motivation of our research efforts.

However, we understand your concern with respect to the fact that the mentioned treatments work in principle consistently in the same direction in a wide range of species. Therefore, we have shortened the paragraph considerably and concentrated on what we consider a more stringent argument:

“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). […]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.”

– The authors state that breeders and non-breeders have "massively diverging aging rates", but I think it is important to keep in mind that differences in lifespan do not necessarily imply differences in aging rate. Especially going from the long-lived state to the short-lived state. Perhaps there is good evidence that functional and molecular declines and diseases/pathologies of aging are accelerated in the non-breeders and delayed in the breeders, which would support this assertion, but this is unclear from the manuscript.

Sorry for being insufficiently clear in this respect in the submitted manuscript. We agree of course that differences in lifespan do not necessarily imply differences in aging rates, as organisms can die for various reasons, many of which have little to do with physiological aging, and as a matter of fact, we have not actually measured aging rates in our study species. However, the survival data for Fukomys non-breeders and breeders (https://pubmed.ncbi.nlm.nih.gov/16488857/, https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3076438/) have been obtained in protected captive environments where both groups lived essentially under the same conditions, and where starvation, predation or dispersal could be excluded as mortality causes. Accidental death cases were very rare. Kaplan-Meier-survival curves and analyses of mortality rates indicate that both groups show actuarial senescence, but with obviously earlier onset of aging and significantly higher annual mortality hazards in non-breeders across all age classes (https://pubmed.ncbi.nlm.nih.gov/16488857/, https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3076438/, https://cloud.leibniz-fli.de/index.php/s/d2Wob9jFAeyRHQZ [accepted at Philos. Trans. R. Soc. B]). Hence, we are confident to say that the longer mean and absolute lifespans of breeders result from a combination of delayed onset of aging and henceforth lower age-specific mortality rates in Fukomys breeders as compared to non-breeders. As, in the absence of extrinsic mortality causes, both parameters are parameters of aging, we consider it justified to speak of “slower aging breeders” and/or “faster aging non-breeders” whenever appropriate.

However, we modified the text to avoid misunderstandings:

“… compared naturally occurring cohorts of species with massively diverging longevities.”

plus numerous other passages where the phrase “aging-rate” has been substituted likewise.

– The phrase "unilaterally described as harmful" in the conclusion is simply not true. GH/IGF-1 signaling limits lifespan in worms, flies, and mice but none of the papers cited unilaterally claim that it is harmful. In fact, some of those same papers note that high GH/IGF-1 signaling often confers a selective advantage in terms of faster maturation and reproduction. So, at the species level, this is beneficial. Even at the individual level, high GH/IGF-1 may be associated with better outcomes in the wild where predation is a factor. High GF/IGF-1 in these species is detrimental (only?) in the context of aging/longevity in a relatively safe laboratory environment at the individual level.

Yes, we phrased that wrongly. In our discussion, we even deal with the possible evolutionary advantages of increased IGF1/GH signalling in non-breeders in such a tradeoff view.

“In particular, this applies to the role of the GH/IGF1 axis, which has in recent years been predominantly described as having a negative impact on lifespan (Lopez-Otin et al., 2013; Pitt and Kaeberlein, 2015; Bartke, 2017).”

Reviewer #2:

Sahm et al. have analyzed gene expression of 22 samples across 15 tissues of Fukomys mole rat, and with these resources, they try to explain why breeders and non-breeders have different lifespan in those species. While I was intrigued by the idea and design of the research, I found the bioinformatic approaches, as well as the overall result of the study, fell far short of being acceptable for publication. Below I highlight key aspects of the approach (and related conclusions, or lack thereof) that I believe represent serious issues. One positive note, again, the questions and data of this paper that I believe are highly interesting and important if the author can pursue it in a more focused and solid way.

1) First, in my opinion, the discussion of the paper is not well synthesized, and the content does not help the reader for the aim of study and of what is discussed in the Title (e.g., HPA) and Abstract. While the entire discussion fail to build logics between HPA stress and aging in mole rats, it says story of intervention and hypothesis testing instead, with the last of results immediately jumps into comparisons positive selection genes and DEG (It is unclear why it is relevant even one gene is found under positive selection).

In response to the comment, we now discuss the hypotheses, results and verifications concerning the HPA axis in greater detail. Furthermore, we have shortened those parts of the manuscript describing the up-regulation of various anabolic components such as the IGF1/GH axis, which are predominantly associated with a shortening of lifespan (see answer 1.10 to Reviewer #1). You also assert that we jump “immediately” into a comparison of positively selected genes and DEGs. This is very surprising to us, since this issue is addressed only once in the entire manuscript: in the last of about fifteen paragraphs of our Discussion – by far not the main pillar of our manuscript.

2) Given hundreds of samples have been sequenced in the paper, there is no extensive examination of batch effect, which could ultimately, in its present form, put all the results (e.g., DEG analysis, pathways enrichment, multifactors analysis) and conclusions at high risk. Particularly, the PCA analysis has indicated the species, tissue, and the combination of both variables accounted for 98.4 % of the total variance in the data set, it becomes more important to know how the author have organized the sequencing strategy.

We have added the rational one for our sequencing strategy as well as a new supplement table that details the organization of sequencing experiments together. In brief, we have taken great care to distribute samples of breeders and non-breeders across the sequencing runs to avoid such systematic errors.

“ We have taken care to concentrate the samples of a tissue on a few sequencing runs and to create a balance of breeders and non-breeders within each run (Supplementary file 1J).”

3) Unsophisticated use of enrichment analysis on pathway likely lead to misleading conclusions – a large portion of the analyses presented use an unreasonable approach to identify genes with expression shift across tissues and species to make conclusions that I believe are largely misleading about genes important for longevity. The authors identify pathways that have an overall (gene-wide) changes in gene expression as p value looks like an potential indication of expression shift. While this approach MAY be lucky enough to catch a few of these “true positives", the VAST majority of what it will identify will be pathways with global accumulated changes (which is what the point would be), but rather small pathways or pathways fully of with wired p values or undetectable expression (and thus perhaps some of the least important genes). This overly coarse approach seems to me unreasonable and publishable for this purpose. Then, the extension of these approaches to pathways further muddies the waters of any discussions to the extent that I believe it is largely nonsense that overwhelms the results. The approaches for detecting affected global pathways have been the subject of a very large body of literature, and there are well developed hierarchical/empirical models for testing these hypotheses that the authors have completely and inexcusably ignored. And, of course, it is strong encouraged the author could formulate an new model/index to re-evaluate pathways analysis as the data and experimental design is unique to other study.

$ℱ=-2*\sum _{i=1}^m{log}_e(p_i)$

Since we have reason to believe that this could be a methodological misunderstanding (see below), we would like to take the opportunity to explain the basic aspects of our approach again with references to the relevant literature. In order to combine different p-values in a threshold-free way into a test statistic, we have used the Fisher’s method.In fact, Fisher's method is one of the most widely used methods in statistics (https://academic.oup.com/biomet/article/105/1/239/4788722). The provided reference also illustrates the use of weights. This is an appropriate approach to account for different levels of gene expression in the individual tissues in the context of cross-tissue analysis. In addition, the method is also specifically recommended for the identification of differentially expressed pathways in gene expression analyses (e.g. https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0012693,

https://academic.oup.com/biomet/article/105/1/239/4788722, https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0041018, https://www.tandfonline.com/doi/abs/10.1080/19466315.2014.888013, https://bmcbioinformatics.biomedcentral.com/articles/10.1186/1471-2105-14-368, https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-8-96, https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3632109/).

“Simulation results demonstrated that overall Fisher's method and the global model with random effects have the highest power for a wide range of scenarios” (first linked article)

Unfortunately, you did not specify any instance of the "large body of literature" that you consider better suited for threshold-free, weighted enrichment analyses of gene expression comparisons across several tissues.

Strategies used to determine the combined p-value based on Fisher's method, are described, for example, in the first, third and seventh linked article. In general Fisher’s method assumes independence of m individual test statistics and uniformity of the associated p-values under the null. Consequently, the test statistic obtained with Fisher’s method is approximately chi square distributed with 2m degrees of freedom. In gene expression analyses, however, the assumption of independence may be violated (https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0012693, https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3632109/). Therefore, a Monte-Carlo approach is often used instead to determine the p-values with empirically determined null distributions.

We have taken your comments as an opportunity to describe the procedure in more detail:

“If the underlying test statistics are independent, a combined p-value may be determined by Fishers’s method using an χ²-distribution. […] The combined p-values for our hypothesis tests can now be obtained by determining the probability of a test statistic in the empirically estimated null-distribution.”

Regarding the assertions about validity and integrity of our results we would like to refer to our cross-tissue-level analysis. We are able to falsify the hypothesis that our method preferentially identifies poorly expressed genes/pathways or small pathways:

“We checked whether this method may have a detection bias with regard to expression level of genes/pathways or size of pathways at the cross-tissue level. For this, we looked at the correlations of the respective p-values with gene expression, pathway expression, and pathway size (Figure 4—figure supplement 8-10). In summary, we find a slightly enhanced sensitivity for highly expressed pathways at the cross-tissue level.”

Furthermore, a preference for pathways in terms of size can already be excluded on the basis of the empirical determination of the null distribution described above. There, the number of elements in the respective pathway is taken into account (see above).

4) Last but not the least, small KEGG set, i.e., KEGG that with very few gene could be also confusing such whole categories analysis. The author should check this potential bias using random sampling “pseudo-KEGG set” with same gene number and/or identify a threshold to filer the each true KEGG categories that considered.

Confusingly, the description of this “pseudo-KEGG set” approach, which you suggest as a possible solution, matches exactly the method used in our manuscript. Given this last statement, we believe that the fundamental criticism our method may have been the result of a misunderstanding. We sincerely hope to have resolved this with our explanations.

Reviewer #3:

The manuscript by Sahm et al. describes the transcriptomic comparison of breeders and non-breeders in two species of mole rats from genus Fukomys. The remarkable aspect of Fukomys mole rats is that the breeders live significantly longer than workers. The authors produced new breeder couples by pairing animals from different family groups and then compared their transcriptomes to non-breeders of the same age from the original colonies.

There were very few differences identified between breeders and non-breeders. Most transcriptomic changes were confined to gonads and endocrine glands. This is somewhat unsurprising and these organs become active for breeding. The pathways that became activated were related to ribosome biogenesis, protein translation, and MYC signaling, which all reflect physiological activation of these tissues in breeders.

Interestingly, some activation of these signaling pathways was also observed in non-gonadal tissues, which contradicts the common dogma that downregulation of these pathways is associated with lifespan extension. This is a remarkable result that suggests that short-lived model organisms may not correctly reflect signaling effects required for longevity in long-lived species. The authors also point out to higher glucocorticoid levels in workers and speculate that may be showing Cushing's like syndrome.

Overall, this is an important study of high interest.

Items to address:

1) It is not clear how long were the animals maintained in a breeder status prior to analysis.

Yes, we added the detailed information for every animal in the study to the supplement and summarized it in the manuscript as follows:

“In total, females/males selected as breeders spent on average 513/353 (F. mechowii) and 1095/997 (F. micklemi) days in this state before sampling (Supplementary file 1D,E).”

2) Figure 3A, the rational for the comparison of changes in young breeder versus non breeder animals to changes that occur with age over time is unclear.

We thank you for the question and introduced the two corresponding paragraphs with the rationale:

“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.”

“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…”

3) Figure 4, were the changes driven by gonads mainly? Hoe many were non-gonadal? For example, in Figure 2, muscle showed 0 DEGs by status, but in pathway analysis we see upregulation of pathways. Please explain.

The enrichment analysis is not based on the differentially expressed genes only but was performed in a threshold-free manner. This was done in order not to lose the ability to detect differences at the pathway level for which the statistical power at the single gene level is insufficient due to the available sample size. An enrichment analysis that would focus only on DEGs weights all genes with a p-value below a certain threshold (e.g. 0.05) with 1 and all genes above this threshold with 0. Instead, we have adopted an approach that takes each gene of a pathway into account in the test statistics – weighted according to its p-value, providing test sensitivity in tissues with no or few detected DEGs. We have addressed this approach in the Results section of our manuscript as follows:

“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. We used an established method that combines all p-values of genes in a given pathway in a threshold-free manner. 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 section for details).”

To avoid confusion when viewing the figure in isolation or when the paragraph above has been overlooked, we have also added the following to the figure description:

“This enrichment analysis was carried out threshold-free, which is why tissues without differentially expressed individual genes (see Figure 2B) can also show differentially expressed pathways.”

4) Table 2: TOR is downregulated in breeders while ribosome processed are upregulated, please explain.

It should be noted that the mTOR downregulation, as shown in Table 2, is limited to the adrenal gland only. Nevertheless, some interesting questions do indeed arise with regard to mTOR regulation due to the scenario on the one hand (massive life span differences) and the concrete results on the other hand (mTOR largely unchanged, but protein biosynthesis up-regulated). Although we cannot conclusively answer these questions, we will discuss them in a new paragraph:

“An important positive regulator of ribosome biogenesis and protein synthesis is mTOR (Johnson et al., 2013). Based on this, one could expect an upregulation 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 downregulated 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.”

5) The finding of Cushing's syndrome needs more support. Please consider comparing Cushing's related transcriptome changes as a whole to the changes observed in mole rats, otherwise this conclusion may need to be toned down.

Many thanks for this really good suggestion! We were able to find one single data set that performs an RNA-seq comparison of Cushing vs. controls. We have correlated this data on a global level with our non-breeder/breeder comparison and have extended the Materials and methods and the last section of the Results accordingly. We discuss the new results and the limitations of this comparison together with our other hypothesis tests about the HPA axis and Cushing in a new section. The part corresponding to your suggestion is this:

“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^-36). 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^‑40). 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.”

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

We are glad to see that you smoothed your HPA stress hypothesis according to our suggestions, but further revisions are still necessary to address our previous concerns about data quality and robustness of results have and remove the overinterpretation from the title.

Specifically you need to:

1) Control the batch effect using established packages, such as "BatchQC" or "BEclear";

As suggested, we analyzed potential batch effects systematically using BatchQC, added a respective source data set 2h with all reports for full transparency, summarized them in a paragraph in the Control analyses section and referenced the issue in the first paragraph of the Results.

“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).”

“Control analyses: For practical reasons, samples collected for this project over several years were sequenced progressively in different batches. Therefore, we analyzed our data and sample scheme across all tissues and both species using BatchQC (Manimaran et al., 2016). In summary, the extensive search did not reveal any possible batch effects in the dataset from F. micklemi, although they may not be ruled out completely for F. mechowii (Source data 2h). A source for the minor batch effects in F. mechowii may be attributable to the repeated sequencing of the same sample(s) in different sequencing experiments – an approach avoided during generation of the F. micklemi dataset. However, because we followed a robust multifactorial analysis, in which differences must occur consistently across both species to be considered (see “Differentially expressed genes analysis” in Materials and methods), false-positive results due to batch effects are substantially reduced.”

2) Detect DEG using limma with voom to account for large variances and compare the results with DESeq2's output (see PMID: 23497356);

As suggested, we searched for DEGs using limma/voom. The respective p- and FDR-values obtained with this method are juxtaposed to the respective DESeq2 values in the tables in Source data 1e-g. We compare DESeq2 and limma/voom results across single tissues in the new Figure 2—figure supplement 2 and summarized this cross-validation in the following paragraph as part of the new Control analyses section.

“Control analyses: In addition to detecting DEGs with DESeq2 (Love et al., 2014), we also examined the central contrast, i.e., the juxtaposition of workers and breeders, using the combination of limma (Smyth, 2004) and voom (Law et al., 2014). The absolute numbers of DEGs found in each tissue, and thus the highly uneven distribution, were nearly identical for all tissues except thyroid (Figure 2—figure supplement 2). For the latter, however, a closer look shows that 85% of the DESeq2 DEGs (threshold FDR <0.05) are confirmed if the FDR threshold for limma/voom is raised from 0.05 to 0.1. If an FDR threshold of 0.05 is applied for both approaches across all tissues, a total of 58% of the DESeq2-DEGs are confirmed by limma/voom. Notably, all genes mentioned in the manuscript, including the candidates from Table 1, are among this set. Using a less strict threshold of 0.1 for limma/voom, 91% of DESeq2-DEGs are confirmed.”

Discussion: “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 manuscript were confirmed by this validation. This specifically includes all genes supporting our HPA hypothesis, e.g., 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).”

3) Redo the pathway enrichment with changes of expression values or other static indicators (see PMID: 23625889 and PMID: 32515539 for examples) but not p-value; to remove the bias and validate the HPA hypothesis.

The requested analysis was added to the Control analyses section and illustrated in a new supplemental figure. We refer to this additional check in our Results section. However, we would like to emphasize at the aggregation of p-values for the purpose of pathway enrichment in the context of gene expression studies has been used or even recommended in numerous previous studies (e.g. https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0012693, https://www.tandfonline.com/doi/abs/10.1080/19466315.2014.888013, https://bmcbioinformatics.biomedcentral.com/articles/10.1186/1471-2105-14-368, https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-8-96, https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3632109/).

Results: “As standard approach, we used an established method (Fridley et al., 2010; Evangelou et al., 2012; Poole et al., 2016) combining all p-values of genes in a given pathway in a threshold-free manner (Figure 4). […] 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 section).”

Discussion (see issue 2, above).

Control analyses: “In order to avoid potential biases by using solely p-values as indicator of differential gene expression, we additionally resorted to fold-changes as a test statistic for our pathway enrichment analyses. This exercise roughly confirms half of the pathways/hallmarks identified as differentially expressed at the cross-tissue level. Notably the set includes, consistent with the HPA hypothesis, steroid hormone biosynthesis and ribosomes as a central part of up-regulated anabolism in breeders (Figure 4—figure supplement 11). Furthermore, myogenesis, the P53 pathway, and coagulation were confirmed as cross-tissue differentially expressed by this approach. On the other hand, proteasome and oxidative phosphorylation gene sets, were not identified by the fold-change method.”

We have extended the corresponding Materials and methods section with the respective formulas. Also, other analyses pertaining to pathway enrichments were moved to the “Control analyses” section.

4) Reconsider the title, which is still an over-interpretation of the data and therefore needs to be fixed as follows. The words "slower aging" imply that the same biological aging process is occurring in the breeders and non-breeders, which is not actually clear from the data. Indeed, it seems quite plausible that the non-breeders die prematurely due to too much stress (as you speculate). Whether the rate of aging is different or not remains unclear. Replacing "slower aging" with "Increased longevity" or something like that would address this issue.

As suggested, we changed the title to “Increased longevity due to sexual activity in mole-rats is associated with transcriptional changes in the HPA stress axis”.

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,295

  Page views

• 240

  Downloads

• 14

  Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.