BATCH-GE: Batch analysis of Next-Generation Sequencing data for genome editing assessment

Bmpr1aa modulates the severity of the skeletal phenotype in an fkbp10-deficient Bruck syndrome zebrafish model

Rare monogenic disorders often exhibit significant phenotypic variability among individuals sharing identical genetic mutations. Bruck syndrome (BS), a prime example, is characterized by bone fragility and congenital contractures, although with a pronounced variability among family members. BS arises from recessive biallelic mutations in FKBP10 or PLOD2. FKBP65, the protein encoded by FKBP10, collaborates with the LH2 enzyme (PLOD2) in type I collagen telopeptide lysine hydroxylation, crucial for collagen cross-linking. To identify potential modifier genes and to investigate the mechanistic role of FKBP10 in BS pathogenesis, we established an fkbp10a knockout zebrafish model. Mass-spectrometry analysis in fkbp10a-/- mutants revealed a generally decreased type I collagen lysyl hydroxylation, paralleled by a wide skeletal variability similar to human patients. Ultrastructural examination of the skeleton in severely affected mutants showed enlarged type I collagen fibrils and disturbed elastin layers. Whole-exome sequencing of 7 mildly and 7 severely affected mutant zebrafish siblings, followed by single nucleotide polymorphism-based linkage analysis, indicated a linked region on chromosome 13, which segregates with phenotypic severity. Transcriptome analysis identified 6 differentially expressed genes (DEGs) between mildly and severely affected mutants. The convergence of genes within the linked region and DEGs highlighted bmpr1aa as a potential modifier gene, as its reduced expression correlates with increased skeletal severity. In summary, our study provides deeper insights into the role of FKBP10 in BS pathogenesis. Additionally, we identified a pivotal gene that influences phenotypic severity in a zebrafish model of BS. These findings hold promise for novel treatments in the field of bone diseases.

Deletion upstream of MAB21L2 highlights the importance of evolutionarily conserved non-coding sequences for eye development

Anophthalmia, microphthalmia and coloboma (AMC) comprise a spectrum of developmental eye disorders, accounting for approximately 20% of childhood visual impairment. While non-coding regulatory sequences are increasingly recognised as contributing to disease burden, characterising their impact on gene function and phenotype remains challenging. Furthermore, little is known of the nature and extent of their contribution to AMC phenotypes. We report two families with variants in or near MAB21L2, a gene where genetic variants are known to cause AMC in humans and animal models. The first proband, presenting with microphthalmia and coloboma, has a likely pathogenic missense variant (c.338 G > C; p.[Trp113Ser]), segregating within the family. The second individual, presenting with microphthalmia, carries an ~ 113.5 kb homozygous deletion 19.38 kb upstream of MAB21L2. Modelling of the deletion results in transient small lens and coloboma as well as midbrain anomalies in zebrafish, and microphthalmia and coloboma in Xenopus tropicalis. Using conservation analysis, we identify 15 non-coding conserved elements (CEs) within the deleted region, while ChIP-seq data from mouse embryonic stem cells demonstrates that two of these (CE13 and 14) bind Otx2, a protein with an established role in eye development. Targeted disruption of CE14 in Xenopus tropicalis recapitulates an ocular coloboma phenotype, supporting its role in eye development. Together, our data provides insights into regulatory mechanisms underlying eye development and highlights the importance of non-coding sequences as a source of genetic diagnoses in AMC.

Quantifying allele-specific CRISPR editing activity with CRISPECTOR2.0

Off-target effects present a significant impediment to the safe and efficient use of CRISPR-Cas genome editing. Since off-target activity is influenced by the genomic sequence, the presence of sequence variants leads to varying on- and off-target profiles among different alleles or individuals. However, a reliable tool that quantifies genome editing activity in an allelic context is not available. Here, we introduce CRISPECTOR2.0, an extended version of our previously published software tool CRISPECTOR, with an allele-specific editing activity quantification option. CRISPECTOR2.0 enables reference-free, allele-aware, precise quantification of on- and off-target activity, by using de novo sample-specific single nucleotide variant (SNV) detection and statistical-based allele-calling algorithms. We demonstrate CRISPECTOR2.0 efficacy in analyzing samples containing multiple alleles and quantifying allele-specific editing activity, using data from diverse cell types, including primary human cells, plants, and an original extensive human cell line database. We identified instances where an SNV induced changes in the protospacer adjacent motif sequence, resulting in allele-specific editing. Intriguingly, differential allelic editing was also observed in regions carrying distal SNVs, hinting at the involvement of additional epigenetic factors. Our findings highlight the importance of allele-specific editing measurement as a milestone in the adaptation of efficient, accurate, and safe personalized genome editing.

Evolutionary origin of Hoxc13-dependent skin appendages in amphibians

Cornified skin appendages, such as hair and nails, are major evolutionary innovations of terrestrial vertebrates. Human hair and nails consist largely of special intermediate filament proteins, known as hair keratins, which are expressed under the control of the transcription factor Hoxc13. Here, we show that the cornified claws of Xenopus frogs contain homologs of hair keratins and the genes encoding these keratins are flanked by promoters in which binding sites of Hoxc13 are conserved. Furthermore, these keratins and Hoxc13 are co-expressed in the claw-forming epithelium of frog toe tips. Upon deletion of hoxc13, the expression of hair keratin homologs is abolished and the development of cornified claws is abrogated in X. tropicalis. These results indicate that Hoxc13-dependent expression of hair keratin homologs evolved already in stem tetrapods, presumably as a mechanism for protecting toe tips, and that this ancestral genetic program was coopted to the growth of hair in mammals.

Mutation Knock-in Methods Using Single-Stranded DNA and Gene Editing Tools in Zebrafish

Introduction or knock-in of precise genomic modifications remains one of the most important applications of CRISPR/Cas9 in all model systems including zebrafish. The most widely used type of donor template containing the desired modification is single-stranded DNA (ssDNA), either in the form of single-stranded oligodeoxynucleotides (ssODN) (<150 nucleotides (nt)) or as long ssDNA (lssDNA) molecules (up to about 2000 nt). Despite the challenges posed by DNA repair after DNA double-strand breaks, knock-in of precise mutations is relatively straightforward in zebrafish. Knock-in efficiency can be enhanced by careful donor template design, using lssDNA as template or tethering the donor template DNA to the Cas9-guide RNA complex. Other point mutation methods such as base editing and prime editing are starting to be applied in zebrafish and many other model systems. However, these methods may not always be sufficiently accessible or may have limited capacity to perform all desired mutation knock-ins which are possible with ssDNA-based knock-in methods. Thus, it is likely that there will be complementarity in the technologies used for generating precise mutants. Here, we review and describe a suite of CRISPR/Cas9 knock-in procedures utilizing ssDNA as the donor template in zebrafish, point out the potential challenges and suggest possible approaches for their solution ultimately leading to successful generation of precise mutant lines.

CRISPR/Cas9-gene editing approaches in plant breeding

CRISPR/Cas9 gene editing system is recently developed robust genome editing technology for accelerating plant breeding. Various modifications of this editing system have been established for adaptability in plant varieties as well as for its improved efficiency and portability. This review provides an in-depth look at the various strategies for synthesizing gRNAs for efficient delivery in plant cells, including chemical synthesis and in vitro transcription. It also covers traditional analytical tools and emerging developments in detection methods to analyze CRISPR/Cas9 mediated mutation in plant breeding. Additionally, the review outlines the various analytical tools which are used to detect and analyze CRISPR/Cas9 mediated mutations, such as next-generation sequencing, restriction enzyme analysis, and southern blotting. Finally, the review discusses emerging detection methods, including digital PCR and qPCR. Hence, CRISPR/Cas9 has great potential for transforming agriculture and opening avenues for new advancements in the system for gene editing in plants.

Mutations in the histone methyltransferase Ezh2 drive context-dependent leukemia in Xenopus tropicalis

CRISPR-mediated simultaneous targeting of candidate tumor suppressor genes in Xenopus tropicalis allows fast functional assessment of co-driver genes for various solid tumors. Genotyping of tumors that emerge in the mosaic mutant animals rapidly exposes the gene mutations under positive selection for tumor establishment. However, applying this simple approach to the blood lineage has not been attempted. Multiple hematologic malignancies have mutations in EZH2, encoding the catalytic subunit of the Polycomb Repressive Complex 2. Interestingly, EZH2 can act as an oncogene or a tumor suppressor, depending on cellular context and disease stage. We show here that mosaic CRISPR/Cas9 mediated ezh2 disruption in the blood lineage resulted in early and penetrant acute myeloid leukemia (AML) induction. While animals were co-targeted with an sgRNA that induces notch1 gain-of-function mutations, sequencing of leukemias revealed positive selection towards biallelic ezh2 mutations regardless of notch1 mutational status. Co-targeting dnm2, recurrently mutated in T/ETP-ALL, induced a switch from myeloid towards acute T-cell leukemia. Both myeloid and T-cell leukemias engrafted in immunocompromised hosts. These data underline the potential of Xenopus tropicalis for modeling human leukemia, where mosaic gene disruption, combined with deep amplicon sequencing of the targeted genomic regions, can rapidly and efficiently expose co-operating driver gene mutations.

