Luz Orozco

A high-resolution association mapping panel for the dissection of complex traits in mice.

Systems genetics relies on common genetic variants to elucidate biologic networks contributing to complex disease-related phenotypes. Mice are ideal model organisms for such approaches, but linkage analysis has been only modestly successful due to low mapping resolution. Association analysis in mice has the potential of much better resolution, but it is confounded by population structure and inadequate power to map traits that explain less than 10% of the variance, typical of mouse quantitative trait loci (QTL). We report a novel strategy for association mapping that combines classic inbred strains for mapping resolution and recombinant inbred strains for mapping power. Using a mixed model algorithm to correct for population structure, we validate the approach by mapping over 2500 cis-expression QTL with a resolution an order of magnitude narrower than traditional QTL analysis. We also report the fine mapping of metabolic traits such as plasma lipids. This resource, termed the Hybrid Mouse Diversity Panel, makes possible the integration of multiple data sets and should prove useful for systems-based approaches to complex traits and studies of gene-by-environment interactions.

Heme Oxygenase-1 Expression in Macrophages Plays a Beneficial Role in Atherosclerosis

Heme oxygenase (HO-1) is the rate-limiting enzyme in the catabolism of heme, which leads to the generation of biliverdin, iron, and carbon monoxide. It has been shown to have important antioxidant and antiinflammatory properties that result in a vascular antiatherogenic effect. To determine whether HO-1 expression in macrophages constitutes a significant component of the protective role in atherosclerosis, we evaluated the effect of decreased or absent HO-1 expression in peritoneal macrophages on oxidative stress and inflammation in vitro, and the effect of complete deficiency of HO-1 expression in macrophages in atherosclerotic lesion formation in vivo. We found that compared with HO-1+/+ controls, peritoneal macrophages from HO-1−/− and HO-1+/− mice exhibited (1) increased reactive oxygen species (ROS) generation, (2) increased proinflammatory cytokines such as monocyte chemotactic protein 1 (MCP-1) and interleukin 6 (IL-6), and (3) increased foam cell formation when treated with oxLDL, attributable in part to increased expression of scavenger receptor A (SR-A). Bone marrow transplantation experiments performed in lethally irradiated LDL-R null female mice, reconstituted with bone marrow from HO-1−/− versus HO-1+/+ mice, revealed that HO-1−/− reconstituted animals exhibited atherosclerotic lesions with a greater macrophage content as evaluated by immunohistochemistry and planimetric assessment. We conclude that HO-1 expression in macrophages constitutes an important component of the antiatherogenic effect by increasing antioxidant protection and decreasing the inflammatory component of atherosclerotic lesions.

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.

Mouse Genome-Wide Association and Systems Genetics Identify Asxl2 As a Regulator of Bone Mineral Density and Osteoclastogenesis

Significant advances have been made in the discovery of genes affecting bone mineral density (BMD); however, our understanding of its genetic basis remains incomplete. In the current study, genome-wide association (GWA) and co-expression network analysis were used in the recently described Hybrid Mouse Diversity Panel (HMDP) to identify and functionally characterize novel BMD genes. In the HMDP, a GWA of total body, spinal, and femoral BMD revealed four significant associations (−log10P>5.39) affecting at least one BMD trait on chromosomes (Chrs.) 7, 11, 12, and 17. The associations implicated a total of 163 genes with each association harboring between 14 and 112 genes. This list was reduced to 26 functional candidates by identifying those genes that were regulated by local eQTL in bone or harbored potentially functional non-synonymous (NS) SNPs. This analysis revealed that the most significant BMD SNP on Chr. 12 was a NS SNP in the additional sex combs like-2 (Asxl2) gene that was predicted to be functional. The involvement of Asxl2 in the regulation of bone mass was confirmed by the observation that Asxl2 knockout mice had reduced BMD. To begin to unravel the mechanism through which Asxl2 influenced BMD, a gene co-expression network was created using cortical bone gene expression microarray data from the HMDP strains. Asxl2 was identified as a member of a co-expression module enriched for genes involved in the differentiation of myeloid cells. In bone, osteoclasts are bone-resorbing cells of myeloid origin, suggesting that Asxl2 may play a role in osteoclast differentiation. In agreement, the knockdown of Asxl2 in bone marrow macrophages impaired their ability to form osteoclasts. This study identifies a new regulator of BMD and osteoclastogenesis and highlights the power of GWA and systems genetics in the mouse for dissecting complex genetic traits.

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.

