Conserved aging-related signatures of senescence and inflammation in different tissues and species

Increasing evidence indicates that chronic inflammation and senescence are the cause of many severe age-related diseases, with both biological processes highly upregulated during aging. However, until now, it has remained unknown whether specific inflammation- or senescence-related genes exist that are common between different species or tissues. These potential markers of aging could help to identify possible targets for therapeutic interventions of aging-associated afflictions and might also deepen our understanding of the principal mechanisms of aging. With the objective of identifying such signatures of aging and tissue-specific aging markers, we analyzed a multitude of cross-sectional RNA-Seq data from four evolutionarily distinct species (human, mouse and two fish) and four different tissues (blood, brain, liver and skin). In at least three different species and three different tissues, we identified several genes that displayed similar expression patterns that might serve as potential aging markers. Additionally, we show that genes involved in aging-related processes tend to be tighter controlled in long-lived than in average-lived individuals. These observations hint at a general genetic level that affect an individual’s life span. Altogether, this descriptive study contributes to a better understanding of common aging signatures as well as tissue-specific aging patterns and supplies the basis for further investigative age-related studies.


INTRODUCTION
The phenomenon of biological aging is still far from satisfactorily explained due to its highly complex nature. Aging is often described as a decline in cellular functions over time and the cause of several severe diseases, such as cardiovascular diseases, neurodegenerative diseases or cancer [1]. Because aging is incompletely understood, a decades-old and still ongoing debate exists on the true source of aging, giving rise to a variety of competing theories. Senescence and inflammatory processes are two of the most common themes within these discussions of the molecular driving forces of aging. Cellular senescence is a state in which permanent replication is halted, and thus cells are unable to further proliferate, and their overall function is strongly diminished [2]. Senescence is meant to be a protective mechanism, stopping further proliferation if cells are on the verge to turn into malignant tumor cell due to severe DNA damage, for example because of telomere AGING shortening. During an organism's lifetime, most cells continuously undergo proliferation and cell division and reach the state of senescence at their own pace, when their telomeres have reached a certain shortening threshold. Over time, this process leads to an accumulation of senescent cells accompanied by loss of function and integrity of the respective tissues, which reflects the close connection of senescence with aging [3]. Recently, cellular senescence was described as the "nexus of aging" by Bhatia-Dey et al., suggesting it as the main driver of the aging process [4].
Another systemic process, that is observed in most tissues with age is the increased release of proinflammatory messenger substances. As a consequence, low-grade chronic inflammatory processes slowly but irresistibly begin to damage organs and are viewed as the cause of other age-related chronic diseases, such as Alzheimer's disease, osteoporosis or diabetes [5,6]. This state of chronic age-dependent inflammation is also suspected as one of the main causes of biological aging and was described as "inflammaging" [7]. Miquel et al. proposed an integrative oxidation-inflammation theory of aging, arguing that chronic oxidative stress originating from mitochondria leads to senescence in cells of the regulatory systems, such as the immune system [8].
Unarguably, senescence and inflammation processes are strongly connected and contribute to an organism's aging phenotype as well as its rate of aging. Therefore, understanding how those systemic processes are regulated during aging in different species and tissues could supply many answers related to biological aging itself.
Here we present a descriptive transcriptomic study on the age-dependent genetic changes of senescence and inflammation in the four evolutionarily distinct species of Homo sapiens, Mus musculus, Danio rerio and the short-lived fish Nothobranchius furzeri and up to four different tissues (brain, blood, liver, and skin). Our aim was to identify potential markers of aging across species and tissues by comparing a young mature time point against an aged and old-aged time point. We report 26 different genes that showed consistent upregulation or downregulation towards old age in multiple tissues and discuss their roles in aging. Furthermore, we identified several genes that were similarly regulated during aging among the investigated species in a tissue-specific manner. Additionally, we observed a stricter control of gene expression of agingrelated processes in the rather old-aged individuals compared with the aged individuals. We conclude that certain of the identified genes that show a conserved age-dependent expression pattern are potentially interesting targets for therapeutic developments designed to achieve healthier aging.

