Systems genetics of metabolic health in the BXD mouse genetic reference population.

Susceptibility to metabolic syndrome (MetS) is dependent on genetics, environment, and gene-by-environment interactions, rendering the study of underlying mechanisms challenging. The majority of experiments in model organisms do not incorporate genetic variation and lack specific evaluation criteria for MetS. Here, we derived a continuous metric, the metabolic health score (MHS), based on standard clinical parameters and defined its molecular signatures in the liver and circulation. In human UK Biobank, the MHS associated with MetS status and was predictive of future disease incidence, even in individuals without MetS. Using quantitative trait locus analyses in mice, we found two MHS-associated genetic loci and replicated them in unrelated mouse populations. Through a prioritization scheme in mice and human genetic data, we identified TNKS and MCPH1 as candidates mediating differences in the MHS. Our findings provide insights into the molecular mechanisms sustaining metabolic health across species and uncover likely regulators. A record of this paper's transparent peer review process is included in the supplemental information.

Hippo-YAP/TAZ signalling coordinates adipose plasticity and energy balance by uncoupling leptin expression from fat mass.

Adipose tissues serve as an energy reservoir and endocrine organ, yet the mechanisms that coordinate these functions remain elusive. Here, we show that the transcriptional coregulators, YAP and TAZ, uncouple fat mass from leptin levels and regulate adipocyte plasticity to maintain metabolic homeostasis. Activating YAP/TAZ signalling in adipocytes by deletion of the upstream regulators Lats1 and Lats2 results in a profound reduction in fat mass by converting mature adipocytes into delipidated progenitor-like cells, but does not cause lipodystrophy-related metabolic dysfunction, due to a paradoxical increase in circulating leptin levels. Mechanistically, we demonstrate that YAP/TAZ-TEAD signalling upregulates leptin expression by directly binding to an upstream enhancer site of the leptin gene. We further show that YAP/TAZ activity is associated with, and functionally required for, leptin regulation during fasting and refeeding. These results suggest that adipocyte Hippo-YAP/TAZ signalling constitutes a nexus for coordinating adipose tissue lipid storage capacity and systemic energy balance through the regulation of adipocyte plasticity and leptin gene transcription.

Genetic and dietary modulators of the inflammatory response in the gastrointestinal tract of the BXD mouse genetic reference population.

Inflammatory gut disorders, including inflammatory bowel disease (IBD), can be impacted by dietary, environmental, and genetic factors. While the incidence of IBD is increasing worldwide, we still lack a complete understanding of the gene-by-environment interactions underlying inflammation and IBD. Here, we profiled the colon transcriptome of 52 BXD mouse strains fed with a chow or high-fat diet (HFD) and identified a subset of BXD strains that exhibit an IBD-like transcriptome signature on HFD, indicating that an interplay of genetics and diet can significantly affect intestinal inflammation. Using gene co-expression analyses, we identified modules that are enriched for IBD-dysregulated genes and found that these IBD-related modules share cis-regulatory elements that are responsive to the STAT2, SMAD3, and REL transcription factors. We used module quantitative trait locus analyses to identify genetic loci associated with the expression of these modules. Through a prioritization scheme involving systems genetics in the mouse and integration with external human datasets, we identified Muc4 and Epha6 as the top candidates mediating differences in HFD-driven intestinal inflammation. This work provides insights into the contribution of genetics and diet to IBD risk and identifies two candidate genes, MUC4 and EPHA6, that may mediate IBD susceptibility in humans.

The genetic background shapes the susceptibility to mitochondrial dysfunction and NASH progression.

Non-alcoholic steatohepatitis (NASH) is a global health concern without treatment. The challenge in finding effective therapies is due to the lack of good mouse models and the complexity of the disease, characterized by gene-environment interactions. We tested the susceptibility of seven mouse strains to develop NASH. The severity of the clinical phenotypes observed varied widely across strains. PWK/PhJ mice were the most prone to develop hepatic inflammation and the only strain to progress to NASH with extensive fibrosis, while CAST/EiJ mice were completely resistant. Levels of mitochondrial transcripts and proteins as well as mitochondrial function were robustly reduced specifically in the liver of PWK/PhJ mice, suggesting a central role of mitochondrial dysfunction in NASH progression. Importantly, the NASH gene expression profile of PWK/PhJ mice had the highest overlap with the human NASH signature. Our study exposes the limitations of using a single mouse genetic background in metabolic studies and describes a novel NASH mouse model with features of the human NASH.

COX7A2L genetic variants determine cardiorespiratory fitness in mice and human.

