Deep genome sequencing and variation analysis of 13 inbred mouse strains defines candidate phenotypic alleles, private variation and homozygous truncating mutations

Computational identification of disease models through cross-species phenotype comparison

The use of standardised phenotyping screens to identify abnormal phenotypes in mouse knockouts, together with the use of ontologies to describe such phenotypic features, allows the implementation of an automated and unbiased pipeline to identify new models of disease by performing phenotype comparisons across species. Using data from the International Mouse Phenotyping Consortium (IMPC), approximately half of mouse mutants are able to mimic, at least partially, the human ortholog disease phenotypes as computed by the PhenoDigm algorithm. We found the number of phenotypic abnormalities in the mouse and the corresponding Mendelian disorder, the pleiotropy and severity of the disease, and the viability and zygosity status of the mouse knockout to be associated with the ability of mouse models to recapitulate the human disorder. An analysis of the IMPC impact on disease gene discovery through a publication-tracking system revealed that the resource has been implicated in at least 109 validated rare disease-gene associations over the last decade.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation - a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

Insertion of short L1 sequences generates inter-strain histone acetylation differences in the mouse

BACKGROUND

Gene expression divergence between populations and between individuals can emerge from genetic variations within the genes and/or in the cis regulatory elements. Since epigenetic modifications regulate gene expression, it is conceivable that epigenetic variations in cis regulatory elements can also be a source of gene expression divergence.

RESULTS

In this study, we compared histone acetylation (namely, H3K9ac) profiles in two mouse strains of different subspecies origin, C57BL/6 J (B6) and MSM/Ms (MSM), as well as their F1 hybrids. This identified 319 regions of strain-specific acetylation, about half of which were observed between the alleles of F1 hybrids. While the allele-specific presence of the interferon regulatory factor 3 (IRF3) binding sequence was associated with allele-specific histone acetylation, we also revealed that B6-specific insertions of a short 3' fragment of LINE-1 (L1) retrotransposon occur within or proximal to MSM-specific acetylated regions. Furthermore, even in hyperacetylated domains, flanking regions of non-polymorphic 3' L1 fragments were hypoacetylated, suggesting a general activity of the 3' L1 fragment to induce hypoacetylation. Indeed, we confirmed the binding of the 3' region of L1 by three Krüppel-associated box domain-containing zinc finger proteins (KZFPs), which interact with histone deacetylases. These results suggest that even a short insertion of L1 would be excluded from gene- and acetylation-rich regions by natural selection. Finally, mRNA-seq analysis for F1 hybrids was carried out, which disclosed a link between allele-specific promoter/enhancer acetylation and gene expression.

CONCLUSIONS

This study disclosed a number of genetic changes that have changed the histone acetylation levels during the evolution of mouse subspecies, a part of which is associated with gene expression changes. Insertions of even a very short L1 fragment can decrease the acetylation level in their neighboring regions and thereby have been counter-selected in gene-rich regions, which may explain a long-standing mystery of discrete genomic distribution of LINEs and SINEs.

Single-mitosis dissection of acute and chronic DNA mutagenesis and repair

How chronic mutational processes and punctuated bursts of DNA damage drive evolution of the cancer genome is poorly understood. Here, we demonstrate a strategy to disentangle and quantify distinct mechanisms underlying genome evolution in single cells, during single mitoses and at single-strand resolution. To distinguish between chronic (reactive oxygen species (ROS)) and acute (ultraviolet light (UV)) mutagenesis, we microfluidically separate pairs of sister cells from the first mitosis following burst UV damage. Strikingly, UV mutations manifest as sister-specific events, revealing mirror-image mutation phasing genome-wide. In contrast, ROS mutagenesis in transcribed regions is reduced strand agnostically. Successive rounds of genome replication over persisting UV damage drives multiallelic variation at CC dinucleotides. Finally, we show that mutation phasing can be resolved to single strands across the entire genome of liver tumors from F1 mice. This strategy can be broadly used to distinguish the contributions of overlapping cancer relevant mutational processes.

Whole‑exome sequencing reveals Lewis lung carcinoma is a hypermutated Kras/Nras-mutant cancer with extensive regional mutation clusters in its genome

Lewis lung carcinoma (LLC), as a widely used preclinical cancer model, has still not been genetically and genomically characterized. Here, we performed a whole-exome sequencing analysis on the LLC cell line to elucidate its molecular characteristics and etiologies. Our data showed that LLC originated from a male mouse belonging to C57BL/6L (a transitional strain between C57BL/6J and C57BL/6N) and contains substantial somatic SNV and InDel mutations (> 20,000). Extensive regional mutation clusters are present in its genome, which were caused mainly by the mutational processes underlying the SBS1, SBS5, SBS15, SBS17a, and SBS21 signatures during frequent structural rearrangements. Thirty three deleterious mutations are present in 30 cancer genes including Kras, Nras, Trp53, Dcc, and Cacna1d. Cdkn2a and Cdkn2b are biallelically deleted from the genome. Five pathways (RTK/RAS, p53, cell cycle, TGFB, and Hippo) are oncogenically deregulated or affected. The major mutational processes in LLC include chromosomal instability, exposure to metabolic mutagens, spontaneous 5-methylcytosine deamination, defective DNA mismatch repair, and reactive oxygen species. Our data also suggest that LLC is a lung cancer similar to human lung adenocarcinoma. This study lays a molecular basis for the more targeted application of LLC in preclinical research.

Evaluation of genetic demultiplexing of single-cell sequencing data from model species

Single-cell sequencing (sc-seq) provides a species agnostic tool to study cellular processes. However, these technologies are expensive and require sufficient cell quantities and biological replicates to avoid artifactual results. An option to address these problems is pooling cells from multiple individuals into one sc-seq library. In humans, genotype-based computational separation (i.e., demultiplexing) of pooled sc-seq samples is common. This approach would be instrumental for studying non-isogenic model organisms. We set out to determine whether genotype-based demultiplexing could be more broadly applied among species ranging from zebrafish to non-human primates. Using such non-isogenic species, we benchmark genotype-based demultiplexing of pooled sc-seq datasets against various ground truths. We demonstrate that genotype-based demultiplexing of pooled sc-seq samples can be used with confidence in several non-isogenic model organisms and uncover limitations of this method. Importantly, the only genomic resource required for this approach is sc-seq data and a de novo transcriptome. The incorporation of pooling into sc-seq study designs will decrease cost while simultaneously increasing the reproducibility and experimental options in non-isogenic model organisms.

Resolution of structural variation in diverse mouse genomes reveals chromatin remodeling due to transposable elements

