INDEL detection, the 'Achilles heel' of precise genome editing: a survey of methods for accurate profiling of gene editing induced indels.

Advances in genome editing technologies have enabled manipulation of genomes at the single base level. These technologies are based on programmable nucleases (PNs) that include meganucleases, zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated 9 (Cas9) nucleases and have given researchers the ability to delete, insert or replace genomic DNA in cells, tissues and whole organisms. The great flexibility in re-designing the genomic target specificity of PNs has vastly expanded the scope of gene editing applications in life science, and shows great promise for development of the next generation gene therapies. PN technologies share the principle of inducing a DNA double-strand break (DSB) at a user-specified site in the genome, followed by cellular repair of the induced DSB. PN-elicited DSBs are mainly repaired by the non-homologous end joining (NHEJ) and the microhomology-mediated end joining (MMEJ) pathways, which can elicit a variety of small insertion or deletion (indel) mutations. If indels are elicited in a protein coding sequence and shift the reading frame, targeted gene knock out (KO) can readily be achieved using either of the available PNs. Despite the ease by which gene inactivation in principle can be achieved, in practice, successful KO is not only determined by the efficiency of NHEJ and MMEJ repair; it also depends on the design and properties of the PN utilized, delivery format chosen, the preferred indel repair outcomes at the targeted site, the chromatin state of the target site and the relative activities of the repair pathways in the edited cells. These variables preclude accurate prediction of the nature and frequency of PN induced indels. A key step of any gene KO experiment therefore becomes the detection, characterization and quantification of the indel(s) induced at the targeted genomic site in cells, tissues or whole organisms. In this survey, we briefly review naturally occurring indels and their detection. Next, we review the methods that have been developed for detection of PN-induced indels. We briefly outline the experimental steps and describe the pros and cons of the various methods to help users decide a suitable method for their editing application. We highlight recent advances that enable accurate and sensitive quantification of indel events in cells regardless of their genome complexity, turning a complex pool of different indel events into informative indel profiles. Finally, we review what has been learned about PN-elicited indel formation through the use of the new methods and how this insight is helping to further advance the genome editing field.

Low incidence of SNVs and indels in trio genomes of Cas9-mediated multiplex edited sheep.

BACKGROUND: The simplicity of the CRISPR/Cas9 system has enabled its widespread applications in generating animal models, functional genomic screening and in treating genetic and infectious diseases. However, unintended mutations produced by off-target CRISPR/Cas9 nuclease activity may lead to negative consequences. Especially, a very recent study found that gene editing can introduce hundreds of unintended mutations into the genome, and have attracted wide attention. RESULTS: To address the off-target concerns, urgent characterization of the CRISPR/Cas9-mediated off-target mutagenesis is highly anticipated. Here we took advantage of our previously generated gene-edited sheep and performed family trio-based whole genome sequencing which is capable of discriminating variants in the edited progenies that are inherited, naturally generated, or induced by genetic modification. Three family trios were re-sequenced at a high average depth of genomic coverage (~ 25.8×). After developing a pipeline to comprehensively analyze the sequence data for de novo single nucleotide variants, indels and structural variations from the genome; we only found a single unintended event in the form of a 2.4 kb inversion induced by site-specific double-strand breaks between two sgRNA targeting sites at the MSTN locus with a low incidence. CONCLUSIONS: We provide the first report on the fidelity of CRISPR-based modification for sheep genomes targeted simultaneously for gene breaks at three coding sequence locations. The trio-based sequencing approach revealed almost negligible off-target modifications, providing timely evidences of the safe application of genome editing in vivo with CRISPR/Cas9.

Generation of gene-edited sheep with a defined Booroola fecundity gene (FecBB) mutation in bone morphogenetic protein receptor type 1B (BMPR1B) via clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) 9.

