An Integrated Systems Genetics and Omics Toolkit to Probe Gene Function

Scrutinizing neurodegenerative diseases: decoding the complex genetic architectures through a multi-omics lens

Neurodegenerative diseases present complex genetic architectures, reflecting a continuum from monogenic to oligogenic and polygenic models. Recent advances in multi-omics data, coupled with systems genetics, have significantly refined our understanding of how these data impact neurodegenerative disease mechanisms. To contextualize these genetic discoveries, we provide a comprehensive critical overview of genetic architecture concepts, from Mendelian inheritance to the latest insights from oligogenic and omnigenic models. We explore the roles of common and rare genetic variants, gene-gene and gene-environment interactions, and epigenetic influences in shaping disease phenotypes. Additionally, we emphasize the importance of multi-omics layers including genomic, transcriptomic, proteomic, epigenetic, and metabolomic data in elucidating the molecular mechanisms underlying neurodegeneration. Special attention is given to missing heritability and the contribution of rare variants, particularly in the context of pleiotropy and network pleiotropy. We examine the application of single-cell omics technologies, transcriptome-wide association studies, and epigenome-wide association studies as key approaches for dissecting disease mechanisms at tissue- and cell-type levels. Our review introduces the OmicPeak Disease Trajectory Model, a conceptual framework for understanding the genetic architecture of neurodegenerative disease progression, which integrates multi-omics data across biological layers and time points. This review highlights the critical importance of adopting a systems genetics approach to unravel the complex genetic architecture of neurodegenerative diseases. Finally, this emerging holistic understanding of multi-omics data and the exploration of the intricate genetic landscape aim to provide a foundation for establishing more refined genetic architectures of these diseases, enhancing diagnostic precision, predicting disease progression, elucidating pathogenic mechanisms, and refining therapeutic strategies for neurodegenerative conditions.

Unlocking metabolic insights with mouse genetic diversity

As part of EMBO Journal’s 2024 metabolism methods series, this commentary revisits the impact of genetics on metabolic studies, enabling dissection of novel mechanisms and phenotypes.

High-content phenotypic analysis of a C. elegans recombinant inbred population identifies genetic and molecular regulators of lifespan

Lifespan is influenced by complex interactions between genetic and environmental factors. Studying those factors in model organisms of a single genetic background limits their translational value for humans. Here, we mapped lifespan determinants in 85 C. elegans recombinant inbred advanced intercross lines (RIAILs). We assessed molecular profiles-transcriptome, proteome, and lipidome-and life-history traits, including lifespan, development, growth dynamics, and reproduction. RIAILs exhibited large variations in lifespan, which correlated positively with developmental time. We validated three longevity modulators, including rict-1, gfm-1, and mltn-1, among the top candidates obtained from multiomics data integration and quantitative trait locus (QTL) mapping. We translated their relevance to humans using UK Biobank data and showed that variants in GFM1 are associated with an elevated risk of age-related heart failure. We organized our dataset as a resource that allows interactive explorations for new longevity targets.

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.

TEAD1, MYO7A and NDUFC2 are novel functional genes associated with glucose metabolism in BXD recombinant inbred population

AIM

The liver is an important metabolic organ that governs glucolipid metabolism, and its dysfunction may cause non-alcoholic fatty liver disease, type 2 diabetes mellitus, dyslipidaemia, etc. We aimed to systematic investigate the key factors related to hepatic glucose metabolism, which may be beneficial for understanding the underlying pathogenic mechanisms for obesity and diabetes mellitus.

MATERIALS AND METHODS

Oral glucose tolerance test (OGTT) phenotypes and liver transcriptomes of BXD mice under chow and high-fat diet conditions were collected from GeneNetwork. QTL mapping was conducted to pinpoint genomic regions associated with glucose homeostasis. Candidate genes were further nominated using a multi-criteria approach and validated to confirm their functional relevance in vitro.

RESULTS

Our results demonstrated that plasma glucose levels in OGTT were significantly affected by both diet and genetic background, with six genetic regulating loci were mapped on chromosomes 1, 4, and 7. Moreover, TEAD1, MYO7A and NDUFC2 were identified as the candidate genes. Functionally, siRNA-mediated TEAD1, MYO7A and NDUFC2 knockdown significantly decreased the glucose uptake and inhibited the transcription of genes related to insulin and glucose metabolism pathways.

