Systems genetics applications in metabolism research

Combinatorial in vivo genome editing identifies widespread epistasis during lung tumorigenesis

Lung adenocarcinoma, the most common subtype of lung cancer, is genomically complex, with tumors containing tens to hundreds of non-synonymous mutations. However, little is understood about how genes interact with each other to enable tumorigenesis in vivo , largely due to a lack of methods for investigating genetic interactions in a high-throughput and multiplexed manner. Here, we employed a novel platform to generate tumors with all pairwise inactivation of ten tumor suppressor genes within an autochthonous mouse model of oncogenic KRAS-driven lung cancer. By quantifying the fitness of tumors with every single and double mutant genotype, we show that most tumor suppressor genetic interactions exhibited negative epistasis, with diminishing returns on tumor fitness. In contrast, Apc inactivation showed positive epistasis with the inactivation of several other genes, including dramatically synergistic effects on tumor fitness in combination with Lkb1 or Nf1 inactivation. This approach has the potential to expand the scope of genetic interactions that may be functionally characterized in vivo , which could lead to a better understanding of how complex tumor genotypes impact each step of carcinogenesis.

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.

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.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

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.

When diet meets genetics

Gene expression profiling of a diverse mouse population helps to decipher how a fat-rich diet contributes to inflammatory bowel disease.

An atlas of genome-wide gene expression and metabolite associations and possible mediation effects towards body mass index

Investigating the cross talk of different omics layers is crucial to understand molecular pathomechanisms of metabolic diseases like obesity. Here, we present a large-scale association meta-analysis of genome-wide whole blood and peripheral blood mononuclear cell (PBMC) gene expressions profiled with Illumina HT12v4 microarrays and metabolite measurements from dried blood spots (DBS) characterized by targeted liquid chromatography tandem mass spectrometry (LC-MS/MS) in three large German cohort studies with up to 7706 samples. We found 37,295 associations comprising 72 amino acids (AA) and acylcarnitine (AC) metabolites (including ratios) and 8579 transcripts. We applied this catalogue of associations to investigate the impact of associating transcript-metabolite pairs on body mass index (BMI) as an example metabolic trait. This is achieved by conducting a comprehensive mediation analysis considering metabolites as mediators of gene expression effects and vice versa. We discovered large mediation networks comprising 27,023 potential mediation effects within 20,507 transcript-metabolite pairs. Resulting networks of highly connected (hub) transcripts and metabolites were leveraged to gain mechanistic insights into metabolic signaling pathways. In conclusion, here, we present the largest available multi-omics integration of genome-wide transcriptome data and metabolite data of amino acid and fatty acid metabolism and further leverage these findings to characterize potential mediation effects towards BMI proposing candidate mechanisms of obesity and related metabolic diseases. KEY MESSAGES: Thousands of associations of 72 amino acid and acylcarnitine metabolites and 8579 genes expand the knowledge of metabolome-transcriptome associations. A mediation analysis of effects on body mass index revealed large mediation networks of thousands of obesity-related gene-metabolite pairs. Highly connected, potentially mediating hub genes and metabolites enabled insight into obesity and related metabolic disease pathomechanisms.

Leveraging genetic diversity to identify small molecules that reverse mouse skeletal muscle insulin resistance

Systems genetics has begun to tackle the complexity of insulin resistance by capitalising on computational advances to study high-diversity populations. 'Diversity Outbred in Australia (DOz)' is a population of genetically unique mice with profound metabolic heterogeneity. We leveraged this variance to explore skeletal muscle's contribution to whole-body insulin action through metabolic phenotyping and skeletal muscle proteomics of 215 DOz mice. Linear modelling identified 553 proteins that associated with whole-body insulin sensitivity (Matsuda Index) including regulators of endocytosis and muscle proteostasis. To enrich for causality, we refined this network by focusing on negatively associated, genetically regulated proteins, resulting in a 76-protein fingerprint of insulin resistance. We sought to perturb this network and restore insulin action with small molecules by integrating the Broad Institute Connectivity Map platform and in vitro assays of insulin action using the Prestwick chemical library. These complementary approaches identified the antibiotic thiostrepton as an insulin resistance reversal agent. Subsequent validation in ex vivo insulin-resistant mouse muscle and palmitate-induced insulin-resistant myotubes demonstrated potent insulin action restoration, potentially via upregulation of glycolysis. This work demonstrates the value of a drug-centric framework to validate systems-level analysis by identifying potential therapeutics for insulin resistance.

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.

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.

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.

Atherosclerosis: Recent developments

