Mutation discovery in mice by whole exome sequencing

CAGI, the Critical Assessment of Genome Interpretation, establishes progress and prospects for computational genetic variant interpretation methods

BACKGROUND

The Critical Assessment of Genome Interpretation (CAGI) aims to advance the state-of-the-art for computational prediction of genetic variant impact, particularly where relevant to disease. The five complete editions of the CAGI community experiment comprised 50 challenges, in which participants made blind predictions of phenotypes from genetic data, and these were evaluated by independent assessors.

RESULTS

Performance was particularly strong for clinical pathogenic variants, including some difficult-to-diagnose cases, and extends to interpretation of cancer-related variants. Missense variant interpretation methods were able to estimate biochemical effects with increasing accuracy. Assessment of methods for regulatory variants and complex trait disease risk was less definitive and indicates performance potentially suitable for auxiliary use in the clinic.

CONCLUSIONS

Results show that while current methods are imperfect, they have major utility for research and clinical applications. Emerging methods and increasingly large, robust datasets for training and assessment promise further progress ahead.

A mouse model of cone photoreceptor function loss (cpfl9) with degeneration due to a mutation in Gucy2e

During routine screening of mouse strains and stocks by the Eye Mutant Resource at The Jackson Laboratory for genetic mouse models of human ocular disorders, we identified cpfl9, a mouse model with cone photoreceptor function loss. The mice exhibited an early-onset phenotype that was easily recognized by the absence of a cone-mediated b-wave electroretinography response and by a reduction in rod-mediated photoresponses at four weeks of age. By genetic mapping and high-throughput sequencing of a whole exome capture library of cpfl9, a homozygous 25 bp deletion within exon 11 of the Gucy2e gene was identified, which is predicted to result in a frame shift leading to premature termination. The corresponding protein in human, retinal guanylate cyclase 1 (GUCY2D), plays an important role in rod and cone photoreceptor cell function. Loss-of-function mutations in human GUCY2D cause LCA1, one of the most common forms of Leber congenital amaurosis, which results in blindness at birth or in early childhood. The early loss of cone and reduced rod photoreceptor cell function in the cpfl9 mutant is accompanied by a later, progressive loss of cone and rod photoreceptor cells, which may be relevant to understanding disease pathology in a subset of LCA1 patients and in individuals with cone-rod dystrophy caused by recessive GUCY2D variants. cpfl9 mice will be useful for studying the role of Gucy2e in the retina.

A Dpagt1 Missense Variant Causes Degenerative Retinopathy without Myasthenic Syndrome in Mice

Congenital disorders of glycosylation (CDG) are a heterogenous group of primarily autosomal recessive mendelian diseases caused by disruptions in the synthesis of lipid-linked oligosaccharides and their transfer to proteins. CDGs usually affect multiple organ systems and vary in presentation, even within families. There is currently no cure, and treatment is aimed at ameliorating symptoms and improving quality of life. Here, we describe a chemically induced mouse mutant, tvrm76, with early-onset photoreceptor degeneration. The recessive mutation was mapped to Chromosome 9 and associated with a missense mutation in the Dpagt1 gene encoding UDP-N-acetyl-D-glucosamine:dolichyl-phosphate N-acetyl-D-glucosaminephosphotransferase (EC 2.7.8.15). The mutation is predicted to cause a substitution of aspartic acid with glycine at residue 166 of DPAGT1. This represents the first viable animal model of a Dpagt1 mutation and a novel phenotype for a CDG. The increased expression of Ddit3, and elevated levels of HSPA5 (BiP) suggest the presence of early-onset endoplasmic reticulum (ER) stress. These changes were associated with the induction of photoreceptor apoptosis in tvrm76 retinas. Mutations in human DPAGT1 cause myasthenic syndrome-13 and severe forms of a congenital disorder of glycosylation Type Ij. In contrast, Dpagt1tvrm76 homozygous mice present with congenital photoreceptor degeneration without overt muscle or muscular junction involvement. Our results suggest the possibility of DPAGT1 mutations in human patients that present primarily with retinitis pigmentosa, with little or no muscle disease. Variants in DPAGT1 should be considered when evaluating cases of non-syndromic retinal degeneration.

Identification of Arhgef12 and Prkci as genetic modifiers of retinal dysplasia in the Crb1rd8 mouse model

Mutations in the apicobasal polarity gene CRB1 lead to diverse retinal diseases, such as Leber congenital amaurosis, cone-rod dystrophy, retinitis pigmentosa (with and without Coats-like vasculopathy), foveal retinoschisis, macular dystrophy, and pigmented paravenous chorioretinal atrophy. Limited correlation between disease phenotypes and CRB1 alleles, and evidence that patients sharing the same alleles often present with different disease features, suggest that genetic modifiers contribute to clinical variation. Similarly, the retinal phenotype of mice bearing the Crb1 retinal degeneration 8 (rd8) allele varies with genetic background. Here, we initiated a sensitized chemical mutagenesis screen in B6.Cg-Crb1rd8/Pjn, a strain with a mild clinical presentation, to identify genetic modifiers that cause a more severe disease phenotype. Two models from this screen, Tvrm266 and Tvrm323, exhibited increased retinal dysplasia. Genetic mapping with high-throughput exome and candidate-gene sequencing identified causative mutations in Arhgef12 and Prkci, respectively. Epistasis analysis of both strains indicated that the increased dysplastic phenotype required homozygosity of the Crb1rd8 allele. Retinal dysplastic lesions in Tvrm266 mice were smaller and caused less photoreceptor degeneration than those in Tvrm323 mice, which developed an early, large diffuse lesion phenotype. At one month of age, Müller glia and microglia mislocalization at dysplastic lesions in both modifier strains was similar to that in B6.Cg-Crb1rd8/Pjn mice but photoreceptor cell mislocalization was more extensive. External limiting membrane disruption was comparable in Tvrm266 and B6.Cg-Crb1rd8/Pjn mice but milder in Tvrm323 mice. Immunohistological analysis of mice at postnatal day 0 indicated a normal distribution of mitotic cells in Tvrm266 and Tvrm323 mice, suggesting normal early development. Aberrant electroretinography responses were observed in both models but functional decline was significant only in Tvrm323 mice. These results identify Arhgef12 and Prkci as modifier genes that differentially shape Crb1-associated retinal disease, which may be relevant to understanding clinical variability and underlying disease mechanisms in humans.

Error-prone, stress-induced 3' flap-based Okazaki fragment maturation supports cell survival