CONCLUSIONS

Our study contributes novel insights to the understanding of hepatic glucose metabolism, demonstrating the impact of TEAD1, MYO7A and NDUFC2 on mitochondrial function in the liver and their regulatory role in maintaining in glucose homeostasis.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

The Combination of Quantitative Proteomics and Systems Genetics Analysis Reveals that PTN Is Associated with Sleep-Loss-Induced Cognitive Impairment

Sleep loss is associated with cognitive dysfunction. However, the detailed mechanisms remain unclear. In this study, we established a para-chlorophenylalanine (PCPA)-induced insomniac mouse model with impaired cognitive function. Mass-spectrometry-based proteomics showed that the expression of 164 proteins was significantly altered in the hippocampus of the PCPA mice. To identify critical regulators among the potential markers, a transcriptome-wide association screening was performed in the BXD mice panel. Among the candidates, the expression of pleiotrophin (Ptn) was significantly associated with cognitive functions, indicating that Ptn-mediates sleep-loss-induced cognitive impairment. Gene co-expression analysis further revealed the potential mechanism by which Ptn mediates insomnia-induced cognitive impairment via the MAPK signaling pathway; that is, the decreased secretion of Ptn induced by insomnia leads to reduced binding to Ptprz1 on the postsynaptic membrane with the activation of the MAPK pathway via Fos and Nr4a1, further leading to the apoptosis of neurons. In addition, Ptn is genetically trans-regulated in the mouse hippocampus and implicated in neurodegenerative diseases in human genome-wide association studies. Our study provides a novel biomarker for insomnia-induced cognitive impairment and a new strategy for seeking neurological biomarkers by the integration of proteomics and systems genetics.

Identification and experimental validation of mitochondria-related genes biomarkers associated with immune infiltration for sepsis

Background

Sepsis remains a complex condition with incomplete understanding of its pathogenesis. Further research is needed to identify prognostic factors, risk stratification tools, and effective diagnostic and therapeutic targets.

Methods

Three GEO datasets (GSE54514, GSE65682, and GSE95233) were used to explore the potential role of mitochondria-related genes (MiRGs) in sepsis. WGCNA and two machine learning algorithms (RF and LASSO) were used to identify the feature of MiRGs. Consensus clustering was subsequently carried out to determine the molecular subtypes for sepsis. CIBERSORT algorithm was conducted to assess the immune cell infiltration of samples. A nomogram was also established to evaluate the diagnostic ability of feature biomarkers via "rms" package.

Results

Three different expressed MiRGs (DE-MiRGs) were identified as sepsis biomarkers. A significant difference in the immune microenvironment landscape was observed between healthy controls and sepsis patients. Among the DE-MiRGs, NDUFB3 was selected to be a potential therapeutic target and its significant elevated expression level was confirmed in sepsis using in vitro experiments and confocal microscopy, indicating its significant contribution to the mitochondrial quality imbalance in the LPS-simulated sepsis model.

Conclusion

By digging the role of these pivotal genes in immune cell infiltration, we gained a better understanding of the molecular immune mechanism in sepsis and identified potential intervention and treatment strategies.

GPS2-mediated regulation of the adipocyte secretome modulates adipose tissue remodeling at the onset of diet-induced obesity

OBJECTIVE

Dysfunctional, unhealthy expansion of white adipose tissue due to excess dietary intake is a process at the root of obesity and Type 2 Diabetes development. The objective of this study is to contribute to a better understanding of the underlying mechanism(s) regulating the early stages of adipose tissue expansion and adaptation to dietary stress due to an acute, high-fat diet (HFD) challenge, with a focus on the communication between adipocytes and other stromal cells.

METHODS

We profiled the early response to high-fat diet exposure in wildtype and adipocyte-specific GPS2-KO (GPS2-AKO) mice at the cellular, tissue and organismal level. A multi-pronged approach was employed to disentangle the complex cellular interactions dictating tissue remodeling, via single-cell RNA sequencing and FACS profiling of the stromal fraction, and semi-quantitative proteomics of the adipocyte-derived exosomal cargo after 5 weeks of HFD feeding.

RESULTS