A Digital PCR Method Based on Highly Specific Taq for Detecting Gene Editing and Mutations

Digital PCR (dPCR) has great potential for assessing gene editing or gene mutation due to its ability to independently inspect each DNA template in parallel. However, current dPCR methods use a fluorescence-labeled probe to detect gene variation events, and their ability to distinguish variated sequences from the wild-type sequence is limited by the probe's tolerance to mismatch. To address this, we have developed a novel dPCR method that uses a primer instead of a probe to sense gene variation. The enhanced Taq DNA polymerase in the PCR system has a high mismatch sensitivity, which enables our dPCR method to distinguish gene mutations from wild-type sequences. Compared to current dPCR methods, our method shows superior precision in assessing gene editing efficiency and single-base DNA mutation. This presents a promising opportunity to advance gene editing research and rare gene mutation detection.

Retained chromosomal integrity following CRISPR-Cas9-based mutational correction in human embryos

Human germline gene correction by targeted nucleases holds great promise for reducing mutation transmission. However, recent studies have reported concerning observations in CRISPR-Cas9-targeted human embryos, including mosaicism and loss of heterozygosity (LOH). The latter has been associated with either gene conversion or (partial) chromosome loss events. In this study, we aimed to correct a heterozygous basepair substitution in PLCZ1, related to infertility. In 36% of the targeted embryos that originated from mutant sperm, only wild-type alleles were observed. By performing genome-wide double-digest restriction site-associated DNA sequencing, integrity of the targeted chromosome (i.e., no deletions larger than 3 Mb or chromosome loss) was confirmed in all seven targeted GENType-analyzed embryos (mutant editing and absence of mutation), while short-range LOH events (shorter than 10 Mb) were clearly observed by single-nucleotide polymorphism assessment in two of these embryos. These results fuel the currently ongoing discussion on double-strand break repair in early human embryos, making a case for the occurrence of gene conversion events or partial template-based homology-directed repair.

Poly (A)-specific ribonuclease deficiency impacts oogenesis in zebrafish

Poly (A)-specific ribonuclease (PARN) is the most important 3'-5'exonuclease involved in the process of deadenylation, the removal of poly (A) tails of mRNAs. Although PARN is primarily known for its role in mRNA stability, recent studies suggest several other functions of PARN including a role in telomere biology, non-coding RNA maturation, trimming of miRNAs, ribosome biogenesis and TP53 function. Moreover, PARN expression is de-regulated in many cancers, including solid tumours and hematopoietic malignancies. To better understand the in vivo role of PARN, we used a zebrafish model to study the physiological consequences of Parn loss-of-function. Exon 19 of the gene, which partially codes for the RNA binding domain of the protein, was targeted for CRISPR-Cas9-directed genome editing. Contrary to the expectations, no developmental defects were observed in the zebrafish with a parn nonsense mutation. Intriguingly, the parn null mutants were viable and fertile, but turned out to only develop into males. Histological analysis of the gonads in the mutants and their wild type siblings revealed a defective maturation of gonadal cells in the parn null mutants. The results of this study highlight yet another emerging function of Parn, i.e., its role in oogenesis.

Multi-faceted CRISPR/Cas technological innovation aspects in the framework of 3P medicine

Since 2009, the European Association for Predictive, Preventive and Personalised Medicine (EPMA, Brussels) promotes the paradigm change from reactive approach to predictive, preventive, and personalized medicine (PPPM/3PM) to protect individuals in sub-optimal health conditions from the health-to-disease transition, to increase life-quality of the affected patient cohorts improving, therefore, ethical standards and cost-efficacy of healthcare to great benefits of the society at large. The gene-editing technology utilizing CRISPR/Cas gene-editing approach has demonstrated its enormous value as a powerful tool in a broad spectrum of bio/medical research areas. Further, CRISPR/Cas gene-editing system is considered applicable to primary and secondary healthcare, in order to prevent disease spread and to treat clinically manifested disorders, involving diagnostics of SARS-Cov-2 infection and experimental treatment of COVID-19. Although the principle of the proposed gene editing is simple and elegant, there are a lot of technological challenges and ethical considerations to be solved prior to its broadly scaled clinical implementation. This article highlights technological innovation beyond the state of the art, exemplifies current achievements, discusses unsolved technological and ethical problems, and provides clinically relevant outlook in the framework of 3PM.

CRISPR-GRANT: a cross-platform graphical analysis tool for high-throughput CRISPR-based genome editing evaluation

BACKGROUD

CRISPR/Cas is an efficient genome editing system that has been widely used for functional genetic studies and exhibits high potential in biomedical translational applications. Indel analysis has thus become one of the most common practices in the lab to evaluate DNA editing events generated by CRISPR/Cas. Several indel analysis tools have been reported, however, it is often required that users have certain bioinformatics training and basic command-line processing capability.

RESULTS

Here, we developed CRISPR-GRANT, a stand-alone graphical CRISPR indel analysis tool, which could be easily installed for multi-platforms, including Linux, Windows, and macOS. CRISPR-GRANT offered a straightforward GUI by simple click-and-run for genome editing analysis of single or pooled amplicons and one-step analysis for whole-genome sequencing without the need of data pre-processing, making it ideal for novice lab scientists. Moreover, it also exhibited shorter run-time compared with tools currently available.

CONCLUSION

Therefore, CRISPR-GRANT is a valuable addition to the current CRISPR toolkits that significantly lower the barrier for wet-lab researchers to conduct indel analysis from large NGS datasets. CRISPR-GRANT binaries are freely available for Linux (above Ubuntu 16.04), macOS (above High Sierra 10.13) and Windows (above Windows 7) at https://github.com/fuhuancheng/CRISPR-GRANT . CRISPR-GRANT source code is licensed under the GPLv3 license and free to download and use.

The Wnt/TCF7L1 transcriptional repressor axis drives primitive endoderm formation by antagonizing naive and formative pluripotency

Early during preimplantation development and in heterogeneous mouse embryonic stem cells (mESC) culture, pluripotent cells are specified towards either the primed epiblast or the primitive endoderm (PE) lineage. Canonical Wnt signaling is crucial for safeguarding naive pluripotency and embryo implantation, yet the role and relevance of canonical Wnt inhibition during early mammalian development remains unknown. Here, we demonstrate that transcriptional repression exerted by Wnt/TCF7L1 promotes PE differentiation of mESCs and in preimplantation inner cell mass. Time-series RNA sequencing and promoter occupancy data reveal that TCF7L1 binds and represses genes encoding essential naive pluripotency factors and indispensable regulators of the formative pluripotency program, including Otx2 and Lef1. Consequently, TCF7L1 promotes pluripotency exit and suppresses epiblast lineage formation, thereby driving cells into PE specification. Conversely, TCF7L1 is required for PE specification as deletion of Tcf7l1 abrogates PE differentiation without restraining epiblast priming. Taken together, our study underscores the importance of transcriptional Wnt inhibition in regulating lineage specification in ESCs and preimplantation embryo development as well as identifies TCF7L1 as key regulator of this process.

A tapt1 knock-out zebrafish line with aberrant lens development and impaired vision models human early-onset cataract

Bi-allelic mutations in the gene coding for human trans-membrane anterior-posterior transformation protein 1 (TAPT1) result in a broad phenotypic spectrum, ranging from syndromic disease with severe skeletal and congenital abnormalities to isolated early-onset cataract. We present here the first patient with a frameshift mutation in the TAPT1 gene, resulting in both bilateral early-onset cataract and skeletal abnormalities, in addition to several dysmorphic features, in this way further expanding the phenotypic spectrum associated with TAPT1 mutations. A tapt1a/tapt1b double knock-out (KO) zebrafish model generated by CRISPR/Cas9 gene editing revealed an early larval phenotype with eye malformations, loss of vision, increased photokinetics and hyperpigmentation, without visible skeletal involvement. Ultrastructural analysis of the eyes showed a smaller condensed lens, loss of integrity of the lens capsule with formation of a secondary lens and hyperplasia of the cells in the ganglion and inner plexiform layers of the retina. Transcriptomic analysis pointed to an impaired lens development with aberrant expression of many of the crystallin and other lens-specific genes. Furthermore, the phototransduction and visual perception pathways were found to be significantly disturbed. Differences in light perception are likely the cause of the increased dark photokinetics and generalized hyperpigmentation observed in this zebrafish model. In conclusion, this study validates TAPT1 as a new gene for early-onset cataract and sheds light on its ultrastructural and molecular characteristics.

Improved selection of zebrafish CRISPR editing by early next-generation sequencing based genotyping

Despite numerous prior attempts to improve knock-in (KI) efficiency, the introduction of precise base pair substitutions by the CRISPR-Cas9 technique in zebrafish remains challenging. In our efforts to generate KI zebrafish models of human CACNA1C mutations, we have tested the effect of several CRISPR determinants on KI efficiency across two sites in a single gene and developed a novel method for early selection to ameliorate KI efficiency. We identified optimal KI conditions for Cas9 protein and non-target asymmetric PAM-distal single stranded deoxynucleotide repair templates at both cacna1c sites. An effect of distance to the cut site on the KI efficiency was only observed for a single repair template conformation at one of the two sites. By combining minimally invasive early genotyping with the zebrafish embryo genotyper (ZEG) device and next-generation sequencing, we were able to obtain an almost 17-fold increase in somatic editing efficiency. The added benefit of the early selection procedure was particularly evident for alleles with lower somatic editing efficiencies. We further explored the potential of the ZEG selection procedure for the improvement of germline transmission by demonstrating germline transmission events in three groups of pre-selected embryos.

