The Hybrid Mouse Diversity Panel: a resource for systems genetics analyses of metabolic and cardiovascular traits

Gene-environment interactions within a precision environmental health framework

Understanding the complex interplay of genetic and environmental factors in disease etiology and the role of gene-environment interactions (GEIs) across human development stages is important. We review the state of GEI research, including challenges in measuring environmental factors and advantages of GEI analysis in understanding disease mechanisms. We discuss the evolution of GEI studies from candidate gene-environment studies to genome-wide interaction studies (GWISs) and the role of multi-omics in mediating GEI effects. We review advancements in GEI analysis methods and the importance of large-scale datasets. We also address the translation of GEI findings into precision environmental health (PEH), showcasing real-world applications in healthcare and disease prevention. Additionally, we highlight societal considerations in GEI research, including environmental justice, the return of results to participants, and data privacy. Overall, we underscore the significance of GEI for disease prediction and prevention and advocate for integrating the exposome into PEH omics studies.

Genetic and environmental interactions contribute to immune variation in rewilded mice

The relative and synergistic contributions of genetics and environment to interindividual immune response variation remain unclear, despite implications in evolutionary biology and medicine. Here we quantify interactive effects of genotype and environment on immune traits by investigating C57BL/6, 129S1 and PWK/PhJ inbred mice, rewilded in an outdoor enclosure and infected with the parasite Trichuris muris. Whereas cellular composition was shaped by interactions between genotype and environment, cytokine response heterogeneity including IFNγ concentrations was primarily driven by genotype with consequence on worm burden. In addition, we show that other traits, such as expression of CD44, were explained mostly by genetics on T cells, whereas expression of CD44 on B cells was explained more by environment across all strains. Notably, genetic differences under laboratory conditions were decreased following rewilding. These results indicate that nonheritable influences interact with genetic factors to shape immune variation and parasite burden.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

The genetic and dietary landscape of the muscle insulin signalling network

Metabolic disease is caused by a combination of genetic and environmental factors, yet few studies have examined how these factors influence signal transduction, a key mediator of metabolism. Using mass spectrometry-based phosphoproteomics, we quantified 23,126 phosphosites in skeletal muscle of five genetically distinct mouse strains in two dietary environments, with and without acute in vivo insulin stimulation. Almost half of the insulin-regulated phosphoproteome was modified by genetic background on an ordinary diet, and high-fat high-sugar feeding affected insulin signalling in a strain-dependent manner. Our data revealed coregulated subnetworks within the insulin signalling pathway, expanding our understanding of the pathway's organisation. Furthermore, associating diverse signalling responses with insulin-stimulated glucose uptake uncovered regulators of muscle insulin responsiveness, including the regulatory phosphosite S469 on Pfkfb2, a key activator of glycolysis. Finally, we confirmed the role of glycolysis in modulating insulin action in insulin resistance. Our results underscore the significance of genetics in shaping global signalling responses and their adaptability to environmental changes, emphasising the utility of studying biological diversity with phosphoproteomics to discover key regulatory mechanisms of complex traits.

Leveraging inter-individual transcriptional correlation structure to infer discrete signaling mechanisms across metabolic tissues

Inter-organ communication is a vital process to maintain physiologic homeostasis, and its dysregulation contributes to many human diseases. Given that circulating bioactive factors are stable in serum, occur naturally, and are easily assayed from blood, they present obvious focal molecules for therapeutic intervention and biomarker development. Recently, studies have shown that secreted proteins mediating inter-tissue signaling could be identified by 'brute force' surveys of all genes within RNA-sequencing measures across tissues within a population. Expanding on this intuition, we reasoned that parallel strategies could be used to understand how individual genes mediate signaling across metabolic tissues through correlative analyses of gene variation between individuals. Thus, comparison of quantitative levels of gene expression relationships between organs in a population could aid in understanding cross-organ signaling. Here, we surveyed gene-gene correlation structure across 18 metabolic tissues in 310 human individuals and 7 tissues in 103 diverse strains of mice fed a normal chow or high-fat/high-sucrose (HFHS) diet. Variation of genes such as FGF21, ADIPOQ, GCG, and IL6 showed enrichments which recapitulate experimental observations. Further, similar analyses were applied to explore both within-tissue signaling mechanisms (liver PCSK9) and genes encoding enzymes producing metabolites (adipose PNPLA2), where inter-individual correlation structure aligned with known roles for these critical metabolic pathways. Examination of sex hormone receptor correlations in mice highlighted the difference of tissue-specific variation in relationships with metabolic traits. We refer to this resource as gene-derived correlations across tissues (GD-CAT) where all tools and data are built into a web portal enabling users to perform these analyses without a single line of code (gdcat.org). This resource enables querying of any gene in any tissue to find correlated patterns of genes, cell types, pathways, and network architectures across metabolic organs.

The role of gonadal hormones and sex chromosomes in sex-dependent effects of early nutrition on metabolic health

Early-life conditions such as prenatal nutrition can have long-term effects on metabolic health, and these effects may differ between males and females. Understanding the biological mechanisms underlying sex differences in the response to early-life environment will improve interventions, but few such mechanisms have been identified, and there is no overall framework for understanding sex differences. Biological sex differences may be due to chromosomal sex, gonadal sex, or interactions between the two. This review describes approaches to distinguish between the roles of chromosomal and gonadal sex, and summarizes findings regarding sex differences in metabolism. The Four Core Genotypes (FCG) mouse model allows dissociation of the sex chromosome genotype from gonadal type, whereas the XY* mouse model can be used to distinguish effects of X chromosome dosage vs the presence of the Y chromosome. Gonadectomy can be used to distinguish between organizational (permanent) and activational (reversible) effects of sex hormones. Baseline sex differences in a variety of metabolic traits are influenced by both activational and organizational effects of gonadal hormones, as well as sex chromosome complement. Thus far, these approaches have not been widely applied to examine sex-dependent effects of prenatal conditions, although a number of studies have found activational effects of estradiol to be protective against the development of hypertension following early-life adversity. Genes that escape X chromosome inactivation (XCI), such as Kdm5c, contribute to baseline sex-differences in metabolism, while Ogt, another XCI escapee, leads to sex-dependent responses to prenatal maternal stress. Genome-wide approaches to the study of sex differences include mapping genetic loci influencing metabolic traits in a sex-dependent manner. Seeking enrichment for binding sites of hormone receptors among genes showing sexually-dimorphic expression can elucidate the relative roles of hormones. Using the approaches described herein to identify mechanisms underlying sex-dependent effects of early nutrition on metabolic health may enable the identification of fundamental mechanisms and potential interventions.

