The Genome of C57BL/6J "Eve", the Mother of the Laboratory Mouse Genome Reference Strain.

Isogenic laboratory mouse strains enhance reproducibility because individual animals are genetically identical. For the most widely used isogenic strain, C57BL/6, there exists a wealth of genetic, phenotypic, and genomic data, including a high-quality reference genome (GRCm38.p6). Now 20 years after the first release of the mouse reference genome, C57BL/6J mice are at least 26 inbreeding generations removed from GRCm38 and the strain is now maintained with periodic reintroduction of cryorecovered mice derived from a single breeder pair, aptly named Adam and Eve. To provide an update to the mouse reference genome that more accurately represents the genome of today's C57BL/6J mice, we took advantage of long read, short read, and optical mapping technologies to generate a de novo assembly of the C57BL/6J Eve genome (B6Eve). Using these data, we have addressed recurring variants observed in previous mouse genomic studies. We have also identified structural variations, closed gaps in the mouse reference assembly, and revealed previously unannotated coding sequences. This B6Eve assembly explains discrepant observations that have been associated with GRCm38-based analyses, and will inform a reference genome that is more representative of the C57BL/6J mice that are in use today.

Mutations in CTNNA1 cause butterfly-shaped pigment dystrophy and perturbed retinal pigment epithelium integrity.

Butterfly-shaped pigment dystrophy is an eye disease characterized by lesions in the macula that can resemble the wings of a butterfly. Here we report the identification of heterozygous missense mutations in the CTNNA1 gene (encoding α-catenin 1) in three families with butterfly-shaped pigment dystrophy. In addition, we identified a Ctnna1 missense mutation in a chemically induced mouse mutant, tvrm5. Parallel clinical phenotypes were observed in the retinal pigment epithelium (RPE) of individuals with butterfly-shaped pigment dystrophy and in tvrm5 mice, including pigmentary abnormalities, focal thickening and elevated lesions, and decreased light-activated responses. Morphological studies in tvrm5 mice demonstrated increased cell shedding and the presence of large multinucleated RPE cells, suggesting defects in intercellular adhesion and cytokinesis. This study identifies CTNNA1 gene variants as a cause of macular dystrophy, indicates that CTNNA1 is involved in maintaining RPE integrity and suggests that other components that participate in intercellular adhesion may be implicated in macular disease.

Genetic Analysis of Substrain Divergence in Non-Obese Diabetic (NOD) Mice.

The non-obese diabetic (NOD) mouse is a polygenic model for type 1 diabetes that is characterized by insulitis, a leukocytic infiltration of the pancreatic islets. During ~35 years since the original inbred strain was developed in Japan, NOD substrains have been established at different laboratories around the world. Although environmental differences among NOD colonies capable of impacting diabetes incidence have been recognized, differences arising from genetic divergence have not been analyzed previously. We use both mouse diversity array and whole-exome capture sequencing platforms to identify genetic differences distinguishing five NOD substrains. We describe 64 single-nucleotide polymorphisms, and two short indels that differ in coding regions of the five NOD substrains. A 100-kb deletion on Chromosome 3 distinguishes NOD/ShiLtJ and NOD/ShiLtDvs from three other substrains, whereas a 111-kb deletion in the Icam2 gene on Chromosome 11 is unique to the NOD/ShiLtDvs genome. The extent of genetic divergence for NOD substrains is compared with similar studies for C57BL6 and BALB/c substrains. As mutations are fixed to homozygosity by continued inbreeding, significant differences in substrain phenotypes are to be expected. These results emphasize the importance of using embryo freezing methods to minimize genetic drift within substrains and of applying appropriate genetic nomenclature to permit substrain recognition when one is used.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2.

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

Development and validation of the JAX Cancer Treatment Profile™ for detection of clinically actionable mutations in solid tumors.

