Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 (Abhd2) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 (Abhd2), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2. The Abhd2KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Mining cholesterol genes from thousands of mouse livers identifies aldolase C as a regulator of cholesterol biosynthesis.

The availability of genome-wide transcriptomic and proteomic datasets is ever-increasing and often not used beyond initial publication. Here, we applied module-based coexpression network analysis to a comprehensive catalog of 35 mouse genome-wide liver expression datasets (encompassing more than 3800 mice) with the goal of identifying and validating unknown genes involved in cholesterol metabolism. From these 35 datasets, we identified a conserved module of genes enriched with cholesterol biosynthetic genes. Using a systematic approach across the 35 datasets, we identified three genes (Rdh11, Echdc1, and Aldoc) with no known role in cholesterol metabolism. We then performed functional validation studies and show that each gene is capable of regulating cholesterol metabolism. For the glycolytic gene, Aldoc, we demonstrate that it contributes to de novo cholesterol biosynthesis and regulates cholesterol and triglyceride levels in mice. As Aldoc is located within a genome-wide significant genome-wide association studies locus for human plasma cholesterol levels, our studies establish Aldoc as a causal gene within this locus. Through our work, we develop a framework for leveraging mouse genome-wide liver datasets for identifying and validating genes involved in cholesterol metabolism.

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.

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.

Cross-species data integration to prioritize causal genes in lipid metabolism.

PURPOSE OF REVIEW: More than one hundred loci have been identified from human genome-wide association studies (GWAS) for blood lipids. Despite the success of GWAS in identifying loci, subsequent prioritization of causal genes related to these loci remains a challenge. To address this challenge, recent work suggests that candidate causal genes within loci can be prioritized through cross-species integration using genome-wide data from the mouse. RECENT FINDINGS: Mouse model systems provide unparalleled access to primary tissues, like the liver, that are not readily available for human studies. Given the key role the liver plays in controlling blood lipid levels and the wealth of liver genome-wide transcript and protein data available in the mouse, these data can be leveraged. Using coexpression network analysis approaches with mouse genome-wide data, coupled with cross-species analysis of human lipid GWAS, causal genes within lipid loci can be prioritized. Prioritization through both mouse and human along with biochemical validation provide a systematic and valuable method to discover lipid metabolism genes. SUMMARY: The prioritization of causal lipid genes within GWAS loci is a challenging process requiring a multidisciplinary approach. Integration of data types across species, such as the mouse, can aid in causal gene prioritization.

Genetic Regulation of Adipose Gene Expression and Cardio-Metabolic Traits.

Subcutaneous adipose tissue stores excess lipids and maintains energy balance. We performed expression quantitative trait locus (eQTL) analyses by using abdominal subcutaneous adipose tissue of 770 extensively phenotyped participants of the METSIM study. We identified cis-eQTLs for 12,400 genes at a 1% false-discovery rate. Among an approximately 680 known genome-wide association study (GWAS) loci for cardio-metabolic traits, we identified 140 coincident cis-eQTLs at 109 GWAS loci, including 93 eQTLs not previously described. At 49 of these 140 eQTLs, gene expression was nominally associated (p < 0.05) with levels of the GWAS trait. The size of our dataset enabled identification of five loci associated (p < 5 × 10-8) with at least five genes located >5 Mb away. These trans-eQTL signals confirmed and extended the previously reported KLF14-mediated network to 55 target genes, validated the CIITA regulation of class II MHC genes, and identified ZNF800 as a candidate master regulator. Finally, we observed similar expression-clinical trait correlations of genes associated with GWAS loci in both humans and a panel of genetically diverse mice. These results provide candidate genes for further investigation of their potential roles in adipose biology and in regulating cardio-metabolic traits.

Regulator of Calcineurin 1 helps coordinate whole-body metabolism and thermogenesis.