Our results indicate that loss of GPS2 in mature adipocytes leads to impaired adaptation to the metabolic stress imposed by HFD feeding. GPS2-AKO mice are significantly more inflamed, insulin resistant, and obese, compared to the WT counterparts. At the cellular level, lack of GPS2 in adipocytes impacts upon other stromal populations, with both the eWAT and scWAT depots exhibiting changes in the immune and non-immune compartments that contribute to an increase in inflammatory and anti-adipogenic cell types. Our studies also revealed that adipocyte to stromal cell communication is facilitated by exosomes, and that transcriptional rewiring of the exosomal cargo is crucial for tissue remodeling. Loss of GPS2 results in increased expression of secreted factors promoting a TGFβ-driven fibrotic microenvironment favoring unhealthy tissue remodeling and expansion.

CONCLUSIONS

Adipocytes serve as an intercellular signaling hub, communicating with the stromal compartment via paracrine signaling. Our study highlights the importance of proper regulation of the 'secretome' released by energetically stressed adipocytes at the onset of obesity. Altered transcriptional regulation of factors secreted via adipocyte-derived exosomes (AdExos), in the absence of GPS2, contributes to the establishment of an anti-adipogenic, pro-fibrotic adipose tissue environment, and to hastened progression towards a metabolically dysfunctional phenotype.

Transcriptome-wide association study reveals cholesterol metabolism gene Lpl is a key regulator of cognitive dysfunction

Cholesterol metabolism in the brain plays a crucial role in normal physiological function, and its aberrations are associated with cognitive dysfunction. The present study aimed to determine which cholesterol-related genes play a vital role in cognitive dysfunction and to dissect its underlying molecular mechanisms using a systems genetics approach in the BXD mice family. We first systematically analyzed the association of expression of 280 hippocampal genes related to cholesterol metabolism with cognition-related traits and identified lipoprotein lipase (Lpl) as a critical regulator. This was further confirmed by phenome-wide association studies that indicate Lpl associated with hippocampus volume residuals and anxiety-related traits. By performing expression quantitative trait locus mapping, we demonstrate that Lpl is strongly cis-regulated in the BXD hippocampus. We also identified ∼3,300 genes significantly (p < 0.05) correlated with the Lpl expression. Those genes are mainly involved in the regulation of neuron-related traits through the MAPK signaling pathway, axon guidance, synaptic vesicle cycle, and NF-kappa B signaling pathway. Furthermore, a protein-protein interaction network analysis identified several direct interactors of Lpl, including Rab3a, Akt1, Igf1, Crp, and Lrp1, which indicates that Lpl involves in the regulation of cognitive dysfunction through Rab3a-mediated synaptic vesicle cycle and Akt1/Igf1/Crp/Lrp1-mediated MAPK signaling pathway. Our findings demonstrate the importance of the Lpl, among the cholesterol-related genes, in regulating cognitive dysfunction and highlighting the potential signaling pathways, which may serve as novel therapeutic targets for the treatment of cognitive dysfunction.

Behavioral phenotypes revealed during reversal learning are linked with novel genetic loci in diversity outbred mice

Impulsive behavior and impulsivity are heritable phenotypes that are strongly associated with risk for substance use disorders. Identifying the neurogenetic mechanisms that influence impulsivity may also reveal novel biological insights into addiction vulnerability. Our past studies using the BXD and Collaborative Cross (CC) recombinant inbred mouse panels have revealed that behavioral indicators of impulsivity measured in a reversal-learning task are heritable and are genetically correlated with aspects of intravenous cocaine self-administration. Genome-wide linkage studies in the BXD panel revealed a quantitative trait locus (QTL) on chromosome 10, but we expect to identify additional QTL by testing in a population with more genetic diversity. To this end, we turned to Diversity Outbred (DO) mice; 392 DO mice (156 males, 236 females) were phenotyped using the same reversal learning test utilized previously. Our primary indicator of impulsive responding, a measure that isolates the relative difficulty mice have with reaching performance criteria under reversal conditions, revealed a genome-wide significant QTL on chromosome 7 (max LOD score = 8.73, genome-wide corrected p<0.05). A measure of premature responding akin to that implemented in the 5-choice serial reaction time task yielded a suggestive QTL on chromosome 17 (max LOD score = 9.14, genome-wide corrected <0.1). Candidate genes were prioritized (2900076A07Rik, Wdr73 and Zscan2) based upon expression QTL data we collected in DO and CC mice and analyses using publicly available gene expression and phenotype databases. These findings may advance understanding of the genetics that drive impulsive behavior and enhance risk for substance use disorders.

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.