Gene expression discriminates among tissues more strongly than ages
To evaluate the homogeneity of the investigated transcriptomic data, we performed t-SNE dimensionality reduction of the RNA-Seq libraries based on the measured expression strengths of our 464 preselected senescence-and inflammation-associated genes (see Figure 1). We observed a homogeneous clustering of the samples with respect to species and tissue, with almost no segregation of the three different ages within these clusters. This observation was already made and reported in similar studies based on different sets of genes and indicates that aging is a relative subtle process, at least on the transcriptional level [9]. Nevertheless, separations between the youngest and oldest time point can be observed within the skin clusters of all four species but to samples were clustered based on the expression patterns of the selected senescence and inflammation related genes, utilizing the t-SNE approach. All tissues form distinctive species-specific and non-overlapping clusters with exception to very few single outliers. Additionally, larger species-independent tissue clusters were drawn to improve the visualization of the data. However, the three different time points did not generally separate in independent clusters of their own. A weak segregation can only be observed among the mature and old-aged skin samples of all four species. For more information, see Supplementary Data 10. AGING a much weaker extent in the human samples. This observation might suggest a pronounced and conserved difference in the activity of senescence and inflammation processes during skin aging of evolutionarily distinct species. This observation and the fact that the skin is also relatively easily accessible makes it especially interesting for interspecies comparisons of age-related senescence and inflammation. However, it remains to be further investigated whether this observation can be made in similar transcriptomic studies.
Another observation is that although the different tissue samples form species-specific clusters, they are still well separated from each other, i.e., forming tissue-specific clusters, with the exception to the Homo sapiens and Mus musculus blood samples. This observation hints that the similarities of the inflammaging process among the same tissues of different species prevails over any systematic effect within one individual.

Expression variance is more controlled in long-lived individuals
It is proposed that individuals that have reached a comparatively high age have somehow counteracted the effects of inflammaging by unknown anti-inflammation processes [5,6]. As one approach to validate this hypothesis, we analyzed how the mean variance in gene expression changes during aging. Interestingly, we found a significantly lower relative standard deviation of the old-age time points compared with the aged time points in almost all species and tissues and even lower than the mature time points in selected cases (see Figure 2 and Supplementary Data 4). This observation indicates that the selected genes associated with inflammation and senescence have a more stable expression, suggesting that both age-driving processes are in fact more controlled in the long-lived individuals. Similar observations were recently made in brains of long-lived mice [10]. Exceptions are the skin samples of all species other than humans, which show the opposite pattern of a significantly increased variance in gene expression in the old-age time point. The skin is the most diverse of the tissues studied in this work, with mice and humans having fur and hair, whereas fish have scales. In addition, it is one of the most exposed tissues and therefore faces a potentially higher diversity of factors that affect ageing in skin [11]. As a result, skin cells are subjected to more and different stresses, and the rate of senescence is increased compared with other tissues, which also affects the rate of inflammation within the skin tissues due to the pro-inflammatory secretory phenotype that senescent cells can develop [12]. In contrast, human skin showed a decreased variance in gene expression, indicating less pronounced or at least tighter controlled inflammatory and senescence processes. Such a distinct behavior could be the result of various intrinsic or extrinsic physiological factors like the thinning of the epidermis, the degeneration of the extra cellular matrix, or a decline of the microvasculature due to a reduced angiogenic capacity [24]. Besides, additional factors like external applications of antioxidants on skin might also play role in it [11].
The only other observed exception is the brain of Nothobranchius furzeri, which displayed a strongly increased gene expression variance at the old-age time point compared with both younger time points. This pattern is similar to that in fish skin and is already reflected in the sample clustering as mentioned above (see Figure 1). This observation could indicate a weaker protection of the Nothobranchius furzeri brain against inflammaging, but it remains to be further studied.
To verify that the more controlled gene expression observed in long-lived individuals is not biased due to our preselected genes, we repeated the same analysis with all expressed genes in every species and tissue. We confirmed our observation with all tissues showing the same relative variance changes as before, in principal. The overall variance within each age group was decreased, which was expected because many more genes were included in the analysis that show no change in expression during aging (see Supplementary Data 4). Nevertheless, in order to confirm that longlived individuals really benefited from controlled gene expression, further investigations, such as longitudinal studies, are necessary.
Oxidative stress response tends to be the prevailing process towards old age To examine whether certain biological functions predominantly drive the inflammaging processes in a tissue-or species-dependent manner, we analyzed the identified differentially expressed genes (DEGs) with respect to their molecular function. We found that most of the DEGs within the brain and liver samples belong to immune and inflammatory response processes, whereas changes in the expression of senescence-related genes are predominately related to skin aging (see Figure 3, Supplementary Data 6 and Supplementary Data 11). Only few DEGs could be identified in blood age comparisons for humans and mice, making it difficult to interpret these results in a meaningful way.
It is a known and common observation of aged skin that it is more susceptible to infections, physical damage and reduced epidermal barrier integrity as well as other agerelated deficiencies [14]. Epidermal stem cells maintain the tissue's homeostasis and loss of those stem cells due to premature senescence is the main cause of aging within the skin [15]. Our data confirm this observation, showing mainly senescence-related DEGs within the skin of all four species already in the early and the late aging comparisons. However, whether the source of epidermal stem cell senescence is primarily intrinsic or extrinsic factors is controversially discussed [16,17]. Apoptosis appears to play a more important role during aging of the liver than in any of the other investigated tissues. This is a confirmed observation, because liver homeostasis is mainly regulated through apoptotic processes and many liver dysfunctions and diseases are related to apoptosis [18,19].
Most interestingly, we observed that DEGs associated with the oxidative stress response and oxidative processes occur predominantly within the late age and longevity comparisons (i.e., comparisons with the longlived individuals) of the brain, liver and skin samples. This observation suggests that those individuals that grew older than their species average survival age have either a generally more active or better regulated response mechanism against reactive oxygen species (ROS). Free radicals such as ROS and their great impact on the aging process are much discussed [8,20]. During respiration of oxygen, free radicals are produced as a harmful byproduct that can oxidize macromolecules such as DNA or proteins and damage them in this process. The main point of aging theories involving oxidative stress is that with age, the rate of harmful ROS accelerates because ROS first target the mitochondria (because they are mainly generated in that location), leading to an accumulation of damage in these cell organelles. The more degenerated the mitochondria become, the higher the rate of released ROS becomes, injuring genomes and membranes of neighboring cells and, as a final after-effect, resulting in their senescence or apoptosis and therefore aging [21]. Different antioxidant mechanisms exist within the cells to prevent the release of free radicals and protect the cells from such deleterious consequences [22]. It was observed that long-lived species appear to have higher protection against ROS than short-lived species [23,24]. Our data suggest that the same is true for long-lived individuals within one species, confirming the importance of oxidative stress response for longevity.

