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.

Comparative analysis of proteome and transcriptome variation in mouse.

The relationships between the levels of transcripts and the levels of the proteins they encode have not been examined comprehensively in mammals, although previous work in plants and yeast suggest a surprisingly modest correlation. We have examined this issue using a genetic approach in which natural variations were used to perturb both transcript levels and protein levels among inbred strains of mice. We quantified over 5,000 peptides and over 22,000 transcripts in livers of 97 inbred and recombinant inbred strains and focused on the 7,185 most heritable transcripts and 486 most reliable proteins. The transcript levels were quantified by microarray analysis in three replicates and the proteins were quantified by Liquid Chromatography-Mass Spectrometry using O(18)-reference-based isotope labeling approach. We show that the levels of transcripts and proteins correlate significantly for only about half of the genes tested, with an average correlation of 0.27, and the correlations of transcripts and proteins varied depending on the cellular location and biological function of the gene. We examined technical and biological factors that could contribute to the modest correlation. For example, differential splicing clearly affects the analyses for certain genes; but, based on deep sequencing, this does not substantially contribute to the overall estimate of the correlation. We also employed genome-wide association analyses to map loci controlling both transcript and protein levels. Surprisingly, little overlap was observed between the protein- and transcript-mapped loci. We have typed numerous clinically relevant traits among the strains, including adiposity, lipoprotein levels, and tissue parameters. Using correlation analysis, we found that a low number of clinical trait relationships are preserved between the protein and mRNA gene products and that the majority of such relationships are specific to either the protein levels or transcript levels. Surprisingly, transcript levels were more strongly correlated with clinical traits than protein levels. In light of the widespread use of high-throughput technologies in both clinical and basic research, the results presented have practical as well as basic implications.

Systems genetics of susceptibility to obesity-induced diabetes in mice.

Inbred strains of mice are strikingly different in susceptibility to obesity-driven diabetes. For instance, deficiency in leptin receptor (db/db) leads to hyperphagia and obesity in both C57BL/6 and DBA/2 mice, but only on the DBA/2 background do the mice develop beta-cell loss leading to severe diabetes, while C57BL/6 mice are relatively resistant. To further investigate the genetic factors predisposing to diabetes, we have studied leptin receptor-deficient offspring of an F2 cross between C57BL/6J (db/+) males and DBA/2J females. The results show that the genetics of diabetes susceptibility are enormously complex and a number of quantitative trait loci (QTL) contributing to diabetes-related traits were identified, notably on chromosomes 4, 6, 7, 9, 10, 11, 12, and 19. The Chr. 4 locus is likely due to a disruption of the Zfp69 gene in C57BL/6J mice. To identify candidate genes and to model coexpression networks, we performed global expression array analysis in livers of the F2 mice. Expression QTL (eQTL) were identified and used to prioritize candidate genes at clinical trait QTL. In several cases, clusters of eQTLs colocalized with clinical trait QTLs, suggesting a common genetic basis. We constructed coexpression networks for both 5 and 12 wk old mice and identified several modules significantly associated with clinical traits. One module in 12 wk old mice was associated with several measures of hepatic fat content as well as with other lipid- and diabetes-related traits. These results add to the understanding of the complex genetic interactions contributing to obesity-induced diabetes.

High-resolution mapping of gene expression using association in an outbred mouse stock.

Quantitative trait locus (QTL) analysis is a powerful tool for mapping genes for complex traits in mice, but its utility is limited by poor resolution. A promising mapping approach is association analysis in outbred stocks or different inbred strains. As a proof of concept for the association approach, we applied whole-genome association analysis to hepatic gene expression traits in an outbred mouse population, the MF1 stock, and replicated expression QTL (eQTL) identified in previous studies of F2 intercross mice. We found that the mapping resolution of these eQTL was significantly greater in the outbred population. Through an example, we also showed how this precise mapping can be used to resolve previously identified loci (in intercross studies), which affect many different transcript levels (known as eQTL "hotspots"), into distinct regions. Our results also highlight the importance of correcting for population structure in whole-genome association studies in the outbred stock.