CRISPR genome editing using computational approaches: A survey

Clustered regularly interspaced short palindromic repeats (CRISPR)-based gene editing has been widely used in various cell types and organisms. To make genome editing with Clustered regularly interspaced short palindromic repeats far more precise and practical, we must concentrate on the design of optimal gRNA and the selection of appropriate Cas enzymes. Numerous computational tools have been created in recent years to help researchers design the best gRNA for Clustered regularly interspaced short palindromic repeats researches. There are two approaches for designing an appropriate gRNA sequence (which targets our desired sites with high precision): experimental and predicting-based approaches. It is essential to reduce off-target sites when designing an optimal gRNA. Here we review both traditional and machine learning-based approaches for designing an appropriate gRNA sequence and predicting off-target sites. In this review, we summarize the key characteristics of all available tools (as far as possible) and compare them together. Machine learning-based tools and web servers are believed to become the most effective and reliable methods for predicting on-target and off-target activities of Clustered regularly interspaced short palindromic repeats in the future. However, these predictions are not so precise now and the performance of these algorithms -especially deep learning one's-depends on the amount of data used during training phase. So, as more features are discovered and incorporated into these models, predictions become more in line with experimental observations. We must concentrate on the creation of ideal gRNA and the choice of suitable Cas enzymes in order to make genome editing with Clustered regularly interspaced short palindromic repeats far more accurate and feasible.

Evaluation of genome and base editing tools in maize protoplasts

INTRODUCTION

Despite its rapid worldwide adoption as an efficient mutagenesis tool, plant genome editing remains a labor-intensive process requiring often several months of in vitro culture to obtain mutant plantlets. To avoid a waste in time and money and to test, in only a few days, the efficiency of molecular constructs or novel Cas9 variants (clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9) prior to stable transformation, rapid analysis tools are helpful.

METHODS

To this end, a streamlined maize protoplast system for transient expression of CRISPR/Cas9 tools coupled to NGS (next generation sequencing) analysis and a novel bioinformatics pipeline was established.

RESULTS AND DISCUSSION

Mutation types found with high frequency in maize leaf protoplasts had a trend to be the ones observed after stable transformation of immature maize embryos. The protoplast system also allowed to conclude that modifications of the sgRNA (single guide RNA) scaffold leave little room for improvement, that relaxed PAM (protospacer adjacent motif) sites increase the choice of target sites for genome editing, albeit with decreased frequency, and that efficient base editing in maize could be achieved for certain but not all target sites. Phenotypic analysis of base edited mutant maize plants demonstrated that the introduction of a stop codon but not the mutation of a serine predicted to be phosphorylated in the bHLH (basic helix loop helix) transcription factor ZmICEa (INDUCER OF CBF EXPRESSIONa) caused abnormal stomata, pale leaves and eventual plant death two months after sowing.

Engraftment of Allotransplanted Tumor Cells in Adult rag2 Mutant Xenopus tropicalis

Simple Summary

While the mouse is without doubt the most studied animal for experimental cancer research, aquatic vertebrates such as zebrafish have also contributed to the field. More recently, thanks to the Nobel-prize winning technology of CRISPR/Cas mediated genomic engineering, the frog Xenopus tropicalis has emerged as an additional powerful model for studying human cancer. Via CRISPR-mediated genome editing, several models for different human cancers have been obtained in this animal. However, what has been lacking in Xenopus is the possibility to transplant tumor cells between different frogs. This is important to allow better characterization of the tumor cells and exploration of therapeutic opportunities. In this paper, we describe the generation of a genetic mutant in Xenopus tropicalis that has a compromised immune system, thereby allowing the grafting and expansion of tumors obtained in this species. In addition, an optimized protocol is provided for the irradiation of wild-type Xenopus frogs that subsequently are temporarily immunocompromised and during that period allow tumor engraftment. This work will expand the toolbox for modeling human cancer in Xenopus tropicalis, thereby further establishing it as a powerful experimental cancer model.

Abstract

Modeling human genetic diseases and cancer in lab animals has been greatly aided by the emergence of genetic engineering tools such as TALENs and CRISPR/Cas9. We have previously demonstrated the ease with which genetically engineered Xenopus models (GEXM) can be generated via injection of early embryos with Cas9 recombinant protein loaded with sgRNAs targeting single or multiple tumor suppressor genes. What has been lacking so far is the possibility to propagate and characterize the induced cancers via transplantation. Here, we describe the generation of a rag2 knockout line in Xenopus tropicalis that is deficient in functional T and B cells. This line was validated by means of allografting experiments with primary tp53^−/− and apc^+/−/tp53^−/− donor tumors. In addition, we optimized available protocols for the sub-lethal irradiation of wild-type X. tropicalis froglets. Irradiated animals also allowed the stable, albeit transient, engraftment of transplanted X. tropicalis tumor cells. The novel rag2^−/− line and the irradiated wild-type froglets will further expand the experimental toolbox in the diploid amphibian X. tropicalis and help to establish it as a versatile and relevant model for exploring human cancer.

TEAD4 regulates trophectoderm differentiation upstream of CDX2 in a GATA3-independent manner in the human preimplantation embryo

STUDY QUESTION

What is the role of transcriptional-enhanced associate (TEA) domain family member 4 (TEAD4) in trophectoderm (TE) differentiation during human embryo preimplantation development in comparison to mouse?

SUMMARY ANSWER

TEAD4 regulates TE lineage differentiation in the human preimplantation embryo acting upstream of caudal-type homeobox protein 2 (CDX2), but in contrast to the mouse in a GATA-binding protein 3 (GATA3)-independent manner.

WHAT IS KNOWN ALREADY

Tead4 is one of the earliest transcription factors expressed during mouse embryo preimplantation development and is required for the expression of TE-associated genes. Functional knock-out studies in mouse, inactivating Tead4 by site-specific recombination, have shown that Tead4-targeted embryos have compromised development and expression of the TE-specific Cdx2 and Gata3 is downregulated. Cdx2 and Gata3 act in parallel pathways downstream of Tead4 to induce successful TE differentiation. Downstream loss of Cdx2 expression, compromises TE differentiation and subsequent blastocoel formation and leads to the ectopic expression of inner cell mass (ICM) genes, including POU Class 5 homeobox 1 (Pou5f1) and SRY-box transcription factor (Sox2). Cdx2 is a more potent regulator of TE fate in mouse as loss of Cdx2 expression induces more severe phenotypes compared with loss of Gata3 expression. The role of TEAD4 and its downstream effectors during human preimplantation embryo development has not been investigated yet.

STUDY DESIGN, SIZE, DURATION

The clustered regularly interspaced short palindromic repeats-clustered regularly interspaced short palindromic repeats (CRISPR)-associated genes (CRISPR-Cas9) system was first introduced in pronuclei (PN)-stage mouse zygotes aiming to identify a guide RNA (gRNA), yielding high editing efficiency and effective disruption of the Tead4 locus. Three guides were tested (gRNA1-3), each time targeting a distinct region of Exon 2 of Tead4. The effects of targeting on developmental capacity were studied in Tead4-targeted embryos (n = 164-summarized data from gRNA1-3) and were compared with two control groups; sham-injected embryos (n = 26) and non-injected media-control embryos (n = 51). The editing efficiency was determined by next-generation sequencing (NGS). In total, n = 55 (summarized data from gRNA1-3) targeted mouse embryos were analysed by NGS. Immunofluorescence analysis to confirm successful targeting by gRNA1 was performed in Tead4-targeted embryos, and non-injected media-control embryos. The downregulation of secondary TE-associated markers Cdx2 and Gata3 was used as an indirect confirmation of successful Tead4-targeting (previously shown to be expressed downstream of Tead4). Additional groups of gRNA1 Tead4-targeted (n = 45) and media control (n = 36) embryos were cultured for an extended period of 8.5 days, to further assess the developmental capacity of the Tead4-targeted group to develop beyond implantation stages. Following the mouse investigation, human metaphase-II (MII) oocytes obtained by IVM were microinjected with gRNA-Cas9 during ICSI (n = 74) to target TEAD4 or used as media-control (n = 33). The editing efficiency was successfully assessed in n = 25 TEAD4-targeted human embryos. Finally, immunofluorescence analysis for TEAD4, CDX2, GATA3 and the ICM marker SOX2 was performed in TEAD4-targeted (n = 10) and non-injected media-control embryos (n = 29).

PARTICIPANTS/MATERIALS, SETTING, METHODS

A ribonucleoprotein complex consisting of a gRNA-Cas9 mixture, designed to target Exon 2 of Tead4/TEAD4, was microinjected in mouse PN stage zygotes or human IVM MII oocytes along with sperm. Generated embryos were cultured in vitro for 4 days in mouse or 6.5 days in human. In mouse, an additional group of Tead4-targeted and media-control embryos was cultured in vitro for an extended period of 8.5 days. Embryonic development and morphology were assessed daily, during culture in vitro of mouse and human embryos and was followed by a detailed scoring at late blastocyst stage. Targeting efficiency following gRNA-Cas9 introduction was assessed via immunostaining and NGS analysis.

