Skip to main content
Download PDF
  • Methodology article
  • Open access
  • Published:

Engineered Cas12a-Plus nuclease enables gene editing with enhanced activity and specificity

BMC Biology volume 20, Article number: 91 (2022) Cite this article

Abstract

Background

The CRISPR-Cas12a (formerly Cpf1) system is a versatile gene-editing tool with properties distinct from the broadly used Cas9 system. Features such as recognition of T-rich protospacer-adjacent motif (PAM) and generation of sticky breaks, as well as amenability for multiplex editing in a single crRNA and lower off-target nuclease activity, broaden the targeting scope of available tools and enable more accurate genome editing. However, the widespread use of the nuclease for gene editing, especially in clinical applications, is hindered by insufficient activity and specificity despite previous efforts to improve the system. Currently reported Cas12a variants achieve high activity with a compromise of specificity. Here, we used structure-guided protein engineering to improve both editing efficiency and targeting accuracy of Acidaminococcus sp. Cas12a (AsCas12a) and Lachnospiraceae bacterium Cas12a (LbCas12a).

Results

We created new AsCas12a variant termed “AsCas12a-Plus” with increased activity (1.5~2.0-fold improvement) and specificity (reducing off-targets from 29 to 23 and specificity index increased from 92% to 94% with 33 sgRNAs), and this property was retained in multiplex editing and transcriptional activation. When used to disrupt the oncogenic BRAFV600E mutant, AsCas12a-Plus showed less off-target activity while maintaining comparable editing efficiency and BRAFV600E cancer cell killing. By introducing the corresponding substitutions into LbCas12a, we also generated LbCas12a-Plus (activity improved ~1.1-fold and off-targets decreased from 20 to 12 while specificity index increased from 78% to 89% with 15 sgRNAs), suggesting this strategy may be generally applicable across Cas12a orthologs. We compared Cas12a-Plus, other variants described in this study, and the reported enCas12a-HF, enCas12a, and Cas12a-ultra, and found that Cas12a-Plus outperformed other variants with a good balance for enhanced activity and improved specificity.

Conclusions

Our discoveries provide alternative AsCas12a and LbCas12a variants with high specificity and activity, which expand the gene-editing toolbox and can be more suitable for clinical applications.

Background

The clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR associated protein (Cas) system, an adaptive immunity system in bacteria and archaea, is a promising genome editing tool that has been widely used in a broad range of areas [1]. However, off-target cleavage of Cas-nucleases is routinely observed and remains an obstacle for clinical applications [2,3,4]. Therefore, the improvement of their targeting accuracy is essential for CRISPR-Cas tools in genome editing research, particularly in therapeutic applications [1]. So far, several strategies have been developed to improve this technique, and they can generally be divided into two categories: the sgRNA modification [5,6,7,8,9,10,11,12,13,14,15,16] and the Cas-nuclease protein engineering [17,18,19,20,21,22,23,24,25,26,27]. In particular, protein engineering is an efficient and widely used approach for the development of high-fidelity Cas-nuclease variants, which has been well-proved in SpCas9, such as the unbiased engineered variants HiFi-Cas9 [17], evoCas9 [18], xCas9 [19], Sniper-Cas9 [20], LZ3-Cas9 [21], etc., and the structure-guided engineered mutants SpCas9-HF [22], eSpCas9 [23], HeFSpCas9 [24], HypaCas9 [25], etc. These novel high-specific Cas-nuclease variants broaden the repertoire of CRISPR-Cas9 tools in gene, epigenome, and base editing applications [28].

As- (Acidaminococcus sp.) and Lb- (Lachnospiraceae bacterium) Cas12a are the two commonly used Cas-nuclease in the CRISPR-Cas12a system, which is also a promising genome editing tool in addition to the extensively investigated CRISPR-Cas9 [29, 30]. Several unique features make Cas12a distinguished from Cas9. First, Cas12a recognizes T-rich PAMs and generates sticky break ends [29], which makes it a complement to Cas9 in genome editing and broadens the genomic targeting scope. Second, Cas12a is a single crRNA-guided endonuclease and has the ribonuclease activity to process its pre-crRNA into mature crRNA [29, 31], which enables multiplex editing in a single crRNA transcript [32, 33]. Third, rather than using both RuvC and HNH domains in Cas9 [34], Cas12a cuts target DNA with a single RuvC domain. Fourth, Cas12a possesses the ability to trans-cleave single-stranded DNA (ssDNA) [35], making it a powerful platform for nucleic acid detection [35, 36]. Finally, Cas12a displays less off-target nuclease activity than does Cas9 [37, 38], enabling more precise genome editing for therapeutic applications [39].

Although Cas12a nuclease has shown powerful potentials in gene editing [40], insufficient efficiency and specificity remain a major obstacle for its broad application [37, 38]. Similar to SpCas9, approaches including Cas-protein engineering and sgRNA-modification had been employed in CRISPR-AsCas12a/LbCas12a system to improve their cutting efficiency or targeting accuracy [6,7,8,9,10,11,12,13,14,15,16, 26, 27, 41,42,43,44,45,46,47,48,49,50,51,52,53]. However, none of them achieved both improved specificity and enhanced activity (Additional file 1: Table S1) [6,7,8,9,10,11,12,13,14,15,16, 26, 27, 41,42,43,44,45,46,47,48,49,50,51,52,53], calling for the development of new variants to improve this system. Moreover, other Cas12a orthologs also have been explored and engineered for the same purpose (Additional file 1: Table S2) [42, 54,55,56,57,58,59,60,61,62,63], reflecting the urgent needs in this field. Here, we developed high-active and high-specific Cas12a nuclease variants by structure-guided protein engineering, which can expand the CRISPR-Cas toolbox and provide new genome editing tools for the applications in fundamental research and translational medicine.

Results

Structure-guided protein engineering for high-fidelity AsCas12a variants

To generate the high-fidelity AsCas12a mutant, we employed the energy hypothesis [22, 23], because numerous high specific Cas-nuclease variants had been successfully developed based on this hypothesis, including Cas9-HF [22] and eSpCas9 [23]. According to the crystal structure of the AsCas12a-crRNA-target DNA complex [64], we identified three positively charged amino acid residues (K780, R951, and R955) forming hydrogen bond contact with the backbone of the target DNA strand (Fig. 1a). R951 and R955 are located in the Bridge helix and K780 is in the WED domain [64], all of which are likely to be involved in stabilizing the interaction between AsCas12a and the targeted strand and thus leading to cleave mismatched off-target sites [22, 23, 64]. We then mutated these three amino acids and combined them for constructing five different AsCas12a variants bearing single, double, or triple substitutions (Fig. 1b), and tested whether these mutants possessed a relatively higher specificity. Western blotting showed that these mutants were expressed equivalently to the wild-type (WT) AsCas12a (Additional file 2: Figure S1a). Then, we performed the editing assay by using mismatched sgRNAs targeting DNMT1-site3, since this site had been well-studied in wild-type AsCas12a with different mismatched crRNAs [38]. Using deep-sequencing, T7 endonuclease I (T7E1), and polyacrylamide gel electrophoresis (PAGE)-based methods, we found that the AsCas12a-KK and AsCas12a-KA mutants retained comparable on-target activities but fewer cleavages with mismatched sgRNAs (Additional file 2: Figure S1b-d), indicating they may have a higher specificity. Based on this, we chose AsCas12a-KK and AsCas12a-KA variants for further study.

Fig. 1
figure 1

Generation of AsCas12a variants with increased specificity by weakening non-specific DNA contacts. a Schematic of wild-type AsCas12a interaction with the target DNA-sgRNA duplex. b The mutation sites of the AsCas12a variants. c Tag-seq-based comparative analyses of wild-type AsCas12a (WT), AsCas12a variant KA (KA), and AsCas12a variant KK (KK) with seventeen sgRNAs targeting nine genes (also see Additional file 2: Figure S1e). The sgRNA was shown on the top and the on-target and the off-target cleavages were displayed without or with mismatches to the sgRNA reference by color highlighting. Sequencing read counts were shown to the right of each site. d Total number of off-target sites detected with the seventeen sgRNAs. e Specificity Index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads within the seventeen sites). f Normalization of on-target activity of KK and KA to wild-type AsCas12a

Full size image

Genome-wide specificity of AsCas12a-KK and AsCas12a-KA

