Introduction

Aging is complex biological phenomenon that results in functional decline across tissues and organs systems (1). Age is one of the most significant contributors to disease risk (2), and the aging process is influenced by a variety of genetic and environmental factors at the tissue and organismal level (3–5). The application of omics technologies in recent years has led to unprecedented insights into the complex biology of aging (6). Bulk and single-cell transcriptomic analyses of mouse tissues across the lifespan have shown that aging tempo and trajectory, as indicated by tissue transcriptomic signatures, varies widely across tissues (7–9). Remarkably, many of the transcriptomic signatures of aging across tissues can be significantly reversed by caloric restriction or rejuvenated by the transfusion of young blood (10, 11). Large-scale proteomics analyses of the human plasma proteome also reveal distinct waves of changes across the lifespan, which are associated with distinct biological pathways and affect age-related phenotypic traits and diseases (12, 13).

Despite enormous progress, there is no consensus regarding the mechanisms underlying aging at the cellular and tissue level. Although many non-mutually exclusive hypotheses have been put forth to explain the root cause of aging—oxidative damage, genomic instability, epigenetic changes, loss of proteostasis, mitochondrial dysfunction, DNA damage, telomere shortening, cellular senescence, stem cell exhaustion—it remains a major challenge to distinguish between the driver and passenger mechanisms of aging (14). Nevertheless, efforts to understand the proximal and ultimate cause of aging will facilitate development of therapeutics aimed at improving healthy aging (15, 16).

In the present study, we focused on aging from a mitochondrial perspective, as this organelle is known to play an important role in the aging process (17–20). Mitochondria supply the bulk of the energy needed to maintain tissue health and repair tissue damage, and their function tends to decline with age. Over time, damage accumulates in mitochondrial DNA, proteins, and lipids, which compromises their functional integrity and leads to dysregulated metabolism and increased oxidative stress. Recent transcriptomic analyses have highlighted major reductions in electron transport chain genes across the lifespan (7), and these changes can be significantly reversed by the transfusion of young blood into an old mouse (10).

Aging-associated reduction in mitochondrial OXPHOS genes appear to be conserved between human, mouse, fly, and worm (21). Accordingly, mitochondrial dysfunction has been implicated in various age-related diseases, including neurodegenerative disorders, cardiovascular diseases, and metabolic syndromes (22). Boosting mitochondrial health has been shown to delay age-related decline in organ function (23–26).

Given the central role of mitochondria in tissue health, we aimed to address the extent and magnitude of aging-induced changes in mitochondrial function across tissues and organ systems. Although many studies have examined mitochondrial respiratory capacity in various tissues, the scale was limited in that only a very small number of tissues could be interrogated at the same time. This is largely due to the inherent low-throughput method of assessing respiration which requires freshly isolated mitochondria or cells from tissues (27). Consequently, it was not previously feasible to have a comprehensive and systems-level analysis of mitochondrial function across many tissues and the lifespan.

This barrier, however, has been recently overcome. An innovative method by Acin-Perez and coworkers has made it possible now to assess mitochondrial function in previously frozen tissues (28). This new method circumvents the requirement for the isolation of mitochondria at the time of tissue harvest, allowing a large number of tissues (e..g., from a mouse) to be collected, frozen, and assayed at a later time. We adopted the new method in a standardized workflow to profile mitochondrial activity in 33 tissues from young and old mice of both sexes. The dataset consists of a total of 1320 tissue samples from 40 mice and 3960 high-resolution respirometry assays encompassing three technical replicates. Our study represents the largest and the most comprehensive tissue respirometry analysis to date. Our data provide an unprecedented view on the variations and changes in mitochondrial functional capacity across tissues, sex, and age, thus informing ongoing studies on the causes and consequences of aging.

Results

Pan-tissue mitochondrial respiration atlas overview, work-flow, and analysis pipeline

To assay mitochondrial respiration and its maintenance across age and sex, we collected 33 tissues from young (2.5 months; ∼18-year-equivalent in human) or old (20 months; ∼65-year-equivalent in human) male or female mice (n = 10 mice per age and sex). The tissues collected include different brain regions (hippocampus, cortex, cerebellum, and hypothalamus), different sections of the GI tract (stomach, duodenum, ileum, jejunum, cecum, proximal colon, and distal colon), various fat depots (gonadal, inguinal, mesenteric), different skeletal muscle groups (tongue, diaphragm, quadriceps complex, hamstrings, gastrocnemius, plantaris, and soleus), reproductive organs (testis and fallopian tubes), as well as liver, pancreas, heart atria and ventricle, spleen, kidney cortex and medulla, eye, and skin (Figure 1A). Our aim was to provide a comprehensive systems-level view of mitochondrial respiration across tissues, sex, and age.

