Ryan J. Buus

Subspecific origin and haplotype diversity in the laboratory mouse

Here we provide a genome-wide, high-resolution map of the phylogenetic origin of the genome of most extant laboratory mouse inbred strains. Our analysis is based on the genotypes of wild-caught mice from three subspecies of Mus musculus. We show that classical laboratory strains are derived from a few fancy mice with limited haplotype diversity. Their genomes are overwhelmingly Mus musculus domesticus in origin, and the remainder is mostly of Japanese origin. We generated genome-wide haplotype maps based on identity by descent from fancy mice and show that classical inbred strains have limited and non-randomly distributed genetic diversity. In contrast, wild-derived laboratory strains represent a broad sampling of diversity within M. musculus. Intersubspecific introgression is pervasive in these strains, and contamination by laboratory stocks has played a role in this process. The subspecific origin, haplotype diversity and identity by descent maps can be visualized using the Mouse Phylogeny Viewer (see URLs).

Genetic analysis of complex traits in the emerging Collaborative Cross

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.

Modeling Host Genetic Regulation of Influenza Pathogenesis in the Collaborative Cross

Genetic variation contributes to host responses and outcomes following infection by influenza A virus or other viral infections. Yet narrow windows of disease symptoms and confounding environmental factors have made it difficult to identify polymorphic genes that contribute to differential disease outcomes in human populations. Therefore, to control for these confounding environmental variables in a system that models the levels of genetic diversity found in outbred populations such as humans, we used incipient lines of the highly genetically diverse Collaborative Cross (CC) recombinant inbred (RI) panel (the pre-CC population) to study how genetic variation impacts influenza associated disease across a genetically diverse population. A wide range of variation in influenza disease related phenotypes including virus replication, virus-induced inflammation, and weight loss was observed. Many of the disease associated phenotypes were correlated, with viral replication and virus-induced inflammation being predictors of virus-induced weight loss. Despite these correlations, pre-CC mice with unique and novel disease phenotype combinations were observed. We also identified sets of transcripts (modules) that were correlated with aspects of disease. In order to identify how host genetic polymorphisms contribute to the observed variation in disease, we conducted quantitative trait loci (QTL) mapping. We identified several QTL contributing to specific aspects of the host response including virus-induced weight loss, titer, pulmonary edema, neutrophil recruitment to the airways, and transcriptional expression. Existing whole-genome sequence data was applied to identify high priority candidate genes within QTL regions. A key host response QTL was located at the site of the known anti-influenza Mx1 gene. We sequenced the coding regions of Mx1 in the eight CC founder strains, and identified a novel Mx1 allele that showed reduced ability to inhibit viral replication, while maintaining protection from weight loss.

Analyses of allele-specific gene expression in highly divergent mouse crosses identifies pervasive allelic imbalance

Complex human traits are influenced by variation in regulatory DNA through mechanisms that are not fully understood. Because regulatory elements are conserved between humans and mice, a thorough annotation of cis regulatory variants in mice could aid in further characterizing these mechanisms. Here we provide a detailed portrait of mouse gene expression across multiple tissues in a three-way diallel. Greater than 80% of mouse genes have cis regulatory variation. Effects from these variants influence complex traits and usually extend to the human ortholog. Further, we estimate that at least one in every thousand SNPs creates a cis regulatory effect. We also observe two types of parent-of-origin effects, including classical imprinting and a new global allelic imbalance in expression favoring the paternal allele. We conclude that, as with humans, pervasive regulatory variation influences complex genetic traits in mice and provide a new resource toward understanding the genetic control of transcription in mammals.

Nerve growth factor sequestering therapy attenuates non-malignant skeletal pain following fracture

Abstract Current therapies to treat skeletal fracture pain are extremely limited. Some non‐steroidal anti‐inflammatory drugs have been shown to inhibit bone healing and opiates induce cognitive dysfunction and respiratory depression which are especially problematic in the elderly suffering from osteoporotic fractures. In the present report, we developed a closed femur fracture pain model in the mouse where skeletal pain behaviors such as flinching and guarding of the fractured limb are reversed by 10 mg/kg morphine. Using this model we showed that the administration of a monoclonal antibody against nerve growth factor (anti‐NGF) reduced fracture‐induced pain‐related behaviors by over 50%. Treatment with anti‐NGF reduced c‐Fos and dynorphin up‐regulation in the spinal cord at day 2 post‐fracture. However, anti‐NGF treatment did not reduce p‐ERK and c‐Fos expression at 20 and 90 min, respectively, following fracture. This suggests NGF is involved in maintenance but not the acute generation of fracture pain. Anti‐NGF therapy did not inhibit bone healing as measured by callus formation, bridging of the fracture site or mechanical strength of the bone. As the anti‐NGF antibody does not appreciably cross the blood–brain barrier, the present data suggest that the anti‐hyperalgesic action of anti‐NGF therapy results from blockade of activation and/or sensitization of the CGRP/trkA positive fibers that normally constitute the majority of sensory fibers that innervate the bone. These results demonstrate that NGF plays a significant role in driving fracture pain and that NGF sequestering therapies may be efficacious in attenuating this pain.

