Rapid identification of a disease allele in mouse through whole genome sequencing and bulk segregation analysis

C3H/HeNSlc mouse with low phospholipid transfer protein expression showed dyslipidemia

High serum levels of triglycerides (TG) and low levels of high-density lipoprotein cholesterol (HDL-C) increase the risk of coronary heart disease in humans. Herein, we first reported that the C3H/HeNSlc (C3H-S) mouse, a C3H/HeN-derived substrain, is a novel model for dyslipidemia. C3H-S showed hypertriglyceridemia and low total cholesterol (TC), HDL-C, and phospholipid (PL) concentrations. To identify the gene locus causing dyslipidemia in C3H-S, we performed genetic analysis. In F2 intercrosses between C3H-S mice and strains with normal serum lipids, the locus associated with serum lipids was identified as 163-168 Mb on chromosome 2. The phospholipid transfer protein (Pltp) gene was a candidate gene within this locus. Pltp expression and serum PLTP activity were markedly lower in C3H-S mice. Pltp expression was negatively correlated with serum TG and positively correlated with serum TC and HDL-C in F2 mice. Genome sequencing analysis revealed that an endogenous retrovirus (ERV) sequence called intracisternal A particle was inserted into intron 12 of Pltp in C3H-S. These results suggest that ERV insertion within Pltp causes aberrant splicing, leading to reduced Pltp expression in C3H-S. This study demonstrated the contribution of C3H-S to our understanding of the relationship between TG, TC, and PL metabolism via PLTP.

COL4A4 variant recently identified: lessons learned in variant interpretation-a case report

Background

Alport syndrome is a hereditary kidney disease characterized by hematuria and proteinuria. Although there have been reports of autosomal dominant COL4A4 variants, this is likely an underdiagnosed condition. Improved access to affordable genetic testing has increased the diagnosis of Alport syndrome. As genetic testing becomes ubiquitous, it is imperative that clinical nephrologists understand the benefits and challenges associated with clinical genetic testing.

Case Presentation

We present a family of Mexican descent with a heterozygous COL4A4 variant (c.5007delC, ClinVar accession numbers: SCV001580980.2, SCV001993731.1) not previously discussed in detail in the literature. The proband received a biopsy diagnosis suggestive of Fabry disease 18 years after she first developed hematuria and progressed to chronic kidney disease stage III. One year later, the proband was provisionally diagnosed with Alport syndrome after a variant of uncertain significance in the COL4A4 gene was identified following targeted family variant testing of her daughter. Upon review of the medical histories of the proband’s children and niece, all but one had the same variant. Of the four with the variant, three display clinical symptoms of hematuria, and/or proteinuria. The youngest of the four, only months old, has yet to exhibit clinical symptoms. Despite these findings there was a considerable delay in synthesizing this data, as patients were tested in different commercial genetic testing laboratories. Subsequently, understanding this family’s inheritance pattern, family history, and clinical symptoms, as well as the location of the COL4A4 variant resulted in the upgrade of the variant’s classification. Although the classification of this variant varied among different clinical genetic testing laboratories, the consensus was that this variant is likely pathogenic.

Conclusions

This COL4A4 variant (c.5007delC) not yet discussed in detail in the literature is associated with Alport syndrome. The inheritance pattern is suggestive of autosomal dominant inheritance. This report highlights the intricacies of variant interpretation and classification, the siloed nature of commercial genetic testing laboratories, and the importance of a thorough family history for proper variant interpretation. Additionally, the cases demonstrate the varied clinical presentations of Alport syndrome and suggest the utility of early screening, diagnosis, monitoring, and treatment.

Molecular Basis, Diagnostic Challenges and Therapeutic Approaches of Alport Syndrome: A Primer for Clinicians

Alport syndrome is a genetic and hereditary disease, caused by mutations in the type IV collagen genes COL4A3, COL4A4 and COL4A5, that affects the glomerular basement membrane of the kidney. It is a rare disease with an underestimated prevalence. Genetic analysis of population cohorts has revealed that it is the second most common inherited kidney disease after polycystic kidney disease. Renal involvement is the main manifestation, although it may have associated extrarenal manifestations such as hearing loss or ocular problems. The degree of expression of the disease changes according to the gene affected and other factors, known or yet to be known. The pathophysiology is not yet fully understood, although some receptors, pathways or molecules are known to be linked to the disease. There is also no specific treatment for Alport syndrome; the most commonly used are renin–angiotensin–aldosterone system inhibitors. In recent years, diagnosis has come a long way, thanks to advances in DNA sequencing technologies such as next-generation sequencing (NGS). Further research at the genetic and molecular levels in the future will complete the partial vision of the pathophysiological mechanism that we have, and will allow us to better understand what is happening and how to solve it.

Modification of an aggressive model of Alport Syndrome reveals early differences in disease pathogenesis due to genetic background

The link between mutations in collagen genes and the development of Alport Syndrome has been clearly established and a number of animal models, including knock-out mouse lines, have been developed that mirror disease observed in patients. However, it is clear from both patients and animal models that the progression of disease can vary greatly and can be modified genetically. We have identified a point mutation in Col4a4 in mice where disease is modified by strain background, providing further evidence of the genetic modification of disease symptoms. Our results indicate that C57BL/6J is a protective background and postpones end stage renal failure from 7 weeks, as seen on a C3H background, to several months. We have identified early differences in disease progression, including expression of podocyte-specific genes and podocyte morphology. In C57BL/6J mice podocyte effacement is delayed, prolonging normal renal function. The slower disease progression has allowed us to begin dissecting the pathogenesis of murine Alport Syndrome in detail. We find that there is evidence of differential gene expression during disease on the two genetic backgrounds, and that disease diverges by 4 weeks of age. We also show that an inflammatory response with increasing MCP-1 and KIM-1 levels precedes loss of renal function.

Trait mapping in diverse arthropods by bulked segregant analysis

Bulked segregant analysis (BSA) is a cross-based method for genetic mapping in sexually reproducing organisms. The method's use of bulked (pooled) samples markedly reduces the genotyping effort associated with traditional linkage mapping studies. Further, it can be applied to species with life histories or physical attributes (as for micro-insects) that render genetic mapping with other methods impractical. Recent studies in both insects and mites have revealed that advanced BSA experimental designs can resolve causal loci to narrow genomic intervals, facilitating follow-up investigations. As high-quality genomes become more widely available, BSA methods are poised to become an increasingly important tool for the rapid mapping of both monogenic and polygenic traits in diverse arthropod species.

