Nam Che

Genetics of gene expression surveyed in maize, mouse and man

Treating messenger RNA transcript abundances as quantitative traits and mapping gene expression quantitative trait loci for these traits has been pursued in gene-specific ways. Transcript abundances often serve as a surrogate for classical quantitative traits in that the levels of expression are significantly correlated with the classical traits across members of a segregating population. The correlation structure between transcript abundances and classical traits has been used to identify susceptibility loci for complex diseases such as diabetes and allergic asthma. One study recently completed the first comprehensive dissection of transcriptional regulation in budding yeast, giving a detailed glimpse of a genome-wide survey of the genetics of gene expression. Unlike classical quantitative traits, which often represent gross clinical measurements that may be far removed from the biological processes giving rise to them, the genetic linkages associated with transcript abundance affords a closer look at cellular biochemical processes. Here we describe comprehensive genetic screens of mouse, plant and human transcriptomes by considering gene expression values as quantitative traits. We identify a gene expression pattern strongly associated with obesity in a murine cross, and observe two distinct obesity subtypes. Furthermore, we find that these obesity subtypes are under the control of different loci.

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.

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.

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. Conclusion—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.

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.

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. Developmental Dynamics 237:403–410, 2008.

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.

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. DOI: http://dx.doi.org/10.7554/eLife.05607.001