Diverse inbred mouse strains are important biomedical research models, yet genome characterization of many strains is fundamentally lacking in comparison with humans. In particular, catalogs of structural variants (SVs) (variants ≥ 50 bp) are incomplete, limiting the discovery of causative alleles for phenotypic variation. Here, we resolve genome-wide SVs in 20 genetically distinct inbred mice with long-read sequencing. We report 413,758 site-specific SVs affecting 13% (356 Mbp) of the mouse reference assembly, including 510 previously unannotated coding variants. We substantially improve the Mus musculus transposable element (TE) callset, and we find that TEs comprise 39% of SVs and account for 75% of altered bases. We further utilize this callset to investigate how TE heterogeneity affects mouse embryonic stem cells and find multiple TE classes that influence chromatin accessibility. Our work provides a comprehensive analysis of SVs found in diverse mouse genomes and illustrates the role of TEs in epigenetic differences.

Knockout mice are an important tool for human monogenic heart disease studies

Mouse models are relevant to studying the functionality of genes involved in human diseases; however, translation of phenotypes can be challenging. Here, we investigated genes related to monogenic forms of cardiovascular disease based on the Genomics England PanelApp and aligned them to International Mouse Phenotyping Consortium (IMPC) data. We found 153 genes associated with cardiomyopathy, cardiac arrhythmias or congenital heart disease in humans, of which 151 have one-to-one mouse orthologues. For 37.7% (57/151), viability and heart data captured by electrocardiography, transthoracic echocardiography, morphology and pathology from embryos and young adult mice are available. In knockout mice, 75.4% (43/57) of these genes showed non-viable phenotypes, whereas records of prenatal, neonatal or infant death in humans were found for 35.1% (20/57). Multisystem phenotypes are common, with 58.8% (20/34) of heterozygous (homozygous lethal) and 78.6% (11/14) of homozygous (viable) mice showing cardiovascular, metabolic/homeostasis, musculoskeletal, hematopoietic, nervous system and/or growth abnormalities mimicking the clinical manifestations observed in patients. These IMPC data are critical beyond cardiac diagnostics given their multisystemic nature, allowing detection of abnormalities across physiological systems and providing a valuable resource to understand pleiotropic effects.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Pisani L, Yoo B, Bejerano G, Peltz G. Analysis of structural variation among inbred mouse strains

BACKGROUND

'Long read' sequencing methods have been used to identify previously uncharacterized structural variants that cause human genetic diseases. Therefore, we investigated whether long read sequencing could facilitate genetic analysis of murine models for human diseases.

RESULTS

The genomes of six inbred strains (BTBR T + Itpr3tf/J, 129Sv1/J, C57BL/6/J, Balb/c/J, A/J, SJL/J) were analyzed using long read sequencing. Our results revealed that (i) Structural variants are very abundant within the genome of inbred strains (4.8 per gene) and (ii) that we cannot accurately infer whether structural variants are present using conventional short read genomic sequence data, even when nearby SNP alleles are known. The advantage of having a more complete map was demonstrated by analyzing the genomic sequence of BTBR mice. Based upon this analysis, knockin mice were generated and used to characterize a BTBR-unique 8-bp deletion within Draxin that contributes to the BTBR neuroanatomic abnormalities, which resemble human autism spectrum disorder.

CONCLUSION

A more complete map of the pattern of genetic variation among inbred strains, which is produced by long read genomic sequencing of the genomes of additional inbred strains, could facilitate genetic discovery when murine models of human diseases are analyzed.

Quantitative trait loci on chromosomes 9 and 19 modulate AII amacrine cell number in the mouse retina

Sequence variants modulating gene function or expression affect various heritable traits, including the number of neurons within a population. The present study employed a forward-genetic approach to identify candidate causal genes and their sequence variants controlling the number of one type of retinal neuron, the AII amacrine cell. Data from twenty-six recombinant inbred (RI) strains of mice derived from the parental C57BL/6J (B6/J) and A/J laboratory strains were used to identify genomic loci regulating cell number. Large variation in cell number is present across the RI strains, from a low of ∼57,000 cells to a high of ∼87,000 cells. Quantitative trait locus (QTL) analysis revealed three prospective controlling genomic loci, on Chromosomes (Chrs) 9, 11, and 19, each contributing additive effects that together approach the range of variation observed. Composite interval mapping validated two of these loci, and chromosome substitution strains, in which the A/J genome for Chr 9 or 19 was introgressed on a B6/J genetic background, showed increased numbers of AII amacrine cells as predicted by those two QTL effects. Analysis of the respective genomic loci identified candidate controlling genes defined by their retinal expression, their established biological functions, and by the presence of sequence variants expected to modulate gene function or expression. Two candidate genes, Dtx4 on Chr 19, being a regulator of Notch signaling, and Dixdc1 on Chr 9, a modulator of the WNT-β-catenin signaling pathway, were explored in further detail. Postnatal overexpression of Dtx4 was found to reduce the frequency of amacrine cells, while Dixdc1 knockout retinas contained an excess of AII amacrine cells. Sequence variants in each gene were identified, being the likely sources of variation in gene expression, ultimately contributing to the final number of AII amacrine cells.

E96V Mutation in the Kdelr3 Gene Is Associated with Type 2 Diabetes Susceptibility in Obese NZO Mice

Type 2 diabetes (T2D) represents a multifactorial metabolic disease with a strong genetic predisposition. Despite elaborate efforts in identifying the genetic variants determining individual susceptibility towards T2D, the majority of genetic factors driving disease development remain poorly understood. With the aim to identify novel T2D risk genes we previously generated an N2 outcross population using the two inbred mouse strains New Zealand obese (NZO) and C3HeB/FeJ (C3H). A linkage study performed in this population led to the identification of the novel T2D-associated quantitative trait locus (QTL) Nbg15 (NZO blood glucose on chromosome 15, Logarithm of odds (LOD) 6.6). In this study we used a combined approach of positional cloning, gene expression analyses and in silico predictions of DNA polymorphism on gene/protein function to dissect the genetic variants linking Nbg15 to the development of T2D. Moreover, we have generated congenic strains that associated the distal sublocus of Nbg15 to mechanisms altering pancreatic beta cell function. In this sublocus, Cbx6, Fam135b and Kdelr3 were nominated as potential causative genes associated with the Nbg15 driven effects. Moreover, a putative mutation in the Kdelr3 gene from NZO was identified, negatively influencing adaptive responses associated with pancreatic beta cell death and induction of endoplasmic reticulum stress. Importantly, knockdown of Kdelr3 in cultured Min6 beta cells altered insulin granules maturation and pro-insulin levels, pointing towards a crucial role of this gene in islets function and T2D susceptibility.

Neurog2 regulates Isl1 to modulate horizontal cell number

The population sizes of different retinal cell types vary between different strains of mice, and that variation can be mapped to genomic loci in order to identify its polygenic origin. In some cases, controlling genes act independently, whereas in other instances, they exhibit epistasis. Here, we identify an epistatic interaction revealed through the mapping of quantitative trait loci from a panel of recombinant inbred strains of mice. The population of retinal horizontal cells exhibits a twofold variation in number, mapping to quantitative trait loci on chromosomes 3 and 13, where these loci are shown to interact epistatically. We identify a prospective genetic interaction underlying this, mediated by the bHLH transcription factor Neurog2, at the chromosome 3 locus, functioning to repress the LIM homeodomain transcription factor Isl1, at the chromosome 13 locus. Using single and double conditional knockout mice, we confirm the countervailing actions of each gene, and validate in vitro a crucial role for two single nucleotide polymorphisms in the 5'UTR of Isl1, one of which yields a novel E-box, mediating the repressive action of Neurog2.