HCFC2 is needed for IRF1- and IRF2-dependent Tlr3 transcription and for survival during viral infections

Sun et al. show that host cell factor C2 (HCFC2) is necessary for basal and induced Tlr3 transcription; deficiency of HCFC2 compromises survival during influenza virus and herpes simplex virus 1 infections in mice.

Transcriptional regulation of numerous interferon-regulated genes, including Toll-like receptor 3 (Tlr3), which encodes an innate immune sensor of viral double-stranded RNA, depends on the interferon regulatory factor 1 (IRF1) and IRF2 transcription factors. We detected specific abrogation of macrophage responses to polyinosinic-polycytidylic acid (poly(I:C)) resulting from three independent N-ethyl-N-nitrosourea–induced mutations in host cell factor C2 (Hcfc2). Hcfc2 mutations compromised survival during influenza virus and herpes simplex virus 1 infections. HCFC2 promoted the binding of IRF1 and IRF2 to the Tlr3 promoter, without which inflammatory cytokine and type I IFN responses to the double-stranded RNA analogue poly(I:C) are reduced in mouse macrophages. HCFC2 was also necessary for the transcription of a large subset of other IRF2-dependent interferon-regulated genes. Deleterious mutations of Hcfc2 may therefore increase susceptibility to diverse infectious diseases.

Innate immunity and the new forward genetics

As it is a hard-wired system for responses to microbes, innate immunity is particularly susceptible to classical genetic analysis. Mutations led the way to the discovery of many of the molecular elements of innate immune sensing and signaling pathways. In turn, the need for a faster way to find the molecular causes of mutation-induced phenotypes triggered a huge transformation in forward genetics. During the 1980s and 1990s, many heritable phenotypes were ascribed to mutations through positional cloning. In mice, this required three steps. First, a genetic mapping step was used to show that a given phenotype emanated from a circumscribed region of the genome. Second, a physical mapping step was undertaken, in which all of the region was cloned and its gene content determined. Finally, a concerted search for the mutation was performed. Such projects usually lasted for several years, but could produce breakthroughs in our understanding of biological processes. Publication of the annotated mouse genome sequence in 2002 made physical mapping unnecessary. More recently we devised a new technology for automated genetic mapping, which eliminated both genetic mapping and the search for mutations among candidate genes. The cause of phenotype can now be determined instantaneously. We have created more than 100,000 coding/splicing mutations. And by screening for defects of innate and adaptive immunity we have discovered many "new" proteins needed for innate immune function.

Genetics of enteric neuropathies

Abnormal development or disturbed functioning of the enteric nervous system (ENS), the intrinsic innervation of the gastrointestinal tract, is associated with the development of neuropathic gastrointestinal motility disorders. Here, we review the underlying molecular basis of these disorders and hypothesize that many of them have a common defective biological mechanism. Genetic burden and environmental components affecting this common mechanism are ultimately responsible for disease severity and symptom heterogeneity. We believe that they act together as the fulcrum in a seesaw balanced with harmful and protective factors, and are responsible for a continuum of symptoms ranging from neuronal hyperplasia to absence of neurons.

Variant mapping and mutation discovery in inbred mice using next-generation sequencing

Background

The development of powerful new methods for DNA sequencing enable the discovery of sequence variants, their utilization for the mapping of mutant loci, and the identification of causal variants in a single step. We have applied this approach for the analysis of ENU-mutagenized mice maintained on an inbred background.

Results

We ascertained ENU-induced variants in four different phenotypically mutant lines. These were then used as informative markers for positional cloning of the mutated genes. We tested both whole genome (WGS) and whole exome (WES) datasets.

Conclusion

Both approaches were successful as a means to localize a region of homozygosity, as well as identifying mutations of candidate genes, which could be individually assessed. As expected, the WGS strategy was more reliable, since many more ENU-induced variants were ascertained.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2173-1) contains supplementary material, which is available to authorized users.

Current strategies for mutation detection in phenotype-driven screens utilising next generation sequencing

Mutagenesis-based screens in mice are a powerful discovery platform to identify novel genes or gene functions associated with disease phenotypes. An N-ethyl-N-nitrosourea (ENU) mutagenesis screen induces single nucleotide variants randomly in the mouse genome. Subsequent phenotyping of mutant and wildtype mice enables the identification of mutated pathways resulting in phenotypes associated with a particular ENU lesion. This unbiased approach to gene discovery conducts the phenotyping with no prior knowledge of the functional mutations. Before the advent of affordable next generation sequencing (NGS), ENU variant identification was a limiting step in gene characterization, akin to ‘finding a needle in a haystack’. The emergence of a reliable reference genome alongside advances in NGS has propelled ENU mutation discovery from an arduous, time-consuming exercise to an effective and rapid form of mutation discovery. This has permitted large mouse facilities worldwide to use ENU for novel mutation discovery in a high-throughput manner, helping to accelerate basic science at the mechanistic level. Here, we describe three different strategies used to identify ENU variants from NGS data and some of the subsequent steps for mutation characterisation.

Applications of comparative evolution to human disease genetics

Direct comparison of human diseases with model phenotypes allows exploration of key areas of human biology which are often inaccessible for practical or ethical reasons. We review recent developments in comparative evolutionary approaches for finding models for genetic disease, including high-throughput generation of gene/phenotype relationship data, the linking of orthologous genes and phenotypes across species, and statistical methods for linking human diseases to model phenotypes.

Getting started in mapping-by-sequencing

Next-generation sequencing (NGS) technologies allow the cost-effective sequencing of whole genomes and have expanded the scope of genomics to novel applications, such as the genome-wide characterization of intraspecific polymorphisms and the rapid mapping and identification of point mutations. Next-generation sequencing platforms, such as the Illumina HiSeq2000 platform, are now commercially available at affordable prices and routinely produce an enormous amount of sequence data, but their wide use is often hindered by a lack of knowledge on how to manipulate and process the information produced. In this review, we focus on the strategies that are available to geneticists who wish to incorporate these novel approaches into their research but who are not familiar with the necessary bioinformatic concepts and computational tools. In particular, we comprehensively summarize case studies where the use of NGS technologies has led to the identification of point mutations, a strategy that has been dubbed "mapping-by-sequencing", and review examples from plants and other model species such as Caenorhabditis elegans, Saccharomyces cerevisiae, and Drosophila melanogaster. As these technologies are becoming cheaper and more powerful, their use is also expanding to allow mutation identification in species with larger genomes, such as many crop plants.