Atherosclerosis is an inflammatory disease of the large arteries that is the major cause of cardiovascular disease (CVD) and stroke. Here, we review the current understanding of the molecular, cellular, genetic, and environmental contributions to atherosclerosis, from both individual pathway and systems perspectives. We place an emphasis on recent developments, some of which have yielded unexpected biology, including previously unknown heterogeneity of inflammatory and smooth muscle cells in atherosclerotic lesions, roles for senescence and clonal hematopoiesis, and links to the gut microbiome.

Genetic variation of putative myokine signaling is dominated by biological sex and sex hormones

Skeletal muscle plays an integral role in coordinating physiological homeostasis, where signaling to other tissues via myokines allows for coordination of complex processes. Here, we aimed to leverage natural genetic correlation structure of gene expression both within and across tissues to understand how muscle interacts with metabolic tissues. Specifically, we performed a survey of genetic correlations focused on myokine gene regulation, muscle cell composition, cross-tissue signaling, and interactions with genetic sex in humans. While expression levels of a majority of myokines and cell proportions within skeletal muscle showed little relative differences between males and females, nearly all significant cross-tissue enrichments operated in a sex-specific or hormone-dependent fashion; in particular, with estradiol. These sex- and hormone-specific effects were consistent across key metabolic tissues: liver, pancreas, hypothalamus, intestine, heart, visceral, and subcutaneous adipose tissue. To characterize the role of estradiol receptor signaling on myokine expression, we generated male and female mice which lack estrogen receptor α specifically in skeletal muscle (MERKO) and integrated with human data. These analyses highlighted potential mechanisms of sex-dependent myokine signaling conserved between species, such as myostatin enriched for divergent substrate utilization pathways between sexes. Several other putative sex-dependent mechanisms of myokine signaling were uncovered, such as muscle-derived tumor necrosis factor alpha (TNFA) enriched for stronger inflammatory signaling in females compared to males and GPX3 as a male-specific link between glycolytic fiber abundance and hepatic inflammation. Collectively, we provide a population genetics framework for inferring muscle signaling to metabolic tissues in humans. We further highlight sex and estradiol receptor signaling as critical variables when assaying myokine functions and how changes in cell composition are predicted to impact other metabolic organs.

Cancer Metabolism: The Role of Immune Cells Epigenetic Alteration in Tumorigenesis, Progression, and Metastasis of Glioma

Glioma is a type of brain and spinal cord tumor that begins in glial cells that support the nervous system neurons functions. Age, radiation exposure, and family background of glioma constitute are risk factors of glioma initiation. Gliomas are categorized on a scale of four grades according to their growth rate. Grades one and two grow slowly, while grades three and four grow faster. Glioblastoma is a grade four gliomas and the deadliest due to its aggressive nature (accelerated proliferation, invasion, and migration). As such, multiple therapeutic approaches are required to improve treatment outcomes. Recently, studies have implicated the significant roles of immune cells in tumorigenesis and the progression of glioma. The energy demands of gliomas alter their microenvironment quality, thereby inducing heterogeneity and plasticity change of stromal and immune cells via the PI3K/AKT/mTOR pathway, which ultimately results in epigenetic modifications that facilitates tumor growth. PI3K is utilized by many intracellular signaling pathways ensuring the proper functioning of the cell. The activation of PI3K/AKT/mTOR regulates the plasma membrane activities, contributing to the phosphorylation reaction necessary for transcription factors activities and oncogenes hyperactivation. The pleiotropic nature of PI3K/AKT/mTOR makes its activity unpredictable during altered cellular functions. Modification of cancer cell microenvironment affects many cell types, including immune cells that are the frontline cells involved in inflammatory cascades caused by cancer cells via high cytokines synthesis. Typically, the evasion of immunosurveillance by gliomas and their resistance to treatment has been attributed to epigenetic reprogramming of immune cells in the tumor microenvironment, which results from cancer metabolism. Hence, it is speculative that impeding cancer metabolism and/or circumventing the epigenetic alteration of immune cell functions in the tumor microenvironment might enhance treatment outcomes. Herein, from an oncological and immunological perspective, this review discusses the underlying pathomechanism of cell-cell interactions enhancing glioma initiation and metabolism activation and tumor microenvironment changes that affect epigenetic modifications in immune cells. Finally, prospects for therapeutic intervention were highlighted.

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.

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.

Sex-specific genetic regulation of adipose mitochondria and metabolic syndrome by Ndufv2

We have previously suggested a central role for mitochondria in the observed sex differences in metabolic traits. However, the mechanisms by which sex differences affect adipose mitochondrial function and metabolic syndrome are unclear. Here we show that in both mice and humans, adipose mitochondrial functions are elevated in females and are strongly associated with adiposity, insulin resistance and plasma lipids. Using a panel of diverse inbred strains of mice, we identify a genetic locus on mouse chromosome 17 that controls mitochondrial mass and function in adipose tissue in a sex- and tissue-specific manner. This locus contains Ndufv2 and regulates the expression of at least 89 mitochondrial genes in females, including oxidative phosphorylation genes and those related to mitochondrial DNA content. Overexpression studies indicate that Ndufv2 mediates these effects by regulating supercomplex assembly and elevating mitochondrial reactive oxygen species production, which generates a signal that increases mitochondrial biogenesis.