Spatiotemporal-social association predicts immunological similarity in rewilded mice

Environmental influences on immune phenotypes are well-documented, but our understanding of which elements of the environment affect immune systems, and how, remains vague. Behaviors, including socializing with others, are central to an individual's interaction with its environment. We therefore tracked behavior of rewilded laboratory mice of three inbred strains in outdoor enclosures and examined contributions of behavior, including associations measured from spatiotemporal co-occurrences, to immune phenotypes. We found extensive variation in individual and social behavior among and within mouse strains upon rewilding. In addition, we found that the more associated two individuals were, the more similar their immune phenotypes were. Spatiotemporal association was particularly predictive of similar memory T and B cell profiles and was more influential than sibling relationships or shared infection status. These results highlight the importance of shared spatiotemporal activity patterns and/or social networks for immune phenotype and suggest potential immunological correlates of social life.

Sex difference in liver diseases: How preclinical models help to dissect the sex-related mechanisms sustaining NAFLD and hepatocellular carcinoma

Only a few preclinical findings are confirmed in the clinic, posing a critical issue for clinical development. Therefore, identifying the best preclinical models can help to dissect molecular and mechanistic insights into liver disease pathogenesis while being clinically relevant. In this context, the sex relevance of most preclinical models has been only partially considered. This is particularly significant in NAFLD and HCC, which have a higher prevalence in men when compared to pre-menopause women but not to those in post-menopausal status, suggesting a role for sex hormones in the pathogenesis of the diseases. This review gathers the sex-relevant findings and the available preclinical models focusing on both in vitro and in vivo studies and discusses the potential implications and perspectives of introducing the sex effect in the selection of the best preclinical model. This is a critical aspect that would help to tailor personalized therapies based on sex.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

Gene Regulatory Networks in Coronary Artery Disease

PURPOSE OF REVIEW

Coronary artery disease is a complex disorder and the leading cause of mortality worldwide. As technologies for the generation of high-throughput multiomics data have advanced, gene regulatory network modeling has become an increasingly powerful tool in understanding coronary artery disease. This review summarizes recent and novel gene regulatory network tools for bulk tissue and single cell data, existing databases for network construction, and applications of gene regulatory networks in coronary artery disease.

RECENT FINDINGS

New gene regulatory network tools can integrate multiomics data to elucidate complex disease mechanisms at unprecedented cellular and spatial resolutions. At the same time, updates to coronary artery disease expression data in existing databases have enabled researchers to build gene regulatory networks to study novel disease mechanisms. Gene regulatory networks have proven extremely useful in understanding CAD heritability beyond what is explained by GWAS loci and in identifying mechanisms and key driver genes underlying disease onset and progression. Gene regulatory networks can holistically and comprehensively address the complex nature of coronary artery disease. In this review, we discuss key algorithmic approaches to construct gene regulatory networks and highlight state-of-the-art methods that model specific modes of gene regulation. We also explore recent applications of these tools in coronary artery disease patient data repositories to understand disease heritability and shared and distinct disease mechanisms and key driver genes across tissues, between sexes, and between species.

Mechanisms Behind NAFLD: a System Genetics Perspective

PURPOSE OF REVIEW

To summarize the key factors contributing to the onset and progress of nonalcoholic fatty liver disease (NAFLD) and put them in a system genetics context. We particularly focus on how genetic regulation of hepatic lipids contributes to NAFLD.

RECENT FINDINGS

NAFLD is characterized by excessive accumulation of fat in the liver. This can progress to steatohepatitis (inflammation and hepatocyte injury) and eventually, cirrhosis. The severity of NAFLD is determined by a combination of factors including obesity, insulin resistance, and lipotoxic lipids, along with genetic susceptibility. Numerous studies have been conducted on large human cohorts and mouse panels, to identify key determinants in the genome, transcriptome, proteome, lipidome, microbiome and different environmental conditions contributing to NAFLD. We review common factors contributing to NAFLD and put them in a systems genetics context. In particular, we describe how genetic regulation of liver lipids contributes to NAFLD. The combination of an unhealthy lifestyle and genetic predisposition increases the likelihood of accumulating lipotoxic specie lipids that may be one of the driving forces behind developing severe forms of NAFLD.

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.

Genetic pathways regulating the longitudinal acquisition of cocaine self-administration in a panel of inbred and recombinant inbred mice

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

Integrative analysis of hepatic transcriptional profiles reveals genetic regulation of atherosclerosis in hyperlipidemic Diversity Outbred-F1 mice