Mitochondrial respiratory complexes form superassembled structures called supercomplexes. COX7A2L is a supercomplex-specific assembly factor in mammals, although its implication for supercomplex formation and cellular metabolism remains controversial. Here we identify a role for COX7A2L for mitochondrial supercomplex formation in humans. By using human cis-expression quantitative trait loci data, we highlight genetic variants in the COX7A2L gene that affect its skeletal muscle expression specifically. The most significant cis-expression quantitative trait locus is a 10-bp insertion in the COX7A2L 3' untranslated region that increases messenger RNA stability and expression. Human myotubes harboring this insertion have more supercomplexes and increased respiration. Notably, increased COX7A2L expression in the muscle is associated with lower body fat and improved cardiorespiratory fitness in humans. Accordingly, specific reconstitution of Cox7a2l expression in C57BL/6J mice leads to higher maximal oxygen consumption, increased lean mass and increased energy expenditure. Furthermore, Cox7a2l expression in mice is induced specifically in the muscle upon exercise. These findings elucidate the genetic basis of mitochondrial supercomplex formation and function in humans and show that COX7A2L plays an important role in cardiorespiratory fitness, which could have broad therapeutic implications in reducing cardiovascular mortality.

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.

Sex- and age-dependent genetics of longevity in a heterogeneous mouse population.

DNA variants that modulate life span provide insight into determinants of health, disease, and aging. Through analyses in the UM-HET3 mice of the Interventions Testing Program (ITP), we detected a sex-independent quantitative trait locus (QTL) on chromosome 12 and identified sex-specific QTLs, some of which we detected only in older mice. Similar relations between life history and longevity were uncovered in mice and humans, underscoring the importance of early access to nutrients and early growth. We identified common age- and sex-specific genetic effects on gene expression that we integrated with model organism and human data to create a hypothesis-building interactive resource of prioritized longevity and body weight genes. Finally, we validated Hipk1, Ddost, Hspg2, Fgd6, and Pdk1 as conserved longevity genes using Caenorhabditis elegans life-span experiments.

Multi-omics analysis identifies essential regulators of mitochondrial stress response in two wild-type C. elegans strains.

The mitochondrial unfolded protein response (UPRmt) is a promising pharmacological target for aging and age-related diseases. However, the integrative analysis of the impact of UPRmt activation on different signaling layers in animals with different genetic backgrounds is lacking. Here, we applied systems approaches to investigate the effect of UPRmt induced by doxycycline (Dox) on transcriptome, proteome, and lipidome in two genetically divergent worm strains, named N2 and CB4856. From the integrated omics datasets, we found that Dox prolongs lifespan of both worm strains through shared and strain-specific mechanisms. Specifically, Dox strongly impacts mitochondria, upregulates defense response, and lipid metabolism, while decreasing triglycerides. We further validated that lipid genes acs-2/20 and fat-7/6 were required for Dox-induced UPRmt and longevity in N2 and CB4856 worms, respectively. Our data have translational value as they indicate that the beneficial effects of Dox-induced UPRmt on lifespan are consistent across different genetic backgrounds through different regulators.

Sphingolipids accumulate in aged muscle, and their reduction counteracts sarcopenia.

Age-related muscle dysfunction and sarcopenia are major causes of physical incapacitation in older adults and currently lack viable treatment strategies. Here we find that sphingolipids accumulate in mouse skeletal muscle upon aging and that both genetic and pharmacological inhibition of sphingolipid synthesis prevent age-related decline in muscle mass while enhancing strength and exercise capacity. Inhibition of sphingolipid synthesis confers increased myogenic potential and promotes protein synthesis. Within the sphingolipid pathway, we show that accumulation of dihydroceramides is the culprit disturbing myofibrillar homeostasis. The relevance of sphingolipid pathways in human aging is demonstrated in two cohorts, the UK Biobank and Helsinki Birth Cohort Study in which gene expression-reducing variants of SPTLC1 and DEGS1 are associated with improved and reduced fitness of older individuals, respectively. These findings identify sphingolipid synthesis inhibition as an attractive therapeutic strategy for age-related sarcopenia and co-occurring pathologies.

The Gene-Regulatory Footprint of Aging Highlights Conserved Central Regulators.