Increasing non-shivering thermogenesis (NST), which expends calories as heat rather than storing them as fat, is championed as an effective way to combat obesity and metabolic disease. Innate mechanisms constraining the capacity for NST present a fundamental limitation to this approach, yet are not well understood. Here, we provide evidence that Regulator of Calcineurin 1 (RCAN1), a feedback inhibitor of the calcium-activated protein phosphatase calcineurin (CN), acts to suppress two distinctly different mechanisms of non-shivering thermogenesis (NST): one involving the activation of UCP1 expression in white adipose tissue, the other mediated by sarcolipin (SLN) in skeletal muscle. UCP1 generates heat at the expense of reducing ATP production, whereas SLN increases ATP consumption to generate heat. Gene expression profiles demonstrate a high correlation between Rcan1 expression and metabolic syndrome. On an evolutionary timescale, in the context of limited food resources, systemic suppression of prolonged NST by RCAN1 might have been beneficial; however, in the face of caloric abundance, RCAN1-mediated suppression of these adaptive avenues of energy expenditure may now contribute to the growing epidemic of obesity.

Genetic, dietary, and sex-specific regulation of hepatic ceramides and the relationship between hepatic ceramides and IR.

Elevated hepatic ceramide levels have been implicated in both insulin resistance (IR) and hepatic steatosis. To understand the factors contributing to hepatic ceramide levels in mice of both sexes, we have quantitated ceramides in a reference population of mice, the Hybrid Mouse Diversity Panel that has been previously characterized for a variety of metabolic syndrome traits. We observed significant positive correlations between Cer(d18:1/16:0) and IR/hepatic steatosis, consistent with previous findings, although the relationship broke down between sexes, as females were less insulin resistant, but had higher Cer(d18:1/16:0) levels than males. The sex difference was due in part to testosterone-mediated repression of ceramide synthase 6. One ceramide species, Cer(d18:1/20:0), was present at higher levels in males and was associated with IR only in males. Clear evidence of gene-by-sex and gene-by-diet interactions was observed, including sex-specific genome-wide association study results. Thus, our studies show clear differences in how hepatic ceramides are regulated between the sexes, which again suggests that the physiological roles of certain hepatic ceramides differ between the sexes.

Genetic and environmental control of host-gut microbiota interactions.

Genetics provides a potentially powerful approach to dissect host-gut microbiota interactions. Toward this end, we profiled gut microbiota using 16s rRNA gene sequencing in a panel of 110 diverse inbred strains of mice. This panel has previously been studied for a wide range of metabolic traits and can be used for high-resolution association mapping. Using a SNP-based approach with a linear mixed model, we estimated the heritability of microbiota composition. We conclude that, in a controlled environment, the genetic background accounts for a substantial fraction of abundance of most common microbiota. The mice were previously studied for response to a high-fat, high-sucrose diet, and we hypothesized that the dietary response was determined in part by gut microbiota composition. We tested this using a cross-fostering strategy in which a strain showing a modest response, SWR, was seeded with microbiota from a strain showing a strong response, A×B19. Consistent with a role of microbiota in dietary response, the cross-fostered SWR pups exhibited a significantly increased response in weight gain. To examine specific microbiota contributing to the response, we identified various genera whose abundance correlated with dietary response. Among these, we chose Akkermansia muciniphila, a common anaerobe previously associated with metabolic effects. When administered to strain A×B19 by gavage, the dietary response was significantly blunted for obesity, plasma lipids, and insulin resistance. In an effort to further understand host-microbiota interactions, we mapped loci controlling microbiota composition and prioritized candidate genes. Our publicly available data provide a resource for future studies.

Genetic and hormonal control of hepatic steatosis in female and male mice.

The etiology of nonalcoholic fatty liver disease is complex and influenced by factors such as obesity, insulin resistance, hyperlipidemia, and sex. We now report a study on sex difference in hepatic steatosis in the context of genetic variation using a population of inbred strains of mice. While male mice generally exhibited higher concentration of hepatic TG levels on a high-fat high-sucrose diet, sex differences showed extensive interaction with genetic variation. Differences in percentage body fat were the best predictor of hepatic steatosis among the strains and explained about 30% of the variation in both sexes. The difference in percent gonadal fat and HDL explained 9.6% and 6.7% of the difference in hepatic TGs between the sexes, respectively. Genome-wide association mapping of hepatic TG revealed some striking differences in genetic control of hepatic steatosis between females and males. Gonadectomy increased the hepatic TG to body fat percentage ratio among male, but not female, mice. Our data suggest that the difference between the sexes in hepatic TG can be partly explained by differences in body fat distribution, plasma HDL, and genetic regulation. Future studies are required to understand the molecular interactions between sex, genetics, and the environment.