How cells with DNA replication defects acquire mutations that allow them to escape apoptosis under environmental stress is a long-standing question. Here, we report that an error-prone Okazaki fragment maturation (OFM) pathway is activated at restrictive temperatures in rad27Δ yeast cells. Restrictive temperature stress activated Dun1, facilitating transformation of unprocessed 5′ flaps into 3′ flaps, which were removed by 3′ nucleases, including DNA polymerase δ (Polδ). However, at certain regions, 3′ flaps formed secondary structures that facilitated 3′ end extension rather than degradation, producing alternative duplications with short spacer sequences, such as pol3 internal tandem duplications. Consequently, little 5′ flap was formed, suppressing rad27Δ-induced lethality at restrictive temperatures. We define a stress-induced, error-prone OFM pathway that generates mutations that counteract replication defects and drive cellular evolution and survival.

Somatic Genetic Mosaicism in the Apolipoprotein E-null Mouse Aorta

In addition to genetic and epigenetic inheritance, somatic variation may contribute to cardiovascular disease (CVD) risk. CVD-associated somatic mutations have been reported in human clonal hematopoiesis, but evidence in the atheroma is lacking. To probe for somatic variation in atherosclerosis, we sought single-nucleotide private variants (PVs) in whole-exome sequencing (WES) data of aorta, liver, and skeletal muscle of two C57BL/6J coisogenic male ApoE null/wild-type (WT) sibling pairs, and RNA-seq data of one of the two pairs. Relative to the C57BL/6 reference genome, we identified 9 and 11 ApoE null aorta- and liver-specific PVs that were shared by all WES and RNA-seq datasets. Corresponding PVs in WT sibling aorta and liver were 1 and 0, respectively, and not overlapping with ApoE null PVs. Pyrosequencing analysis of 4 representative PVs in 17 ApoE null aortas and livers confirmed tissue-specific shifts toward the alternative allele, in addition to significant deviations from mendelian allele ratios. Notably, all aorta and liver PVs were present in the dbSNP database and were predominantly transition mutations within atherosclerosis-related genes. The majority of PVs were in discrete clusters approximately 3 Mb and 65 to 73 Mb away from hypermutable immunoglobin loci in chromosome 6. These features were largely shared with previously reported CVD-associated somatic mutations in human clonal hematopoiesis. The observation that SNPs exhibit tissue-specific somatic DNA mosaicism in ApoE null mice is potentially relevant for genetic association study design. The proximity of PVs to hypermutable loci suggests testable mechanistic hypotheses.

A missense mutation in Pitx2 leads to early-onset glaucoma via NRF2-YAP1 axis

Glaucoma is a leading cause of blindness, affecting 70 million people worldwide. Owing to the similarity in anatomy and physiology between human and mouse eyes and the ability to genetically manipulate mice, mouse models are an invaluable resource for studying mechanisms underlying disease phenotypes and for developing therapeutic strategies. Here, we report the discovery of a new mouse model of early-onset glaucoma that bears a transversion substitution c. G344T, which results in a missense mutation, p. R115L in PITX2. The mutation causes an elevation in intraocular pressure (IOP) and progressive death of retinal ganglion cells (RGC). These ocular phenotypes recapitulate features of pathologies observed in human glaucoma. Increased oxidative stress was evident in the inner retina. We demonstrate that the mutant PITX2 protein was not capable of binding to Nuclear factor-like 2 (NRF2), which regulates Pitx2 expression and nuclear localization, and to YAP1, which is necessary for co-initiation of transcription of downstream targets. PITX2-mediated transcription of several antioxidant genes were also impaired. Treatment with N-Acetyl-L-cysteine exerted a profound neuroprotective effect on glaucoma-associated neuropathies, presumably through inhibition of oxidative stress. Our study demonstrates that a disruption of PITX2 leads to glaucoma optic pathogenesis and provides a novel early-onset glaucoma model that will enable elucidation of mechanisms underlying the disease as well as to serve as a resource to test new therapeutic strategies.

Evolutionary and functional lessons from human-specific amino acid substitution matrices

Human genetic variation in coding regions is fundamental to the study of protein structure and function. Most methods for interpreting missense variants consider substitution measures derived from homologous proteins across different species. In this study, we introduce human-specific amino acid (AA) substitution matrices that are based on genetic variations in the modern human population. We analyzed the frequencies of >4.8M single nucleotide variants (SNVs) at codon and AA resolution and compiled human-centric substitution matrices that are fundamentally different from classic cross-species matrices (e.g. BLOSUM, PAM). Our matrices are asymmetric, with some AA replacements showing significant directional preference. Moreover, these AA matrices are only partly predicted by nucleotide substitution rates. We further test the utility of our matrices in exposing functional signals of experimentally-validated protein annotations. A significant reduction in AA transition frequencies was observed across nine post-translational modification (PTM) types and four ion-binding sites. Our results propose a purifying selection signal in the human proteome across a diverse set of functional protein annotations and provide an empirical baseline for interpreting human genetic variation in coding regions.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

Museum genomics reveals the rapid decline and extinction of Australian rodents since European settlement

Australia has the highest historically recorded rate of mammalian extinction in the world, with 34 terrestrial species declared extinct since European colonization in 1788. Among Australian mammals, rodents have been the most severely affected by these recent extinctions; however, given a sparse historical record, the scale and timing of their decline remain unresolved. Using museum specimens up to 184 y old, we generate genomic-scale data from across the entire assemblage of Australian hydromyine rodents (i.e., eight extinct species and their 42 living relatives). We reconstruct a phylogenomic tree for these species spanning ∼5.2 million years, revealing a cumulative total of 10 million years (>10%) of unique evolutionary history lost to extinction within the past ∼150 y. We find no evidence for reduced genetic diversity in extinct species just prior to or during decline, indicating that their extinction was extremely rapid. This suggests that populations of extinct Australian rodents were large prior to European colonization, and that genetic diversity does not necessarily protect species from catastrophic extinction. In addition, comparative analyses suggest that body size and biome interact to predict extinction and decline, with larger species more likely to go extinct. Finally, we taxonomically resurrect a species from extinction, Gould's mouse (Pseudomys gouldii Waterhouse, 1839), which survives as an island population in Shark Bay, Western Australia (currently classified as Pseudomys fieldi Waite, 1896). With unprecedented sampling across a radiation of extinct and living species, we unlock a previously inaccessible historical perspective on extinction in Australia. Our results highlight the capacity of collections-based research to inform conservation and management of persisting species.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Rat PRDM9 shapes recombination landscapes, duration of meiosis, gametogenesis, and age of fertility

BACKGROUND

