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. Since regulatory elements are conserved between humans and mice, a thorough annotation of cis regulatory variants in mice could aid in this process. 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. These effects 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 novel, global allelic imbalance in favor of 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

Genetic diversity between mouse strains allows identification of the CC027/GeniUnc strain as an orally reactive model of peanut allergy

Background: Improved animal models are needed to understand the genetic and environmental factors that contribute to food allergy. Objective: We sought to assess food allergy phenotypes in a genetically diverse collection of mice. Methods: We selected 16 Collaborative Cross (CC) mouse strains, as well as the classic inbred C57BL/6J, C3H/HeJ, and BALB/cJ strains, for screening. Female mice were sensitized to peanut intragastrically with or without cholera toxin and then challenged with peanut by means of oral gavage or intraperitoneal injection and assessed for anaphylaxis. Peanut-specific immunoglobulins, T-cell cytokines, regulatory T cells, mast cells, and basophils were quantified. Results: Eleven of the 16 CC strains had allergic reactions to intraperitoneal peanut challenge, whereas only CC027/GeniUnc mice reproducibly experienced severe symptoms after oral food challenge (OFC). CC027/GeniUnc, C3H/HeJ, and C57BL/6J mice all mounted a T H 2 response against peanut, leading to production of IL-4 and IgE, but only the CC027/GeniUnc mice reacted to OFC. Orally induced anaphylaxis in CC027/GeniUnc mice was correlated with serum levels of Ara h 2 in circulation but not with allergen-specific IgE or mucosal mast cell protease 1 levels, indicating systemic allergen absorption is important for anaphylaxis through the gastrointestinal tract. Furthermore, CC027/GeniUnc, but not C3H/HeJ or BALB/cJ, mice can be sensitized in the absence of cholera toxin and react on OFC to peanut. Conclusions: We have identified and characterized CC027/GeniUnc mice as a strain that is genetically susceptible to peanut allergy and prone to severe reactions after OFC. More broadly, these findings demonstrate the untapped potential of the CC population in developing novel models for allergy research

Transcriptome Atlases of Mouse Brain Reveals Differential Expression Across Brain Regions and Genetic Backgrounds

Mouse models play a crucial role in the study of human behavioral traits and diseases. Variation of gene expression in brain may play a critical role in behavioral phenotypes, and thus it is of great importance to understand regulation of transcription in mouse brain. In this study, we analyzed the role of two important factors influencing steady-state transcriptional variation in mouse brain. First we considered the effect of assessing whole brain vs. discrete regions of the brain. Second, we investigated the genetic basis of strain effects on gene expression. We examined the transcriptome of three brain regions using Affymetrix expression arrays: whole brain, forebrain, and hindbrain in adult mice from two common inbred strains (C57BL/6J vs. NOD/ShiLtJ) with eight replicates for each brain region and strain combination. We observed significant differences between the transcriptomes of forebrain and hindbrain. In contrast, the transcriptomes of whole brain and forebrain were very similar. Using 4.3 million single-nucleotide polymorphisms identified through whole-genome sequencing of C57BL/6J and NOD/ShiLtJ strains, we investigated the relationship between strain effect in gene expression and DNA sequence similarity. We found that cis-regulatory effects play an important role in gene expression differences between strains and that the cis-regulatory elements are more often located in 5′ and/or 3′ transcript boundaries, with no apparent preference on either 5′ or 3′ ends

The Evolutionary Fates of a Large Segmental Duplication in Mouse

Gene duplication and loss are major sources of genetic polymorphism in populations, and are important forces shaping the evolution of genome content and organization. We have reconstructed the origin and history of a 127-kbp segmental duplication, R2d, in the house mouse (Mus musculus). R2d contains a single protein-coding gene, Cwc22. De novo assembly of both the ancestral (R2d1) and the derived (R2d2) copies reveals that they have been subject to nonallelic gene conversion events spanning tens of kilobases. R2d2 is also a hotspot for structural variation: its diploid copy number ranges from zero in the mouse reference genome to >80 in wild mice sampled from around the globe. Hemizygosity for high copy-number alleles of R2d2 is associated in cis with meiotic drive; suppression of meiotic crossovers; and copy-number instability, with a mutation rate in excess of 1 per 100 transmissions in some laboratory populations. Our results provide a striking example of allelic diversity generated by duplication and demonstrate the value of de novo assembly in a phylogenetic context for understanding the mutational processes affecting duplicate genes