To globally evaluate the editing specificity of AsCas12a-KK and AsCas12a-KA, we performed Tag-seq experiments [65] to assess seventeen different sgRNAs targeting different sites in the endogenous human EMX1, DNMT1, RUNX1, PD1, CTLA4, CD47, SIRPa, CCR5, and CXCR4 genes (Fig. 1c and Additional file 2: Figure S1e), as these sites had been well-studied or were of clinical relevance. As a result, Tag-seq showed that the off-target cleavage was significantly decreased for AsCas12a-KK and AsCas12a-KA variants at most of the tested sites, with reducing the total off-target sites from 31 with AsCas12a-WT to 13 with AsCas12a-KK and to 8 with AsCas12a-KA (Fig. 1d). As expected, the specificity of the two mutants was increased, with the specificity index of 0.85, 0.82, and 0.71 for KA, KK, and WT, respectively (Fig. 1e). However, AsCas12a-KA was less active with only an average of 82% editing efficiency compared to the WT AsCas12a (Fig. 1f).

Improvement of the AsCas12a mutants for genome editing

Inspired by the enhanced AsCas12a (enAsCas12a), a highly active AsCas12a mutant previously reported [27], we next constructed AsCas12a-RKA mutant by introducing the E174R substitution (Fig. 2a), because this site was proximal to PAM DNA [64] and the charged arginine residue mutation could alter or form novel PAM proximal DNA contacts, which had been proved to increase the editing activity of the AsCas12a nuclease [27]. Western blotting showed that this substitution did not affect protein expression (Additional file 2: Figure S2a). Then, we performed specificity comparative analyses among AsCas12a-WT, AsCas12a-RKA, and AsCas12a-HF (enAsCas12a-HF, the high-fidelity version of enAsCas12a) by targeting RUNX1 and Site 6, two well-studied sites for specificity assessment of the CRISPR-Cas12a system [38]. Tag-seq indicated that the AsCas12a-RKA increased the editing efficiency with ~3-fold promotion (on-target reads was 2922 in RKA while WT is 1059) at site RUNX1 and at least equivalent editing level at Site 6 (2054 reads in RKA versus 1951 in WT) (Additional file 2: Figure S2b, c). More importantly, unlike the HF mutant, the activity improvement of the RKA mutant did not compromise but slightly improved the specificity (Additional file 2: Figure S2b-e).

Fig. 2
figure 2

Generation of AsCas12a variants with increased efficiency by introducing the high-activity substitution. a The mutation sites of the high active AsCas12a variants. b Tag-seq-based comparative analyses of wild-type AsCas12a (WT), AsCas12a variant RKA (RKA), AsCas12a variant RKK (RKK), and AsCas12a-HF (HF, the reported high-fidelity variant of enAsCas12a) with twenty-two sgRNAs targeting twelve genes (also see Additional file 2: Figure S3). c Normalization of on-target activity of RKA, RKK, and HF to wild-type AsCas12a. d Total number of off-target sites detected with the twenty-two sgRNAs. e Specificity index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads within the twenty-two sites). f Detection of the editing abilities for the non-canonical PAM with AsCas12a-RKA. Mean values are presented with SEM, n=4 independent experiments. Indel was revealed by Deep-seq

Full size image

Next, to determine whether this strategy could be applied to another variant, KK, we constructed AsCas12a-RKK as well and examined the protein expression level (Fig. 2a and Additional file 2: Figure S2a). Then, we utilized Tag-seq with twenty-two sgRNAs targeting twelve genes to comprehensively assess the editing abilities among AsCas12a-WT, AsCas12a-RKA, AsCas12a-RKK, and AsCas12a-HF (Fig. 2b and Additional file 2: Figure S3). Expectedly, with the introduction of the E174R, the average efficiency of RKA was increased about 1.64-fold compared to the WT AsCas12a (Fig. 2c). Moreover, among the four tested enzymes, RKA exhibited the highest specificity with the least off-target sites (Fig. 2d, e). Notably, the improvement of activity was also observed in mutant RKK (Fig. 2c); however, its specificity was slightly affected (Fig. 2d, e). Consistent with the previous report [27], enAsCas12a-HF displayed a robust efficiency which was ~2.23-fold improvement compared to WT AsCas12a (Fig. 2c); however, it induced 32 additional off-target cleavages in twenty-two tested sgRNAs and exhibited the lowest specificity (Fig. 2d, e). These data demonstrated that with the combination of E174R substitution, AsCas12a-RKA exhibited improved activity and slightly increased specificity, indicating it was a high-active and high-specific Cas12a nuclease. As the high-active site E174R displayed extended targeting range for non-canonical PAMs (such as ATTA, CTTA, GTTA, and TCTA) [27], we next tested whether AsCas12a-RKA possessed this ability. As shown in Fig. 2f, AsCas12a-RKA induced indels to an extent between enAsCas12a and WT AsCas12a, indicating a slightly expanded PAM recognition.

Apart from protein engineering, sgRNA modification is also an efficient way of enhancing the efficiency of the CRISPR-Cas12a system [66]. It had been reported that adding a “U4AU6” motif at the end of the crRNA [6] or using a pol-II-driven truncated pre-tRNA [7] to express the crRNA could improve the activity of the CRISPR-Cas12a system (Additional file 2: Figure S4a). Thus, we designed such sgRNAs to target exogenous EGFP gene and endogenous genes, FANCF and Site 6. FACS and Deep-seq results showed that the U4AU6-crRNA combined with some Cas12a variants tended to increase the editing efficiency at EGFP-g1 and site 6 loci, while the truncated pre-tRNA displayed comparable or less efficiency at all the four tested sites (Additional file 2: Figure S4b, c), suggesting that these two methods may work in a site-dependent manner, similar to a previous report of improving SpCas9 specificity by truncated-sgRNA [2, 5].

Multiplex editing of the AsCas12a variants using a single crRNA array

Next, we tested whether the new AsCas12a variants could improve Cas12a-based approaches. An advantage of the Cas12a enzyme over Cas9 is the multiplex editing, in which Cas12a processes individual crRNAs from a single crRNA array to simplify multiplex targeting in cells [32, 33, 47]. To assess this property in the engineered Cas12a variants, we cloned a poly-crRNA transcript including six crRNAs targeting DNMT1, EMX1, CTLA4, CCR5, SIPRa, and RUNX1 (Fig. 3a). Tag-seq experiments showed that all the mutants could mediate gene editing in these six sites with different levels (Fig. 3b), demonstrating the amino acid substitutions did not affect the crRNA self-processing activity. Among all the tested variants, RKA displayed both improved efficiency and specificity, while the enAsCas12a-HF showed the highest activity but with the worst specificity (Fig. 3b–e). These results demonstrated that AsCas12a-RKA was able to do multiplex editing with high activity and slightly improved specificity.

Fig. 3
figure 3

Multiplex editing with the AsCas12a variants. a Schematic of the multiplex editing. b Tag-seq-based comparative analyses of wild-type AsCas12a (WT), AsCas12a variants KA, KK, RKA, RKK, and HF with a single crRNA array targeting six sites. c Total number of off-target sites detected with the six sgRNAs within a single crRNA array. d Specificity index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads counts within these six sites). e Normalization of on-target activity of KA, KK, RKA, RKK, and HF to wild-type AsCas12a in multiplex editing

Full size image

Transcriptional activation with the AsCas12a variants

Cas12a has also been used for transcriptional activation of endogenous genes by fusing DNase-inactive Cas12a (dCas12a) to a gene activator [67, 68]. We then examined this application with the engineered AsCas12a variants. We found that using dAsCas12a fused to the synthetic VPR (VP64-p65-Rta) activation domain (dAsCas12a-VPR) (Fig. 4a), the dRKA-, dRKK-, and dHF-VPR systems can activate the transcriptional expression of IL1RN, MOYD, and HBG in human cells and Fgf21 in mouse cells with comparable level to dWT-VPR (Fig. 4b). However, dKK-VPR and dKA-VPR showed much lower capability to activate HBG and even failed to activate Fgf21 (Fig. 4b), which might reflect that the binding ability at these two sites was remarkably attenuated [64, 69]. Next, we performed RNA-seq for specificity comparison among dWT-, dRKA-, and dHF-VPR to activate IL1RN. As a result, the dRKA-VPR system displayed a slight improvement in activation of the endogenous gene IL1RN compared with the WT and HF (Fig. 4c), demonstrating the ability of the RKA in transcriptional activation.

Fig. 4
figure 4