Figure 1.

Our standardized workflow involved thawing frozen tissue samples in buffer, minced, and homogenized. Samples were centrifuged to pellet cell debris and the supernatant was collected for immediate protein and mitochondrial content quantification (via MitoTracker Deep Red, MTDR). Sample respiration rates were then assayed using a Seahorse XFe96 Analyzer (Figure 1B). The basic respiration assay consisted of four sequential steps: first, baseline unstimulated measurements were obtained. Then, NADH was used to assess respiration via mitochondrial complex I (CI), or succinate was used as to assess respiration via mitochondrial complex II (CII) in the presence of rotenone (Rot, a CI inhibitor). Following this, samples were exposed to rotenone and antimycin A (AA, a complex III inhibitor) to silence respiration. Then, TMPD in the presence of ascorbate was used to assess respiration through mitochondrial complex IV (CIV), via donation of electrons to cytochrome c (Figure 1C). Detailed information of the methodology can be found in the methods section, which closely follows the method first described by Acin-Perez et al (28). High resolution respirometry data for each of the 33 tissues were used in all subsequent comparison (Figure 1-figure supplement 1-33).

Mitochondrial function across different organ systems in male and female mice

The first analysis made was within a group (male, female, young, or old) across all tissues. This allowed us to focus on shared and unique mitochondrial properties across different tissues within a single mouse system. Ranking young male or young female tissues by their respiration via CI (NADH-stimulated), CII (succinate-stimulated), or CIV (TMPD and ascorbate-stimulated) showed that both sexes have the greatest oxygen consumption in the heart atria and ventricles, brown adipose tissue (BAT), kidney cortex and medulla, and the lowest respiration in the colon (distal or proximal), plantaris muscle, jejunum, ileum, and mesenteric white adipose tissue (mesWAT) (Figure 2A-B).

Figure 2.

Ranking old male tissues by their respiration via CI, CII, or CIV showed that the heart atria and ventricles, BAT, kidney cortex and medulla have the highest oxygen consumption, while the ileum, pancreas, skin, duodenum, stomach, gonadal white adipose tissue (gWAT), mesWAT, and inguinal white adipose tissue (iWAT) have the lowest (Figure 2C). Ranking all old female tissues by their respiration via CI, CII, or CIV showed that the heart atria and ventricles, BAT, quadriceps muscles, and diaphragm have the highest respiration, while the pancreas, iWAT, mesWAT, ileum, distal colon, eyes, and skin have the lowest (Figure 2D).

Our analyses revealed a wide range in mitochondrial respiratory capacity across nearly all tissues. We observed tissue-types of high similarity, such as the different regions of the brain (Figure 2-figure supplement 1) and kidney cortex and medulla (Figure 2-figure supplement 2), as well as those displaying degrees of functional heterogeneity, such as the heart atria and ventricle (Figure 2-figure supplement 3), different skeletal muscle groups (Figure 2-figure supplement 4) and white adipose tissue depots (Figure 2-figure supplement 5), and different sections of the gastrointestinal tract (Figure 2-figure supplement 6). Together, these data help illustrate the heterogeneity or lack thereof in mitochondrial function across different tissues and organ systems. Remarkably, even when age and sex are controlled for, mitochondria within a single tissue-classification or organ system can have distinct respiration signatures.

Sex differences in mitochondrial function in young mice

