Two Conserved Histone Demethylases Regulate Mitochondrial Stress-Induced Longevity.

Across eukaryotic species, mild mitochondrial stress can have beneficial effects on the lifespan of organisms. Mitochondrial dysfunction activates an unfolded protein response (UPR(mt)), a stress signaling mechanism designed to ensure mitochondrial homeostasis. Perturbation of mitochondria during larval development in C. elegans not only delays aging but also maintains UPR(mt) signaling, suggesting an epigenetic mechanism that modulates both longevity and mitochondrial proteostasis throughout life. We identify the conserved histone lysine demethylases jmjd-1.2/PHF8 and jmjd-3.1/JMJD3 as positive regulators of lifespan in response to mitochondrial dysfunction across species. Reduction of function of the demethylases potently suppresses longevity and UPR(mt) induction, while gain of function is sufficient to extend lifespan in a UPR(mt)-dependent manner. A systems genetics approach in the BXD mouse reference population further indicates conserved roles of the mammalian orthologs in longevity and UPR(mt) signaling. These findings illustrate an evolutionary conserved epigenetic mechanism that determines the rate of aging downstream of mitochondrial perturbations.

The Convergence of Systems and Reductionist Approaches in Complex Trait Analysis.

Research into the genetic and environmental factors behind complex trait variation has traditionally been segregated into distinct scientific camps. The reductionist approach aims to decrypt phenotypic variability bit by bit, founded on the underlying hypothesis that genome-to-phenome relations are largely constructed from the additive effects of their molecular players. In contrast, the systems approach aims to examine large-scale interactions of many components simultaneously, on the premise that interactions in gene networks can be both linear and non-linear. Both approaches are complementary, and they are becoming increasingly intertwined due to developments in gene editing tools, omics technologies, and population resources. Together, these strategies are beginning to drive the next era in complex trait research, paving the way to improve agriculture and toward more personalized medicine.

Multilayered genetic and omics dissection of mitochondrial activity in a mouse reference population.

The manner by which genotype and environment affect complex phenotypes is one of the fundamental questions in biology. In this study, we quantified the transcriptome--a subset of the metabolome--and, using targeted proteomics, quantified a subset of the liver proteome from 40 strains of the BXD mouse genetic reference population on two diverse diets. We discovered dozens of transcript, protein, and metabolite QTLs, several of which linked to metabolic phenotypes. Most prominently, Dhtkd1 was identified as a primary regulator of 2-aminoadipate, explaining variance in fasted glucose and diabetes status in both mice and humans. These integrated molecular profiles also allowed further characterization of complex pathways, particularly the mitochondrial unfolded protein response (UPR(mt)). UPR(mt) shows strikingly variant responses at the transcript and protein level that are remarkably conserved among C. elegans, mice, and humans. Overall, these examples demonstrate the value of an integrated multilayered omics approach to characterize complex metabolic phenotypes.

Systems genetics of metabolism: the use of the BXD murine reference panel for multiscalar integration of traits.

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.

NCoR1 is a conserved physiological modulator of muscle mass and oxidative function.

Transcriptional coregulators control the activity of many transcription factors and are thought to have wide-ranging effects on gene expression patterns. We show here that muscle-specific loss of nuclear receptor corepressor 1 (NCoR1) in mice leads to enhanced exercise endurance due to an increase of both muscle mass and of mitochondrial number and activity. The activation of selected transcription factors that control muscle function, such as MEF2, PPARß/δ, and ERRs, underpins these phenotypic alterations. NCoR1 levels are decreased in conditions that require fat oxidation, resetting transcriptional programs to boost oxidative metabolism. Knockdown of gei-8, the sole C. elegans NCoR homolog, also robustly increased muscle mitochondria and respiration, suggesting conservation of NCoR1 function. Collectively, our data suggest that NCoR1 plays an adaptive role in muscle physiology and that interference with NCoR1 action could be used to improve muscle function.

PCfun: a hybrid computational framework for systematic characterization of protein complex function.

