Unraveling inflammatory responses using systems genetics and gene-environment interactions in macrophages.

Many common diseases have an important inflammatory component mediated in part by macrophages. Here we used a systems genetics strategy to examine the role of common genetic variation in macrophage responses to inflammatory stimuli. We examined genome-wide transcript levels in macrophages from 92 strains of the Hybrid Mouse Diversity Panel. We exposed macrophages to control media, bacterial lipopolysaccharide (LPS), or oxidized phospholipids. We performed association mapping under each condition and identified several thousand expression quantitative trait loci (eQTL), gene-by-environment interactions, and eQTL "hot spots" that specifically control LPS responses. We used siRNA knockdown of candidate genes to validate an eQTL hot spot in chromosome 8 and identified the gene 2310061C15Rik as a regulator of inflammatory responses in macrophages. We have created a public database where the data presented here can be used as a resource for understanding many common inflammatory traits that are modeled in the mouse and for the dissection of regulatory relationships between genes.

Transcription Factor MAFF (MAF Basic Leucine Zipper Transcription Factor F) Regulates an Atherosclerosis Relevant Network Connecting Inflammation and Cholesterol Metabolism.

BACKGROUND: Coronary artery disease (CAD) is a multifactorial condition with both genetic and exogenous causes. The contribution of tissue-specific functional networks to the development of atherosclerosis remains largely unclear. The aim of this study was to identify and characterize central regulators and networks leading to atherosclerosis. METHODS: Based on several hundred genes known to affect atherosclerosis risk in mouse (as demonstrated in knockout models) and human (as shown by genome-wide association studies), liver gene regulatory networks were modeled. The hierarchical order and regulatory directions of genes within the network were based on Bayesian prediction models, as well as experimental studies including chromatin immunoprecipitation DNA-sequencing, chromatin immunoprecipitation mass spectrometry, overexpression, small interfering RNA knockdown in mouse and human liver cells, and knockout mouse experiments. Bioinformatics and correlation analyses were used to clarify associations between central genes and CAD phenotypes in both human and mouse. RESULTS: The transcription factor MAFF (MAF basic leucine zipper transcription factor F) interacted as a key driver of a liver network with 3 human genes at CAD genome-wide association studies loci and 11 atherosclerotic murine genes. Most importantly, expression levels of the low-density lipoprotein receptor (LDLR) gene correlated with MAFF in 600 CAD patients undergoing bypass surgery (STARNET [Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task]) and a hybrid mouse diversity panel involving 105 different inbred mouse strains. Molecular mechanisms of MAFF were tested in noninflammatory conditions and showed positive correlation between MAFF and LDLR in vitro and in vivo. Interestingly, after lipopolysaccharide stimulation (inflammatory conditions), an inverse correlation between MAFF and LDLR in vitro and in vivo was observed. Chromatin immunoprecipitation mass spectrometry revealed that the human CAD genome-wide association studies candidate BACH1 (BTB domain and CNC homolog 1) assists MAFF in the presence of lipopolysaccharide stimulation with respective heterodimers binding at the MAF recognition element of the LDLR promoter to transcriptionally downregulate LDLR expression. CONCLUSIONS: The transcription factor MAFF was identified as a novel central regulator of an atherosclerosis/CAD-relevant liver network. MAFF triggered context-specific expression of LDLR and other genes known to affect CAD risk. Our results suggest that MAFF is a missing link between inflammation, lipid and lipoprotein metabolism, and a possible treatment target.

Regulation of NF-κB signaling by oxidized glycerophospholipid and IL-1ß induced miRs-21-3p and -27a-5p in human aortic endothelial cells.