Real-time resolution of point mutations that cause phenovariance in mice

With the wide availability of massively parallel sequencing technologies, genetic mapping has become the rate limiting step in mammalian forward genetics. Here we introduce a method for real-time identification of N-ethyl-N-nitrosourea-induced mutations that cause phenotypes in mice. All mutations are identified by whole exome G1 progenitor sequencing and their zygosity is established in G2/G3 mice before phenotypic assessment. Quantitative and qualitative traits, including lethal effects, in single or multiple combined pedigrees are then analyzed with Linkage Analyzer, a software program that detects significant linkage between individual mutations and aberrant phenotypic scores and presents processed data as Manhattan plots. As multiple alleles of genes are acquired through mutagenesis, pooled "superpedigrees" are created to analyze the effects. Our method is distinguished from conventional forward genetic methods because it permits (1) unbiased declaration of mappable phenotypes, including those that are incompletely penetrant (2), automated identification of causative mutations concurrent with phenotypic screening, without the need to outcross mutant mice to another strain and backcross them, and (3) exclusion of genes not involved in phenotypes of interest. We validated our approach and Linkage Analyzer for the identification of 47 mutations in 45 previously known genes causative for adaptive immune phenotypes; our analysis also implicated 474 genes not previously associated with immune function. The method described here permits forward genetic analysis in mice, limited only by the rates of mutant production and screening.

Cell Receptor-Basement Membrane Interactions in Health and Disease: A Kidney-Centric View

Cell-extracellular matrix (ECM) interactions are essential for tissue development, homeostasis, and response to injury. Basement membranes (BMs) are specialized ECMs that separate epithelial or endothelial cells from stromal components and interact with cells via cellular receptors, including integrins and discoidin domain receptors. Disruption of cell-BM interactions due to either injury or genetic defects in either the ECM components or cellular receptors often lead to irreversible tissue injury and loss of organ function. Animal models that lack specific BM components or receptors either globally or in selective tissues have been used to help with our understanding of the molecular mechanisms whereby cell-BM interactions regulate organ function in physiological and pathological conditions. We review recently published works on animal models that explore how cell-BM interactions regulate kidney homeostasis in both health and disease.

A Comparison of the Chicken and Turkey Proteomes and Phosphoproteomes in the Development of Poultry-Specific Immuno-Metabolism Kinome Peptide Arrays

The use of species-specific peptide arrays for the study of animal kinomes has a proven track record of success. This technique has been used in a variety of species for the study of host–pathogen interactions and metabolism. Species-specific peptide arrays have been designed previously for use with chicken but a turkey array has never been attempted. In addition, arrays designed around individual cellular functions have been designed and utilized, but cross-function immuno-metabolic arrays have not been considered previously. Antecedent to designing separate chicken and turkey immuno-metabolic kinome peptide arrays, we show that while the chicken and turkey genomes are quite similar, the two species are much more distinct at the proteome and phosphoproteome levels. Despite a genome identity of approximately 90%, we observe that only 83% of chicken and turkey orthologous proteins display sequence matches between the two species. Further, less than 70% of kinase recognition target sequences are exact matches between chicken and turkey. Thus, our analysis shows that, at the proteome and kinome level, these two species must be considered separately in the design of novel peptide arrays. Our ultimate array design covers numerous immune and metabolic processes including innate and adaptive immunity, inflammatory responses, carbohydrate, protein, and fat metabolism, and response to hormones. We have shown the proteomic and phosphoproteomic diversity of chicken and turkey and have designed a valuable research tool for the study of immuno-metabolism within these two species.

ESCRT-II/Vps25 constrains digit number by endosome-mediated selective modulation of FGF-SHH signaling

Sorting and degradation of receptors and associated signaling molecules maintain homeostasis of conserved signaling pathways during cell specification and tissue development. Yet, whether machineries that sort signaling proteins act preferentially on different receptors and ligands in different contexts remains mysterious. Here, we show that Vacuolar protein sorting 25, Vps25, a component of ESCRT-II (Endosomal Sorting Complex Required for Transport II), directs preferential endosome-mediated modulation of FGF signaling in limbs. By ENU-induced mutagenesis, we isolated a polydactylous mouse line carrying a hypomorphic mutation of Vps25 (Vps25(ENU)). Unlike Vps25-null embryos we generated, Vps25(ENU/ENU) mutants survive until late gestation. Their limbs display FGF signaling enhancement and consequent hyperactivation of the FGF-SHH feedback loop causing polydactyly, whereas WNT and BMP signaling remain unperturbed. Notably, Vps25(ENU/ENU) Mouse Embryonic Fibroblasts exhibit aberrant FGFR trafficking and degradation; however, SHH signaling is unperturbed. These studies establish that the ESCRT-II machinery selectively limits FGF signaling in vertebrate skeletal patterning.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

The functional diversity of essential genes required for mammalian cardiac development

Genes required for an organism to develop to maturity (for which no other gene can compensate) are considered essential. The continuing functional annotation of the mouse genome has enabled the identification of many essential genes required for specific developmental processes including cardiac development. Patterns are now emerging regarding the functional nature of genes required at specific points throughout gestation. Essential genes required for development beyond cardiac progenitor cell migration and induction include a small and functionally homogenous group encoding transcription factors, ligands and receptors. Actions of core cardiogenic transcription factors from the Gata, Nkx, Mef, Hand, and Tbx families trigger a marked expansion in the functional diversity of essential genes from midgestation onwards. As the embryo grows in size and complexity, genes required to maintain a functional heartbeat and to provide muscular strength and regulate blood flow are well represented. These essential genes regulate further specialization and polarization of cell types along with proliferative, migratory, adhesive, contractile, and structural processes. The identification of patterns regarding the functional nature of essential genes across numerous developmental systems may aid prediction of further essential genes and those important to development and/or progression of disease. genesis 52:713–737, 2014.