In molecular biology, it is a general assumption that the ensemble of expressed molecules, their activities and interactions determine biological function, cellular states and phenotypes. Stable protein complexes-or macromolecular machines-are, in turn, the key functional entities mediating and modulating most biological processes. Although identifying protein complexes and their subunit composition can now be done inexpensively and at scale, determining their function remains challenging and labor intensive. This study describes Protein Complex Function predictor (PCfun), the first computational framework for the systematic annotation of protein complex functions using Gene Ontology (GO) terms. PCfun is built upon a word embedding using natural language processing techniques based on 1 million open access PubMed Central articles. Specifically, PCfun leverages two approaches for accurately identifying protein complex function, including: (i) an unsupervised approach that obtains the nearest neighbor (NN) GO term word vectors for a protein complex query vector and (ii) a supervised approach using Random Forest (RF) models trained specifically for recovering the GO terms of protein complex queries described in the CORUM protein complex database. PCfun consolidates both approaches by performing a hypergeometric statistical test to enrich the top NN GO terms within the child terms of the GO terms predicted by the RF models. The documentation and implementation of the PCfun package are available at https://github.com/sharmavaruns/PCfun. We anticipate that PCfun will serve as a useful tool and novel paradigm for the large-scale characterization of protein complex function.

Diagnostics and correction of batch effects in large-scale proteomic studies: a tutorial.

Advancements in mass spectrometry-based proteomics have enabled experiments encompassing hundreds of samples. While these large sample sets deliver much-needed statistical power, handling them introduces technical variability known as batch effects. Here, we present a step-by-step protocol for the assessment, normalization, and batch correction of proteomic data. We review established methodologies from related fields and describe solutions specific to proteomic challenges, such as ion intensity drift and missing values in quantitative feature matrices. Finally, we compile a set of techniques that enable control of batch effect adjustment quality. We provide an R package, "proBatch", containing functions required for each step of the protocol. We demonstrate the utility of this methodology on five proteomic datasets each encompassing hundreds of samples and consisting of multiple experimental designs. In conclusion, we provide guidelines and tools to make the extraction of true biological signal from large proteomic studies more robust and transparent, ultimately facilitating reliable and reproducible research in clinical proteomics and systems biology.

Quantifying and Localizing the Mitochondrial Proteome Across Five Tissues in A Mouse Population.

We have used SWATH mass spectrometry to quantify 3648 proteins across 76 proteomes collected from genetically diverse BXD mouse strains in two fractions (mitochondria and total cell) from five tissues: liver, quadriceps, heart, brain, and brown adipose (BAT). Across tissues, expression covariation between genes' proteins and transcripts-measured in the same individuals-broadly aligned. Covariation was however far stronger in certain subsets than others: only 8% of transcripts in the lowest expression and variance quintile covaried with their protein, in contrast to 65% of transcripts in the highest quintiles. Key functional differences among the 3648 genes were also observed across tissues, with electron transport chain (ETC) genes particularly investigated. ETC complex proteins covary and form strong gene networks according to tissue, but their equivalent transcripts do not. Certain physiological consequences, such as the depletion of ATP synthase in BAT, are thus obscured in transcript data. Lastly, we compared the quantitative proteomic measurements between the total cell and mitochondrial fractions for the five tissues. The resulting enrichment score highlighted several hundred proteins which were strongly enriched in mitochondria, which included several dozen proteins were not reported in literature to be mitochondrially localized. Four of these candidates were selected for biochemical validation, where we found MTAP, SOAT2, and IMPDH2 to be localized inside the mitochondria, whereas ABCC6 was in the mitochondria-associated membrane. These findings demonstrate the synergies of a multi-omics approach to study complex metabolic processes, and this provides a resource for further discovery and analysis of proteoforms, modified proteins, and protein localization.

Mild inborn errors of metabolism in commonly used inbred mouse strains.