Forward genetic screen of homeostatic antibody levels in the Collaborative Cross identifies MBD1 as a novel regulator of B cell homeostasis

Variation in immune homeostasis, the state in which the immune system is maintained in the absence of stimulation, is highly variable across populations. This variation is attributed to both genetic and environmental factors. However, the identity and function of specific regulators have been difficult to identify in humans. We evaluated homeostatic antibody levels in the serum of the Collaborative Cross (CC) mouse genetic reference population. We found heritable variation in all antibody isotypes and subtypes measured. We identified 4 quantitative trait loci (QTL) associated with 3 IgG subtypes: IgG1, IgG2b, and IgG2c. While 3 of these QTL map to genome regions of known immunological significance (major histocompatibility and immunoglobulin heavy chain locus), Qih1 (associated with variation in IgG1) mapped to a novel locus on Chromosome 18. We further associated this locus with B cell proportions in the spleen and identify Methyl-CpG binding domain protein 1 under this locus as a novel regulator of homeostatic IgG1 levels in the serum and marginal zone B cells (MZB) in the spleen, consistent with a role in MZB differentiation to antibody secreting cells.

Isometric skeletal muscle contractile properties in common strains of male laboratory mice

Assessing contractile function of skeletal muscle in murine models is a commonly employed laboratory technique that investigators utilize to measure the impact of genetic manipulations, drug efficacy, or other therapeutic interventions. Often overlooked is the potential for the strain of the mouse to influence the functional properties of the skeletal muscle. Thus, we sought to characterize commonly assessed isometric force measures in the hindlimb muscles across a variety of mouse strains. Using 6-8-week-old male mice, we measured isometric force, fatigue susceptibility, relaxation kinetics, muscle mass, myofiber cross-sectional area, and fiber type composition of the extensor digitorum longus (EDL) and soleus muscles in C57BL/6NJ, BALB/cJ, FVB/NJ, C57BL/6J, and C57BL/10 mice. The data demonstrate both unique differences and a number of similarities between both muscles in the various genetic backgrounds of mice. Soleus muscle specific force (i.e., force per unit size) exhibited higher variation across strains while specific force of the EDL muscle exhibited minimal variation. In contrast, absolute force differed only in a few mouse strains whereas analysis of muscle morphology revealed many distinctions when compared across all the groups. Collectively, the data suggest that the strain of the mouse can potentially influence the measured biological outcome and may possibly promote a synergistic effect with any genetic manipulation or therapeutic intervention. Thus, it is critical for the investigator to carefully consider the genetic background of the mouse used in the experimental design and precisely document the strain of mouse employed during publication.

Pierre CL, Williams A, Sohni A, Wilkinson MF, Gymrek M, Sebat J, Palmer AA. SNPs, short tandem repeats, and structural variants are responsible for differential gene expression across C57BL/6 and C57BL/10 substrains

Mouse substrains are an invaluable model for understanding disease. We compared C57BL/6J, which is the most commonly used inbred mouse strain, with eight C57BL/6 and five C57BL/10 closely related inbred substrains. Whole-genome sequencing and RNA-sequencing analysis yielded 352,631 SNPs, 109,096 indels, 150,344 short tandem repeats (STRs), 3,425 structural variants (SVs), and 2,826 differentially expressed genes (DE genes) among these 14 strains; 312,981 SNPs (89%) distinguished the B6 and B10 lineages. These SNPs were clustered into 28 short segments that are likely due to introgressed haplotypes rather than new mutations. Outside of these introgressed regions, we identified 53 SVs, protein-truncating SNPs, and frameshifting indels that were associated with DE genes. Our results can be used for both forward and reverse genetic approaches and illustrate how introgression and mutational processes give rise to differences among these widely used inbred substrains.

Genomic variation in captive deer mouse (Peromyscus maniculatus) populations

BACKGROUND

Deer mice (genus Peromyscus) are the most common rodents in North America. Despite the availability of reference genomes for some species, a comprehensive database of polymorphisms, especially in those maintained as living stocks and distributed to academic investigators, is missing. In the present study we surveyed two populations of P. maniculatus that are maintained at the Peromyscus Genetic Stock Center (PGSC) for polymorphisms across their 2.5 × 109 bp genome.

RESULTS

High density of variation was identified, corresponding to one SNP every 55 bp for the high altitude stock (SM2) or 207 bp for the low altitude stock (BW) using snpEff (v4.3). Indels were detected every 1157 bp for BW or 311 bp for SM2. The average Watterson estimator for the BW and SM2 populations is 248813.70388 and 869071.7671 respectively. Some differences in the distribution of missense, nonsense and silent mutations were identified between the stocks, as well as polymorphisms in genes associated with inflammation (NFATC2), hypoxia (HIF1a) and cholesterol metabolism (INSIG1) and may possess value in modeling pathology.

CONCLUSIONS

This genomic resource, in combination with the availability of P. maniculatus from the PGSC, is expected to promote genetic and genomic studies with this animal model.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Background

Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines.

Results

In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters.

Conclusions

Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07577-3.

Chd8 regulates X chromosome inactivation in mouse through fine-tuning control of Xist expression

Female mammals achieve dosage compensation by inactivating one of their two X chromosomes during development, a process entirely dependent on Xist, an X-linked long non-coding RNA (lncRNA). At the onset of X chromosome inactivation (XCI), Xist is up-regulated and spreads along the future inactive X chromosome. Contextually, it recruits repressive histone and DNA modifiers that transcriptionally silence the X chromosome. Xist regulation is tightly coupled to differentiation and its expression is under the control of both pluripotency and epigenetic factors. Recent evidence has suggested that chromatin remodelers accumulate at the X Inactivation Center (XIC) and here we demonstrate a new role for Chd8 in Xist regulation in differentiating ES cells, linked to its control and prevention of spurious transcription factor interactions occurring within Xist regulatory regions. Our findings have a broader relevance, in the context of complex, developmentally-regulated gene expression.

Andrea Cerase et al. report that the chromatin remodeler Chd8 is a key regulator of mammalian Xist expression and therefore X chromosome inactivation. They find that Chd8 activates Xist expression in embryonic stem cells, while in differentiating cells it acts to prevent spurious Xist expression through blocking transcription factor binding to the Xist promoter.

In silico candidate variant and gene identification using inbred mouse strains

