Long-read sequencing for fast and robust identification of correct genome-edited alleles: PCR-based and Cas9 capture methods.

BACKGROUND: Recent developments in CRISPR/Cas9 genome-editing tools have facilitated the introduction of precise alleles, including genetic intervals spanning several kilobases, directly into the embryo. However, the introduction of donor templates, via homology directed repair, can be erroneous or incomplete and these techniques often produce mosaic founder animals. Thus, newly generated alleles must be verified at the sequence level across the targeted locus. Screening for the presence of the desired mutant allele using traditional sequencing methods can be challenging due to the size of the interval to be sequenced, together with the mosaic nature of founders. METHODOLOGY/PRINCIPAL FINDINGS: In order to help disentangle the genetic complexity of these animals, we tested the application of Oxford Nanopore Technologies long-read sequencing at the targeted locus and found that the achievable depth of sequencing is sufficient to offset the sequencing error rate associated with the technology used to validate targeted regions of interest. We have assembled an analysis workflow that facilitates interrogating the entire length of a targeted segment in a single read, to confirm that the intended mutant sequence is present in both heterozygous animals and mosaic founders. We used this workflow to compare the output of PCR-based and Cas9 capture-based targeted sequencing for validation of edited alleles. CONCLUSION: Targeted long-read sequencing supports in-depth characterisation of all experimental models that aim to produce knock-in or conditional alleles, including those that contain a mix of genome-edited alleles. PCR- or Cas9 capture-based modalities bring different advantages to the analysis.

TMEM161B regulates cerebral cortical gyration, Sonic Hedgehog signaling, and ciliary structure in the developing central nervous system.

Sonic hedgehog signaling regulates processes of embryonic development across multiple tissues, yet factors regulating context-specific Shh signaling remain poorly understood. Exome sequencing of families with polymicrogyria (disordered cortical folding) revealed multiple individuals with biallelic deleterious variants in TMEM161B, which encodes a multi-pass transmembrane protein of unknown function. Tmem161b null mice demonstrated holoprosencephaly, craniofacial midline defects, eye defects, and spinal cord patterning changes consistent with impaired Shh signaling, but were without limb defects, suggesting a CNS-specific role of Tmem161b. Tmem161b depletion impaired the response to Smoothened activation in vitro and disrupted cortical histogenesis in vivo in both mouse and ferret models, including leading to abnormal gyration in the ferret model. Tmem161b localizes non-exclusively to the primary cilium, and scanning electron microscopy revealed shortened, dysmorphic, and ballooned ventricular zone cilia in the Tmem161b null mouse, suggesting that the Shh-related phenotypes may reflect ciliary dysfunction. Our data identify TMEM161B as a regulator of cerebral cortical gyration, as involved in primary ciliary structure, as a regulator of Shh signaling, and further implicate Shh signaling in human gyral development.

Neuroplastin genetically interacts with Cadherin 23 and the encoded isoform Np55 is sufficient for cochlear hair cell function and hearing.

Mammalian hearing involves the mechanoelectrical transduction (MET) of sound-induced fluid waves in the cochlea. Essential to this process are the specialised sensory cochlear cells, the inner (IHCs) and outer hair cells (OHCs). While genetic hearing loss is highly heterogeneous, understanding the requirement of each gene will lead to a better understanding of the molecular basis of hearing and also to therapeutic opportunities for deafness. The Neuroplastin (Nptn) gene, which encodes two protein isoforms Np55 and Np65, is required for hearing, and homozygous loss-of-function mutations that affect both isoforms lead to profound deafness in mice. Here we have utilised several distinct mouse models to elaborate upon the spatial, temporal, and functional requirement of Nptn for hearing. While we demonstrate that both Np55 and Np65 are present in cochlear cells, characterisation of a Np65-specific mouse knockout shows normal hearing thresholds indicating that Np65 is functionally redundant for hearing. In contrast, we find that Nptn-knockout mice have significantly reduced maximal MET currents and MET channel open probabilities in mature OHCs, with both OHCs and IHCs also failing to develop fully mature basolateral currents. Furthermore, comparing the hearing thresholds and IHC synapse structure of Nptn-knockout mice with those of mice that lack Nptn only in IHCs and OHCs shows that the majority of the auditory deficit is explained by hair cell dysfunction, with abnormal afferent synapses contributing only a small proportion of the hearing loss. Finally, we show that continued expression of Neuroplastin in OHCs of adult mice is required for membrane localisation of Plasma Membrane Ca2+ ATPase 2 (PMCA2), which is essential for hearing function. Moreover, Nptn haploinsufficiency phenocopies Atp2b2 (encodes PMCA2) mutations, with heterozygous Nptn-knockout mice exhibiting hearing loss through genetic interaction with the Cdh23ahl allele. Together, our findings provide further insight to the functional requirement of Neuroplastin for mammalian hearing.