Conserved aging expression signatures across tissues and species
To identify a more precise signature of the abovementioned observations (Figure 3), we analyzed the expression patterns of our preselected gene set in greater detail. We were particularly interested in genes that show a constant increase or decrease towards old age that could also be consistently observed between the investigated tissues and species.
We found 16 genes that were consistently higher expressed and 10 that were consistently lower expressed during aging in at least three of the four species and tissues, totaling 26 genes that show a conserved expression pattern with aging (see Figure 4). Certain of these genes have already been reported in the context of aging, and we observed them to be the most conserved genes from our preselected set among species and tissues. Hence, we discuss several of these genes based on their biological function in additional detail. More senescence-and inflammation-related genes were found to have conserved expression patterns but to a lesser extent and can be found in Supplementary Data 7.

Immune/inflammatory response
The gene marco encodes for a scavenger receptor (Marco -macrophage receptor with collagenous structure), which is typically found on macrophages but also other immune cells [25,26], and acts in the innate immune system by binding and clearing pathogens and initiating an inflammatory response [27,28]. We observed marco to be constantly upregulated (mostly by an increase of two-fold or more) towards old-age in all four species and within the blood, liver and skin. Although Marco does not cause inflammation directly, it is important for the activation of other receptors, such as members of the toll-like family [28]. These receptors subsequently activate NF-κB, inducing an inflammatory response. Interestingly, Marco is associated with Alzheimer's disease, because it is also present on microglia and a decreased response of these cells has been observed after binding amyloid beta peptides [29]. The coagulation factor II receptor (PAR1) is encoded by the gene f2r, which similar to marco shows persistent upregulation with age in blood, liver or skin within all four species. The major function of PAR1 is mediation between coagulation and inflammation activity, and therefore it has a key role in inflammatory response activation [30,31]. In addition, PAR1 activity is implicated in aging-associated cardiovascular diseases [32].
Four more genes (cd40, sh2d1b1, ptafr, adam8) coding for membranous receptors have a common upregulated signature through several aging tissues and species. These genes have all been shown to promote proinflammatory signaling, either directly by binding platelet-activating factors (ptafr [33]), activating antigen presenting cells (cd40 [34]), and acting as a regulator of antigen receptor signal transduction (sh2d1b1 [35]) or more indirectly by releasing and degrading other cell surface receptors of leukocytes (adam8 [36]).
Observation of a generally higher expression of ceruloplasmin encoded by the gene cp in all species and tissues with age (except Nothobranchius furzeri brain) is not unexpected due to its role as an acute-phase protein.
These reactants are increasingly released in the blood plasma during an inflammatory response and support the innate immune response [37]. Thus, increased cp expression is most likely a result of the chronic inflammatory nature of aged tissues.
Lef1 is one of the few immune-related genes that showed a decreased expression with age in the blood and skin of human and mice respectively and the brains of both examined fish. The main function of the LEF1 (lymphoid enhancer-binding factor 1) protein is to enhance expression of the T-cell receptor alpha chain [38]. However, a variety of different interaction partners are known, linking LEF1 as an important regulator within the WNT and TGF pathways and implicating a role in apoptosis and cell proliferation of leukocytes [39,40].

