Genome-wide association studies using commercially available outbred mice can detect genes involved in phenotypes of biomedical interest. Useful populations need high-frequency alleles to ensure high power to detect quantitative trait loci (QTLs), low linkage disequilibrium between markers to obtain accurate mapping resolution, and an absence of population structure to prevent false positive associations. We surveyed 66 colonies for inbreeding, genetic diversity, and linkage disequilibrium, and we demonstrate that some have haplotype blocks of less than 100 Kb, enabling gene-level mapping resolution. The same alleles contribute to variation in different colonies, so that when mapping progress stalls in one, another can be used in its stead. Colonies are genetically diverse: 45% of the total genetic variation is attributable to differences between colonies. However, quantitative differences in allele frequencies, rather than the existence of private alleles, are responsible for these population differences. The colonies derive from a limited pool of ancestral haplotypes resembling those found in inbred strains: over 95% of sequence variants segregating in outbred populations are found in inbred strains. Consequently it is possible to impute the sequence of any mouse from a dense SNP map combined with inbred strain sequence data, which opens up the possibility of cataloguing and testing all variants for association, a situation that has so far eluded studies in completely outbred populations. We demonstrate the colonies' potential by identifying a deletion in the promoter of H2-Ea as the molecular change that strongly contributes to setting the ratio of CD4+ and CD8+ lymphocytes.

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.

No abstract available

The cytomegalovirus resistance locus Cmv3 has been linked to an epistatic interaction between two loci: a Natural Killer (NK) cell receptor gene and the major histocompatibility complex class I (MHC-I) locus. To demonstrate the interaction between Cmv3 and H2(k), we generated double congenic mice between MA/My and BALB.K mice and an F(2) cross between FVB/N (H-2(q)) and BALB.K (H2(k)) mice, two strains susceptible to mouse cytomegalovirus (MCMV). Only mice expressing H2(k) in conjunction with Cmv3(MA/My) or Cmv3(FVB) were resistant to MCMV infection. Subsequently, an F(3) cross was carried out between transgenic FVB/H2-D(k) and MHC-I deficient mice in which only the progeny expressing Cmv3(FVB) and a single H2-D(k) class-I molecule completely controlled MCMV viral loads. This phenotype was shown to be NK cell-dependent and associated with subsequent NK cell proliferation. Finally, we demonstrated that a number of H2(q) alleles influence the expression level of H2(q) molecules, but not intrinsic functional properties of NK cells; viral loads, however, were quantitatively proportional to the number of H2(q) alleles. Our results support a model in which H-2(q) molecules convey Ly49-dependent inhibitory signals that interfere with the action of H2-D(k) on NK cell activation against MCMV infection. Thus, the integration of activating and inhibitory signals emanating from various MHC-I/NK cell receptor interactions regulates NK cell-mediated control of viral load.

New systems genetics approaches are needed to rapidly identify host genes and genetic networks that regulate complex disease outcomes. Using genetically diverse animals from incipient lines of the Collaborative Cross mouse panel, we demonstrate a greatly expanded range of phenotypes relative to classical mouse models of SARS-CoV infection including lung pathology, weight loss and viral titer. Genetic mapping revealed several loci contributing to differential disease responses, including an 8.5Mb locus associated with vascular cuffing on chromosome 3 that contained 23 genes and 13 noncoding RNAs. Integrating phenotypic and genetic data narrowed this region to a single gene, Trim55, an E3 ubiquitin ligase with a role in muscle fiber maintenance. Lung pathology and transcriptomic data from mice genetically deficient in Trim55 were used to validate its role in SARS-CoV-induced vascular cuffing and inflammation. These data establish the Collaborative Cross platform as a powerful genetic resource for uncovering genetic contributions of complex traits in microbial disease severity, inflammation and virus replication in models of outbred populations.