Exposure of endothelial cells (ECs) to agents such as oxidized glycerophospholipids (oxGPs) and cytokines, known to accumulate in atherosclerotic lesions, perturbs the expression of hundreds of genes in ECs involved in inflammatory and other biological processes. We hypothesized that microRNAs (miRNAs) are involved in regulating the inflammatory response in human aortic endothelial cells (HAECs) in response to oxGPs and interleukin 1ß (IL-1ß). Using next-generation sequencing and RT-quantitative PCR, we characterized the profile of expressed miRNAs in HAECs pre- and postexposure to oxGPs. Using this data, we identified miR-21-3p and miR-27a-5p to be induced 3- to 4-fold in response to oxGP and IL-1ß treatment compared with control treatment. Transient overexpression of miR-21-3p and miR-27a-5p resulted in the downregulation of 1,253 genes with 922 genes overlapping between the two miRNAs. Gene Ontology functional enrichment analysis predicted that the two miRNAs were involved in the regulation of nuclear factor κB (NF-κB) signaling. Overexpression of these two miRNAs leads to changes in p65 nuclear translocation. Using 3' untranslated region luciferase assay, we identified 20 genes within the NF-κB signaling cascade as putative targets of miRs-21-3p and -27a-5p, implicating these two miRNAs as modulators of NF-κB signaling in ECs.

Flavin containing monooxygenase 3 exerts broad effects on glucose and lipid metabolism and atherosclerosis.

We performed silencing and overexpression studies of flavin containing monooxygenase (FMO) 3 in hyperlipidemic mouse models to examine its effects on trimethylamine N-oxide (TMAO) levels and atherosclerosis. Knockdown of hepatic FMO3 in LDL receptor knockout mice using an antisense oligonucleotide resulted in decreased circulating TMAO levels and atherosclerosis. Surprisingly, we also observed significant decreases in hepatic lipids and in levels of plasma lipids, ketone bodies, glucose, and insulin. FMO3 overexpression in transgenic mice, on the other hand, increased hepatic and plasma lipids. Global gene expression analyses suggested that these effects of FMO3 on lipogenesis and gluconeogenesis may be mediated through the PPARα and Kruppel-like factor 15 pathways. In vivo and in vitro results were consistent with the concept that the effects were mediated directly by FMO3 rather than trimethylamine/TMAO; in particular, overexpression of FMO3 in the human hepatoma cell line, Hep3B, resulted in significantly increased glucose secretion and lipogenesis. Our results indicate a major role for FMO3 in modulating glucose and lipid homeostasis in vivo, and they suggest that pharmacologic inhibition of FMO3 to reduce TMAO levels would be confounded by metabolic interactions.

Genetic regulation of human adipose microRNA expression and its consequences for metabolic traits.

The genetics of messenger RNA (mRNA) expression has been extensively studied in humans and other organisms, but little is known about genetic factors contributing to microRNA (miRNA) expression. We examined natural variation of miRNA expression in adipose tissue in a population of 200 men who have been carefully characterized for metabolic syndrome (MetSyn) phenotypes as part of the Metabolic Syndrome in Men (METSIM) study. We genotyped the subjects using high-density single-nucleotide polymorphism microarrays and quantified the mRNA abundance using genome-wide expression arrays and miRNA abundance using next-generation sequencing. We reliably quantified 356 miRNA species that were expressed in human adipose tissue, a limited number of which made up most of the expressed miRNAs. We mapped the miRNA abundance as an expression quantitative trait and determined cis regulation of expression for nine of the miRNAs and of the processing of one miRNA (miR-28). The degree of genetic variation of miRNA expression was substantially less than that of mRNAs. For the majority of the miRNAs, genetic regulation of expression was independent of the expression of mRNA from which the miRNA is transcribed. We also showed that for 108 miRNAs, mapped reads displayed widespread variation from the canonical sequence. We found a total of 24 miRNAs to be significantly associated with MetSyn traits. We suggest a regulatory role for miR-204-5p which was predicted to inhibit acetyl coenzyme A carboxylase ß, a key fatty acid oxidation enzyme that has been shown to play a role in regulating body fat and insulin resistance in adipose tissue.

Genetic regulation of mouse liver metabolite levels.