Hyplip2, a new gene for combined hyperlipidemia and increased atherosclerosis.

OBJECTIVE: We previously reported the mapping of a quantitative trait locus (QTL) on chromosome 15 contributing to hyperlipidemia in a cross between inbred strains MRL/MpJ (MRL) and BALB/cJ (BALB). Using marker-assisted breeding, we constructed a congenic strain in which chromosome 15 interval from MRL is placed on the genetic background of BALB. The congenic allowed us to confirm the QTL result and to further characterize the properties and location of the underlying gene. METHODS AND RESULTS: On chow and high-fat (atherogenic) diets, the congenic mice exhibited higher levels of plasma triglycerides and cholesterol than BALB mice. In response to the atherogenic diet, the congenic mice but not BALB mice exhibited a dramatic approximately 30-fold increase in atherogenic lesions accompanied by approximately 2-fold decrease in high-density lipoprotein cholesterol levels. With respect to atherosclerotic lesions and some lipid parameters, this chromosome 15 gene, designated Hyplip2, exhibited dominant inheritance. Expression array analyses suggested that Hyplip2 may influence inflammatory and bile acid synthesis pathways. Finally, we demonstrated the usefulness of subcongenic strains to narrow the locus (50 Mbp) with the goal of positionally cloning Hyplip2. CONCLUSIONS: Our data demonstrate that the Hyplip2 gene significantly contributes to combined hyperlipidemia and increased atherosclerosis in mice.

Novel QTLs for HDL levels identified in mice by controlling for Apoa2 allelic effects: confirmation of a chromosome 6 locus in a congenic strain.

Atherosclerosis is a complex disease resulting from the interaction of multiple genes, including those causing dyslipidemia. Relatively few of the causative genes have been identified. Previously, we identified Apoa2 as a major determinant of high-density lipoprotein cholesterol (HDL-C) levels in the mouse model. To identify additional HDL-C level quantitative trait loci (QTLs), while controlling for the effect of the Apoa2 locus, we performed linkage analysis in 179 standard diet-fed F(2) mice derived from strains BALB/cJ and B6.C-H25(c) (a congenic strain carrying the BALB/c Apoa2 allele). Three significant QTLs and one suggestive locus were identified. A female-specific locus mapping to chromosome 6 (Chr 6) also exhibited effects on plasma non-HDL-C, apolipoprotein AII (apoAII), apoB, and apoE levels. A Chr 6 QTL was independently isolated in a related congenic strain (C57BL/6J vs. B6.NODc6: P = 0.003 and P = 0.0001 for HDL-C and non-HDL-C levels, respectively). These data are consistent with polygenic inheritance of HDL-C levels in the mouse model and provide candidate loci for HDL-C and non-HDL-C level determination in humans.

Dissection of multigenic obesity traits in congenic mouse strains.

Previous quantitative trait locus mapping (QTL) identified multigenic obesity (MOB) loci on mouse Chromosome (Chr) 2 that influence the interrelated phenotypes of obesity, insulin resistance, and dyslipidemia. To better localize and characterize the MOB locus, three congenic mouse strains were created. Overlapping genomic intervals from the lean CAST/Ei (CAST) strain were introgressed onto an obesity-susceptible C57BL/6 (BL6) background to create proximal (15 Mb-73 Mb), middle (63 Mb-165 Mb), and distal (83 Mb-182 Mb) congenic strains. The congenic strains showed differences in obesity, insulin, and lipid traits consistent with the original QTL analysis for the locus. Importantly, characterization of the MOB congenics localized the effects of genes that underlie obesity-related traits to an introgressed interval (73-83 Mb) unique to the middle MOB congenic. Conversely, significant differences between the lipid and insulin profiles of the middle and distal MOB congenics implicated the presence of at least two genes that underlie these traits. When fed an atherogenic diet, several traits associated with metabolic syndrome were observed in the distal MOB congenic, while alterations in plasma lipoproteins were observed in the middle MOB congenic strain.