BACKGROUND: The continued development of targeted therapeutics for cancer treatment has required the concomitant development of more expansive methods for the molecular profiling of the patient's tumor. We describe the validation of the JAX Cancer Treatment Profile™ (JAX-CTP™), a next generation sequencing (NGS)-based molecular diagnostic assay that detects actionable mutations in solid tumors to inform the selection of targeted therapeutics for cancer treatment. METHODS: NGS libraries are generated from DNA extracted from formalin fixed paraffin embedded tumors. Using hybrid capture, the genes of interest are enriched and sequenced on the Illumina HiSeq 2500 or MiSeq sequencers followed by variant detection and functional and clinical annotation for the generation of a clinical report. RESULTS: The JAX-CTP™ detects actionable variants, in the form of single nucleotide variations and small insertions and deletions (≤50 bp) in 190 genes in specimens with a neoplastic cell content of ≥10%. The JAX-CTP™ is also validated for the detection of clinically actionable gene amplifications. CONCLUSIONS: There is a lack of consensus in the molecular diagnostics field on the best method for the validation of NGS-based assays in oncology, thus the importance of communicating methods, as contained in this report. The growing number of targeted therapeutics and the complexity of the tumor genome necessitate continued development and refinement of advanced assays for tumor profiling to enable precision cancer treatment.

Discovery of transgene insertion sites by high throughput sequencing of mate pair libraries.

BACKGROUND: Transgenesis by random integration of a transgene into the genome of a zygote has become a reliable and powerful method for the creation of new mouse strains that express exogenous genes, including human disease genes, tissue specific reporter genes or genes that allow for tissue specific recombination. Nearly 6,500 transgenic alleles have been created by random integration in embryos over the last 30 years, but for the vast majority of these strains, the transgene insertion sites remain uncharacterized. RESULTS: To obtain a complete understanding of how insertion sites might contribute to phenotypic outcomes, to more cost effectively manage transgenic strains, and to fully understand mechanisms of instability in transgene expression, we've developed methodology and a scoring scheme for transgene insertion site discovery using high throughput sequencing data. CONCLUSIONS: Similar to other molecular approaches to transgene insertion site discovery, high-throughput sequencing of standard paired-end libraries is hindered by low signal to noise ratios. This problem is exacerbated when the transgene consists of sequences that are also present in the host genome. We've found that high throughput sequencing data from mate-pair libraries are more informative when compared to data from standard paired end libraries. We also show examples of the genomic regions that harbor transgenes, which have in common a preponderance of repetitive sequences.

GPR179 is required for depolarizing bipolar cell function and is mutated in autosomal-recessive complete congenital stationary night blindness.

Complete congenital stationary night blindness (cCSNB) is a clinically and genetically heterogeneous group of retinal disorders characterized by nonprogressive impairment of night vision, absence of the electroretinogram (ERG) b-wave, and variable degrees of involvement of other visual functions. We report here that mutations in GPR179, encoding an orphan G protein receptor, underlie a form of autosomal-recessive cCSNB. The Gpr179(nob5/nob5) mouse model was initially discovered by the absence of the ERG b-wave, a component that reflects depolarizing bipolar cell (DBC) function. We performed genetic mapping, followed by next-generation sequencing of the critical region and detected a large transposon-like DNA insertion in Gpr179. The involvement of GPR179 in DBC function was confirmed in zebrafish and humans. Functional knockdown of gpr179 in zebrafish led to a marked reduction in the amplitude of the ERG b-wave. Candidate gene analysis of GPR179 in DNA extracted from patients with cCSNB identified GPR179-inactivating mutations in two patients. We developed an antibody against mouse GPR179, which robustly labeled DBC dendritic terminals in wild-type mice. This labeling colocalized with the expression of GRM6 and was absent in Gpr179(nob5/nob5) mutant mice. Our results demonstrate that GPR179 plays a critical role in DBC signal transduction and expands our understanding of the mechanisms that mediate normal rod vision.

Mutation discovery in mice by whole exome sequencing.