Atherogenesis is an insipidus but precipitating process leading to serious consequences of many cardiovascular diseases (CVD). Numerous genetic loci contributing to atherosclerosis have been identified in human genome-wide association studies, but these studies have limitations in the ability to control environmental factors and to decipher cause/effect relationships. To assess the power of hyperlipidemic Diversity Outbred (DO) mice in facilitating quantitative trait loci (QTL) analysis of complex traits, we generated a high-resolution genetic panel of atherosclerosis susceptible (DO-F1) mouse cohort by crossing 200 DO females with C57BL/6J males carrying two human genes: encoding apolipoprotein E3-Leiden and cholesterol ester transfer protein. We examined atherosclerotic traits including plasma lipids and glucose in the 235 female and 226 male progeny before and after 16 weeks of a high-fat/cholesterol diet, and aortic plaque size at 24 weeks. We also assessed the liver transcriptome using RNA-sequencing. Our QTL mapping for atherosclerotic traits identified one previously reported female-specific QTL on Chr10 with a narrower interval of 22.73 to 30.80 Mb, and one novel male-specific QTL at 31.89 to 40.25 Mb on Chr19. Liver transcription levels of several genes within each QTL were highly correlated with the atherogenic traits. A majority of these candidates have already known atherogenic potential in humans and/or mice, but integrative QTL, eQTL, and correlation analyses further pointed Ptprk as a major candidate of the Chr10 QTL, while Pten and Cyp2c67 of the Chr19 QTL in our DO-F1 cohort. Finally, through additional analyses of RNA-seq data we identified genetic regulation of hepatic transcription factors, including Nr1h3, contributes to atherogenesis in this cohort. Thus, an integrative approach using DO-F1 mice effectively validates the influence of genetic factors on atherosclerosis in DO mice and suggests an opportunity to discover therapeutics in the setting of hyperlipidemia.

Genetic architecture of heart mitochondrial proteome influencing cardiac hypertrophy

Mitochondria play an important role in both normal heart function and disease etiology. We report analysis of common genetic variations contributing to mitochondrial and heart functions using an integrative proteomics approach in a panel of inbred mouse strains called the Hybrid Mouse Diversity Panel (HMDP). We performed a whole heart proteome study in the HMDP (72 strains, n=2-3 mice) and retrieved 848 mitochondrial proteins (quantified in ≥50 strains). High-resolution association mapping on their relative abundance levels revealed three trans-acting genetic loci on chromosomes (chr) 7, 13 and 17 that regulate distinct classes of mitochondrial proteins as well as cardiac hypertrophy. DAVID enrichment analyses of genes regulated by each of the loci revealed that the chr13 locus was highly enriched for complex-I proteins (24 proteins, P=2.2E-61), the chr17 locus for mitochondrial ribonucleoprotein complex (17 proteins, P=3.1E-25) and the chr7 locus for ubiquinone biosynthesis (3 proteins, P=6.9E-05). Follow-up high resolution regional mapping identified NDUFS4, LRPPRC and COQ7 as the candidate genes for chr13, chr17 and chr7 loci, respectively, and both experimental and statistical analyses supported their causal roles. Furthermore, a large cohort of Diversity Outbred mice was used to corroborate Lrpprc gene as a driver of mitochondrial DNA (mtDNA)-encoded gene regulation, and to show that the chr17 locus is specific to heart. Variations in all three loci were associated with heart mass in at least one of two independent heart stress models, namely, isoproterenol-induced heart failure and diet-induced obesity. These findings suggest that common variations in certain mitochondrial proteins can act in trans to influence tissue-specific mitochondrial functions and contribute to heart hypertrophy, elucidating mechanisms that may underlie genetic susceptibility to heart failure in human populations.

Conserved multi-tissue transcriptomic adaptations to exercise training in humans and mice

Physical activity is associated with beneficial adaptations in human and rodent metabolism. We studied over 50 complex traits before and after exercise intervention in middle-aged men and a panel of 100 diverse strains of female mice. Candidate gene analyses in three brain regions, muscle, liver, heart, and adipose tissue of mice indicate genetic drivers of clinically relevant traits, including volitional exercise volume, muscle metabolism, adiposity, and hepatic lipids. Although ∼33% of genes differentially expressed in skeletal muscle following the exercise intervention are similar in mice and humans independent of BMI, responsiveness of adipose tissue to exercise-stimulated weight loss appears controlled by species and underlying genotype. We leveraged genetic diversity to generate prediction models of metabolic trait responsiveness to volitional activity offering a framework for advancing personalized exercise prescription. The human and mouse data are publicly available via a user-friendly Web-based application to enhance data mining and hypothesis development.

Genetic and Environmental interactions contribute to immune variation in rewilded mice

The relative and synergistic contributions of genetics and environment to inter-individual immune response variation remain unclear, despite its implications for understanding both evolutionary biology and medicine. Here, we quantify interactive effects of genotype and environment on immune traits by investigating three inbred mouse strains rewilded in an outdoor enclosure and infected with the parasite, Trichuris muris. Whereas cytokine response heterogeneity was primarily driven by genotype, cellular composition heterogeneity was shaped by interactions between genotype and environment. Notably, genetic differences under laboratory conditions can be decreased following rewilding, and variation in T cell markers are more driven by genetics, whereas B cell markers are driven more by environment. Importantly, variation in worm burden is associated with measures of immune variation, as well as genetics and environment. These results indicate that nonheritable influences interact with genetic factors to shape immune variation, with synergistic impacts on the deployment and evolution of defense mechanisms.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Isometric skeletal muscle contractile properties in common strains of male laboratory mice

Assessing contractile function of skeletal muscle in murine models is a commonly employed laboratory technique that investigators utilize to measure the impact of genetic manipulations, drug efficacy, or other therapeutic interventions. Often overlooked is the potential for the strain of the mouse to influence the functional properties of the skeletal muscle. Thus, we sought to characterize commonly assessed isometric force measures in the hindlimb muscles across a variety of mouse strains. Using 6-8-week-old male mice, we measured isometric force, fatigue susceptibility, relaxation kinetics, muscle mass, myofiber cross-sectional area, and fiber type composition of the extensor digitorum longus (EDL) and soleus muscles in C57BL/6NJ, BALB/cJ, FVB/NJ, C57BL/6J, and C57BL/10 mice. The data demonstrate both unique differences and a number of similarities between both muscles in the various genetic backgrounds of mice. Soleus muscle specific force (i.e., force per unit size) exhibited higher variation across strains while specific force of the EDL muscle exhibited minimal variation. In contrast, absolute force differed only in a few mouse strains whereas analysis of muscle morphology revealed many distinctions when compared across all the groups. Collectively, the data suggest that the strain of the mouse can potentially influence the measured biological outcome and may possibly promote a synergistic effect with any genetic manipulation or therapeutic intervention. Thus, it is critical for the investigator to carefully consider the genetic background of the mouse used in the experimental design and precisely document the strain of mouse employed during publication.