Oxidative stress response
As previously mentioned, ROS and oxidative stress are widely believed to be the driving force behind the cellular aging process [8]. The source of oxidative stress is likely an imbalanced activity of the respiratory chain and antioxidant processes. We observed a constant agedepended upregulation of the genes cyba, cybb, cyp26b1 and nrros in various combinations of the investigated tissues and species. The first two genes encode for the light and heavy chains of the Cytochrome b-245 protein, a superoxide-producing subunit of the NADPH oxidase [41]. By producing and releasing ROS, phagocytes use Cytochrome b-245, mainly as an antimicrobial strategy during an infection. However, the role of Cytochrome b-245 during aging is still debated [42,43], and it is associated with several degenerative diseases due to uncontrolled production of ROS [44,45]. Cyp26b1 encodes another cytochrome protein, Cytochrome P450, which is involved in many metabolic reactions, especially the oxidation of NADPH [46]. Higher expression of cyp26b1 during aged time points is suggested to aid in degradation of toxic substances that have accumulated with age [47]. The genes nrros and gpx8 both code for proteins that belong to antioxidant mechanisms controlling the production of ROS. NRROS directly interacts with Cytochrome b-245 and mediates its degradation, limiting the rate of ROS production of the associated NADPH oxidase complex [48]. Similarly, the enzyme GPX8 protects cells from ROS-induced damage by its peroxidase activity by reducing free hydrogen peroxide to water [49]. However, whereas gene expression of nrros rises with age, most likely to regulate the increased Cytochrome b-245 activity, the translation of gpx8 is downregulated in many tissues and species. We observed the opposite regulation during aging only in the human skin and mouse liver samples. Nevertheless, gpx8 serves as an interesting potential therapeutic target and it was previously shown in the Caenorhabditis elegans model that the deletion of several GPx family members is the cause of an accelerated aging process that results in shorter life-span [50].

Senescence and apoptosis
Most of the 26 identified genes with conserved expression patterns in aging, are strongly related to the cell cycle control or apoptotic processes. Both replicative senescence and apoptosis are immediate causes for the loss of somatic and stem cells with age and can be similarly triggered by chronic inflammation and oxidative stress [51,52].
We found four of these genes to encode direct inducers of senescence: ccnf (Cyclin-F) [53], scrib [54]. cdkn2b (p15 INK4b ) [55] and rassf4 encoding an inhibitor of Cyclin D1 [56]. We found all of these genes to be expressed in a senescence-promoting fashion towards aging (rassf4 and cdkn2b upregulated; scrib and ccnf downregulated) in the examined tissues and species. Additionally, selected genes interacting with cell cycle regulators also showed persistent expression towards old age, namely, dab2 [57], lgals1 [58], and brat1 [59]. The proteins CDCA8 and Survivin encoded by cdca8 and birc5, respectively, are both associated with apoptosis. By inhibiting caspase activity, Survivin is a repressor of apoptosis, whereas CDCA8 directly interacts and stabilizes Survivin [60,61]. However, both proteins show converse directions of regulation during aging between the same tissues and within individual species, resembling rather weakly conserved expression signatures. AGING