Intrauterine calorie restriction affects placental DNA methylation and gene expression.

Maternal nutrient restriction causes the development of adult onset chronic diseases in the intrauterine growth restricted (IUGR) fetus. Investigations in mice have shown that either protein or calorie restriction during pregnancy leads to glucose intolerance, increased fat mass, and hypercholesterolemia in adult male offspring. Some of these phenotypes are shown to persist in successive generations. The molecular mechanisms underlying IUGR remain unclear. The placenta is a critical organ for mediating changes in the environment and the development of embryos. To shed light on molecular mechanisms that might affect placental responses to differing environments we examined placentas from mice that had been exposed to different diets. We measured gene expression and whole genome DNA methylation in both male and female placentas of mice exposed to either caloric restriction or ad libitum diets. We observed several differentially expressed pathways associated with IUGR phenotypes and, most importantly, a significant decrease in the overall methylation between these groups as well as sex-specific effects that are more pronounced in males. In addition, a set of significantly differentially methylated genes that are enriched for known imprinted genes were identified, suggesting that imprinted loci may be particularly susceptible to diet effects. Lastly, we identified several differentially methylated microRNAs that target genes associated with immunological, metabolic, gastrointestinal, cardiovascular, and neurological chronic diseases, as well as genes responsible for transplacental nutrient transfer and fetal development.

Epigenome-Wide Association of Liver Methylation Patterns and Complex Metabolic Traits in Mice

Heritable epigenetic factors can contribute to complex disease etiology. Here we examine the contribution of DNA methylation to complex traits that are precursors to heart disease, diabetes, and osteoporosis. We profiled DNA methylation in the liver using bisulfite sequencing in 90 mouse inbred strains, genome-wide expression levels, proteomics, metabolomics, and 68 clinical traits and performed epigenome-wide association studies (EWAS). We found associations with numerous clinical traits including bone density, insulin resistance, expression, and protein and metabolite levels. A large proportion of associations were unique to EWAS and were not identified using GWAS. Methylation levels were regulated by genetics largely in cis, but we also found evidence of trans regulation, and we demonstrate that genetic variation in the methionine synthase reductase gene Mtrr affects methylation of hundreds of CpGs throughout the genome. Our results indicate that natural variation in methylation levels contributes to the etiology of complex clinical traits.

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.

Genetic Architecture of Atherosclerosis in Mice: A Systems Genetics Analysis of Common Inbred Strains

Common forms of atherosclerosis involve multiple genetic and environmental factors. While human genome-wide association studies have identified numerous loci contributing to coronary artery disease and its risk factors, these studies are unable to control environmental factors or examine detailed molecular traits in relevant tissues. We now report a study of natural variations contributing to atherosclerosis and related traits in over 100 inbred strains of mice from the Hybrid Mouse Diversity Panel (HMDP). The mice were made hyperlipidemic by transgenic expression of human apolipoprotein E-Leiden (APOE-Leiden) and human cholesteryl ester transfer protein (CETP). The mice were examined for lesion size and morphology as well as plasma lipid, insulin and glucose levels, and blood cell profiles. A subset of mice was studied for plasma levels of metabolites and cytokines. We also measured global transcript levels in aorta and liver. Finally, the uptake of acetylated LDL by macrophages from HMDP mice was quantitatively examined. Loci contributing to the traits were mapped using association analysis, and relationships among traits were examined using correlation and statistical modeling. A number of conclusions emerged. First, relationships among atherosclerosis and the risk factors in mice resemble those found in humans. Second, a number of trait-loci were identified, including some overlapping with previous human and mouse studies. Third, gene expression data enabled enrichment analysis of pathways contributing to atherosclerosis and prioritization of candidate genes at associated loci in both mice and humans. Fourth, the data provided a number of mechanistic inferences; for example, we detected no association between macrophage uptake of acetylated LDL and atherosclerosis. Fifth, broad sense heritability for atherosclerosis was much larger than narrow sense heritability, indicating an important role for gene-by-gene interactions. Sixth, stepwise linear regression showed that the combined variations in plasma metabolites, including LDL/VLDL-cholesterol, trimethylamine N-oxide (TMAO), arginine, glucose and insulin, account for approximately 30 to 40% of the variation in atherosclerotic lesion area. Overall, our data provide a rich resource for studies of complex interactions underlying atherosclerosis.

Effect of 9p21.3 Coronary Artery Disease Locus Neighboring Genes on Atherosclerosis in Mice