Genetics of coronary artery disease in the post-GWAS era

During the past decade, genome-wide association studies (GWAS) have transformed our understanding of many heritable traits. Three recent large-scale GWAS meta-analyses now further markedly expand the knowledge on coronary artery disease (CAD) genetics in doubling the number of loci with genome-wide significant signals. Here, we review the unprecedented discoveries of CAD GWAS on low-frequency variants, underrepresented populations, sex differences and integrated polygenic risk. We present the milestones of CAD GWAS and post-GWAS studies from 2007 to 2021, and the trend in identification of variants with smaller odds ratio by year due to the increasing sample size. We compile the 321 CAD loci discovered thus far and classify candidate genes as well as distinct functional pathways on the road to indepth biological investigation and identification of novel treatment targets. We draw attention to systems genetics in integrating these loci into gene regulatory networks within and across tissues. We review the traits, biomarkers and diseases scrutinized by Mendelian randomization studies for CAD. Finally, we discuss the potentials and concerns of polygenic scores in predicting CAD risk in patient care as well as future directions of GWAS and post-GWAS studies in the field of precision medicine.

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.

CoffeeProt: an online tool for correlation and functional enrichment of systems genetics data

The integration of genomics, transcriptomics, proteomics and phenotypic traits across genetically diverse populations is a powerful approach to discover novel biological regulators. The increasing volume of complex data require new and easy-to-use tools accessible to a variety of scientists for the discovery and visualization of functionally relevant associations. To meet this requirement, we developed CoffeeProt, an open-source tool that analyses genetic variants associated to protein networks, other omics datatypes and phenotypic traits. CoffeeProt uses transcriptomics or proteomics data to perform correlation network analyses and annotates results with protein-protein interactions, subcellular localisations and drug associations. It then integrates genetic variants associated with gene expression (eQTLs) or protein abundance (pQTLs) and includes predictions of the potential consequences of variants on gene function. Finally, genetic variants are co-mapped to molecular or phenotypic traits either provided by the user or retrieved directly from publicly available GWAS results. We demonstrate its utility with the analysis of mouse and human population data enabling the rapid identification of genetic variants associated with druggable proteins and clinical traits. We expect that CoffeeProt will serve the systems genetics and basic science research communities, leading to the discovery of novel biologically relevant associations. CoffeeProt is available at www.coffeeprot.com.

Genetic regulation of liver lipids in a mouse model of insulin resistance and hepatic steatosis

To elucidate the contributions of specific lipid species to metabolic traits, we integrated global hepatic lipid data with other omics measures and genetic data from a cohort of about 100 diverse inbred strains of mice fed a high-fat/high-sucrose diet for 8 weeks. Association mapping, correlation, structure analyses, and network modeling revealed pathways and genes underlying these interactions. In particular, our studies lead to the identification of Ifi203 and Map2k6 as regulators of hepatic phosphatidylcholine homeostasis and triacylglycerol accumulation, respectively. Our analyses highlight mechanisms for how genetic variation in hepatic lipidome can be linked to physiological and molecular phenotypes, such as microbiota composition.

Exploring the Multi-Tissue Crosstalk Relevant to Insulin Resistance Through Network-Based Analysis

Insulin resistance (IR) is a precursor event that occurs in multiple organs and underpins many metabolic disorders. However, due to the lack of effective means to systematically explore and interpret disease-related tissue crosstalk, the tissue communication mechanism in pathogenesis of IR has not been elucidated yet. To solve this issue, we profiled all proteins in white adipose tissue (WAT), liver, and skeletal muscle of a high fat diet induced IR mouse model via proteomics. A network-based approach was proposed to explore IR related tissue communications. The cross-tissue interface was constructed, in which the inter-tissue connections and also their up and downstream processes were particularly inspected. By functional quantification, liver was recognized as the only organ that can output abnormal carbohydrate metabolic signals, clearly highlighting its central role in regulation of glucose homeostasis. Especially, the CD36-PPAR axis in liver and WAT was identified and verified as a potential bridge that links cross-tissue signals with intracellular metabolism, thereby promoting the progression of IR through a PCK1-mediated lipotoxicity mechanism. The cross-tissue mechanism unraveled in this study not only provides novel insights into the pathogenesis of IR, but also is conducive to development of precision therapies against various IR associated diseases. With further improvement, our network-based cross-tissue analytic method would facilitate other disease-related tissue crosstalk study in the near future.