Mice are the most widely used animal model to study genotype to phenotype relationships. Inbred mice are genetically identical, which eliminates genetic heterogeneity and makes them particularly useful for genetic studies. Many different strains have been bred over decades and a vast amount of phenotypic data has been generated. In addition, recently whole genome sequencing-based genome-wide genotype data for many widely used inbred strains has been released. Here, we present an approach for in silico fine-mapping that uses genotypic data of 37 inbred mouse strains together with phenotypic data provided by the user to propose candidate variants and genes for the phenotype under study. Public genome-wide genotype data covering more than 74 million variant sites is queried efficiently in real-time to provide those variants that are compatible with the observed phenotype differences between strains. Variants can be filtered by molecular consequences and by corresponding molecular impact. Candidate gene lists can be generated from variant lists on the fly. Fine-mapping together with annotation or filtering of results is provided in a Bioconductor package called MouseFM. In order to characterize candidate variant lists under various settings, MouseFM was applied to two expression data sets across 20 inbred mouse strains, one from neutrophils and one from CD4⁺ T cells. Fine-mapping was assessed for about 10,000 genes, respectively, and identified candidate variants and haplotypes for many expression quantitative trait loci (eQTLs) reported previously based on these data. For albinism, MouseFM reports only one variant allele of moderate or high molecular impact that only albino mice share: a missense variant in the Tyr gene, reported previously to be causal for this phenotype. Performing in silico fine-mapping for interfrontal bone formation in mice using four strains with and five strains without interfrontal bone results in 12 genes. Of these, three are related to skull shaping abnormality. Finally performing fine-mapping for dystrophic cardiac calcification by comparing 9 strains showing the phenotype with eight strains lacking it, we identify only one moderate impact variant in the known causal gene Abcc6. In summary, this illustrates the benefit of using MouseFM for candidate variant and gene identification.

The genetic background significantly impacts the severity of kidney cystic disease in the Pkd1RC/RC mouse model of autosomal dominant polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD), primarily due to PKD1 or PKD2 mutations, causes progressive kidney cyst development and kidney failure. There is significant intrafamilial variability likely due to the genetic background and environmental/lifestyle factors; variability that can be modeled in PKD mice. Here, we characterized mice homozygous for the PKD1 hypomorphic allele, p.Arg3277Cys (Pkd1^RC/RC), inbred into the BALB/cJ (BC) or the 129S6/SvEvTac (129) strains, plus F1 progeny bred with the previously characterized C57BL/6J (B6) model; F1(BC/B6) or F1(129/B6). By one-month cystic disease in both the BC and 129 Pkd1^RC/RC mice was more severe than in B6 and continued with more rapid progression to six to nine months. Thereafter, the expansive disease stage plateaued/declined, coinciding with increased fibrosis and a clear decline in kidney function. Greater severity correlated with more inter-animal and inter-kidney disease variability, especially in the 129-line. Both F1 combinations had intermediate disease severity, more similar to B6 but progressive from one-month of age. Mild biliary dysgenesis, and an early switch from proximal tubule to collecting duct cysts, was seen in all backgrounds. Preclinical testing with a positive control, tolvaptan, employed the F1(129/B6)-Pkd1^RC/RC line, which has moderately progressive disease and limited isogenic variability. Magnetic resonance imaging was utilized to randomize animals and provide total kidney volume endpoints; complementing more traditional data. Thus, we show how genetic background can tailor the Pkd1^RC/RC model to address different aspects of pathogenesis and disease modification, and describe a possible standardized protocol for preclinical testing.

Content and Performance of the MiniMUGA Genotyping Array: A New Tool To Improve Rigor and Reproducibility in Mouse Research

The laboratory mouse is the most widely used animal model for biomedical research, due in part to its well-annotated genome, wealth of genetic resources, and the ability to precisely manipulate its genome. Despite the importance of genetics for mouse research, genetic quality control (QC) is not standardized, in part due to the lack of cost-effective, informative, and robust platforms. Genotyping arrays are standard tools for mouse research and remain an attractive alternative even in the era of high-throughput whole-genome sequencing. Here, we describe the content and performance of a new iteration of the Mouse Universal Genotyping Array (MUGA), MiniMUGA, an array-based genetic QC platform with over 11,000 probes. In addition to robust discrimination between most classical and wild-derived laboratory strains, MiniMUGA was designed to contain features not available in other platforms: (1) chromosomal sex determination, (2) discrimination between substrains from multiple commercial vendors, (3) diagnostic SNPs for popular laboratory strains, (4) detection of constructs used in genetically engineered mice, and (5) an easy-to-interpret report summarizing these results. In-depth annotation of all probes should facilitate custom analyses by individual researchers. To determine the performance of MiniMUGA, we genotyped 6899 samples from a wide variety of genetic backgrounds. The performance of MiniMUGA compares favorably with three previous iterations of the MUGA family of arrays, both in discrimination capabilities and robustness. We have generated publicly available consensus genotypes for 241 inbred strains including classical, wild-derived, and recombinant inbred lines. Here, we also report the detection of a substantial number of XO and XXY individuals across a variety of sample types, new markers that expand the utility of reduced complexity crosses to genetic backgrounds other than C57BL/6, and the robust detection of 17 genetic constructs. We provide preliminary evidence that the array can be used to identify both partial sex chromosome duplication and mosaicism, and that diagnostic SNPs can be used to determine how long inbred mice have been bred independently from the relevant main stock. We conclude that MiniMUGA is a valuable platform for genetic QC, and an important new tool to increase the rigor and reproducibility of mouse research.

Variation under domestication in animal models: the case of the Mexican axolotl

BACKGROUND

Species adaptation to laboratory conditions is a special case of domestication that has modified model organisms phenotypically and genetically. The characterisation of these changes is crucial to understand how this variation can affect the outcome of biological experiments. Yet despite the wide use of laboratory animals in biological research, knowledge of the genetic diversity within and between different strains and populations of some animal models is still scarce. This is particularly the case of the Mexican axolotl, which has been bred in captivity since 1864.

RESULTS

Using gene expression data from nine different projects, nucleotide sequence variants were characterised, and distinctive genetic background of the experimental specimens was uncovered. This study provides a catalogue of thousands of nucleotide variants along predicted protein-coding genes, while identifying genome-wide differences between pigment phenotypes in laboratory populations.

CONCLUSIONS

Awareness of the genetic variation could guide a better experimental design while helping to develop molecular tools for monitoring genetic diversity and studying gene functions in laboratory axolotls. Overall, this study highlights the cross-taxa utility that transcriptomic data might have to assess the genetic variation of the experimental specimens, which might help to shorten the journey towards reproducible research.

Functional signatures of evolutionarily young CTCF binding sites

Background

The introduction of novel CTCF binding sites in gene regulatory regions in the rodent lineage is partly the effect of transposable element expansion, particularly in the murine lineage. The exact mechanism and functional impact of evolutionarily novel CTCF binding sites are not yet fully understood. We investigated the impact of novel subspecies-specific CTCF binding sites in two Mus genus subspecies, Mus musculus domesticus and Mus musculus castaneus, that diverged 0.5 million years ago.

Results