Liver-heart cross-talk mediated by coagulation factor XI protects against heart failure

Tissue-tissue communication by endocrine factors is a vital mechanism for physiologic homeostasis. A systems genetics analysis of transcriptomic and functional data from a cohort of diverse, inbred strains of mice predicted that coagulation factor XI (FXI), a liver-derived protein, protects against diastolic dysfunction, a key trait of heart failure with preserved ejection fraction. This was confirmed using gain- and loss-of-function studies, and FXI was found to activate the bone morphogenetic protein (BMP)-SMAD1/5 pathway in the heart. The proteolytic activity of FXI is required for the cleavage and activation of extracellular matrix-associated BMP7 in the heart, thus inhibiting genes involved in inflammation and fibrosis. Our results reveal a protective role of FXI in heart injury that is distinct from its role in coagulation.

Evaluation and characterization of expression quantitative trait analysis methods in the Hybrid Rat Diversity Panel

The Hybrid Rat Diversity Panel (HRDP) is a stable and well-characterized set of more than 90 inbred rat strains that can be leveraged for systems genetics approaches to understanding the genetic and genomic variation associated with complex disease. The HRDP exhibits substantial between-strain diversity while retaining substantial within-strain isogenicity, allowing for the precise mapping of genetic variation associated with complex phenotypes and providing statistical power to identify associated variants. In order to robustly identify associated genetic variants, it is important to account for the population structure induced by inbreeding. To this end, we investigate the performance of four plausible approaches towards modeling quantitative traits in the HRDP and quantify their operating characteristics. In particular, we investigate three approaches based on genome-wide mixed model analysis, and one approach based on ordinary least squares linear regression. Towards facilitating study planning and design, we conduct extensive simulations to investigate the power of genetic association analyses in the HRDP, and characterize the impressive attained power. In simulation of eQTL data in the HRDP, we find that a mixed model approach that leverages leave-one-chromosome-out kinship estimation attains the highest power while controlling type I error.

Identification of the Transcription Factor ATF3 as a Direct and Indirect Regulator of the LDLR

Coronary artery disease (CAD) is a complex, multifactorial disease caused, in particular, by inflammation and cholesterol metabolism. At the molecular level, the role of tissue-specific signaling pathways leading to CAD is still largely unexplored. This study relied on two main resources: (1) genes with impact on atherosclerosis/CAD, and (2) liver-specific transcriptome analyses from human and mouse studies. The transcription factor activating transcription factor 3 (ATF3) was identified as a key regulator of a liver network relevant to atherosclerosis and linked to inflammation and cholesterol metabolism. ATF3 was predicted to be a direct and indirect (via MAF BZIP Transcription Factor F (MAFF)) regulator of low-density lipoprotein receptor (LDLR). Chromatin immunoprecipitation DNA sequencing (ChIP-seq) data from human liver cells revealed an ATF3 binding motif in the promoter regions of MAFF and LDLR. siRNA knockdown of ATF3 in human Hep3B liver cells significantly upregulated LDLR expression (p < 0.01). Inflammation induced by lipopolysaccharide (LPS) stimulation resulted in significant upregulation of ATF3 (p < 0.01) and subsequent downregulation of LDLR (p < 0.001). Liver-specific expression data from human CAD patients undergoing coronary artery bypass grafting (CABG) surgery (STARNET) and mouse models (HMDP) confirmed the regulatory role of ATF3 in the homeostasis of cholesterol metabolism. This study suggests that ATF3 might be a promising treatment candidate for lowering LDL cholesterol and reducing cardiovascular risk.

Cochlear ribbon synapse maturation requires Nlgn1 and Nlgn3

Hearing depends on precise synaptic transmission between cochlear inner hair cells and spiral ganglion neurons through afferent ribbon synapses. Neuroligins (Nlgns) facilitate synapse maturation in the brain, but they have gone unstudied in the cochlea. We report Nlgn3 and Nlgn1 knockout (KO) cochleae have fewer ribbon synapses and have impaired hearing. Nlgn3 KO is more vulnerable to noise trauma with limited activity at high frequencies one day after noise. Furthermore, Nlgn3 KO cochleae have a 5-fold reduction in synapse number compared to wild type after two weeks of recovery. Double KO cochlear phenotypes are more prominent than the KOs, for example, 5-fold smaller synapses, 25% reduction in synapse density, and 30% less synaptic output. These observations indicate Nlgn3 and Nlgn1 are essential to cochlear ribbon synapse maturation and function.

Mouse models of atherosclerosis in translational research

Atherosclerotic cardiovascular disease (CVD), the major cause of premature human mortality, is a chronic and progressive metabolic and inflammatory disease in large- and medium-sized arteries. Mouse models are widely used to gain mechanistic insights into the pathogenesis of atherosclerosis and have facilitated the discovery of anti-atherosclerotic drugs. Despite promising preclinical studies, many drug candidates have not translated to clinical use because of the complexity of disease patho-mechanisms including lipid metabolic traits and inflammatory, genetic, and hemodynamic factors. We review the current preclinical utility and translation potential of traditional [apolipoprotein E (APOE)- and low-density lipoprotein (LDL) receptor (LDLR)-deficient mice] and emerging mouse models that include partial carotid ligation and AAV8-Pcsk9-D377Y injection in atherosclerosis research and drug discovery. This article represents an important resource in atherosclerosis research.

Genetic Background Matters: Population-Based Studies in Model Organisms for Translational Research

We are all similar but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly, and yeast panels with insightful applications for translational research.

Sex-specific changes in metabolism during the transition from chow to high-fat diet feeding are abolished in response to dieting in C57BL/6J mice

Background/Objective