We profiled and analyzed 283 metabolites representing eight major classes of molecules including Lipids, Carbohydrates, Amino Acids, Peptides, Xenobiotics, Vitamins and Cofactors, Energy Metabolism, and Nucleotides in mouse liver of 104 inbred and recombinant inbred strains. We find that metabolites exhibit a wide range of variation, as has been previously observed with metabolites in blood serum. Using genome-wide association analysis, we mapped 40% of the quantified metabolites to at least one locus in the genome and for 75% of the loci mapped we identified at least one candidate gene by local expression QTL analysis of the transcripts. Moreover, we validated 2 of 3 of the significant loci examined by adenoviral overexpression of the genes in mice. In our GWAS results, we find that at significant loci the peak markers explained on average between 20 and 40% of variation in the metabolites. Moreover, 39% of loci found to be regulating liver metabolites in mice were also found in human GWAS results for serum metabolites, providing support for similarity in genetic regulation of metabolites between mice and human. We also integrated the metabolomic data with transcriptomic and clinical phenotypic data to evaluate the extent of co-variation across various biological scales.

Network for activation of human endothelial cells by oxidized phospholipids: a critical role of heme oxygenase 1.

RATIONALE: Oxidized palmitoyl arachidonyl phosphatidylcholine (Ox-PAPC) accumulates in atherosclerotic lesions, is proatherogenic, and influences the expression of more than 1000 genes in endothelial cells. OBJECTIVE: To elucidate the major pathways involved in Ox-PAPC action, we conducted a systems analysis of endothelial cell gene expression after exposure to Ox-PAPC. METHODS AND RESULTS: We used the variable responses of primary endothelial cells from 149 individuals exposed to Ox-PAPC to construct a network that consisted of 11 groups of genes, or modules. Modules were enriched for a broad range of Gene Ontology pathways, some of which have not been identified previously as major Ox-PAPC targets. Further validating our method of network construction, modules were consistent with relationships established by cell biology studies of Ox-PAPC effects on endothelial cells. This network provides novel hypotheses about molecular interactions, as well as candidate molecular regulators of inflammation and atherosclerosis. We validated several hypotheses based on network connections and genomic association. Our network analysis predicted that the hub gene CHAC1 (cation transport regulator homolog 1) was regulated by the ATF4 (activating transcription factor 4) arm of the unfolded protein response pathway, and here we showed that ATF4 directly activates an element in the CHAC1 promoter. We showed that variation in basal levels of heme oxygenase 1 (HMOX1) contribute to the response to Ox-PAPC, consistent with its position as a hub in our network. We also identified G-protein-coupled receptor 39 (GPR39) as a regulator of HMOX1 levels and showed that it modulates the promoter activity of HMOX1. We further showed that OKL38/OSGN1 (oxidative stress-induced growth inhibitor), the hub gene in the blue module, is a key regulator of both inflammatory and antiinflammatory molecules. CONCLUSIONS: Our systems genetics approach has provided a broad view of the pathways involved in the response of endothelial cells to Ox-PAPC and also identified novel regulatory mechanisms.

NF-E2-related factor 2 promotes atherosclerosis by effects on plasma lipoproteins and cholesterol transport that overshadow antioxidant protection.

OBJECTIVE: To test the hypothesis that NF-E2-related factor 2 (Nrf2) expression plays an antiatherogenic role by its vascular antioxidant and anti-inflammatory properties. METHODS AND RESULTS: Nrf2 is an important transcription factor that regulates the expression of phase 2 detoxifying enzymes and antioxidant genes. Its expression in vascular cells appears to be an important factor in the protection against vascular oxidative stress and inflammation. We developed Nrf2 heterozygous (HET) and homozygous knockout (KO) mice on an apolipoprotein (apo) E-null background by sequential breeding, resulting in Nrf2(-/-), apoE(-/-) (KO), Nrf2(-/+), apoE(-/-) (HET) and Nrf2(+/+), and apoE(-/-) wild-type littermates. KO mice exhibited decreased levels of antioxidant genes with evidence of increased reactive oxygen species generation compared with wild-type controls. Surprisingly, KO males exhibited 47% and 53% reductions in the degree of aortic atherosclerosis compared with HET or wild-type littermates, respectively. Decreased atherosclerosis in KO mice correlated with lower plasma total cholesterol in a sex-dependent manner. KO mice also had a decreased hepatic cholesterol content and a lower expression of lipogenic genes, suggesting that hepatic lipogenesis could be reduced. In addition, KO mice exhibited atherosclerotic plaques characterized by a lesser macrophage component and decreased foam cell formation in an in vitro lipid-loading assay. This was associated with a lower rate of cholesterol influx, mediated in part by decreased expression of the scavenger receptor CD36. CONCLUSIONS: Nrf2 expression unexpectedly promotes atherosclerotic lesion formation in a sex-dependent manner, most likely by a combination of systemic metabolic and local vascular effects.