Inbred mouse strains are a cornerstone of translational research but paradoxically many strains carry mild inborn errors of metabolism. For example, α-aminoadipic acidemia and branched-chain ketoacid dehydrogenase deficiency are known in C57BL/6J mice. Using RNA sequencing, we now reveal the causal variants in Dhtkd1 and Bckdhb, and the molecular mechanism underlying these metabolic defects. C57BL/6J mice have decreased Dhtkd1 mRNA expression due to a solitary long terminal repeat (LTR) in intron 4 of Dhtkd1. This LTR harbors an alternate splice donor site leading to a partial splicing defect and as a consequence decreased total and functional Dhtkd1 mRNA, decreased DHTKD1 protein and α-aminoadipic acidemia. Similarly, C57BL/6J mice have decreased Bckdhb mRNA expression due to an LTR retrotransposon in intron 1 of Bckdhb. This transposable element encodes an alternative exon 1 causing aberrant splicing, decreased total and functional Bckdhb mRNA and decreased BCKDHB protein. Using a targeted metabolomics screen, we also reveal elevated plasma C5-carnitine in 129 substrains. This biochemical phenotype resembles isovaleric acidemia and is caused by an exonic splice mutation in Ivd leading to partial skipping of exon 10 and IVD protein deficiency. In summary, this study identifies three causal variants underlying mild inborn errors of metabolism in commonly used inbred mouse strains.

Genetic cartography of longevity in humans and mice: Current landscape and horizons.

Aging is a complex and highly variable process. Heritability of longevity among humans and other species is low, and this finding has given rise to the idea that it may be futile to search for DNA variants that modulate aging. We argue that the problem in mapping longevity genes is mainly one of low power and the genetic and environmental complexity of aging. In this review we highlight progress made in mapping genes and molecular networks associated with longevity, paying special attention to work in mice and humans. We summarize 40 years of linkage studies using murine cohorts and 15 years of studies in human populations that have exploited candidate gene and genome-wide association methods. A small but growing number of gene variants contribute to known longevity mechanisms, but a much larger set have unknown functions. We outline these and other challenges and suggest some possible solutions, including more intense collaboration between research communities that use model organisms and human cohorts. Once hundreds of gene variants have been linked to differences in longevity in mammals, it will become feasible to systematically explore gene-by-environmental interactions, dissect mechanisms with more assurance, and evaluate the roles of epistasis and epigenetics in aging. A deeper understanding of complex networks-genetic, cellular, physiological, and social-should position us well to improve healthspan.

An evolutionarily conserved role for the aryl hydrocarbon receptor in the regulation of movement.

The BXD genetic reference population is a recombinant inbred panel descended from crosses between the C57BL/6 (B6) and DBA/2 (D2) strains of mice, which segregate for about 5 million sequence variants. Recently, some of these variants have been established with effects on general metabolic phenotypes such as glucose response and bone strength. Here we phenotype 43 BXD strains and observe they have large variation (-5-fold) in their spontaneous activity during waking hours. QTL analyses indicate that -40% of this variance is attributable to a narrow locus containing the aryl hydrocarbon receptor (Ahr), a basic helix-loop-helix transcription factor with well-established roles in development and xenobiotic metabolism. Strains with the D2 allele of Ahr have reduced gene expression compared to those with the B6 allele, and have significantly higher spontaneous activity. This effect was also observed in B6 mice with a congenic D2 Ahr interval, and in B6 mice with a humanized AHR allele which, like the D2 allele, is expressed much less and has less enzymatic activity than the B6 allele. Ahr is highly conserved in invertebrates, and strikingly inhibition of its orthologs in D. melanogaster and C. elegans (spineless and ahr-1) leads to marked increases in basal activity. In mammals, Ahr has numerous ligands, but most are either non-selective (e.g. resveratrol) or highly toxic (e.g., 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)). Thus, we chose to examine a major environmental influence--long term feeding with high fat diet (HFD)--to see if the effects of Ahr are dependent on major metabolic differences. Interestingly, while HFD robustly halved movement across all strains, the QTL position and effects of Ahr remained unchanged, indicating that the effects are independent. The highly consistent effects of Ahr on movement indicate that changes in its constitutive activity have a role on spontaneous movement and may influence human behavior.

