Next-generation sequencing of experimental mouse strains

Single-housing-induced islet epigenomic changes are related to polymorphisms in diabetic KK mice

A lack of social relationships is increasingly recognized as a type 2 diabetes (T2D) risk. To investigate the underlying mechanism, we used male KK mice, an inbred strain with spontaneous diabetes. Given the association between living alone and T2D risk in humans, we divided the non-diabetic mice into singly housed (KK-SH) and group-housed control mice. Around the onset of diabetes in KK-SH mice, we compared H3K27ac ChIP-Seq with RNA-Seq using pancreatic islets derived from each experimental group, revealing a positive correlation between single-housing-induced changes in H3K27ac and gene expression levels. In particular, single-housing-induced H3K27ac decreases revealed a significant association with islet cell functions and GWAS loci for T2D and related diseases, with significant enrichment of binding motifs for transcription factors representative of human diabetes. Although these H3K27ac regions were preferentially localized to a polymorphic genomic background, SNVs and indels did not cause sequence disruption of enriched transcription factor motifs in most of these elements. These results suggest alternative roles of genetic variants in environment-dependent epigenomic changes and provide insights into the complex mode of disease inheritance.

Excessive self-grooming, gene dysregulation and imbalance between the striosome and matrix compartments in the striatum of Shank3 mutant mice

Autism is characterized by atypical social communication and stereotyped behaviors. Mutations in the gene encoding the synaptic scaffolding protein SHANK3 are detected in 1–2% of patients with autism and intellectual disability, but the mechanisms underpinning the symptoms remain largely unknown. Here, we characterized the behavior of Shank3Δ11/Δ11 mice from 3 to 12 months of age. We observed decreased locomotor activity, increased stereotyped self-grooming and modification of socio-sexual interaction compared to wild-type littermates. We then used RNAseq on four brain regions of the same animals to identify differentially expressed genes (DEGs). DEGs were identified mainly in the striatum and were associated with synaptic transmission (e.g., Grm2, Dlgap1), G-protein-signaling pathways (e.g., Gnal, Prkcg1, and Camk2g), as well as excitation/inhibition balance (e.g., Gad2). Downregulated and upregulated genes were enriched in the gene clusters of medium-sized spiny neurons expressing the dopamine 1 (D1-MSN) and the dopamine 2 receptor (D2-MSN), respectively. Several DEGs (Cnr1, Gnal, Gad2, and Drd4) were reported as striosome markers. By studying the distribution of the glutamate decarboxylase GAD65, encoded by Gad2, we showed that the striosome compartment of Shank3Δ11/Δ11 mice was enlarged and displayed much higher expression of GAD65 compared to wild-type mice. Altogether, these results indicate altered gene expression in the striatum of Shank3-deficient mice and strongly suggest, for the first time, that the excessive self-grooming of these mice is related to an imbalance in the striatal striosome and matrix compartments.

Genetic variant in 3' untranslated region of the mouse pycard gene regulates inflammasome activity

Quantitative trait locus mapping for interleukin-1β release after inflammasome priming and activation was performed on bone-marrow-derived macrophages (BMDM) from an AKRxDBA/2 mouse strain intercross. The strongest associated locus mapped very close to the Pycard gene on chromosome 7, which codes for the inflammasome adaptor protein apoptosis-associated speck-like protein containing a CARD (ASC). The DBA/2 and AKR Pycard genes only differ at a single-nucleotide polymorphism (SNP) in their 3' untranslated region (UTR). DBA/2 vs. AKR BMDM had increased levels of Pycard mRNA expression and ASC protein, and increased inflammasome speck formation, which was associated with increased Pycard mRNA stability without an increased transcription rate. CRISPR/Cas9 gene editing was performed on DBA/2 embryonic stem cells to change the Pycard 3'UTR SNP from the DBA/2 to the AKR allele. This single base change significantly reduced Pycard expression and inflammasome activity after cells were differentiated into macrophages due to reduced Pycard mRNA stability.

High-throughput genotyping of high-homology mutant mouse strains by next-generation sequencing

Genotyping of knockout alleles in mice is commonly performed by end-point PCR or gene-specific/universal cassette qPCR. Both have advantages and limitations in terms of assay design and interpretation of results. As an alternative method for high-throughput genotyping, we investigated next generation sequencing (NGS) of PCR amplicons, with a focus on CRISPR-mediated exon deletions where antibiotic selection markers are not present. By multiplexing the wild type and mutant-specific PCR reactions, the genotype can be called by the relative sequence counts of each product. The system is highly scalable and can be applied to a variety of different allele types, including those produced by the International Mouse Phenotyping Consortium and associated projects. One potential challenge with any assay design is locating unique areas of the genome, especially when working with gene families or regions of high homology. These can result in misleading or ambiguous genotypes for either qPCR or end-point assays. Here, we show that genotyping by NGS can negate these issues by simple, automated filtering of undesired sequences. Analysis and genotype calls can also be fully automated, using FASTQ or FASTA input files and an in-house Perl script and SQL database.

Diploid genome architecture revealed by multi-omic data of hybrid mice

Although mammalian genomes are diploid, previous studies extensively investigated the average chromatin architectures without considering the differences between homologous chromosomes. We generated Hi-C, ChIP-seq, and RNA-seq data sets from CD4 T cells of B6, Cast, and hybrid mice, to investigate the diploid chromatin organization and epigenetic regulation. Our data indicate that inter-chromosomal interaction patterns between homologous chromosomes are similar, and the similarity is highly correlated with their allelic coexpression levels. Reconstruction of the 3D nucleus revealed that distances of the homologous chromosomes to the center of nucleus are almost the same. The inter-chromosomal interactions at centromere ends are significantly weaker than those at telomere ends, suggesting that they are located in different regions within the chromosome territories. The majority of A|B compartments or topologically associated domains (TADs) are consistent between B6 and Cast. We found 58% of the haploids in hybrids maintain their parental compartment status at B6/Cast divergent compartments owing to cis effect. About 95% of the trans-effected B6/Cast divergent compartments converge to the same compartment status potentially because of a shared cellular environment. We showed the differentially expressed genes between the two haploids in hybrid were associated with either genetic or epigenetic effects. In summary, our multi-omics data from the hybrid mice provided haploid-specific information on the 3D nuclear architecture and a rich resource for further understanding the epigenetic regulation of haploid-specific gene expression.

Transcriptional activity and strain-specific history of mouse pseudogenes

Pseudogenes are ideal markers of genome remodelling. In turn, the mouse is an ideal platform for studying them, particularly with the recent availability of strain-sequencing and transcriptional data. Here, combining both manual curation and automatic pipelines, we present a genome-wide annotation of the pseudogenes in the mouse reference genome and 18 inbred mouse strains (available via the mouse.pseudogene.org resource). We also annotate 165 unitary pseudogenes in mouse, and 303, in human. The overall pseudogene repertoire in mouse is similar to that in human in terms of size, biotype distribution, and family composition (e.g. with GAPDH and ribosomal proteins being the largest families). Notable differences arise in the pseudogene age distribution, with multiple retro-transpositional bursts in mouse evolutionary history and only one in human. Furthermore, in each strain about a fifth of all pseudogenes are unique, reflecting strain-specific evolution. Finally, we find that ~15% of the mouse pseudogenes are transcribed, and that highly transcribed parent genes tend to give rise to many processed pseudogenes.