Atrio-ventricular conduction disease is a common feature in Mendelian rhythm disorders associated with sudden cardiac death and is characterized by prolongation of the PR interval on the surface electrocardiogram (ECG). Prolongation of the PR interval is also a strong predictor of atrial fibrillation, the most prevalent sustained cardiac arrhythmia. Despite the significant genetic component in PR duration variability, the genes regulating PR interval duration remain largely elusive. We here aimed to dissect the quantitative trait locus (QTL) for PR interval duration that we previously mapped in murine F2 progeny of a sensitized 129P2 and FVBN/J cross. To determine the underlying gene responsible for this QTL, genome-wide transcriptional profiling was carried out on myocardial tissue from 109 F2 mice. Expression QTLs (eQTLs) were mapped and the PR interval QTL was inspected for the co-incidence of eQTLs. We further determined the correlation of each of these transcripts to the PR interval. Tnni3k was the only eQTL, mapping to the PR-QTL, with an established abundant cardiac-specific expression pattern and a significant correlation to PR interval duration. Genotype inspection in various inbred mouse strains revealed the presence of at least three independent haplotypes at the Tnni3k locus. Measurement of PR interval duration and Tnni3k mRNA expression levels in six inbred lines identified a positive correlation between the level of Tnni3k mRNA and PR interval duration. Furthermore, in DBA/2J mice overexpressing hTNNI3K, and in DBA.AKR.hrtfm2 congenic mice, which harbor the AKR/J "high-Tnni3k expression" haplotype in the DBA/2J genetic background, PR interval duration was prolonged as compared to DBA/2J wild-type mice ("low-Tnni3k expression" haplotype). Our data provide the first evidence for a role of Tnni3k in controlling the electrocardiographic PR interval indicating a function of Tnni3k in atrio-ventricular conduction.

In the United States, roughly 10% of the population is exposed daily to hazardous levels of noise in the workplace. Twin studies estimate heritability for noise-induced hearing loss (NIHL) of approximately 36%, and strain specific variation in sensitivity has been demonstrated in mice. Based upon the difficulties inherent to the study of NIHL in humans, we have turned to the study of this complex trait in mice. We exposed 5 week-old mice from the Hybrid Mouse Diversity Panel (HMDP) to a 10 kHz octave band noise at 108 dB for 2 hours and assessed the permanent threshold shift 2 weeks post exposure using frequency specific stimuli. These data were then used in a genome-wide association study (GWAS) using the Efficient Mixed Model Analysis (EMMA) to control for population structure. In this manuscript we describe our GWAS, with an emphasis on a significant peak for susceptibility to NIHL on chromosome 17 within a haplotype block containing NADPH oxidase-3 (Nox3). Our peak was detected after an 8 kHz tone burst stimulus. Nox3 mutants and heterozygotes were then tested to validate our GWAS. The mutants and heterozygotes demonstrated a greater susceptibility to NIHL specifically at 8 kHz both on measures of distortion product otoacoustic emissions (DPOAE) and on auditory brainstem response (ABR). We demonstrate that this sensitivity resides within the synaptic ribbons of the cochlea in the mutant animals specifically at 8 kHz. Our work is the first GWAS for NIHL in mice and elucidates the power of our approach to identify tonotopic genetic susceptibility to NIHL.

During ischemic stroke, occlusion of the cerebrovasculature causes neuronal cell death (infarction), but naturally occurring genetic factors modulating infarction have been difficult to identify in human populations. In a surgically induced mouse model of ischemic stroke, we have previously mapped Civq1 to distal chromosome 7 as a quantitative trait locus determining infarct volume. In this study, genome-wide association mapping using 32 inbred mouse strains and an additional linkage scan for infarct volume confirmed that the size of the infarct is determined by ancestral alleles of the causative gene(s). The genetically isolated Civq1 locus in reciprocal recombinant congenic mice refined the critical interval and demonstrated that infarct size is determined by both vascular (collateral vessel anatomy) and non-vascular (neuroprotection) effects. Through the use of interval-specific SNP haplotype analysis, we further refined the Civq1 locus and identified integrin alpha L (Itgal) as one of the causative genes for Civq1. Itgal is the only gene that exhibits both strain-specific amino acid substitutions and expression differences. Coding SNPs, a 5-bp insertion in exon 30b, and increased mRNA and protein expression of a splice variant of the gene (Itgal-003, ENSMUST00000120857), all segregate with infarct volume. Mice lacking Itgal show increased neuronal cell death in both ex vivo brain slice and in vivo focal cerebral ischemia. Our data demonstrate that sequence variation in Itgal modulates ischemic brain injury, and that infarct volume is determined by both vascular and non-vascular mechanisms.

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.