We report the development and optimization of reagents for in-solution, hybridization-based capture of the mouse exome. By validating this approach in a multiple inbred strains and in novel mutant strains, we show that whole exome sequencing is a robust approach for discovery of putative mutations, irrespective of strain background. We found strong candidate mutations for the majority of mutant exomes sequenced, including new models of orofacial clefting, urogenital dysmorphology, kyphosis and autoimmune hepatitis.

A mouse TRAPP-related protein is involved in pigmentation.

We identified a new spontaneous recessive mutation in the mouse, mhyp (mosaic hypopigmentation), in a screen for novel proviral integration sites in a multiple ecotropic provirus mapping stock. Integration of an 8.4-kb retrovirus results in mosaic loss of coat pigment in mhyp homozygotes. Patchy loss of pigmentation in the retinal pigmented epithelial layer of the eye with abnormal melanosomes is also evident. We mapped mhyp to mouse chromosome 7 and cloned the underlying gene. mhyp is a defect in the Trappc6a gene. Expression of Trappc6a is markedly diminished in mhyp homozygotes. The normal protein, TRAPPC6A, is a subunit of the TRAPP (transport protein particle) I and II complexes. While TRAPP complexes are essential for ER-to-Golgi and intra-Golgi vesicle trafficking in yeast, TRAPP subunits participate in additional, including post-Golgi, transport events in mammals. The data implicate mammalian TRAPPC6A in vesicle trafficking during melanosome biogenesis.

The development of a highly informative mouse Simple Sequence Length Polymorphism (SSLP) marker set and construction of a mouse family tree using parsimony analysis.

To identify highly informative markers for a large number of commonly employed murine crosses, we selected a subset of the extant mouse simple sequence length polymorphism (SSLP) marker set for further development. Primer pairs for 314 SSLP markers were designed and typed against 54 inbred mouse strains. We designed new PCR primer sequences for the markers selected for multiplexing using the fluorescent dyes FAM, VIC, NED, and ROX. The number of informative markers for C57BL/6J x DBA/2J is 217, with an average spacing of 6.8 centiMorgans (cM). For all other pairs of strains, the mean number of informative markers per cross is 197.0 (SD 37.8) with a mean distance between markers of 6.8 cM (SD 1.1). To confirm map positions of the 224 markers in our set that are polymorphic between Mus musculus and Mus spretus, we used The Jackson Laboratory (TJL) interspecific backcross mapping panel (TJL BSS); 168 (75%) of these markers had not been previously mapped in this cross by other investigators, adding new information to this community map resource. With this large data set, we sought to reconstruct a phylogenetic history of the laboratory mouse using Wagner parsimony analysis. Our results are largely congruent with the known history of inbred mouse strains.

The comprehensive mouse radiation hybrid map densely cross-referenced to the recombination map: a tool to support the sequence assemblies.

We have developed a unique comprehensive mouse radiation hybrid (RH) map of nearly 23,000 markers integrating data from three international genome centers and over 400 independent laboratories. We have cross-referenced this map to the 0.5-cM resolution recombination-based Jackson Laboratory (TJL) backcross panel map, building a complete set of RH framework chromosome maps based on a high density of known-ordered anchor markers. We have systematically typed markers to improve coverage and resolve discrepancies, and have reanalyzed data sets as needed. The cross-linking of the RH and recombination maps has resulted in a highly accurate genome-wide map with consistent marker order. We have compared these linked framework maps to the Ensemble mouse genome sequence assembly, and show that they are a useful medium resolution tool for both validating sequence assembly and elucidating chromosome biology.

Crossover interference in the mouse.

We present an analysis of crossover interference in the mouse genome, on the basis of high-density genotype data from two reciprocal interspecific backcrosses, comprising 188 meioses. Overwhelming evidence was found for strong positive crossover interference with average strength greater than that implied by the Carter-Falconer map function. There was some evidence for interchromosomal variation in the level of interference, with smaller chromosomes exhibiting stronger interference. We further compared the observed numbers of crossovers to previous cytological observations on the numbers of chiasmata and evaluated evidence for the obligate chiasma hypothesis.