Evolutionary analysis of genes coding for Cysteine-RIch Secretory Proteins (CRISPs) in mammals

BACKGROUND

Cysteine-RIch Secretory Proteins (CRISP) are expressed in the reproductive tract of mammalian males and are involved in fertilization and related processes. Due to their important role in sperm performance and sperm-egg interaction, these genes are likely to be exposed to strong selective pressures, including postcopulatory sexual selection and/or male-female coevolution. We here perform a comparative evolutionary analysis of Crisp genes in mammals. Currently, the nomenclature of CRISP genes is confusing, as a consequence of discrepancies between assignments of orthologs, particularly due to numbering of CRISP genes. This may generate problems when performing comparative evolutionary analyses of mammalian clades and species. To avoid such problems, we first carried out a study of possible orthologous relationships and putative origins of the known CRISP gene sequences. Furthermore, and with the aim to facilitate analyses, we here propose a different nomenclature for CRISP genes (EVAC1-4, "EVolutionarily-analyzed CRISP") to be used in an evolutionary context.

RESULTS

We found differing selective pressures among Crisp genes. CRISP1/4 (EVAC1) and CRISP2 (EVAC2) orthologs are found across eutherian mammals and seem to be conserved in general, but show signs of positive selection in primate CRISP1/4 (EVAC1). Rodent Crisp1 (Evac3a) seems to evolve under a comparatively more relaxed constraint with positive selection on codon sites. Finally, murine Crisp3 (Evac4), which appears to be specific to the genus Mus, shows signs of possible positive selection. We further provide evidence for sexual selection on the sequence of one of these genes (Crisp1/4) that, unlike others, is thought to be exclusively expressed in male reproductive tissues.

CONCLUSIONS

We found differing selective pressures among CRISP genes and sexual selection as a contributing factor in CRISP1/4 gene sequence evolution. Our evolutionary analysis of this unique set of genes contributes to a better understanding of Crisp function in particular and the influence of sexual selection on reproductive mechanisms in general.

De novo leaf transcriptome assembly of Bougainvillea spectabilis for the identification of genes involves in the secondary metabolite pathways

Bougainvillea spectabilis is known as a vital medicinal, ornamental as well as an essential oil producing plant. It is also a rich source of important secondary metabolites with several therapeutic properties. Various studies on its pharmacological and toxicological aspects have been published but there is no genomic or transcriptomic resource available in the public databases. To address this important issue, the de-novo transcriptome assembly of B. spectabilis leaf tissue has been done for the identification of genes involved in various important secondary metabolites, Single nucleotide polymorphism (SNPs) and Simple sequence repeats (SSRs). The transcriptome sequencing of B. spectabilis leaf tissue generated 79,811,024 raw reads with GC value 42.77%. The transcriptomic assembly was performed by Trinity software which generated 100,374 transcripts and 99,793 unigenes with minimum and maximum length of 201 bp and 13,237 bp and N50 value of 1470 and 1472 respectively. Annotation of these unigenes was performed using seven databases including NR, PFAM, GO and KEGG. Approximately, 44,302 unigenes were annotated in GO database. The KEGG pathway analysis revealed 23,102 unigenes in which 19,054 genes were assigned to five groups in KEGG and 130 biochemical pathways. The highest group among the five groups was Metabolism with 9230 unigenes. Moreover, about 63,226 SNPs and 30,333 SSRs in the leaf transcriptome of B. spectabilis were identified. To the best of my understanding it will be the first comprehensive transcriptome analysis of B. spectabilis from family Nyctaginaceae which will help as a reference line for further genomic and transcriptomic studies.

Comparative genomic analysis of inbred rat strains reveals the existence of ancestral polymorphisms

In an alignment of closely related genomic sequences, the existence of discordant mutation sites, which do not reflect the phylogenetic relationship of the genomes, is often observed. Although these discordant mutation sites are thought to have emerged by ancestral polymorphism or gene flow, their frequency and distribution in the genome have not yet been analyzed in detail. Using the genome sequences of all protein coding genes of 25 inbred rat strains, we analyzed the frequency and genome-wide distribution of the discordant mutation sites. From the comparison of different substrains, it was found that these loci are not substrain specific, but are common among different groups of substrains, suggesting that the discordant sites might have mainly emerged through ancestral polymorphism. It was also revealed that the discordant sites are not uniformly distributed along chromosomes, but are concentrated at certain genomic loci, such as RT1, major histocompatibility complex of rats, and olfactory receptors, indicating that genes known to be highly polymorphic tend to have more discordant sites. Our results also showed that loci with a high density of discordant sites are also rich in heterozygous variants, even though these are inbred strains.

A Microbe Associated with Sleep Revealed by a Novel Systems Genetic Analysis of the Microbiome in Collaborative Cross Mice

The microbiome influences health and disease through complex networks of host genetics, genomics, microbes, and environment. Identifying the mechanisms of these interactions has remained challenging. Systems genetics in laboratory mice (Mus musculus) enables data-driven discovery of biological network components and mechanisms of host-microbial interactions underlying disease phenotypes. To examine the interplay among the whole host genome, transcriptome, and microbiome, we mapped QTL and correlated the abundance of cecal messenger RNA, luminal microflora, physiology, and behavior in a highly diverse Collaborative Cross breeding population. One such relationship, regulated by a variant on chromosome 7, was the association of Odoribacter (Bacteroidales) abundance and sleep phenotypes. In a test of this association in the BKS.Cg-Dock7m +/+ Leprdb/J mouse model of obesity and diabetes, known to have abnormal sleep and colonization by Odoribacter, treatment with antibiotics altered sleep in a genotype-dependent fashion. The many other relationships extracted from this study can be used to interrogate other diseases, microbes, and mechanisms.

Microarray-Based Quality Assessment as a Supporting Criterion for de novo Transcriptome Assembly Selection

RNA-Sequencing and de novo assembly have enabled the analysis of species with non-available reference transcriptomes, although intrinsic features (biological and technical) induce errors in the reconstruction. A strategy to resolve these errors consists of varying assembling process parameters to generate multiple reconstructions. However, the best assembly selection remains a challenge. Quantitative metrics for quality assessment have been inconsistent when compared with pertinent references. In this paper, a criterion for supporting assembly selection based on mapping DNA microarray hybridized probes to assembly sets is proposed. Mouse and fruit fly RNA-Seq datasets were assembled with standard de novo procedures. Quality assessment was estimated using quantitative metrics and the proposed criterion. The assembly that best mapped to the available reference transcriptomes of these model species provided the highest quality assembly. The hybridized probes identified the best assemblies, whereas quantitative metrics remained inconsistent. For example, subtle probe mapping difference of 0.25 percent, but statistically significant (ANOVA, p < 0.05), enabled the assembly selection that led to identify 3,719 more contigs and led to 1,049 further mapped contigs to the mouse reference transcriptome. The microarray data availability for non-model species makes the proposed criterion suitable for quality assessment of multiple de novo assembly strategies.