MAIN RESULTS AND THE ROLE OF CHANCE

NGS analysis of the Tead4-targeted locus revealed very high editing efficiencies for all three guides, with 100% of the mouse embryos (55 out of 55) carrying genetic modifications resulting from CRISPR-Cas9 genome editing. More specifically, 65.22% (15 out 23) of the PN zygotes microinjected with gRNA1-Cas9, which exhibited the highest efficiency, carried exclusively mutated alleles. The developmental capacity of targeted embryos was significantly reduced (data from gRNA1), as 44.17% of the embryos arrested at the morula stage (2.5 days post coitum), coincident with the initiation of TE lineage differentiation, compared with 8.51% in control and 12.50% in sham control groups. High-quality blastocyst formation rates (Grade 3) were 8.97% in the gRNA1-targeted group, compared with 87.23% in the media-control and 87.50% in the sham group. Immunofluorescence analysis in targeted embryos confirmed downregulation of Tead4, Cdx2, and Gata3 expression, which resulted from successful targeting of the Tead4 locus. Tead4-targeted mouse embryos stained positive for the ICM markers Pou5f1 and Sox2, indicating that expression of ICM lineage markers is not affected. Tead4-targeted embryos were able to cavitate and form a blastocoel without being able to hatch. Extended embryo culture following zona pellucida removal, revealed that the targeted embryos can attach and form egg-cylinder-like structures in the absence of trophoblast giant cells. In human embryos, Exon 2 of TEAD4 was successfully targeted by CRISPR-Cas9 (n = 74). In total, 25 embryos from various developmental stages were analysed by NGS and 96.00% (24 out of 25) of the embryos carried genetic modifications because of gRNA-Cas9 editing. In the subgroup of the 24 edited embryos, 17 (70.83%) carried only mutant alleles and 11 out of these 17 (64.70%) carried exclusively frameshift mutations. Six out of 11 embryos reached the blastocyst stage. In contrast to mice, human-targeted embryos formed blastocysts at a rate (25.00%) that did not differ significantly from the control group (23.81%). However, blastocyst morphology and TE quality were significantly compromised following TEAD4-targeting, showing grade C TE scores, with TE containing very few cells. Immunofluorescence analysis of TEAD4-targeted embryos (n = 10) confirmed successful editing by the complete absence of TEAD4 and its downstream TE marker CDX2, but the embryos generated retained expression of GATA3, which is in contrast to what we have observed and has previously been reported in mouse. In this regard, our results indicate that GATA3 acts in parallel with TEAD4/CDX2 towards TE differentiation in human.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

CRISPR-Cas9 germline genome editing, in some cases, induces mosaic genotypes. These genotypes are a result of inefficient and delayed editing, and complicate the phenotypic analysis and developmental assessment of the injected embryos. We cannot exclude the possibility that the observed differences between mouse and human are the result of variable effects triggered by the culture conditions, which were however similar for both mouse and human embryos in this study. Furthermore, this study utilized human oocytes obtained by IVM, which may not fully recapitulate the developmental behaviour of in vivo matured oocytes.

WIDER IMPLICATIONS OF THE FINDINGS

Elucidation of the evolutionary conservation of molecular mechanisms that regulate the differentiation and formation of the trophoblast lineage can give us fundamental insights into early implantation failure, which accounts for ∼15% of human conceptions.

STUDY FUNDING/COMPETING INTEREST(S)

The research was funded by the FWO-Vlaanderen (Flemish fund for scientific research, Grant no. G051516N), and Hercules funding (FWO.HMZ.2016.00.02.01) and Ghent University (BOF.BAS.2018.0018.01). G.C. is supported by FWO-Vlaanderen (Flemish fund for scientific research, Grant no. 11L8822N). A.B. is supported by FWO-Vlaanderen (Flemish fund for scientific research, Grant no. 1298722 N). We further thank Ferring Pharmaceuticals (Aalst, Belgium) for their unrestricted educational grant. The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Computational tools and resources for CRISPR/Cas genome editing

The past decade has witnessed a rapid evolution in identifying more versatile clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) nucleases and their functional variants as well as in developing precise CRISPR/Cas-derived genome editors. The programmable and robust features of the genomic editors provide an effective RNA-guided platform for fundamental life science research and subsequent applications in diverse scenarios, including biomedical innovation and targeted crop improvement. One of the most essential principles is to guide alterations in genomic sequences or genes in the intended manner without undesired off-target impacts, which strongly depends on the efficiency and specificity of single guide RNA (sgRNA)-directed recognition of targeted DNA sequences. Recent advances in empirical scoring algorithms and machine learning models have facilitated sgRNA design and off-target prediction. In this review, we first briefly introduced the different features of CRISPR/Cas tools that should be taken into consideration to achieve specific purposes. Secondly, we focused on the computer-assisted tools and resources that are widely used in designing sgRNAs and analyzing CRISPR/Cas-induced on- and off-target mutations. Thirdly, we provide insights on the limitations of available computational tools that surely help researchers of this field for further optimization. Lastly, we suggested a simple but effective workflow for choosing and applying web-based resources and tools for CRISPR/Cas genome editing.

A Bioinformatic Workflow for InDel Analysis in the Wheat Multi-Copy α-Gliadin Gene Family Engineered with CRISPR/Cas9

The α-gliadins of wheat, along with other gluten components, are responsible for bread viscoelastic properties. However, they are also related to human pathologies as celiac disease or non-celiac wheat sensitivity. CRISPR/Cas was successfully used to knockout α-gliadin genes in bread and durum wheat, therefore, obtaining low gluten wheat lines. Nevertheless, the mutation analysis of these genes is complex as they present multiple and high homology copies arranged in tandem in A, B, and D subgenomes. In this work, we present a bioinformatic pipeline based on NGS amplicon sequencing for the analysis of insertions and deletions (InDels) in α-gliadin genes targeted with two single guides RNA (sgRNA). This approach allows the identification of mutated amplicons and the analysis of InDels through comparison to the most similar wild type parental sequence. TMM normalization was performed for inter-sample comparisons; being able to study the abundance of each InDel throughout generations and observe the effects of the segregation of Cas9 coding sequence in different lines. The usefulness of the workflow is relevant to identify possible genomic rearrangements such as large deletions due to Cas9 cleavage activity. This pipeline enables a fast characterization of mutations in multiple samples for a multi-copy gene family.

CRISPR-SID: Identifying EZH2 as a druggable target for desmoid tumors via in vivo dependency mapping

Cancer precision medicine implies identification of tumor-specific vulnerabilities associated with defined oncogenic pathways. Desmoid tumors are soft-tissue neoplasms strictly driven by Wnt signaling network hyperactivation. Despite this clearly defined genetic etiology and the strict and unique implication of the Wnt/β-catenin pathway, no specific molecular targets for these tumors have been identified. To address this caveat, we developed fast, efficient, and penetrant genetic Xenopus tropicalis desmoid tumor models to identify and characterize drug targets. We used multiplexed CRISPR/Cas9 genome editing in these models to simultaneously target a tumor suppressor gene (apc) and candidate dependency genes. Our methodology CRISPR/Cas9 selection-mediated identification of dependencies (CRISPR-SID) uses calculated deviations between experimentally observed gene editing outcomes and deep-learning-predicted double-strand break repair patterns to identify genes under negative selection during tumorigenesis. This revealed EZH2 and SUZ12, both encoding polycomb repressive complex 2 components, and the transcription factor CREB3L1 as genetic dependencies for desmoid tumors. In vivo EZH2 inhibition by Tazemetostat induced partial regression of established autochthonous tumors. In vitro models of patient desmoid tumor cells revealed a direct effect of Tazemetostat on Wnt pathway activity. CRISPR-SID represents a potent approach for in vivo mapping of tumor vulnerabilities and drug target identification.

Computational resources to define alleles and altered regulatory motifs at genomically edited candidate response elements

Unequivocal functional assessment of candidate genomic regulatory regions, such as transcriptional response elements, requires genetic alteration at their native chromosomal loci. Targeted DNA cleavage by Cas9 or other programmable nucleases enables analysis at virtually any genomic region, and diverse alleles generated by editing can be defined by deep sequencing for functional analysis. Interpretation of disrupted response elements, however, presents a special challenge, as these regions typically comprise clustered DNA binding motifs for multiple transcriptional regulatory factors (TFs); DNA sequence differences, natural or engineered, that affect binding by one TF can confer loss or gain of binding sites for other TFs. To address these and other analytical complexities, we created three computational tools that together integrate, in a single experiment, allele definition and TF binding motif evaluation for up to 9216 clones isolated, sequenced and propagated from Cas9-treated cell populations. We demonstrate 1) the capacity to functionally assess edited TF binding sites to query response element function, and 2) the efficacy and utility of these tools, by analyzing cell populations targeted by Cas9 for disruption of example glucocorticoid receptor (GR) binding motifs near FKBP5, a GR-regulated gene in the human adenocarcinoma cell line A549.

