The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

Hematological parameters, including red and white blood cell counts and hemoglobin concentration, are widely used clinical indicators of health and disease. These traits are tightly regulated in healthy individuals and are under genetic control. Mutations in key genes that affect hematological parameters have important phenotypic consequences, including multiple variants that affect susceptibility to malarial disease. However, most variation in hematological traits is continuous and is presumably influenced by multiple loci and variants with small phenotypic effects. We used a newly developed mouse resource population, the Collaborative Cross (CC), to identify genetic determinants of hematological parameters. We surveyed the eight founder strains of the CC and performed a mapping study using 131 incipient lines of the CC. Genome scans identified quantitative trait loci for several hematological parameters, including mean red cell volume (Chr 7 and Chr 14), white blood cell count (Chr 18), percent neutrophils/lymphocytes (Chr 11), and monocyte number (Chr 1). We used evolutionary principles and unique bioinformatics resources to reduce the size of candidate intervals and to view functional variation in the context of phylogeny. Many quantitative trait loci regions could be narrowed sufficiently to identify a small number of promising candidate genes. This approach not only expands our knowledge about hematological traits but also demonstrates the unique ability of the CC to elucidate the genetic architecture of complex traits.

A diverse suite of effector immune responses provide protection against various pathogens. However, the array of effector responses must be immunologically regulated to limit pathogen- and immune-associated damage. CD4+Foxp3+ regulatory T cells (Treg) calibrate immune responses; however, how Treg cells adapt to control different effector responses is unclear. To investigate the molecular mechanism of Treg diversity we used whole genome expression profiling and next generation small RNA sequencing of Treg cells isolated from type-1 or type-2 inflamed tissue following Leishmania major or Schistosoma mansoni infection, respectively. In-silico analyses identified two miRNA “regulatory hubs” miR-10a and miR-182 as critical miRNAs in Th1- or Th2-associated Treg cells, respectively. Functionally and mechanistically, in-vitro and in-vivo systems identified that an IL-12/IFNγ axis regulated miR-10a and its putative transcription factor, Creb. Importantly, reduced miR-10a in Th1-associated Treg cells was critical for Treg function and controlled a suite of genes preventing IFNγ production. In contrast, IL-4 regulated miR-182 and cMaf in Th2-associed Treg cells, which mitigated IL-2 secretion, in part through repression of IL2-promoting genes. Together, this study indicates that CD4+Foxp3+ cells can be shaped by local environmental factors, which orchestrate distinct miRNA pathways preserving Treg stability and suppressor function.

Airway allergen exposure induces inflammation among individuals with atopy that is characterized by altered airway gene expression, elevated levels of T helper type 2 cytokines, mucus hypersecretion, and airflow obstruction. To identify the genetic determinants of the airway allergen response, we employed a systems genetics approach. We applied a house dust mite mouse model of allergic airway disease to 151 incipient lines of the Collaborative Cross, a new mouse genetic reference population, and measured serum IgE, airway eosinophilia, and gene expression in the lung. Allergen-induced serum IgE and airway eosinophilia were not correlated. We detected quantitative trait loci (QTL) for airway eosinophilia on chromosome (Chr) 11 (71.802-87.098 megabases [Mb]) and allergen-induced IgE on Chr 4 (13.950-31.660 Mb). More than 4,500 genes expressed in the lung had gene expression QTL (eQTL), the majority of which were located near the gene itself. However, we also detected approximately 1,700 trans-eQTL, and many of these trans-eQTL clustered into two regions on Chr 2. We show that one of these loci (at 147.6 Mb) is associated with the expression of more than 100 genes, and, using bioinformatics resources, fine-map this locus to a 53 kb-long interval. We also use the gene expression and eQTL data to identify a candidate gene, Tlcd2, for the eosinophil QTL. Our results demonstrate that hallmark allergic airway disease phenotypes are associated with distinct genetic loci on Chrs 4 and 11, and that gene expression in the allergically inflamed lung is controlled by both cis and trans regulatory factors.

Allergic asthma is a complex disease characterized in part by granulocytic inflammation of the airways. In addition to eosinophils, neutrophils (PMN) are also present, particularly in cases of severe asthma. We sought to identify the genetic determinants of neutrophilic inflammation in a mouse model of house dust mite (HDM)-induced asthma. We applied an HDM model of allergic asthma to the eight founder strains of the Collaborative Cross (CC) and 151 incipient lines of the CC (preCC). Lung lavage fluid was analyzed for PMN count and the concentration of CXCL1, a hallmark PMN chemokine. PMN and CXCL1 were strongly correlated in preCC mice. We used quantitative trait locus (QTL) mapping to identify three variants affecting PMN, one of which colocalized with a QTL for CXCL1 on chromosome (Chr) 7. We used lung eQTL data to implicate a variant in the gene Zfp30 in the CXCL1/PMN response. This genetic variant regulates both CXCL1 and PMN by altering Zfp30 expression, and we model the relationships between the QTL and these three endophenotypes. We show that Zfp30 is expressed in airway epithelia in the normal mouse lung and that altering Zfp30 expression in vitro affects CXCL1 responses to an immune stimulus. Our results provide strong evidence that Zfp30 is a novel regulator of neutrophilic airway inflammation.