Estrogen receptor (ER)α-regulated lipocalin 2 expression in adipose tissue links obesity with breast cancer progression.

Obesity is associated with increased breast cancer (BrCA) incidence. Considering that inactivation of estrogen receptor (ER)α promotes obesity and metabolic dysfunction in women and female mice, understanding the mechanisms and tissue-specific sites of ERα action to combat metabolic-related disease, including BrCA, is of clinical importance. To study the role of ERα in adipose tissue we generated fat-specific ERα knock-out (FERKO) mice. Herein we show that ERα deletion increased adipocyte size, fat pad weight, and tissue expression and circulating levels of the secreted glycoprotein, lipocalin 2 (Lcn2), an adipokine previously associated with BrCA development. Chromatin immunoprecipitation and luciferase reporter studies showed that ERα binds the Lcn2 promoter to repress its expression. Because adipocytes constitute an important cell type of the breast microenvironment, we examined the impact of adipocyte ERα deletion on cancer cell behavior. Conditioned medium from ERα-null adipocytes and medium containing pure Lcn2 increased proliferation and migration of a subset of BrCA cells in culture. The proliferative and promigratory effects of ERα-deficient adipocyte-conditioned medium on BrCA cells was reversed by Lcn2 deletion. BrCA cell responsiveness to exogenous Lcn2 was heightened in cell types where endogenous Lcn2 expression was minimal, but components of the Lcn2 signaling pathway were enriched, i.e. SLC22A17 and 3-hydroxybutyrate dehydrogenase (BDH2). In breast tumor biopsies from women diagnosed with BrCA we found that BDH2 expression was positively associated with adiposity and circulating Lcn2 levels. Collectively these data suggest that reduction of ERα expression in adipose tissue promotes adiposity and is linked with the progression and severity of BrCA via increased adipocyte-specific Lcn2 production and enhanced tumor cell Lcn2 sensitivity.

Open chromatin profiling in mice livers reveals unique chromatin variations induced by high fat diet.

Metabolic diseases result from multiple genetic and environmental factors. We report here that one manner in which environmental factors can contribute to metabolic disease progression is through modification to chromatin. We demonstrate that high fat diet leads to chromatin remodeling in the livers of C57BL/6J mice, as compared with mice fed a control diet, and that these chromatin changes are associated with changes in gene expression. We further show that the regions of greatest variation in chromatin accessibility are targeted by liver transcription factors, including HNF4α, CCAAT/enhancer-binding protein α (CEBP/α), and FOXA1. Repeating the chromatin and gene expression profiling in another mouse strain, DBA/2J, revealed that the regions of greatest chromatin change are largely strain-specific and that integration of chromatin, gene expression, and genetic data can be used to characterize regulatory regions. Our data indicate dramatic changes in the epigenome due to diet and demonstrate strain-specific dynamics in chromatin remodeling.

CD36, but not G2A, modulates efferocytosis, inflammation, and fibrosis following bleomycin-induced lung injury.

Macrophage G2A and CD36 lipid receptors are thought to mediate efferocytosis following tissue injury and thereby prevent excessive inflammation that could compromise tissue repair. To test this, we subjected mice lacking G2A or CD36 receptor to bleomycin-induced lung injury and measured efferocytosis, inflammation, and fibrosis. Loss of CD36 (but not G2A) delayed clearance of apoptotic alveolar cells (mean 78% increase in apoptotic cells 7 days postinjury), potentiated inflammation (mean 56% increase in lung neutrophils and 75% increase in lung KC levels 7 days postinjury, 51% increase in lung macrophages 14 days postinjury), and reduced lung fibrosis (mean 41% and 29% reduction 14 and 21 days postinjury, respectively). Reduced fibrosis in CD36(-/-) mice was associated with lower levels of profibrotic TH2 cytokines (IL-9, IL-13, IL-4), decreased expression of the M2 macrophage marker Arginase-1, and reduced interstitial myofibroblasts. G2A, on the other hand, was required for optimal clearance of apoptotic neutrophils during zymosan-induced peritoneal inflammation (50.3% increase in apoptotic neutrophils and 30.6% increase in total neutrophils 24 h following zymosan administration in G2A(-/-) mice). Thus, CD36 is required for timely removal of apoptotic cells in the context of lung injury and modulates subsequent inflammatory and fibrotic processes relevant to fibrotic lung disease.