LAMA: automated image analysis for the developmental phenotyping of mouse embryos.

Advanced 3D imaging modalities, such as micro-computed tomography (micro-CT), have been incorporated into the high-throughput embryo pipeline of the International Mouse Phenotyping Consortium (IMPC). This project generates large volumes of raw data that cannot be immediately exploited without significant resources of personnel and expertise. Thus, rapid automated annotation is crucial to ensure that 3D imaging data can be integrated with other multi-dimensional phenotyping data. We present an automated computational mouse embryo phenotyping pipeline that harnesses the large amount of wild-type control data available in the IMPC embryo pipeline in order to address issues of low mutant sample number as well as incomplete penetrance and variable expressivity. We also investigate the effect of developmental substage on automated phenotyping results. Designed primarily for developmental biologists, our software performs image pre-processing, registration, statistical analysis and segmentation of embryo images. We also present a novel anatomical E14.5 embryo atlas average and, using it with LAMA, show that we can uncover known and novel dysmorphology from two IMPC knockout lines.

Cleft palate and minor metabolic disturbances in a mouse global Arl15 gene knockout.

ARL15, a small GTPase protein, was linked to metabolic traits in association studies. We aimed to test the Arl15 gene as a functional candidate for metabolic traits in the mouse. CRISPR/Cas9 germline knockout (KO) of Arl15 showed that homozygotes were postnatal lethal and exhibited a complete cleft palate (CP). Also, decreased cell migration was observed from Arl15 KO mouse embryonic fibroblasts (MEFs). Metabolic phenotyping of heterozygotes showed that females had reduced fat mass on a chow diet from 14 weeks of age. Mild body composition phenotypes were also observed in heterozygous mice on a high-fat diet (HFD)/low-fat diet (LFD). Females on a HFD showed reduced body weight, gonadal fat depot weight and brown adipose tissue (BAT) weight. In contrast, in the LFD group, females showed increased bone mineral density (BMD), while males showed a trend toward reduced BMD. Clinical biochemistry analysis of plasma on HFD showed transient lower adiponectin at 20 weeks of age in females. Urinary and plasma Mg2+ concentrations were not significantly different. Our phenotyping data showed that Arl15 is essential for craniofacial development. Adult metabolic phenotyping revealed potential roles in brown adipose tissue and bone development.

Ap2s1 mutation causes hypercalcaemia in mice and impairs interaction between calcium-sensing receptor and adaptor protein-2.

Adaptor protein 2 (AP2), a heterotetrameric complex comprising AP2α, AP2ß2, AP2µ2 and AP2σ2 subunits, is ubiquitously expressed and involved in endocytosis and trafficking of membrane proteins, such as the calcium-sensing receptor (CaSR), a G-protein coupled receptor that signals via Gα11. Mutations of CaSR, Gα11 and AP2σ2, encoded by AP2S1, cause familial hypocalciuric hypercalcaemia types 1-3 (FHH1-3), respectively. FHH3 patients have heterozygous AP2S1 missense Arg15 mutations (p.Arg15Cys, p.Arg15His or p.Arg15Leu) with hypercalcaemia, which may be marked and symptomatic, and occasional hypophosphataemia and osteomalacia. To further characterize the phenotypic spectrum and calcitropic pathophysiology of FHH3, we used CRISPR/Cas9 genome editing to generate mice harboring the AP2S1 p.Arg15Leu mutation, which causes the most severe FHH3 phenotype. Heterozygous (Ap2s1+/L15) mice were viable, and had marked hypercalcaemia, hypermagnesaemia, hypophosphataemia, and increases in alkaline phosphatase activity and fibroblast growth factor-23. Plasma 1,25-dihydroxyvitamin D was normal, and no alterations in bone mineral density or bone turnover were noted. Homozygous (Ap2s1L15/L15) mice invariably died perinatally. Co-immunoprecipitation studies showed that the AP2S1 p.Arg15Leu mutation impaired protein-protein interactions between AP2σ2 and the other AP2 subunits, and also with the CaSR. Cinacalcet, a CaSR positive allosteric modulator, decreased plasma calcium and parathyroid hormone concentrations in Ap2s1+/L15 mice, but had no effect on the diminished AP2σ2-CaSR interaction in vitro. Thus, our studies have established a mouse model that is representative for FHH3 in humans, and demonstrated that the AP2S1 p.Arg15Leu mutation causes a predominantly calcitropic phenotype, which can be ameliorated by treatment with cinacalcet.