Transcriptional activation with the AsCas12a variants. a Schematic of the gene activation system based on catalytic-dead AsCas12a variants. VPR, synthetic VP64-p65-Rta activation domain. b qPCR analyses of the transcriptional activation levels with the AsCas12a variants guided by a single sgRNA targeting each promoter region of IL1RN, MYOD, and HBG in human HEK293T cells and Fgf21 in mouse B16 cells, respectively. Mean values are presented with SEM, n=3 independent experiments. **p<0.01. ns, no significance. Student’s t-test, dWT-VPR sample versus corresponding dAsCas12a-VPR variant samples. c RNA-seq assessing the specificity of the dWT-, dRKA-, and dHF-VPR systems to activate IL1RN gene. d Schematic of the orthogonal gene editing and activation system based on catalytic active AsCas12a variants. e Relative mRNA expression (Activation) and genome editing efficiency (Editing) of the AsCas12a-VPR systems at MYOD and IL1RN sites with indicated lengthed sgRNAs. 15-bp sgRNAs for gene activation and 23-bp sgRNAs for genome editing. Mean values are presented with SEM, n=3 independent experiments. Indel was revealed by Deep-seq

Full size image

Simultaneous orthogonal gene activation and genome editing for multiplex genes with catalytical active SpCas9 or AsCas12a fused to a gene activator has been reported [32, 70, 71]. Therefore, we compared orthogonally activation and editing ability among catalytically active AsCas12a-WT, AsCas12a-RKA, and AsCas12a-HF fused to VPR (termed WT-VPR, RKA-VPR, and HF-VPR system) by using short 15-bp sgRNA for gene activation and long 23-bp sgRNA for gene editing (Fig. 4d). To this end, we designed two experiments, one for MYOD activation and IL1RN editing, the other for IL1RN activation and MYOD editing. Deep-seq assays showed that all the three AsCas12a nucleases cleaved genomic DNA at the MYOD promoter region with comparable levels when using a 23-bp sgRNA, and failed to induce indels when using a 15-bp sgRNA. For activation, we observed that MYOD expression can be activated to a similar extent with the WT-, RKA-, and HF-VPR systems when transfected with 15-bp sgRNAs (Fig. 4e). And similar results were observed at the IL1RN site (Fig. 4e). However, we also noticed that 23-bp guides could activate transcription, although to a less extent than 15-bp guides (Fig. 4e), which was consistent with a previous report [71]. Cas12a usually cleaves DNA at around 18–23 bps distant from its PAM [29, 72], and DNA repair after cleavage could form a new imperfect matched sgRNA targeting site with mismatches at the sgRNA distal end. Since 15-bp sgRNA is able to activate gene expression, distal-mismatched 23-bp sgRNA might guide Cas12a to the repaired site using the proximal matched sequence and activate gene expression. Together, these data indicated that the RKA-VPR system could mediate gene activation and editing with different lengthed sgRNAs.

AsCas12a-RKA holds editing safety in disruption oncogenic BRAFV600E

BRAF-V600E (1799T>A) is one of the most frequently reported driver mutations in multiple types of cancers, and patients with such mutations could benefit from disrupting this mutant allele [39, 73]. However, a major concern for implementing CRISPR/Cas9 for gene therapy is the relatively high frequency of off-target effects. Therefore, we sought to examine the therapeutic potential of the high-fidelity Cas12a-RKA for editing this mutation. Melanoma cell line A375 is a homozygous genotype with BRAF-V600E [74] (Additional file 2: Figure S5a). By using a mut-sgRNA, Tag-PCR assay [65] roughly displayed that the AsCas12a nucleases (WT, RKA, and HF) retained high editing selectivity at this site, as they did not cut wild-type but mutated sequence, whereas Cas9 recognized and cut both wild-type and mutated alleles (Additional file 2: Figure S5b, c). Next, to more accurately assess the editing selectivity of the Cas-nucleases (SpCas9-WT, AsCas12a-WT, and all the AsCas12a mutants in this study), Tag-seq experiments were performed by using the Cas9- and Cas12a-sgRNA (both contained WT- and mut-sgRNAs) in both BRAF+/+ HEK293T and BRAFV600E/V600E A375 cells (Fig. 5a). Consistent with the Tag-PCR results, SpCas9 edited both wild-type and mutant BRAF with the mut-sgRNA and induced abundant off-targets editing (Fig. 5a), indicating a low selectivity for this mutation editing. In contrast, wild-type AsCas12a and the engineered variants displayed higher specificity with only few cleavages in wild-type BRAF and no off-targets detection in mutant BRAF when applied with the mut-sgRNA (Fig. 5a). However, among all the tested AsCas12a nucleases, variant RKA showed high-specificity and high-activity at this site (Fig. 5a). Further, by disruption BRAFV600E, AsCas12a-RKA induced A375 cell apoptosis with comparable level to AsCas12a-WT and AsCas12a-HF (Fig. 5b), demonstrating the therapeutic potential of AsCas12a-RKA to treat BRAFV600E tumors.

Fig. 5
figure 5

Disruption of the oncogenic BRAFV600E mutation with AsCas12a variants in human A375 melanoma cells. a Tag-seq assessing the editing specificity of AsCas12a-WT and AsCas12a variants for BRAF with WT- and mut-sgRNAs in HEK-293T Cell (BRAF +/+) and A375 Cell (BRAF V600E/V600E). b Flow cytometry detecting the cells apoptosis induced by disrupting the BRAFV600E mutation using AsCas12a-WT, AsCas12a-RKA, and AsCas12a-HF. Blank, A375 cells without treatment. NC, AsCas12a-WT, AsCas12a-RKA, and AsCas12a-HF mix + sgRNA-EGFP. WT/RKA/HF, AsCas12a-WT/AsCas12a-RKA/AsCas12a-HF + mut-BRAF-sgRNA. Mean values are presented with SEM, n=3 independent experiments

Full size image

Engineering of high-performance LbCas12a variants via analogous substitutions to AsCas12a-RKA

Encouraged by AsCas12a-RKA, we next examined whether these analogous positions in AsCas12a-RKA could be deployed in LbCas12a, another commonly used Cas12a nuclease, to generate LbCas12a mutants with high performance as well. Via amino acid sequence alignment between AsCas12a and LbCas12a, we identified the conserved amino acid residues and constructed four LbCas12a variants, KK, KA, RKA, and RKK (Fig. 6a, b). Western blotting showed similar protein expression levels of these variants (Additional file 2: Figure S6). Then, we compared the activity and specificity of the four LbCas12a variants to WT LbCas12a by Tag-seq with fifteen sgRNAs targeting nine human endogenous genes (Fig. 6c). Consistently, the total off-targets of LbCas12a-KK, LbCas12a-KA, LbCas12a-RKA, and LbCas12a-RKK reduced from 27 to 2, 2, 4, and 3, respectively (Fig. 6d). The specificity of these four mutants was increased as well (Fig. 6e). Notably, similar to AsCas12a-RKA, LbCas12a-RKA showed high efficiency and specificity (Fig. 6c–f).

Fig. 6
figure 6

Generation of high-performance LbCas12a variants. a The conserved sites in AsCas12a and LbCas12a, underlined with red. b The mutation sites of LbCas12a variants. c Tag-seq-based comparative analyses of wild-type LbCas12a (WT), LbCas12a variant KA (KA), LbCas12a variant KK (KK), LbCas12a variant RKA (RKA), and LbCas12a variant RKK (RKK) with fifteen sgRNAs targeting nine genes. d Total number of off-target sites detected with the fifteen sgRNAs. e Specificity index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads within the fifteen sites). f Normalization of on-target activity of LbCas12a-KA, LbCas12a-KK, LbCas12a-RKA, and LbCas12a-RKK to wild-type LbCas12a

Full size image

Collectively, all the above results demonstrated that the engineered Cas12a-RKA variant behaved as a high-active and high-specific nuclease; we hence termed it as “Cas12a-Plus”.

Systematical comparison of available high-active AsCas12a variants

When our manuscript was underwriting, a new AsCas12a variant, AsCas12a-ultra, was reported with significantly enhanced activity [52]. To systematically compare the performance among these high-active Cas12a variants, we constructed the AsCas12a-ultra, and the LbCas12a-ultra that was created by amino acid sequence conservation (Fig. 7a). Correspondingly, more variants were generated by combining with the high-fidelity mutant sites KK and KA, or by introducing the high-active substitutions, or by combining the RKA and RKK mutations (Fig. 7b). After confirming the comparable protein expression level by Western blotting (Additional file 2: Figures S7a and S8a), we comprehensively analyzed their performance using Tag-seq with twenty-eight sgRNAs targeting nineteen genes among AsCas12a variants (Additional file 2: Figure S7b), and with fifteen sgRNAs targeting nine genes among LbCas12a variants (Additional file 2: Figure S8b), respectively. As a result, although the mutant enAsCas12a and its high-fidelity version enAsCas12-HF exhibited the highest editing abilities, they induced numbers of extra off-targets, particularly the enAsCas12a (Fig. 7c–e). Surprisingly, the AsCas12a-ultra showed a slightly improved activity and an obvious decreased specificity (Fig. 7c–e), and disruption of mNeonGreen expression in HEK293T-KI reporter cell line leads to similar results (Additional file 2: Figure S9), which was different from the previous report [52]. We speculated that the possible reason was the difference in delivery approach. RNP delivery and plasmid transfection were employed by the previous report and the current study, respectively. And these two methods had been demonstrated to result in different efficiency of gene editing [75, 76]. When introduced the high-active mutation E174R into the AsCas12a-ultra (termed RU), efficiency was increased while specificity was significantly compromised (Fig. 7c–e). All the variants containing the KA and KK mutations exhibited improved specificity. Notably, among the ten tested variants of AsCas12a, AsCas12a-Plus displayed the best balance in editing performance with moderately enhanced activity and specificity (Fig. 7c–e). And similar results were observed in the LbCas12a versions (Fig. 7f–h, Additional file 2: Figures S8 and S10).