Role of gut microbe-derived metabolites in cardiometabolic diseases: Systems based approach

BACKGROUND

The gut microbiome influences host physiology and cardiometabolic diseases by interacting directly with intestinal cells or by producing molecules that enter the host circulation. Given the large number of microbial species present in the gut and the numerous factors that influence gut bacterial composition, it has been challenging to understand the underlying biological mechanisms that modulate risk of cardiometabolic disease.

SCOPE OF THE REVIEW

Here we discuss a systems-based approach that involves simultaneously examining individuals in populations for gut microbiome composition, molecular traits using "omics" technologies, such as circulating metabolites quantified by mass spectrometry, and clinical traits. We summarize findings from landmark studies using this approach and discuss future applications.

MAJOR CONCLUSIONS

Population-based integrative approaches have identified a large number of microbe-derived or microbe-modified metabolites that are associated with cardiometabolic traits. The knowledge gained from these studies provide new opportunities for understanding the mechanisms involved in gut microbiome-host interactions and may have potentially important implications for developing novel therapeutic approaches.

Genetic Variants Associated With Resilience in Human and Animal Studies

Resilience is broadly defined as the ability to maintain or regain functioning in the face of adversity and is influenced by both environmental and genetic factors. The identification of specific genetic factors and their biological pathways underpinning resilient functioning can help in the identification of common key factors, but heterogeneities in the operationalisation of resilience have hampered advances. We conducted a systematic review of genetic variants associated with resilience to enable the identification of general resilience mechanisms. We adopted broad inclusion criteria for the definition of resilience to capture both human and animal model studies, which use a wide range of resilience definitions and measure very different outcomes. Analyzing 158 studies, we found 71 candidate genes associated with resilience. OPRM1 (Opioid receptor mu 1), NPY (neuropeptide Y), CACNA1C (calcium voltage-gated channel subunit alpha1 C), DCC (deleted in colorectal carcinoma), and FKBP5 (FKBP prolyl isomerase 5) had both animal and human variants associated with resilience, supporting the idea of shared biological pathways. Further, for OPRM1, OXTR (oxytocin receptor), CRHR1 (corticotropin-releasing hormone receptor 1), COMT (catechol-O-methyltransferase), BDNF (brain-derived neurotrophic factor), APOE (apolipoprotein E), and SLC6A4 (solute carrier family 6 member 4), the same allele was associated with resilience across divergent resilience definitions, which suggests these genes may therefore provide a starting point for further research examining commonality in resilience pathways.

Extreme phenotypic diversity in operant response to intravenous cocaine or saline infusion in the hybrid mouse diversity panel

Cocaine self-administration is a complexly determined trait, with a substantial proportion of individual differences being determined by genetic variation. However, the relevant genetic variants that drive heritable differences in cocaine use remain undiscovered. Cocaine intravenous self-administration (IVSA) procedures in laboratory animals provide opportunities to prospectively investigate neurogenetic influences on the acquisition of voluntary cocaine use. Here, we provide information on cocaine (or saline-as a control) IVSA in 84 members of the hybrid mouse diversity panel (HMDP), an array of genetically distinct classical or recombinant inbred strains. We found cocaine IVSA to be substantially heritable in this population, with strain-level intake ranging for near 0 to >25 mg/kg/session. Though saline IVSA was also found to be heritable, a modest genetic correlation between cocaine and saline IVSA indicates that operant responding for the cocaine reinforcer was influenced, at least in part, by unique genetic variants. Genome-wide association studies (GWAS) of infusions earned in cocaine and saline groups revealed significant quantitative trait loci (QTL) on Chromosomes 3 and 14 for cocaine, but not saline, IVSA. Positional candidates were further prioritized through use of bulk RNA sequencing data that revealed genes with cis-eQTL and genetic correlation to number of infusions. Additionally, these data identify reference strains with extreme cocaine IVSA phenotypes, revealing them as polygenic models of risk and resilience to cocaine reinforcement. This work is part of an ongoing effort to characterize genetic variation that moderates cocaine IVSA that may, in turn, provide a more comprehensive understanding of cocaine risk genetics and neurobiology.

Translational approaches to understanding resilience to Alzheimer's disease