Anticipating and Identifying Collateral Damage in Genome Editing.

Genome editing has powerful applications in research, healthcare, and agriculture. However, the range of possible molecular events resulting from genome editing has been underestimated and the technology remains unpredictable on, and away from, the target locus. This has considerable impact in providing a safe approach for therapeutic genome editing, agriculture, and other applications. This opinion article discusses how to anticipate and detect those editing events by a combination of assays to capture all possible genomic changes. It also discusses strategies for preventing unwanted effects, critical to appraise the benefit or risk associated with the use of the technology. Anticipating and verifying the result of genome editing are essential for the success for all applications.

Corrigendum: High-throughput discovery of novel developmental phenotypes.

This corrects the article DOI: 10.1038/nature19356.

Mouse mutant phenotyping at scale reveals novel genes controlling bone mineral density.

The genetic landscape of diseases associated with changes in bone mineral density (BMD), such as osteoporosis, is only partially understood. Here, we explored data from 3,823 mutant mouse strains for BMD, a measure that is frequently altered in a range of bone pathologies, including osteoporosis. A total of 200 genes were found to significantly affect BMD. This pool of BMD genes comprised 141 genes with previously unknown functions in bone biology and was complementary to pools derived from recent human studies. Nineteen of the 141 genes also caused skeletal abnormalities. Examination of the BMD genes in osteoclasts and osteoblasts underscored BMD pathways, including vesicle transport, in these cells and together with in silico bone turnover studies resulted in the prioritization of candidate genes for further investigation. Overall, the results add novel pathophysiological and molecular insight into bone health and disease.

Loss of O-GlcNAcase catalytic activity leads to defects in mouse embryogenesis.

O-GlcNAcylation is an essential post-translational modification that has been implicated in neurodevelopmental and neurodegenerative disorders. O-GlcNAcase (OGA), the sole enzyme catalyzing the removal of O-GlcNAc from proteins, has emerged as a potential drug target. OGA consists of an N-terminal OGA catalytic domain and a C-terminal pseudo histone acetyltransferase (HAT) domain with unknown function. To investigate phenotypes specific to loss of OGA catalytic activity and dissect the role of the HAT domain, we generated a constitutive knock-in mouse line, carrying a mutation of a catalytic aspartic acid to alanine. These mice showed perinatal lethality and abnormal embryonic growth with skewed Mendelian ratios after day E18.5. We observed tissue-specific changes in O-GlcNAc homeostasis regulation to compensate for loss of OGA activity. Using X-ray microcomputed tomography on late gestation embryos, we identified defects in the kidney, brain, liver, and stomach. Taken together, our data suggest that developmental defects during gestation may arise upon prolonged OGA inhibition specifically because of loss of OGA catalytic activity and independent of the function of the HAT domain.

Mice with endogenous TDP-43 mutations exhibit gain of splicing function and characteristics of amyotrophic lateral sclerosis.