Female mice are often excluded from diet-induced obesity studies as they are more resistant to the obesifying effects of a high-fat diet (HFD). However, the underlying mechanisms behind this sex disparity may actually have important implications for the development and management of obesity in humans. Therefore, we systematically investigated the immediate sex-specific effects of transitioning to a HFD in C57BL/6J mice as well as monitored whether these effects are altered after sustained HFD feeding and whether sex affects the response to a return to chow, representative of dieting.

Methods

Dual X-ray absorptiometry (DXA) analysis of body composition, indirect calorimetry measurements, and qPCR analysis of hypothalamic and brainstem regions were performed on male and female C57BL/6J mice.

Results

HFD had immediate and dramatic effects in males, increasing fat mass by 58% in the first 3 days. The resistance to the obesifying effect of HFD in females was linked both to an ability to maintain activity levels as well as to an immediate and significantly enhanced reduction in respiratory quotient (RQ), suggesting a greater ability to utilise fat in the diet as a source of fuel. Mechanistically, this sex disparity may be at least partially due to inherent sex differences in the catabolic (POMC/CART) versus anabolic (NPY/AgRP) neurological signalling pathways. Interestingly, the reintroduction of chow following HFD had immediate and consistent responses between the sexes with body composition and most metabolic parameters normalised within 3 days. However, both sexes displayed elevated hypothalamic Npy levels reminiscent of starvation. The difference in RQ seen between the sexes on HFD was immediately abolished suggesting similar abilities to burn fat reserves for fuel.

Conclusions

C57BL/6J mice have markedly different sex-specific behavioural and metabolic responses to the introduction as well as the sustained intake of a HFD, but consistent responses to a dieting situation.

Sex differences in heart mitochondria regulate diastolic dysfunction

Heart failure with preserved ejection fraction (HFpEF) exhibits a sex bias, being more common in women than men, and we hypothesize that mitochondrial sex differences might underlie this bias. As part of genetic studies of heart failure in mice, we observe that heart mitochondrial DNA levels and function tend to be reduced in females as compared to males. We also observe that expression of genes encoding mitochondrial proteins are higher in males than females in human cohorts. We test our hypothesis in a panel of genetically diverse inbred strains of mice, termed the Hybrid Mouse Diversity Panel (HMDP). Indeed, we find that mitochondrial gene expression is highly correlated with diastolic function, a key trait in HFpEF. Consistent with this, studies of a "two-hit" mouse model of HFpEF confirm that mitochondrial function differs between sexes and is strongly associated with a number of HFpEF traits. By integrating data from human heart failure and the mouse HMDP cohort, we identify the mitochondrial gene Acsl6 as a genetic determinant of diastolic function. We validate its role in HFpEF using adenoviral over-expression in the heart. We conclude that sex differences in mitochondrial function underlie, in part, the sex bias in diastolic function.

Illuminating the Mechanisms Underlying Sex Differences in Cardiovascular Disease

Sex is a key risk factor for many types of cardiovascular disease. It is imperative to understand the mechanisms underlying sex differences to devise optimal preventive and therapeutic approaches for all individuals. Both biological sex (determined by sex chromosomes and gonadal hormones) and gender (social and cultural behaviors associated with femininity or masculinity) influence differences between men and women in disease susceptibility and pathology. Here, we focus on the application of experimental mouse models that elucidate the influence of 2 components of biological sex-sex chromosome complement (XX or XY) and gonad type (ovaries or testes). These models have revealed that in addition to well-known effects of gonadal hormones, sex chromosome complement influences cardiovascular risk factors, such as plasma cholesterol levels and adiposity, as well as the development of atherosclerosis and pulmonary hypertension. One mechanism by which sex chromosome dosage influences cardiometabolic traits is through sex-biased expression of X chromosome genes that escape X inactivation. These include chromatin-modifying enzymes that regulate gene expression throughout the genome. The identification of factors that determine sex-biased gene expression and cardiometabolic traits will expand our mechanistic understanding of cardiovascular disease processes and provide insight into sex differences that remain throughout the lifespan as gonadal hormone levels alter with age.

Hearing loss and tinnitus: association studies for complex-hearing disorders in mouse and man

Genome-wide association studies (GWAS) provide an unbiased first look at genetic loci involved in aging and noise-induced sensorineural hearing loss and tinnitus. The hearing phenotype, whether audiogram-based or self-report, is regressed against genotyped information at representative single nucleotide polymorphisms (SNPs) across the genome. Findings include the fact that both hearing loss and tinnitus are polygenic disorders, with up to thousands of genes, each of effect size of < 0.02. Smaller human GWAS' were able to use objective measures and identified a few loci; however, hundreds of thousands of participants have been required for the statistical power to identify significant variants, and GWAS is unable to assess rare variants with mean allele frequency < 1%. Animal studies are required as well because of inability to access the human cochlea. Mouse GWAS builds on linkage techniques and the known phenotypic differences in auditory function between inbred strains. With the advantage that the laboratory environment can be controlled for noise and aging, the Hybrid Mouse Diversity Panel (HDMP) combines 100 strains sequenced at high resolution. Lift-over regions between mice and humans have identified over 17,000 homologous genes. Since most significant SNPs are either intergenic or in introns, and binding sites between species are poorly preserved between species, expression quantitative trait locus information is required to bring humans and mice into agreement. Transcriptome-wide analysis studies (TWAS) can prioritize putative causal genes and tissues. Diverse species, each making a distinct contribution, carry a synergistic advantage in the quest for treatment and ultimate cure of sensorineural hearing difficulties.

Atherosclerosis Induced by Adeno-Associated Virus Encoding Gain-of-Function PCSK9

Induction of atherosclerosis in mice with one or more genetic alterations (e.g., conditional deletion of a gene of interest) has traditionally required crossbreeding with Apoe or Ldlr deficient mice to achieve sufficient hypercholesterolemia. However, this procedure is time consuming and generates a surplus of mice with genotypes that are irrelevant for experiments. Several alternative methods exist that obviate the need to work in mice with germline-encoded hypercholesterolemia. In this chapter, we detail an efficient and increasingly used method to induce hypercholesterolemia in mice through adeno-associated virus-mediated transfer of the proprotein convertase subtilisin/kexin type 9 (PCSK9) gene.

Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.