CTCF binding site evolution is influenced by the action of the B2-B4 family of transposable elements independently in both lineages, leading to the proliferation of novel CTCF binding sites. A subset of evolutionarily young sites may harbour transcriptional functionality as evidenced by the stability of their binding across multiple tissues in M. musculus domesticus (BL6), while overall the distance of subspecies-specific CTCF binding to the nearest transcription start sites and/or topologically associated domains (TADs) is largely similar to musculus-common CTCF sites. Remarkably, we discovered a recurrent regulatory architecture consisting of a CTCF binding site and an interferon gene that appears to have been tandemly duplicated to create a 15-gene cluster on chromosome 4, thus forming a novel BL6 specific immune locus in which CTCF may play a regulatory role.

Conclusions

Our results demonstrate that thousands of CTCF binding sites show multiple functional signatures rapidly after incorporation into the genome.

Transgenerational epigenetic effects from male exposure to endocrine-disrupting compounds: a systematic review on research in mammals

Assessing long-term health effects from a potentially harmful environment is challenging. Endocrine-disrupting compounds (EDCs) have become omnipresent in our environment. Individuals may or may not experience clinical health issues from being exposed to the increasing environmental pollution in daily life, but an issue of high concern is that also the non-exposed progeny may encounter consequences of these ancestral exposures. Progress in understanding epigenetic mechanisms opens new perspectives to estimate the risk of man-made EDCs. However, the field of epigenetic toxicology is new and its application in public health or in the understanding of disease etiology is almost non-existent, especially if it concerns future generations. In this review, we investigate the literature on transgenerational inheritance of diseases, published in the past 10 years. We question whether persistent epigenetic changes occur in the male germ line after exposure to synthesized EDCs. Our systematic search led to an inclusion of 43 articles, exploring the effects of commonly used synthetic EDCs, such as plasticizers (phthalates and bisphenol A), pesticides (dichlorodiphenyltrichloroethane, atrazine, vinclozin, methoxychlor), dioxins, and polycyclic aromatic hydrocarbons (PAHs, such as benzo(a)pyrene). Most studies found transgenerational epigenetic effects, often linked to puberty- or adult-onset diseases, such as testicular or prostate abnormalities, metabolic disorders, behavioral anomalies, and tumor development. The affected epigenetic mechanisms included changes in DNA methylation patterns, transcriptome, and expression of DNA methyltransferases. Studies involved experiments in animal models and none were based on human data. In the future, human studies are needed to confirm animal findings. If not transgenerational, at least intergenerational human studies and studies on EDC-induced epigenetic effects on germ cells could help to understand early processes of inheritance. Next, toxicity tests of new chemicals need a more comprehensive approach before they are introduced on the market. We further point to the relevance of epigenetic toxicity tests in regard to public health of the current population but also of future generations. Finally, this review sheds a light on how the interplay of genetics and epigenetics may explain the current knowledge gap on transgenerational inheritance.

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.

An Automated Method To Predict Mouse Gene and Protein Sequences Using Variant Data

With recent advances in sequencing technologies, the scientific community has begun to probe the potential genetic bases behind complex phenotypes in humans and model organisms. In many cases, the genomes of genetically distinct strains of model organisms, such as the mouse (Mus musculus), have not been fully sequenced. Here, we report on a tool designed to use single-nucleotide polymorphism (SNP) and insertion-deletion (indel) data to predict gene, mRNA, and protein sequences for up to 36 genetically distinct mouse strains. By automated querying of freely accessible databases through a graphical interface, the software requires no data and little computational experience. As a proof of concept, we predicted the gene and amino acid sequence of the aryl hydrocarbon receptor (Ahr) for all inbred mouse strains of which variant data were currently available through Mouse Genome Project. Predicted sequences were compared with fully sequenced genomes to show that the tool is effective in predicting gene and protein sequences.

In situ CRISPR-Cas9 base editing for the development of genetically engineered mouse models of breast cancer

Genetically engineered mouse models (GEMMs) of cancer have proven to be of great value for basic and translational research. Although CRISPR-based gene disruption offers a fast-track approach for perturbing gene function and circumvents certain limitations of standard GEMM development, it does not provide a flexible platform for recapitulating clinically relevant missense mutations in vivo. To this end, we generated knock-in mice with Cre-conditional expression of a cytidine base editor and tested their utility for precise somatic engineering of missense mutations in key cancer drivers. Upon intraductal delivery of sgRNA-encoding vectors, we could install point mutations with high efficiency in one or multiple endogenous genes in situ and assess the effect of defined allelic variants on mammary tumorigenesis. While the system also produces bystander insertions and deletions that can stochastically be selected for when targeting a tumor suppressor gene, we could effectively recapitulate oncogenic nonsense mutations. We successfully applied this system in a model of triple-negative breast cancer, providing the proof of concept for extending this flexible somatic base editing platform to other tissues and tumor types.

Population Genomics of the House Mouse and the Brown Rat

Mice (Mus musculus) and rats (Rattus norvegicus) have long served as model systems for biomedical research. However, they are also excellent models for studying the evolution of populations, subspecies, and species. Within the past million years, they have spread in various waves across large parts of the globe, with the most recent spread in the wake of human civilization. They have developed into commensal species, but have also been able to colonize extreme environments on islands free of human civilization. Given that ample genomic and genetic resources are available for these species, they have thus also become ideal mammalian systems for evolutionary studies on adaptation and speciation, particularly in the combination with the rapid developments in population genomics. The chapter provides an overview of the systems and their history, as well as of available resources.

Immune Relevant and Immune Deficient Mice: Options and Opportunities in Translational Research

In 1989 ILAR published a list and description of immunodeficient rodents used in research. Since then, advances in understanding of molecular mechanisms; recognition of genetic, epigenetic microbial, and other influences on immunity; and capabilities in manipulating genomes and microbiomes have increased options and opportunities for selecting mice and designing studies to answer important mechanistic and therapeutic questions. Despite numerous scientific breakthroughs that have benefitted from research in mice, there is debate about the relevance and predictive or translational value of research in mice. Reproducibility of results obtained from mice and other research models also is a well-publicized concern. This review summarizes resources to inform the selection and use of immune relevant mouse strains and stocks, aiming to improve the utility, validity, and reproducibility of research in mice. Immune sufficient genetic variations, immune relevant spontaneous mutations, immunodeficient and autoimmune phenotypes, and selected induced conditions are emphasized.

Mouse protein coding diversity: What's left to discover?

For over a century, mice have been used to model human disease, leading to many fundamental discoveries about mammalian biology and the development of new therapies. Mouse genetics research has been further catalysed by a plethora of genomic resources developed in the last 20 years, including the genome sequence of C57BL/6J and more recently the first draft reference genomes for 16 additional laboratory strains. Collectively, the comparison of these genomes highlights the extreme diversity that exists at loci associated with the immune system, pathogen response, and key sensory functions, which form the foundation for dissecting phenotypic traits in vivo. We review the current status of the mouse genome across the diversity of the mouse lineage and discuss the value of mice to understanding human disease.

Drug-induced chromatin accessibility changes associate with sensitivity to liver tumor promotion