Gene network based analysis identifies a coexpression module involved in regulating plasma lipids with high-fat diet response.

Plasma lipids are modulated by gene variants and many environmental factors, including diet-associated weight gain. However, understanding how these factors jointly interact to influence molecular networks that regulate plasma lipid levels is limited. Here, we took advantage of the BXD recombinant inbred family of mice to query weight gain as an environmental stressor on plasma lipids. Coexpression networks were examined in both nonobese and obese livers, and a network was identified that specifically responded to the obesogenic diet. This obesity-associated module was significantly associated with plasma lipid levels and enriched with genes known to have functions related to inflammation and lipid homeostasis. We identified key drivers of the module, including Cidec, Cidea, Pparg, Cd36, and Apoa4. The Pparg emerged as a potential master regulator of the module as it can directly target 19 of the top 30 hub genes. Importantly, activation of this module is causally linked to lipid metabolism in humans, as illustrated by correlation analysis and inverse-variance weighed Mendelian randomization. Our findings provide novel insights into gene-by-environment interactions for plasma lipid metabolism that may ultimately contribute to new biomarkers, better diagnostics, and improved approaches to prevent or treat dyslipidemia in patients.

Multiomic profiling of the liver across diets and age in a diverse mouse population.

We profiled the liver transcriptome, proteome, and metabolome in 347 individuals from 58 isogenic strains of the BXD mouse population across age (7 to 24 months) and diet (low or high fat) to link molecular variations to metabolic traits. Several hundred genes are affected by diet and/or age at the transcript and protein levels. Orthologs of two aging-associated genes, St7 and Ctsd, were knocked down in C. elegans, reducing longevity in wild-type and mutant long-lived strains. The multiomics data were analyzed as segregating gene networks according to each independent variable, providing causal insight into dietary and aging effects. Candidates were cross-examined in an independent diversity outbred mouse liver dataset segregating for similar diets, with ∼80%-90% of diet-related candidate genes found in common across datasets. Together, we have developed a large multiomics resource for multivariate analysis of complex traits and demonstrate a methodology for moving from observational associations to causal connections.

GCN5 maintains muscle integrity by acetylating YY1 to promote dystrophin expression.

Protein lysine acetylation is a post-translational modification that regulates protein structure and function. It is targeted to proteins by lysine acetyltransferases (KATs) or removed by lysine deacetylases. This work identifies a role for the KAT enzyme general control of amino acid synthesis protein 5 (GCN5; KAT2A) in regulating muscle integrity by inhibiting DNA binding of the transcription factor/repressor Yin Yang 1 (YY1). Here we report that a muscle-specific mouse knockout of GCN5 (Gcn5skm-/-) reduces the expression of key structural muscle proteins, including dystrophin, resulting in myopathy. GCN5 was found to acetylate YY1 at two residues (K392 and K393), disrupting the interaction between the YY1 zinc finger region and DNA. These findings were supported by human data, including an observed negative correlation between YY1 gene expression and muscle fiber diameter. Collectively, GCN5 positively regulates muscle integrity through maintenance of structural protein expression via acetylation-dependent inhibition of YY1. This work implicates the role of protein acetylation in the regulation of muscle health and for consideration in the design of novel therapeutic strategies to support healthy muscle during myopathy or aging.

Integrative systems analysis identifies genetic and dietary modulators of bile acid homeostasis.