Effects of a Monoclonal Antibody Raised Against Nerve Growth Factor on Skeletal Pain and Bone Healing After Fracture of the C57BL/6J Mouse Femur†

A closed femur fracture pain model was developed in the C57BL/6J mouse. One day after fracture, a monoclonal antibody raised against nerve growth factor (anti‐NGF) was delivered intraperitoneally and resulted in a reduction in fracture pain‐related behaviors of ∼50%. Anti‐NGF therapy did not interfere with bone healing as assessed by mechanical testing and histomorphometric analysis.

Genetic Analysis of Hematological Parameters in Incipient Lines of the Collaborative Cross

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.

Genome Wide Identification of SARS-CoV Susceptibility Loci Using the Collaborative Cross.

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.

A Novel Intronic Single Nucleotide Polymorphism in the Myosin heavy polypeptide 4 Gene Is Responsible for the Mini-Muscle Phenotype Characterized by Major Reduction in Hind-Limb Muscle Mass in Mice

Replicated artificial selection for high levels of voluntary wheel running in an outbred strain of mice favored an autosomal recessive allele whose primary phenotypic effect is a 50% reduction in hind-limb muscle mass. Within the High Runner (HR) lines of mice, the numerous pleiotropic effects (e.g., larger hearts, reduced total body mass and fat mass, longer hind-limb bones) of this hypothesized adaptive allele include functional characteristics that facilitate high levels of voluntary wheel running (e.g., doubling of mass-specific muscle aerobic capacity, increased fatigue resistance of isolated muscles, longer hind-limb bones). Previously, we created a backcross population suitable for mapping the responsible locus. We phenotypically characterized the population and mapped the Minimsc locus to a 2.6-Mb interval on MMU11, a region containing ∼100 known or predicted genes. Here, we present a novel strategy to identify the genetic variant causing the mini-muscle phenotype. Using high-density genotyping and whole-genome sequencing of key backcross individuals and HR mice with and without the mini-muscle mutation, from both recent and historical generations of the HR lines, we show that a SNP representing a C-to-T transition located in a 709-bp intron between exons 11 and 12 of the Myosin heavy polypeptide 4 (Myh4) skeletal muscle gene (position 67,244,850 on MMU11; assembly, December 2011, GRCm38/mm10; ENSMUSG00000057003) is responsible for the mini-muscle phenotype, Myh4Minimsc. Using next-generation sequencing, our approach can be extended to identify causative mutations arising in mouse inbred lines and thus offers a great avenue to overcome one of the most challenging steps in quantitative genetics.

Identification of quantitative trait loci influencing skeletal architecture in mice: emergence of Cdh11 as a primary candidate gene regulating femoral morphology.

Bone strength is influenced by many properties intrinsic to bone, including its mass, geometry, and mineralization. To further advance our understanding of the genetic basis of bone‐strength‐related traits, we used a large (n = 815), moderately (G4) advanced intercross line (AIL) of mice derived from a high‐runner selection line (HR) and the C57BL/6J inbred strain. In total, 16 quantitative trait loci (QTLs) were identified that affected areal bone mineral density (aBMD) and femoral length and width. Four significant (p < .05) and one suggestive (p < .10) QTLs were identified for three aBMD measurements: total body, vertebral, and femoral. A QTL on chromosome (Chr.) 3 influenced all three aBMD measures, whereas the other four QTLs were unique to a single measure. A total of 10 significant and one suggestive QTLs were identified for femoral length (FL) and two measures of femoral width, anteroposterior (AP) and mediolateral (ML). FL QTLs were distinct from loci affecting AP and ML width, and of the 7 AP QTLs, only three affected ML. A QTL on Chr. 8 that explained 7.1% and 4.0% of the variance in AP and ML, respectively, was mapped to a 6‐Mb region harboring 12 protein‐coding genes. The pattern of haplotype diversity across the QTL region and expression profiles of QTL genes suggested that of the 12, cadherin 11 (Cdh11) was most likely the causal gene. These findings, when combined with existing data from gene knockouts, identify Cdh11 as a strong candidate gene within which genetic variation may affect bone morphology.