Since its emergence, the clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR-associated (Cas) 9 system has been increasingly used to generate animals for economically important traits. However, most CRISPR/Cas9 applications have been focused on non-homologous end joining, which results in base deletions and insertions, leading to a functional knockout of the targeted gene. The Booroola fecundity gene (FecBB) mutation (p.Q249R) in bone morphogenetic protein receptor type 1B (BMPR1B) has been demonstrated to exert a profound effect on fecundity in many breeds of sheep. In the present study, we successfully obtained lambs with defined point mutations resulting in a p.249Q>R substitution through the coinjection of Cas9 mRNA, a single guide RNA and single-stranded DNA oligonucleotides into Tan sheep zygotes. In the newborn lambs, the observed efficiency of the single nucleotide exchange was as high as 23.8%. We believe that our findings will contribute to improved reproduction traits in sheep, as well as to the generation of defined point mutations in other large animals.

Efficient generation of goats with defined point mutation (I397V) in GDF9 through CRISPR/Cas9.

The recent emergence of the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) 9 system has attracted significant attention for its potential to improve traits of agricultural importance. However, most applications in livestock species to date have depended on aberrant DNA repair to generate frameshifting indels. Whether this genomic engineering technique involving homology-dependent repair (HDR) can be used to introduce defined point mutations has been less explored. Previously, we reported a G→A point mutation (g.231A>G, p.Val397Ile) in the growth differentiation factor 9 (GDF9) gene that has a large effect on the litter size of cashmere goats. In the present study we report that by co-injecting synthesised RNAs and single-stranded oligo deoxynucleotide (ssODN) donor sequences into goat zygotes, we successfully introduced defined point mutations resulting in single amino acid substitutions in the proteins as expected. The efficiency of this precise single-nucleotide substitution in newborn kids was as high as 24% (4/17), indicating that ssODN-directed HDR via zygote injection is efficient at introducing point mutations in the goat genome. The findings of the present study further highlight the complex genome modifications facilitated by the CRISPR/Cas9 system, which is able to introduce defined point mutations. This represents a significant development for the improvement of reproduction traits in goats, as well as for validating the roles of specific nucleotides in functional genetic elements in large animals.

Tß4-overexpression based on the piggyBac transposon system in cashmere goats alters hair fiber characteristics.

Increasing cashmere yield is one of the vital aims of cashmere goats breeding. Compared to traditional breeding methods, transgenic technology is more efficient and the piggyBac (PB) transposon system has been widely applied to generate transgenic animals. For the present study, donor fibroblasts were stably transfected via a PB donor vector containing the coding sequence of cashmere goat thymosin beta-4 (Tß4) and driven by a hair follicle-specific promoter, the keratin-associated protein 6.1 (KAP6.1) promoter. To obtain genetically modified cells as nuclear donors, we co-transfected donor vectors into fetal fibroblasts of cashmere goats. Five transgenic cashmere goats were generated following somatic cell nuclear transfer (SCNT). Via determination of the copy numbers and integration sites, the Tß4 gene was successfully inserted into the goat genome. Histological examination of skin tissue revealed that Tß4-overexpressing, transgenic goats had a higher secondary to primary hair follicle (S/P) ratio compared to wild type goats. This indicates that Tß4-overexpressing goats possess increased numbers of secondary hair follicles (SHF). Our results indicate that Tß4-overexpression in cashmere goats could be a feasible strategy to increase cashmere yield.

Intraspecific predator interference promotes biodiversity in ecosystems.

Explaining biodiversity is a fundamental issue in ecology. A long-standing puzzle lies in the paradox of the plankton: many species of plankton feeding on a limited variety of resources coexist, apparently flouting the competitive exclusion principle (CEP), which holds that the number of predator (consumer) species cannot exceed that of the resources at a steady state. Here, we present a mechanistic model and demonstrate that intraspecific interference among the consumers enables a plethora of consumer species to coexist at constant population densities with only one or a handful of resource species. This facilitated biodiversity is resistant to stochasticity, either with the stochastic simulation algorithm or individual-based modeling. Our model naturally explains the classical experiments that invalidate the CEP, quantitatively illustrates the universal S-shaped pattern of the rank-abundance curves across a wide range of ecological communities, and can be broadly used to resolve the mystery of biodiversity in many natural ecosystems.