Vertebrate meiotic recombination events are concentrated in regions (hotspots) that display open chromatin marks, such as trimethylation of lysines 4 and 36 of histone 3 (H3K4me3 and H3K36me3). Mouse and human PRDM9 proteins catalyze H3K4me3 and H3K36me3 and determine hotspot positions, whereas other vertebrates lacking PRDM9 recombine in regions with chromatin already opened for another function, such as gene promoters. While these other vertebrate species lacking PRDM9 remain fertile, inactivation of the mouse Prdm9 gene, which shifts the hotspots to the functional regions (including promoters), typically causes gross fertility reduction; and the reasons for these species differences are not clear.

RESULTS

We introduced Prdm9 deletions into the Rattus norvegicus genome and generated the first rat genome-wide maps of recombination-initiating double-strand break hotspots. Rat strains carrying the same wild-type Prdm9 allele shared 88% hotspots but strains with different Prdm9 alleles only 3%. After Prdm9 deletion, rat hotspots relocated to functional regions, about 40% to positions corresponding to Prdm9-independent mouse hotspots, including promoters. Despite the hotspot relocation and decreased fertility, Prdm9-deficient rats of the SHR/OlaIpcv strain produced healthy offspring. The percentage of normal pachytene spermatocytes in SHR-Prdm9 mutants was almost double than in the PWD male mouse oligospermic sterile mutants. We previously found a correlation between the crossover rate and sperm presence in mouse Prdm9 mutants. The crossover rate of SHR is more similar to sperm-carrying mutant mice, but it did not fully explain the fertility of the SHR mutants. Besides mild meiotic arrests at rat tubular stages IV (mid-pachytene) and XIV (metaphase), we also detected postmeiotic apoptosis of round spermatids. We found delayed meiosis and age-dependent fertility in both sexes of the SHR mutants.

CONCLUSIONS

We hypothesize that the relative increased fertility of rat versus mouse Prdm9 mutants could be ascribed to extended duration of meiotic prophase I. While rat PRDM9 shapes meiotic recombination landscapes, it is unnecessary for recombination. We suggest that PRDM9 has additional roles in spermatogenesis and speciation-spermatid development and reproductive age-that may help to explain male-specific hybrid sterility.

Precise detection of a murine germline mutation of the Notch3 gene associated with kyphosis and developmental disorders

Objective:

Humpback (hpbk) mice harbor a pathogenic mutation in the Notch3 gene and can serve as a beneficial animal model for investigating human myopathy, kyphosis, and developmental disorders, including lateral meningocele syndrome. Detection of the point mutation in hpbk mice is important for maintaining strains and scrutinizing genetic rescues, especially considering that homozygous mice are infertile and indistinguishable from their littermates at a young age. This study aimed for the development of a novel, precise, and time-saving genotyping method to identify the mutation in hpbk mice.

Materials and Methods:

In order to study the hpbk mouse line, we describe how we applied several tools, including quantitative polymerase chain reaction (qPCR), multiplex tetra-primer amplification-refractory mutation system (ARMS-PCR) and Sanger sequencing, toward the recognition of heterozygous and homozygous mice.

Results:

The Notch3 mutation was clearly identified using qPCR and ARMS assays, but the latter was a more precise and cost-effective approach. The lengths of the ARMS-PCR amplicons are 210 bp and 164 bp for the wild-type and hpbk alleles, respectively. Moreover, the genotyping results for each mouse were corroborated by Sanger DNA sequencing.

Conclusion:

Our newly developed PCR-based ARMS system affords a swift and precise way to genotype the hpbk mice. ARMS-PCR does not rely on any advanced equipment and is useful as a genotyping method for other model organisms that harbor a pathogenic variant.

Extracellular Matrix Remodeling During Palate Development

The morphogenesis of the mammalian secondary plate is a series of highly dynamic developmental process, including the palate shelves vertical outgrowth, elevation to the horizontal plane and complete fusion in the midline. Extracellular matrix (ECM) proteins not only form the basic infrastructure for palatal mesenchymal cells to adhere via integrins but also interact with cells to regulate their functions such as proliferation and differentiation. ECM remodeling is essential for palatal outgrowth, expansion, elevation, and fusion. Multiple signaling pathways important for palatogenesis such as FGF, TGF β, BMP, and SHH remodels ECM dynamics. Dysregulation of ECM such as HA synthesis or ECM breakdown enzymes MMPs or ADAMTS causes cleft palate in mouse models. A better understanding of ECM remodeling will contribute to revealing the pathogenesis of cleft palate.

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 ENU-induced mutation in Twist1 transactivation domain causes hindlimb polydactyly with complete penetrance and dominant-negatively impairs E2A-dependent transcription

Twist1 encodes a basic helix-loop-helix transcription factor (TF), which forms homodimer or heterodimer with other TFs, like E2A, to regulate target genes' expression. Mutations in TWIST1 are associated with Saethre-Chotzen syndrome (SCS), a rare congenital disorder characterized with osteogenesis abnormalities. However, how dysfunction of TWIST1 leads to SCS is still largely unknown. Here, using an unbiased ENU-induced mutagenesis screening, we identified a novel Twist1 mutation and the mutant mouse phenocopies some features of SCS in a dominant manner. Physically, our mutation p.F191S lies at the edge of a predicted α-helix in Twist1 transactivation (TA) domain. Adjacent to F191, a consecutive three-residue (AFS) has been hit by 3 human and 2 mouse disease-associated mutations, including ours. Unlike previously reported mouse null and p.S192P alleles that lead to hindlimb polydactyly with incomplete penetrance but a severe craniofacial malformation, our p.F191S causes the polydactyly (84.2% bilateral and 15.8% unilateral) with complete penetrance but a mild craniofacial malformation. Consistent with the higher penetrance, p.F191S has stronger impairment on E2A-dependent transcription than p.S192P. Although human p.A186T and mouse p.S192P disease mutations are adjacent to ours, these three mutations function differently to impair the E2A-dependent transcription. Unlike p.A186T and p.S192S that disturb local protein conformation and unstabilize the mutant proteins, p.F191S keeps the mutant protein stable and its interaction with E2A entire. Therefore, we argue that p.F191S we identified acts in a dominant-negative manner to impair E2A-dependent transcription and to cause the biological consequences. In addition, the mutant mouse we provided here could be an additional and valuable model for better understanding the disease mechanisms underlying SCS caused by TWIST1 dysfunction.

Mice with a Brd4 Mutation Represent a New Model of Nephrocalcinosis