Copy number variation influences gene expression and metabolic traits in mice.

Copy number variants (CNVs) are genomic segments which are duplicated or deleted among different individuals. CNVs have been implicated in both Mendelian and complex traits, including immune and behavioral disorders, but the study of the mechanisms by which CNVs influence gene expression and clinical phenotypes in humans is complicated by the limited access to tissues and by population heterogeneity. We now report studies of the effect of 19 CNVs on gene expression and metabolic traits in a mouse intercross between strains C57BL/6J and C3H/HeJ. We found that 83% of genes predicted to occur within CNVs were differentially expressed. The expression of most CNV genes was correlated with copy number, but we also observed evidence that gene expression was altered in genes flanking CNVs, suggesting that CNVs may contain regulatory elements for these genes. Several CNVs mapped to hotspots, genomic regions influencing expression of tens or hundreds of genes. Several metabolic traits including cholesterol, triglycerides, glucose and body weight mapped to three CNVs in the genome, in mouse chromosomes 1, 4 and 17. Predicted CNV genes, such as Itlna, Defcr-1, Trim12 and Trim34 were highly correlated with these traits. Our results suggest that CNVs have a significant impact on gene expression and that CNVs may be playing a role in the mechanisms underlying metabolic traits in mice.

Integration of Multi-omics Data from Mouse Diversity Panel Highlights Mitochondrial Dysfunction in Non-alcoholic Fatty Liver Disease.

The etiology of non-alcoholic fatty liver disease (NAFLD), the most common form of chronic liver disease, is poorly understood. To understand the causal mechanisms underlying NAFLD, we conducted a multi-omics, multi-tissue integrative study using the Hybrid Mouse Diversity Panel, consisting of ∼100 strains of mice with various degrees of NAFLD. We identified both tissue-specific biological processes and processes that were shared between adipose and liver tissues. We then used gene network modeling to predict candidate regulatory genes of these NAFLD processes, including Fasn, Thrsp, Pklr, and Chchd6. In vivo knockdown experiments of the candidate genes improved both steatosis and insulin resistance. Further in vitro testing demonstrated that downregulation of both Pklr and Chchd6 lowered mitochondrial respiration and led to a shift toward glycolytic metabolism, thus highlighting mitochondria dysfunction as a key mechanistic driver of NAFLD.

A Strategy for Discovery of Endocrine Interactions with Application to Whole-Body Metabolism.

Inter-tissue communication via secreted proteins has been established as a vital mechanism for proper physiologic homeostasis. Here, we report a bioinformatics framework using a mouse reference population, the Hybrid Mouse Diversity Panel (HMDP), which integrates global multi-tissue expression data and publicly available resources to identify and functionally annotate novel circuits of tissue-tissue communication. We validate this method by showing that we can identify known as well as novel endocrine factors responsible for communication between tissues. We further show the utility of this approach by identification and mechanistic characterization of two new endocrine factors. Adipose-derived Lipocalin-5 is shown to enhance skeletal muscle mitochondrial function, and liver-secreted Notum promotes browning of white adipose tissue, also known as "beiging." We demonstrate the general applicability of the method by providing in vivo evidence for three additional novel molecules mediating tissue-tissue interactions.

Author Correction: Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition.

In the version of this article originally published, minus signs were missing from the three ß-values for BMI given in Table 1. The errors have been corrected in the HTML and PDF versions of the article.

Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition.