The surface waters of the ocean are teeming with microscopic creatures known as plankton, which get carried across vast distances by the currents. In a single ecosystem, thousands of plankton species may coexist, all competing for very few types of food sources. According to the principle of competitive exclusion, this should not be the case. Indeed, this theory states that the population levels of two species competing for the same resource cannot remain steady over time ­ or more generally, that the number of consumer species in an ecosystem cannot be higher than the number of resource types on which they rely. And yet, the Earth abounds with examples where a limited variety of resources supports a large number of competing yet coexisting consumer species. This is known as the paradox of the plankton. Many models have been proposed to explain how the limitations set by the competitive exclusion principle can be overcome, yet it is still unknown how to resolve the paradox of the plankton in a steady environment. In response, Kang et al. set out to test whether a phenomenon known as predator interference, which emerges when two or more individuals of the same species compete for the same resources, could help address the paradox of the plankton. To test this idea, Kang et al. developed a mathematical model of predator interference for multiple species of plankton feeding on a limited variety of food sources. The model put predators of the same species into encountering pairs to simulate predator interference. In this scenario, a wide range of predator species were able to live alongside each other with the numbers of each type of predator remaining steady over time. These results can be understood as follows: as a species becomes more successful at extracting resources from its environment, its population grows ­ and with it, the number of individuals engaged in intraspecific interference. Locked in interference, these species become less effective at getting food, creating a negative feedback loop that slows down the expansion of the species, allowing others to occupy the same niche. These findings may benefit ecologists and conservationists by offering insights into how to maintain biodiversity in ecosystems and protect endangered species. Further work is needed to test how well the model applies to other ecosystems.

Multiparametric and accurate functional analysis of genetic sequence variants using CRISPR-Select.

Determining the functional role of thousands of genetic sequence variants (mutations) associated with genetic diseases is a major challenge. Here we present clustered regularly interspaced short palindromic repeat (CRISPR)-SelectTIME, CRISPR-SelectSPACE and CRISPR-SelectSTATE, a set of flexible knock-in assays that introduce a genetic variant in a cell population and track its absolute frequencies relative to an internal, neutral control mutation as a function of time, space or a cell state measurable by flow cytometry. Phenotypically, CRISPR-Select can thereby determine, for example, pathogenicity, drug responsiveness/resistance or in vivo tumor promotion by a specific variant. Mechanistically, CRISPR-Select can dissect how the variant elicits the phenotype by causally linking the variant to motility/invasiveness or any cell state or biochemical process with a flow cytometry marker. The method is applicable to organoids, nontransformed or cancer cell lines. It is accurate, quantitative, fast and simple and works in single-well or 96-well higher throughput format. CRISPR-Select provides a versatile functional variant assay for research, diagnostics and drug development for genetic disorders.

CRISPR/Cas9-mediated VDR knockout plays an essential role in the growth of dermal papilla cells through enhanced relative genes.

BACKGROUND: Hair follicles in cashmere goats are divided into primary and secondary hair follicles (HFs). HF development, which determines the morphological structure, is regulated by a large number of vital genes; however, the key functional genes and their interaction networks are still unclear. Although the vitamin D receptor (VDR) is related to cashmere goat HF formation, its precise effects are largely unknown. In the present study, we verified the functions of key genes identified in previous studies using hair dermal papilla (DP) cells as an experimental model. Furthermore, we used CRISPR/Cas9 technology to modify the VDR in DP cells to dissect the molecular mechanism underlying HF formation in cashmere goats. RESULTS: The VDR expression levels in nine tissues of Shaanbei white cashmere goats differed significantly between embryonic day 60 (E60) and embryonic day 120 (E120). At E120, VDR expression was highest in the skin. At the newborn and E120 stages, the VDR protein was highly expressed in the root sheath and hair ball region of Shaanbei cashmere goats. We cloned the complete CDS of VDR in the Shaanbei white cashmere goat and constructed a VDR-deficient DP cell model by CRISPR/Cas9. Heterozygous and homozygous mutant DP cells were produced. The growth rate of mutant DP cells was significantly lower than that of wild-type DP cells (P < 0.05) and VDR mRNA levels in DP cells decreased significantly after VDR knockdown (P < 0.05). Further, the expression levels of VGF, Noggin, Lef1, and ß-catenin were significantly downregulated (P < 0.05). CONCLUSIONS: Our results indicated that VDR has a vital role in DP cells, and that its effects are mediated by Wnt and BMP4 signaling.