Liver cancer susceptibility varies amongst humans and between experimental animal models because of multiple genetic and epigenetic factors. The molecular characterization of such susceptibilities has the potential to enhance cancer risk assessment of xenobiotic exposures and disease prevention strategies. Here, using DNase I hypersensitivity mapping coupled with transcriptomic profiling, we investigate perturbations in cis-acting gene regulatory elements associated with the early stages of phenobarbital (PB)-mediated liver tumor promotion in susceptible versus resistant mouse strains (B6C3F1 versus C57BL/6J). Integrated computational analyses of strain-selective changes in liver chromatin accessibility underlying PB response reveal differential epigenetic regulation of molecular pathways associated with PB-mediated tumor promotion, including Wnt/β-catenin signaling. Complementary transcription factor motif analyses reveal mouse strain-selective gene regulatory networks and a novel role for Stat, Smad, and Fox transcription factors in the early stages of PB-mediated tumor promotion. Mapping perturbations in cis-acting gene regulatory elements provides novel insights into the molecular basis for susceptibility to xenobiotic-induced rodent liver tumor promotion and has the potential to enhance mechanism-based cancer risk assessments of xenobiotic exposures.

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.

Gene-diet interactions associated with complex trait variation in an advanced intercross outbred mouse line

Phenotypic variation of quantitative traits is orchestrated by a complex interplay between the environment (e.g. diet) and genetics. However, the impact of gene-environment interactions on phenotypic traits mostly remains elusive. To address this, we feed 1154 mice of an autoimmunity-prone intercross line (AIL) three different diets. We find that diet substantially contributes to the variability of complex traits and unmasks additional genetic susceptibility quantitative trait loci (QTL). By performing whole-genome sequencing of the AIL founder strains, we resolve these QTLs to few or single candidate genes. To address whether diet can also modulate genetic predisposition towards a given trait, we set NZM2410/J mice on similar dietary regimens as AIL mice. Our data suggest that diet modifies genetic susceptibility to lupus and shifts intestinal bacterial and fungal community composition, which precedes clinical disease manifestation. Collectively, our study underlines the importance of including environmental factors in genetic association studies.

Complex traits associate with genetic variation and environment and their interaction. Here, the authors study the influence of different diets on trait variability in 1154 outbred mice from an advanced intercross line and find gene-diet interactions associated with spontaneous autoimmunity development in these animals.

Different effects of titanium dioxide nanoparticles instillation in young and adult mice on DNA methylation related with lung inflammation and fibrosis

Wide use of titanium dioxide nanoparticles (TiO₂ NPs) as white pigments induces unintentionally release in environment which increases concerns about their adverse health effects on respiratory system. So it is crucial to get a deep understanding of the disease process and molecular mechanism. Epigenetic mechanisms, such as DNA methylation, have been found to play a role in the development of lung diseases by affecting expression of key genes. In addition, there could be potential different toxic effects of TiO₂ NPs between young and adult. Thus, the comparative toxicity of TiO₂ NPs in 5-week (young) and 10-week (adult) old NIH mice is investigated in this study following nasal inhalation of TiO₂ NPs at dose of 20 mg/kg (body weight)/day for 30 days. Global DNA methylation and hydroxymethylation in lung were measured. Promoter methylation of inflammatory genes (IFN-γ and TNF-α) and tissue fibrosis gene (Thy-1) were determined. Additional, RNA-sequencing runs were performed on the pulmonic libraries. We found the induced pulmonary inflammation and fibrosis were more severe in young mice. Decreased global methylation and hydroxymethylation were only found in the young group. The altered methylation in promoter of TNF-α and Thy-1 were found to play a role in the inflammatory response and fibration. RNA-sequencing showed that in pathways in cancer expression of 197 genes was up-regulated in the young mice more that in the adult mice. All these results suggested that the young ages are more sensitive to TiO₂ NP exposure and the potential of abnormal DNA methylation might be used as biomarkers of both exposure and disease development.

Whole Genome Sequencing and Progress Toward Full Inbreeding of the Mouse Collaborative Cross Population

Two key features of recombinant inbred panels are well-characterized genomes and reproducibility. Here we report on the sequenced genomes of six additional Collaborative Cross (CC) strains and on inbreeding progress of 72 CC strains. We have previously reported on the sequences of 69 CC strains that were publicly available, bringing the total of CC strains with whole genome sequence up to 75. The sequencing of these six CC strains updates the efforts toward inbreeding undertaken by the UNC Systems Genetics Core. The timing reflects our competing mandates to release to the public as many CC strains as possible while achieving an acceptable level of inbreeding. The new six strains have a higher than average founder contribution from non-domesticus strains than the previously released CC strains. Five of the six strains also have high residual heterozygosity (>14%), which may be related to non-domesticus founder contributions. Finally, we report on updated estimates on residual heterozygosity across the entire CC population using a novel, simple and cost effective genotyping platform on three mice from each strain. We observe a reduction in residual heterozygosity across all previously released CC strains. We discuss the optimal use of different genetic resources available for the CC population.

Frequency of off-targeting in genome edited pigs produced via direct injection of the CRISPR/Cas9 system into developing embryos

Background

The CRISPR/Cas9 system can effectively introduce site-specific modifications to the genome. The efficiency is high enough to induce targeted genome modifications during embryogenesis, thus increasing the efficiency of producing genetically modified animal models and having potential clinical applications as an assisted reproductive technology. Because most of the CRISPR/Cas9 systems introduce site-specific double-stranded breaks (DSBs) to induce site-specific modifications, a major concern is its potential off-targeting activity, which may hinder the application of the technology in clinics. In this study, we investigated off-targeting events in genome edited pigs/fetuses that were generated through direct injection of the CRISPR/Cas9 system into developing embryos; off-targeting activity of four different sgRNAs targeting RAG2, IL2RG, SCD5, and Ig Heavy chain were examined.

Results

First, bioinformatics analysis was applied to identify 27 potential off-targeting genes from the sgRNAs. Then, PCR amplification followed by sequencing analysis was used to verify the presence of off-targeting events. Off-targeting events were only identified from the sgRNA used to disrupt Ig Heavy chain in pigs; frequency of off-targeting was 80 and 70% on AR and RBFOX1 locus respectively. A potential PAM sequence was present in both of the off-targeting genes adjacent to probable sgRNA binding sites. Mismatches against sgRNA were present only on the 5′ side of AR, suggesting that off-targeting activities are systematic events. However, the mismatches on RBFOX1 were not limited to the 5′ side, indicating unpredictability of the events.

Conclusions

The prevalence of off-targeting is low via direct injection of CRISPR/Cas9 system into developing embryos, but the events cannot be accurately predicted. Off-targeting frequency of each CRISPR/Cas9 system should be deliberately assessed prior to its application in clinics.

Electronic supplementary material

The online version of this article (10.1186/s12896-019-0517-7) contains supplementary material, which is available to authorized users.

Structural and diffusion weighted MRI demonstrates responses to ibrutinib in a mouse model of follicular helper (Tfh) T-cell lymphoma