Fig. 7
figure 7

Performance comparison of Cas12a variants. a The conserved sites of ultra variant in AsCas12a and LbCas12a, underlined with red. b The mutation sites of AsCas12a and LbCas12a variants. RKA=Plus, UKA=ultra+KA, UKK=ultra+KK, RU=high-active substitution E174R/D156R+ultra, RKAU=RKA+ultra, RKKU=RKK+ultra. c–e Tag-seq-based comparative analyses of wild-type AsCas12a (WT), and AsCas12a variants with twenty-eight sgRNAs targeting nineteen genes (see Additional file 2: Figure S7b). c Total number of off-target sites detected with the twenty-eight sgRNAs. d Specificity index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads within the twenty-eight sites). e Normalization of on-target activity of AsCas12a variants to wild-type AsCas12a. f–h Tag-seq-based comparative analyses of wild-type LbCas12a (WT), and LbCas12a variants with fifteen sgRNAs targeting nine genes (see Additional file 2: Figure S8b). f Total number of off-target sites detected with the fifteen sgRNAs. g Specificity Index (value was calculated by the ratio of total on-target reads to the on-target reads plus the off-target reads within the fifteen sites). h Normalization of on-target activity of LbCas12a variants to wild-type LbCas12a

Full size image

Discussion

The off-target effect of the CRISPR-Cas genome editing tools is a major concern for therapeutic applications. It has been reported that Cas12a exhibits a higher specificity over the widely used SpCas9; however, the relatively low activity restricts their broad use [26, 27, 29, 30, 44]. Given the advantageous properties of the Cas12a nuclease, such as the higher specificity and distinct PAM preference [29, 31], it represents a powerful alternative for gene editing. Here, we constructed novel Cas12a variants termed “Cas12a-Plus” by rational structure-guided engineering to enable more potent and more specific gene editing.

Since Cas12a and Cas9 are complementary to each other as genome editing tools, like Cas9, Cas12a has also attracted plenty of effort for protein engineering to expand editing range, enhance activity, and improve specificity (for detailed information, please refer to Additional file 1: Table S1 and Table S2). For instance, AsCas12a-RVR and AsCas12a-RR have been engineered to expand the PAM recognition range and their high-fidelity derivative mutants AsCas12a-RVRA and AsCas12a-RRA can improve editing specificity but with compromised activity [26, 27, 41]. More recently, enhanced AsCas12a (enAsCas12a/enAsCas12a-HF) [27], AsCas12a-ultra [52], and the impLbCas12a [45] have been reported to significantly enhance editing efficiency and broaden editing ranges. However, our data (Figs. 2, 3, and 7) and the previous reports [27, 45, 52] demonstrate that these variants induce extra off-targets cleavages, indicating a compromise of specificity. In this study, we created alternative high-fidelity Cas12a-KA and KK, and high-active and high-specific Cas12a-Plus (Figs. 1, 2, 3, and 7), which expanded the Cas12a toolbox. Therefore, we recommend to use Cas12a variants with different properties according to the intended applications. Our study provides Cas12a-Plus as the first variant with enhanced activity and increased specificity, which holds great potential for broad applications, especially for clinical disease therapy.

With unique features, Cas12a outperforms Cas9 in some applications. For example, Cas12a has been reported to be better than Cas9 for one-step generation of modular CAR-T cells [77]. With the crRNA self-processing activity, Cas12a could be used for combinatorial genetic screening [50, 51]. In this study, we also found that, with higher specificity, Cas12a outperformed Cas9 to disrupt BRAFV600E mutated allele and thus triggered cancer cell death with much less genome editing in normal cells (Fig. 5). With enhanced specificity and activity, the AsCas12a-Plus variant further improved the performance (Fig. 5). Since it maintained the enhanced-activity and improved-specificity in multiplex editing and transcriptional activation (Figs. 3 and 4), we believe that AsCas12a-Plus could outperform wild-type AsCas12a in combinatorial genetic screening as well as other Cas12a-based applications.

According to the energy hypothesis [22, 23], the Watson-Crick base pairing between gRNA and the target DNA strand as well as the binding between Cas protein and the PAM bases provided specific energy (A, T, C, and G base-dependent), while the binding between Cas protein and the backbone of target DNA strand and PAM DNA as well as the binding between Cas protein and non-target DNA strand provided non-specific energy (base-independent). Both specific and non-specific energy contributed to recognition and cleavage. Besides, it has been reported that Bridge helix arginines (Rs) play a critical role in sensitivity to mismatched sequences [69]. Based on these, we finally focused on R951 and R955 by analyzing the crystal structure of the AsCas12a-crRNA-target DNA complex [64]. Because they were located at Bridge helix and seemed to provide non-specific contacts [64]. Since R951A mutation reduces Cas12a activity [64], we thus used R951K as an alternative, since lysine (K) is highly similar to arginine (R) in structure and is of less possibility to form hydrogen bonds between Cas12a and the target DNA strand as predicted by the complex structure. We did demonstrate that R951K/R955A (KA) and R951K/R955K (KK) mutation could improve editing specificity (Fig. 1). By combining with a known activity-enhancing E174R mutation [27], we created AsCas12a-Plus with high-activity and high-fidelity, which was also applicable with LbCas12a (Figs. 2 and 6). Other variants combined with the RKA mutation, such as the AsCas12a-ultra, could improve activity and specificity (the RU and RKAU variants, Fig. 7). In theory, this strategy could be introduced into other Cas12a variants, such as the PAM-less-restricted RVR and RR mutations [26, 27, 41], or combined with the chemical modification [10, 15], to improve performance. Amino acid residues other than the tested E174/K780/R951/R955 could contribute to the activity and specificity of AsCas12a, such as the residues contacting the non-target DNA strand or the crRNA, and thus mutating these residues might improve AsCas12a performance. All of these hypotheses need to be tested in the future.

Conclusions

In summary, we created novel AsCas12a and LbCas12a variants with both high-activity and high-fidelity, expanding the Cas12a toolbox, and thus, these variants could enhance the performance of Cas12a in a wide spectrum of applications.

Methods

Plasmid construction

AsCas12a and LbCas12a variants expression plasmids bearing amino acid substitutions were generated by standard PCR and molecular cloning into a plasmid contained a CAG promoter, HA, P2A-mcherry cassette via Gibson Assembly. sgRNA expression plasmids were constructed by ligating oligonucleotide duplexes into EcoR V and Hind III cut pBlueScript backbone with a human U6 promoter and an As- or Lb-crRNA sequence. All the plasmids were confirmed by Sanger sequencing, and all the sgRNAs used in this study are shown in Additional file 1: Table S3.

Cell culture and cell transfection

HEK293T, B16, and A375 cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM, Life Technologies) at 37°C in a 5% CO2 humidified incubator. All growth media were supplemented with 2 mM L-glutamine (Life Technologies), 100 U/mL penicillin, 100 μg/mL streptomycin (Life Technologies), and 10% fetal bovine serum. All the cell lines in this study were cultured no more than 10 passages.

Cells were transfected with PEI reagent (Polysciences, Inc., PA, USA) according to the manufacturer’s instructions. Briefly, 250 ng of pCAG-Cas12a-mcherry and 250 ng of sgRNA-encoding plasmids were transfected per well in a 24-well plate. Cells were harvested 2–3 days after transfection, then the genomic DNA or the total RNA were extracted for the following assays.

Tag-seq method