Multiplex Gene Editing via CRISPR/Cas9 System in Sheep.

Sheep is a major large animal model for studying development and disease in biomedical research. We utilized CRISPR/Cas9 system successfully to modify multiple genes in sheep. Here we provide a detailed protocol for one-cell-stage embryo manipulation by co-injecting Cas9 mRNA and RNA guides targeting three genes (MSTN, ASIP, and BCO2) to create genetic-modified sheep. Procedure described sgRNA design, construction of gRNA-Cas9 plasmid, efficient detection in fibroblast, embryos and sheep, and some manipulative technologies. Our findings suggested that the CRISPR/Cas9 method can be exploited as a powerful tool for livestock improvement by targeting multiple genes that are in charge of economically significant traits simultaneously.

Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats.

Cashmere goat skin contains two types of hair follicles (HF): primary hair follicles (PHF) and secondary hair follicles (SHF). Although multiple genetic determinants associated with HF formation have been identified, the molecules that determine the independent morphogenesis of HF in cashmere goats remain elusive. The growth and development of SHF directly influence the quantity and quality of cashmere production. Here, we report the transcriptome profiling analysis of nine skin samples from cashmere goats using 60- and 120-day-old embryos (E60 and E120, respectively), as well as newborns (NB), through RNA-sequencing (RNA-seq). HF morphological changes indicated that PHF were initiated at E60, with maturation from E120, while differentiation of SHF was identified at E120 until formation of cashmere occurred after birth (NB). The RNA-sequencing analysis generated over 20.6 million clean reads from each mRNA library. The number of differentially expressed genes (DEGs) in E60 vs. E120, E120 vs. NB, and E60 vs. NB were 1,024, 0 and 1,801, respectively, indicating that no significant differences were found at transcriptomic levels between E120 and NB. Key genes including B4GALT4, TNC, a-integrin, and FGFR1, were up-regulated and expressed in HF initiation from E60 to E120, while regulatory genes such as GPRC5D, PAD3, HOXC13, PRR9, VSIG8, LRRC15, LHX2, MSX-2, and FOXN1 were up-regulated and expressed in HF keratinisation and hair shaft differentiation from E120 and NB to E60. Several genes belonging to the KRT and KRTAP gene families were detected throughout the three HF developmental stages. The transcriptional trajectory analyses of all DEGs indicated that immune privilege, glycosaminoglycan biosynthesis, extracellular matrix receptor interaction, and growth factor receptors all played dominant roles in the epithelial-mesenchymal interface and HF formation. We found that the Wnt, transforming growth factor-beta/bone morphogenetic protein, and Notch family members played vital roles in HF differentiation and maturation. The DEGs we found could be attributed to the generation and development of HF, and thus will be critically important for improving the quantity and quality of fleece production in animals for fibres.

Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits.

The goat (Capra hircus) is one of the first farm animals that have undergone domestication and extensive natural and artificial selection by adapting to various environments, which in turn has resulted in its high level of phenotypic diversity. Here, we generated medium-coverage (9-13×) sequences from eight domesticated goat breeds, representing morphologically or geographically specific populations, to identify genomic regions representing selection signatures. We discovered ~10 million single nucleotide polymorphisms (SNPs) for each breed. By combining two approaches, ZHp and di values, we identified 22 genomic regions that may have contributed to the phenotypes in coat color patterns, body size, cashmere traits, as well as high altitude adaptation in goat populations. Candidate genes underlying strong selection signatures including coloration (ASIP, KITLG, HTT, GNA11, and OSTM1), body size (TBX15, DGCR8, CDC25A, and RDH16), cashmere traits (LHX2, FGF9, and WNT2), and hypoxia adaptation (CDK2, SOCS2, NOXA1, and ENPEP) were identified. We also identified candidate functional SNPs within selected genes that may be important for each trait. Our results demonstrated the potential of using sequence data in identifying genomic regions that are responsible for agriculturally significant phenotypes in goats, which in turn can be used in the selection of goat breeds for environmental adaptation and domestication.