Nondisjunction and transmission ratio distortion ofChromosome 2 in a (2.8) Robertsonian translocation mouse strain

Aneuploidy results from nondisjunction of chromosomes in meiosis and is the leading cause of developmental disabilities and mental retardation in humans. Therefore, understanding aspects of chromosome segregation in a genetic model is of value. Mice heterozygous for a (2.8) Robertsonian translocation were intercrossed with chromosomally normal mice and Chromosome 2 was genotyped for number and parental origin in 836 individuals at 8.5 dpc. The frequency of nondisjunction of this Robertsonian chromosome is 1.58%. Trisomy of Chromosome 2 with two maternally derived chromosomes is the most developmentally successful aneuploid karyotype at 8.5 dpc. Trisomy of Chromosome 2 with two paternally derived chromosomes is developmentally delayed and less frequent than the converse. Individuals with maternal or paternal uniparental disomy of Chromosome 2 were not detected at 8.5 dpc. Nondisjunction events were distributed randomly across litters, i.e., no evidence for clustering was found. Transmission ratio distortion is frequently observed in Robertsonian chromosomes and a bias against the transmission of the (2.8) Chromosome was detected. Interestingly, this was observed for female and male transmitting parents

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

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

IsoDOT Detects Differential RNA-isoform Expression/Usage with respect to a Categorical or Continuous Covariate with High Sensitivity and Specificity

We have developed a statistical method named IsoDOT to assess differential isoform expression (DIE) and differential isoform usage (DIU) using RNA-seq data. Here isoform usage refers to relative isoform expression given the total expression of the corresponding gene. IsoDOT performs two tasks that cannot be accomplished by existing methods: to test DIE/DIU with respect to a continuous covariate, and to test DIE/DIU for one case versus one control. The latter task is not an uncommon situation in practice, e.g., comparing paternal and maternal allele of one individual or comparing tumor and normal sample of one cancer patient. Simulation studies demonstrate the high sensitivity and specificity of IsoDOT. We apply IsoDOT to study the effects of haloperidol treatment on mouse transcriptome and identify a group of genes whose isoform usages respond to haloperidol treatment

Erratum to: Status and access to the Collaborative Cross population (Mammalian Genome (2012) 23, (706-712) DOI: 10.1007/s00335-012-9410-6)

The Collaborative Cross (CC) is a panel of recombinant inbred lines derived from eight genetically diverse laboratory inbred strains. Recently, the genetic architecture of the CC population was reported based on the genotype of a single male per line, and other publications reported incompletely inbred CC mice that have been used to map a variety of traits. The three breeding sites, in the US, Israel, and Australia, are actively collaborating to accelerate the inbreeding process through marker-assisted inbreeding and to expedite community access of CC lines deemed to have reached defined thresholds of inbreeding. Plans are now being developed to provide access to this novel genetic reference population through distribution centers. Here we provide a description of the distribution efforts by the University of North Carolina Systems Genetics Core, Tel Aviv University, Israel and the University of Western Australia

C57BL/6N Mutation in Cytoplasmic FMRP interacting protein 2 Regulates Cocaine Response

The inbred mouse C57BL/6J is the reference strain for genome sequence and for most behavioral and physiological phenotypes. However the International Knockout Mouse Consortium uses an embryonic stem cell line derived from a related C57BL/6N substrain. We found that C57BL/6N has lower acute and sensitized response to cocaine and methamphetamine. We mapped a single causative locus and identified a non-synonymous mutation of serine to phenylalanine (S968F) in Cytoplasmic FMR interacting protein 2 (Cyfip2) as the causative variant. The S968F mutation destabilizes CYFIP2 and deletion of the C57BL/6N mutant allele leads to acute and sensitized cocaine response phenotypes. We propose CYFIP2 is a key regulator of cocaine response in mammals and present a framework to utilize mouse substrains to discover novel genes and alleles regulating behavior