Recent analyses of the genetics of peripheral T-cell lymphoma (PTCL) have shown that a large proportion of cases are derived from normal follicular helper (Tfh) T-cells. The sanroque mouse strain bears a mutation that increases Tfh cell number and heterozygous animals (Roquin^san/+) develop lymphomas similar to human Tfh lymphoma. Here we demonstrate the usefulness of Roquin^san/+ animals as a pre-clinical model of Tfh lymphoma. Long latency of development and incomplete penetrance in this strain suggests the lymphomas are genetically diverse. We carried out preliminary genetic characterisation by whole exome sequencing and detected tumor specific mutations in Hsp90ab1, Ccnb3 and RhoA. Interleukin-2-inducible kinase (ITK) is expressed in Tfh lymphoma and is a potential therapeutic agent. A preclinical study of ibrutinib, a small molecule inhibitor of mouse and human ITK, in established lymphoma was carried out and showed lymphoma regression in 8/12 (67%) mice. Using T2-weighted MRI to assess lymph node volume and diffusion weighted MRI scanning as a measure of function, we showed that treatment increased mean apparent diffusion coefficient (ADC) suggesting cell death, and that change in ADC following treatment correlated with change in lymphoma volume. We suggest that heterozygous sanroque mice are a useful model of Tfh cell derived lymphomas in an immunocompetent animal.

Maternal-biased H3K27me3 correlates with paternal-specific gene expression in the human morula

In this study from Zhang et al., the authors found that H3K27me3 is strongly biased toward the maternal allele with some associated with DNA methylation–independent paternally expressed genes (PEGs) in human morulae. Their findings uncover the H3K27me3 landscape and establish a correlation between maternal-biased H3K27me3 and PEGs in human morulae.

Genomic imprinting is an epigenetic mechanism by which genes are expressed in a parental origin–dependent manner. We recently discovered that, like DNA methylation, oocyte-inherited H3K27me3 can also serve as an imprinting mark in mouse preimplantation embryos. In this study, we found H3K27me3 is strongly biased toward the maternal allele with some associated with DNA methylation–independent paternally expressed genes (PEGs) in human morulae. The H3K27me3 domains largely overlap with DNA partially methylated domains (PMDs) and occupy developmental gene promoters. Thus, our study not only reveals the H3K27me3 landscape but also establishes a correlation between maternal-biased H3K27me3 and PEGs in human morulae.

Streamlined computational pipeline for genetic background characterization of genetically engineered mice based on next generation sequencing data

BACKGROUND

Genetically engineered mice (GEM) are essential tools for understanding gene function and disease modeling. Historically, gene targeting was first done in embryonic stem cells (ESCs) derived from the 129 family of inbred strains, leading to a mixed background or congenic mice when crossed with C57BL/6 mice. Depending on the number of backcrosses and breeding strategies, genomic segments from 129-derived ESCs can be introgressed into the C57BL/6 genome, establishing a unique genetic makeup that needs characterization in order to obtain valid conclusions from experiments using GEM lines. Currently, SNP genotyping is used to detect the extent of 129-derived ESC genome introgression into C57BL/6 recipients; however, it fails to detect novel/rare variants.

RESULTS

Here, we present a computational pipeline implemented in the Galaxy platform and in BASH/R script to determine genetic introgression of GEM using next generation sequencing data (NGS), such as whole genome sequencing (WGS), whole exome sequencing (WES) and RNA-Seq. The pipeline includes strategies to uncover variants linked to a targeted locus, genome-wide variant visualization, and the identification of potential modifier genes. Although these methods apply to congenic mice, they can also be used to describe variants fixed by genetic drift. As a proof of principle, we analyzed publicly available RNA-Seq data from five congenic knockout (KO) lines and our own RNA-Seq data from the Sall2 KO line. Additionally, we performed target validation using several genetics approaches.

CONCLUSIONS

We revealed the impact of the 129-derived ESC genome introgression on gene expression, predicted potential modifier genes, and identified potential phenotypic interference in KO lines. Our results demonstrate that our new approach is an effective method to determine genetic introgression of GEM.

Genetic Models Reveal cis and trans Immune-Regulatory Activities for lincRNA-Cox2

An inducible gene expression program is a hallmark of the host inflammatory response. Recently, long intergenic non-coding RNAs (lincRNAs) have been shown to regulate the magnitude, duration, and resolution of these responses. Among these is lincRNA-Cox2, a dynamically regulated gene that broadly controls immune gene expression. To evaluate the in vivo functions of this lincRNA, we characterized multiple models of lincRNA-Cox2-deficient mice. LincRNA-Cox2-deficient macrophages and murine tissues had altered expression of inflammatory genes. Transcriptomic studies from various tissues revealed that deletion of the lincRNA-Cox2 locus also strongly impaired the basal and inducible expression of the neighboring gene prostaglandin-endoperoxide synthase (Ptgs2), encoding cyclooxygenase-2, a key enzyme in the prostaglandin biosynthesis pathway. By utilizing different genetic manipulations in vitro and in vivo, we found that lincRNA-Cox2 functions through an enhancer RNA mechanism to regulate Ptgs2. More importantly, lincRNA-Cox2 also functions in trans, independently of Ptgs2, to regulate critical innate immune genes in vivo.

Sixteen diverse laboratory mouse reference genomes define strain-specific haplotypes and novel functional loci

We report full-length draft de novo genome assemblies for 16 widely used inbred mouse strains and find extensive strain-specific haplotype variation. We identify and characterize 2,567 regions on the current mouse reference genome exhibiting the greatest sequence diversity. These regions are enriched for genes involved in pathogen defence and immunity and exhibit enrichment of transposable elements and signatures of recent retrotransposition events. Combinations of alleles and genes unique to an individual strain are commonly observed at these loci, reflecting distinct strain phenotypes. We used these genomes to improve the mouse reference genome, resulting in the completion of 10 new gene structures. Also, 62 new coding loci were added to the reference genome annotation. These genomes identified a large, previously unannotated, gene (Efcab3-like) encoding 5,874 amino acids. Mutant Efcab3-like mice display anomalies in multiple brain regions, suggesting a possible role for this gene in the regulation of brain development.

Dotting the eyes: mouse strain dependency of the lens epithelium to low dose radiation-induced DNA damage

PURPOSE

Epidemiological evidence regarding the radiosensitivity of the lens of the eye and radiation cataract development has led to changes in the EU Basic Safety Standards for protection of the lens against ionizing radiation. However, mechanistic details of lens radiation response pathways and their significance for cataractogenesis remain unclear. Radiation-induced DNA damage and the potential impairment of repair pathways within the lens epithelium, a cell monolayer that covers the anterior hemisphere of the lens, are likely to be involved.

MATERIALS AND METHODS

In this work, the lens epithelium has been analyzed for its DNA double-strand break (DSB) repair response to ionizing radiation. The responses of epithelial cells located at the anterior pole (central region) have been compared to at the very periphery of the monolayer (germinative and transitional zones). Described here are the different responses in the two regions and across four strains (C57BL/6, 129S2, BALB/c and CBA/Ca) over a low dose (0-25 mGy) in-vivo whole body X-irradiation range up to 24 hours post exposure.