Disruption of FGF5 in Cashmere Goats Using CRISPR/Cas9 Results in More Secondary Hair Follicles and Longer Fibers.

Precision genetic engineering accelerates the genetic improvement of livestock for agriculture and biomedicine. We have recently reported our success in producing gene-modified goats using the CRISPR/Cas9 system through microinjection of Cas9 mRNA and sgRNAs targeting the MSTN and FGF5 genes in goat embryos. By investigating the influence of gene modification on the phenotypes of Cas9-mediated goats, we herein demonstrate that the utility of this approach involving the disruption of FGF5 results in increased number of second hair follicles and enhanced fiber length in Cas9-mediated goats, suggesting more cashmere will be produced. The effects of genome modifications were characterized using H&E and immunohistochemistry staining, quantitative PCR, and western blotting techniques. These results indicated that the gene modifications induced by the disruption of FGF5 had occurred at the morphological and genetic levels. We further show that the knockout alleles were likely capable of germline transmission, which is essential for goat population expansion. These results provide sufficient evidences of the merit of using the CRISPR/Cas9 approach for the generation of gene-modified goats displaying the corresponding mutant phenotypes.

Correction: Disruption of FGF5 in Cashmere Goats Using CRISPR/Cas9 Results in More Secondary Hair Follicles and Longer Fibers.

[This corrects the article DOI: 10.1371/journal.pone.0164640.].

Multiplex gene editing via CRISPR/Cas9 exhibits desirable muscle hypertrophy without detectable off-target effects in sheep.

The CRISPR/Cas9 system provides a flexible approach for genome engineering of genetic loci. Here, we successfully achieved precise gene targeting in sheep by co-injecting one-cell-stage embryos with Cas9 mRNA and RNA guides targeting three genes (MSTN, ASIP, and BCO2). We carefully examined the sgRNAs:Cas9-mediated targeting effects in injected embryos, somatic tissues, as well as gonads via cloning and sequencing. The targeting efficiencies in these three genes were within the range of 27-33% in generated lambs, and that of simultaneously targeting the three genes was 5.6%, which demonstrated that micro-injection of zygotes is an efficient approach for generating gene-modified sheep. Interestingly, we observed that disruption of the MSTN gene resulted in the desired muscle hypertrophy that is characterized by enlarged myofibers, thereby providing the first detailed evidence supporting that gene modifications had occurred at both the genetic and morphological levels. In addition, prescreening for the off-target effect of sgRNAs was performed on fibroblasts before microinjection, to ensure that no detectable off-target mutations from founder animals existed. Our findings suggested that the CRISPR/Cas9 method can be exploited as a powerful tool for livestock improvement by simultaneously targeting multiple genes that are responsible for economically significant traits.

Generation of gene-modified goats targeting MSTN and FGF5 via zygote injection of CRISPR/Cas9 system.

Recent advances in the study of the CRISPR/Cas9 system have provided a precise and versatile approach for genome editing in various species. However, the applicability and efficiency of this method in large animal models, such as the goat, have not been extensively studied. Here, by co-injection of one-cell stage embryos with Cas9 mRNA and sgRNAs targeting two functional genes (MSTN and FGF5), we successfully produced gene-modified goats with either one or both genes disrupted. The targeting efficiency of MSTN and FGF5 in cultured primary fibroblasts was as high as 60%, while the efficiency of disrupting MSTN and FGF5 in 98 tested animals was 15% and 21% respectively, and 10% for double gene modifications. The on- and off-target mutations of the target genes in fibroblasts, as well as in somatic tissues and testis of founder and dead animals, were carefully analyzed. The results showed that simultaneous editing of several sites was achieved in large animals, demonstrating that the CRISPR/Cas9 system has the potential to become a robust and efficient gene engineering tool in farm animals, and therefore will be critically important and applicable for breeding.