A strategy to identify dominant point mutant modifiers of a quantitative trait

A central goal in the analysis of complex traits is to identify genes that modify a phenotype. Modifiers of a cancer phenotype may act either intrinsically or extrinsically on the salient cell lineage. Germline point mutagenesis by ethylnitrosourea can provide alleles for a gene of interest that include loss-, gain-, or alteration-of-function. Unlike strain polymorphisms, point mutations with heterozygous quantitative phenotypes are detectable in both essential and nonessential genes and are unlinked from other variants that might confound their identification and analysis. This report analyzes strategies seeking quantitative mutational modifiers of Apc^Min in the mouse. To identify a quantitative modifier of a phenotype of interest, a cluster of test progeny is needed. The cluster size can be increased as necessary for statistical significance if the founder is a male whose sperm is cryopreserved. A second critical element in this identification is a mapping panel free of polymorphic modifiers of the phenotype, to enable low-resolution mapping followed by targeted resequencing to identify the causative mutation. Here, we describe the development of a panel of six “isogenic mapping partner lines” for C57BL/6J, carrying single-nucleotide markers introduced by mutagenesis. One such derivative, B6.SNVg, shown to be phenotypically neutral in combination with Apc^Min, is an appropriate mapping partner to locate induced mutant modifiers of the Apc^Min phenotype. The evolved strategy can complement four current major initiatives in the genetic analysis of complex systems: the Genome-wide Association Study; the Collaborative Cross; the Knockout Mouse Project; and The Cancer Genome Atlas.

A mouse Col4a4 mutation causing Alport glomerulosclerosis with abnormal collagen α3α4α5(IV) trimers

A spontaneous mutation termed bilateral wasting kidneys (bwk) was identified in a colony of NONcNZO recombinant inbred mice. These mice exhibit a rapid increase of urinary albumin at an early age associated with glomerulosclerosis, interstitial nephritis, and tubular atrophy. The mutation was mapped to a location on chromosome 1 containing the Col4a3 and Col4a4 genes, for which mutations in the human orthologs cause the hereditary nephritis Alport syndrome. DNA sequencing identified a G-to-A mutation in the conserved GT splice donor of Col4a4 intron 30, resulting in skipping of exon 30 but maintaining the mRNA reading frame. Protein analyses showed that mutant collagen α3α4α5(IV) trimers were secreted and incorporated into the glomerular basement membrane (GBM), but levels were low, and GBM lesions typical of Alport syndrome were observed. Moving the mutation into the more renal damage-prone DBA/2J and 129S1/SvImJ backgrounds revealed differences in albuminuria and its rate of increase, suggesting an interaction between the Col4a4 mutation and modifier genes. This novel mouse model of Alport syndrome is the only one shown to accumulate abnormal collagen α3α4α5(IV) in the GBM, as also found in a subset of Alport patients. These mice will be valuable for testing potential therapies, for understanding abnormal collagen IV structure and assembly, and for gaining better insights into the mechanisms leading to Alport syndrome, and to the variability in the age of onset and associated phenotypes.

A nonsense mutation in TMEM95 encoding a nondescript transmembrane protein causes idiopathic male subfertility in cattle

Genetic variants underlying reduced male reproductive performance have been identified in humans and model organisms, most of them compromising semen quality. Occasionally, male fertility is severely compromised although semen analysis remains without any apparent pathological findings (i.e., idiopathic subfertility). Artificial insemination (AI) in most cattle populations requires close examination of all ejaculates before insemination. Although anomalous ejaculates are rejected, insemination success varies considerably among AI bulls. In an attempt to identify genetic causes of such variation, we undertook a genome-wide association study (GWAS). Imputed genotypes of 652,856 SNPs were available for 7962 AI bulls of the Fleckvieh (FV) population. Male reproductive ability (MRA) was assessed based on 15.3 million artificial inseminations. The GWAS uncovered a strong association signal on bovine chromosome 19 (P = 4.08×10⁻⁵⁹). Subsequent autozygosity mapping revealed a common 1386 kb segment of extended homozygosity in 40 bulls with exceptionally poor reproductive performance. Only 1.7% of 35,671 inseminations with semen samples of those bulls were successful. None of the bulls with normal reproductive performance was homozygous, indicating recessive inheritance. Exploiting whole-genome re-sequencing data of 43 animals revealed a candidate causal nonsense mutation (rs378652941, c.483C>A, p.Cys161X) in the transmembrane protein 95 encoding gene TMEM95 which was subsequently validated in 1990 AI bulls. Immunohistochemical investigations evidenced that TMEM95 is located at the surface of spermatozoa of fertile animals whereas it is absent in spermatozoa of subfertile animals. These findings imply that integrity of TMEM95 is required for an undisturbed fertilisation. Our results demonstrate that deficiency of TMEM95 severely compromises male reproductive performance in cattle and reveal for the first time a phenotypic effect associated with genomic variation in TMEM95.

Impaired male fertility is a prevalent condition in many species and is often explained by aberrant semen quality. In some cases, male fertility is severely compromised although semen quality is without any apparent pathological findings (i.e., idiopathic male subfertility). The genetic mechanisms underlying idiopathic male subfertility often remain unexplained. In the present paper, we report a recessively inherited variant of idiopathic male subfertility in a cattle population. We use 650,000 genome-wide SNP markers genotyped in >7900 artificial insemination bulls to pinpoint the underlying genomic region. We take advantage of whole-genome re-sequencing data of 43 animals to identify a causal loss-of-function mutation in TMEM95 encoding a nondescript transmembrane protein. We demonstrate that transmembrane protein 95 is located at the plasma membrane of spermatozoa of fertile animals whereas it is absent in spermatozoa of subfertile animals. Our results indicate that integrity of transmembrane protein 95 is required for an undisturbed fertilisation. This is the first report to reveal a phenotypic effect associated with genomic variation in TMEM95 in any organism.

Ion Torrent sequencing for conducting genome-wide scans for mutation mapping analysis