Lrp5 Mutant and Crispant Zebrafish Faithfully Model Human Osteoporosis, Establishing the Zebrafish as a Platform for CRISPR-Based Functional Screening of Osteoporosis Candidate Genes

Genomewide association studies (GWAS) have improved our understanding of the genetic architecture of common complex diseases such as osteoporosis. Nevertheless, to attribute functional skeletal contributions of candidate genes to osteoporosis-related traits, there is a need for efficient and cost-effective in vivo functional testing. This can be achieved through CRISPR-based reverse genetic screens, where phenotyping is traditionally performed in stable germline knockout (KO) mutants. Recently it was shown that first-generation (F0) mosaic mutant zebrafish (so-called crispants) recapitulate the phenotype of germline KOs. To demonstrate feasibility of functional validation of osteoporosis candidate genes through crispant screening, we compared a crispant to a stable KO zebrafish model for the lrp5 gene. In humans, recessive loss-of-function mutations in LRP5, a co-receptor in the Wnt signaling pathway, cause osteoporosis-pseudoglioma syndrome. In addition, several GWAS studies identified LRP5 as a major risk locus for osteoporosis-related phenotypes. In this study, we showed that early stage lrp5 KO larvae display decreased notochord mineralization and malformations of the head cartilage. Quantitative micro-computed tomography (micro-CT) scanning and mass-spectrometry element analysis of the adult skeleton revealed decreased vertebral bone volume and bone mineralization, hallmark features of osteoporosis. Furthermore, regenerating fin tissue displayed reduced Wnt signaling activity in lrp5 KO adults. We next compared lrp5 mutants with crispants. Next-generation sequencing analysis of adult crispant tissue revealed a mean out-of-frame mutation rate of 76%, resulting in strongly reduced levels of Lrp5 protein. These crispants generally showed a milder but nonetheless highly comparable skeletal phenotype and a similarly reduced Wnt pathway response compared with lrp5 KO mutants. In conclusion, we show through faithful modeling of LRP5-related primary osteoporosis that crispant screening in zebrafish is a promising approach for rapid functional screening of osteoporosis candidate genes. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Loss of zebrafish atp6v1e1b, encoding a subunit of vacuolar ATPase, recapitulates human ARCL type 2C syndrome and identifies multiple pathobiological signatures

The inability to maintain a strictly regulated endo(lyso)somal acidic pH through the proton-pumping action of the vacuolar-ATPases (v-ATPases) has been associated with various human diseases including heritable connective tissue disorders. Autosomal recessive (AR) cutis laxa (CL) type 2C syndrome is associated with genetic defects in the ATP6V1E1 gene and is characterized by skin wrinkles or loose redundant skin folds with pleiotropic systemic manifestations. The underlying pathological mechanisms leading to the clinical presentations remain largely unknown. Here, we show that loss of atp6v1e1b in zebrafish leads to early mortality, associated with craniofacial dysmorphisms, vascular anomalies, cardiac dysfunction, N-glycosylation defects, hypotonia, and epidermal structural defects. These features are reminiscent of the phenotypic manifestations in ARCL type 2C patients. Our data demonstrates that loss of atp6v1e1b alters endo(lyso)somal protein levels, and interferes with non-canonical v-ATPase pathways in vivo. In order to gain further insights into the processes affected by loss of atp6v1e1b, we performed an untargeted analysis of the transcriptome, metabolome, and lipidome in early atp6v1e1b-deficient larvae. We report multiple affected pathways including but not limited to oxidative phosphorylation, sphingolipid, fatty acid, and energy metabolism together with profound defects on mitochondrial respiration. Taken together, our results identify complex pathobiological effects due to loss of atp6v1e1b in vivo.

Cutis laxa syndromes are pleiotropic disorders of the connective tissue, characterized by skin redundancy and variable systemic manifestations. Cutis laxa syndromes are caused by pathogenic variants in genes encoding structural and regulatory components of the extracellular matrix or in genes encoding components of cellular trafficking, metabolism, and mitochondrial function. Pathogenic variants in genes coding for vacuolar-ATPases, a multisubunit complex responsible for the acidification of multiple intracellular vesicles, cause type 2 cutis laxa syndromes, a group of cutis laxa subtypes further characterized by neurological, skeletal, and rarely cardiopulmonary manifestations. To investigate the pathomechanisms of vacuolar-ATPase dysfunction, we generated zebrafish models that lack a crucial subunit of the vacuolar-ATPases. The mutant zebrafish models show morphological and functional features reminiscent of the phenotypic manifestations in cutis laxa patients carrying pathogenic variants in ATP6V1E1. In-depth analysis at multiple -omic levels identified biological signatures that indicate impairment of signaling pathways, lipid metabolism, and mitochondrial respiration. We anticipate that these data will contribute to a better understanding of the pathogenesis of cutis laxa syndromes and other disorders involving defective v-ATPase function, which may eventually improve patient treatment and management.

Bi-allelic premature truncating variants in LTBP1 cause cutis laxa syndrome

Latent transforming growth factor β (TGFβ)-binding proteins (LTBPs) are microfibril-associated proteins essential for anchoring TGFβ in the extracellular matrix (ECM) as well as for correct assembly of ECM components. Variants in LTBP2, LTBP3, and LTBP4 have been identified in several autosomal recessive Mendelian disorders with skeletal abnormalities with or without impaired development of elastin-rich tissues. Thus far, the human phenotype associated with LTBP1 deficiency has remained enigmatic. In this study, we report homozygous premature truncating LTBP1 variants in eight affected individuals from four unrelated consanguineous families. Affected individuals present with connective tissue features (cutis laxa and inguinal hernia), craniofacial dysmorphology, variable heart defects, and prominent skeletal features (craniosynostosis, short stature, brachydactyly, and syndactyly). In vitro studies on proband-derived dermal fibroblasts indicate distinct molecular mechanisms depending on the position of the variant in LTBP1. C-terminal variants lead to an altered LTBP1 loosely anchored in the microfibrillar network and cause increased ECM deposition in cultured fibroblasts associated with excessive TGFβ growth factor activation and signaling. In contrast, N-terminal truncation results in a loss of LTBP1 that does not alter TGFβ levels or ECM assembly. In vivo validation with two independent zebrafish lines carrying mutations in ltbp1 induce abnormal collagen fibrillogenesis in skin and intervertebral ligaments and ectopic bone formation on the vertebrae. In addition, one of the mutant zebrafish lines shows voluminous and hypo-mineralized vertebrae. Overall, our findings in humans and zebrafish show that LTBP1 function is crucial for skin and bone ECM assembly and homeostasis.

Comparative analysis of mouse and human preimplantation development following POU5F1 CRISPR/Cas9 targeting reveals interspecies differences

STUDY QUESTION

What is the role of POU class 5 homeobox 1 (POU5F1) in human preimplantation development and how does it compare with the mouse model?

SUMMARY ANSWER

POU5F1 is required for successful development of mouse and human embryos to the blastocyst stage as knockout embryos exhibited a significantly lower blastocyst formation rate, accompanied by lack of inner cell mass (ICM) formation.

WHAT IS KNOWN ALREADY

Clustered regularly interspaced short palindromic repeats-CRISPR associated genes (CRISPR-Cas9) has previously been used to examine the role of POU5F1 during human preimplantation development. The reported POU5F1-targeted blastocysts always retained POU5F1 expression in at least one cell, because of incomplete CRISPR-Cas9 editing. The question remains of whether the inability to obtain fully edited POU5F1-targeted blastocysts in human results from incomplete editing or the actual inability of these embryos to reach the blastocyst stage.

STUDY DESIGN, SIZE, DURATION

The efficiency of CRISPR-Cas9 to induce targeted gene mutations was first optimized in the mouse model. Two CRISPR-Cas9 delivery methods were compared in the B6D2F1 strain: S-phase injection (zygote stage) (n = 135) versus metaphase II-phase (M-phase) injection (oocyte stage) (n = 23). Four control groups were included: non-injected media-control zygotes (n = 43)/oocytes (n = 48); sham-injected zygotes (n = 45)/oocytes (n = 47); Cas9-protein injected zygotes (n = 23); and Cas9 protein and scrambled guide RNA (gRNA)-injected zygotes (n = 27). Immunofluorescence analysis was performed in Pou5f1-targeted zygotes (n = 37), media control zygotes (n = 19), and sham-injected zygotes (n = 15). To assess the capacity of Pou5f1-null embryos to develop further in vitro, additional groups of Pou5f1-targeted zygotes (n = 29) and media control zygotes (n = 30) were cultured to postimplantation stages (8.5 dpf). Aiming to identify differences in developmental capacity of Pou5f1-null embryos attributed to strain variation, zygotes from a second mouse strain-B6CBA (n = 52) were targeted. Overall, the optimized methodology was applied in human oocytes following IVM (metaphase II stage) (n = 101). The control group consisted of intracytoplasmically sperm injected (ICSI) IVM oocytes (n = 33). Immunofluorescence analysis was performed in human CRISPR-injected (n = 10) and media control (n = 9) human embryos.

PARTICIPANTS/MATERIALS, SETTING, METHODS