Bile acids (BAs) are complex and incompletely understood enterohepatic-derived hormones that control whole-body metabolism. Here, we profiled postprandial BAs in the liver, feces, and plasma of 360 chow- or high-fat-diet-fed BXD male mice and demonstrated that both genetics and diet strongly influence BA abundance, composition, and correlation with metabolic traits. Through an integrated systems approach, we mapped hundreds of quantitative trait loci that modulate BAs and identified both known and unknown regulators of BA homeostasis. In particular, we discovered carboxylesterase 1c (Ces1c) as a genetic determinant of plasma tauroursodeoxycholic acid (TUDCA), a BA species with established disease-preventing actions. The association between Ces1c and plasma TUDCA was validated using data from independent mouse cohorts and a Ces1c knockout mouse model. Collectively, our data are a unique resource to dissect the physiological importance of BAs as determinants of metabolic traits, as underscored by the identification of CES1C as a master regulator of plasma TUDCA levels.

A platform for experimental precision medicine: The extended BXD mouse family.

The challenge of precision medicine is to model complex interactions among DNA variants, phenotypes, development, environments, and treatments. We address this challenge by expanding the BXD family of mice to 140 fully isogenic strains, creating a uniquely powerful model for precision medicine. This family segregates for 6 million common DNA variants-a level that exceeds many human populations. Because each member can be replicated, heritable traits can be mapped with high power and precision. Current BXD phenomes are unsurpassed in coverage and include much omics data and thousands of quantitative traits. BXDs can be extended by a single-generation cross to as many as 19,460 isogenic F1 progeny, and this extended BXD family is an effective platform for testing causal modeling and for predictive validation. BXDs are a unique core resource for the field of experimental precision medicine.

Molecular pathways behind acquired obesity: Adipose tissue and skeletal muscle multiomics in monozygotic twin pairs discordant for BMI.

Tissue-specific mechanisms prompting obesity-related development complications in humans remain unclear. We apply multiomics analyses of subcutaneous adipose tissue and skeletal muscle to examine the effects of acquired obesity among 49 BMI-discordant monozygotic twin pairs. Overall, adipose tissue appears to be more affected by excess body weight than skeletal muscle. In heavier co-twins, we observe a transcriptional pattern of downregulated mitochondrial pathways in both tissues and upregulated inflammatory pathways in adipose tissue. In adipose tissue, heavier co-twins exhibit lower creatine levels; in skeletal muscle, glycolysis- and redox stress-related protein and metabolite levels remain higher. Furthermore, metabolomics analyses in both tissues reveal that several proinflammatory lipids are higher and six of the same lipid derivatives are lower in acquired obesity. Finally, in adipose tissue, but not in skeletal muscle, mitochondrial downregulation and upregulated inflammation are associated with a fatty liver, insulin resistance, and dyslipidemia, suggesting that adipose tissue dominates in acquired obesity.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family.

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Body weight and high-fat diet are associated with epigenetic aging in female members of the BXD murine family.

DNA methylation (DNAm) is shaped by genetic and environmental factors and modulated by aging. Here, we examine interrelations between epigenetic aging, body weight (BW), and life span in 12 isogenic strains from the BXD family of mice that exhibit over twofold variation in longevity. Genome-wide DNAm was assayed in 70 liver specimens from predominantly female cases, 6-25 months old, that were maintained on normal chow or high-fat diet (HFD). We defined subsets of CpG regions associated with age, BW at young adulthood, and strain-by-diet-dependent life span. These age-associated differentially methylated CpG regions (age-DMRs) featured distinct genomic characteristics, with DNAm gains over time occurring in sites such as promoters and exons that have high CpG density and low average methylation. CpG regions associated with BW were enriched in introns, tended to have lower methylation in mice with higher BW, and were inversely correlated with gene expression (i.e., higher mRNA levels in mice with higher BW). CpG regions associated with life span were linked to genes involved in life span modulation, including the telomerase reverse transcriptase gene, Tert, which had both lower methylation and higher expression in long-lived strains. An epigenetic clock defined from age-DMRs revealed accelerated aging in mice belonging to strains with shorter life spans. Both higher BW and the HFD were associated with accelerated epigenetic aging. Our results highlight the age-accelerating effect of heavier BW. Furthermore, we demonstrate that the measure of epigenetic aging derived from age-DMRs can predict genotype and diet-induced differences in life span among female BXD members.