Mutation mapping in mice can be readily accomplished by genome wide segregation analysis of polymorphic DNA markers. In this study, we showed the efficacy of Ion Torrent next generation sequencing for conducting genome-wide scans to map and identify a mutation causing congenital heart disease in a mouse mutant, Bishu, recovered from a mouse mutagenesis screen. The Bishu mutant line generated in a C57BL/6J (B6) background was intercrossed with another inbred strain, C57BL/10J (B10), and the resulting B6/B10 hybrid offspring were intercrossed to generate mutants used for the mapping analysis. For each mutant sample, a panel of 123 B6/B10 polymorphic SNPs distributed throughout the mouse genome was PCR amplified, bar coded, and then pooled to generate a single library used for Ion Torrent sequencing. Sequencing carried out using the 314 chip yielded >600,000 usable reads. These were aligned and mapped using a custom bioinformatics pipeline. Each SNP was sequenced to a depth >500×, allowing accurate automated calling of the B6/B10 genotypes. This analysis mapped the mutation in Bishu to an interval on the proximal region of mouse chromosome 4. This was confirmed by parallel capillary sequencing of the 123 polymorphic SNPs. Further analysis of genes in the map interval identified a splicing mutation in Dnaic1(c.204+1G>A), an intermediate chain dynein, as the disease causing mutation in Bishu. Overall, our experience shows Ion Torrent amplicon sequencing is high throughput and cost effective for conducting genome-wide mapping analysis and is easily scalable for other high volume genotyping analyses.

MutMap+: genetic mapping and mutant identification without crossing in rice

Advances in genome sequencing technologies have enabled researchers and breeders to rapidly associate phenotypic variation to genome sequence differences. We recently took advantage of next-generation sequencing technology to develop MutMap, a method that allows rapid identification of causal nucleotide changes of rice mutants by whole genome resequencing of pooled DNA of mutant F2 progeny derived from crosses made between candidate mutants and the parental line. Here we describe MutMap+, a versatile extension of MutMap, that identifies causal mutations by comparing SNP frequencies of bulked DNA of mutant and wild-type progeny of M3 generation derived from selfing of an M2 heterozygous individual. Notably, MutMap+ does not necessitate artificial crossing between mutants and the wild-type parental line. This method is therefore suitable for identifying mutations that cause early development lethality, sterility, or generally hamper crossing. Furthermore, MutMap+ is potentially useful for gene isolation in crops that are recalcitrant to artificial crosses.

Massively parallel sequencing of the mouse exome to accurately identify rare, induced mutations: an immediate source for thousands of new mouse models

Accurate identification of sparse heterozygous single-nucleotide variants (SNVs) is a critical challenge for identifying the causative mutations in mouse genetic screens, human genetic diseases and cancer. When seeking to identify causal DNA variants that occur at such low rates, they are overwhelmed by false-positive calls that arise from a range of technical and biological sources. We describe a strategy using whole-exome capture, massively parallel DNA sequencing and computational analysis, which identifies with a low false-positive rate the majority of heterozygous and homozygous SNVs arising de novo with a frequency of one nucleotide substitution per megabase in progeny of N-ethyl-N-nitrosourea (ENU)-mutated C57BL/6j mice. We found that by applying a strategy of filtering raw SNV calls against known and platform-specific variants we could call true SNVs with a false-positive rate of 19.4 per cent and an estimated false-negative rate of 21.3 per cent. These error rates are small enough to enable calling a causative mutation from both homozygous and heterozygous candidate mutation lists with little or no further experimental validation. The efficacy of this approach is demonstrated by identifying the causative mutation in the Ptprc gene in a lymphocyte-deficient strain and in 11 other strains with immune disorders or obesity, without the need for meiotic mapping. Exome sequencing of first-generation mutant mice revealed hundreds of unphenotyped protein-changing mutations, 52 per cent of which are predicted to be deleterious, which now become available for breeding and experimental analysis. We show that exome sequencing data alone are sufficient to identify induced mutations. This approach transforms genetic screens in mice, establishes a general strategy for analysing rare DNA variants and opens up a large new source for experimental models of human disease.

Unlocking the bottleneck in forward genetics using whole-genome sequencing and identity by descent to isolate causative mutations

Forward genetics screens with N-ethyl-N-nitrosourea (ENU) provide a powerful way to illuminate gene function and generate mouse models of human disease; however, the identification of causative mutations remains a limiting step. Current strategies depend on conventional mapping, so the propagation of affected mice requires non-lethal screens; accurate tracking of phenotypes through pedigrees is complex and uncertain; out-crossing can introduce unexpected modifiers; and Sanger sequencing of candidate genes is inefficient. Here we show how these problems can be efficiently overcome using whole-genome sequencing (WGS) to detect the ENU mutations and then identify regions that are identical by descent (IBD) in multiple affected mice. In this strategy, we use a modification of the Lander-Green algorithm to isolate causative recessive and dominant mutations, even at low coverage, on a pure strain background. Analysis of the IBD regions also allows us to calculate the ENU mutation rate (1.54 mutations per Mb) and to model future strategies for genetic screens in mice. The introduction of this approach will accelerate the discovery of causal variants, permit broader and more informative lethal screens to be used, reduce animal costs, and herald a new era for ENU mutagenesis.

MMAPPR: mutation mapping analysis pipeline for pooled RNA-seq

Forward genetic screens in model organisms are vital for identifying novel genes essential for developmental or disease processes. One drawback of these screens is the labor-intensive and sometimes inconclusive process of mapping the causative mutation. To leverage high-throughput techniques to improve this mapping process, we have developed a Mutation Mapping Analysis Pipeline for Pooled RNA-seq (MMAPPR) that works without parental strain information or requiring a preexisting SNP map of the organism, and adapts to differential recombination frequencies across the genome. MMAPPR accommodates the considerable amount of noise in RNA-seq data sets, calculates allelic frequency by Euclidean distance followed by Loess regression analysis, identifies the region where the mutation lies, and generates a list of putative coding region mutations in the linked genomic segment. MMAPPR can exploit RNA-seq data sets from isolated tissues or whole organisms that are used for gene expression and transcriptome analysis in novel mutants. We tested MMAPPR on two known mutant lines in zebrafish, nkx2.5 and tbx1, and used it to map two novel ENU-induced cardiovascular mutants, with mutations found in the ctr9 and cds2 genes. MMAPPR can be directly applied to other model organisms, such as Drosophila and Caenorhabditis elegans, that are amenable to both forward genetic screens and pooled RNA-seq experiments. Thus, MMAPPR is a rapid, cost-efficient, and highly automated pipeline, available to perform mutant mapping in any organism with a well-assembled genome.