Potential common tissue-specific marker genes of aging
We also focused on senescence-and inflammation-related genes that change their expression during aging consistently and similarly in a tissue-specific manner. These genes should present conserved age-related changes and might be the driving factors of tissuespecific aspects of the inflammaging process, probably closely linked to the specialized function of the respective tissue. In a first comparison of the identified DEGs of the different species, we observed only a few common genes that occur as differentially expressed in all species of a single tissue (see Figure 5 and Supplementary Data 9). However, the blood and liver samples revealed no DEGs to be common in all of the respective species, and two were identified within the brain (cybb, cd68) and another two within the skin (Gpx7, Tnfrsf25). Cd68 is a member of the cell surface receptors of macrophages and other monocytes and is a well-known marker for macrophage activation. The Cytochrome b-245 subunit encoding gene cybb and its potentially harmful role during aging were already introduced in the previous section. Tnfrsf25 (a member of the TNF receptor superfamily) stimulates the proliferation of T cells or can initiate apoptotic signaling, leading either to survival or cell death, and is therefore an important regulator of T cell development [62]. Gpx7 encodes for another glutathione peroxidase and has the same antioxidant function as its above described homolog gpx8.
However, none of the four genes exhibits features that serve as a tissue-specific marker because they do not display a consistent up-or downregulation during aging within each investigated species or share the same expression pattern with other tissues. Therefore, we applied a similar but more stringent filtering process to reveal genes displaying the characteristics of potential tissue-specific markers. With this approach, we identified several genes that showed a constant age-related transcriptional increase exclusively in the liver or skin samples of the investigated species.

Liver
Four constantly upregulated genes were identified in the liver samples of Mus musculus, Danio rerio and Nothobranchius furzeri: jag2, anax1, ralb, and sfrb. The protein encoded by jag2 (Jagged-2) is a known activating ligand of Notch2, which is mainly involved in many different developmental processes regulating cell fate decisions [63]. Most interestingly, jag2 overexpression is not only shown to play a critical role in the formation of plasma cell myeloma and progression of other tumors [64,65] but also induces the secretion of interleukin-6 [66]. Interleukin-6 is a potent stimulator of immune and inflammatory responses and acts in the development and progression of many age-associated diseases, such as Alzheimer's disease [67], atherosclerosis [68], diabetes [69] and various cancers, making it an attractive therapeutic target [70,71]. Because the liver shows strong signs of inflammation during aging (Figure 3), overexpression of jag2 might be one conserved driving factor of inflammatory processes within hepatic cells and could be of interest as an additional therapeutic target. In contrast, Annexin A1, the protein encoded by anax1, has distinct anti-inflammatory and protective properties, by inhibiting NF-κB signal transduction and counterregulating pro-inflammatory signals in a variety of immune cells [72,73]. The age-related upregulation of Annexin A1 in the liver could be a component of a tissue-specific mechanism attempting to cope with the chronic state of inflammation. SFRP2 (encoded by sfrp2) has an oncogenic character and is associated with cancer formation by acting as an inhibitor of the canonical WNT/βcatenin pathway [74,75]. The last of the four conserved liver-specific upregulated genes, ralb, encodes for one of two Ral protein paralogs (Ral-B) and acts as a major modulator of a multitude of cellular processes [76]. Although many functions are shared between the paralogs Ral-A and Ral-B, the latter is more specifically involved in the activation of apoptosis [77], which we observed to occur more predominantly during liver aging (Figure 3). Additionally, Ral-B has clinical significance because it promotes tumor progression of several cancers and activation of the innate immune response [76,78] and hence might contribute to inflammatory stress in the liver.