The first pan-tissue comparison of mitochondrial function across groups was made between young male and female mice. In this comparison, we focused on the effects of sex on mitochondrial respiration in young mice. A systems-level view of mitochondrial respiration via CI, CII, or CIV across 32 tissues (omitting reproductive tissues) showed that young male and female mice have similar distributions in oxygen consumption rate (OCR; Figure 3A-C, left panel). Individual tissue-level analysis of respiration through CI showed no significant difference between young male and female mice (Figure 3A, right panel). However, individual tissue comparisons of respiration through CII showed that young males have higher mitochondrial respiration in the stomach and gastrocnemius muscles, and reduced mitochondrial activity in gWAT, duodenum, ileum, and distal colon, relative to young females (Figure 3B, right panel). Respiration at CIV showed that young males have higher mitochondrial activity in the BAT, kidney medulla, hippocampus, and stomach, and reduced mitochondrial activity in the gWAT, duodenum, and distal colon, when compared to young females (Figure 3C, right panel). A relative and global view of the sex-specific differences across young male and female mice showed that males have higher mitochondrial respiration in the stomach, kidney, and skeletal muscle, whereas females have higher mitochondrial respiration in the gastrointestinal tract and different fat depots (gWAT, and mesWAT (Figure 3D).

Figure 3.

Effects of age on male mitochondrial function

The second tissue-by-tissue comparison of mitochondrial function was made across age, where we focused on the effects of age on mitochondrial respiration in young and old male mice. A systems-level view of respiration via CI, CII, or CIV across all 33 tissues showed that old and young male mice have similar distributions of OCR (Figure 4A-C, left panel). However, an individual tissue view showed many differences in response to age. Old males showed a striking increase in mitochondrial activity via CI in the heart atria and ventricles, skeletal muscles (quadriceps, diaphragm, soleus, tongue, plantaris), and mesWAT, and reduced mitochondrial respiration in the brain cortex, relative to young males (Figure 4A, right panel). Viewing respiration through CII also showed that old males have markedly higher mitochondrial activity in the heart, skeletal muscles (quadriceps, diaphragm, plantaris), jejunum, ileum, and mesWAT when compared to young males (Figure 4B, right panel). In contrast, old males had reduced mitochondrial activity via CII in the brain (cerebellum, cortex, hippocampus), BAT, stomach, and eyes when compared to young males. Respiration at CIV showed that old males have markedly elevated mitochondrial activity in the skeletal muscle (quadriceps, diaphragm, soleus, plantaris) and cecum, and reduced mitochondrial activity in the brain (cortex, cerebellum, hippocampus), BAT, liver, spleen, stomach, and iWAT, relative to young males (Figure 4C, right panel). A relative and global view of the age-specific differences across male mice showed that old males have significantly elevated mitochondrial activity in the skeletal muscles, mesWAT, GI tract, and heart, with a concomitant reduction in mitochondrial activity in the stomach, iWAT, gWAT, eyes, and nearly all brain regions (Figure 4D). Together, these data indicate age has a strong effect in modulating mitochondrial respiration in male mice.

Figure 4.

Effects of age on female mitochondrial function

The third tissue-by-tissue comparison of mitochondrial function was made between old and young female mice. Systems-level analysis of respiration via CI, CII, or CIV across all 33 tissues showed no significant differences in OCR distribution between old and young female mice (Figure 5A-C, left panel). Individual tissue comparisons, however, showed many significant differences in response to age. Respiration through CI showed that old females have elevated mitochondrial activity in the heart ventricles, skeletal muscles (quadriceps, diaphragm, soleus, tongue, and plantaris), stomach, and mesWAT, and a concomitant reduction in mitochondrial activity in the cerebellum when compared to young females (Figure 5A, right panel). Respiration through CII also showed that old females have increased mitochondrial respiration in the heart ventricles, skeletal muscles (quadriceps, diaphragm, plantaris), and mesWAT relative to young females (Figure 5B, right panel). In contrast, old females had reduced mitochondrial respiration via CII in the liver, brain (cortex and cerebellum), gWAT, fallopian tubes, duodenum, eyes, and pancreas when compared to young females. Respiration at CIV likewise showed that old females have elevated mitochondrial activity in the skeletal muscle (quadriceps, diaphragm, plantaris) and mesWAT relative to young females; concomitantly, old females had reduced mitochondrial activity in the brain (cortex and cerebellum), liver, gWAT, stomach, lung, eyes, duodenum, and skin when compared to young females (Figure 5C, right panel). A relative and global view of the age-specific differences across female mice showed that old females have elevated mitochondrial respiration in the skeletal muscles, mesenteric fat, and heart ventricle, and a concomitant reduction in mitochondrial activity in the eyes, skin, duodenum, gWAT, iWAT, and nearly all brain regions when compared to young females. Similar to males, these data also indicate age has a strong effect in modulating mitochondrial activity in female mice.

Figure 5.

Age differentially affects mitochondrial function in male and female mice

The fourth tissue-by-tissue comparison of mitochondrial function was made between old male and female mice, where we focused on the effects of sex on mitochondrial respiration in old mice. A systems-level view of respiration through CI, CII, or CIV across all 33 tissues showed no differences in OCR distribution across sex in old mice (Figure 6A-C, left panel). At the individual tissue level, however, old male mice had greater mitochondrial activity via CI in the heart atria, kidney medulla, and soleus muscle as compared to old females. Respiration through CII showed that old males have higher mitochondrial activity in the kidney (cortex and medulla) and pancreas, and reduced mitochondrial activity in the stomach, relative to old females (Figure 6B, right panel). Respiration at CIV showed that old males have higher mitochondrial respiration in the kidney medulla, soleus muscle, and lungs, and reduced mitochondrial activity in the cerebellum and mesWAT, when compared to old females (Figure 6C, right panel). A global and relative view of the sex-specific differences across age showed that old males have elevated mitochondrial activity in the distal colon, lung, eyes, pancreas, kidney, and heart when compared to old females, whereas old females have higher mitochondrial activity in all adipose tissue depots, stomach, and nearly all brain regions when compared to old males (Figure 6D).

Figure 6.

Age has a much larger effect than sex on mitochondrial function

To assess the relative contributions of sex and age to mitochondrial respiration, we first counted the total number of significant differences from all cross-group comparisons made in Figures 3-6, which was 98 in total. Clustering significant differences by the type of respiration affected showed CIV to have the most with 39 differences, followed by CII with 38 differences, and CI with 21 differences (Figure 7A).

Figure 7.

Grouping significant differences by the cross-group comparison in which they occurred showed that young male (YM)-by-old male (OM) have the most significant differences with 37 in total, followed by young female (YF)-by-old female (OF) with 35, YM-by-YF with 14, and OM-by-OF with 12 significant differences (Figure 7B). Clustering the significant difference counts by type, sex (originating from YM-by-YF or OM-by-OF) or age (originating from YM-by-OM or YF-by-OF), we saw that although both sex and age exerted an effect on mitochondrial respiration, it was age that resulted in the highest number of significant changes (Figure 7B). We next sought to quantify the number of tissues significantly affected in at least one mitochondrial parameter per cross-group comparison. There were eight tissues with a difference when comparing young male and female, 10 when comparing old male and female, and 21 tissues when comparing young and old male, or young and old female (Figure 7C). Together, these data indicate that the majority of significant differences observed are the result of age.

Next, we calculated the absolute difference of means and summed the resultant values based on their effect-type, sex (originating from YM-by-YF or OM-by-OF) or age (originating from YM-by-YF or OM-by-OF) (Figure 7D). This analysis offered an added magnitude perspective to our view of changes in mitochondrial respiration across sex and age. It allowed us to see that respiration via CI, CII, and CIV all have a greater magnitude of difference as a result of age as compared to sex. Each graph in Figure 7D shows ranked values (lowest/left to highest/right) within the age and sex groupings with the top tissues colored and labeled with their percent contribution to the total. The highest contributors to the differences in respiration via CI as a result of age were heart (ventricles and atria), BAT, and quadriceps muscles, while those for sex are heart (ventricles and atria) and BAT. The highest contributors to the differences in respiration through CII as a result of age were heart ventricles, BAT, quadriceps, and diaphragm, while those for sex are heart (ventricles and atria), kidney (cortex and medulla), and BAT. The highest contributors to the differences in respiration at CIV as a result of age were heart (ventricles and atria), BAT, quadriceps, and diaphragm, while those for sex were BAT and heart (ventricles and atria). It is important to note that these rankings do not indicate directionality of a difference, only total magnitude. Additionally, these changes are represented as absolute magnitude of difference. Therefore, tissues with lower absolute oxygen consumption values would skew towards lower ranking, even though the relative change occurring within that tissue may be large.

Principle component (PC) analysis of each mitochondrial respiration dataset showed a primary and consistent separation across PC1 as a result of age (Figure 7E). A small percentage of the variance can be explained through PC2 as a result of sex, although this effect was not as consistent as that seen with age. Pearson correlation heat maps further highlighted age as a key factor relating samples across groups (Figure 7F). Male and female samples of the same age positively correlated with one another, while young and old samples regardless of sex showed a negative correlation. These combined data indicate that while sex does affect mitochondrial respiration, age is the dominant factor.

The sex-specific impact of age on mitochondrial function

To investigate the sex-dependent effects on mitochondrial function across age, we first calculated all of the relative change values for males and females from young to old. Summation of relative change values of all tissues within a specific mitochondrial parameter—respiration through CI, CII, or CIV—allowed us to view sex-specific systems-level effects across age (Figure 8A-B). Both the male and female systems showed a relative net increase in CI and decrease in CIV activity with age. Interestingly, with respect to CII activity, the male and female systems had divergent response. The net relative CII activity increased in males and decreased in females with age. These data suggest that: 1) independent of sex, CI and CIV are uniquely regulated, showing net opposite responses with age, and 2) CII activity is the most sex-affected respiratory component across age and tissues.