Exon capture and bulk segregant analysis: rapid discovery of causative mutations using high-throughput sequencing

Background

Exome sequencing has transformed human genetic analysis and may do the same for other vertebrate model systems. However, a major challenge is sifting through the large number of sequence variants to identify the causative mutation for a given phenotype. In models like Xenopus tropicalis, an incomplete and occasionally incorrect genome assembly compounds this problem. To facilitate cloning of X. tropicalis mutants identified in forward genetic screens, we sought to combine bulk segregant analysis and exome sequencing into a single step.

Results

Here we report the first use of exon capture sequencing to identify mutations in a non-mammalian, vertebrate model. We demonstrate that bulk segregant analysis coupled with exon capture sequencing is not only able to identify causative mutations but can also generate linkage information, facilitate the assembly of scaffolds, identify misassembles, and discover thousands of SNPs for fine mapping.

Conclusion

Exon capture sequencing and bulk segregant analysis is a rapid, inexpensive method to clone mutants identified in forward genetic screens. With sufficient meioses, this method can be generalized to any model system with a genome assembly, polished or unpolished, and in the latter case, it also provides many critical genomic resources.

Dissecting immunity by germline mutagenesis

The last decades have seen numerous approaches being used to decipher biological phenomena, notably the strategies we employ to defend ourselves against pathogenic attacks. From microarrays to genetics to computing technologies, all have supported a better but not yet comprehensive understanding of the pathways regulating our immune system. Limitations are notably exemplified by cases of immune deficiencies in humans that often result in high susceptibility to infections or even death, without the genetic cause being evident. To provide further insight into the mechanisms by which pathogen detection and eradication occur, several in vivo strategies can be used. The current review focuses on one of them, namely germline mutagenesis in the mouse. After describing the main technical aspects of this forward genetic approach, we will discuss particular germline mutants that have all been instrumental in deciphering innate or adaptive immune responses. Mutations in previously uncharacterized genes in the mouse, like Unc93B or Themis, have demonstrated the impartiality of forward genetics and led to the identification of new crucial immunity actors. Some mutants, like PanR1, have informed us on particular protein domains and their specific functions. Finally, certain mutations identified by this non-hypothesis-driven method have revealed previously unknown gene functions, as recently illustrated by memi, which links a particular nucleoside salvage enzyme to cell proliferation and apoptosis.

ENU-induced phenovariance in mice: inferences from 587 mutations

Background

We present a compendium of N-ethyl-N-nitrosourea (ENU)-induced mouse mutations, identified in our laboratory over a period of 10 years either on the basis of phenotype or whole genome and/or whole exome sequencing, and archived in the Mutagenetix database. Our purpose is threefold: 1) to formally describe many point mutations, including those that were not previously disclosed in peer-reviewed publications; 2) to assess the characteristics of these mutations; and 3) to estimate the likelihood that a missense mutation induced by ENU will create a detectable phenotype.

Findings

In the context of an ENU mutagenesis program for C57BL/6J mice, a total of 185 phenotypes were tracked to mutations in 129 genes. In addition, 402 incidental mutations were identified and predicted to affect 390 genes. As previously reported, ENU shows strand asymmetry in its induction of mutations, particularly favoring T to A rather than A to T in the sense strand of coding regions and splice junctions. Some amino acid substitutions are far more likely to be damaging than others, and some are far more likely to be observed. Indeed, from among a total of 494 non-synonymous coding mutations, ENU was observed to create only 114 of the 182 possible amino acid substitutions that single base changes can achieve. Based on differences in overt null allele frequencies observed in phenotypic vs. non-phenotypic mutation sets, we infer that ENU-induced missense mutations create detectable phenotype only about 1 in 4.7 times. While the remaining mutations may not be functionally neutral, they are, on average, beneath the limits of detection of the phenotypic assays we applied.

Conclusions

Collectively, these mutations add to our understanding of the chemical specificity of ENU, the types of amino acid substitutions it creates, and its efficiency in causing phenovariance. Our data support the validity of computational algorithms for the prediction of damage caused by amino acid substitutions, and may lead to refined predictions as to whether specific amino acid changes are responsible for observed phenotypes. These data form the basis for closer in silico estimations of the number of genes mutated to a state of phenovariance by ENU within a population of G3 mice.

The Mammalian Phenotype Ontology as a unifying standard for experimental and high-throughput phenotyping data

The Mammalian Phenotype Ontology (MP) is a structured vocabulary for describing mammalian phenotypes and serves as a critical tool for efficient annotation and comprehensive retrieval of phenotype data. Importantly, the ontology contains broad and specific terms, facilitating annotation of data from initial observations or screens and detailed data from subsequent experimental research. Using the ontology structure, data are retrieved inclusively, i.e., data annotated to chosen terms and to terms subordinate in the hierarchy. Thus, searching for "abnormal craniofacial morphology" also returns annotations to "megacephaly" and "microcephaly," more specific terms in the hierarchy path. The development and refinement of the MP is ongoing, with new terms and modifications to its organization undergoing continuous assessment as users and expert reviewers propose expansions and revisions. A wealth of phenotype data on mouse mutations and variants annotated to the MP already exists in the Mouse Genome Informatics database. These data, along with data curated to the MP by many mouse mutagenesis programs and mouse repositories, provide a platform for comparative analyses and correlative discoveries. The MP provides a standard underpinning to mouse phenotype descriptions for existing and future experimental and large-scale phenotyping projects. In this review we describe the MP as it presently exists, its application to phenotype annotations, the relationship of the MP to other ontologies, and the integration of the MP within large-scale phenotyping projects. Finally we discuss future application of the MP in providing standard descriptors of the phenotype pipeline test results from the International Mouse Phenotype Consortium projects.

High throughput sequencing approaches to mutation discovery in the mouse