Skin
If including all four species, we did not observe any gene displaying the characteristics of a potential conserved marker gene for skin aging. Considering the diverse nature of skin among the studied species we decided to examine the skin in two separate groups of, fishes and mammals.
Between human and mouse, only pecam1 was identified to be consistently upregulated towards old age. The gene encodes the protein Pecam-1, also known as CD31, which is a common immunohistochemistry marker used to evaluate tumor growth and has a major role in the removal of old neutrophils, which could explain its higher expression with age. However, this gene is already known as a biomarker of inflammatory processes [79].
The two fishes Danio rerio and Nothobranchius furzeri share 6 genes that showed a persistent increased expression with aging exclusively in their skin and their encoded proteins are all involved in mitotic progression of the cell cycle: ckap2 [80], cdca8 [61], aurkc [81], cep55 [82], mis12 [83] and spc25 [84]. This observation might indicate that fish skin cells are not as affected by senescence as the skin of humans and mice, which is in line with the observation that the skin of zebrafish shows a high regenerative capacity even at higher age [85].

CONCLUSIONS
In the current descriptive study, we have investigated the inflammation-and senescence-related gene expression during aging in multiple tissues of four evolutionarily distinct species. By analyzing the gene expression profiles, we were able to identify common signatures of aging, i.e., genes that illustrate steady increased or decreased expression with age. We identified 26 genes that shared an age-dependent expression pattern in at least three of the four investigated species and at least three different tissues. These genes represent interesting targets for further study, because they hint at general molecular mechanisms of aging that occur similarly not only within various organs of one organism but also between different species. It is important to note that we focused our study on a specific set of senescence-and inflammation-related genes and other processes might exist that share similar expression changes during aging across distinct species.
Several of the identified genes are not only directly involved in the initiation of an inflammatory response but also play a major role in sustaining the state of inflammation. Modulation of the expression level of certain these genes, such as marco and f2r, or addressing the respective proteins could be one potential approach to controlling chronic inflammation and thus might result in reduced inflammaging-related cellular stress with aging. In total, six different genes coding for macrophage cell surface receptors (marco, f2r, cd40, sh2d1b1, ptafr, adam8) have been found to show conserved upregulated expression across different tissues and species. Accumulation and higher activity of macrophages has already been linked to aging and removal of macrophages has shown beneficial effects for age-related disorders, such as neuro-degeneration, or atherosclerosis [86][87][88].
Another strongly aging-associated intrinsic stress factor is the production of cytotoxic ROS as a byproduct of the respiratory chain. Although phagocytes use Cytochrome b-245 (a heterodimer encoded by the genes cyba and cybb) to produce ROS to kill microbes, excessive generation of ROS causes premature replicative senescence of cells due to DNA damage [89]. We observed a conserved age-dependent upregulation of selected genes involved in ROS production (cyba, cybb, cyp26b1) together with a similarly upregulated but apparently inefficient (nrros) or even downregulated oxidative stress response (gpx8). Most of the identified genes with conserved expression patterns were related to senescence and apoptosis and showed an increase in both processes with age in general. As noted by many different theories on aging, each of the mentioned cellular processes is suspected to be the main driving force behind biological aging. However, it has become increasingly apparent that these processes should not be treated as separate and seemingly competing sources of aging, because they are strongly interconnected [89]. During a state of inflammatory stress, excessive production of ROS is enhanced due to dysfunction of the respiratory chain in mitochondria, leading to an accumulation of senescent cells, which can develop a pro-inflammatory secretory phenotype [12]. In addition, inflammation is a known trigger of apoptotic processes and during the process of apoptosis, certain cellular components are released that can further activate inflammatory processes, adding to the overall cellular stress [90]. As a consequence, an environment of chronic AGING inflammatory, oxidative, senescent and apoptotic stress is established, with these processes mutually triggering each other and boosting the cycle of self-harming development [90] over time. Additionally, constant stimulation of cellular stress response mechanisms promotes genetic deregulation, which is reflected in the observed gene expression changes.
Based on the annotated functions of the identified DEGs, we observed that the individual tissues express different aging-related processes more strongly than others. Although the impact on age-related expression changes of immune and inflammatory response processes was stronger in the blood, brain and liver samples of all investigated species, the skin samples displayed greater modulation of cell cycle and senescence-associated genes. This observation is most likely due to the specialized functions of these different organs. Only few tissue-specific and conserved regulated genes could be observed in the liver and skin and none in the blood and brain samples. Nevertheless, the identified genes potentially represent molecular sources for why certain aging-related processes appear more tissue-specific than others.
Additionally, we report that gene transcription in longlived individuals is generally more controlled compared with average-lived individuals, showing significantly lower variance in gene expression in all tissues and species, except the skin. Maintaining a more stable transcriptional activity, not only for senescence and inflammatory response processes, appears to have a significant life-prolonging effect [91]. It is also considered that oxidative stress response processes were observed to be more regulated in the old-age comparisons, and a possible reason could be better management of oxidative processes and more efficient antioxidant mechanisms. However, this conclusion remains speculative until further experimental studies can prove this observation.
Overall, in our work, we describe previously unknown conserved transcriptional changes across different species and tissues as well as tissue-specific changes with age, supplying a complementary overview of how changes in gene expression relate to processes of aging. In addition, our findings could serve as a basis for new strategies in the development of therapies against agingrelated diseases.