TDP-43 (encoded by the gene TARDBP) is an RNA binding protein central to the pathogenesis of amyotrophic lateral sclerosis (ALS). However, how TARDBP mutations trigger pathogenesis remains unknown. Here, we use novel mouse mutants carrying point mutations in endogenous Tardbp to dissect TDP-43 function at physiological levels both in vitro and in vivo Interestingly, we find that mutations within the C-terminal domain of TDP-43 lead to a gain of splicing function. Using two different strains, we are able to separate TDP-43 loss- and gain-of-function effects. TDP-43 gain-of-function effects in these mice reveal a novel category of splicing events controlled by TDP-43, referred to as "skiptic" exons, in which skipping of constitutive exons causes changes in gene expression. In vivo, this gain-of-function mutation in endogenous Tardbp causes an adult-onset neuromuscular phenotype accompanied by motor neuron loss and neurodegenerative changes. Furthermore, we have validated the splicing gain-of-function and skiptic exons in ALS patient-derived cells. Our findings provide a novel pathogenic mechanism and highlight how TDP-43 gain of function and loss of function affect RNA processing differently, suggesting they may act at different disease stages.

Universal Southern blot protocol with cold or radioactive probes for the validation of alleles obtained by homologous recombination.

The widespread availability of recombineered vectors and gene targeted embryonic stem cells from large-scale repositories facilitates the generation of mouse models for functional genetic studies. Southern blotting validates the structure of these targeted alleles produced by homologous recombination, as well as indicating any additional integrations of the vector into the genome. Traditionally this technique employs radioactively-labelled probes; however, there are many laboratories that are restricted in their use of radioactivity. Here, we present a widely applicable protocol for Southern blot analysis using cold probes and alternative procedures employing radioactive probes. Furthermore, the probes are designed to recognise standardised regions of gene-targeting cassettes and so represent universally applicable reagents for assessing allelic integrity.

High-throughput discovery of novel developmental phenotypes.

Approximately one-third of all mammalian genes are essential for life. Phenotypes resulting from knockouts of these genes in mice have provided tremendous insight into gene function and congenital disorders. As part of the International Mouse Phenotyping Consortium effort to generate and phenotypically characterize 5,000 knockout mouse lines, here we identify 410 lethal genes during the production of the first 1,751 unique gene knockouts. Using a standardized phenotyping platform that incorporates high-resolution 3D imaging, we identify phenotypes at multiple time points for previously uncharacterized genes and additional phenotypes for genes with previously reported mutant phenotypes. Unexpectedly, our analysis reveals that incomplete penetrance and variable expressivity are common even on a defined genetic background. In addition, we show that human disease genes are enriched for essential genes, thus providing a dataset that facilitates the prioritization and validation of mutations identified in clinical sequencing efforts.

A novel knockout mouse for the small EDRK-rich factor 2 (Serf2) showing developmental and other deficits.

The small EDRK-rich factor 2 (SERF2) is a highly conserved protein that modifies amyloid fibre assembly in vitro and promotes protein misfolding. However, the role of SERF2 in regulating age-related proteotoxicity remains largely unexplored due to a lack of in vivo models. Here, we report the generation of Serf2 knockout mice using an ES cell targeting approach, with Serf2 knockout alleles being bred onto different defined genetic backgrounds. We highlight phenotyping data from heterozygous Serf2+/- mice, including unexpected male-specific phenotypes in startle response and pre-pulse inhibition. We report embryonic lethality in Serf2-/- null animals when bred onto a C57BL/6 N background. However, homozygous null animals were viable on a mixed genetic background and, remarkably, developed without obvious abnormalities. The Serf2 knockout mice provide a powerful tool to further investigate the role of SERF2 protein in previously unexplored pathophysiological pathways in the context of a whole organism.

Variability in Genome Editing Outcomes: Challenges for Research Reproducibility and Clinical Safety.

Genome editing tools have already revolutionized biomedical research and are also expected to have an important impact in the clinic. However, their extensive use in research has revealed much unpredictability, both off and on target, in the outcome of their application. We discuss the challenges associated with this unpredictability, both for research and in the clinic. For the former, an extensive validation of the model is essential. For the latter, potential unpredicted activity does not preclude the use of these tools but requires that molecular evidence to underpin the relevant risk:benefit evaluation is available. Safe and successful clinical application will also depend on the mode of delivery and the cellular context.

Microhomologies are prevalent at Cas9-induced larger deletions.

The CRISPR system is widely used in genome editing for biomedical research. Here, using either dual paired Cas9D10A nickases or paired Cas9 nuclease we characterize unintended larger deletions at on-target sites that frequently evade common genotyping practices. We found that unintended larger deletions are prevalent at multiple distinct loci on different chromosomes, in cultured cells and mouse embryos alike. We observed a high frequency of microhomologies at larger deletion breakpoint junctions, suggesting the involvement of microhomology-mediated end joining in their generation. In populations of edited cells, the distribution of larger deletion sizes is dependent on proximity to sgRNAs and cannot be predicted by microhomology sequences alone.