A gRNA-Cas9 protein mixture targeting exon 2 of Pou5f1/POU5F1 was microinjected in mouse oocytes/zygotes or human IVM oocytes. Reconstructed embryos were cultured for 4 days (mouse) or 6.5 days (human) in sequential culture media. An additional group of mouse-targeted zygotes was cultured to postimplantation stages. Embryonic development was assessed daily, with detailed scoring at late blastocyst stage. Genomic editing was assessed by immunofluorescence analysis and next-generation sequencing.

MAIN RESULTS AND THE ROLE OF CHANCE

Genomic analysis in mouse revealed very high editing efficiencies with 95% of the S-Phase and 100% of the M-Phase embryos containing genetic modifications, of which 89.47% in the S-Phase and 84.21% in the M-Phase group were fully edited. The developmental capacity was significantly compromised as only 46.88% embryos in the S-Phase and 19.05% in the M-Phase group reached the blastocyst stage, compared to 86.36% in control M-Phase and 90.24% in control S-Phase groups, respectively. Immunofluorescence analysis confirmed the loss of Pou5f1 expression and downregulation of the primitive marker SRY-Box transcription factor (Sox17). Our experiments confirmed the requirement of Pou5f1 expression for blastocyst development in the second B6CBA strain. Altogether, our data obtained in mouse reveal that Pou5f1 expression is essential for development to the blastocyst stage. M-Phase injection in human IVM oocytes (n = 101) similarly resulted in 88.37% of the POU5F1-targeted embryos being successfully edited. The developmental capacity of generated embryos was compromised from the eight-cell stage onwards. Only 4.55% of the microinjected embryos reached the late blastocyst stage and the embryos exhibited complete absence of ICM and an irregular trophectoderm cell layer. Loss of POU5F1 expression resulted in absence of SOX17 expression, as in mouse. Interestingly, genetic mosaicism was eliminated in a subset of targeted human embryos (9 out of 38), three of which developed into blastocysts.

LIMITATIONS, REASONS FOR CAUTION

One of the major hurdles of CRISPR-Cas9 germline genome editing is the occurrence of mosaicism, which may complicate phenotypic analysis and interpretation of developmental behavior of the injected embryos. Furthermore, in this study, spare IVM human oocytes were used, which may not recapitulate the developmental behavior of in vivo matured oocytes.

WIDER IMPLICATIONS OF THE FINDINGS

Comparison of developmental competency following CRISPR-Cas-mediated gene targeting in mouse and human may be influenced by the selected mouse strain. Gene targeting by CRISPR-Cas9 is subject to variable targeting efficiencies. Therefore, striving to reduce mosaicism can provide novel molecular insights into mouse and human embryogenesis.

STUDY FUNDING/COMPETING INTEREST(S)

The research was funded by the Ghent University Hospital and Ghent University and supported by the FWO-Vlaanderen (Flemish fund for scientific research, Grant no. G051516N), and Hercules funding (FWO.HMZ.2016.00.02.01). The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

CRISPECTOR provides accurate estimation of genome editing translocation and off-target activity from comparative NGS data

Controlling off-target editing activity is one of the central challenges in making CRISPR technology accurate and applicable in medical practice. Current algorithms for analyzing off-target activity do not provide statistical quantification, are not sufficiently sensitive in separating signal from noise in experiments with low editing rates, and do not address the detection of translocations. Here we present CRISPECTOR, a software tool that supports the detection and quantification of on- and off-target genome-editing activity from NGS data using paired treatment/control CRISPR experiments. In particular, CRISPECTOR facilitates the statistical analysis of NGS data from multiplex-PCR comparative experiments to detect and quantify adverse translocation events. We validate the observed results and show independent evidence of the occurrence of translocations in human cell lines, after genome editing. Our methodology is based on a statistical model comparison approach leading to better false-negative rates in sites with weak yet significant off-target activity.

In silico Method in CRISPR/Cas System: An Expedite and Powerful Booster

The CRISPR/Cas system has stood in the center of attention in the last few years as a revolutionary gene editing tool with a wide application to investigate gene functions. However, the labor-intensive workflow requires a sophisticated pre-experimental and post-experimental analysis, thus becoming one of the hindrances for the further popularization of practical applications. Recently, the increasing emergence and advancement of the in silico methods play a formidable role to support and boost experimental work. However, various tools based on distinctive design principles and frameworks harbor unique characteristics that are likely to confuse users about how to choose the most appropriate one for their purpose. In this review, we will present a comprehensive overview and comparisons on the in silico methods from the aspects of CRISPR/Cas system identification, guide RNA design, and post-experimental assistance. Furthermore, we establish the hypotheses in light of the new trends around the technical optimization and hope to provide significant clues for future tools development.

The Diversity of Genetic Outcomes from CRISPR/Cas Gene Editing is Regulated by the Length of the Symmetrical Donor DNA Template

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas gene editing systems have enabled molecular geneticists to manipulate prokaryotic and eukaryotic genomes with greater efficiency and precision. CRISPR/Cas provides adaptive immunity in bacterial cells by degrading invading viral genomes. By democratizing this activity into human cells, it is possible to knock out specific genes to disable their function and repair errors. The latter of these activities requires the participation of a single-stranded donor DNA template that provides the genetic information to execute correction in a process referred to as homology directed repair (HDR). Here, we utilized an established cell-free extract system to determine the influence that the donor DNA template length has on the diversity of products from CRISPR-directed gene editing. This model system enables us to view all outcomes of this reaction and reveals that donor template length can influence the efficiency of the reaction and the categories of error-prone products that accompany it. A careful measurement of the products revealed a category of error-prone events that contained the corrected template along with insertions and deletions (indels). Our data provides foundational information for those whose aim is to translate CRISPR/Cas from bench to bedside.

Maximizing CRISPR/Cas9 phenotype penetrance applying predictive modeling of editing outcomes in Xenopus and zebrafish embryos

CRISPR/Cas9 genome editing has revolutionized functional genomics in vertebrates. However, CRISPR/Cas9 edited F0 animals too often demonstrate variable phenotypic penetrance due to the mosaic nature of editing outcomes after double strand break (DSB) repair. Even with high efficiency levels of genome editing, phenotypes may be obscured by proportional presence of in-frame mutations that still produce functional protein. Recently, studies in cell culture systems have shown that the nature of CRISPR/Cas9-mediated mutations can be dependent on local sequence context and can be predicted by computational methods. Here, we demonstrate that similar approaches can be used to forecast CRISPR/Cas9 gene editing outcomes in Xenopus tropicalis, Xenopus laevis, and zebrafish. We show that a publicly available neural network previously trained in mouse embryonic stem cell cultures (InDelphi-mESC) is able to accurately predict CRISPR/Cas9 gene editing outcomes in early vertebrate embryos. Our observations can have direct implications for experiment design, allowing the selection of guide RNAs with predicted repair outcome signatures enriched towards frameshift mutations, allowing maximization of CRISPR/Cas9 phenotype penetrance in the F0 generation.

Technologies and Computational Analysis Strategies for CRISPR Applications

The CRISPR-Cas system offers a programmable platform for eukaryotic genome and epigenome editing. The ability to perform targeted genetic and epigenetic perturbations enables researchers to perform a variety of tasks, ranging from investigating questions in basic biology to potentially developing novel therapeutics for the treatment of disease. While CRISPR systems have been engineered to target DNA and RNA with increased precision, efficiency, and flexibility, assays to identify off-target editing are becoming more comprehensive and sensitive. Furthermore, techniques to perform high-throughput genome and epigenome editing can be paired with a variety of readouts and are uncovering important cellular functions and mechanisms. These technological advances drive and are driven by accompanying computational approaches. Here, we briefly present available CRISPR technologies and review key computational advances and considerations for various CRISPR applications. In particular, we focus on the analysis of on- and off-target editing and CRISPR pooled screen data.

Design and analysis of CRISPR-Cas experiments

A large and ever-expanding set of CRISPR-Cas systems now enables the rapid and flexible manipulation of genomes in both targeted and large-scale experiments. Numerous software tools and analytical methods have been developed for the design and analysis of CRISPR-Cas experiments, including resources to design optimal guide RNAs for various modes of manipulation and to analyze the results of such experiments. A major recent focus has been the development of comprehensive tools for use on data from large-scale CRISPR-based genetic screens. As this field continues to progress, a clear ongoing challenge is not only to innovate, but to actively maintain and improve existing tools so that researchers across disciplines can rely on a stable set of excellent computational resources for CRISPR-Cas experiments.

gRNA Sequence Heterology Tolerance Catalyzed by CRISPR/Cas in an In Vitro Homology-Directed Repair Reaction

CRISPR and associated Cas nucleases are genetic engineering tools revolutionizing innovative approaches to cancer and inherited diseases. CRISPR-directed gene editing relies heavily on proper DNA sequence alignment between the guide RNA (gRNA)/CRISPR complex and its genomic target. Accurate hybridization of complementary DNA initiates gene editing in human cells, but inherent gRNA sequence variation that could influence the gene editing reaction has been clearly established among diverse genetic populations. As this technology advances toward clinical implementation, it will be essential to assess what degree of gRNA variation generates unwanted and erroneous CRISPR activity. With the use of a system in which a cell-free extract catalyzes nonhomologous end joining (NHEJ) and homology-directed repair (HDR), it is possible to observe a more representative population of all forms of gene editing outcomes. In this manuscript, we demonstrate CRISPR/Cas complexation at heterologous binding sites that facilitate precise and error-prone HDR. The tolerance of mispairing between the gRNA and target site of the DNA to enable HDR is surprisingly high and greatly influenced by polarity of the donor DNA strand in the reaction. These results suggest that some collateral genomic activity could occur at unintended sites in CRISPR-directed gene editing in human cells.