High-throughput transcriptomic data processing
For each species we investigated up to four different tissues at three different ages: young mature (M), aged (A) and old-aged (OA) (see Figure 6).  [92] was used to trim and filter the RNA-Seq libraries, by clipping reads on either ends to achieve a minimum base quality score of 20 and discarding all reads with a length of less than 15 nt or containing more than two ambiguous N bases. Second, quality-trimmed RNA-Seq libraries were mapped onto the according recent genomes (as from Ensembl release version 92 [93] and [94,9]) using TopHat2 (v2.1.1) [95] under consideration of the default parameters. This procedure allowed the identification of spliced reads, and the mapping of single reads to multiple best-fit locations. Third, read counting was performed by Featurecounts (v1.5.3) [96], with the counts normalized to transcripts per million (TPM) [97] to discard all genes with a value of ≤ 1 in every sample from subsequent expression analyses. Finally, differentially expressed genes were identified via the DESeq2 (v1.10.0) [98] Bioconductor package by comparing the three different time points of each species and tissue. False discovery rate adjustment of the resulting gene's p-values was performed according to previous work [99]. Details of all DEG results, together with the raw and normalized count values are given in detail in Supplementary Data 1.

Collection of genes of interest and t-SNE dimensionality reduction
An initial list of 769 genes relevant in the context of senescence and inflammaging was compiled from a detailed literature search (see Supplementary Data 2). From those 769 genes, we removed all genes that were either missing an ortholog in any of the four species or were not differentially expressed (FDR adjusted p-value < 0.05) at least once in any species and tissue, ending with the final list of 464 genes (see Supplementary Data 3). Based on the measured expression of these 464 genes the t-SNE algorithm [100] was used for dimensionality reduction in order to investigate the clustering behavior of the analyzed RNA-Seq libraries (see Figure 1).

Gene expression coefficient of variation analysis
To estimate the general gene expression variation among the different time points of each species and tissue, we calculated the coefficient of variation for every expressed gene (i.e., a TPM value greater than 1) as the ratio of its standard deviation to its mean. The significance of the variation differences among time points was determined using the twosided t-test. We repeated this analysis using only the 464 genes of interest (for details, Supplementary Data 4).

Gene function annotation and expression pattern analysis
The exact biological functions of our preselected gene set of interest were obtained from the functional annotation database David (version 6.8) [101] (Supplementary Data 5). For every age comparison ( Figure 6) we determined the dominant biological processes by the majority of the annotated functions of the respective DEGs (for details, see Supplementary Data 11). To identify common signatures of senescenceor inflammation-related processes, we filtered for genes that were differentially expressed at least once within one of the age comparisons and showed the same direction of up-or downregulation by a minimum threshold of 10% in every other age comparison. This rather low effect-size threshold can be justified, because it was previously shown that expression changes related to aging tend to be more subtle [102]. This process was conducted for every species and tissue individually. Based on this gene subset, we further searched for common expression patterns of single genes among the different tissues and species. When searching for tissuespecific signatures, we applied the same filtering strategy but raised the minimum threshold of expression changes to 25%.

Additional information
All relevant data can be found within the paper and its supporting information files. The online supplement is integrated into the open science framework (https://osf.io/) and available under https://osf.io/kzq5y/ or https://doi.org/10.17605/OSF.IO/KZQ5Y. As soon as this manuscript gets published, the Supplementary Data will be made publicly available and the corresponding link and doi will be provided.

ACKNOWLEDGMENTS
We thank U. Teschner for the coordination of human sample material acquisition. We also thank I. Goerlich and I. Heinze for their skillful and diligent technical assistance in NGS data acquisition, M. Guenther for tissue-sampling and isolating RNA from mice and M.