Individuals who maintain cognitive function despite high levels of Alzheimer's disease (AD)-associated pathology are said to be 'resilient' to AD. Identifying mechanisms underlying resilience represents an exciting therapeutic opportunity. Human studies have identified a number of molecular and genetic factors associated with resilience, but the complexity of these cohorts prohibits a complete understanding of which factors are causal or simply correlated with resilience. Genetically and phenotypically diverse mouse models of AD provide new and translationally relevant opportunities to identify and prioritize new resilience mechanisms for further cross-species investigation. This review will discuss insights into resilience gained from both human and animal studies and highlight future approaches that may help translate these insights into therapeutics designed to prevent or delay AD-related dementia.

Maturation and application of phenome-wide association studies

In the past 10 years since its introduction, phenome-wide association studies (PheWAS) have uncovered novel genotype-phenotype relationships. Along the way, PheWAS have evolved in many aspects as a study design with the expanded availability of large data repositories with genome-wide data linked to detailed phenotypic data. Advancement in methods, including algorithms, software, and publicly available integrated resources, makes it feasible to more fully realize the potential of PheWAS, overcoming the previous computational and analytical limitations. We review here the most recent improvements and notable applications of PheWAS since the second half of the decade from its inception. We also note the challenges that remain embedded along the entire PheWAS analytical pipeline that necessitate further development of tools and resources to further advance the understanding of the complex genetic architecture underlying human diseases and traits.

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.

•

Dox extends lifespan of N2 and CB4856 via shared and strain-specific mechanisms

•

Dox controls mitochondria, defense responses, and lipid metabolism in both strains

•

Dox-mediated longevity requires acs-2/20 in N2 and fat-7/6 in CB4856 worms

Inhibiting serotonin signaling through HTR2B in visceral adipose tissue improves obesity-related insulin resistance

Insulin resistance is a cornerstone of obesity-related complications such as type 2 diabetes, metabolic syndrome, and nonalcoholic fatty liver disease. A high rate of lipolysis is known to be associated with insulin resistance, and inhibiting adipose tissue lipolysis improves obesity-related insulin resistance. Here, we demonstrate that inhibition of serotonin (5-hydroxytryptamine [5-HT]) signaling through serotonin receptor 2B (HTR2B) in adipose tissues ameliorates insulin resistance by reducing lipolysis in visceral adipocytes. Chronic high-fat diet (HFD) feeding increased Htr2b expression in epididymal white adipose tissue, resulting in increased HTR2B signaling in visceral white adipose tissue. Moreover, HTR2B expression in white adipose tissue was increased in obese humans and positively correlated with metabolic parameters. We further found that adipocyte-specific Htr2b-knockout mice are resistant to HFD-induced insulin resistance, visceral adipose tissue inflammation, and hepatic steatosis. Enhanced 5-HT signaling through HTR2B directly activated lipolysis through phosphorylation of hormone-sensitive lipase in visceral adipocytes. Moreover, treatment with a selective HTR2B antagonist attenuated HFD-induced insulin resistance, visceral adipose tissue inflammation, and hepatic steatosis. Thus, adipose HTR2B signaling could be a potential therapeutic target for treatment of obesity-related insulin resistance.

Combining power of different methods to detect associations in large data sets

Exploring the relationship between factors of interest is a fundamental step for further analysis on various scientific problems such as understanding the genetic mechanism underlying specific disease, brain functional connectivity analysis. There are many methods proposed for association analysis and each has its own advantages, but none of them is suitable for all kinds of situations. This brings difficulties and confusions to practitioner on which one to use when facing a real problem. In this paper, we propose to combine power of different methods to detect associations in large data sets. It goes as combining the weaker to be stronger. Numerical results from simulation study and real data applications show that our new framework is powerful. Importantly, the framework can also be applied to other problems. Availability: The R script is available at https://jiangdata.github.io/resources/DM.zip.

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

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

Do gene expression findings from mouse models of cocaine use recapitulate human cocaine use disorder in reward circuitry?