RESULTS

DNA damage and repair as visualized through 53BP1 staining was present across the lens epithelium, although repair kinetics appeared non-uniform. Epithelial cells in the central region have significantly more 53BP1 foci. The sensitivities of different mouse strains have also been compared.

CONCLUSIONS

129S2 and BALB/c showed higher levels of DNA damage, with BALB/c showing significantly less inter-individual variability and appearing to be a more robust model for future DNA damage and repair studies. As a result of this study, BALB/c was identified as a suitable radiosensitive lens strain to detect and quantify early low dose ionizing radiation DNA damage effects in the mouse eye lens specifically, as an indicator of cataract formation.

A severe atherosclerosis mouse model on the resistant NOD background

Atherosclerosis is a complex disease affecting arterial blood vessels and blood flow that could result in a variety of life-threatening consequences. Disease models with diverged genomes are necessary for understanding the genetic architecture of this complex disease. Non-obese diabetic (NOD) mice are highly polymorphic and widely used for studies of type 1 diabetes and autoimmunity. Understanding atherosclerosis development in the NOD strain is of particular interest as human atherosclerosis on the diabetic and autoimmune background has not been successfully modeled. In this study, we used CRISPR/Cas9 genome editing to genetically disrupt apolipoprotein E (ApoE) and low-density lipoprotein receptor (LDLR) expression on the pure NOD background, and compared phenotype between single-gene-deleted mice and double-knockout mutants with reference to ApoE-deficient C57BL/6 mice. We found that genetic ablation of Ldlr or Apoe in NOD mice was not sufficient to establish an atherosclerosis model, in contrast to ApoE-deficient C57BL/6 mice fed a high-fat diet (HFD) for over 12 weeks. We further obtained NOD mice deficient in both LDLR and ApoE, and assessed the severity of atherosclerosis and immune response to hyperlipidemia in comparison to ApoE-deficient C57BL/6 mice. Strikingly, the double-knockout NOD mice treated with a HFD developed severe atherosclerosis with aorta narrowed by over 60% by plaques, accompanied by destruction of pancreatic islets and an inflammatory response to hyperlipidemia. Therefore, we succeeded in obtaining a genetic model with severe atherosclerosis on the NOD background, which is highly resistant to the disease. This model is useful for the study of atherosclerosis in the setting of autoimmunity.

Sequence and Structural Diversity of Mouse Y Chromosomes

Over the 180 My since their origin, the sex chromosomes of mammals have evolved a gene repertoire highly specialized for function in the male germline. The mouse Y chromosome is unique among mammalian Y chromosomes characterized to date in that it is large, gene-rich and euchromatic. Yet, little is known about its diversity in natural populations. Here, we take advantage of published whole-genome sequencing data to survey the diversity of sequence and copy number of sex-linked genes in three subspecies of house mice. Copy number of genes on the repetitive long arm of both sex chromosomes is highly variable, but sequence diversity in nonrepetitive regions is decreased relative to expectations based on autosomes. We use simulations and theory to show that this reduction in sex-linked diversity is incompatible with neutral demographic processes alone, but is consistent with recent positive selection on genes active during spermatogenesis. Our results support the hypothesis that the mouse sex chromosomes are engaged in ongoing intragenomic conflict.

Top-Down Proteomics Enables Comparative Analysis of Brain Proteoforms Between Mouse Strains

Over the past decade, advances in mass spectrometry-based proteomics have accelerated brain proteome research aimed at studying the expression, dynamic modification, interaction and function of proteins in the nervous system that are associated with physiological and behavioral processes. With the latest hardware and software improvements in top-down mass spectrometry, the technology has expanded from mere protein profiling to high-throughput identification and quantification of intact proteoforms. Murine systems are broadly used as models to study human diseases. Neuroscientists specifically study the mouse brain from inbred strains to help understand how strain-specific genotype and phenotype affect development, functioning, and disease progression. This work describes the first application of label-free quantitative top-down proteomics to the analysis of the mouse brain proteome. Operating in discovery mode, we determined physiochemical differences in brain tissue from four healthy inbred strains, C57BL/6J, DBA/2J, FVB/NJ, and BALB/cByJ, after probing their intact proteome in the 3.5-30 kDa mass range. We also disseminate these findings using a new tool for top-down proteomics, TDViewer and cataloged them in a newly established Mouse Brain Proteoform Atlas. The analysis of brain tissues from the four strains identified 131 gene products leading to the full characterization of 343 of the 593 proteoforms identified. Within the results, singly and doubly phosphorylated ARPP-21 proteoforms, known to inhibit calmodulin, were differentially expressed across the four strains. Gene ontology (GO) analysis for detected differentially expressed proteoforms also helps to illuminate the similarities and dissimilarities in phenotypes among these inbred strains.

Incorporating population-level genetic variability within laboratory models in toxicology: From the individual to the population

Humans respond to chemical exposures differently due to many factors, such as previous and concurrent stressors, age, sex, and genetic background. The vast majority of laboratory-based toxicology studies, however, have not considered the impact of population-level variability within dose-response relationships. The lack of data dealing with the influence of genetic diversity on the response to chemical exposure provides a difficult challenge for risk assessment as individuals within the population will display a wide-range of responses following toxicant challenge. Notably, the genetic background of individuals plays a major role in the variability seen in a population-level response to a drug or chemical and, thus, there is growing interest in including genetic diversity into laboratory-models. Here we outline several laboratory-based models that can be used to assay the influence of genetic variability on an individual's response to chemicals: 1) genetically-diverse cell lines, 2) human primary cells, 3) and genetically-diverse mouse panels. We also provide a succinct review for several seminal studies to highlight the capability, feasibility, and power of each of these models. This article is intended to highlight the need to include population-level genetic diversity into toxicological study designs via laboratory-based models with the goal to provide and supplement evidence in assessing the risk posed by chemicals to the human population. As such, incorporation of genetic variability will positively impact human-based risk assessment and provide empirical data to aid and influence decision-making processes in relation to chemical exposures.

Genetics of metabolic syndrome: potential clues from wild-derived inbred mouse strains

The metabolic syndrome (MetS) is a complex constellation of metabolic abnormalities including obesity, abnormal glucose metabolism, dyslipidemia, and elevated blood pressure that together substantially increase risk for cardiovascular disease and Type 2 diabetes. Both genetic and environmental factors contribute to the development of MetS, but this process is still far from understood. Human studies have revealed only part of the underlying basis. Studies in mice offer many strengths that can complement human studies to help elucidate the etiology and pathophysiology of MetS. Here we review the ways mice can contribute to MetS research. In particular, we focus on the information that can be obtained from studies of the inbred strains, with specific focus on the phenotypes of the wild-derived inbred strains. These are newly derived inbred strains that were created from wild-caught mice. They contain substantial genetic variation that is not present in the classical inbred strains, have phenotypes of relevance for MetS, and various mouse strain resources have been created to facilitate the mining of this new genetic variation. Thus studies using wild-derived inbred strains hold great promise for increasing our understanding of MetS.