Background— The human 9p21.3 chromosome locus has been shown to be an independent risk factor for atherosclerosis in multiple large-scale genome-wide association studies, but the underlying mechanism remains unknown. We set out to investigate the potential role of the 9p21.3 locus neighboring genes, including Mtap, the 2 isoforms of Cdkn2a, p16Ink4a and p19Arf, and Cdkn2b, in atherosclerosis using knockout mice models. Methods and Results— Gene-targeted mice for neighboring genes, including Mtap, Cdkn2a, p19Arf, and Cdkn2b, were each bred to mice carrying the human APO*E3 Leiden transgene that sensitizes the mice for atherosclerotic lesions through elevated plasma cholesterol. We found that the mice heterozygous for Mtap developed larger lesions compared with wild-type mice (49623±21650 versus 18899±9604 &mgr;m2 per section [mean±SD]; P=0.01), with morphology similar to that of wild-type mice. The Mtap heterozygous mice demonstrated changes in metabolic and methylation profiles and CD4+ cell counts. The Cdkn2a knockout mice had smaller lesions compared with wild-type and heterozygous mice, and there were no significant differences in lesion size in p19Arf and Cdkn2b mutants compared with wild type. We observed extensive, tissue-specific compensatory regulation of the Cdkn2a and Cdkn2b genes among the various knockout mice, making the effects on atherosclerosis difficult to interpret. Conclusions— Mtap plays a protective role against atherosclerosis, whereas Cdkn2a appears to be modestly proatherogenic. However, no relation was found between the 9p21 genotype and the transcription of 9p21 neighboring genes in primary human aortic vascular cells in vitro. There is extensive compensatory regulation in the highly conserved 9p21 orthologous region in mice.The 9p21.3 region of the genome has been identified as the locus with strongest association to coronary artery disease (CAD) and myocardial infarction (MI) in multiple independent large scale genome-wide association studies (GWAS).1-3 The locus is within a 58kb region that is devoid of protein coding genes, suggestive of a regulatory function (Figure 1). Interestingly, the neighboring genes in the region include well-known tumor suppressor genes, including CDKN2A and CDKN2B.4-6 The CDKN2A locus encodes a cyclin-dependent protein kinase (CDK) inhibitory protein (CKI), p16INK4A, and a p53-regulatory protein, p19ARF. The CDKN2B gene encodes another CKI, p15INK4B. Another gene in the region is methylthioadenosine phosphorylase (MTAP), which encodes a ubiquitously expressed metabolic enzyme S-methyl-5′-thioadenosine phosphorylase7 that processes the polyamine biosynthesis byproduct in the methionine salvage pathway. Loss or inactivation of MTAP has frequently been observed in a number of different human tumors, and it has been shown to have a tumor suppressive role in a mice model.8 Figure 1 The landscape of the 9p21.3 region Multiple studies demonstrated a potential role for cell cycle regulatory mechanisms in atherosclerosis progression. Previously, the master tumor suppressor gene p53 has been implicated in the development of atherosclerosis in apolipoprotein E (ApoE)-null mice9, 10, affecting both cell proliferation and apoptosis within the atheroma. Another tumor suppressor gene, p21Waf1, was also shown to increase the atheroma size in ApoE-null mice11, whereas the tumor suppressor p27Kip1 was shown to protect against atherosclerosis.12 Correlations of the 9p21 locus SNP genotype to differential expression of the neighboring genes have been observed in several studies with inconsistent findings.13-15 A knockout (KO) mouse model involving the entire region orthologous to the 9p21.3 CAD locus showed significant decreases in the expressed levels of Cdkn2a and Cdkn2b, and increased proliferation of primary smooth muscle cells (SMC) and mouse embryonic fibroblasts (MEF), although an effect on atherosclerosis in vivo was not demonstrated.16 Mice deficient in the p19Arf gene were found to have increased atherosclerotic lesions in an ApoE null background with significant attenuation of apoptosis in lesions as well as in cultured primary macrophages and vascular smooth muscle cells.17 However, to date no observation regarding atherosclerotic phenotype has been made involving the other neighboring genes. We set out to survey the 9p21.3 orthologous region using knockout mice models to systematically address the role of the neighboring protein-coding genes in atherosclerosis. We chose the APOE*3 Leiden sensitizing model because it is dominant, simplifying the construction of the models, and also because it exhibits relatively modest elevations of cholesterol, more realistically modeling the human disease than other widely used models.