Many genes and pathways have been linked to aging, yet our understanding of underlying molecular mechanisms is still lacking. Here, we measure changes in the transcriptome, histone modifications, and DNA methylome in three metabolic tissues of adult and aged mice. Transcriptome and methylome changes dominate the liver aging footprint, whereas heart and muscle globally increase chromatin accessibility, especially in aging pathways. In mouse and human data from multiple tissues and regulatory layers, age-related transcription factor expression changes and binding site enrichment converge on putative aging modulators, including ZIC1, CXXC1, HMGA1, MECP2, SREBF1, SREBF2, ETS2, ZBTB7A, and ZNF518B. Using Mendelian randomization, we establish possible epidemiological links between expression of some of these transcription factors or their targets, including CXXC1, ZNF518B, and BBC3, and longevity. We conclude that conserved modulators are at the core of the molecular footprint of aging, and variation in tissue-specific expression of some may affect human longevity.

Bile Acids Signal via TGR5 to Activate Intestinal Stem Cells and Epithelial Regeneration.

BACKGROUND & AIMS: Renewal and patterning of the intestinal epithelium is coordinated by intestinal stem cells (ISCs); dietary and metabolic factors provide signals to the niche that control ISC activity. Bile acids (BAs), metabolites in the gut, signal nutrient availability by activating the G protein-coupled bile acid receptor 1 (GPBAR1, also called TGR5). TGR5 is expressed in the intestinal epithelium, but it is not clear how its activation affects ISCs and regeneration of the intestinal epithelium. We studied the role of BAs and TGR5 in intestinal renewal, and regulation of ISC function in mice and intestinal organoids. METHODS: We derived intestinal organoids from wild-type mice and Tgr5-/- mice, incubated them with BAs or the TGR5 agonist INT-777, and monitored ISC function by morphologic analyses and colony-forming assays. We disrupted Tgr5 specifically in Lgr5-positive ISCs in mice (Tgr5ISC-/- mice) and analyzed ISC number, proliferation, and differentiation by flow cytometry, immunofluorescence, and organoid assays. Tgr5ISC-/- mice were given cholecystokinin; we measured the effects of BA release into the intestinal lumen and on cell renewal. We induced colitis in Tgr5ISC-/- mice by administration of dextran sulfate sodium; disease severity was determined based on body weight, colon length, and histopathology analysis of colon biopsies. RESULTS: BAs and TGR5 agonists promoted growth of intestinal organoids. Administration of cholecystokinin to mice resulted in acute release of BAs into the intestinal lumen and increased proliferation of the intestinal epithelium. BAs and Tgr5 expression in ISCs were required for homeostatic intestinal epithelial renewal and fate specification, and for regeneration after colitis induction. Tgr5ISC-/- mice developed more severe colitis than mice without Tgr5 disruption in ISCs. ISCs incubated with INT-777 increased activation of yes-associated protein 1 (YAP1) and of its upstream regulator SRC. Inhibitors of YAP1 and SRC prevented organoid growth induced by TGR5 activation. CONCLUSIONS: BAs promote regeneration of the intestinal epithelium via activation of TGR5 in ISCs, resulting in activation of SRC and YAP and activation of their target genes. Release of endogenous BAs in the intestinal lumen is sufficient to promote ISC renewal and drives regeneration in response to injury.

cis-regulatory variation modulates susceptibility to enteric infection in the Drosophila genetic reference panel.

BACKGROUND: Resistance to enteric pathogens is a complex trait at the crossroads of multiple biological processes. We have previously shown in the Drosophila Genetic Reference Panel (DGRP) that resistance to infection is highly heritable, but our understanding of how the effects of genetic variants affect different molecular mechanisms to determine gut immunocompetence is still limited. RESULTS: To address this, we perform a systems genetics analysis of the gut transcriptomes from 38 DGRP lines that were orally infected with Pseudomonas entomophila. We identify a large number of condition-specific, expression quantitative trait loci (local-eQTLs) with infection-specific ones located in regions enriched for FOX transcription factor motifs. By assessing the allelic imbalance in the transcriptomes of 19 F1 hybrid lines from a large round robin design, we independently attribute a robust cis-regulatory effect to only 10% of these detected local-eQTLs. However, additional analyses indicate that many local-eQTLs may act in trans instead. Comparison of the transcriptomes of DGRP lines that were either susceptible or resistant to Pseudomonas entomophila infection reveals nutcracker as the only differentially expressed gene. Interestingly, we find that nutcracker is linked to infection-specific eQTLs that correlate with its expression level and to enteric infection susceptibility. Further regulatory analysis reveals one particular eQTL that significantly decreases the binding affinity for the repressor Broad, driving differential allele-specific nutcracker expression. CONCLUSIONS: Our collective findings point to a large number of infection-specific cis- and trans-acting eQTLs in the DGRP, including one common non-coding variant that lowers enteric infection susceptibility.