A simple and practical workflow for genotyping of CRISPR-Cas9-based knockout phenotypes using multiplexed amplicon sequencing

Founder animals carrying high proportions of somatic mutation induced by CRISPR-Cas9 enable a rapid and scalable strategy for the functional screening of numerous target genes in vivo. In this functional screening, genotyping using pooled amplicons with next-generation sequencing is the most suitable approach for large-scale management of multiple samples and accurate evaluation of the efficiency of Cas9-induced somatic mutations at target sites. Here, we present a simple workflow for genotyping of multiple CRISPR-Cas9-based knockout founders by pooled amplicon sequencing. Using custom barcoded primers, pooled amplicons from multiple individuals can be run in a single-indexed library on the Illumina MiSeq platform. Additionally, a user-friendly web tool, CLiCKAR, is available to simultaneously perform demultiplexing of pooled sequence data and evaluation of somatic mutation in each phenotype. CLiCKAR provides users with practical reports regarding the positions of insertions/deletions, as well as the frameshift ratio and tables containing mutation sequences, and read counts of each phenotype, with just a few clicks by the implementation of demultiplexing for pooled sample data and calculation of the frameshift ratio. This genotyping workflow can be harnessed to evaluate genotype-phenotype correlations in CRISPR-Cas9-based loss-of-function screening of numerous target genes in various organisms.

Investigation of the role of VHL-HIF signaling in DNA repair and apoptosis in zebrafish

pVHL is a tumor suppressor. The lack of its function leads to various tumors, among which ccRCC (clear cell renal cell carcinoma) has the most serious outcome due to its resistance to chemotherapies and radiotherapies. Although HIF promotes the progression of ccRCC, the precise mechanism by which the loss of VHL leads to tumor initiation remains unclear. We exploited two zebrafish vhl mutants, vhl and vll, and Tg (phd3:: EGFP)i144 fish to identify crucial functions of Vhl in tumor initiation. Through the mutant analysis, we found that the role of pVHL in DNA repair is conserved in zebrafish Vll. Interestingly, we also discovered that Hif activation strongly suppressed genotoxic stress induced DNA repair defects and apoptosis in vll and brca2 mutants and in embryos lacking ATM activity. These results suggest the potential of HIF as a clinical modulator that can protect cells from accumulating DNA damage and apoptosis which can lead to cancers and neurodegenerative disorders.

RBL1 (p107) functions as tumor suppressor in glioblastoma and small-cell pancreatic neuroendocrine carcinoma in Xenopus tropicalis

Alterations of the retinoblastoma and/or the p53 signaling network are associated with specific cancers such as high-grade astrocytoma/glioblastoma, small-cell lung cancer (SCLC), choroid plexus tumors, and small-cell pancreatic neuroendocrine carcinoma (SC-PaNEC). However, the intricate functional redundancy between RB1 and the related pocket proteins RBL1/p107 and RBL2/p130 in suppressing tumorigenesis remains poorly understood. Here we performed lineage-restricted parallel inactivation of rb1 and rbl1 by multiplex CRISPR/Cas9 genome editing in the true diploid Xenopus tropicalis to gain insight into this in vivo redundancy. We show that while rb1 inactivation is sufficient to induce choroid plexus papilloma, combined rb1 and rbl1 inactivation is required and sufficient to drive SC-PaNEC, retinoblastoma and astrocytoma. Further, using a novel Li-Fraumeni syndrome-mimicking tp53 mutant X. tropicalis line, we demonstrate increased malignancy of rb1/rbl1-mutant glioma towards glioblastoma upon concomitant inactivation of tp53. Interestingly, although clinical SC-PaNEC samples are characterized by abnormal p53 expression or localization, in the current experimental models, the tp53 status had little effect on the establishment and growth of SC-PaNEC, but may rather be essential for maintaining chromosomal stability. SCLC was only rarely observed in our experimental setup, indicating requirement of additional or alternative oncogenic insults. In conclusion, we used CRISPR/Cas9 to delineate the tumor suppressor properties of Rbl1, generating new insights in the functional redundancy within the retinoblastoma protein family in suppressing neuroendocrine pancreatic cancer and glioma/glioblastoma.

Zebrafish: A Resourceful Vertebrate Model to Investigate Skeletal Disorders

Animal models are essential tools for addressing fundamental scientific questions about skeletal diseases and for the development of new therapeutic approaches. Traditionally, mice have been the most common model organism in biomedical research, but their use is hampered by several limitations including complex generation, demanding investigation of early developmental stages, regulatory restrictions on breeding, and high maintenance cost. The zebrafish has been used as an efficient alternative vertebrate model for the study of human skeletal diseases, thanks to its easy genetic manipulation, high fecundity, external fertilization, transparency of rapidly developing embryos, and low maintenance cost. Furthermore, zebrafish share similar skeletal cells and ossification types with mammals. In the last decades, the use of both forward and new reverse genetics techniques has resulted in the generation of many mutant lines carrying skeletal phenotypes associated with human diseases. In addition, transgenic lines expressing fluorescent proteins under bone cell- or pathway- specific promoters enable in vivo imaging of differentiation and signaling at the cellular level. Despite the small size of the zebrafish, many traditional techniques for skeletal phenotyping, such as x-ray and microCT imaging and histological approaches, can be applied using the appropriate equipment and custom protocols. The ability of adult zebrafish to remodel skeletal tissues can be exploited as a unique tool to investigate bone formation and repair. Finally, the permeability of embryos to chemicals dissolved in water, together with the availability of large numbers of small-sized animals makes zebrafish a perfect model for high-throughput bone anabolic drug screening. This review aims to discuss the techniques that make zebrafish a powerful model to investigate the molecular and physiological basis of skeletal disorders.

A qPCR method for genome editing efficiency determination and single-cell clone screening in human cells

CRISPR/Cas9 technology has been widely used for targeted genome modification both in vivo and in vitro. However, an effective method for evaluating genome editing efficiency and screening single-cell clones for desired modification is still lacking. Here, we developed this real time PCR method based on the sensitivity of Taq DNA polymerase to nucleotide mismatch at primer 3' end during initiating DNA replication. Applications to CRISPR gRNAs targeting EMX1, DYRK1A and HOXB13 genes in Lenti-X 293 T cells exhibited comprehensive advantages. Just in one-round qPCR analysis using genomic DNA from cells underwent CRISPR/Cas9 or BE4 treatments, the genome editing efficiency could be determined accurately and quickly, for indel, HDR as well as base editing. When applied to single-cell clone screening, the genotype of each cell colony could also be determined accurately. This method defined a rigorous and practical way in quantify genome editing events.

Homozygous Null TBX4 Mutations Lead to Posterior Amelia with Pelvic and Pulmonary Hypoplasia

The development of hindlimbs in tetrapod species relies specifically on the transcription factor TBX4. In humans, heterozygous loss-of-function TBX4 mutations cause dominant small patella syndrome (SPS) due to haploinsufficiency. Here, we characterize a striking clinical entity in four fetuses with complete posterior amelia with pelvis and pulmonary hypoplasia (PAPPA). Through exome sequencing, we find that PAPPA syndrome is caused by homozygous TBX4 inactivating mutations during embryogenesis in humans. In two consanguineous couples, we uncover distinct germline TBX4 coding mutations, p.Tyr113^∗ and p.Tyr127Asn, that segregated with SPS in heterozygous parents and with posterior amelia with pelvis and pulmonary hypoplasia syndrome (PAPPAS) in one available homozygous fetus. A complete absence of TBX4 transcripts in this proband with biallelic p.Tyr113^∗ stop-gain mutations revealed nonsense-mediated decay of the endogenous mRNA. CRISPR/Cas9-mediated TBX4 deletion in Xenopus embryos confirmed its restricted role during leg development. We conclude that SPS and PAPPAS are allelic diseases of TBX4 deficiency and that TBX4 is an essential transcription factor for organogenesis of the lungs, pelvis, and hindlimbs in humans.

Optimized knock-in of point mutations in zebrafish using CRISPR/Cas9

We have optimized point mutation knock-ins into zebrafish genomic sites using clustered regularly interspaced palindromic repeats (CRISPR)/Cas9 reagents and single-stranded oligodeoxynucleotides. The efficiency of knock-ins was assessed by a novel application of allele-specific polymerase chain reaction and confirmed by high-throughput sequencing. Anti-sense asymmetric oligo design was found to be the most successful optimization strategy. However, cut site proximity to the mutation and phosphorothioate oligo modifications also greatly improved knock-in efficiency. A previously unrecognized risk of off-target trans knock-ins was identified that we obviated through the development of a workflow for correct knock-in detection. Together these strategies greatly facilitate the study of human genetic diseases in zebrafish, with additional applicability to enhance CRISPR-based approaches in other animal model systems.