Characterization of genetic subclonal evolution in pancreatic cancer mouse models

The KPC mouse model, driven by the Kras and Trp53 transgenes, is well regarded for faithful recapitulation of human pancreatic cancer biology. However, the extent that this model recapitulates the subclonal evolution of this tumor type is unknown. Here we report evidence of continuing subclonal evolution after tumor initiation that largely reflect copy number alterations that target cellular processes of established significance in human pancreatic cancer. The evolutionary trajectories of the mouse tumors show both linear and branching patterns as well as clonal mixing. We propose the KPC model and derivatives have unexplored utility as a functional system to model the mechanisms and modifiers of tumor evolution.

In pancreatic cancer the Kras and Trp53 transgene driven KPC mouse model is used to experimentally study disease processes. Here, the authors analyse tumour evolution within the KPC model, finding both linear and branched evolution and highlighting the utility of this model in mechanistic research of tumour evolution.

Impact of mouse contamination in genomic profiling of patient-derived models and best practice for robust analysis

BACKGROUND

Patient-derived xenograft and cell line models are popular models for clinical cancer research. However, the inevitable inclusion of a mouse genome in a patient-derived model is a remaining concern in the analysis. Although multiple tools and filtering strategies have been developed to account for this, research has yet to demonstrate the exact impact of the mouse genome and the optimal use of these tools and filtering strategies in an analysis pipeline.

RESULTS

We construct a benchmark dataset of 5 liver tissues from 3 mouse strains using human whole-exome sequencing kit. Next-generation sequencing reads from mouse tissues are mappable to 49% of the human genome and 409 cancer genes. In total, 1,207,556 mouse-specific alleles are aligned to the human genome reference, including 467,232 (38.7%) alleles with high sensitivity to contamination, which are pervasive causes of false cancer mutations in public databases and are signatures for predicting global contamination. Next, we assess the performance of 8 filtering methods in terms of mouse read filtration and reduction of mouse-specific alleles. All filtering tools generally perform well, although differences in algorithm strictness and efficiency of mouse allele removal are observed. Therefore, we develop a best practice pipeline that contains the estimation of contamination level, mouse read filtration, and variant filtration.

CONCLUSIONS

The inclusion of mouse cells in patient-derived models hinders genomic analysis and should be addressed carefully. Our suggested guidelines improve the robustness and maximize the utility of genomic analysis of these models.

Whole Genome Sequencing of the Mutamouse Model Reveals Strain- and Colony-Level Variation, and Genomic Features of the Transgene Integration Site

The MutaMouse transgenic rodent model is widely used for assessing in vivo mutagenicity. Here, we report the characterization of MutaMouse's whole genome sequence and its genetic variants compared to the C57BL/6 reference genome. High coverage (>50X) next-generation sequencing (NGS) of whole genomes from multiple MutaMouse animals from the Health Canada (HC) colony showed ~5 million SNVs per genome, ~20% of which are putatively novel. Sequencing of two animals from a geographically separated colony at Covance indicated that, over the course of 23 years, each colony accumulated 47,847 (HC) and 17,677 (Covance) non-parental homozygous single nucleotide variants. We found no novel nonsense or missense mutations that impair the MutaMouse response to genotoxic agents. Pairing sequencing data with array comparative genomic hybridization (aCGH) improved the accuracy and resolution of copy number variants (CNVs) calls and identified 300 genomic regions with CNVs. We also used long-read sequence technology (PacBio) to show that the transgene integration site involved a large deletion event with multiple inversions and rearrangements near a retrotransposon. The MutaMouse genome gives important genetic context to studies using this model, offers insight on the mechanisms of structural variant formation, and contributes a framework to analyze aCGH results alongside NGS data.

Tissue-Specific Trans Regulation of the Mouse Epigenome

The epigenetic landscape varies greatly among cell types. Although a variety of writers, readers, and erasers of epigenetic features are known, we have little information about the underlying regulatory systems controlling the establishment and maintenance of these features. Here, we have explored how natural genetic variation affects the epigenome in mice. Studying levels of H3K4me3, a histone modification at sites such as promoters, enhancers, and recombination hotspots, we found tissue-specific trans-regulation of H3K4me3 levels in four highly diverse cell types: male germ cells, embryonic stem cells, hepatocytes, and cardiomyocytes. To identify the genetic loci involved, we measured H3K4me3 levels in male germ cells in a mapping population of 59 BXD recombinant inbred lines. We found extensive trans-regulation of H3K4me3 peaks, including six major histone quantitative trait loci (QTL). These chromatin regulatory loci act dominantly to suppress H3K4me3, which at hotspots reduces the likelihood of subsequent DNA double-strand breaks. QTL locations do not correspond with genes encoding enzymes known to metabolize chromatin features. Instead their locations match clusters of zinc finger genes, making these possible candidates that explain the dominant suppression of H3K4me3. Collectively, these data describe an extensive, set of chromatin regulatory loci that control the epigenetic landscape.

Murine Models for Viral Hemorrhagic Fever

Hemorrhagic fever (HF) viruses, such as Lassa, Ebola, and dengue viruses, represent major human health risks due to their highly contagious nature, the severity of the clinical manifestations induced, the lack of vaccines, and the very limited therapeutic options currently available. Appropriate animal models are obviously critical to study disease pathogenesis and develop efficient therapies. We recently reported that the clone 13 (Cl13) variant of the lymphocytic choriomeningitis virus (LCMV-Cl13), a prototype arenavirus closely related to Lassa virus, causes in some mouse strains endothelial damage, vascular leakage, platelet loss, and death, mimicking pathological aspects typically observed in Lassa and other HF syndromes. This model has the advantage that the mice used are fully immunocompetent, allowing studies on the contribution of the immune response to disease progression. Moreover, LCMV is very well characterized and exhibits limited pathogenicity in humans, allowing handling in convenient BSL-2 facilities. In this chapter we outline protocols for the induction and analysis of arenavirus-mediated pathogenesis in the NZB/LCMV model, including mouse infection, virus titer determination, platelet counting, phenotypic analysis of virus-specific T cells, and assessment of vascular permeability.

The tyrosine kinase receptor Tyro3 enhances lifespan and neuropeptide Y (Npy) neuron survival in the mouse anorexia (anx) mutation