Placental genomics mediates genetic associations with complex health traits and disease

As the master regulator in utero, the placenta is core to the Developmental Origins of Health and Disease (DOHaD) hypothesis but is historically understudied. To identify placental gene-trait associations (GTAs) across the life course, we perform distal mediator-enriched transcriptome-wide association studies (TWAS) for 40 traits, integrating placental multi-omics from the Extremely Low Gestational Age Newborn Study. At $$P \; < \; 2.5\times {10}^{-6}$$ P < 2.5 × 10 − 6 , we detect 248 GTAs, mostly for neonatal and metabolic traits, across 176 genes, enriched for cell growth and immunological pathways. In aggregate, genetic effects mediated by placental expression significantly explain 4 early-life traits but no later-in-life traits. 89 GTAs show significant mediation through distal genetic variants, identifying hypotheses for distal regulation of GTAs. Investigation of one hypothesis in human placenta-derived choriocarcinoma cells reveal that knockdown of mediator gene EPS15 upregulates predicted targets SPATA13 and FAM214A, both associated with waist-hip ratio in TWAS, and multiple genes involved in metabolic pathways. These results suggest profound health impacts of placental genomic regulation in developmental programming across the life course.

Stem Cell-Derived β Cells: A Versatile Research Platform to Interrogate the Genetic Basis of β Cell Dysfunction

Pancreatic β cell dysfunction is a central component of diabetes progression. During the last decades, the genetic basis of several monogenic forms of diabetes has been recognized. Genome-wide association studies (GWAS) have also facilitated the identification of common genetic variants associated with an increased risk of diabetes. These studies highlight the importance of impaired β cell function in all forms of diabetes. However, how most of these risk variants confer disease risk, remains unanswered. Understanding the specific contribution of genetic variants and the precise role of their molecular effectors is the next step toward developing treatments that target β cell dysfunction in the era of personalized medicine. Protocols that allow derivation of β cells from pluripotent stem cells, represent a powerful research tool that allows modeling of human development and versatile experimental designs that can be used to shed some light on diabetes pathophysiology. This article reviews different models to study the genetic basis of β cell dysfunction, focusing on the recent advances made possible by stem cell applications in the field of diabetes research.

A mechanistic framework for cardiometabolic and coronary artery diseases

Coronary atherosclerosis results from the delicate interplay of genetic and exogenous risk factors, principally taking place in metabolic organs and the arterial wall. Here we show that 224 gene-regulatory coexpression networks (GRNs) identified by integrating genetic and clinical data from patients with (n = 600) and without (n = 250) coronary artery disease (CAD) with RNA-seq data from seven disease-relevant tissues in the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) study largely capture this delicate interplay, explaining >54% of CAD heritability. Within 89 cross-tissue GRNs associated with clinical severity of CAD, 374 endocrine factors facilitated inter-organ interactions, primarily along an axis from adipose tissue to the liver (n = 152). This axis was independently replicated in genetically diverse mouse strains and by injection of recombinant forms of adipose endocrine factors (EPDR1, FCN2, FSTL3 and LBP) that markedly altered blood lipid and glucose levels in mice. Altogether, the STARNET database and the associated GRN browser (http://starnet.mssm.edu) provide a multiorgan framework for exploration of the molecular interplay between cardiometabolic disorders and CAD.

Proteome-wide systems genetics identifies UFMylation as a regulator of skeletal muscle function

Improving muscle function has great potential to improve the quality of life. To identify novel regulators of skeletal muscle metabolism and function, we performed a proteomic analysis of gastrocnemius muscle from 73 genetically distinct inbred mouse strains, and integrated the data with previously acquired genomics and >300 molecular/phenotypic traits via quantitative trait loci mapping and correlation network analysis. These data identified thousands of associations between protein abundance and phenotypes and can be accessed online (https://muscle.coffeeprot.com/) to identify regulators of muscle function. We used this resource to prioritize targets for a functional genomic screen in human bioengineered skeletal muscle. This identified several negative regulators of muscle function including UFC1, an E2 ligase for protein UFMylation. We show UFMylation is up-regulated in a mouse model of amyotrophic lateral sclerosis, a disease that involves muscle atrophy. Furthermore, in vivo knockdown of UFMylation increased contraction force, implicating its role as a negative regulator of skeletal muscle function.

Transcriptome-wide association study of coronary artery disease identifies novel susceptibility genes

The majority of risk loci identified by genome-wide association studies (GWAS) are in non-coding regions, hampering their functional interpretation. Instead, transcriptome-wide association studies (TWAS) identify gene-trait associations, which can be used to prioritize candidate genes in disease-relevant tissue(s). Here, we aimed to systematically identify susceptibility genes for coronary artery disease (CAD) by TWAS. We trained prediction models of nine CAD-relevant tissues using EpiXcan based on two genetics-of-gene-expression panels, the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) and the Genotype-Tissue Expression (GTEx). Based on these prediction models, we imputed gene expression of respective tissues from individual-level genotype data on 37,997 CAD cases and 42,854 controls for the subsequent gene-trait association analysis. Transcriptome-wide significant association (i.e. P < 3.85e-6) was observed for 114 genes. Of these, 96 resided within previously identified GWAS risk loci and 18 were novel. Stepwise analyses were performed to study their plausibility, biological function, and pathogenicity in CAD, including analyses for colocalization, damaging mutations, pathway enrichment, phenome-wide associations with human data and expression-traits correlations using mouse data. Finally, CRISPR/Cas9-based gene knockdown of two newly identified TWAS genes, RGS19 and KPTN, in a human hepatocyte cell line resulted in reduced secretion of APOB100 and lipids in the cell culture medium. Our CAD TWAS work (i) prioritized candidate causal genes at known GWAS loci, (ii) identified 18 novel genes to be associated with CAD, and iii) suggested potential tissues and pathways of action for these TWAS CAD genes.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Thymus-Derived CD4+CD8+ Cells Reside in Mediastinal Adipose Tissue and the Aortic Arch