Enteric infection induces Lark-mediated intron retention at the 5' end of Drosophila genes.

BACKGROUND: RNA splicing is a key post-transcriptional mechanism that generates protein diversity and contributes to the fine-tuning of gene expression, which may facilitate adaptation to environmental challenges. Here, we employ a systems approach to study alternative splicing changes upon enteric infection in females from classical Drosophila melanogaster strains as well as 38 inbred lines. RESULTS: We find that infection leads to extensive differences in isoform ratios, which results in a more diverse transcriptome with longer 5' untranslated regions (5'UTRs). We establish a role for genetic variation in mediating inter-individual splicing differences, with local splicing quantitative trait loci (local-sQTLs) being preferentially located at the 5' end of transcripts and directly upstream of splice donor sites. Moreover, local-sQTLs are more numerous in the infected state, indicating that acute stress unmasks a substantial number of silent genetic variants. We observe a general increase in intron retention concentrated at the 5' end of transcripts across multiple strains, whose prevalence scales with the degree of pathogen virulence. The length, GC content, and RNA polymerase II occupancy of these introns with increased retention suggest that they have exon-like characteristics. We further uncover that retained intron sequences are enriched for the Lark/RBM4 RNA binding motif. Interestingly, we find that lark is induced by infection in wild-type flies, its overexpression and knockdown alter survival, and tissue-specific overexpression mimics infection-induced intron retention. CONCLUSION: Our collective findings point to pervasive and consistent RNA splicing changes, partly mediated by Lark/RBM4, as being an important aspect of the gut response to infection.

Identifying gene function and module connections by the integration of multispecies expression compendia.

The functions of many eukaryotic genes are still poorly understood. Here, we developed and validated a new method, termed GeneBridge, which is based on two linked approaches to impute gene function and bridge genes with biological processes. First, Gene-Module Association Determination (G-MAD) allows the annotation of gene function. Second, Module-Module Association Determination (M-MAD) allows predicting connectivity among modules. We applied the GeneBridge tools to large-scale multispecies expression compendia-1700 data sets with over 300,000 samples from human, mouse, rat, fly, worm, and yeast-collected in this study. G-MAD identifies novel functions of genes-for example, DDT in mitochondrial respiration and WDFY4 in T cell activation-and also suggests novel components for modules, such as for cholesterol biosynthesis. By applying G-MAD on data sets from respective tissues, tissue-specific functions of genes were identified-for instance, the roles of EHHADH in liver and kidney, as well as SLC6A1 in brain and liver. Using M-MAD, we identified a list of module-module associations, such as those between mitochondria and proteasome, mitochondria and histone demethylation, as well as ribosomes and lipid biosynthesis. The GeneBridge tools together with the expression compendia are available as an open resource, which will facilitate the identification of connections linking genes, modules, phenotypes, and diseases.

Commensal Gut Bacteria Buffer the Impact of Host Genetic Variants on Drosophila Developmental Traits under Nutritional Stress.

Eukaryotic genomes encode several buffering mechanisms that robustly maintain invariant phenotypic outcome despite fluctuating environmental conditions. Here we show that the Drosophila gut-associated commensals, represented by a single facultative symbiont, Lactobacillus plantarum (LpWJL), constitutes a so far unexpected buffer that masks the contribution of the host's cryptic genetic variation (CGV) to developmental traits while the host is under nutritional stress. During chronic under-nutrition, LpWJL consistently reduces variation in different host phenotypic traits and ensures robust organ patterning during development; LpWJL also decreases genotype-dependent expression variation, particularly for development-associated genes. We further provide evidence that LpWJL buffers via reactive oxygen species (ROS) signaling whose inhibition impairs microbiota-mediated phenotypic robustness. We thus identified a hitherto unappreciated contribution of the gut facultative symbionts to host fitness that, beyond supporting growth rates and maturation timing, confers developmental robustness and phenotypic homogeneity in times of nutritional stress.

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.

Microfluidics-enabled phenotyping of a whole population of C. elegans worms over their embryonic and post-embryonic development at single-organism resolution.