Individual risk of type 2 diabetes (T2D) is modified by perturbations to the mass, distribution and function of adipose tissue. To investigate the mechanisms underlying these associations, we explored the molecular, cellular and whole-body effects of T2D-associated alleles near KLF14. We show that KLF14 diabetes-risk alleles act in adipose tissue to reduce KLF14 expression and modulate, in trans, the expression of 385 genes. We demonstrate, in human cellular studies, that reduced KLF14 expression increases pre-adipocyte proliferation but disrupts lipogenesis, and in mice, that adipose tissue-specific deletion of Klf14 partially recapitulates the human phenotype of insulin resistance, dyslipidemia and T2D. We show that carriers of the KLF14 T2D risk allele shift body fat from gynoid stores to abdominal stores and display a marked increase in adipocyte cell size, and that these effects on fat distribution, and the T2D association, are female specific. The metabolic risk associated with variation at this imprinted locus depends on the sex both of the subject and of the parent from whom the risk allele derives.

The TMAO-Producing Enzyme Flavin-Containing Monooxygenase 3 Regulates Obesity and the Beiging of White Adipose Tissue.

Emerging evidence suggests that microbes resident in the human intestine represent a key environmental factor contributing to obesity-associated disorders. Here, we demonstrate that the gut microbiota-initiated trimethylamine N-oxide (TMAO)-generating pathway is linked to obesity and energy metabolism. In multiple clinical cohorts, systemic levels of TMAO were observed to strongly associate with type 2 diabetes. In addition, circulating TMAO levels were associated with obesity traits in the different inbred strains represented in the Hybrid Mouse Diversity Panel. Further, antisense oligonucleotide-mediated knockdown or genetic deletion of the TMAO-producing enzyme flavin-containing monooxygenase 3 (FMO3) conferred protection against obesity in mice. Complimentary mouse and human studies indicate a negative regulatory role for FMO3 in the beiging of white adipose tissue. Collectively, our studies reveal a link between the TMAO-producing enzyme FMO3 and obesity and the beiging of white adipose tissue.

The genetic architecture of NAFLD among inbred strains of mice.

To identify genetic and environmental factors contributing to the pathogenesis of non-alcoholic fatty liver disease, we examined liver steatosis and related clinical and molecular traits in more than 100 unique inbred mouse strains, which were fed a diet rich in fat and carbohydrates. A >30-fold variation in hepatic TG accumulation was observed among the strains. Genome-wide association studies revealed three loci associated with hepatic TG accumulation. Utilizing transcriptomic data from the liver and adipose tissue, we identified several high-confidence candidate genes for hepatic steatosis, including Gde1, a glycerophosphodiester phosphodiesterase not previously implicated in triglyceride metabolism. We confirmed the role of Gde1 by in vivo hepatic over-expression and shRNA knockdown studies. We hypothesize that Gde1 expression increases TG production by contributing to the production of glycerol-3-phosphate. Our multi-level data, including transcript levels, metabolite levels, and gut microbiota composition, provide a framework for understanding genetic and environmental interactions underlying hepatic steatosis.

Intergenerational genomic DNA methylation patterns in mouse hybrid strains.

BACKGROUND: DNA methylation is a contributing factor to both rare and common human diseases, and plays a major role in development and gene silencing. While the variation of DNA methylation among individuals has been partially characterized, the degree to which methylation patterns are preserved across generations is still poorly understood. To determine the extent of methylation differences between two generations of mice we examined DNA methylation patterns in the livers of eight parental and F1 mice from C57BL/6J and DBA/2J mouse strains using bisulfite sequencing. RESULTS: We find a large proportion of reproducible methylation differences between C57BL/6J and DBA/2J chromosomes in CpGs, which are highly heritable between parent and F1 mice. We also find sex differences in methylation levels in 396 genes, and 11% of these are differentially expressed between females and males. Using a recently developed approach to identify allelically methylated regions independently of genotypic differences, we identify 112 novel putative imprinted genes and microRNAs, and validate imprinting at the RNA level in 10 of these genes. CONCLUSIONS: The majority of DNA methylation differences among individuals are associated with genetic differences, and a much smaller proportion of these epigenetic differences are due to sex, imprinting or stochastic intergenerational effects. Epigenetic differences can be a determining factor in heritable traits and should be considered in association studies for molecular and clinical traits, as we observed that methylation differences in the mouse model are highly heritable and can have functional consequences on molecular traits such as gene expression.