Genetic variation in the major histocompatibility complex (MHC) affects CD4∶CD8 lineage commitment and MHC expression. However, the contribution of specific genes in this gene-dense region has not yet been resolved. Nor has it been established whether the same genes regulate MHC expression and T cell selection. Here, we assessed the impact of natural genetic variation on MHC expression and CD4∶CD8 lineage commitment using two genetic models in the rat. First, we mapped Quantitative Trait Loci (QTLs) associated with variation in MHC class I and II protein expression and the CD4∶CD8 T cell ratio in outbred Heterogeneous Stock rats. We identified 10 QTLs across the genome and found that QTLs for the individual traits colocalized within a region spanning the MHC. To identify the genes underlying these overlapping QTLs, we generated a large panel of MHC-recombinant congenic strains, and refined the QTLs to two adjacent intervals of ∼0.25 Mb in the MHC-I and II regions, respectively. An interaction between these intervals affected MHC class I expression as well as negative selection and lineage commitment of CD8 single-positive (SP) thymocytes. We mapped this effect to the transporter associated with antigen processing 2 (Tap2) in the MHC-II region and the classical MHC class I gene(s) (RT1-A) in the MHC-I region. This interaction was revealed by a recombination between RT1-A and Tap2, which occurred in 0.2% of the rats. Variants of Tap2 have previously been shown to influence the antigenicity of MHC class I molecules by altering the MHC class I ligandome. Our results show that a restricted peptide repertoire on MHC class I molecules leads to reduced negative selection of CD8SP cells. To our knowledge, this is the first study showing how a recombination between natural alleles of genes in the MHC influences lineage commitment of T cells.

Experimental autoimmune orchitis (EAO), the principal model of non-infectious testicular inflammatory disease, can be induced in susceptible mouse strains by immunization with autologous testicular homogenate and appropriate adjuvants. As previously established, the genome of DBA/2J mice encodes genes that are capable of conferring dominant resistance to EAO, while the genome of BALB/cByJ mice does not and they are therefore susceptible to EAO. In a genome scan, we previously identified Orch3 as the major quantitative trait locus controlling dominant resistance to EAO and mapped it to chromosome 11. Here, by utilizing a forward genetic approach, we identified kinesin family member 1C (Kif1c) as a positional candidate for Orch3 and, using a transgenic approach, demonstrated that Kif1c is Orch3. Mechanistically, we showed that the resistant Kif1c(D2) allele leads to a reduced antigen-specific T cell proliferative response as a consequence of decreased MHC class II expression by antigen presenting cells, and that the L(578) → P(578) and S(1027) → P(1027) polymorphisms distinguishing the BALB/cByJ and DBA/2J alleles, respectively, can play a role in transcriptional regulation. These findings may provide mechanistic insight into how polymorphism in other kinesins such as KIF21B and KIF5A influence susceptibility and resistance to human autoimmune diseases.

Osteoporosis is a genetic disease characterized by progressive reductions in bone mineral density (BMD) leading to an increased risk of fracture. Over the last decade, genome-wide association studies (GWASs) have identified over 1000 associations for BMD. However, as a phenotype BMD is challenging as bone is a multicellular tissue affected by both local and systemic physiology. Here, we focused on a single component of BMD, osteoblast-mediated bone formation in mice, and identified associations influencing osteoblast activity on mouse Chromosomes (Chrs) 1, 4, and 17. The locus on Chr. 4 was in an intergenic region between Wnt4 and Zbtb40, homologous to a locus for BMD in humans. We tested both Wnt4 and Zbtb40 for a role in osteoblast activity and BMD. Knockdown of Zbtb40, but not Wnt4, in osteoblasts drastically reduced mineralization. Additionally, loss-of-function mouse models for both genes exhibited reduced BMD. Our results highlight that investigating the genetic basis of in vitro osteoblast mineralization can be used to identify genes impacting bone formation and BMD.