Phenotype-driven approaches in mice are powerful strategies for the discovery of genes and gene functions and for unravelling complex biological mechanisms. Traditional methods for mutation discovery are reliable and robust, but they can also be laborious and time consuming. Recently, high-throughput sequencing (HTS) technologies have revolutionised the process of forward genetics in mice by paving the way to rapid mutation discovery. However, successful application of HTS for mutation discovery relies heavily on the sequencing approach employed and strategies for data analysis. Here we review current HTS applications and resources for mutation discovery and provide an overview of the practical considerations for HTS implementation and data analysis.

A forward genetic screen reveals roles for Nfkbid, Zeb1, and Ruvbl2 in humoral immunity

Using chemical germ-line mutagenesis, we screened mice for defects in the humoral immune response to a type II T-independent immunogen and an experimental alphavirus vector. A total of 26 mutations that impair humoral immunity were recovered, and 19 of these mutations have been positionally cloned. Among the phenovariants were bumble, cellophane, and Worker ascribed to mutations in Nfkbid, Zeb1, and Ruvbl2, respectively. We show that IκBNS, the nuclear IκB-like protein encoded by Nfkbid, is required for the development of marginal zone and peritoneal B-1 B cells and additionally required for extrafollicular antibody responses to T-independent and -dependent immunogens. Zeb1 is also required for marginal zone and peritoneal B-1 B-cell development as well as T-cell development, germinal center formation, and memory B-cell responses. Finally, Ruvbl2 is required for T-cell development and maximal T-dependent antibody responses. Collectively, the mutations that we identified give us insight into the points at which disruption of an antibody response can occur. All of the mutations identified to date directly affect lymphocyte development or function; none have an exclusive effect on cells of the innate immune system.

Discovery Genetics - The History and Future of Spontaneous Mutation Research

Historically, spontaneous mutations in mice have served as valuable models of heritable human diseases, contributing substantially to our understanding of both disease mechanisms and basic biological pathways. While advances in molecular technologies have improved our ability to create mouse models of human disease through targeted mutagenesis and transgenesis, spontaneous mutations continue to provide valuable research tools for discovery of novel genes and functions. In addition, the genetic defects caused by spontaneous mutations are molecularly similar to mutations in the human genome and, therefore often produce phenotypes that more closely resemble those characteristic of human disease than do genetically engineered mutations. Due to the rarity with which spontaneous mutations arise and the animal intensive nature of their genetic analysis, large-scale spontaneous mutation analysis has traditionally been limited to large mammalian genetics institutes. More recently, ENU mutagenesis and new screening methods have increased the rate of mutant strain discovery, and high-throughput DNA sequencing has enabled rapid identification of the underlying genes and their causative mutations. Here, we discuss the continued value of spontaneous mutations for biomedical research.

Mutation mapping and identification by whole-genome sequencing

Genetic mapping of mutations in model systems has facilitated the identification of genes contributing to fundamental biological processes including human diseases. However, this approach has historically required the prior characterization of informative markers. Here we report a fast and cost-effective method for genetic mapping using next-generation sequencing that combines single nucleotide polymorphism discovery, mutation localization, and potential identification of causal sequence variants. In contrast to prior approaches, we have developed a hidden Markov model to narrowly define the mutation area by inferring recombination breakpoints of chromosomes in the mutant pool. In addition, we created an interactive online software resource to facilitate automated analysis of sequencing data and demonstrate its utility in the zebrafish and mouse models. Our novel methodology and online tools will make next-generation sequencing an easily applicable resource for mutation mapping in all model systems.

Modifier genes and the plasticity of genetic networks in mice

Modifier genes are an integral part of the genetic landscape in both humans and experimental organisms, but have been less well explored in mammals than other systems. A growing number of modifier genes in mouse models of disease nonetheless illustrate the potential for novel findings, while new technical advances promise many more to come. Modifier genes in mouse models include induced mutations and spontaneous or wild-derived variations captured in inbred strains. Identification of modifiers among wild-derived variants in particular should detect disease modifiers that have been shaped by selection and might therefore be compatible with high fitness and function. Here we review selected examples and argue that modifier genes derived from natural variation may provide a bias for nodes in genetic networks that have greater intrinsic plasticity and whose therapeutic manipulation may therefore be more resilient to side effects than conventional targets.

Fast homozygosity mapping and identification of a zebrafish ENU-induced mutation by whole-genome sequencing

Forward genetics using zebrafish is a powerful tool for studying vertebrate development through large-scale mutagenesis. Nonetheless, the identification of the molecular lesion is still laborious and involves time-consuming genetic mapping. Here, we show that high-throughput sequencing of the whole zebrafish genome can directly locate the interval carrying the causative mutation and at the same time pinpoint the molecular lesion. The feasibility of this approach was validated by sequencing the m1045 mutant line that displays a severe hypoplasia of the exocrine pancreas. We generated 13 Gb of sequence, equivalent to an eightfold genomic coverage, from a pool of 50 mutant embryos obtained from a map-cross between the AB mutant carrier and the WIK polymorphic strain. The chromosomal region carrying the causal mutation was localized based on its unique property to display high levels of homozygosity among sequence reads as it derives exclusively from the initial AB mutated allele. We developed an algorithm identifying such a region by calculating a homozygosity score along all chromosomes. This highlighted an 8-Mb window on chromosome 5 with a score close to 1 in the m1045 mutants. The sequence analysis of all genes within this interval revealed a nonsense mutation in the snapc4 gene. Knockdown experiments confirmed the assertion that snapc4 is the gene whose mutation leads to exocrine pancreas hypoplasia. In conclusion, this study constitutes a proof-of-concept that whole-genome sequencing is a fast and effective alternative to the classical positional cloning strategies in zebrafish.

Genetic Modifier Screens in Mice

ENU mutagenesis is a forward genetics strategy in which random mutagenesis and phenotypic screening is used to identify genes based on the phenotype induced when they are mutated. A modifier screen is a type of screen in which mice with a pre-existing phenotype are utilized to identify mutations that can enhance or suppress this phenotype. This approach has the potential to uncover missing pathway members, reveal novel genetic interactions, and pinpoint new drug targets. Considerations when planning a suppressor screen include current knowledge, genomic footprint, penetrance, variance, robustness, latency of the starting phenotype, viability, fertility, genetic background and ENU tolerance of starting strain, screening assay, mouse numbers required, and mutation identification strategy. Practical advice on each of these is provided in this review. Curr. Protoc. Mouse Biol. 2:75-87 © 2012 by John Wiley & Sons, Inc.