Nephrolithiasis (NL) and nephrocalcinosis (NC), which comprise renal calcification of the collecting system and parenchyma, respectively, have a multifactorial etiology with environmental and genetic determinants and affect ∼10% of adults by age 70 years. Studies of families with hereditary NL and NC have identified >30 causative genes that have increased our understanding of extracellular calcium homeostasis and renal tubular transport of calcium. However, these account for <20% of the likely genes that are involved, and to identify novel genes for renal calcification disorders, we investigated 1745 12-month-old progeny from a male mouse that had been treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) for radiological renal opacities. This identified a male mouse with renal calcification that was inherited as an autosomal dominant trait with >80% penetrance in 152 progeny. The calcification consisted of calcium phosphate deposits in the renal papillae and was associated with the presence of the urinary macromolecules osteopontin and Tamm-Horsfall protein, which are features found in Randall's plaques of patients with NC. Genome-wide mapping located the disease locus to a ∼30 Mbp region on chromosome 17A3.3-B3 and whole-exome sequence analysis identified a heterozygous mutation, resulting in a missense substitution (Met149Thr, M149T), in the bromodomain-containing protein 4 (BRD4). The mutant heterozygous (Brd4^+/M149T ) mice, when compared with wild-type (Brd4^+/+ ) mice, were normocalcemic and normophosphatemic, with normal urinary excretions of calcium and phosphate, and had normal bone turnover markers. BRD4 plays a critical role in histone modification and gene transcription, and cDNA expression profiling, using kidneys from Brd4^+/M149T and Brd4^+/+ mice, revealed differential expression of genes involved in vitamin D metabolism, cell differentiation, and apoptosis. Kidneys from Brd4^+/M149T mice also had increased apoptosis at sites of calcification within the renal papillae. Thus, our studies have established a mouse model, due to a Brd4 Met149Thr mutation, for inherited NC. © 2019 American Society for Bone and Mineral Research.

Histone methyltransferase PRDM9 is not essential for meiosis in male mice

A hallmark of meiosis is the rearrangement of parental alleles to ensure genetic diversity in the gametes. These chromosome rearrangements are mediated by the repair of programmed DNA double-strand breaks (DSBs) as genetic crossovers between parental homologs. In mice, humans, and many other mammals, meiotic DSBs occur primarily at hotspots, determined by sequence-specific binding of the PRDM9 protein. Without PRDM9, meiotic DSBs occur near gene promoters and other functional sites. Studies in a limited number of mouse strains showed that functional PRDM9 is required to complete meiosis, but despite its apparent importance, Prdm9 has been repeatedly lost across many animal lineages. Both the reason for mouse sterility in the absence of PRDM9 and the mechanism by which Prdm9 can be lost remain unclear. Here, we explore whether mice can tolerate the loss of Prdm9 By generating Prdm9 functional knockouts in an array of genetic backgrounds, we observe a wide range of fertility phenotypes and ultimately demonstrate that PRDM9 is not required for completion of male meiosis. Although DSBs still form at a common subset of functional sites in all mice lacking PRDM9, meiotic outcomes differ substantially. We speculate that DSBs at functional sites are difficult to repair as a crossover and that by increasing the efficiency of crossover formation at these sites, genetic modifiers of recombination rates can allow for meiotic progression. This model implies that species with a sufficiently high recombination rate may lose Prdm9 yet remain fertile.

Spontaneous Posterior Segment Vascular Disease Phenotype of a Mouse Model, rnv3, Is Dependent on the Crb1rd8 Allele

Purpose

To determine the molecular basis of lesion development in a murine model of spontaneous retinal vascularization, rnv3 (retinal vascularization 3, aka JR5558).

Methods

Disease progression of rnv3 was examined in longitudinal studies by clinical evaluation, electroretinography (ERG) and light microscopy analyses. The chromosomal position for the recessive rnv3 mutation was determined by DNA pooling and genome-wide linkage analysis. The causative mutation was discovered by comparison of whole exome sequences of rnv3 mutant and wild-type (WT) controls. In order to confirm the causative mutation, transcription activator-like effector nuclease (TALEN)-mediated oligonucleotide directed repair (ODR) was utilized to correct the mutant allele. Phenotypic correction was assessed by fundus imaging and optical coherence tomography of live mice.

Results

rnv3 exhibits early-onset, multifocal depigmented retinal lesions observable by fundus examination starting at 18 days of age. The retinal lesions are associated with fluorescein leakage around 25 days of age, with peak leakage at about 4 weeks of age. ERG responses deteriorate as rnv3 mutants age, concomitant with progressive photoreceptor disruption and loss that is observable by histology. Genetic analysis localized rnv3 to mouse chromosome (Chr) 1. By high throughput sequencing of a whole exome capture library of a rnv3/rnv3 mutant and subsequent sequence analysis, a single base deletion (del) in the Crb1 [crumbs family member 1] gene, which was previously reported to cause retinal degeneration 8, was identified. The TALEN-mediated ODR rescued the posterior segment vascularization phenotype; heterozygous Crb1^rd8+em1Boc/Crb1^rd8 and homozygous Crb1^rd8+em1Boc/Crb1^rd8+em1Boc mice showed a normal retinal phenotype. Additionally, six novel disruptions of Crb1 that were generated through aberrant non-homologous end joining induced by TALEN exhibited variable levels of vascularization, suggesting allelic effects.

Conclusions

The rnv3 model and the models of six novel disruptions of Crb1 are all reliable, novel mouse models for the study of both early and late events associated with posterior segment vascularization and can also be used to test the effects of pharmacological targets for treating human ocular vascular disorders. Further study of these models may provide a greater understanding about how different Crb1 alleles result in aberrant angiogenesis.

Mouse NC/Jic strain provides novel insights into host genetic factors for malaria research

Malaria is caused by Plasmodium parasites and is one of the most life-threatening infectious diseases in humans. Infection can result in severe complications such as cerebral malaria, acute lung injury/acute respiratory distress syndrome, and acute renal injury. These complications are mainly caused by P. falciparum infection and are major causes of death associated with malaria. There are a few species of rodent-infective malaria parasites, and mice infected with such parasites are now widely used for screening candidate drugs and vaccines and for studying host immune responses and pathogenesis associated with disease-related complications. We found that mice of the NC/Jic strain infected with rodent malarial parasites exhibit distinctive disease-related complications such as cerebral malaria and nephrotic syndrome, in addition to a rapid increase in parasitemia. Here, we focus on the analysis of host genetic factors that affect malarial pathogenesis and describe the characteristic features, utility, and future prospects for exploitation of the NC/Jic strain as a novel mouse model for malaria research.

Mutation in NADPH oxidase 3 (NOX3) impairs SHH signaling and increases cerebellar neural stem/progenitor cell proliferation