Tag-seq experiments were performed and analyzed as previously described [65]. Briefly, HEK293T cells were transfected by PEI with 20 nM Tag, 1000 ng of Cas nuclease, and 1000 ng single sgRNA or a pool sgRNAs (30–50 ng/sgRNA) per well in a six-well plate. A375 cells were transfected by Amaxa Cell Line Nucleofector Kit V (VCA-1003, LONZA, Switzerland) following the manufacturer’s instructions (2D) with 20 nM Tag, 1200 ng of Cas nuclease, and 800 ng WT/Mut-BRAF-sgRNA. All cells were harvested 3 days after transfection and genomic DNA was extracted for one-step libraries preparation by the Fragmentation, End Preparation, and dA-Tailing Module and Adapter Ligation Module kit (Vazyme Biotech Co., Ltd., Nanjing, China). The R and L libraries were constructed by PCR with library preparation primers, which were followed by sequencing (Hiseq/NovaSeq platform, Novogene, Beijing, China) and analysis with a Tag-seq bioinformatics pipeline. Tag-seq experiments were performed with the same input gDNA and an equal sequencing depth. The analysis pipeline is available at https://github.com/zhoujj2013/Tag-seq and https://doi.org/10.5281/zenodo.4679460.

Activity and specificity scoring

For the comparisons of performance among Cas12a variants, Tag-seq reads were used for calculating the editing activity and targeting specificity. Activity scores were calculated as the mean ratio of the on target reads across all the tested sites, normalized to the WT Cas12a nuclease. Specificity scores were calculated as the ratio of the on target reads to the on-target reads plus the off-target reads across all the tested sites.

Tag-PCR assay

Tag-PCR was used to roughly determine the editing events of the CRISPR-Cas systems, which reflected the insertion efficiency of the Tag at the editing sites [65]. Briefly, cell transfection was the same as the Tag-seq method. After extraction of genomic DNA with integrated Tag sequence at break sites, PCR was performed by using the Tag-specific primer and a locus-specific R primer, then PCR products were assessed by running on an agarose gel.

Deep-seq analysis

Deep-seq was used to determine the indel frequency. Briefly, the primers with forward and reverse indexes were used to amplify the genomic regions in the first-round PCR. Then, equal amounts of the first PCR products were mixed and subjected to a second round of PCR with the P5- and P7-containing primers to generate the sequencing libraries. Paired-end sequencing was performed using the Hiseq/NovaSeq system (Novogene, Beijing, China). Indel frequency was calculated as the ratio of (read counts with indel sequence)/(total sequencing read counts). The deep-seq primers and the samples’ index information were listed in Additional file 1: Tables S4 and S5.

T7EI and polyacrylamide gel electrophoresis (PAGE) assay

For T7EI analysis, the amplicons were purified, denatured at 95°C for 5 min and annealed in NEB Buffer 2 with a slow ramp down (approximately −2°C/min) to 4°C, then subjected to T7 endonuclease I (NEB, UK) digestion for 3 h at 37°C before loading on a 2% agarose gel. For PAGE assay, genomic DNA was isolated using sarkosyl lysis buffer (10 mM Tris pH7.6, 0.5% Sarkosyl, 10 mM NaCl, 10 mM EDTA, 0.1 mg/ml proteinase K) and the target sites were amplified by PCR. The purified amplicons were reannealed to form heteroduplexes and then subjected to 5% polyacrylamide gel electrophoresis. All T7EI and the PAGE primers in this study are listed in Additional file 1: Table S4.

Quantitative real-time PCR

Total RNA from the transfected cells was isolated using Trizol Reagent (Thermo Fisher, USA) following the manufacturer’s instructions. Total RNA (1 μg) was reverse transcribed into cDNA and then quantitative real-time PCR (SYBR Premix Ex Taq II, TAKARA, China) was performed using a LightCycler 96 System (Roche, Switzerland). Relative gene expression was calculated using the 2−ΔΔCt method after normalizing to GAPDH expression. All the qPCR primers are listed in Additional file 1: Table S4.

Western blotting