Figure 8.

Heat maps of the individual tissue-level data presented in Figure 8A-B further highlight the tissue-level similarities and differences across age, sex, and mitochondrial respiration type (i.e., through CI, CII, or CIV). For example, respiration through CI, CII, or CIV showed that the plantaris muscles (PL), diaphragm muscle (DI), quadriceps muscles (QD), jejunum (JE), and mesenteric white adipose tissue (MW) have the greatest positive relative change in response to age, while the eye (EY), gonadal white adipose tissue (GW), skin (SK), and stomach (ST) were among the tissues that showed the greatest negative relative change with age (Figure 8C).

To determine the tissues with the strongest sex-specific age effects, we organized male and female data by the largest (left) to smallest (right) difference of relative change. This multi-dimensional view highlights differences within a tissue across sex in response to age (Figure 8D). The largest sex-specific effects of age on respiration via CI were found in the stomach (ST), jenujum (JE), distal colon (DC), quadriceps muscles (QD), and pancreas (PN); in contrast, the kidney cortex (KC), cecum (CC), cerebellum (CE), soleus muscle (SL), and diaphragm muscle (DI) displayed the most similarity of change (Figure 8D, top panel). The largest sex-specific effects of age on respiration through CII were found in the stomach (ST), duodenum (DU), distal colon (DC), pancreas (PN), and jejunum (JE); in contrast, the liver (LV), skin (SK), hamstring muscle (HS), kidney medulla (KM), and tongue (TN) displayed the smallest sex-differences (Figure 8D, middle panel). Finally, the most sex-specific effects of age on respiration through CIV were found in mesenteric white adipose tissue (MW), jejunum (JE), duodenum (DU), distal colon (DC), and skin (SK); in contrast, the tongue (TN), heart ventricles (HV), gastrocnemius muscle (GS), distal colon (DI), and brain cortex (CO) displayed the greatest similarity of change with age (Figure 8D, bottom panel).