Muscleblind-like 2 (Mbnl2) -deficient mice as a model for myotonic dystrophy.

Myotonic dystrophy (DM), the most common adult-onset muscular dystrophy, is caused by CTG or CCTG microsatellite repeat expansions. Expanded DM mRNA microsatellite repeats are thought to accumulate in the nucleus, sequester Muscleblind proteins, and interfere with alternative mRNA splicing. Muscleblind2 (Mbnl2) is a member of the family of Muscleblind RNA binding proteins (that also include Mbnl1 and Mbnl3) that are known to bind CTG/CCTG RNA repeats. Recently, it was demonstrated that Mbnl1-deficient mice have characteristic features of human DM, including myotonia and defective chloride channel expression. Here, we demonstrate that Mbnl2-deficient mice also develop myotonia and have skeletal muscle pathology consistent with human DM. We also find defective expression and mRNA splicing of the chloride channel (Clcn1) in skeletal muscle that likely contributes to the myotonia phenotype. Our results support the hypothesis that Muscleblind proteins and specifically MBNL2 contribute to the pathogenesis of human DM.

Myocardin-related transcription factor B is required for normal mouse vascular development and smooth muscle gene expression.

Smooth muscle gene expression is required for the proper development and function of multiple organ systems. Expression of smooth muscle genes is critical for contractile function and tissue architectural integrity. One critical transcription factor for smooth muscle gene expression is the Serum Response Factor (SRF). SRF is expressed ubiquitously, but tissue-specific transcriptional regulation is conferred by its binding to cofactors such as myocardin. Myocardin-related transcription factor B (MRTF-B) is a member of a family of genes (Myocardin(Myocd),Myocardin-related transcription factor A(MRTF-A),MRTF-B) that provides tissue-specificity and potentiate SRF-dependent transcription. Unlike myocardin, which is expressed specifically in smooth and cardiac muscle, MRTF-B is expressed in a wide variety of tissues. To examine the function of MRTF-B, we generated mice containing an insertional mutation of MRTF-B. MRTF-B homozygous mutants die in late gestation with vascular defects and liver hemorrhage. At E9.5, MRTF-B is expressed strongly in the septum transversum mesoderm critical for development of the vitelline system that produces the liver sinusoids and portal venous system. MRTF-B deficiency results in defective smooth muscle gene expression in the liver sinusoids, vitelline veins, and yolk sac, which contributes significantly to the lethal phenotype. These data support our hypothesis that MRTF-B has a unique role in regulating smooth muscle genes important for liver, yolk sac, and portal vascular development.

Identification of Abcc6 as the major causal gene for dystrophic cardiac calcification in mice through integrative genomics.

The genetic factors contributing to the complex disorder of myocardial calcification are largely unknown. Using a mouse model, we fine-mapped the major locus (Dyscalc1) contributing to the dystrophic cardiac calcification (DCC) to an 840-kb interval containing 38 genes. We then identified the causal gene by using an approach integrating genetic segregation and expression array analyses to identify, on a global scale, cis-acting DNA variations that perturb gene expression. By studying two intercrosses, in which the DCC trait segregates, a single candidate gene (encoding the ATP-binding cassette transporter ABCC6) was identified. Transgenic complementation confirmed Abcc6 as the underlying causal gene for Dyscalc1. We demonstrate that in the cross, the expression of Abcc6 is highly correlated with the local mineralization regulatory system and the BMP2-Wnt signaling pathway known to be involved in the systemic regulation of calcification, suggesting potential pathways for the action of Abcc6 in DCC. Our results demonstrate the power of the integrative genomics in discovering causal genes and pathways underlying complex traits.