Abnormalities in cerebellar structure and function may cause ataxia, a neurological dysfunction of motor coordination. In the course of the present study, we characterized a mutant mouse lineage with an ataxia-like phenotype. We localized the mutation on chromosome 17 and mapped it to position 1534 of the Nox3 gene, resulting in p.Asn64Tyr change. The primary defect observed in Nox3^eqlb mice was increased proliferation of cerebellar granule cell precursors (GCPs). cDNA microarray comparing Nox3^eqlb and BALB/c neonatal cerebellum revealed changes in the expression of genes involved in the control of cell proliferation. Nox3^eqlb GCPs and NSC produce higher amounts of reactive oxygen species (ROS) and upregulate the expression of SHH target genes, such as Gli1-3 and Ccnd1 (CyclinD1). We hypothesize that this new mutation is responsible for an increase in proliferation via stimulation of the SHH pathway. We suggest this mutant mouse lineage as a new model to investigate the role of ROS in neuronal precursor cell proliferation.

An exonic splicing enhancer mutation in DUOX2 causes aberrant alternative splicing and severe congenital hypothyroidism in Bama pigs

Pigs share many similarities with humans in terms of anatomy, physiology and genetics, and have long been recognized as important experimental animals in biomedical research. Using an N-ethyl-N-nitrosourea (ENU) mutagenesis screen, we previously identified a large number of pig mutants, which could be further established as human disease models. However, the identification of causative mutations in large animals with great heterogeneity remains a challenging endeavor. Here, we select one pig mutant, showing congenital nude skin and thyroid deficiency in a recessive inheritance pattern. We were able to efficiently map the causative mutation using family-based genome-wide association studies combined with whole-exome sequencing and a small sample size. A loss-of-function variant (c.1226 A>G) that resulted in a highly conserved amino acid substitution (D409G) was identified in the DUOX2 gene. This mutation, located within an exonic splicing enhancer motif, caused aberrant splicing of DUOX2 transcripts and resulted in lower H₂O₂ production, which might cause a severe defect in thyroid hormone production. Our findings suggest that exome sequencing is an efficient way to map causative mutations and that DUOX2^D409G/D409G mutant pigs could be a potential large animal model for human congenital hypothyroidism.

Summary: Here, we show that an exonic splicing enhancer variant in DUOX2 (c.1226 A>G) causes aberrant splicing of DUOX2 transcripts, resulting in lower H₂O₂ production, to cause severe congenital hypothyroidism in Bama pigs.

Whole exome sequencing of ENU-induced thrombosis modifier mutations in the mouse

Although the Factor V Leiden (FVL) gene variant is the most prevalent genetic risk factor for venous thrombosis, only 10% of FVL carriers will experience such an event in their lifetime. To identify potential FVL modifier genes contributing to this incomplete penetrance, we took advantage of a perinatal synthetic lethal thrombosis phenotype in mice homozygous for FVL (F5L/L) and haploinsufficient for tissue factor pathway inhibitor (Tfpi+/-) to perform a sensitized dominant ENU mutagenesis screen. Linkage analysis conducted in the 3 largest pedigrees generated from the surviving F5L/L Tfpi+/- mice ('rescues') using ENU-induced coding variants as genetic markers was unsuccessful in identifying major suppressor loci. Whole exome sequencing was applied to DNA from 107 rescue mice to identify candidate genes enriched for ENU mutations. A total of 3,481 potentially deleterious candidate ENU variants were identified in 2,984 genes. After correcting for gene size and multiple testing, Arl6ip5 was identified as the most enriched gene, though not reaching genome-wide significance. Evaluation of CRISPR/Cas9 induced loss of function in the top 6 genes failed to demonstrate a clear rescue phenotype. However, a maternally inherited (not ENU-induced) de novo mutation (Plcb4R335Q) exhibited significant co-segregation with the rescue phenotype (p = 0.003) in the corresponding pedigree. Thrombosis suppression by heterozygous Plcb4 loss of function was confirmed through analysis of an independent, CRISPR/Cas9-induced Plcb4 mutation (p = 0.01).

Let's think again about using mammalian temperature-sensitive mutants to investigate functional molecules-The perspectives from the studies on three mutants showing chromosome instability

This review evaluates the use of temperature-sensitive (ts) mutants to investigate functional molecules in mammalian cells. A series of studies were performed in which mammalian cells expressing functional molecules were isolated from ts mutants using complementation by the introduction and expression of the responsible protein tagged with the green fluorescent protein. The results showed that chromosome instability and cell-cycle arrest were caused by ts defects in the following three molecules: the largest subunit of RNA polymerase II, a protein involved in splicing, and ubiquitin-activating enzyme. The cells expressing functional protein were then isolated by introducing the responsible gene tagged with the green fluorescent protein to complement the ts phenotype. These cells proved to be useful in analyzing the dynamics of RNA polymerase II in living cells. Analyses of the functional interaction between proteins involved in splicing were also useful in the investigation of ts mutants and their derivatives. In addition, these cells demonstrated the functional localization of ubiquitin-activating enzyme in the nucleus. Mammalian ts mutants continue to show great potential to aid in understanding the functions of the essential molecules in cells. Therefore, it is highly important that studies on the identification and characterization of the genes responsible for the phenotype of a mutant are carried out.

Defective immuno- and thymoproteasome assembly causes severe immunodeficiency

By N-ethyl-N-nitrosourea (ENU) mutagenesis, we generated the mutant mouse line TUB6 that is characterised by severe combined immunodeficiency (SCID) and systemic sterile autoinflammation in homozygotes, and a selective T cell defect in heterozygotes. The causative missense point mutation results in the single amino acid exchange G170W in multicatalytic endopeptidase complex subunit-1 (MECL-1), the β2i-subunit of the immuno- and thymoproteasome. Yeast mutagenesis and crystallographic data suggest that the severe TUB6-phenotype compared to the MECL-1 knockout mouse is caused by structural changes in the C-terminal appendage of β2i that prevent the biogenesis of immuno- and thymoproteasomes. Proteasomes are essential for cell survival, and defective proteasome assembly causes selective death of cells expressing the mutant MECL-1, leading to the severe immunological phenotype. In contrast to the immunosubunits β1i (LMP2) and β5i (LMP7), mutations in the gene encoding MECL-1 have not yet been assigned to human disorders. The TUB6 mutant mouse line exemplifies the involvement of MECL-1 in immunopathogenesis and provides the first mouse model for primary immuno- and thymoproteasome-associated immunodeficiency that may also be relevant in humans.

The developmental biology of genetic Notch disorders

Notch signaling regulates a vast array of crucial developmental processes. It is therefore not surprising that mutations in genes encoding Notch receptors or ligands lead to a variety of congenital disorders in humans. For example, loss of function of Notch results in Adams-Oliver syndrome, Alagille syndrome, spondylocostal dysostosis and congenital heart disorders, while Notch gain of function results in Hajdu-Cheney syndrome, serpentine fibula polycystic kidney syndrome, infantile myofibromatosis and lateral meningocele syndrome. Furthermore, structure-abrogating mutations in NOTCH3 result in CADASIL. Here, we discuss these human congenital disorders in the context of known roles for Notch signaling during development. Drawing on recent analyses by the exome aggregation consortium (EXAC) and on recent studies of Notch signaling in model organisms, we further highlight additional Notch receptors or ligands that are likely to be involved in human genetic diseases.