Metalloproteinase processing of HBEGF is a proximal event in the response of human aortic endothelial cells to oxidized phospholipids.

OBJECTIVE: Atherosclerosis is a chronic inflammatory disease initiated by monocyte recruitment and retention in the vessel wall. An important mediator of monocyte endothelial interaction is the chemokine interleukin (IL)-8. The oxidation products of phospholipids, including oxidized 1-palmitoyl-2-arachidonyl-sn-glycerol-3-phosphocholine (Ox-PAPC), accumulate in atherosclerotic lesions and strongly induce IL-8 in human aortic endothelial cells (HAECs). The goal of this study was to identify the proximal events leading to induction of IL-8 by Ox-PAPC in vascular endothelial cells. METHODS AND RESULTS: In a systems genetics analysis of HAECs isolated from 96 different human donors, we showed that heparin-binding EGF-like growth factor (HBEGF) transcript levels are strongly correlated to IL-8 induction by Ox-PAPC. The silencing and overexpression of HBEGF in HAECs confirmed the role of HBEGF in regulating IL-8 expression. HBEGF has been shown to be stored in an inactive form and activation is dependent on processing by a dysintegrin and metalloproteinases (ADAM) to a form that can activate the epidermal growth factor (EGF) receptor. Ox-PAPC was shown to rapidly induce HBEGF processing and EGF receptor activation in HAECs. Using siRNA we identified 3 ADAMs that regulate IL-8 induction and directly demonstrated that Ox-PAPC increases ADAM activity in the cells using a substrate cleavage assay. We provide evidence for one mechanism of Ox-PAPC activation of ADAM involving covalent binding of Ox-PAPC to cysteine on ADAM. Free thiol cysteine analogs showed inhibition of IL-8 induction by Ox-PAPC, and both a cysteine analog and a cell surface thiol blocker strongly inhibited ADAM activity induction by Ox-PAPC. Using microarray analyses, we determined that this ADAM pathway may regulate at least 30% of genes induced by Ox-PAPC in HAECs. CONCLUSIONS: This study is the first report demonstrating a role for the ADAM-HBEGF-EGF receptor axis in Ox-PAPC induction of IL-8 in HAECs. These studies highlight a role for specific ADAMs as initiators of Ox-PAPC action and provide evidence for a role of covalent interaction of Ox-PAPC in activation of ADAMs.

"Humanizing" mouse environments: Humidity, diurnal cycles and thermoneutrality.

Thermoneutral housing has been shown to promote more accurate and robust development of several pathologies in mice. Raising animal housing temperatures a few degrees may create a relatively straightforward opportunity to improve translatability of mouse models. In this commentary, we discuss the changes of physiology induced in mice housed at thermoneutrality, and review techniques for measuring systemic thermogenesis, specifically those affecting storage and mobilization of lipids in adipose depots. Environmental cues are a component of the information integrated by the brain to calculate food consumption and calorie deposition. We show that relative humidity is one of those cues, inducing a rapid sensory response that is converted to a more chronic susceptibility to obesity. Given high inter-institutional variability in the regulation of relative humidity, study reproducibility may be improved by consideration of this factor. We evaluate a "humanized" environmental cycling protocol, where mice sleep in warm temperature housing, and are cool during the wake cycle. We show that this protocol suppresses adaptation to cool exposure, with consequence for adipose-associated lipid storage. To evaluate systemic cues in mice housed at thermoneutral temperatures, we characterized the circulating lipidome, and show that sera are highly depleted in some HDL-associated phospholipids, specifically phospholipids containing the essential fatty acid, 18:2 linoleic acid, and its derivative, arachidonic acid (20:4) and related ether-phospholipids. Given the role of these fatty acids in inflammatory responses, we propose they may underlie the differences in disease progression observed at thermoneutrality.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/ß hydrolase domain 2 ( Abhd2 ) as an enzyme that metabolizes phosphatidylcholine and cardiolipin.

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/ß-hydrolase domain 2 ( Abhd2 ), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2 . The Abhd2 KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2 KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

Hybrid mouse diversity panel: a panel of inbred mouse strains suitable for analysis of complex genetic traits.