Animal models of drug use have investigated possible mechanisms governing human substance use traits for over 100 years. Most cross-species research on drug use/addiction examines behavioral overlap, but studies assessing neuromolecular (e.g. RNA) correspondence are lacking. Our study utilized transcriptome-wide data from the hippocampus and ventral tegmental area (VTA)/midbrain from a total of 35 human males with cocaine use disorder/controls and 49 male C57BL/6J cocaine/saline administering/exposed mice. We hypothesized differential expressed genes and systems of co-expressed genes (gene networks) would show appreciable overlap across mouse cocaine self-administration and human cocaine use disorder. We found modest, but significant relationships between differentially expressed genes associated with cocaine self-administration (short access) and cocaine use disorder within reward circuitry. Differentially expressed genes underlying models of acute cocaine exposure (cocaine), context re-exposure and cocaine + context re-exposure were not consistently associated with human CUD across brain regions. Investigating systems of co-expressed genes, we found several validated gene networks with weak to moderate conservation between cocaine/saline self-administering mice and disordered cocaine users/controls. The most conserved hippocampal and VTA gene networks demonstrated substantial overlap (2029 common genes) and included both novel and previously implicated targets for cocaine use/addiction. Lastly, we conducted (expression-based) phenome-wide association studies of the nine common hub genes across conserved gene networks. Common hub genes were associated with dopamine/serotonin function, cocaine self-administration and other relevant mouse traits. Overall, our study pinpointed and characterized conserved brain-related RNA patterns across mouse cocaine self-administration and human cocaine use disorder. We offer recommendations for future research and add to the dialogue surrounding pre-clinical animal research for human disease.

Mouse Systems Genetics as a Prelude to Precision Medicine

Mouse models have been instrumental in understanding human disease biology and proposing possible new treatments. The precise control of the environment and genetic composition of mice allows more rigorous observations, but limits the generalizability and translatability of the results into human applications. In the era of precision medicine, strategies using mouse models have to be revisited to effectively emulate human populations. Systems genetics is one promising paradigm that may promote the transition to novel precision medicine strategies. Here, we review the state-of-the-art resources and discuss how mouse systems genetics helps to understand human diseases and to advance the development of precision medicine, with an emphasis on the existing resources and strategies.

Unknown to Known: Advancing Knowledge of Coral Gene Function

Given the catastrophic changes befalling coral reefs, understanding coral gene function is essential to advance reef conservation. This has proved challenging due to the paucity of genomic data and genetic tools available for corals. Recently, CRISPR/Cas9 gene editing was applied to these species; however, a major bottleneck is the identification and prioritization of candidate genes for manipulation. This issue is exacerbated by the many unknown ('dark') coral genes that may play key roles in the stress response. We review the use of gene coexpression networks that incorporate both known and unknown genes to identify targets for reverse genetic analysis. This approach also provides a framework for the annotation of dark genes in established interaction networks to improve our fundamental knowledge of coral gene function.

Genetic Dissection of Femoral and Tibial Microarchitecture

Our understanding of the genetic control of bone strength has relied mainly on estimates of bone mineral density. Here we have mapped genetic factors that influence femoral and tibial microarchitecture using high-resolution x-ray computed tomography (8-μm isotropic voxels) across a family of 61 BXD strains of mice, roughly 10 isogenic cases per strain and balanced by sex. We computed heritabilities for 25 cortical and trabecular traits. Males and females have well-matched heritabilities, ranging from 0.25 to 0.75. We mapped 16 genetic loci most of which were detected only in females. There is also a bias in favor of loci that control cortical rather than trabecular bone. To evaluate candidate genes, we combined well-established gene ontologies with bone transcriptome data to compute bone-enrichment scores for all protein-coding genes. We aligned candidates with those of human genome-wide association studies. A subset of 50 strong candidates fell into three categories: (1) experimentally validated genes already known to modulate bone function (Adamts4, Ddr2, Darc, Adam12, Fkbp10, E2f6, Adam17, Grem2, Ifi204); (2) candidates without any experimentally validated function in bone (eg, Greb1, Ifi202b), but linked to skeletal phenotypes in human cohorts; and (3) candidates that have high bone-enrichment scores, but for which there is not yet any functional link to bone biology or skeletal system disease (including Ifi202b, Ly9, Ifi205, Mgmt, F2rl1, Iqgap2). Our results highlight contrasting genetic architecture between sexes and among major bone compartments. The alignment of murine and human data facilitates function analysis and should prove of value for preclinical testing of molecular control of bone structure. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

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.

Systems genetics applications in metabolism research