Phylogenomic Insights into Mouse Evolution Using a Pseudoreference Approach

Comparative genomic studies are now possible across a broad range of evolutionary timescales, but the generation and analysis of genomic data across many different species still present a number of challenges. The most sophisticated genotyping and down-stream analytical frameworks are still predominantly based on comparisons to high-quality reference genomes. However, established genomic resources are often limited within a given group of species, necessitating comparisons to divergent reference genomes that could restrict or bias comparisons across a phylogenetic sample. Here, we develop a scalable pseudoreference approach to iteratively incorporate sample-specific variation into a genome reference and reduce the effects of systematic mapping bias in downstream analyses. To characterize this framework, we used targeted capture to sequence whole exomes (∼54 Mbp) in 12 lineages (ten species) of mice spanning the Mus radiation. We generated whole exome pseudoreferences for all species and show that this iterative reference-based approach improved basic genomic analyses that depend on mapping accuracy while preserving the associated annotations of the mouse reference genome. We then use these pseudoreferences to resolve evolutionary relationships among these lineages while accounting for phylogenetic discordance across the genome, contributing an important resource for comparative studies in the mouse system. We also describe patterns of genomic introgression among lineages and compare our results to previous studies. Our general approach can be applied to whole or partitioned genomic data and is easily portable to any system with sufficient genomic resources, providing a useful framework for phylogenomic studies in mice and other taxa.

Identification of a Novel ENU-Induced Mutation in Mouse Tbx1 Linked to Human DiGeorge Syndrome

The patients with DiGeorge syndrome (DGS), caused by deletion containing dozens of genes in chromosome 22, often carry cardiovascular problem and hearing loss associated with chronic otitis media. Inside the deletion region, a transcription factor TBX1 was highly suspected. Furthermore, similar DGS phenotypes were found in the Tbx1 heterozygous knockout mice. Using ENU-induced mutagenesis and G1 dominant screening strategy, here we identified a nonsynonymous mutation p.W118R in T-box of TBX1, the DNA binding domain for transcription activity. The mutant mice showed deficiency of inner ear functions, including head tossing and circling, plus increased hearing threshold determined by audiometry. Therefore, our result further confirms the pathogenic basis of Tbx1 in DGS, points out the crucial role of DNA binding activity of TBX1 for the ear function, and provides additional animal model for studying the DGS disease mechanisms.

Design and application of a target capture sequencing of exons and conserved non-coding sequences for the rat

Background

Target capture sequencing is an efficient approach to directly identify the causative mutations of genetic disorders. To apply this strategy to laboratory rats exhibiting various phenotypes, we developed a novel target capture probe set, TargetEC (target capture for exons and conserved non-coding sequences), which can identify mutations not only in exonic regions but also in conserved non-coding sequences and thus can detect regulatory mutations.

Results

TargetEC covers 1,078,129 regions spanning 146.8 Mb of the genome. We applied TargetEC to four inbred rat strains (WTC/Kyo, WTC-swh/Kyo, PVG/Seac, and KFRS4/Kyo) maintained by the National BioResource Project for the Rat in Japan, and successfully identified mutations associated with these phenotypes, including one mutation detected in a conserved non-coding sequence.

Conclusions

The method developed in this study can be used to efficiently identify regulatory mutations, which cannot be detected using conventional exome sequencing, and will help to deepen our understanding of the relationships between regulatory mutations and associated phenotypes.

Electronic supplementary material

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

Retinal Pigment Epithelium Atrophy 1 (rpea1): A New Mouse Model With Retinal Detachment Caused by a Disruption of Protein Kinase C, θ

Purpose

Retinal detachments (RDs), a separation of the light-sensitive tissue of the retina from its supporting layers in the posterior eye, isolate retinal cells from their normal supply of nourishment and can lead to their deterioration and death. We identified a new, spontaneous murine model of exudative retinal detachment, nm3342 (new mutant 3342, also referred to as rpea1: retinal pigment epithelium atrophy 1), which we characterize herein.

Methods

The chromosomal position for the recessive nm3342 mutation was determined by DNA pooling, and the causative mutation was discovered by comparison of whole exome sequences of mutant and wild-type controls. The effects of the mutation were examined in longitudinal studies by clinical evaluation, electroretinography (ERG), light microscopy, and marker and Western blot analyses.

Results

New mutant 3342, nm3342, also referred to as rpea1, causes an early-onset, complete RD on the ABJ/LeJ strain background, and central exudative RD and late-onset RPE atrophy on the C57BL/6J background. The ERG responses were normal at 2 months of age but deteriorate as mice age, concomitant with progressive pan-retinal photoreceptor loss. Genetic analysis localized rpea1 to mouse chromosome 2. By high-throughput sequencing of a whole exome capture library of an rpea1/rpea1 mutant and subsequent sequence analysis, a splice donor site mutation in the Prkcq (protein kinase C, θ) gene, was identified, leading to a skipping of exon 6, frame shift and premature termination. Homozygotes with a Prkcq-targeted null allele (Prkcq^tm1Litt) have similar retinal phenotypes as homozygous rpea1 mice. We determined that the PKCθ protein is abundant in the lateral surfaces of RPE cells and colocalizes with both tight and adherens junction proteins. Phalloidin-stained RPE whole mounts showed abnormal RPE cell morphology with aberrant actin ring formation.

Conclusions

The homozygous Prkcq^rpea1 and the null Prkcq^tm1Litt mutants are reliable novel mouse models of RD and can also be used to study the effects of the disruption of PRKCQ (PKCθ) signaling in RPE cells.

Cilia gene mutations cause atrioventricular septal defects by multiple mechanisms

Atrioventricular septal defects (AVSDs) are a common severe form of congenital heart disease (CHD). In this study we identified deleterious non-synonymous mutations in two cilia genes, Dnah11 and Mks1, in independent N-ethyl-N-nitrosourea-induced mouse mutant lines with heritable recessive AVSDs by whole-exome sequencing. Cilia are required for left/right body axis determination and second heart field (SHF) Hedgehog (Hh) signaling, and we find that cilia mutations affect these requirements differentially. Dnah11^avc4 did not disrupt SHF Hh signaling and caused AVSDs only concurrently with heterotaxy, a left/right axis abnormality. In contrast, Mks1^avc6 disrupted SHF Hh signaling and caused AVSDs without heterotaxy. We performed unbiased whole-genome SHF transcriptional profiling and found that cilia motility genes were not expressed in the SHF whereas cilia structural and signaling genes were highly expressed. SHF cilia gene expression predicted the phenotypic concordance between AVSDs and heterotaxy in mice and humans with cilia gene mutations. A two-step model of cilia action accurately predicted the AVSD/heterotaxyu phenotypic expression pattern caused by cilia gene mutations. We speculate that cilia gene mutations contribute to both syndromic and non-syndromic AVSDs in humans and provide a model that predicts the phenotypic consequences of specific cilia gene mutations.