Double-positive CD4⁺CD8αβ⁺ (DP) cells are thought to reside as T cell progenitors exclusively within the thymus. We recently discovered an unexpected CD4⁺ and CD8αβ⁺ immune cell population in healthy and atherosclerotic mice by single-cell RNA sequencing. Transcriptomically, these cells resembled thymic DPs. Flow cytometry and three-dimensional whole-mount imaging confirmed DPs in thymus, mediastinal adipose tissue, and aortic adventitia, but nowhere else. Deep transcriptional profiling revealed differences between DP cells isolated from the three locations. All DPs were dependent on RAG2 expression and the presence of the thymus. Mediastinal adipose tissue DPs resided in close vicinity to invariant NKT cells, which they could activate in vitro. Thymus transplantation failed to reconstitute extrathymic DPs, and frequencies of extrathymic DPs were unaltered by pharmacologic inhibition of S1P1, suggesting that their migration may be locally confined. Our results define two new, transcriptionally distinct subsets of extrathymic DPs that may play a role in aortic vascular homeostasis.

Matricellular Protein Cilp1 Promotes Myocardial Fibrosis in Response to Myocardial Infarction

[Figure: see text].

Integrative Analyses Reveal Tstd1 as a Potential Modulator of HDL Cholesterol and Mitochondrial Function in Mice

High-density lipoprotein (HDL) cholesterol levels are closely associated with human health and diseases. To identify genes modulating plasma HDL levels, we integrated HDL measurements and multi-omics data collected from diverse mouse cohorts and combined a list of systems genetics methods, including quantitative trait loci (QTL) mapping analysis, mediation analysis, transcriptome-wide association analysis (TWAS), and correlation analysis. We confirmed a significant and conserved QTL for plasma HDL on chromosome 1 and identified that Tstd1 liver transcript correlates with plasma HDL in several independent mouse cohorts, suggesting Tstd1 may be a potential modulator of plasma HDL levels. Correlation analysis using over 70 transcriptomics datasets in humans and mice revealed consistent correlations between Tstd1 and genes known to be involved in cholesterol and HDL regulation. Consistent with strong enrichment in gene sets related to cholesterol and lipoproteins in the liver, mouse strains with high Tstd1 exhibited higher plasma levels of HDL, total cholesterol and other lipid markers. GeneBridge using large-scale expression datasets identified conserved and positive associations between TSTD1/Tstd1 and mitochondrial pathways, as well as cholesterol and lipid pathways in human, mouse and rat. In summary, we identified Tstd1 as a new modulator of plasma HDL and mitochondrial function through integrative systems analyses, and proposed a new mechanism of HDL modulation and a potential therapeutic target for relevant diseases. This study highlights the value of such integrative approaches in revealing molecular mechanisms of complex traits or diseases.

Pcpe2, a Novel Extracellular Matrix Protein, Regulates Adipocyte SR-BI-Mediated High-Density Lipoprotein Uptake

Objective

To investigate the role of adipocyte Pcpe2 (procollagen C-endopeptidase enhancer 2) in SR-BI (scavenger receptor class BI)-mediated HDL-C (high-density lipoprotein cholesterol) uptake and contributions to adipose lipid storage.

Approach and Results

Pcpe2, a glycoprotein devoid of intrinsic proteolytic activity, is believed to participate in extracellular protein-protein interactions, supporting SR-BI- mediated HDL-C uptake. In published studies, Pcpe2 deficiency increased the development of atherosclerosis by reducing SR-BI-mediated HDL-C catabolism, but the biological impact of this deficiency on adipocyte SR-BI-mediated HDL-C uptake is unknown. Differentiated cells from Ldlr-/-/Pcpe2-/- (Pcpe2-/-) mouse adipose tissue showed elevated SR-BI protein levels, but significantly reduced HDL-C uptake compared to Ldlr-/- (control) adipose tissue. SR-BI-mediated HDL-C uptake was restored by preincubation of cells with exogenous Pcpe2. In diet-fed mice lacking Pcpe2, significant reductions in visceral, subcutaneous, and brown adipose tissue mass were observed, despite elevations in plasma triglyceride and cholesterol concentrations. Significant positive correlations exist between adipose mass and Pcpe2 expression in both mice and humans.

Conclusions

Overall, these findings reveal a novel and unexpected function for Pcpe2 in modulating SR-BI expression and function as it relates to adipose tissue expansion and cholesterol balance in both mice and humans.

Gene-Environment Interactions for Cardiovascular Disease

PURPOSE OF REVIEW

We provide an overview of recent findings with respect to gene-environment (GxE) interactions for cardiovascular disease (CVD) risk and discuss future opportunities for advancing the field.

RECENT FINDINGS

Over the last several years, GxE interactions for CVD have mostly been identified for smoking and coronary artery disease (CAD) or related risk factors. By comparison, there is more limited evidence for GxE interactions between CVD outcomes and other exposures, such as physical activity, air pollution, diet, and sex. The establishment of large consortia and population-based cohorts, in combination with new computational tools and mouse genetics platforms, can potentially overcome some of the limitations that have hindered human GxE interaction studies and reveal additional association signals for CVD-related traits. The identification of novel GxE interactions is likely to provide a better understanding of the pathogenesis and genetic liability of CVD, with significant implications for healthy lifestyles and therapeutic strategies.

Integration of feeding behavior by the liver circadian clock reveals network dependency of metabolic rhythms