The organism Caenorhabditis elegans is a performant model system for studying human biological processes and diseases, but until now all phenome data are produced as population-averaged read-outs. Monitoring of individual responses to drug treatments would however be more informative. Here, a new strategy to track different phenotypic traits of individual C. elegans nematodes throughout their full life-cycle-i.e., embryonic and post-embryonic development, until adulthood onset, differently from life-span-is presented. In an automated fashion, single worms were synchronized, isolated, and cultured from egg to adulthood in a microfluidic device, where their identity was preserved during their whole development. Several phenotypes were monitored and quantified for each animal, resulting in high-content phenome data. Specifically, the method was validated by analyzing the response of C. elegans to doxycycline, an antibiotic fairly well-known to prolong the development and activate mitochondrial stress-response pathways in different species. Interestingly, the obtained extensive single-worm phenome not only confirmed the dramatic doxycycline effect on the worm developmental delay, but more importantly revealed subtle yet severe treatment-dependent phenotypes that are representative of minority subgroups and would have otherwise stayed hidden in an averaged dataset. Such heterogeneous response started during the embryonic development, which makes essential having a dedicated chip that allows including this early developmental stage in the drug assay. Our approach would therefore allow elucidating pharmaceutical or therapeutic responses that so far were still being overlooked.

An Integrated Systems Genetics and Omics Toolkit to Probe Gene Function.

Identifying genetic and environmental factors that impact complex traits and common diseases is a high biomedical priority. Here, we developed, validated, and implemented a series of multi-layered systems approaches, including (expression-based) phenome-wide association, transcriptome-/proteome-wide association, and (reverse-) mediation analysis, in an open-access web server (systems-genetics.org) to expedite the systems dissection of gene function. We applied these approaches to multi-omics datasets from the BXD mouse genetic reference population, and identified and validated associations between genes and clinical and molecular phenotypes, including previously unreported links between Rpl26 and body weight, and Cpt1a and lipid metabolism. Furthermore, through mediation and reverse-mediation analysis we established regulatory relations between genes, such as the co-regulation of BCKDHA and BCKDHB protein levels, and identified targets of transcription factors E2F6, ZFP277, and ZKSCAN1. Our multifaceted toolkit enabled the identification of gene-gene and gene-phenotype links that are robust and that translate well across populations and species, and can be universally applied to any populations with multi-omics datasets.

Genetic, molecular and physiological basis of variation in Drosophila gut immunocompetence.

Gut immunocompetence involves immune, stress and regenerative processes. To investigate the determinants underlying inter-individual variation in gut immunocompetence, we perform enteric infection of 140 Drosophila lines with the entomopathogenic bacterium Pseudomonas entomophila and observe extensive variation in survival. Using genome-wide association analysis, we identify several novel immune modulators. Transcriptional profiling further shows that the intestinal molecular state differs between resistant and susceptible lines, already before infection, with one transcriptional module involving genes linked to reactive oxygen species (ROS) metabolism contributing to this difference. This genetic and molecular variation is physiologically manifested in lower ROS activity, lower susceptibility to ROS-inducing agent, faster pathogen clearance and higher stem cell activity in resistant versus susceptible lines. This study provides novel insights into the determinants underlying population-level variability in gut immunocompetence, revealing how relatively minor, but systematic genetic and transcriptional variation can mediate overt physiological differences that determine enteric infection susceptibility.

Randomized codon mutagenesis reveals that the HIV Rev arginine-rich motif is robust to substitutions and that double substitution of two critical residues alters specificity.

The binding of the arginine-rich motif (ARM) of HIV Rev protein to its high-affinity site in stem IIB in the Rev response element (RRE) initiates assembly of a ribonucleoprotein complex that mediates the export of essential, incompletely spliced viral transcripts. Many biochemical, genetic, and structural studies of Rev-RRE IIB have been published, yet the roles of many peptide residues in Rev ARM are unconfirmed by mutagenesis. Rev aptamer I (RAI) is an optimized RRE IIB that binds Rev with higher affinity and for which mutational data are sparse. Randomized-codon libraries of Rev ARM were assayed for their ability to bind RRE IIB and RAI using a bacterial reporter system based on bacteriophage λ N-nut antitermination. Most Rev ARM residues tolerated substitutions without strong loss of binding to RRE IIB, and all except arginine 39 tolerated substitution without strong loss of binding to RAI. The pattern of critical Rev residues is not the same for RRE IIB and RAI, suggesting important differences between the interactions. The results support and aid the interpretation of existing structural models. Observed clinical variation is consistent with additional constraints on Rev mutation. By chance, we found double mutants of two highly critical residues, arginine 35 (to glycine) and asparagine 40 (to valine or lysine), that bind RRE IIB well, but not RAI. That an apparently distinct binding mode occurs with only two mutations highlights the ability of ARMs to evolve new recognition strategies and supports the application of neutral theories of evolution to protein-RNA recognition.