We have developed an association-based approach using classical inbred strains of mice in which we correct for population structure, which is very extensive in mice, using an efficient mixed-model algorithm. Our approach includes inbred parental strains as well as recombinant inbred strains in order to capture loci with effect sizes typical of complex traits in mice (in the range of 5% of total trait variance). Over the last few years, we have typed the hybrid mouse diversity panel (HMDP) strains for a variety of clinical traits as well as intermediate phenotypes and have shown that the HMDP has sufficient power to map genes for highly complex traits with resolution that is in most cases less than a megabase. In this essay, we review our experience with the HMDP, describe various ongoing projects, and discuss how the HMDP may fit into the larger picture of common diseases and different approaches.

Autoimmune-mediated reduction of high-density lipoprotein-cholesterol and paraoxonase 1 activity in systemic lupus erythematosus-prone gld mice.

OBJECTIVE: To characterize modifications of high-density lipoprotein (HDL) in autoimmune gld mice that may be relevant to premature atherosclerosis in systemic lupus erythematosus, and to assess their relationship to specific aspects of autoimmune disease. METHODS: HDL cholesterol (HDL-C), apolipoprotein A-I (Apo A-I), paraoxonase 1 (PON1) activity, hepatic gene expression, and HDL biogenesis were measured in aging female gld and wild-type congenic mice. Autoantibodies, lymphoid organs, and cytokines were analyzed by enzyme-linked immunosorbent assay, flow cytometry, and multiplex assay, respectively. RESULTS: Plasma HDL-C, HDL Apo A-I, and HDL-associated PON1 activity were reduced in aging gld mice in association with the development of autoimmunity, independent of changes in hepatic Apo A-I and PON1 expression or HDL biogenesis. Hepatic induction of the acute-phase reactant serum amyloid A1 resulted in its incorporation into HDL in gld mice. Deletion of the lipid-sensitive receptor G2A in gld mice (G2A-/- gld) attenuated reductions in HDL-C and PON1 activity without altering hepatic Apo A-I and PON1 expression, HDL biogenesis, or levels of acute-phase proinflammatory cytokines. Plasma anti-Apo A-I autoantibodies were elevated in aging gld mice commensurate with detectable increases in Apo A-I immune complexes. Autoantibody levels were lower in aging G2A-/- gld mice compared with gld mice, and anti-Apo A-I autoantibody levels were significantly related to HDL-C concentrations (r=-0.645, P<0.00004) and PON1 activity (r=-0.555, P<0.0007) among autoimmune gld and G2A-/- gld mice. CONCLUSION: Autoantibodies against Apo A-I contribute to reducing HDL-C and PON1 activity in autoimmune gld mice independently of hepatic HDL biogenesis, suggesting that functional impairment and premature clearance of HDL immune complexes may be principal mechanisms involved.

The Genetic Architecture of Carbon Tetrachloride-Induced Liver Fibrosis in Mice.

BACKGROUND & AIMS: Liver fibrosis is a multifactorial trait that develops in response to chronic liver injury. Our aim was to characterize the genetic architecture of carbon tetrachloride (CCl4)-induced liver fibrosis using the Hybrid Mouse Diversity Panel, a panel of more than 100 genetically distinct mouse strains optimized for genome-wide association studies and systems genetics. METHODS: Chronic liver injury was induced by CCl4 injections twice weekly for 6 weeks. Four hundred thirty-seven mice received CCl4 and 256 received vehicle, after which animals were euthanized for liver histology and gene expression. Using automated digital image analysis, we quantified fibrosis as the collagen proportionate area of the whole section, excluding normal collagen. RESULTS: We discovered broad variation in fibrosis among the Hybrid Mouse Diversity Panel strains, demonstrating a significant genetic influence. Genome-wide association analyses revealed significant and suggestive loci underlying susceptibility to fibrosis, some of which overlapped with loci identified in mouse crosses and human population studies. Liver global gene expression was assessed by RNA sequencing across the strains, and candidate genes were identified using differential expression and expression quantitative trait locus analyses. Gene set enrichment analyses identified the underlying pathways, of which stellate cell involvement was prominent, and coexpression network modeling identified modules associated with fibrosis. CONCLUSIONS: Our results provide a rich resource for the design of experiments to understand mechanisms underlying fibrosis and for rational strain selection when testing antifibrotic drugs.