All sexually divergent tissues are marked specifically within each graph per mitochondrial respiration-type. Respiration via CI showed 5 tissues with sexual divergence in response to aging, and these were stomach (ST), distal colon (DC), pancreas (PN), lung (LN), and kidney medulla (KM) (Figure 8D, top panel). Respiration through CII showed 5 tissues with sexual divergence in response to aging, and these were duodenum (DU), pancreas (PN), lung (LN), heart atria (HA), and tongue (TN) (Figure 8D, second panel). Respiration at CIV also showed 5 tissues with sexual divergence in response to aging, and these were jejunum (JE), distal colon (DC), hamstring (HS), heart atria (HA), and kidney medulla (KM) (Figure 8D, bottom). Together, these data show the relative magnitude of change across sex and age, and highlight the tissue- and sex-specific regulation of mitochondrial respiration with aging.

We further grouped the above data by directionality to visualize male and female tissue responses to aging (Figure 8E). From this, we can quickly visualize tissues displaying a shared direction of change in respiration via CI, CII, or CIV across sex and age. Additionally, some tissues responded in an age-and-sex-specific manner (i.e., display divergence in Figure 8D). This data clustering summarizes the shared (age-specific) or opposite (sex-specific) directionality of mitochondrial function across tissues in respond to age.

It should be noted, however, that the mitochondrial activity of a particular tissue could increase or decrease in a sex-specific way even though the directionality is the same. For example, respiration through CII for the stomach (ST) showed male and female values decreasing with age, yet the male had a larger relative decrease when compared to female (Figure 8D, middle panel). Another consideration should be given to the divergent group with varying degrees of divergence. For example, CI-associated respiration of the stomach (Figure 8D, top panel) showed a strong divergence given that male had a strong decline in function while the female increased over time. On the other hand, the CI-associated respiration of the kidney medulla showed a weak divergence; although the directionality of change was opposite, the magnitude was small.