Danforth's short tail mutant (Sd) mouse, first described in 1930, is a classic spontaneous mutant exhibiting defects of the axial skeleton, hindgut, and urogenital system. We used meiotic mapping in 1,497 segregants to localize the mutation to a 42.8-kb intergenic segment on chromosome 2. Resequencing of this region identified an 8.5-kb early retrotransposon (ETn) insertion within the highly conserved regulatory sequences upstream of Pancreas Specific Transcription Factor, 1a (Ptf1a). This mutation resulted in up to tenfold increased expression of Ptf1a as compared to wild-type embryos at E9.5 but no detectable changes in the expression levels of other neighboring genes. At E9.5, Sd mutants exhibit ectopic Ptf1a expression in embryonic progenitors of every organ that will manifest a developmental defect: the notochord, the hindgut, and the mesonephric ducts. Moreover, at E 8.5, Sd mutant mice exhibit ectopic Ptf1a expression in the lateral plate mesoderm, tail bud mesenchyme, and in the notochord, preceding the onset of visible defects such as notochord degeneration. The Sd heterozygote phenotype was not ameliorated by Ptf1a haploinsufficiency, further suggesting that the developmental defects result from ectopic expression of Ptf1a. These data identify disruption of the spatio-temporal pattern of Ptf1a expression as the unifying mechanism underlying the multiple congenital defects in Danforth's short tail mouse. This striking example of an enhancer mutation resulting in profound developmental defects suggests that disruption of conserved regulatory elements may also contribute to human malformation syndromes.

Hybrid sterility (HS) belongs to reproductive isolation barriers that safeguard the integrity of species in statu nascendi. Although hybrid sterility occurs almost universally among animal and plant species, most of our current knowledge comes from the classical genetic studies on Drosophila interspecific crosses or introgressions. With the house mouse subspecies Mus m. musculus and Mus m. domesticus as a model, new research tools have become available for studies of the molecular mechanisms and genetic networks underlying HS. Here we used QTL analysis and intersubspecific chromosome substitution strains to identify a 4.7 Mb critical region on Chromosome X (Chr X) harboring the Hstx2 HS locus, which causes asymmetrical spermatogenic arrest in reciprocal intersubspecific F1 hybrids. Subsequently, we mapped autosomal loci on Chrs 3, 9 and 13 that can abolish this asymmetry. Combination of immunofluorescent visualization of the proteins of synaptonemal complexes with whole-chromosome DNA FISH on pachytene spreads revealed that heterosubspecific, unlike consubspecific, homologous chromosomes are predisposed to asynapsis in F1 hybrid male and female meiosis. The asynapsis is under the trans- control of Hstx2 and Hst1/Prdm9 hybrid sterility genes in pachynemas of male but not female hybrids. The finding concurred with the fertility of intersubpecific F1 hybrid females homozygous for the Hstx2(Mmm) allele and resolved the apparent conflict with the dominance theory of Haldane's rule. We propose that meiotic asynapsis in intersubspecific hybrids is a consequence of cis-acting mismatch between homologous chromosomes modulated by the trans-acting Hstx2 and Prdm9 hybrid male sterility genes.

Understanding mammalian genetic systems is predicated on the determination of the relationship between genetic variation and phenotype. Several international programmes are under way to deliver mutations in every gene in the mouse genome. The challenge for mouse geneticists is to develop approaches that will provide comprehensive phenotype datasets for these mouse mutant libraries. Several factors are critical to success in this endeavour. It will be important to catalogue assay and environment and where possible to adopt standardised procedures for phenotyping tests along with common environmental conditions to ensure comparable datasets of phenotypes. Moreover, the scale of the task underlines the need to invest in technological development improving both the speed and cost of phenotyping platforms. In addition, it will be necessary to develop new informatics standards that capture the phenotype assay as well as other factors, genetic and environmental, that impinge upon phenotype outcome.