To detect the expression of the AsCas12a and LbCas12a variants, the transfected cells were lysed in a 2×SDS loading buffer and boiled for 10 min. Lysates were resolved through SDS/PAGE and transferred onto a nitrocellulose membrane which was blocked using 5% non-fat milk and sequentially incubated with primary antibodies (anti-HA, sigma, USA, anti-GADPH, Proteintech, China) and an HRP-conjugated horse anti-mouse IgG secondary antibody (CST, USA, CAT# 7076S). All the probed proteins were finally detected through chemiluminescence following the manufacturer’s instructions (Pierce, USA).

RNA-seq

RNA-seq experiments were performed and analyzed as previously described [78]. Briefly, total RNA was extracted by Trizol reagent (Invitrogen, Carlsbad, CA, USA), then mRNA was used for the standard RNA libraries’ preparation, and libraries were sequenced by 150 bp paired-end Novaseq device. For data analysis, Hisat2 v2.0.52 was used to build the index of the reference genome and align the paired-end clean reads with the reference genome. Then, StringTie v2.23 was used to count the read numbers mapped to each gene. Fragments per kilobase per million (FPKM) of each gene was calculated based on the length of the gene and the reads count mapped to this gene. Differential expression was defined by a Benjamini-Hochberg adjusted p-value (q value | FDR) of <0.05 and a fold change of >2 or <0.5. All figures were plotted using R package ggplot2.

FACS analysis

All flow cytometry analyses were performed using FlowJo software (TreeStar, USA). Cells were harvested 48 h post-transfection, and the cleavage efficiency of AsCas12a variants was determined as the proportion of GFP negative cells within the AsCas12a-transfected cells (mCherry-positive). To detect the apoptosis of A375, cells were first transduced with the lentivirus encoding the AsCas12a-WT, AsCas12a-RKA, and AsCas12a-HF protein, then cells were co-transfected with the mut-BRAF or EGFP (as a negative control) sgRNA and a pCMV-mCherry reporter plasmid. After 7 days, cells were gated out using mCherry, followed by the standard procedures of the Annexin V-FITC Apoptosis Detection Kit (BestBio, China) according to the manufacturer’s instructions.

Statistics analysis and reproducibility

Student’s t-test and one-way ANOVA were used in this study for the statistical analysis. The reproducibility was shown by performing two-four independent biological replicate experiments.

Availability of data and materials

All data generated or analyzed during this study are included in this published article, the sequencing data are deposited at NCBI SRA (Bioproject PRJNA755186), and the supporting data values for the figures are provided at Additional file 3.

Abbreviations

CRISPR:

Clustered regularly interspaced short palindromic repeats

Cas:

CRISPR associated protein

HF:

High-fidelity

ssDNA:

Single-stranded DNA

T7E1:

T7 endonuclease I

PAGE:

Polyacrylamide gel electrophoresis

WT:

Wild-type

enAsCas12a:

Enhanced AsCas12a

dCas12a:

DNase-inactive Cas12a

PAM:

Protospacer-adjacent motif

References

  1. Pickar-Oliver A, Gersbach CA. The next generation of CRISPR-Cas technologies and applications. Nat Rev Mol Cell Biol. 2019;20(8):490–507. https://doi.org/10.1038/s41580-019-0131-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Tsai SQ, Zheng Z, Nguyen NT, Liebers M, Topkar VV, Thapar V, et al. GUIDE-seq enables genome-wide profiling of off-target cleavage by CRISPR-Cas nucleases. Nat Biotechnol. 2015;33(2):187–97. https://doi.org/10.1038/nbt.3117.

    Article  CAS  PubMed  Google Scholar 

  3. Wienert B, Wyman SK, Richardson CD, Yeh CD, Akcakaya P, Porritt MJ, et al. Unbiased detection of CRISPR off-targets in vivo using DISCOVER-Seq. Science. 2019;364(6437):286–9. https://doi.org/10.1126/science.aav9023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Lazzarotto CR, Malinin NL, Li Y, Zhang R, Yang Y, Lee G, et al. CHANGE-seq reveals genetic and epigenetic effects on CRISPR-Cas9 genome-wide activity. Nat Biotechnol. 2020;38(11):1317–27. https://doi.org/10.1038/s41587-020-0555-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Fu Y, Sander JD, Reyon D, Cascio VM, Joung JK. Improving CRISPR-Cas nuclease specificity using truncated guide RNAs. Nat Biotechnol. 2014;32(3):279–84. https://doi.org/10.1038/nbt.2808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Bin Moon S, Lee JM, Kang JG, Lee NE, Ha DI, Kim DY, et al. Highly efficient genome editing by CRISPR-Cpf1 using CRISPR RNA with a uridinylate-rich 3'-overhang. Nat Commun. 2018;9(1):3651. https://doi.org/10.1038/s41467-018-06129-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Zhang X, Xu L, Fan R, Gao Q, Song Y, Lyu X, et al. Genetic editing and interrogation with Cpf1 and caged truncated pre-tRNA-like crRNA in mammalian cells. Cell Discov. 2018;4:36. https://doi.org/10.1038/s41421-018-0035-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Park HM, Liu H, Wu J, Chong A, Mackley V, Fellmann C, et al. Extension of the crRNA enhances Cpf1 gene editing in vitro and in vivo. Nat Commun. 2018;9(1):3313. https://doi.org/10.1038/s41467-018-05641-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Gao Z, Herrera-Carrillo E, Berkhout B. Improvement of the CRISPR-Cpf1 system with ribozyme-processed crRNA. RNA Biol. 2018;15(12):1458–67. https://doi.org/10.1080/15476286.2018.1551703.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Li B, Zhao W, Luo X, Zhang X, Li C, Zeng C, et al. Engineering CRISPR-Cpf1 crRNAs and mRNAs to maximize genome editing efficiency. Nat. Biomed Eng. 2017;1(5). https://doi.org/10.1038/s41551-017-0066.

  11. Kweon J, Jang AH, Kim DE, Yang JW, Yoon M, Rim Shin H, et al. Fusion guide RNAs for orthogonal gene manipulation with Cas9 and Cpf1. Nat Commun. 2017;8(1):1723. https://doi.org/10.1038/s41467-017-01650-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ding D, Chen K, Chen Y, Li H, Xie K. Engineering introns to express RNA guides for Cas9- and Cpf1-mediated multiplex genome editing. Mol Plant. 2018;11(4):542–52. https://doi.org/10.1016/j.molp.2018.02.005.

    Article  CAS  PubMed  Google Scholar 

  13. Liu P, Luk K, Shin M, Idrizi F, Kwok S, Roscoe B, et al. Enhanced Cas12a editing in mammalian cells and zebrafish. Nucleic Acids Res. 2019;47(8):4169–80. https://doi.org/10.1093/nar/gkz184.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Kim H, Lee WJ, Oh Y, Kang SH, Hur JK, Lee H, et al. Enhancement of target specificity of CRISPR-Cas12a by using a chimeric DNA-RNA guide. Nucleic Acids Res. 2020;48(15):8601–16. https://doi.org/10.1093/nar/gkaa605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Ling X, Chang L, Chen H, Gao X, Yin J, Zuo Y, et al. Improving the efficiency of CRISPR-Cas12a-based genome editing with site-specific covalent Cas12a-crRNA conjugates. Mol Cell. 2021;81(22):4747–56 e7. https://doi.org/10.1016/j.molcel.2021.09.021.

    Article  CAS  PubMed  Google Scholar 

  16. Li X, Wang Y, Liu Y, Yang B, Wang X, Wei J, et al. Base editing with a Cpf1-cytidine deaminase fusion. Nat Biotechnol. 2018;36(4):324–7. https://doi.org/10.1038/nbt.4102.

    Article  CAS  PubMed  Google Scholar 

  17. Vakulskas CA, Dever DP, Rettig GR, Turk R, Jacobi AM, Collingwood MA, et al. A high-fidelity Cas9 mutant delivered as a ribonucleoprotein complex enables efficient gene editing in human hematopoietic stem and progenitor cells. Nat Med. 2018;24(8):1216–24. https://doi.org/10.1038/s41591-018-0137-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Casini A, Olivieri M, Petris G, Montagna C, Reginato G, Maule G, et al. A highly specific SpCas9 variant is identified by in vivo screening in yeast. Nat Biotechnol. 2018;36(3):265–71. https://doi.org/10.1038/nbt.4066.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hu JH, Miller SM, Geurts MH, Tang W, Chen L, Sun N, et al. Evolved Cas9 variants with broad PAM compatibility and high DNA specificity. Nature. 2018;556(7699):57–63. https://doi.org/10.1038/nature26155.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lee JK, Jeong E, Lee J, Jung M, Shin E, Kim YH, et al. Directed evolution of CRISPR-Cas9 to increase its specificity. Nat Commun. 2018;9(1):3048. https://doi.org/10.1038/s41467-018-05477-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Schmid-Burgk JL, Gao L, Li D, Gardner Z, Strecker J, Lash B, et al. Highly parallel profiling of Cas9 variant specificity. Mol Cell. 2020;78(4):794–800 e8. https://doi.org/10.1016/j.molcel.2020.02.023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Kleinstiver BP, Pattanayak V, Prew MS, Tsai SQ, Nguyen NT, Zheng Z, et al. High-fidelity CRISPR-Cas9 nucleases with no detectable genome-wide off-target effects. Nature. 2016;529(7587):490–5. https://doi.org/10.1038/nature16526.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Slaymaker IM, Gao L, Zetsche B, Scott DA, Yan WX, Zhang F. Rationally engineered Cas9 nucleases with improved specificity. Science. 2016;351(6268):84–8. https://doi.org/10.1126/science.aad5227.

    Article  CAS  PubMed  Google Scholar 

  24. Kulcsar PI, Talas A, Huszar K, Ligeti Z, Toth E, Weinhardt N, et al. Crossing enhanced and high fidelity SpCas9 nucleases to optimize specificity and cleavage. Genome Biol. 2017;18(1):190. https://doi.org/10.1186/s13059-017-1318-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Chen JS, Dagdas YS, Kleinstiver BP, Welch MM, Sousa AA, Harrington LB, et al. Enhanced proofreading governs CRISPR-Cas9 targeting accuracy. Nature. 2017;550(7676):407–10. https://doi.org/10.1038/nature24268.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Gao L, Cox DBT, Yan WX, Manteiga JC, Schneider MW, Yamano T, et al. Engineered Cpf1 variants with altered PAM specificities. Nat Biotechnol. 2017;35(8):789–92. https://doi.org/10.1038/nbt.3900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Kleinstiver BP, Sousa AA, Walton RT, Tak YE, Hsu JY, Clement K, et al. Engineered CRISPR-Cas12a variants with increased activities and improved targeting ranges for gene, epigenetic and base editing. Nat Biotechnol. 2019;37(3):276–82. https://doi.org/10.1038/s41587-018-0011-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Anzalone AV, Koblan LW, Liu DR. Genome editing with CRISPR-Cas nucleases, base editors, transposases and prime editors. Nat Biotechnol. 2020;38(7):824–44. https://doi.org/10.1038/s41587-020-0561-9.

    Article  CAS  PubMed  Google Scholar 

  29. Zetsche B, Gootenberg JS, Abudayyeh OO, Slaymaker IM, Makarova KS, Essletzbichler P, et al. Cpf1 is a single RNA-guided endonuclease of a class 2 CRISPR-Cas system. Cell. 2015;163(3):759–71. https://doi.org/10.1016/j.cell.2015.09.038.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zetsche B, Abudayyeh OO, Gootenberg JS, Scott DA, Zhang F. A survey of genome editing activity for 16 Cas12a orthologs. Keio J Med. 2020;69(3):59–65. https://doi.org/10.2302/kjm.2019-0009-OA.

    Article  CAS  PubMed  Google Scholar 

  31. Fonfara I, Richter H, Bratovic M, Le Rhun A, Charpentier E. The CRISPR-associated DNA-cleaving enzyme Cpf1 also processes precursor CRISPR RNA. Nature. 2016;532(7600):517–21. https://doi.org/10.1038/nature17945.

    Article  CAS  PubMed  Google Scholar 

  32. Campa CC, Weisbach NR, Santinha AJ, Incarnato D, Platt RJ. Multiplexed genome engineering by Cas12a and CRISPR arrays encoded on single transcripts. Nat Methods. 2019;16(9):887–93. https://doi.org/10.1038/s41592-019-0508-6.

    Article  CAS  PubMed  Google Scholar 

  33. Zetsche B, Heidenreich M, Mohanraju P, Fedorova I, Kneppers J, DeGennaro EM, et al. Multiplex gene editing by CRISPR-Cpf1 using a single crRNA array. Nat Biotechnol. 2017;35(1):31–4. https://doi.org/10.1038/nbt.3737.

    Article  CAS  PubMed  Google Scholar 

  34. Nishimasu H, Ran FA, Hsu PD, Konermann S, Shehata SI, Dohmae N, et al. Crystal structure of Cas9 in complex with guide RNA and target DNA. Cell. 2014;156(5):935–49. https://doi.org/10.1016/j.cell.2014.02.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Chen JS, Ma E, Harrington LB, Da Costa M, Tian X, Palefsky JM, et al. CRISPR-Cas12a target binding unleashes indiscriminate single-stranded DNase activity. Science. 2018;360(6387):436–9. https://doi.org/10.1126/science.aar6245.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Teng F, Guo L, Cui T, Wang XG, Xu K, Gao Q, et al. CDetection: CRISPR-Cas12b-based DNA detection with sub-attomolar sensitivity and single-base specificity. Genome Biol. 2019;20(1):132. https://doi.org/10.1186/s13059-019-1742-z.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kim D, Kim J, Hur JK, Been KW, Yoon SH, Kim JS. Genome-wide analysis reveals specificities of Cpf1 endonucleases in human cells. Nat Biotechnol. 2016;34(8):863–8. https://doi.org/10.1038/nbt.3609.

    Article  CAS  PubMed  Google Scholar 

  38. Kleinstiver BP, Tsai SQ, Prew MS, Nguyen NT, Welch MM, Lopez JM, et al. Genome-wide specificities of CRISPR-Cas Cpf1 nucleases in human cells. Nat Biotechnol. 2016;34(8):869–74. https://doi.org/10.1038/nbt.3620.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Yang M, Wei H, Wang Y, Deng J, Tang Y, Zhou L, et al. Targeted disruption of V600E-mutant BRAF gene by CRISPR-Cpf1. Mol Ther Nucleic Acids. 2017;8:450–8. https://doi.org/10.1016/j.omtn.2017.05.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Safari F, Zare K, Negahdaripour M, Barekati-Mowahed M, Ghasemi Y. CRISPR Cpf1 proteins: structure, function and implications for genome editing. Cell Biosci. 2019;9:36. https://doi.org/10.1186/s13578-019-0298-7.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Li S, Zhang X, Wang W, Guo X, Wu Z, Du W, et al. Expanding the scope of CRISPR/Cpf1-mediated genome editing in rice. Mol Plant. 2018;11(7):995–8. https://doi.org/10.1016/j.molp.2018.03.009.

    Article  CAS  PubMed  Google Scholar 

  42. Zhong Z, Zhang Y, You Q, Tang X, Ren Q, Liu S, et al. Plant genome editing using FnCpf1 and LbCpf1 nucleases at redefined and altered PAM sites. Mol Plant. 2018;11(7):999–1002. https://doi.org/10.1016/j.molp.2018.03.008.

    Article  CAS  PubMed  Google Scholar 

  43. Xu R, Qin R, Li H, Li J, Yang J, Wei P. Enhanced genome editing in rice using single transcript unit CRISPR-LbCpf1 systems. Plant Biotechnol J. 2019;17(3):553–5. https://doi.org/10.1111/pbi.13028.

    Article  PubMed  Google Scholar 

  44. Liu RM, Liang LL, Freed E, Chang H, Oh E, Liu ZY, et al. Synthetic chimeric nucleases function for efficient genome editing. Nat Commun. 2019;10(1):5524. https://doi.org/10.1038/s41467-019-13500-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Toth E, Varga E, Kulcsar PI, Kocsis-Jutka V, Krausz SL, Nyeste A, et al. Improved LbCas12a variants with altered PAM specificities further broaden the genome targeting range of Cas12a nucleases. Nucleic Acids Res. 2020;48(7):3722–33. https://doi.org/10.1093/nar/gkaa110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Schindele P, Puchta H. Engineering CRISPR/LbCas12a for highly efficient, temperature-tolerant plant gene editing. Plant Biotechnol J. 2020;18(5):1118–20. https://doi.org/10.1111/pbi.13275.

    Article  PubMed  Google Scholar 

  47. Port F, Starostecka M, Boutros M. Multiplexed conditional genome editing with Cas12a in Drosophila. Proc Natl Acad Sci U S A. 2020;117(37):22890–9. https://doi.org/10.1073/pnas.2004655117.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kim D, Lim K, Kim DE, Kim JS. Genome-wide specificity of dCpf1 cytidine base editors. Nat Commun. 2020;11(1):4072. https://doi.org/10.1038/s41467-020-17889-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Shin U, Brondani V. Analysis of wild type LbCpf1 protein, and PAM recognition variants, in a cellular context. Front Genet. 2020;11:571591. https://doi.org/10.3389/fgene.2020.571591.

    Article  CAS  PubMed  Google Scholar 

  50. Gier RA, Budinich KA, Evitt NH, Cao Z, Freilich ES, Chen Q, et al. High-performance CRISPR-Cas12a genome editing for combinatorial genetic screening. Nat Commun. 2020;11(1):3455. https://doi.org/10.1038/s41467-020-17209-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. DeWeirdt PC, Sanson KR, Sangree AK, Hegde M, Hanna RE, Feeley MN, et al. Optimization of AsCas12a for combinatorial genetic screens in human cells. Nat Biotechnol. 2021;39(1):94–104. https://doi.org/10.1038/s41587-020-0600-6.

    Article  CAS  PubMed  Google Scholar 

  52. Zhang L, Zuris JA, Viswanathan R, Edelstein JN, Turk R, Thommandru B, et al. AsCas12a ultra nuclease facilitates the rapid generation of therapeutic cell medicines. Nat Commun. 2021;12(1):3908. https://doi.org/10.1038/s41467-021-24017-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Zhou J, Chen P, Wang H, Liu H, Li Y, Zhang Y, et al. Cas12a variants designed for lower genome-wide off-target effect through stringent PAM recognition. Mol Ther. 2022;30(1):244–55. https://doi.org/10.1016/j.ymthe.2021.10.010.

    Article  CAS  PubMed  Google Scholar 

  54. Liu X, Liu X, Zhou C, Lv J, He X, Liu Y, et al. Engineered FnCas12a with enhanced activity through directional evolution in human cells. J Biol Chem. 2021:100394. https://doi.org/10.1016/j.jbc.2021.100394.

  55. Toth E, Czene BC, Kulcsar PI, Krausz SL, Talas A, Nyeste A, et al. Mb- and FnCpf1 nucleases are active in mammalian cells: activities and PAM preferences of four wild-type Cpf1 nucleases and of their altered PAM specificity variants. Nucleic Acids Res. 2018;46(19):10272–85. https://doi.org/10.1093/nar/gky815.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Teng F, Cui T, Feng G, Guo L, Xu K, Gao Q, et al. Repurposing CRISPR-Cas12b for mammalian genome engineering. Cell Discov. 2018;4:63. https://doi.org/10.1038/s41421-018-0069-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Chen P, Zhou J, Wan Y, Liu H, Li Y, Liu Z, et al. A Cas12a ortholog with stringent PAM recognition followed by low off-target editing rates for genome editing. Genome Biol. 2020;21(1):78. https://doi.org/10.1186/s13059-020-01989-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Jacobsen T, Ttofali F, Liao C, Manchalu S, Gray BN, Beisel CL. Characterization of Cas12a nucleases reveals diverse PAM profiles between closely-related orthologs. Nucleic Acids Res. 2020;48(10):5624–38. https://doi.org/10.1093/nar/gkaa272.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Wu Y, Yuan Q, Zhu Y, Gao X, Song J, Yin Z. Improving FnCas12a genome editing by exonuclease fusion. CRISPR J. 2020;3(6):503–11. https://doi.org/10.1089/crispr.2020.0073.

    Article  CAS  PubMed  Google Scholar 

  60. Tran MH, Park H, Nobles CL, Karunadharma P, Pan L, Zhong G, et al. A more efficient CRISPR-Cas12a variant derived from Lachnospiraceae bacterium MA2020. Mol Ther Nucleic Acids. 2021;24:40–53. https://doi.org/10.1016/j.omtn.2021.02.012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Zhu D, Wang J, Yang D, Xi J, Li J. High-throughput profiling of Cas12a orthologues and engineered variants for enhanced genome editing activity. Int J Mol Sci. 2021;22(24):13301. https://doi.org/10.3390/ijms222413301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Zhang Y, Ren Q, Tang X, Liu S, Malzahn AA, Zhou J, et al. Expanding the scope of plant genome engineering with Cas12a orthologs and highly multiplexable editing systems. Nat Commun. 2021;12(1):1944. https://doi.org/10.1038/s41467-021-22330-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Teng F, Li J, Cui T, Xu K, Guo L, Gao Q, et al. Enhanced mammalian genome editing by new Cas12a orthologs with optimized crRNA scaffolds. Genome Biol. 2019;20(1):15. https://doi.org/10.1186/s13059-019-1620-8.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Yamano T, Nishimasu H, Zetsche B, Hirano H, Slaymaker IM, Li Y, et al. Crystal structure of Cpf1 in complex with guide RNA and target DNA. Cell. 2016;165(4):949–62. https://doi.org/10.1016/j.cell.2016.04.003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Huang H, Hu Y, Huang G, Ma S, Feng J, Wang D, et al. Tag-seq: a convenient and scalable method for genome-wide specificity assessment of CRISPR/Cas nucleases. Commun Biol. 2021;4(1):830. https://doi.org/10.1038/s42003-021-02351-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Kim D, Luk K, Wolfe SA, Kim JS. Evaluating and enhancing target specificity of gene-editing nucleases and deaminases. Annu Rev Biochem. 2019;88:191–220. https://doi.org/10.1146/annurev-biochem-013118-111730.

    Article  CAS  PubMed  Google Scholar 

  67. Tak YE, Kleinstiver BP, Nunez JK, Hsu JY, Horng JE, Gong J, et al. Inducible and multiplex gene regulation using CRISPR-Cpf1-based transcription factors. Nat Methods. 2017;14(12):1163–6. https://doi.org/10.1038/nmeth.4483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Zhang X, Wang W, Shan L, Han L, Ma S, Zhang Y, et al. Gene activation in human cells using CRISPR/Cpf1-p300 and CRISPR/Cpf1-SunTag systems. Protein & Cell. 2018;9(4):380–3. https://doi.org/10.1007/s13238-017-0491-6.

    Article  CAS  Google Scholar 

  69. Bratovic M, Fonfara I, Chylinski K, Galvez EJC, Sullivan TJ, Boerno S, et al. Bridge helix arginines play a critical role in Cas9 sensitivity to mismatches. Nat Chem Biol. 2020;16(5):587–95. https://doi.org/10.1038/s41589-020-0490-4.

    Article  CAS  PubMed  Google Scholar 

  70. Dahlman JE, Abudayyeh OO, Joung J, Gootenberg JS, Zhang F, Konermann S. Orthogonal gene knockout and activation with a catalytically active Cas9 nuclease. Nat Biotechnol. 2015;33(11):1159–61. https://doi.org/10.1038/nbt.3390.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Breinig M, Schweitzer AY, Herianto AM, Revia S, Schaefer L, Wendler L, et al. Multiplexed orthogonal genome editing and transcriptional activation by Cas12a. Nat Methods. 2019;16(1):51–4. https://doi.org/10.1038/s41592-018-0262-1.

    Article  CAS  PubMed  Google Scholar 

  72. Singh D, Mallon J, Poddar A, Wang Y, Tippana R, Yang O, et al. Real-time observation of DNA target interrogation and product release by the RNA-guided endonuclease CRISPR Cpf1 (Cas12a). Proc Natl Acad Sci U S A. 2018;115(21):5444–9. https://doi.org/10.1073/pnas.1718686115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Holderfield M, Deuker MM, McCormick F, McMahon M. Targeting RAF kinases for cancer therapy: BRAF-mutated melanoma and beyond. Nat Rev Cancer. 2014;14(7):455–67. https://doi.org/10.1038/nrc3760.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Lu H, Liu S, Zhang G, Kwong LN, Zhu Y, Miller JP, et al. Oncogenic BRAF-mediated melanoma cell invasion. Cell Rep. 2016;15(9):2012–24. https://doi.org/10.1016/j.celrep.2016.04.073.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Wu Z, Zhang Y, Yu H, Pan D, Wang Y, Wang Y, et al. Programmed genome editing by a miniature CRISPR-Cas12f nuclease. Nat Chem Biol. 2021;17(11):1132–8. https://doi.org/10.1038/s41589-021-00868-6.

    Article  CAS  PubMed  Google Scholar 

  76. Kim S, Kim D, Cho SW, Kim J, Kim JS. Highly efficient RNA-guided genome editing in human cells via delivery of purified Cas9 ribonucleoproteins. Genome Res. 2014;24(6):1012–9. https://doi.org/10.1101/gr.171322.113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Dai X, Park JJ, Du Y, Kim HR, Wang G, Errami Y, et al. One-step generation of modular CAR-T cells with AAV-Cpf1. Nat Methods. 2019;16(3):247–54. https://doi.org/10.1038/s41592-019-0329-7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Zhang X, Lv S, Luo Z, Hu Y, Peng X, Lv J, et al. MiniCAFE, a CRISPR/Cas9-based compact and potent transcriptional activator, elicits gene expression in vivo. Nucleic Acids Res. 2021;49(7):4171–85. https://doi.org/10.1093/nar/gkab174.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to all members of the Rong and Lin labs for helpful comments and discussions on the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2017YFA0105001) to Dr. Xiaoyang Zhao, the National Natural Science Foundation of China (82072329) to Dr. Ying Lin and (82070002) and (81872511) to Dr. Zhili Rong, and the Fellowship of China Postdoctoral Science Foundation (2021M701615) to Hongxin Huang, (2021M701634) to Xin Zhang, and (2021M691473) to Shufeng Ma.

Author information

Authors and Affiliations

  1. Dermatology Hospital, Southern Medical University, Guangzhou, 510091, China

    Hongxin Huang, Yongfei Hu, Jiajian Zhou, Weiqi Lv, Tao Huang, Dong Wang & Zhili Rong

  2. Cancer Research Institute, School of Basic Medical Sciences, State Key Laboratory of Organ Failure Research, National Clinical Research Center of Kidney Disease, Key Laboratory of Organ Failure Research (Ministry of Education), Southern Medical University, Guangzhou, 510515, China

    Guanjie Huang, Zhihong Tan, Lin Shan, Xin Zhang, Shufeng Ma, Tao Huang, Yuchen Liu, Ying Lin & Zhili Rong

  3. Department of Bioinformatics, School of Basic Medical Sciences, Southern Medical University, Guangzhou, 510515, China

    Yongfei Hu & Dong Wang

  4. Department of Development, School of Basic Medical Sciences, Southern Medical University, Guangzhou, China

    Xiaoyang Zhao

  5. Experimental Education/Administration Center, School of Basic Medical Sciences, Southern Medical University, Guangzhou, 510515, China

    Ying Lin & Zhili Rong

Authors
  1. Hongxin Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Guanjie Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Zhihong Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yongfei Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Lin Shan
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jiajian Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Xin Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Shufeng Ma
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Weiqi Lv
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Tao Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Yuchen Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Dong Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Xiaoyang Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Ying Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Zhili Rong
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

All authors read and approved the final manuscript. Z.R., Y.L., and H.H. conceived the study, designed the experiments, analyzed the data, and wrote the manuscript. H.H. performed most experiments. G.H., L.S., W.L, S.M., H.T., and X.Z.1 (Xin Zhang) helped H.H. with some experiments. Z.T., Y.H., J.Z., Y.L., and D.W. helped to perform bioinformatics analyses. X.Z.2 (Xiaoyang Zhao) edited the manuscript and provided critical suggestions.

Corresponding authors

Correspondence to Ying Lin or Zhili Rong.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors read and approved the final manuscript.

Competing interests

Z.R., Y.L., and H.H. have filed a patent application on Cas12a-Plus variants. The other authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary Information, including supplementary tables and DNA and amino acid sequences. Table S1. The engineering of the CRISPR-AsCas12a/LbCas12a systems. Table S2. The engineering of the CRISPR-Cas12a (excluded AsCas12a and LbCas12a) systems. Table S3. The sgRNAs used in this study. Table S4. Primers for this study (including qPCR, T7E1/PAGE, Deep-seq). Table S5. Index information for deep-seq. Supplementary sequences. DNA and Amino acid sequences used in this study.

Additional file 2: Figure S1.

Construction of AsCas12a variants and assessment of their specificity. Figure S2. Specificity assessment of AsCas12a-RKA. Figure S3. Specificity comparison of AsCas12a-WT, -RKA, -RKK, and -HF by Tag-seq. Figure S4. The Engineering of crRNA to improve efficiency. Figure S5. Editing of BRAFV600E with SpCas9 and AsCas12a nucleases. Figure S6. Detection the expression levels of LbCas12a nucleases. Figure S7. Specificity assessment of AsCas12a-WT, -Plus, -HF, enAsCas12a, -ultra, -UKA, -UKK, -RU, -RKAU, and -RKKU by Tag-seq. Figure S8. Specificity assessment of LbCas12a-WT, -Plus, -ultra, -UKA, -UKK, -RU, -RKAU, -RKKU by Tag-seq. Figure S9. Editing efficiency analyses of AsCas12a-WT, -Plus, -HF, enAsCas12a, -ultra, -UKA, -UKK, and -RU. Figure S10. Editing efficiency analyses of LbCas12a-WT, -Plus, -ultra, -UKA, -UKK, -RU, -RKAU, and -RKKU.

Additional file 3.

The supporting data values for the figures, including Figs. 1e, f, 2c, e, f, 3d, e, 4b, e, 5b, 6e, f, 7d, e, g, h, and Figs. S1c, S4b, c, S9, S10.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visithttp://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Huang, H., Huang, G., Tan, Z. et al. Engineered Cas12a-Plus nuclease enables gene editing with enhanced activity and specificity. BMC Biol 20, 91 (2022). https://doi.org/10.1186/s12915-022-01296-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12915-022-01296-1

Keywords

BMC Biology

ISSN: 1741-7007

Contact us