Lastly, we grouped all the tissues with a consistent shared (across sex) relative increase or decrease in CI-, CII-, and CIV-associated respiration in response to age (Figure 8F). The tissues with a shared (across sex) increased response to aging were plantaris (PL), quadriceps (QD), diaphragm (DI), soleus (SL), heart ventricle (HV), ileum (IL), cecum (CC), proximal colon (PC), and mesenteric fat (MW). Interestingly, PL, QD, DI, and SL are all skeletal muscles and IL, CC, and PC are all part of the digestive tract. Also of note, the gastrocnemius (GT), hamstring (HS), and tongue (TN) muscles were not uniformly increasing with age. In fact, these muscles showed reduced mitochondrial respiration with age, in striking contrast to other skeletal muscle types assayed. The tissues with a shared (across sex) decrease in response to aging were brown fat (BT), gonadal fat (GW), inguinal fat (IW), brain cortex (CO), hippocampus (HI), hypothalamus (HY), cerebellum (CE), liver (LV), kidney cortex (KC), eye (EY), skin (SK), and spleen (SP). Different fat depots and many of the brain regions showed a consistent decline in respiratory function with age, regardless of sex. The kidney cortex (KC) oddly showed a respiratory effect distinct from the medulla, which was almost completely resistant to age-related respiratory changes. The eyes (EY) not only showed a consistent decline in respiration with age, but also had one of the largest combined relative effects. Together, the shared and divergent analyses presented (Figure 8E and F) illuminate the modulation of mitochondrial function across age.

Discussion

Here, we provide the most comprehensive catalog of sex- and age-associated mitochondrial respiration signatures across diverse tissues and organs system. One important aspect of our study to emphasize is the experimental uniformity. All the tissue samples were collected from the same set of mice, and processed and analyzed in a uniform and standardized manner on the same respirometry platform in a single laboratory. Combined with the ten biological replicates, our experimental setup and workflow markedly enhanced the rigor and reduced the variations in data quality. From the mitochondrial respiration compendium, we can clearly see that mitochondrial function varies widely across tissues and that aging has a much larger effect than sex on mitochondrial activity across tissues.

The transcriptomic dataset generated by the Tabula Muris Consortium clearly indicates organ-specific temporal signatures across the lifespan (7–9), and this is assumed to track closely with tissue function even though no equivalent systematic interrogation of tissue function across the lifespan has been performed. In this regard, our mitochondrial activity atlas aligned and corroborated with the conclusion from large-scale transcriptomic analyses that the aging process and trajectory can vary substantially across tissues (7, 8). This reinforces the notion that while an organism ages, the aging-associated decline in organ function is not uniform across tissues and organ systems. Given that genetics and environmental factors all contribute to tissue health and their functional breakdown with age, this further amplifies the heterogenous process of aging across a population.

Although sex is a significant biological variable influencing many aspects of physiology and pathophysiology (29, 30), we were surprised that sex has a relatively modest impact on respiratory capacity in mitochondria across tissues and age. The modest impact of sex on mitochondrial respiration across tissues is also consistent with the relatively minor impact of sex on tissue transcriptome across the lifespan (7, 11). The caveat, however, is that our study is based on a single strain of mice with a uniform genetic background. The contribution of biological sex may be more pronounced in an outbred population where sex hormones interact with other genetic and environmental determinants to influence mitochondrial activity across tissues and age.

Aging is a dynamic process that unfolds across the lifespan, and this involves changes at the molecular, cellular, biochemical, and tissue level. Transcriptomics analyses have consistently shown that aging induces a reduction in mitochondrial OXPHOS genes (7, 8), and this likely contributes to the reduced mitochondrial activity we observed across many tissues in old mice. Changes in cell function, however, are often counter measured by homeostatic responses. While many tissues showed a significant reduction in mitochondrial function, some tissues also showed a concomitant increase in mitochondrial activity with age. With the assumption that organ function declines with age, our data suggests a potential organ-selective compensatory response to declining efficiency in cellular function which entails a greater demand on mitochondrial ATP production. This hypothesis awaits future experimental confirmation.