The common forms of metabolic diseases are highly complex, involving hundreds of genes, environmental and lifestyle factors, age-related changes, sex differences and gut-microbiome interactions. Systems genetics is a population-based approach to address this complexity. In contrast to commonly used 'reductionist' approaches, such as gain or loss of function, that examine one element at a time, systems genetics uses high-throughput 'omics' technologies to quantitatively assess the many molecular differences among individuals in a population and then to relate these to physiologic functions or disease states. Unlike genome-wide association studies, systems genetics seeks to go beyond the identification of disease-causing genes to understand higher-order interactions at the molecular level. The purpose of this review is to introduce the systems genetics applications in the areas of metabolic and cardiovascular disease. Here, we explain how large clinical and omics-level data and databases from both human and animal populations are available to help researchers place genes in the context of pathways and networks and formulate hypotheses that can then be experimentally examined. We provide lists of such databases and examples of the integration of reductionist and systems genetics data. Among the important applications emerging is the development of improved nutritional and pharmacological strategies to address the rise of metabolic diseases.

Diet modulates cecum bacterial diversity and physiological phenotypes across the BXD mouse genetic reference population

The BXD family has become one of the preeminent genetic reference populations to understand the genetic and environmental control of phenotypic variation. Here we evaluate the responses to different levels of fat in the diet using both chow diet (CD, 13–18% fat) and a high-fat diet (HFD, 45–60% fat). We studied cohorts of BXD strains, both inbred parents C57BL/6J and DBA/2J (commonly known as B6 and D2, respectively), as well as B6D2 and D2B6 reciprocal F1 hybrids. The comparative impact of genetic and dietary factors was analyzed by profiling a range of phenotypes, most prominently their cecum bacterial composition. The parents of the BXDs and F1 hybrids express limited differences in terms of weight and body fat gain on CD. In contrast, the strain differences on HFD are substantial for percent body fat, with DBA/2J accumulating 12.5% more fat than C57BL/6J (P < 0.0001). The F1 hybrids born to DBA/2J dams (D2B6F1) have 10.6% more body fat (P < 0.001) than those born to C57BL/6J dams. Sequence analysis of the cecum microbiota reveals important differences in bacterial composition among BXD family members with a substantial shift in composition caused by HFD. Relative to CD, the HFD induces a decline in diversity at the phylum level with a substantial increase in Firmicutes (+13.8%) and a reduction in Actinobacteria (-7.9%). In the majority of BXD strains, the HFD also increases cecal sIgA (P < 0.0001)—an important component of the adaptive immunity response against microbial pathogens. Host genetics modulates variation in cecum bacterial composition at the genus level in CD, with significant quantitative trait loci (QTLs) for Oscillibacter mapped to Chr 3 (18.7–19.2 Mb, LRS = 21.4) and for Bifidobacterium mapped to Chr 6 (89.21–89.37 Mb, LRS = 19.4). Introduction of HFD served as an environmental suppressor of these QTLs due to a reduction in the contribution of both genera (P < 0.001). Relations among liver metabolites and cecum bacterial composition were predominant in CD cohort, but these correlations do not persist following the shift to HFD. Overall, these findings demonstrate the important impact of environmental/dietary manipulation on the relationships between host genetics, gastrointestinal bacterial composition, immunological parameters, and metabolites—knowledge that will help in the understanding of the causal sources of metabolic disorders.

Cross-species functional modules link proteostasis to human normal aging

The evolutionarily conserved nature of the few well-known anti-aging interventions that affect lifespan, such as caloric restriction, suggests that aging-related research in model organisms is directly relevant to human aging. Since human lifespan is a complex trait, a systems-level approach will contribute to a more comprehensive understanding of the underlying aging landscape. Here, we integrate evolutionary and functional information of normal aging across human and model organisms at three levels: gene-level, process-level, and network-level. We identify evolutionarily conserved modules of normal aging across diverse taxa, and notably show proteostasis to be conserved in normal aging. Additionally, we find that mechanisms related to protein quality control network are enriched for genes harboring genetic variants associated with 22 age-related human traits and associated to caloric restriction. These results demonstrate that a systems-level approach, combined with evolutionary conservation, allows the detection of candidate aging genes and pathways relevant to human normal aging.

Discovery of early life stress interacting and sex-specific quantitative trait loci impacting cocaine responsiveness

BACKGROUND AND PURPOSE

Addiction vulnerability involves complex gene X environment interactions leading to a pathological response to drugs. Identification of the genes involved in these interactions is an important step in understanding the underlying neurobiology and rarely have such analyses examined sex-specific influences. To dissect this interaction, we examined the impact of prenatal stress (PNS) on cocaine responsiveness in male and female mice of the BXD recombinant inbred panel.