Asthma is etiologically and clinically heterogeneous, making the genomic basis of asthma difficult to identify. We exploited the straindependence of a murine model of allergic airway disease to identify different genomic responses in the lung. BALB/cJ and C57BL/6J mice were sensitized with the immunodominant allergen from the Dermatophagoides pteronyssinus species of house dust mite (Der p 1), without exogenous adjuvant, and the mice then underwent a single challenge with Der p 1. Allergic inflammation, serum antibody titers, mucous metaplasia, and airway hyperresponsiveness were evaluated 72 hours after airway challenge. Whole-lung gene expression analyses were conducted to identify genomic responses to allergen challenge. Der p 1-challenged BALB/cJ mice produced all the key features of allergic airway disease. In comparison, C57BL/6J mice produced exaggerated Th2-biased responses and inflammation, but exhibited an unexpected decrease in airway hyperresponsiveness compared with control mice. Lung gene expression analysis revealed genes that were shared by both strains and a set of down-regulated genes unique to C57BL/6J mice, including several G-protein-coupled receptors involved in airway smooth muscle contraction, most notably the M2 muscarinic receptor, which we show is expressed in airway smooth muscle and was decreased at the protein level after challenge with Der p 1. Murine strain-dependent genomic responses in the lung offer insights into the different biological pathways that develop after allergen challenge. This study of two different murine strains demonstrates that inflammation and airway hyperresponsiveness can be decoupled, and suggests that the downmodulation of expression of G-protein-coupled receptors involved in regulating airway smooth muscle contraction may contribute to this dissociation.

Genome-wide association studies (GWAS) have strongly implicated MIR137 (the gene encoding the microRNA miR-137) in schizophrenia. A parsimonious hypothesis is that a pathway regulated by miR-137 is important in the etiology of schizophrenia. Full evaluation of this hypothesis requires more definitive knowledge about biological targets of miR-137, which is currently lacking. Our goals were to expand knowledge of the biology of miR-137 by identifying its empirical targets, and to test whether the resulting lists of direct and indirect targets were enriched for genes and pathways involved in risk for schizophrenia. We overexpressed miR-137 in a human neural stem cell line and analyzed gene expression changes at 24 and 48 h using RNA sequencing. Following miR-137 overexpression, 202 and 428 genes were differentially expressed after 24 and 48 h. Genes differentially expressed at 24 h were enriched for transcription factors and cell cycle genes, and differential expression at 48 h affected a wider variety of pathways. Pathways implicated in schizophrenia were upregulated in the 48 h findings (major histocompatibility complex, synapses, FMRP interacting RNAs and calcium channels). Critically, differentially expressed genes at 48 h were enriched for smaller association P-values in the largest published schizophrenia GWAS. This work provides empirical support for a role of miR-137 in the etiology of schizophrenia.

Mucus hyper-secretion is a hallmark feature of asthma and other muco-obstructive airway diseases. The mucin proteins MUC5AC and MUC5B are the major glycoprotein components of mucus and have critical roles in airway defense. Despite the biomedical importance of these two proteins, the loci that regulate them in the context of natural genetic variation have not been studied. To identify genes that underlie variation in airway mucin levels, we performed genetic analyses in founder strains and incipient lines of the Collaborative Cross (CC) in a house dust mite mouse model of asthma. CC founder strains exhibited significant differences in MUC5AC and MUC5B, providing evidence of heritability. Analysis of gene and protein expression of and in incipient CC lines ( = 154) suggested that post-transcriptional events were important regulators of mucin protein content in the airways. Quantitative trait locus (QTL) mapping identified distinct, protein QTL for MUC5AC (chromosome 13) and MUC5B (chromosome 2). These two QTL explained 18 and 20% of phenotypic variance, respectively. Examination of the MUC5B QTL allele effects and subsequent phylogenetic analysis allowed us to narrow the MUC5B QTL and identify as a candidate gene. mRNA and protein expression were upregulated in parallel to MUC5B after allergen challenge, and knockout mice exhibited higher MUC5B expression. Thus, BPIFB1 is a novel regulator of MUC5B.