Every respirometry method has its advantages and limitations (27). Unlike the mitochondrial respiration analysis in intact cells or isolated mitochondria, our respirometry analysis was performed in mitochondria-enriched lysate derived from frozen tissues (28). In an intact mitochondrion, the rate of mitochondrial respiration is regulated by both substrate supply across the mitochondrial membrane and the rate of NADH and FADH production via the TCA cycle and fatty acid β-oxidation. In the absence of an intact mitochondrial membrane, our assays measure maximal mitochondrial respiration through CI, CII, or CIV. Thus, our mitochondrial activity reflects maximal respiration across tissues. In an in vivo milieu, mitochondrial activity is likely regulated and may not reach its maximal capacity as reflected in our assay. However, the advantage of the frozen-tissue method makes it feasible to carry out large-scale respirometry analyses, which outweigh its limitations.

In summary, our mitochondrial functional signatures demonstrate similar and divergent responses to aging across tissues. Whether changes in tissue mitochondrial respiration reflect the cause, consequence, or both, of aging remains to be determined. We anticipate that the integration of this knowledge with other large-scale omics data will help uncover potential genetic, epigenetic, and biochemical mechanisms linking mitochondrial health and organismal aging.

Materials and methods

Mouse model

All wild-type C57BL/6J male and female mice were purchased from the Jackson Laboratory and fed a standard chow (Envigo; 2018SX). In total, the young group was comprised of 10 male and 10 female 8-week-old mice, and the old group was comprised of 10 male and 10 female 80-week-old mice. Mice were housed in polycarbonate cages on a 12h:12h light-dark photocycle with ad libitum access to water and food. All mice were fasted for 2 h prior to euthanasia and dissection. Tissues were collected, snap-frozen in liquid nitrogen, and kept at −80°C until analysis. All mouse protocols were approved by the Institutional Animal Care and Use Committee of the Johns Hopkins University School of Medicine (animal protocol # MO22M367). All animal experiments were conducted in accordance with the National Institute of Health guidelines and followed the standards established by the Animal Welfare Acts.

Comprehensive multi-organ dissection

Each mouse was dissected cleanly, swiftly, and in a concerted manner by three people. Each mouse dissection took approximately 8-10 min, with 33 tissues collected per dissection. After euthanasia, blood was collected via decapitation. The head was then immediately given to dissector one to collect brain regions (hypothalamus, cerebellum, hippocampus, and cortex), eyes, and tongue. Simultaneously, the visceral cavity was opened by dissector two. Inguinal white adipose tissue (iWAT) was collected immediately after opening the abdominal skin. After that the abdominal muscle was cut and gonadal white adipose tissue (gWAT) and testes or fallopian tubes were collected. Following this, the visceral organs were partitioned in two groups. Dissector two promptly dissected liver, stomach, kidneys (further separated the cortex and medulla), spleen, diaphragm, heart (further divided into atria and ventricles), lungs, brown adipose tissue (BAT), and skin (cleared of hair using Topical Nair lotion, cleaned, then collected). At the same time, dissector three was collecting the pancreas, mesenteric white adipose tissue (mesWAT), small intestine (further split into duodenum, jejunum, and ileum), cecum, and large intestine (further split into proximal and distal colon). As soon as tissue collection in the head was finished, the mouse carcass was cut transversely at the lumbar spine and handed to dissector one for muscle dissection. Dissector one then rapidly and precisely anatomized (bilaterally) the quadriceps (entire complex – rectus femoris, vastus lateralis, vastus intermedius, and vastus medialis), hamstrings (biceps femoris), gastrocnemius, plantaris, and soleus muscles. All tissues were washed with sterile 1X PBS to remove residual blood prior to snap freeze. Additionally, the stomach, small intestine, and large intestine were cleared and cleaned of debris with PBS prior to freezing. The Cecum, however, was kept whole containing all fecal/food matter and microorganisms present. All dissected tissues were snap frozen in liquid nitrogen, and stored at −80°C for later analysis.

Respirometry of frozen tissue samples

Respirometry was conducted on frozen tissue samples to assay for mitochondrial activity as described previously (28), using a Seahorse XFe96 Analyzer. Samples were thawed in 1X MAS buffer (70 mM sucrose, 220 mM mannitol, 5 mM KH2PO4, 5 mM MgCl2, 1 mM EGTA, 2 mM HEPES pH 7.4), finely minced with scissors, then homogenized with a glass Dounce homogenizer on ice. The entire sample, bilateral for skeletal muscles, was homogenized and used for respirometry with the exception of liver in which only a piece was homogenized. This was done to provide a true tissue average and avoid regional differences in mitochondria. The resulting homogenate was spun at 1000 g for 10 min at 4°C. The supernatant was collected and immediately used for protein quantification by BCA assay (Thermo Scientific, 23225). Each well of the Seahorse microplate was loaded with the designated amount (in µg) of homogenate protein (Table 1). Each biological replicate is comprised of three technical replicates. Samples from all tissues were treated separately with NADH (1 mM) as a complex I substrate or Succinate (a complex II substrate, 5 mM) in the presence of rotenone (a complex I inhibitor, 2 µM), then with the inhibitors rotenone (2 µM) and Antimycin A (4 µM), followed by TMPD (N,N,N’,N’-Tetramethyl-p-Phenylenediamine; also known as Wurster’s Reagent, 0.45 mM) and Ascorbate (Vitamin C, 1 mM) to activate complex IV, and finally treated with Azide (40 mM) to assess non-mitochondrial respiration.

Table 1.

Quantification of mitochondrial content

Mitochondrial content of homogenates used for respirometry was quantified with MitoTracker Deep Red FM (MTDR, Invitrogen, M22426) staining as described previously (28). Briefly, lysates were incubated with MTDR (1 µM) for 10 min at 37°C, then centrifuged at 2000 g for 5 min at 4°C. The supernatant was carefully removed and replaced with 1X MAS solution and fluorescence was read with excitation and emission wavelengths of 625 nm and 670 nm, respectively. All measurements were read with a BioTek Synergy HTX Multimode Plate Reader (Agilent). All samples from a single tissue-type were measured on the same plate, in duplicate, with an equal gain setting of 90 to ensure comparability across samples. To minimize non-specific background signal contribution, control wells were loaded with MTDR + 1X MAS and subtracted from all sample values.

Acknowledgements

This work was supported by the National Institutes of Health (DK084171 to GWW). D.C.S. was supported by an NIH T32 training grant (HL007534).

Author contributions

DCS and GWW contributed to the experimental design; DCS, MS, FC, and GWW performed the experiments; DCS and GWW analyzed and interpreted the data; DCS and GWW drafted the paper with inputs and edits from all authors.

Competing interests

We declare that none of the authors has a conflict of interest.

Data availability

All the data presented are included in the main text or as supplementary tables and figures.

Abbreviations

• BAT: Brown adipose tissue

• CI: Mitochondrial complex I

• CII: Mitochondrial complex II

• CIII: Mitochondrial complex IIII

• CIV: Mitochondrial complex IV

• FADH: flavin adenine dinucleotide

• gWAT: Gonadal white adipose tissue

• iWAT: Inguinal white adipose tissue NADH nicotinamide adenine dinucleotide OCR Oxygen consumption rate

• Rot: Rotenone

• TMPD: N,N,N’,N’-Tetramethyl-p-Phenylenediamine

Supplemental figure legends

Figure 1—figure supplement 1

Figure 1—figure supplement 2

Figure 1—figure supplement 3

Figure 1—figure supplement 4

Figure 1—figure supplement 5

Figure 1—figure supplement 6

Figure 1—figure supplement 7

Figure 1—figure supplement 8

Figure 1—figure supplement 9

Figure 1—figure supplement 10

Figure 1—figure supplement 11

Figure 1—figure supplement 12

Figure 1—figure supplement 13

Figure 1—figure supplement 14

Figure 1—figure supplement 15

Figure 1—figure supplement 16

Figure 1—figure supplement 17

Figure 1—figure supplement 18

Figure 1—figure supplement 19

Figure 1—figure supplement 20

Figure 1—figure supplement 21

Figure 1—figure supplement 22

Figure 1—figure supplement 23

Figure 1—figure supplement 24

Figure 1—figure supplement 25

Figure 1—figure supplement 26

Figure 1—figure supplement 27

Figure 1—figure supplement 28

Figure 1—figure supplement 29

Figure 1—figure supplement 30

Figure 1—figure supplement 31

Figure 1—figure supplement 32

Figure 1—figure supplement 33

Figure 2-figure supplement 1

Figure 2-figure supplement 2

Figure 2-figure supplement 3

Figure 2-figure supplement 4

Figure 2-figure supplement 5

Figure 2-figure supplement 6