Severe appetite and weight loss define the eating disorder anorexia nervosa, and can also accompany the progression of some neurodegenerative disorders such as amyotrophic lateral sclerosis (ALS). Although acute loss of hypothalamic neurons that produce appetite-stimulating neuropeptide Y (Npy) and agouti-related peptide (Agrp) in adult mice or in mice homozygous for the anorexia (anx) mutation causes aphagia, our understanding of the factors that help maintain appetite regulatory circuitry is limited. Here we identify a mutation (C19T) that converts an arginine to a tryptophan (R7W) in the TYRO3 protein tyrosine kinase 3 (Tyro3) gene, which resides within the anx critical interval, as contributing to the severity of anx phenotypes. Our observation that, like Tyro3-/- mice, anx/anx mice exhibit abnormal secondary platelet aggregation suggested that the C19T Tyro3 variant might have functional consequences. Tyro3 is expressed in the hypothalamus and other brain regions affected by the anx mutation, and its mRNA localization appeared abnormal in anx/anx brains by postnatal day 19 (P19). The presence of wild-type Tyro3 transgenes, but not an R7W-Tyro3 transgene, doubled the weight and lifespans of anx/anx mice and near-normal numbers of hypothalamic Npy-expressing neurons were present in Tyro3-transgenic anx/anx mice at P19. Although no differences in R7W-Tyro3 signal sequence function or protein localization were discernible in vitro, distribution of R7W-Tyro3 protein differed from that of Tyro3 protein in the cerebellum of transgenic wild-type mice. Thus, R7W-Tyro3 protein localization deficits are only detectable in vivo Further analyses revealed that the C19T Tyro3 mutation is present in a few other mouse strains, and hence is not the causative anx mutation, but rather an anx modifier. Our work shows that Tyro3 has prosurvival roles in the appetite regulatory circuitry and could also provide useful insights towards the development of interventions targeting detrimental weight loss.

Complex Genetics of Cardiovascular Traits in Mice: F2-Mapping of QTLs and Their Underlying Genes

In this chapter, we will use the example of the identification of Tnni3k as a modulator of cardiac conduction to introduce you to the use of a murine F2-generation intercross as a powerful method for the identification of novel genes relevant for cardiovascular traits. Murine F2-progeny is a genetically diverse panel of mice with differences in phenotype manifestations, e.g. cardiovascular traits such as cardiomyopathy and ECG parameters. This chapter discusses the best strategies for using F2-mice for genetic mapping. Moreover, we provide an example of the feasibility of identification of new genes modulating cardiac function utilizing the technique of mapping quantitative trait loci (QTLs) and a systems genetics integration of available genetic, gene expression, and phenotypic data.

Temporal dynamics of the developing lung transcriptome in three common inbred strains of laboratory mice reveals multiple stages of postnatal alveolar development