Multiplex Chromosomal Exome Sequencing Accelerates Identification of ENU-Induced Mutations in the Mouse

Forward genetic screens in Mus musculus have proved powerfully informative by revealing unsuspected mechanisms governing basic biological processes. This approach uses potent chemical mutagens, such as N-ethyl-N-nitrosourea (ENU), to randomly induce mutations in mice, which are then bred and phenotypically screened to identify lines that disrupt a specific biological process of interest. Although identifying a mutation using the rich resources of mouse genetics is straightforward, it is unfortunately neither fast nor cheap. Here we show that detecting newly induced causal variants in a forward genetic screen can be accelerated dramatically using a methodology that combines multiplex chromosome-specific exome capture, next-generation sequencing, rapid mapping, sequence annotation, and variation filtering. The key innovation of our method is multiplex capture and sequence that allows the simultaneous survey of both mutant, parental, and background strains in a single experiment. By comparing variants identified in mutant offspring with those found in dbSNP, the unmutagenized background strains, and parental lines, induced causative mutations can be distinguished immediately from preexisting variation or experimental artifact. Here we demonstrate this approach to find the causative mutations induced in four novel ENU lines identified from a recent ENU screen. In all four cases, after applying our method, we found six or fewer putative mutations (and sometimes only a single one). Determining the causative variant was then easily achieved through standard segregation approaches. We have developed this process into a community resource that will speed up individual labs’ ability to identify the genetic lesion in mutant mouse lines; all of our reagents and software tools are open source and available to the broader scientific community.

Efficient mapping and cloning of mutations in zebrafish by low-coverage whole-genome sequencing

The generation and analysis of mutants in zebrafish has been instrumental in defining the genetic regulation of vertebrate development, physiology, and disease. However, identifying the genetic changes that underlie mutant phenotypes remains a significant bottleneck in the analysis of mutants. Whole-genome sequencing has recently emerged as a fast and efficient approach for identifying mutations in nonvertebrate model organisms. However, this approach has not been applied to zebrafish due to the complicating factors of having a large genome and lack of fully inbred lines. Here we provide a method for efficiently mapping and detecting mutations in zebrafish using these new parallel sequencing technologies. This method utilizes an extensive reference SNP database to define regions of homozygosity-by-descent by low coverage, whole-genome sequencing of pooled DNA from only a limited number of mutant F(2) fish. With this approach we mapped each of the five different zebrafish mutants we sequenced and identified likely causative nonsense mutations in two and candidate mutations in the remainder. Furthermore, we provide evidence that one of the identified mutations, a nonsense mutation in bmp1a, underlies the welded mutant phenotype.

Increased susceptibility to DNA virus infection in mice with a GCN2 mutation

The downregulation of translation through eIF2α phosphorylation is a cellular response to diverse stresses, including viral infection, and is mediated by the GCN2 kinase, protein kinase R (PKR), protein kinase-like endoplasmic reticulum kinase (PERK), and heme-regulated inhibitor kinase (HRI). Although PKR plays a major role in defense against viruses, other eIF2α kinases also may respond to viral infection and contribute to the shutdown of protein synthesis. Here we describe the recessive, loss-of-function mutation atchoum (atc) in Eif2ak4, encoding GCN2, which increased susceptibility to infection by the double-stranded DNA viruses mouse cytomegalovirus (MCMV) and human adenovirus. This mutation was identified by screening macrophages isolated from mice carrying N-ethyl-N-nitrosourea (ENU)-induced mutations. Cells from Eif2ak4(atc/atc) mice failed to phosphorylate eIF2α in response to MCMV. Importantly, homozygous Eif2ak4(atc) mice showed a modest increase in susceptibility to MCMV infection, demonstrating that translational arrest dependent on GCN2 contributes to the antiviral response in vivo.

Mutation of the gastric hydrogen-potassium ATPase alpha subunit causes iron-deficiency anemia in mice

Iron is an essential component of heme and hemoglobin, and therefore restriction of iron availability directly limits erythropoiesis. In the present study, we report a defect in iron absorption that results in iron-deficiency anemia, as revealed by an N-ethyl-N-nitrosourea-induced mouse phenotype called sublytic. Homozygous sublytic mice develop hypochromic microcytic anemia with reduced osmotic fragility of RBCs. The sublytic phenotype stems from impaired gastrointestinal iron absorption caused by a point mutation of the gastric hydrogen-potassium ATPase α subunit encoded by Atp4a, which results in achlorhydria. The anemia of sublytic homozygotes can be corrected by feeding with a high-iron diet or by parenteral injection of iron dextran; rescue can also be achieved by providing acidified drinking water to sublytic homozygotes. These findings establish the necessity of the gastric proton pump for iron absorption and effective erythropoiesis.

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.

Second-generation high-throughput forward genetic screen in mice to isolate subtle behavioral mutants

Forward genetic screens have been highly successful in revealing roles of genes and pathways in complex biological events. Traditionally these screens have focused on isolating mutants with the greatest phenotypic deviance, with the hopes of discovering genes that are central to the biological event being investigated. Behavioral screens in mice typically use simple activity-based assays as endophenotypes for more complex emotional states of the animal. They generally set the selection threshold for a putative mutant at 3 SDs (z score of 3) from the average behavior of normal animals to minimize false-positive results. Behavioral screens using a high threshold for detection have generally had limited success, with high false-positive rates and subtle phenotypic differences that have made mapping and cloning difficult. In addition, targeted reverse genetic approaches have shown that when genes central to behaviors such as open field behavior, psychostimulant response, and learning and memory tasks are mutated, they produce subtle phenotypes that differ from wild-type animals by 1 to 2 SDs (z scores of 1 to 2). We have conducted a second-generation (G2) dominant N-ethyl-N-nitrosourea (ENU) screen especially designed to detect subtle behavioral mutants for open field activity and psychostimulant response behaviors. We successfully detect mutant lines with only 1 to 2 SD shifts in mean response compared with wild-type control animals and present a robust statistical and methodological framework for conducting such forward genetic screens. Using this methodology we have screened 229 ENU mutant lines and have identified 15 heritable mutant lines. We conclude that for screens in mice that use activity-based endophenotypic measurements for complex behavioral states, this G2 screening approach yields better results.