EXPERIMENTAL APPROACH

BXD strains were subjected to timed mating and assigned to PNS or control groups. PNS dams were subjected to restraint stress (1-hr restraint, three times daily) starting between embryonic day (E) 11 and 14 and continued until parturition. Adult male and female, control and PNS offspring were tested for locomotor response to initial and repeated cocaine injections (sensitization) as well as cocaine-induced conditioned place preference (CPP).

KEY RESULTS

Strain, PNS, and sex interacted to modulate initial and sensitized cocaine-induced locomotion, as well as CPP. Moreover, a quantitative trait locus (QTL) interacting with PNS regulating initial locomotor response to cocaine (chromosome X, 37.91 to 50.95 Mb) was identified. Also PNS-independent, female-specific QTLs regulating CPP (chromosome 11, 65.50 to 81.31 Mb) and sensitized cocaine-induced locomotion (chromosome 16, 95.79 to 98.32 Mb) were identified. Publicly available mRNA expression data were utilized to identify cis-eQTL and transcript covariation with the behavioural phenotype to prioritize candidate genes; including Aifm1.

CONCLUSIONS AND IMPLICATIONS

These QTL encompass genes that may moderate genetic susceptibility to PNS and interact with sex to determine adult responsiveness to cocaine and addiction vulnerability.

LINKED ARTICLES

This article is part of a themed section on The Importance of Sex Differences in Pharmacology Research. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.21/issuetoc.

Tissue-Specific Trans Regulation of the Mouse Epigenome

The epigenetic landscape varies greatly among cell types. Although a variety of writers, readers, and erasers of epigenetic features are known, we have little information about the underlying regulatory systems controlling the establishment and maintenance of these features. Here, we have explored how natural genetic variation affects the epigenome in mice. Studying levels of H3K4me3, a histone modification at sites such as promoters, enhancers, and recombination hotspots, we found tissue-specific trans-regulation of H3K4me3 levels in four highly diverse cell types: male germ cells, embryonic stem cells, hepatocytes, and cardiomyocytes. To identify the genetic loci involved, we measured H3K4me3 levels in male germ cells in a mapping population of 59 BXD recombinant inbred lines. We found extensive trans-regulation of H3K4me3 peaks, including six major histone quantitative trait loci (QTL). These chromatin regulatory loci act dominantly to suppress H3K4me3, which at hotspots reduces the likelihood of subsequent DNA double-strand breaks. QTL locations do not correspond with genes encoding enzymes known to metabolize chromatin features. Instead their locations match clusters of zinc finger genes, making these possible candidates that explain the dominant suppression of H3K4me3. Collectively, these data describe an extensive, set of chromatin regulatory loci that control the epigenetic landscape.

The 40S ribosomal protein uS5 (RPS2) assembles into an extraribosomal complex with human ZNF277 that competes with the PRMT3-uS5 interaction

Ribosomal (r)-proteins are generally viewed as ubiquitous, constitutive proteins that simply function to maintain ribosome integrity. However, findings in the past decade have led to the idea that r-proteins have evolved specialized functions beyond the ribosome. For example, the 40S ribosomal protein uS5 (RPS2) is known to form an extraribosomal complex with the protein arginine methyltransferase PRMT3 that is conserved from fission yeast to humans. However, the full scope of uS5's extraribosomal functions, including whether uS5 interacts with any other proteins, is not known. In this study, we identify the conserved zinc finger protein 277 (ZNF277) as a new uS5-associated protein by using quantitative proteomics approaches in human cells. As previously shown for PRMT3, we found that ZNF277 uses a C2H2-type zinc finger domain to recognize uS5. Analysis of protein-protein interactions in living cells indicated that the ZNF277-uS5 complex is found in the cytoplasm and the nucleolus. Furthermore, we show that ZNF277 and PRMT3 compete for uS5 binding, because overexpression of PRMT3 inhibited the formation of the ZNF277-uS5 complex, whereas depletion of cellular ZNF277 resulted in increased levels of uS5-PRMT3. Notably, our results reveal that ZNF277 recognizes nascent uS5 in the course of mRNA translation, suggesting cotranslational assembly of the ZNF277-uS5 complex. Our findings thus unveil an intricate network of evolutionarily conserved protein-protein interactions involving extraribosomal uS5, suggesting a key role for uS5 beyond the ribosome.