High-resolution µCT of a mouse embryo using a compact laser-driven X-ray betatron source.

In the field of X-ray microcomputed tomography (µCT) there is a growing need to reduce acquisition times at high spatial resolution (approximate micrometers) to facilitate in vivo and high-throughput operations. The state of the art represented by synchrotron light sources is not practical for certain applications, and therefore the development of high-brightness laboratory-scale sources is crucial. We present here imaging of a fixed embryonic mouse sample using a compact laser-plasma-based X-ray light source and compare the results to images obtained using a commercial X-ray µCT scanner. The radiation is generated by the betatron motion of electrons inside a dilute and transient plasma, which circumvents the flux limitations imposed by the solid or liquid anodes used in conventional electron-impact X-ray tubes. This X-ray source is pulsed (duration <30 fs), bright (>1010 photons per pulse), small (diameter <1 µm), and has a critical energy >15 keV. Stable X-ray performance enabled tomographic imaging of equivalent quality to that of the µCT scanner, an important confirmation of the suitability of the laser-driven source for applications. The X-ray flux achievable with this approach scales with the laser repetition rate without compromising the source size, which will allow the recording of high-resolution µCT scans in minutes.

Biallelic Variants in OTUD6B Cause an Intellectual Disability Syndrome Associated with Seizures and Dysmorphic Features.

Ubiquitination is a posttranslational modification that regulates many cellular processes including protein degradation, intracellular trafficking, cell signaling, and protein-protein interactions. Deubiquitinating enzymes (DUBs), which reverse the process of ubiquitination, are important regulators of the ubiquitin system. OTUD6B encodes a member of the ovarian tumor domain (OTU)-containing subfamily of deubiquitinating enzymes. Herein, we report biallelic pathogenic variants in OTUD6B in 12 individuals from 6 independent families with an intellectual disability syndrome associated with seizures and dysmorphic features. In subjects with predicted loss-of-function alleles, additional features include global developmental delay, microcephaly, absent speech, hypotonia, growth retardation with prenatal onset, feeding difficulties, structural brain abnormalities, congenital malformations including congenital heart disease, and musculoskeletal features. Homozygous Otud6b knockout mice were subviable, smaller in size, and had congenital heart defects, consistent with the severity of loss-of-function variants in humans. Analysis of peripheral blood mononuclear cells from an affected subject showed reduced incorporation of 19S subunits into 26S proteasomes, decreased chymotrypsin-like activity, and accumulation of ubiquitin-protein conjugates. Our findings suggest a role for OTUD6B in proteasome function, establish that defective OTUD6B function underlies a multisystemic human disorder, and provide additional evidence for the emerging relationship between the ubiquitin system and human disease.

A bioimage informatics platform for high-throughput embryo phenotyping.

High-throughput phenotyping is a cornerstone of numerous functional genomics projects. In recent years, imaging screens have become increasingly important in understanding gene-phenotype relationships in studies of cells, tissues and whole organisms. Three-dimensional (3D) imaging has risen to prominence in the field of developmental biology for its ability to capture whole embryo morphology and gene expression, as exemplified by the International Mouse Phenotyping Consortium (IMPC). Large volumes of image data are being acquired by multiple institutions around the world that encompass a range of modalities, proprietary software and metadata. To facilitate robust downstream analysis, images and metadata must be standardized to account for these differences. As an open scientific enterprise, making the data readily accessible is essential so that members of biomedical and clinical research communities can study the images for themselves without the need for highly specialized software or technical expertise. In this article, we present a platform of software tools that facilitate the upload, analysis and dissemination of 3D images for the IMPC. Over 750 reconstructions from 80 embryonic lethal and subviable lines have been captured to date, all of which are openly accessible at mousephenotype.org. Although designed for the IMPC, all software is available under an open-source licence for others to use and develop further. Ongoing developments aim to increase throughput and improve the analysis and dissemination of image data. Furthermore, we aim to ensure that images are searchable so that users can locate relevant images associated with genes, phenotypes or human diseases of interest.