Contrasting Levels of Molecular Evolution on the Mouse X Chromosome

The mammalian X chromosome has unusual evolutionary dynamics compared to autosomes. Faster-X evolution of spermatogenic protein-coding genes is known to be most pronounced for genes expressed late in spermatogenesis, but it is unclear if these patterns extend to other forms of molecular divergence. We tested for faster-X evolution in mice spanning three different forms of molecular evolution-divergence in protein sequence, gene expression, and DNA methylation-across different developmental stages of spermatogenesis. We used FACS to isolate individual cell populations and then generated cell-specific transcriptome profiles across different stages of spermatogenesis in two subspecies of house mice (Mus musculus), thereby overcoming a fundamental limitation of previous studies on whole tissues. We found faster-X protein evolution at all stages of spermatogenesis and faster-late protein evolution for both X-linked and autosomal genes. In contrast, there was less expression divergence late in spermatogenesis (slower late) on the X chromosome and for autosomal genes expressed primarily in testis (testis-biased). We argue that slower-late expression divergence reflects strong regulatory constraints imposed during this critical stage of sperm development and that these constraints are particularly acute on the tightly regulated sex chromosomes. We also found slower-X DNA methylation divergence based on genome-wide bisulfite sequencing of sperm from two species of mice (M. musculus and M. spretus), although it is unclear whether slower-X DNA methylation reflects development constraints in sperm or other X-linked phenomena. Our study clarifies key differences in patterns of regulatory and protein evolution across spermatogenesis that are likely to have important consequences for mammalian sex chromosome evolution, male fertility, and speciation.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

2016: A 'Mitochondria' Odyssey

The integration of the many roles of mitochondria in cellular function and the contribution of mitochondrial dysfunction to disease are major areas of research. Within this realm, the roles of mitochondria in immune defense, epigenetics, and stem cell (SC) development have recently come into the spotlight. With new understanding, mitochondria may bring together these seemingly unrelated fields, a crucial process in treatment and prevention for various diseases. In this review we describe novel findings in these three arenas, discussing the significance of the interplay between mitochondria and the cell nucleus in response to environmental cues. While we optimistically anticipate that further research in these areas can have a profound impact on disease management, we also bring forth some of the key questions and challenges that remain.

Toward the most ideal case-control design with related and unrelated dogs in whole-exome sequencing studies

With the recent development of whole-exome sequencing enrichment designs for the dog, a novel tool for disease-association studies became available. The aim of disease-association studies is to identify one or a very limited number of putative causal variants or genes from the large pool of genetic variation. To maximize the efficiency of these studies and to provide some directions of what to expect, we evaluated the effect on variant reduction for various combinations of cases and controls for both dominant and recessive types of inheritance assuming variable degrees of penetrance and detectance. In this study, variant data of 14 dogs (13 Labrador Retrievers and one Dogue de Bordeaux), obtained by whole-exome sequencing, were analyzed. In the filtering process, we found that unrelated dogs from the same breed share up to 70% of their variants, which is likely a consequence of the breeding history of the dog. For the designs tested with unrelated dogs, combining two cases and two controls gave the best result. These results were improved further by adding closely related dogs. Reduced penetrance and/or detectance has a drastic effect on the efficiency and is likely to have a profound effect on the sample size needed to elucidate the causal variant. Overall, we demonstrated that sequencing a small number of dogs results in a marked reduction of variants that are likely sufficient to pinpoint causal variants or genes.

An heuristic filtering tool to identify phenotype-associated genetic variants applied to human intellectual disability and canine coat colors

Background

Identification of one or several disease causing variant(s) from the large collection of variants present in an individual is often achieved by the sequential use of heuristic filters. The recent development of whole exome sequencing enrichment designs for several non-model species created the need for a species-independent, fast and versatile analysis tool, capable of tackling a wide variety of standard and more complex inheritance models. With this aim, we developed “Mendelian”, an R-package that can be used for heuristic variant filtering.

Results

The R-package Mendelian offers fast and convenient filters to analyze putative variants for both recessive and dominant models of inheritance, with variable degrees of penetrance and detectance. Analysis of trios is supported. Filtering against variant databases and annotation of variants is also included. This package is not species specific and supports parallel computation. We validated this package by reanalyzing data from a whole exome sequencing experiment on intellectual disability in humans. In a second example, we identified the mutations responsible for coat color in the dog. This is the first example of whole exome sequencing without prior mapping in the dog.

Conclusion

We developed an R-package that enables the identification of disease-causing variants from the long list of variants called in sequencing experiments. The software and a detailed manual are available at https://github.com/BartBroeckx/Mendelian.

Electronic supplementary material

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

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.

Identification of putative SNPs in progressive retinal atrophy affected Canis lupus familiaris using exome sequencing

Progressive retinal atrophy (PRA) is one of the major causes of retinal photoreceptor cell degeneration in canines. The inheritance pattern of PRA is autosomal recessive and genetically heterogeneous. Here, using targeted sequencing technology, we have performed exome sequencing of 10 PRA-affected (Spitz=7, Cocker Spaniel=1, Lhasa Aphso=1 and Spitz-Labrador cross breed=1) and 6 normal (Spitz=5, Cocker Spaniel=1) dogs. The high-throughput sequencing using 454-Roche Titanium sequencer generated about 2.16 Giga bases of raw data. Initially, we have successfully identified 25,619 single nucleotide polymorphisms (SNPs) that passed the stringent SNP calling parameters. Further, we performed association study on the cohort, and the highly significant (0.001) associations were short-listed and investigated in-depth. Out of the 171 significant SNPs, 113 were previously unreported. Interestingly, six among them were non-synonymous coding (NSC) SNPs, which includes CPPED1 A>G (p.M307V), PITRM1 T>G (p.S715A), APP G>A (p.T266M), RNF213 A>G (p.V1482A), C>A (p.V1456L), and SLC46A3 G>A (p.R168Q). On the other hand, 35 out of 113 unreported SNPs were falling in regulatory regions such as 3'-UTR, 5'-UTR, etc. In-depth bioinformatics analysis revealed that majority of NSC SNPs have damaging effect and alter protein stability. This study highlighted the genetic markers associated with PRA, which will help to develop genetic assay-based screening in effective breeding.