The mammalian circadian clock, expressed throughout the brain and body, controls daily metabolic homeostasis. Clock function in peripheral tissues is required, but not sufficient, for this task. Because of the lack of specialized animal models, it is unclear how tissue clocks interact with extrinsic signals to drive molecular oscillations. Here, we isolated the interaction between feeding and the liver clock by reconstituting Bmal1 exclusively in hepatocytes (Liver-RE), in otherwise clock-less mice, and controlling timing of food intake. We found that the cooperative action of BMAL1 and the transcription factor CEBPB regulates daily liver metabolic transcriptional programs. Functionally, the liver clock and feeding rhythm are sufficient to drive temporal carbohydrate homeostasis. By contrast, liver rhythms tied to redox and lipid metabolism required communication with the skeletal muscle clock, demonstrating peripheral clock cross-talk. Our results highlight how the inner workings of the clock system rely on communicating signals to maintain daily metabolism.

Anterograde regulation of mitochondrial genes and FGF21 signaling by hepatic LSD1

Mitochondrial biogenesis and function are controlled by anterograde regulatory pathways involving more than 1000 nuclear-encoded proteins. Transcriptional networks controlling the nuclear-encoded mitochondrial genes remain to be fully elucidated. Here, we show that histone demethylase LSD1 KO from adult mouse liver (LSD1-LKO) reduces the expression of one-third of all nuclear-encoded mitochondrial genes and decreases mitochondrial biogenesis and function. LSD1-modulated histone methylation epigenetically regulates nuclear-encoded mitochondrial genes. Furthermore, LSD1 regulates gene expression and protein methylation of nicotinamide mononucleotide adenylyltransferase 1 (NMNAT1), which controls the final step of NAD+ synthesis and limits NAD+ availability in the nucleus. Lsd1 KO reduces NAD+-dependent SIRT1 and SIRT7 deacetylase activity, leading to hyperacetylation and hypofunctioning of GABPβ and PGC-1α, the major transcriptional factor/cofactor for nuclear-encoded mitochondrial genes. Despite the reduced mitochondrial function in the liver, LSD1-LKO mice are protected from diet-induced hepatic steatosis and glucose intolerance, partially due to induction of hepatokine FGF21. Thus, LSD1 orchestrates a core regulatory network involving epigenetic modifications and NAD+ synthesis to control mitochondrial function and hepatokine production.

Identification of genetic modifiers of murine hepatic β-glucocerebrosidase activity

The acid β-glucocerebrosidase (GCase) enzyme cleaves glucosylceramide into glucose and ceramide. Loss of function variants in the gene encoding for GCase can lead to Gaucher disease and Parkinson's disease. Therapeutic strategies aimed at increasing GCase activity by targeting a modulating factor are attractive and poorly explored. To identify genetic modifiers, we measured hepatic GCase activity in 27 inbred mouse strains. A genome-wide association study (GWAS) using GCase activity as a trait identified several candidate modifier genes, including Dmrtc2 and Arhgef1 (p=2.1x10⁻⁷), and Grik5 (p=2.1x10⁻⁷). Bayesian integration of the gene mapping with transcriptomics was used to build integrative networks. The analysis uncovered additional candidate GCase regulators, highlighting modules of the acute phase response (p=1.01x10⁻⁸), acute inflammatory response (p=1.01x10⁻⁸), fatty acid beta-oxidation (p=7.43x10⁻⁵), among others. Our study revealed previously unknown candidate modulators of GCase activity, which may facilitate the design of therapies for diseases with GCase dysfunction.

•

Hepatic GCase activity significantly differs among mouse strains.

•

Genome-wide association study revealed putative modifier genes of GCase activity.

•

Bayesian integration of multi-omics identified a regulatory network of GCase activity.

•

This study may facilitate the design of therapies for diseases with GCase dysfunction.

Synergizing Mouse and Human Studies to Understand the Heterogeneity of Obesity

Obesity is routinely considered as a single disease state, which drives a "one-size-fits-all" approach to treatment. We recently convened the first annual University of North Carolina Interdisciplinary Nutrition Sciences Symposium to discuss the heterogeneity of obesity and the need for translational science to advance understanding of this heterogeneity. The symposium aimed to advance scientific rigor in translational studies from animal to human models with the goal of identifying underlying mechanisms and treatments. In this review, we discuss fundamental gaps in knowledge of the heterogeneity of obesity ranging from cellular to population perspectives. We also advocate approaches to overcoming limitations in the field. Examples include the use of contemporary mouse genetic reference population models such as the Collaborative Cross and Diversity Outbred mice that effectively model human genetic diversity and the use of translational models that integrate -omics and computational approaches from pre-clinical to clinical models of obesity. Finally, we suggest best scientific practices to ensure strong rigor that will allow investigators to delineate the sources of heterogeneity in the population with obesity. Collectively, we propose that it is critical to think of obesity as a heterogeneous disease with complex mechanisms and etiologies, requiring unique prevention and treatment strategies tailored to the individual.

Transcriptional circuitry atlas of genetic diverse unstimulated murine and human macrophages define disparity in population-wide innate immunity

Macrophages are ubiquitous custodians of tissues, which play decisive role in maintaining cellular homeostasis through regulatory immune responses. Within tissues, macrophage exhibit extremely heterogeneous population with varying functions orchestrated through regulatory response, which can be further exacerbated in diverse genetic backgrounds. Gene regulatory networks (GRNs) offer comprehensive understanding of cellular regulatory behavior by unfolding the transcription factors (TFs) and regulated target genes. RNA-Seq coupled with ATAC-Seq has revolutionized the regulome landscape influenced by gene expression modeling. Here, we employ an integrative multi-omics systems biology-based analysis and generated GRNs derived from the unstimulated bone marrow-derived macrophages of five inbred genetically defined murine strains, which are reported to be linked with most of the population-wide human genetic variants. Our probabilistic modeling of a basal hemostasis pan regulatory repertoire in diverse macrophages discovered 96 TFs targeting 6279 genes representing 468,291 interactions across five inbred murine strains. Subsequently, we identify core and distinctive GRN sub-networks in unstimulated macrophages to describe the system-wide conservation and dissimilarities, respectively across five murine strains. Our study concludes that discrepancies in unstimulated macrophage-specific regulatory networks not only drives the basal functional plasticity within genetic backgrounds, additionally aid in understanding the complexity of racial disparity among the human population during stress.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.