To characterize temporal patterns of transcriptional activity during normal lung development, we generated genome wide gene expression data for 26 pre- and post-natal time points in three common inbred strains of laboratory mice (C57BL/6J, A/J, and C3H/HeJ). Using Principal Component Analysis and least squares regression modeling, we identified both strain-independent and strain-dependent patterns of gene expression. The 4,683 genes contributing to the strain-independent expression patterns were used to define a murine Developing Lung Characteristic Subtranscriptome (mDLCS). Regression modeling of the Principal Components supported the four canonical stages of mammalian embryonic lung development (embryonic, pseudoglandular, canalicular, saccular) defined previously by morphology and histology. For postnatal alveolar development, the regression model was consistent with four stages of alveolarization characterized by episodic transcriptional activity of genes related to pulmonary vascularization. Genes expressed in a strain-dependent manner were enriched for annotations related to neurogenesis, extracellular matrix organization, and Wnt signaling. Finally, a comparison of mouse and human transcriptomics from pre-natal stages of lung development revealed conservation of pathways associated with cell cycle, axon guidance, immune function, and metabolism as well as organism-specific expression of genes associated with extracellular matrix organization and protein modification. The mouse lung development transcriptome data generated for this study serves as a unique reference set to identify genes and pathways essential for normal mammalian lung development and for investigations into the developmental origins of respiratory disease and cancer. The gene expression data are available from the Gene Expression Omnibus (GEO) archive (GSE74243). Temporal expression patterns of mouse genes can be investigated using a study specific web resource (http://lungdevelopment.jax.org).

Striking Immune Phenotypes in Gene-Targeted Mice Are Driven by a Copy-Number Variant Originating from a Commercially Available C57BL/6 Strain

We describe a homozygous copy-number variant that disrupts the function of Dock2 in a commercially available C57BL/6 mouse strain that is widely used for backcrossing. This Dock2 allele was presumed to have spontaneously arisen in a colony of Irf5 knockout mice. We discovered that this allele has actually been inadvertently backcrossed into multiple mutant mouse lines, including two engineered to be deficient in Siae and Cmah. This particular commercially obtained subline of C57BL/6 mice also exhibits several striking immune phenotypes that have been previously described in the context of Dock2 deficiency. Inadvertent backcrossing of a number of gene-targeted mice into this background has complicated the interpretation of several immunological studies. In light of these findings, published studies involving immune or hematopoietic phenotypes in which these C57BL/6 mice have been used as controls, as experimental animals, or for backcrossing will need to be reinterpreted.

Efficient analysis of mouse genome sequences reveal many nonsense variants

Genetic polymorphisms in coding genes play an important role when using mouse inbred strains as research models. They have been shown to influence research results, explain phenotypical differences between inbred strains, and increase the amount of interesting gene variants present in the many available inbred lines. SPRET/Ei is an inbred strain derived from Mus spretus that has ∼1% sequence difference with the C57BL/6J reference genome. We obtained a listing of all SNPs and insertions/deletions (indels) present in SPRET/Ei from the Mouse Genomes Project (Wellcome Trust Sanger Institute) and processed these data to obtain an overview of all transcripts having nonsynonymous coding sequence variants. We identified 8,883 unique variants affecting 10,096 different transcripts from 6,328 protein-coding genes, which is about 28% of all coding genes. Because only a subset of these variants results in drastic changes in proteins, we focused on variations that are nonsense mutations that ultimately resulted in a gain of a stop codon. These genes were identified by in silico changing the C57BL/6J coding sequences to the SPRET/Ei sequences, converting them to amino acid (AA) sequences, and comparing the AA sequences. All variants and transcripts affected were also stored in a database, which can be browsed using a SPRET/Ei M. spretus variants web tool (www.spretus.org), including a manual. We validated the tool by demonstrating the loss of function of three proteins predicted to be severely truncated, namely Fas, IRAK2, and IFNγR1.

The Mouse Solitary Odorant Receptor Gene Promoters as Models for the Study of Odorant Receptor Gene Choice

BACKGROUND

In vertebrates, several anatomical regions located within the nasal cavity mediate olfaction. Among these, the main olfactory epithelium detects most conventional odorants. Olfactory sensory neurons, provided with cilia exposed to the air, detect volatile chemicals via an extremely large family of seven-transmembrane chemoreceptors named odorant receptors. Their genes are expressed in a monogenic and monoallelic fashion: a single allele of a single odorant receptor gene is transcribed in a given mature neuron, through a still uncharacterized molecular mechanism known as odorant receptor gene choice.

AIM

Odorant receptor genes are typically arranged in genomic clusters, but a few are isolated (we call them solitary) from the others within a region broader than 1 Mb upstream and downstream with respect to their transcript's coordinates. The study of clustered genes is problematic, because of redundancy and ambiguities in their regulatory elements: we propose to use the solitary genes as simplified models to understand odorant receptor gene choice.

PROCEDURES

Here we define number and identity of the solitary genes in the mouse genome (C57BL/6J), and assess the conservation of the solitary status in some mammalian orthologs. Furthermore, we locate their putative promoters, predict their homeodomain binding sites (commonly present in the promoters of odorant receptor genes) and compare candidate promoter sequences with those of wild-caught mice. We also provide expression data from histological sections.

RESULTS

In the mouse genome there are eight intact solitary genes: Olfr19 (M12), Olfr49, Olfr266, Olfr267, Olfr370, Olfr371, Olfr466, Olfr1402; five are conserved as solitary in rat. These genes are all expressed in the main olfactory epithelium of three-day-old mice. The C57BL/6J candidate promoter of Olfr370 has considerably varied compared to its wild-type counterpart. Within the putative promoter for Olfr266 a homeodomain binding site is predicted. As a whole, our findings favor Olfr266 as a model gene to investigate odorant receptor gene choice.

Neurotoxicity to DRG neurons varies between rodent strains treated with cisplatin and bortezomib

Chemotherapy-induced peripheral neuropathy (CIPN) is a major dose limiting side effect that can lead to long-term morbidity. Approximately one-third of patients receiving chemotherapy with taxanes, vinca alkaloids, platinum compounds or proteasome inhibitors develop this toxic side effect. It is not possible to predict who will get CIPN, however, genetic susceptibility may play a role. We explored this hypothesis using an established in vitro dorsal root ganglia neurite outgrowth (DRG-NOG) assay to assess possible genetic influences for cisplatin- and bortezomib-induced neurotoxicity. Almost all previous in vitro studies have used rats or mice. We compared DRG-NOG between four genetically defined, inbred mouse strains (C57BL/6J, DBA/2J, BALB/cJ, and C3H/HeJ) and one rat strain (Sprague Dawley). Our studies found differences in cisplatin and bortezomib-induced neurotoxicity between mouse and rat strains and between the different mouse strains. C57BL/6J and Balb/cJ DRG-NOG was more sensitive to cisplatin than DBA/2J and C3H/HeJ DRG-NOG, and all mouse strains were more sensitive to cisplatin than rat. Bortezomib induced a biphasic dose response in DBA/2J and C3H/H3J mice. C57BL/6J DRG-NOG was most sensitive and Balb/cJ DRG-NOG was least sensitive to bortezomib. Our animal data supports the hypothesis that genetic background may play a role in CIPN and care must be taken when rodent models are used to better understand the contribution of genetics in patient susceptibility to CIPN.

Genetic Architectures of Quantitative Variation in RNA Editing Pathways

RNA editing refers to post-transcriptional processes that alter the base sequence of RNA. Recently, hundreds of new RNA editing targets have been reported. However, the mechanisms that determine the specificity and degree of editing are not well understood. We examined quantitative variation of site-specific editing in a genetically diverse multiparent population, Diversity Outbred mice, and mapped polymorphic loci that alter editing ratios globally for C-to-U editing and at specific sites for A-to-I editing. An allelic series in the C-to-U editing enzyme Apobec1 influences the editing efficiency of Apob and 58 additional C-to-U editing targets. We identified 49 A-to-I editing sites with polymorphisms in the edited transcript that alter editing efficiency. In contrast to the shared genetic control of C-to-U editing, most of the variable A-to-I editing sites were determined by local nucleotide polymorphisms in proximity to the editing site in the RNA secondary structure. Our results indicate that RNA editing is a quantitative trait subject to genetic variation and that evolutionary constraints have given rise to distinct genetic architectures in the two canonical types of RNA editing.

Using the GEMM-ESC strategy to study gene function in mouse models

Preclinical in vivo validation of target genes for therapeutic intervention requires careful selection and characterization of the most suitable animal model in order to assess the role of these genes in a particular process or disease. To this end, genetically engineered mouse models (GEMMs) are typically used. However, the appropriate engineering of these models is often cumbersome and time consuming. Recently, we and others described a modular approach for fast-track modification of existing GEMMs by re-derivation of embryonic stem cells (ESCs) that can be modified by recombinase-mediated transgene insertion and subsequently used for the production of chimeric mice. This 'GEMM-ESC strategy' allows for rapid in vivo analysis of gene function in the chimeras and their offspring. Moreover, this strategy is compatible with CRISPR/Cas9-mediated genome editing. This protocol describes when and how to use the GEMM-ESC strategy effectively, and it provides a detailed procedure for re-deriving and manipulating GEMM-ESCs under feeder- and serum-free conditions. This strategy produces transgenic mice with the desired complex genotype faster than traditional methods: generation of validated GEMM-ESC clones for controlled transgene integration takes 9-12 months, and recombinase-mediated transgene integration and chimeric cohort production takes 2-3 months. The protocol requires skills in embryology, stem cell biology and molecular biology, and it is ideally performed within, or in close collaboration with, a transgenic facility.

Creating reference gene annotation for the mouse C57BL6/J genome assembly

Annotation on the reference genome of the C57BL6/J mouse has been an ongoing project ever since the draft genome was first published. Initially, the principle focus was on the identification of all protein-coding genes, although today the importance of describing long non-coding RNAs, small RNAs, and pseudogenes is recognized. Here, we describe the progress of the GENCODE mouse annotation project, which combines manual annotation from the HAVANA group with Ensembl computational annotation, alongside experimental and in silico validation pipelines from other members of the consortium. We discuss the more recent incorporation of next-generation sequencing datasets into this workflow, including the usage of mass-spectrometry data to potentially identify novel protein-coding genes. Finally, we will outline how the C57BL6/J genebuild can be used to gain insights into the variant sites that distinguish different mouse strains and species.

Comprehensive comparative homeobox gene annotation in human and mouse

Homeobox genes are a group of genes coding for transcription factors with a DNA-binding helix-turn-helix structure called a homeodomain and which play a crucial role in pattern formation during embryogenesis. Many homeobox genes are located in clusters and some of these, most notably the HOX genes, are known to have antisense or opposite strand long non-coding RNA (lncRNA) genes that play a regulatory role. Because automated annotation of both gene clusters and non-coding genes is fraught with difficulty (over-prediction, under-prediction, inaccurate transcript structures), we set out to manually annotate all homeobox genes in the mouse and human genomes. This includes all supported splice variants, pseudogenes and both antisense and flanking lncRNAs. One of the areas where manual annotation has a significant advantage is the annotation of duplicated gene clusters. After comprehensive annotation of all homeobox genes and their antisense genes in human and in mouse, we found some discrepancies with the current gene set in RefSeq regarding exact gene structures and coding versus pseudogene locus biotype. We also identified previously un-annotated pseudogenes in the DUX, Rhox and Obox gene clusters, which helped us re-evaluate and update the gene nomenclature in these regions. We found that human homeobox genes are enriched in antisense lncRNA loci, some of which are known to play a role in gene or gene cluster regulation, compared to their mouse orthologues. Of the annotated set of 241 human protein-coding homeobox genes, 98 have an antisense locus (41%) while of the 277 orthologous mouse genes, only 62 protein coding gene have an antisense locus (22%), based on publicly available transcriptional evidence.

Collaborative Cross and Diversity Outbred data resources in the Mouse Phenome Database

The Mouse Phenome Database was originally conceived as a platform for the integration of phenotype data collected on a defined collection of 40 inbred mouse strains--the "phenome panel." This model provided an impetus for community data sharing, and integration was readily achieved through the reproducible genotypes of the phenome panel strains. Advances in the development of mouse populations lead to an expanded role of the Mouse Phenome Database to encompass new strain panels and inbred strain crosses. The recent introduction of the Collaborative Cross and Diversity Outbred mice, which share an extensive pool of genetic variation from eight founder inbred strains, presents new opportunities and challenges for community data resources. A wide variety of molecular and clinical phenotypes are being collected across genotypes, tissues, ages, environmental exposures, interventions, and treatments. The Mouse Phenome Database provides a framework for retrieval, integration, analysis, and display of these data, enabling them to be evaluated in the context of existing data from standard inbred strains. Primary data in the Mouse Phenome Database are supported by extensive metadata on protocols and procedures. These are centrally curated to ensure accuracy and reproducibility and to provide data in consistent formats. The Mouse Phenome Database represents an established and growing community data resource for mouse phenotype data and encourages submissions from new mouse resources, enabling investigators to integrate existing data into their studies of the phenotypic consequences of genetic variation.

Mouse Genome Database: From sequence to phenotypes and disease models

The Mouse Genome Database (MGD, www.informatics.jax.org) is the international scientific database for genetic, genomic, and biological data on the laboratory mouse to support the research requirements of the biomedical community. To accomplish this goal, MGD provides broad data coverage, serves as the authoritative standard for mouse nomenclature for genes, mutants, and strains, and curates and integrates many types of data from literature and electronic sources. Among the key data sets MGD supports are: the complete catalog of mouse genes and genome features, comparative homology data for mouse and vertebrate genes, the authoritative set of Gene Ontology (GO) annotations for mouse gene functions, a comprehensive catalog of mouse mutations and their phenotypes, and a curated compendium of mouse models of human diseases. Here, we describe the data acquisition process, specifics about MGD's key data areas, methods to access and query MGD data, and outreach and user help facilities. genesis 53:458–473, 2015. © 2015 The Authors. Genesis Published by Wiley Periodicals, Inc.

The pervasiveness and plasticity of circadian oscillations: the coupled circadian-oscillators framework

MOTIVATION

Circadian oscillations have been observed in animals, plants, fungi and cyanobacteria and play a fundamental role in coordinating the homeostasis and behavior of biological systems. Genetically encoded molecular clocks found in nearly every cell, based on negative transcription/translation feedback loops and involving only a dozen genes, play a central role in maintaining these oscillations. However, high-throughput gene expression experiments reveal that in a typical tissue, a much larger fraction ([Formula: see text]) of all transcripts oscillate with the day-night cycle and the oscillating species vary with tissue type suggesting that perhaps a much larger fraction of all transcripts, and perhaps also other molecular species, may bear the potential for circadian oscillations.

RESULTS

To better quantify the pervasiveness and plasticity of circadian oscillations, we conduct the first large-scale analysis aggregating the results of 18 circadian transcriptomic studies and 10 circadian metabolomic studies conducted in mice using different tissues and under different conditions. We find that over half of protein coding genes in the cell can produce transcripts that are circadian in at least one set of conditions and similarly for measured metabolites. Genetic or environmental perturbations can disrupt existing oscillations by changing their amplitudes and phases, suppressing them or giving rise to novel circadian oscillations. The oscillating species and their oscillations provide a characteristic signature of the physiological state of the corresponding cell/tissue. Molecular networks comprise many oscillator loops that have been sculpted by evolution over two trillion day-night cycles to have intrinsic circadian frequency. These oscillating loops are coupled by shared nodes in a large network of coupled circadian oscillators where the clock genes form a major hub. Cells can program and re-program their circadian repertoire through epigenetic and other mechanisms.

AVAILABILITY AND IMPLEMENTATION

High-resolution and tissue/condition specific circadian data and networks available at http://circadiomics.igb.uci.edu.

CONTACT

pfbaldi@ics.uci.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Accessing Data Resources in the Mouse Phenome Database for Genetic Analysis of Murine Life Span and Health Span

Understanding the source of genetic variation in aging and using this variation to define the molecular mechanisms of healthy aging require deep and broad quantification of a host of physiological, morphological, and behavioral endpoints. The murine model is a powerful system in which to understand the relations across age-related phenotypes and to identify research models with variation in life span and health span. The Jackson Laboratory Nathan Shock Center of Excellence in the Basic Biology of Aging has performed broad characterization of aging in genetically diverse laboratory mice and has placed these data, along with data from several other major aging initiatives, into the interactive Mouse Phenome Database. The data may be accessed and analyzed by researchers interested in finding mouse models for specific aging processes, age-related health and disease states, and for genetic analysis of aging variation and trait covariation. We expect that by placing these data in the hands of the aging community that there will be (a) accelerated genetic analyses of aging processes, (b) discovery of genetic loci regulating life span, (c) identification of compelling correlations between life span and susceptibility for age-related disorders, and (d) discovery of concordant genomic loci influencing life span and aging phenotypes between mouse and humans.

A comparative analysis of C57BL/6J and 6N substrains; chemokine/cytokine expression and susceptibility to laser-induced choroidal neovascularization

Age-related macular degeneration (AMD) is the most prevalent cause of blindness in the elderly. To study potential underlying mechanisms of AMD, animal models are utilized, focusing mostly on mice. Recently, genomic and phenotypic differences between the so-called control substrains, C57BL/6J and C57BL/6N, have been described in models of ocular and non-ocular diseases. In particular, the rd8 mutation of the Crb1 gene present in the C57BL/6N has been shown to impact certain ocular phenotypes and appears to augment phenotypes generally associated with inflammation. Here, we investigated angiogenic factor and cytokine expression using pathway arrays as well as the susceptibility to laser-induced choroidal neovascularization (CNV), a model of wet AMD, in the two substrains. Age-matched 3-month-old C57BL/6J and C57BL/6N animals differed in gene expression levels for angiogenic factors and cytokines, with 6N animals expressing higher levels of inflammatory markers than 6Js. Yet laser-induced CNV was comparable in size between the two substrains. This lack of difference in CNV size was correlated with a gene expression profile that was comparable between the two substrains, due to the fact that the degree of change in gene expression of inflammatory markers after CNV was blunted in 6N mice. In summary, significant gene expression differences exist between C57BL/6J and C57BL/6N animals, reinforcing the notion that appropriate litter-mate controls or genetic background controls need to be used. Contrary to our expectation, CNV was not augmented in 6N animals, suggesting that low chronic inflammation in the RPE might provide a level of pre-conditioning and protection against stress.

Parallel universes of Black Six biology

Creation of lethal and synthetic lethal mutations in an experimental organism is a cornerstone of genetic dissection of gene function, and is related to the concept of an essential gene. Common inbred mouse strains carry background mutations, which can act as genetic modifiers, interfering with the assignment of gene essentiality. The inbred strain C57BL/6J, commonly known as "Black Six", stands out, as it carries a spontaneous homozygous deletion in the nicotinamide nucleotide transhydrogenase (Nnt) gene [GenBank: AH009385.2], resulting in impairment of steroidogenic mitochondria of the adrenal gland, and a multitude of indirect modifier effects, coming from alteration of glucocorticoid-regulated processes. Over time, the popular strain has been used, by means of gene targeting technology, to assign "essential" and "redundant" qualifiers to numerous genes, thus creating an internally consistent "parallel universe" of knowledge. It is unrealistic to suggest phasing-out of this strain, given the scope of shared resources built around it, however, continuing on the road of "strain-unawareness" will result in profound waste of effort, particularly where translational research is concerned. The review analyzes the historical roots of this phenomenon and proposes that building of "parallel universes" should be urgently made visible to a critical reader by obligatory use of unambiguous and persistent tags in publications and databases, such as hypertext links, pointing to a vendor's strain description web page, or to a digital object identifier (d.o.i.) of the original publication, so that any research done exclusively in C57BL/6J, could be easily identified.

REVIEWERS

This article was reviewed by Dr. Neil Smalheiser and Dr. Miguel Andrade-Navarro.

Genome-wide patterns of differentiation among house mouse subspecies

One approach to understanding the genetic basis of speciation is to scan the genomes of recently diverged taxa to identify highly differentiated regions. The house mouse, Mus musculus, provides a useful system for the study of speciation. Three subspecies (M. m. castaneus, M. m. domesticus, and M. m. musculus) diverged ∼350 KYA, are distributed parapatrically, show varying degrees of reproductive isolation in laboratory crosses, and hybridize in nature. We sequenced the testes transcriptomes of multiple wild-derived inbred lines from each subspecies to identify highly differentiated regions of the genome, to identify genes showing high expression divergence, and to compare patterns of differentiation among subspecies that have different demographic histories and exhibit different levels of reproductive isolation. Using a sliding-window approach, we found many genomic regions with high levels of sequence differentiation in each of the pairwise comparisons among subspecies. In all comparisons, the X chromosome was more highly differentiated than the autosomes. Sequence differentiation and expression divergence were greater in the M. m. domesticus-M. m. musculus comparison than in either pairwise comparison with M. m. castaneus, which is consistent with laboratory crosses that show the greatest reproductive isolation between M. m. domesticus and M. m. musculus. Coalescent simulations suggest that differences in estimates of effective population size can account for many of the observed patterns. However, there was an excess of highly differentiated regions relative to simulated distributions under a wide range of demographic scenarios. Overlap of some highly differentiated regions with previous results from QTL mapping and hybrid zone studies points to promising candidate regions for reproductive isolation.

Identification of structural variation in mouse genomes

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.

Implications of genomic signatures in the differential vulnerability to fetal alcohol exposure in C57BL/6 and DBA/2 mice

Maternal alcohol consumption inflicts a multitude of phenotypic consequences that range from undetectable changes to severe dysmorphology. Using tightly controlled murine studies that deliver precise amounts of alcohol at discrete developmental stages, our group and other labs demonstrated in prior studies that the C57BL/6 and DBA/2 inbred mouse strains display differential susceptibility to the teratogenic effects of alcohol. Since the phenotypic diversity extends beyond the amount, dosage and timing of alcohol exposure, it is likely that an individual's genetic background contributes to the phenotypic spectrum. To identify the genomic signatures associated with these observed differences in alcohol-induced dysmorphology, we conducted a microarray-based transcriptome study that also interrogated the genomic signatures between these two lines based on genetic background and alcohol exposure. This approach is called a gene x environment (GxE) analysis; one example of a GxE interaction would be a gene whose expression level increases in C57BL/6, but decreases in DBA/2 embryos, following alcohol exposure. We identified 35 candidate genes exhibiting GxE interactions. To identify cis-acting factors that mediated these interactions, we interrogated the proximal promoters of these 35 candidates and found 241 single nucleotide variants (SNVs) in 16 promoters. Further investigation indicated that 186 SNVs (15 promoters) are predicted to alter transcription factor binding. In addition, 62 SNVs created, removed or altered the placement of a CpG dinucleotide in 13 of the proximal promoters, 53 of which overlapped putative transcription factor binding sites. These 53 SNVs are also our top candidates for future studies aimed at examining the effects of alcohol on epigenetic gene regulation.

Going forward with genetics: recent technological advances and forward genetics in mice

Forward genetic analysis is an unbiased approach for identifying genes essential to defined biological phenomena. When applied to mice, it is one of the most powerful methods to facilitate understanding of the genetic basis of human biology and disease. The speed at which disease-causing mutations can be identified in mutagenized mice has been markedly increased by recent advances in DNA sequencing technology. Creating and analyzing mutant phenotypes may therefore become rate-limiting in forward genetic experimentation. We review the forward genetic approach and its future in the context of recent technological advances, in particular massively parallel DNA sequencing, induced pluripotent stem cells, and haploid embryonic stem cells.

Maternal bias and escape from X chromosome imprinting in the midgestation mouse placenta

To investigate the epigenetic landscape at the interface between mother and fetus, we provide a comprehensive analysis of parent-of-origin bias in the mouse placenta. Using F1 interspecies hybrids between mus musculus (C57BL/6J) and mus musculus castaneus, we sequenced RNA from 23 individual midgestation placentas, five late stage placentas, and two yolk sac samples and then used SNPs to determine whether transcripts were preferentially generated from the maternal or paternal allele. In the placenta, we find 103 genes that show significant and reproducible parent-of-origin bias, of which 78 are novel candidates. Most (96%) show a strong maternal bias which we demonstrate, via multiple mathematical models, pyrosequencing, and FISH, is not due to maternal decidual contamination. Analysis of the X chromosome also reveals paternal expression of Xist and several genes that escape inactivation, most significantly Alas2, Fhl1, and Slc38a5. Finally, sequencing individual placentas allowed us to reveal notable expression similarity between littermates. In all, we observe a striking preference for maternal transcription in the midgestation mouse placenta and a dynamic imprinting landscape in extraembryonic tissues, reflecting the complex nature of epigenetic pathways in the placenta.

Improving mRNA 5' coding sequence determination in the mouse genome

The incomplete determination of the mRNA 5' end sequence may lead to the incorrect assignment of the first AUG codon and to errors in the prediction of the encoded protein product. Due to the significance of the mouse as a model organism in biomedical research, we performed a systematic identification of coding regions at the 5' end of all known mouse mRNAs, using an automated expressed sequence tag (EST)-based approach which we have previously described. By parsing almost 4 million BLAT alignments we found 351 mouse loci, out of 20,221 analyzed, in which an extension of the mRNA 5' coding region was identified. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for Apc2 and Mknk2 cDNAs. We also generated a list of 16,330 mouse mRNAs where the presence of an in-frame stop codon upstream of the known start codon indicates completeness of the coding sequence at 5' end in the current form. Systematic searches in the main mouse genome databases and genome browsers showed that 82% of our results are original and have not been identified by their annotation pipelines. Moreover, the same information is not easily derivable from RNA-Seq data, due to short sequence length and laboriousness in building full-length transcript structures. In conclusion, our results improve the determination of full-length 5' coding sequences and might be useful in order to reduce errors when studying mouse gene structure and function in biomedical research.

Mouse phenome database

The Mouse Phenome Database (MPD; phenome.jax.org) was launched in 2001 as the data coordination center for the international Mouse Phenome Project. MPD integrates quantitative phenotype, gene expression and genotype data into a common annotated framework to facilitate query and analysis. MPD contains >3500 phenotype measurements or traits relevant to human health, including cancer, aging, cardiovascular disorders, obesity, infectious disease susceptibility, blood disorders, neurosensory disorders, drug addiction and toxicity. Since our 2012 NAR report, we have added >70 new data sets, including data from Collaborative Cross lines and Diversity Outbred mice. During this time we have completely revamped our homepage, improved search and navigational aspects of the MPD application, developed several web-enabled data analysis and visualization tools, annotated phenotype data to public ontologies, developed an ontology browser and released new single nucleotide polymorphism query functionality with much higher density coverage than before. Here, we summarize recent data acquisitions and describe our latest improvements.

Comparative analysis of context-dependent mutagenesis using human and mouse models

Substitution rates strongly depend on their nucleotide context. One of the most studied examples is the excess of C > T mutations in the CG context in various groups of organisms, including vertebrates. Studies on the molecular mechanisms underlying this mutation regularity have provided insights into evolution, mutagenesis, and cancer development. Recently several other hypermutable motifs were identified in the human genome. There is an increased frequency of T > C mutations in the second position of the words ATTG and ATAG and an increased frequency of A > C mutations in the first position of the word ACAA. For a better understanding of evolution, it is of interest whether these mutation regularities are human specific or present in other vertebrates, as their presence might affect the validity of currently used substitution models and molecular clocks. A comprehensive analysis of mutagenesis in 4 bp mutation contexts requires a vast amount of mutation data. Such data may be derived from the comparisons of individual genomes or from single nucleotide polymorphism (SNP) databases. Using this approach, we performed a systematical comparison of mutation regularities within 2-4 bp contexts in Mus musculus and Homo sapiens and uncovered that even closely related organisms may have notable differences in context-dependent mutation regularities.

QTL analysis of dietary obesity in C57BL/6byj X 129P3/J F2 mice: diet- and sex-dependent effects

Obesity is a heritable trait caused by complex interactions between genes and environment, including diet. Gene-by-diet interactions are difficult to study in humans because the human diet is hard to control. Here, we used mice to study dietary obesity genes, by four methods. First, we bred 213 F₂ mice from strains that are susceptible [C57BL/6ByJ (B6)] or resistant [129P3/J (129)] to dietary obesity. Percent body fat was assessed after mice ate low-energy diet and again after the same mice ate high-energy diet for 8 weeks. Linkage analyses identified QTLs associated with dietary obesity. Three methods were used to filter candidate genes within the QTL regions: (a) association mapping was conducted using >40 strains; (b) differential gene expression and (c) comparison of genomic DNA sequence, using two strains closely related to the progenitor strains from Experiment 1. The QTL effects depended on whether the mice were male or female or which diet they were recently fed. After feeding a low-energy diet, percent body fat was linked to chr 7 (LOD = 3.42). After feeding a high-energy diet, percent body fat was linked to chr 9 (Obq5; LOD = 3.88), chr 12 (Obq34; LOD = 3.88), and chr 17 (LOD = 4.56). The Chr 7 and 12 QTLs were sex dependent and all QTL were diet-dependent. The combination of filtering methods highlighted seven candidate genes within the QTL locus boundaries: Crx, Dmpk, Ahr, Mrpl28, Glo1, Tubb5, and Mut. However, these filtering methods have limitations so gene identification will require alternative strategies, such as the construction of congenics with very small donor regions.

Genomic structure and variation of nuclear factor (erythroid-derived 2)-like 2

High-density mapping of mammalian genomes has enabled a wide range of genetic investigations including the mapping of polygenic traits, determination of quantitative trait loci, and phylogenetic comparison. Genome sequencing analysis of inbred mouse strains has identified high-density single nucleotide polymorphisms (SNPs) for investigation of complex traits, which has become a useful tool for biomedical research of human disease to alleviate ethical and practical problems of experimentation in humans. Nuclear factor (erythroid-derived 2)-like 2 (NRF2) encodes a key host defense transcription factor. This review describes genetic characteristics of human NRF2 and its homologs in other vertebrate species. NRF2 is evolutionally conserved and shares sequence homology among species. Compilation of publically available SNPs and other genetic mutations shows that human NRF2 is highly polymorphic with a mutagenic frequency of 1 per every 72 bp. Functional at-risk alleles and haplotypes have been demonstrated in various human disorders. In addition, other pathogenic alterations including somatic mutations and misregulated epigenetic processes in NRF2 have led to oncogenic cell survival. Comprehensive information from the current review addresses association of NRF2 variation and disease phenotypes and supports the new insights into therapeutic strategies.

Deconstructing Mus gemischus: advances in understanding ancestry, structure, and variation in the genome of the laboratory mouse

The laboratory mouse is an artificial construct with a complex relationship to its natural ancestors. In 2002, the mouse became the first mammalian model organism with a reference genome. Importantly, the mouse genome sequence was assembled from data on a single inbred laboratory strain, C57BL/6. Several large-scale genetic variant discovery efforts have been conducted, resulting in a catalog of tens of millions of SNPs and structural variants. High-density genotyping arrays covering a subset of those variants have been used to produce hundreds of millions of genotypes in laboratory stocks and a small number of wild mice. These landmark resources now enable us to determine relationships among laboratory mice, assign local ancestry at fine scale, resolve important controversies, and identify a new set of challenges-most importantly, the troubling scarcity of genetic data on the very natural populations from which the laboratory mouse was derived. Our aim with this review is to provide the reader with an historical context for the mouse as a model organism and to explain how practical decisions made in the past have influenced both the architecture of the laboratory mouse genome and the design and execution of current large-scale resources. We also provide examples on how the accomplishments of the past decade can be used by researchers to streamline the use of mice in their experiments and correctly interpret results. Finally, we propose future steps that will enable the mouse community to extend its successes in the decade to come.

Association studies in outbred mice in a new era of full-genome sequencing

Thousands of loci that contribute to quantitative traits in outbred crosses of mice have been reported over the last two decades. In this review we discuss how outbred mouse populations can be used to map and identify the genes and sequence variants that give rise to quantitative variation. We discuss heterogeneous stocks, the diversity outbred, and commercially available outbred populations of mice. All of these populations are descended from a small number of progenitor strains. The availability of the complete sequence of laboratory strains means that in many cases it will be possible to reconstruct the genomes of the outbred animals so that in a genetic association study we can detect the effect of all variants, a situation that has so far eluded studies in completely outbred populations. These resources constitute a major advance and make it possible to progress from a quantitative trait locus to a gene at an unprecedented speed.