New technologies for DNA analysis--a review of the READNA Project

The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) project received funding from the European Commission for 41/2 years. The objectives of the project revolved around technological developments in nucleic acid analysis. The project partners have discovered, created and developed a huge body of insights into nucleic acid analysis, ranging from improvements and implementation of current technologies to the most promising sequencing technologies that constitute a 3(rd) and 4(th) generation of sequencing methods with nanopores and in situ sequencing, respectively.

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.

Deleted in Breast Cancer 1 Suppresses B Cell Activation through RelB and Is Regulated by IKKα Phosphorylation

Alternative NF-κB signaling is crucial for B cell activation and Ig production, and it is mainly regulated by the inhibitor of κ B kinase (IKK) regulatory complex. Dysregulation of alternative NF-κB signaling in B cells could therefore lead to hyperactive B cells and Ig overproduction. In our previous, study we found that deleted in breast cancer 1 (DBC1) is a suppressor of the alternative NF-κB pathway to attenuate B cell activation. In this study, we report that loss of DBC1 results in spontaneous overproduction of Ig in mice after 10 mo of age. Using a double mutant genetic model, we confirm that DBC1 suppresses B cell activation through RelB inhibition. At the molecular level, we show that DBC1 interacts with alternative NF-κB members RelB and p52 through its leucine zipper domain. In addition, phosphorylation of DBC1 at its C terminus by IKKα facilitates its interaction with RelB and IKKα, indicating that DBC1-mediated suppression of alternative NF-κB is regulated by IKKα. Our results define the molecular mechanism of DBC1 inhibition of alternative NF-κB activation in suppressing B cell activation.

Multimer Formation Explains Allelic Suppression of PRDM9 Recombination Hotspots

Genetic recombination during meiosis functions to increase genetic diversity, promotes elimination of deleterious alleles, and helps assure proper segregation of chromatids. Mammalian recombination events are concentrated at specialized sites, termed hotspots, whose locations are determined by PRDM9, a zinc finger DNA-binding histone methyltransferase. Prdm9 is highly polymorphic with most alleles activating their own set of hotspots. In populations exhibiting high frequencies of heterozygosity, questions remain about the influences different alleles have in heterozygous individuals where the two variant forms of PRDM9 typically do not activate equivalent populations of hotspots. We now find that, in addition to activating its own hotspots, the presence of one Prdm9 allele can modify the activity of hotspots activated by the other allele. PRDM9 function is also dosage sensitive; Prdm9^+/- heterozygous null mice have reduced numbers and less active hotspots and increased numbers of aberrant germ cells. In mice carrying two Prdm9 alleles, there is allelic competition; the stronger Prdm9 allele can partially or entirely suppress chromatin modification and recombination at hotspots of the weaker allele. In cell cultures, PRDM9 protein variants form functional heteromeric complexes which can bind hotspots sequences. When a heteromeric complex binds at a hotspot of one PRDM9 variant, the other PRDM9 variant, which would otherwise not bind, can still methylate hotspot nucleosomes. We propose that in heterozygous individuals the underlying molecular mechanism of allelic suppression results from formation of PRDM9 heteromers, where the DNA binding activity of one protein variant dominantly directs recombination initiation towards its own hotspots, effectively titrating down recombination by the other protein variant. In natural populations with many heterozygous individuals, allelic competition will influence the recombination landscape.

During formation of sperm and eggs chromosomes exchange DNA in a process known as recombination, creating new combinations responsible for much of the enormous diversity in populations. In some mammals, including humans, the locations of recombination are chosen by a DNA-binding protein named PRDM9. Importantly, there are tens to hundreds of different variations of the Prdm9 gene (termed alleles), many of which are predicted to bind a unique DNA sequence. This high frequency of variation results in many individuals having two different copies of Prdm9, and several lines of evidence indicate that alleles compete to initiate recombination. In seeking to understand the mechanism of this competition we found that Prdm9 activity is sensitive to the number of gene copies present, suggesting that availability of this protein is a limiting factor during recombination. Moreover, we found that variant forms of PRDM9 protein can physically interact suggesting that when this happens one variant can influence which hotspots will become activated. Genetic crosses in mice support these observations; the presence of a dominant Prdm9 allele can completely suppress recombination at some locations. We conclude that allele-dominance of PRDM9 is a consequence of protein-protein interaction and competition for DNA binding in a limited pool of molecules, thus shaping the recombination landscape in natural populations.

Discovery and characterization of spontaneous mouse models of craniofacial dysmorphology

Craniofacial abnormalities are among the most common features of human genetic syndromes and disorders. The etiology of these conditions is often complex, influenced by both genetic context and the environment. Frequently, craniofacial abnormalities present as part of a syndrome with clear comorbid phenotypes, providing additional insight into mechanisms of the causative gene or pathway. The mouse has been a key tool in our understanding of the genetic mechanisms of craniofacial development and disease, and can provide excellent models for human craniofacial abnormalities. While powerful genetic engineering tools in the mouse have contributed significantly our understanding of craniofacial development and dysmorphology, forward genetic approaches provide an unbiased means to identify new genes and pathways. Moreover, spontaneous mutations can occur on any number of genetic backgrounds, potentially revealing critical genes that require a specific genetic context. Here we report discovery and phenotyping of 43 craniofacial mouse models, derived primarily from a screen for spontaneous mutations in production colonies at the Jackson Laboratory. We identify the causative gene for 33 lines, including novel genes in pathways not previously connected to craniofacial development, and novel alleles of known genes that present with unique phenotypes. Together with our detailed characterization, this work provides a valuable gene discovery resource for the craniofacial community, and a rich source of mouse models for further investigation.

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.

Large-scale identification of chemically induced mutations in Drosophila melanogaster

Forward genetic screens using chemical mutagens have been successful in defining the function of thousands of genes in eukaryotic model organisms. The main drawback of this strategy is the time-consuming identification of the molecular lesions causative of the phenotypes of interest. With whole-genome sequencing (WGS), it is now possible to sequence hundreds of strains, but determining which mutations are causative among thousands of polymorphisms remains challenging. We have sequenced 394 mutant strains, generated in a chemical mutagenesis screen, for essential genes on the Drosophila X chromosome and describe strategies to reduce the number of candidate mutations from an average of ∼3500 to 35 single-nucleotide variants per chromosome. By combining WGS with a rough mapping method based on large duplications, we were able to map 274 (∼70%) mutations. We show that these mutations are causative, using small 80-kb duplications that rescue lethality. Hence, our findings demonstrate that combining rough mapping with WGS dramatically expands the toolkit necessary for assigning function to genes.