Xenopus tropicalis: Joining the Armada in the Fight Against Blood Cancer

Aquatic vertebrate organisms such as zebrafish have been used for over a decade to model different types of human cancer, including hematologic malignancies. However, the introduction of gene editing techniques such as CRISPR/Cas9 and TALEN, have now opened the road for other organisms featuring large externally developing embryos that are easily accessible. Thanks to its unique diploid genome that shows a high degree of synteny to the human, combined with its relatively short live cycle, Xenopus tropicalis has now emerged as an additional powerful aquatic model for studying human disease genes. Genome editing techniques are very simple and extremely efficient, permitting the fast and cheap generation of genetic models for human disease. Mosaic disruption of tumor suppressor genes allows the generation of highly penetrant and low latency cancer models. While models for solid human tumors have been recently generated, genetic models for hematologic malignancies are currently lacking for Xenopus. Here we describe our experimental pipeline, based on mosaic genome editing by CRISPR/Cas9, to generate innovative and high-performing leukemia models in X. tropicalis. These add to the existing models in zebrafish and will extend the experimental platform available in aquatic vertebrate organisms to contribute to the field of hematologic malignancies. This will extend our knowledge in the etiology of this cancer and assist the identification of molecular targets for therapeutic intervention.

CRISPR/Cas9-Mediated Knockout of Rb1 in Xenopus tropicalis

At this time, no molecular targeted therapies exist for treatment of retinoblastoma. This can be, in part, attributed to the lack of animal models that allow for both rapid identification of novel therapeutic targets and hypothesis driven drug testing. Within this scope, we have recently reported the first genuine genetic nonmammalian retinoblastoma cancer model within the aquatic model organism Xenopus tropicalis (Naert et al., Sci Rep 6: 35263, 2016). Here we describe the methods to generate rb1 mosaic mutant Xenopus tropicalis by employing the CRISPR/Cas9 technology. In depth, we discuss short guide RNA (sgRNA) design parameters, generation, quality control, quantification, and delivery followed by several methods for assessing genome editing efficiencies. As such the reader should be capable, by minor changes to the methods described here, to (co-) target rb1 or any one or multiple gene(s) within the Xenopus tropicalis genome by multiplex CRISPR/Cas9 methodology.

Tracking CRISPR's Footprints

The programmable clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9) and CRISPR-Cas9-derived gene editing and manipulation tools have revolutionized biomedical research over the past few years. One important category of assisting technologies in CRISPR gene editing is methods used for detecting and quantifying indels (deletions or insertions). These indels are caused by the repair of CRISPR-Cas9-introduced DNA double-stranded breaks (DBSs), known as CRISPR's DNA cleavage footprints. In addition, CRISPR-Cas9 can also leave footprints to the DNA without introducing DSBs, known as CRISPR's DNA-binding footprints. The indel tracking methods have contributed greatly to the improvement of CRISPR-Cas9 activity and specificity. Here, we review and discuss strategies developed over that past few years to track the CRISPR's footprints, their advantages, and limitations.

Hi-TOM: a platform for high-throughput tracking of mutations induced by CRISPR/Cas systems

The CRISPR/Cas system has been extensively applied to make precise genetic modifications in various organisms. Despite its importance and widespread use, large-scale mutation screening remains time-consuming, labour-intensive and costly. Here, we developed Hi-TOM (available at https://doi.org/www.hi-tom.net/hi-tom/ ), an online tool to track the mutations with precise percentage for multiple samples and multiple target sites. We also described a corresponding next-generation sequencing (NGS) library construction strategy by fixing the bridge sequences and barcoding primers. Analysis of the samples from rice, hexaploid wheat and human cells reveals that the Hi-TOM tool has high reliability and sensitivity in tracking various mutations, especially complex chimeric mutations frequently induced by genome editing. Hi-TOM does not require special design of barcode primers, cumbersome parameter configuration or additional data analysis. Thus, the streamlined NGS library construction and comprehensive result output make Hi-TOM particularly suitable for high-throughput identification of all types of mutations induced by CRISPR/Cas systems.

CRISPR/Cas9-mediated homology-directed repair by ssODNs in zebrafish induces complex mutational patterns resulting from genomic integration of repair-template fragments

Targeted genome editing by CRISPR/Cas9 is extremely well fitted to generate gene disruptions, although precise sequence replacement by CRISPR/Cas9-mediated homology-directed repair (HDR) suffers from low efficiency, impeding its use for high-throughput knock-in disease modeling. In this study, we used next-generation sequencing (NGS) analysis to determine the efficiency and reliability of CRISPR/Cas9-mediated HDR using several types of single-stranded oligodeoxynucleotide (ssODN) repair templates for the introduction of disease-relevant point mutations in the zebrafish genome. Our results suggest that HDR rates are strongly determined by repair-template composition, with the most influential factor being homology-arm length. However, we found that repair using ssODNs does not only lead to precise sequence replacement but also induces integration of repair-template fragments at the Cas9 cut site. We observed that error-free repair occurs at a relatively constant rate of 1-4% when using different repair templates, which was sufficient for transmission of point mutations to the F1 generation. On the other hand, erroneous repair mainly accounts for the variability in repair rate between the different repair templates. To further improve error-free HDR rates, elucidating the mechanism behind this erroneous repair is essential. We show that the error-prone nature of ssODN-mediated repair, believed to act via synthesis-dependent strand annealing (SDSA), is most likely due to DNA synthesis errors. In conclusion, caution is warranted when using ssODNs for the generation of knock-in models or for therapeutic applications. We recommend the application of in-depth NGS analysis to examine both the efficiency and error-free nature of HDR events.

This article has an associated First Person interview with the first author of the paper.

Summary: NGS-based analysis reveals that CRISPR/Cas9-induced double-strand-break repair using single-stranded repair templates is error prone in zebrafish, resulting in complex patterns of integrated repair-template fragments.

Functional testing of a human PBX3 variant in zebrafish reveals a potential modifier role in congenital heart defects

Whole-genome and exome sequencing efforts are increasingly identifying candidate genetic variants associated with human disease. However, predicting and testing the pathogenicity of a genetic variant remains challenging. Genome editing allows for the rigorous functional testing of human genetic variants in animal models. Congenital heart defects (CHDs) are a prominent example of a human disorder with complex genetics. An inherited sequence variant in the human PBX3 gene (PBX3 p.A136V) has previously been shown to be enriched in a CHD patient cohort, indicating that the PBX3 p.A136V variant could be a modifier allele for CHDs. Pbx genes encode three-amino-acid loop extension (TALE)-class homeodomain-containing DNA-binding proteins with diverse roles in development and disease, and are required for heart development in mouse and zebrafish. Here, we used CRISPR-Cas9 genome editing to directly test whether this Pbx gene variant acts as a genetic modifier in zebrafish heart development. We used a single-stranded oligodeoxynucleotide to precisely introduce the human PBX3 p.A136V variant in the homologous zebrafish pbx4 gene (pbx4 p.A131V). We observed that zebrafish that are homozygous for pbx4 p.A131V are viable as adults. However, the pbx4 p.A131V variant enhances the embryonic cardiac morphogenesis phenotype caused by loss of the known cardiac specification factor, Hand2. Our study is the first example of using precision genome editing in zebrafish to demonstrate a function for a human disease-associated single nucleotide variant of unknown significance. Our work underscores the importance of testing the roles of inherited variants, not just de novo variants, as genetic modifiers of CHDs. Our study provides a novel approach toward advancing our understanding of the complex genetics of CHDs.

Optimized knock-in of point mutations in zebrafish using CRISPR/Cas9

We have optimized point mutation knock-ins into zebrafish genomic sites using clustered regularly interspaced palindromic repeats (CRISPR)/Cas9 reagents and single-stranded oligodeoxynucleotides. The efficiency of knock-ins was assessed by a novel application of allele-specific polymerase chain reaction and confirmed by high-throughput sequencing. Anti-sense asymmetric oligo design was found to be the most successful optimization strategy. However, cut site proximity to the mutation and phosphorothioate oligo modifications also greatly improved knock-in efficiency. A previously unrecognized risk of off-target trans knock-ins was identified that we obviated through the development of a workflow for correct knock-in detection. Together these strategies greatly facilitate the study of human genetic diseases in zebrafish, with additional applicability to enhance CRISPR-based approaches in other animal model systems.

CRISPR-DAV: CRISPR NGS data analysis and visualization pipeline

Summary

The simplicity and precision of CRISPR/Cas9 system has brought in a new era of gene editing. Screening for desired clones with CRISPR-mediated genomic edits in a large number of samples is made possible by next generation sequencing (NGS) due to its multiplexing. Here we present CRISPR-DAV (CRISPR Data Analysis and Visualization) pipeline to analyze the CRISPR NGS data in a high throughput manner. In the pipeline, Burrows-Wheeler Aligner and Assembly Based ReAlignment are used for small and large indel detection, and results are presented in a comprehensive set of charts and interactive alignment view.

Availability and implementation

CRISPR-DAV is available at GitHub and Docker Hub repositories: https://github.com/pinetree1/crispr-dav.git and https://hub.docker.com/r/pinetree1/crispr-dav/.

Contact

xuning.wang@bms.com.