We report genome sequences of 17 inbred strains of laboratory mice and identify almost ten times more variants than previously known. We use these genomes to explore the phylogenetic history of the laboratory mouse and to examine the functional consequences of allele-specific variation on transcript abundance, revealing that at least 12% of transcripts show a significant tissue-specific expression bias. By identifying candidate functional variants at 718 quantitative trait loci we show that the molecular nature of functional variants and their position relative to genes vary according to the effect size of the locus. These sequences provide a starting point for a new era in the functional analysis of a key model organism.

Background: The mouse inbred line C57BL/6J is widely used in mouse genetics and its genome has been incorporated into many genetic reference populations. More recently large initiatives such as the International Knockout Mouse Consortium (IKMC) are using the C57BL/6N mouse strain to generate null alleles for all mouse genes. Hence both strains are now widely used in mouse genetics studies. Here we perform a comprehensive genomic and phenotypic analysis of the two strains to identify differences that may influence their underlying genetic mechanisms. Results: We undertake genome sequence comparisons of C57BL/6J and C57BL/6N to identify SNPs, indels and structural variants, with a focus on identifying all coding variants. We annotate 34 SNPs and 2 indels that distinguish C57BL/6J and C57BL/6N coding sequences, as well as 15 structural variants that overlap a gene. In parallel we assess the comparative phenotypes of the two inbred lines utilizing the EMPReSSslim phenotyping pipeline, a broad based assessment encompassing diverse biological systems. We perform additional secondary phenotyping assessments to explore other phenotype domains and to elaborate phenotype differences identified in the primary assessment. We uncover significant phenotypic differences between the two lines, replicated across multiple centers, in a number of physiological, biochemical and behavioral systems. Conclusions: Comparison of C57BL/6J and C57BL/6N demonstrates a range of phenotypic differences that have the potential to impact upon penetrance and expressivity of mutational effects in these strains. Moreover, the sequence variants we identify provide a set of candidate genes for the phenotypic differences observed between the two strains.

David Adams, Julia Newton-Bishop, Timothy Bishop, Nicholas Hayward and colleagues identify loss-of-function variants in POT1 in several families with early onset multiple primary melanoma. They further show that these variants disrupt telomere binding by POT1 and are associated with increased telomere length. Deleterious germline variants in CDKN2A account for around 40% of familial melanoma cases1, and rare variants in CDK4, BRCA2, BAP1 and the promoter of TERT have also been linked to the disease2,3,4,5. Here we set out to identify new high-penetrance susceptibility genes by sequencing 184 melanoma cases from 105 pedigrees recruited in the UK, The Netherlands and Australia that were negative for variants in known predisposition genes. We identified families where melanoma cosegregates with loss-of-function variants in the protection of telomeres 1 gene (POT1), with a proportion of family members presenting with an early age of onset and multiple primary tumors. We show that these variants either affect POT1 mRNA splicing or alter key residues in the highly conserved oligonucleotide/oligosaccharide-binding (OB) domains of POT1, disrupting protein-telomere binding and leading to increased telomere length. These findings suggest that POT1 variants predispose to melanoma formation via a direct effect on telomeres.

An in vivo transposon screen in a pancreatic cancer model identifies frequent inactivation of Usp9x; deletion of Usp9x cooperates with KrasG12D to accelerate rapidly pancreatic tumorigenesis in mice, validating their genetic interaction. An insertional-mutagenesis screen for genes that cooperate with oncogenic KrasG12D in a mouse model of pancreatic cancer has identified the disruption of Usp9x in more than half of the almost 200 distinct tumours analysed. This is surprising, because previous work had indicated that the deubiquitinase USP9X promoted tumour-cell survival. Functional studies show that loss of Usp9x protects pancreatic cancer cells against apoptosis through an effect on multiple oncogenic signalling pathways. In human pancreatic cancers, Usp9x expression is often lost and this correlates with more advanced disease. Approaches that can trigger re-expression of Usp9x may therefore be useful in the treatment of pancreatic cancer. Pancreatic ductal adenocarcinoma (PDA) remains a lethal malignancy despite much progress concerning its molecular characterization. PDA tumours harbour four signature somatic mutations1,2,3,4 in addition to numerous lower frequency genetic events of uncertain significance5. Here we use Sleeping Beauty (SB) transposon-mediated insertional mutagenesis6,7 in a mouse model of pancreatic ductal preneoplasia8 to identify genes that cooperate with oncogenic KrasG12D to accelerate tumorigenesis and promote progression. Our screen revealed new candidate genes for PDA and confirmed the importance of many genes and pathways previously implicated in human PDA. The most commonly mutated gene was the X-linked deubiquitinase Usp9x, which was inactivated in over 50% of the tumours. Although previous work had attributed a pro-survival role to USP9X in human neoplasia9, we found instead that loss of Usp9x enhances transformation and protects pancreatic cancer cells from anoikis. Clinically, low USP9X protein and messenger RNA expression in PDA correlates with poor survival after surgery, and USP9X levels are inversely associated with metastatic burden in advanced disease. Furthermore, chromatin modulation with trichostatin A or 5-aza-2′-deoxycytidine elevates USP9X expression in human PDA cell lines, indicating a clinical approach for certain patients. The conditional deletion of Usp9x cooperated with KrasG12D to accelerate pancreatic tumorigenesis in mice, validating their genetic interaction. We propose that USP9X is a major tumour suppressor gene with prognostic and therapeutic relevance in PDA.