- Research
- Open access
- Published:
Creating large-scale genetic diversity in Arabidopsis via base editing-mediated deep artificial evolution
Genome Biology volume 25, Article number: 215 (2024)
Abstract
Background
Base editing is a powerful tool for artificial evolution to create allelic diversity and improve agronomic traits. However, the great evolutionary potential for every sgRNA target has been overlooked. And there is currently no high-throughput method for generating and characterizing as many changes in a single target as possible based on large mutant pools to permit rapid gene directed evolution in plants.
Results
In this study, we establish an efficient germline-specific evolution system to screen beneficial alleles in Arabidopsis which could be applied for crop improvement. This system is based on a strong egg cell-specific cytosine base editor and the large seed production of Arabidopsis, which enables each T1 plant with unedited wild type alleles to produce thousands of independent T2 mutant lines. It has the ability of creating a wide range of mutant lines, including those containing atypical base substitutions, and as well providing a space- and labor-saving way to store and screen the resulting mutant libraries. Using this system, we efficiently generate herbicide-resistant EPSPS, ALS, and HPPD variants that could be used in crop breeding.
Conclusions
Here, we demonstrate the significant potential of base editing-mediated artificial evolution for each sgRNA target and devised an efficient system for conducting deep evolution to harness this potential.
Background
Genetic diversity, e.g., variation in coding sequences, non-coding genes, and expression-regulating elements, is the basis for plant domestication and crop breeding. Natural variation and random mutations induced by physical and chemical treatments are the sources of genetic diversity in conventional breeding [1]. However, with the growing global population and increasing food demand, there is an urgent need for even greater genetic diversity in modern agriculture. Fortunately, genome editing tools offer a promising solution for protein evolution [1]. There are two major ways to achieve directed evolution, insertions/deletions (indels), and base changes [2,3,4,5]. Indels normally disrupt the function of genes and regulatory elements. Base changes generate a wider range of genetic variants than indels, including gain-of-function changes that fine-tune gene function. Prime editor (PE) and base editor (BE) are two major editing tools for introducing targeted base substitutions [6].
PE is a “search-and-replace” technology [7]. It has been used to create a large number of base substitutions at six conserved residues of acetyl-coenzyme A carboxylase 1 (OsACC1), generating 16 herbicide-resistant alleles in rice [5]. However, the editing products of PE are determined by the editing templates in the pegRNAs, and one pegRNA can only generate one editing product [7]. Hence, PE is suitable for modifying a defined base at a critical position rather than generating large mutant pools of altered bases for high-throughput screening for putative critical positions.
Unlike PE, BE can generate multiple types of editing products at each sgRNA target. Thus, BE has shown greater potential than PE for generating ample base variations. Previously, BE-mediated gene artificial evolution focused on designing multiple sgRNAs to cover whole gene or key regions of these genes [2]. However, the enormous evolutionary potential at every single sgRNA site has been overlooked. For example, cytosine base editor (CBE) is capable of converting multiple cytosines within the editing window in a single mutant derivative [8]. In addition to the primary C-to-T conversion, CBE produces atypical base substitutions such as C-to-G and C-to-A, albeit at lower frequencies [9]. As a result, it can theoretically generate a vast number of alleles, with many variants appearing rarely in nature. Previous studies have obtained herbicide-resistant rice by designing multiple sgRNA targets to cover key positions of genes such as acetolactate synthase (ALS) using base editors [2], but this approach requires obtaining a large number of rice mutant plants, which is time and labor-consuming and low-throughput, and this method has a major penalty that the number of nucleotide and protein mutation forms obtained by every single target is thin, which limits the possibilities of evolution. The generation of large-scale genetic variation depends on having a sufficient number of mutant lines, which is hard to create in crops. Moreover, storage and phenotypic screening of large-scale mutant libraries require much space, time, and labor. So, despite the rapid development of genome edit toolbox, there are no such systems capable of large-scale targeted gene mutagenesis and high-throughput screening in crops.
Arabidopsis (Arabidopsis thaliana) is a model plant that is small, has a short life cycle, and produces thousands of small seeds per plant. We can conduct phenotypic screening on young Arabidopsis seedlings in Petri dishes or soil. Furthermore, it shares many essential genes with crop plants, e.g., genes regulating stress tolerance, herbicide resistance, and photosynthetic efficiency [10, 11]. It would be a breakthrough to generate sufficient mutant derivatives of each sgRNA target required for high-throughput directed evolution in Arabidopsis to produce beneficial allelic forms for crop trait improvement. Here, we developed an efficient pipeline for germline-specific evolution that can easily produce and screen thousands to millions of mutant lines. Our studies of large mutant population produced in this way achieved significant improvements in directed evolution of genes and offered genetic resources for herbicide resistant trait and more importantly ensuing successful transmission application in crop plants. We anticipate that this approach will speed up many aspects of directed evolution for crop improvement.
Results
Development of the germline-specific evolution system in Arabidopsis
Base editor has great potential for generating ample base variations. Due to the advantages we mentioned above, Arabidopsis might be an excellent candidate for base editor mediated directed evolution in situ owing to its advantages of high-throughput. However, the diversity of mutation generated by base editor in Arabidopsis still needs to be improved. Thus, we developed a germline-specific evolution system, which combined an efficient germline-specific base editor with next generation screening, aimed to break this hurdle and generate diverse mutant forms during sexual propagation, enabling high-throughput screening for beneficial derivatives which could be further used for crop improvement (Fig. 1a). Base substitutions in some alleles often fail to be generated in the T1 plants probably due to the short time expression after floral dip. However, the same wild-type (WT) alleles continue to be subjected to editing in the germline cells of T1 plants during sexual propagation, thereby generating a variety of base edited lines in the T2 generation (Additional file 1: Fig. S1). A single T1 Arabidopsis plant with an unchanged WT allele has the potential to produce thousands of seeds, approximately 75% or more of which contain T-DNA expressing the base editor. Thus, this system enables a single T1 plant with WT allele to generate thousands of independent T2 mutant lines. Compared with traditional methods that screen the T1 plants, the system built here has the potential to increase the number of mutant lines by several orders of magnitude.
The widely used CBE vector in Arabidopsis, pHEE901, contains the rat APOBEC1-based BE3 (base editor 3) under the control of the EC1.2en-EC1.1p promoter [12]. It has been reported that unedited T1 Arabidopsis plants containing this vector can generated base edited T2 mutants. However, the editing efficiency is less than 10%, and no C-to-A/G mutant has been identified which significantly restricts the diversity of potential allele forms [12]. To address this problem, we replaced rat APOBEC1 with a more efficient cytosine deaminase, human APOBEC3A [13], to create pHEE901-A3A (Fig. 1b). To check whether the editing efficiency of pHEE901-A3A was increased relative to that of the original BE3, we designed five sgRNAs for four target genes including agamous (AG) [14], auxin response factor 2 (ARF2) [15], 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) [16], and hypersensitive to ABA 1 (HAB1) [17] (Additional file 1: Figs. S2, S3 and S4). A total of 817 transgenic lines were generated, and Hi-TOM sequencing was used to assess the genotypes of all T1 plants [18]. pHEE901-A3A proved to be more efficient than pHEE901, with a wide editing window (Fig. 1c, Additional file 2: Tables S1 and S2). Notably, only pHEE901-A3A produced C-to-A and C-to-G conversions (Additional file 2: Table S2).
Although pHEE901-A3A performed base editing efficiently in T1 plants, 50.0–92.3% of these plants were WT, and 2.6–14.0% were heterozygotes containing an intact WT allele (Additional file 2: Table S1). To evaluate the editing efficiency of pHEE901-A3A during the sexual reproduction process from T1 to T2, we randomly selected six T1 WT plants containing sgRNA targeting AG and six targeting EPSPS. We then assessed the genotype of approximately 300 T2 plants from each of the T1 plants. Base editing had occurred in 31.2–76.9% of the T2 plants with sgRNA targeting AG and 14.9–45.8% of the T2 plants with sgRNA targeting EPSPS (Fig. 1d and Additional file 2: Table S3). These results confirmed that thousands of T2 mutant lines could be generated from a single T1 WT plant using this germline-specific evolution system.
Diverse T2 mutants can be generated at each sgRNA target using the germline-specific evolution system
To maximize the mutation depth of every single target, we investigated whether the T2 mutants were more extensive and diverse genetic variations than the T1 mutants at each target. Both protospacers of the AG and EPSPS targets contained seven cytosines (Additional file 1: Fig. S2). In theory, all seven cytosines could be altered; we could expect to obtain lots of mutated allele forms for each target. In the event, we respectively identified 106 and 36 forms of mutated allele for AG and EPSPS in the tested 1800 T2 plants, compared with only 22 and 12 in the 88 and 85 T1 plants (Fig. 2a, b, Additional file 1: Figs. S4 and S5). Our findings provide strong evidences supporting a correlation between the scale of mutants and types of genetic variants. More importantly, in this experiment, we have only tested a very small number of T2 plants, equivalent to only approximately a quarter of seeds produced by one T1 plant, whereas if we characterized and screened the seeds produced by all the T1 transgenic plants, we would obtain significantly larger number of mutation forms.
Analysis of the profiles of the editing products revealed that our system expanded the size of the editing window (Fig. 2c, Additional file 1: Figs. S4 and S5). For instance, the editing window of the EPSPS target expanded to 34 bp (from position − 18 to 16) in the T2 generation compared to 7 bp (positions 5 to 11) in the T1 plants (Fig. 2c). Similarly, base substitutions outside of the protospacer of the AG target were detected only in the T2 plants (Fig. 2c). We also observed that our evolution system yielded a wide range of atypical base substitutions including not only C-to-G and C-to-A but also A-to-T, G-to-A, G-to-C, and T-to-C, suggesting that deamination also occurred on A and G (Fig. 2c). The T-to-C conversion might have been resulted from deamination of A on the complementary strand. These results underscore the significant evolutionary potential of base editing on single sgRNA target.
Although more than half of the T2 altered allele forms were the results of atypical base substitutions, the most frequent deamination events occurred on cytosines in a narrow window, and the most common change was a C-to-T transition (Fig. 2d and Additional file 1: Fig. S5). As a result, most of the altered allele forms, including those resulting from C-to-T conversions in the extended editing window, or involving atypical base substitutions, arose at relatively low frequencies: over 75% occurred at frequencies below 1% and more than 40% at frequencies less than 0.1% (only one plant harboring the corresponding allele in the set of approximately 1800 T2 plants) (Fig. 2e and Additional file 1: Fig. S5). These findings show that our system provides an effective strategy for generating large-scale mutant libraries at every single target needed for achieving high-throughput directed evolution and screening for rare mutant forms. The evolutionary depth and efficiency of atypical base substitution can be further improved by using more potent deaminases and designing multiple target sites.
The germline-specific evolution system facilitates directed evolution of proteins
Since our system increased variation at the DNA level, we expected that it would also generate a wider spectrum of variants of AG and EPSPS at the protein level. Although we observed only 2 types of full-length AG protein variant and 10 types of full-length EPSPS protein variant in T1 plants, we identified 15 and 25 forms of full-length protein, respectively, in the tested T2 plants (Fig. 3a, b).
Next, we analyzed variation at the amino acid level and found that the presence of atypical base substitutions led to extensive amino acid variation. For example, P182 (CCA) of EPSPS was converted not only to Ser (TCA) and Leu (TTA and CTA) but also to Ala (GCA), Val (GTA), Ile (ATA), Thr (ACA), and a stop codon (TGA) (Fig. 3c, Additional file 1: Figs. S5 and S6). T2 plants contained a wider variety of amino acid changes at T152, R153, I154, R155, and S156 of AG and R181 and L183 of EPSPS than T1 plants (Fig. 3c and Additional file 1: Fig. S6). Additionally, due to the broader editing window, amino acid changes at L147, E148, and R149 of AG and G177 and A179 of EPSPS were detected only in T2 plants (Fig. 3c).
Premature termination of protein translation, which often results from C-to-T conversions at CAA, CGA, and CAG codons, poses a significant challenge to CBE-mediated directed evolution of proteins, as it can impact the evolution of not only these specific codons but also adjacent ones. For instance, the protospacer region of the AG target comprises two CGA codons encoding R153 and R155, with cytosines situated at positions 5 and 11 (Fig. 2c). In the T1 generation, only two of 42 mutant lines contained full-length mutated proteins with single amino acid changes (R153G or S156F), while the others resulted in premature termination, which prevented the functional detection of variants of the other amino acids in the protospacer region (Fig. 3d and Additional file 1: Fig. S6). Screening a large T2 mutant library circumvented this limitation: 102 of 1080 T2 mutated alleles encoded 15 different forms of full-length mutated AG protein were created, including amino acid conversions at five different positions (Additional file 1: Fig. S6). We noted that in these full-length AG variants, CGA either underwent atypical C-to-G/A or G-to-A/C changes instead of C-to-T changes or base editing occurred specifically on the surrounding codons but not CGA (Fig. 3e and Additional file 1: Fig. S7). These results showed that our evolution system is well protected against the adverse effects of mutant protein premature stoppage on artificial evolution. Similar to the frequency at the DNA level, over 60% of the different protein variants appeared at a frequency below 1% and more than 40% at a frequency less than 0.1% (Fig. 3f). In sum, the diversity of mutation significantly increased by using this germline-specific evolution system which could enable comprehensive and large-scale directed evolution of proteins.
Screening for EPSPS mutants exhibiting herbicide-resistance
Glyphosate is one of the most widely used herbicides that target EPSPS [16]. As the targeting site of EPSPS sgRNA used above is in close proximity to the coding sequence of key amino acids related to glyphosate resistance, we assessed the potential to generate herbicide resistant mutants. We sowed approximately 50,000 T2 seedlings from 10 T1 WT plants containing this sgRNA to 1/2 MS medium supplemented with 100 mM glyphosate (Fig. 4a). Only three seedlings grew, and all of them produced the same form of EPSPS with a rare combination of amino acid substitutions: T178I/A179V/P182S (Fig. 4b). This variant was the result of three C-to-T substitutions on C-4, C-1, and C8 and was not encountered in the previously sequenced T2 plants (Fig. 4b and Additional file 1: Fig. S5). Moreover, C-4 was not changed in any of the previously sequenced T2 plants (Fig. 2c). These results demonstrate that a large mutant library is critical for isolating rare beneficial mutants. Interestingly, the T178I/A179V/P182S form of EPSPS has been found in a naturally occurring high glyphosate-resistant strain of herbicide-resistant weed, and its equivalent variants have been used in rice and maize breeding [16, 19, 20].
In the above screening, 10 unedited T1 WT plants generated more than 50,000 T2 seeds that could be easily accommodated in 10 Eppendorf tubes. According to the mutation rate calculated above, these T2 seeds should contain approximately 15,000 independent mutant lines, a number that is almost impossible to obtain in other ways. Moreover, we needed only 25 petri dishes to screen for herbicide-resistance for the 50,000 seeds. Thus, the use of the germline-specific evolution system in Arabidopsis is space/labor-saving and effective in generating, storing and high-throughput screening of mutant libraries.
Screening for herbicide-resistant ALS mutants
We also conducted directed evolution of Arabidopsis ALS, the target of five groups of commercial herbicides: imidazolinones (IMI), pyrimidinylthiobenzoates (PTB), sulfonylaminocarbonyltriazolinones (SCT), sulfonylureas (SU), and triazolopyrimidines (TP) [21]. The development of ALS mutants with broad-spectrum resistance (BSR) to several of these herbicides could benefit agricultural breeding focused on weed control. To isolate such a mutant, we designed 12 sgRNAs targeting the coding sequences of conserved ALS regions between Arabidopsis and rice and then created more than 60 transgenic T1 plants for each sgRNA (Fig. 4c, Additional file 1: Figs. S8 and S9). To screen for herbicide-resistant mutants, we grew 120 T2 seedlings from each T1 plant on 1/2 MS plates containing chlorsulfuron (SU herbicide), bispyribac (PTB herbicide), or imazapic (IMI herbicide). In total, we screened over 100,000 T2 seedlings from 834 T1 plants, determined the genotypes of 1244 herbicide-resistant T2 plants by Hi-TOM sequencing, and identified 46 forms of ALS mutant, including 32 that had not been reported previously (Fig. 4d, e, Additional file 1: Fig. S10 and Additional file 2: Table S4).
We focused on alleles with multiple resistance and found that one T2 seedling with a P197V/R198C/R199C mutation exhibited tolerance to imazapic; meanwhile, P197V and P197V/R198C mutations were observed in five chlorsulfuron-resistant seedlings (Additional file 2: Table S4). All P197F, P197F/R198C and P197F/R198C/R199C were resistance to three herbicides, moreover, P197F/R198C and P197F/R198C/R199C mutants showed resistance to all three tested herbicides (Additional file 2: Table S4). To confirm their resistances to different ALS-targeting herbicides, we generated T3 homozygous derivatives. None of these mutations had a significant influence on plant growth in the absence of herbicide treatment (Additional file 1: Fig. S11). We then grew the mutants on 1/2 MS plates with each of five ALS-targeting herbicides (Fig. 4f). P197V offered robust resistance to four tested herbicides, namely imazapic, pyroxsulam, chlorsulfuron, and flucarbazone, whereas P197F showed strong resistance only to bispyribac and chlorsulfuron (Fig. 4f and Additional file 1: Fig. S12). Moreover, the R198C/R199C mutation increased the tolerance of P197F/V to higher doses of several herbicides (Fig. 4f and Additional file 1: Fig. S12).
To check whether P171F and P197V/R198C/R199C could be used as BSR forms of ALS in crops, we generated genome edited rice plants with modified gene sequences encoding the corresponding P171F and P171V/R172C/R173C variants of OsALS1 using prime editing (Additional file 1: Fig. S13). The BSR of both variants in rice was similar to that observed in Arabidopsis, indicating that the germline-specific evolution system works well for crop improvement (Fig. 4g and Additional file 1: Fig. S14). These results demonstrate that the germline-specific evolution system can perform efficient and effective directed evolution of prominent agronomic trait genes, and more importantly, the obtained mutations can be good applied in molecular breeding of other crops.
Screening for herbicide-resistant HPPD mutants
To conduct direct evolution on 4-hydroxyphenylpyruvate dioxygenase (HPPD), another conserved protein and important target of commercial herbicides (Additional file 1: Fig. S15) [22], we generated a Cas9-NG vector [23], pHEE901-A3A-NG, to expand the target scope of our system, and designed 92 sgRNAs to cover its coding sequence (Fig. 4h and Additional file 2: Table S5). Over 3 million T2 seedlings were planted on 1/2 MS medium supplemented with five HPPD inhibitor herbicides, including pyrasulfotole, topramezone, mesotrione, tembotrione, and isoxaflutole (Fig. 4i). We obtained 26 herbicide-resistant HPPD mutation forms caused by 10 sgRNAs (Fig. 4h, Additional file 1: Figs. S16 and S17). These mutant forms contained substitutions on 17 highly conserved amino acid residues, particularly among dicot plants (Additional file 1: Figs. S15 and S16). Among the 26 mutant forms, 14 occurred on the target of sgRNA68, and all of them contained amino acid substitutions on P339, suggesting that the critical roles of P339 and its surrounding amino acids in herbicide binding. In a previous study, six amino acids, including P339, were identified as potentially contributory for the divergence in herbicide-binding activity between AtHPPD and ZmHPPD [24]. Our findings indicate that the diversity on this amino acid may partially account for the HPPD-type herbicide resistance observed in maize. Further testing in T2 generation indicated that mutants on the target of sgRNA68 and sgRNA73 showed higher degree of resistance to these herbicides (Fig. 4j). However, all of these mutants have not provided enough tolerance which can be applied in crop plants. In these assays, we have identified several important regions of HPPD strongly relating to herbicide resistance. Thus, in the future, we will simultaneously express several sgRNAs to cover all of these key regions to perform larger scale evolution for creating high level BSR mutants, which can be applied in crop plants.
Discussion
In this study, we describe the potential of BE-mediated artificial evolution for generating a large allelic diversity at single target which has been overlooked before, probably due to the lack of an efficient tool for producing and screening large mutant libraries. To address this issue, we established a germline-specific evolution system to screen beneficial alleles in Arabidopsis which could be applied for crop improvement. This system allows unedited T1 Arabidopsis plants with germline-specifically expressed A3A-BE to produce thousands of T2 mutant lines, all of which can be stored in a small Eppendorf tube, enabling space-saving and labor-efficient phenotypic screening.
It has been demonstrated that A3A-BE behave broader editing window than APOBEC1-mediated BE due to its higher editing efficiency [13]. Rare base substitutions, which occur at very low frequencies, can only been found in large-scale mutant pool. Thus, the use of this germline-specific evolution system broadens the editing window and permits the detection of atypical base conversions, significantly increasing the available genetic diversity at both the DNA and protein levels. Importantly, this system reduces the inhibitory effect of premature translation termination on the directed evolution of proteins. We believe that the germline-specific evolution system should also permit us to co-evolve key amino acid residues distributed across different regions of a protein.
Several in vivo evolution studies have been conducted in plants, with screenings primarily focusing on crops to identify herbicide-resistant alleles [2, 25,26,27,28,29]. These screenings typically involved selecting calli on plates supplemented with herbicide, from which T0 resistant lines were subsequently generated. However, the labor-intensive and space-demanding nature of callus transformation limited the selection to only a few thousand calli, thus restricting the variety of mutants that could be examined. In contrast, our system conducted screening in the T2 generation of Arabidopsis, facilitating the evaluation of millions of transgenic lines for the same purpose, resulting in deep evolution on each sgRNA and higher yield of herbicide-resistant mutant alleles. For instance, in comparison to the study by Kuang et al., who employed 63 sgRNAs to evolve ALS in rice and identified only five types of herbicide-resistant variants, with P171F being the most potent [2], we discovered 46 variants. Among these, P171F/R172C/R173C and P171V/R172C/R173C exhibited stronger and broader resistance to ALS-targeting herbicides than P171F.
Multiple strategies could be applied to further increase the genetic diversity created by the germline-specific evolution system. For instance, dCas9-AID(x) has shown the capability to induce more C-to-G/A substitutions compared to dCas9-AID(x)-UGI [30]. Therefore, by modifying A3A-CBE through the removal of the UGI or replacing UGI with uracil-DNA glycosylase [31], we could potentially amplify the mutation forms. Furthermore, exploring the substitution of A3A-CBE with different glycosylase base editors to facilitate C-to-N, T-to-N, A-to-N, and G-to-N substitutions is another promising avenue [32, 33]. Additionally, the utilization of adenine and cytosine base editor (ACBE) could enable the simultaneous base substitution of both A and C nucleotides [25, 34,35,36,37,38,39]. These prospective strategies present exciting directions for our research. The use of BE-mediated evolution system is complementary to that of PE, which is well-suited to achieving directed changes at critical positions in a gene rather than high-throughput screening at the whole-gene level. Using BE, our system can identify critical positions for focused PE-mediated directed evolution.
In addition to studying resistance traits, Arabidopsis is an excellent model plant for investigating various aspects such as photosynthesis rate, fertilizer use efficiency, disease resistance, abiotic tolerance, and heavy metal absorption. Many key genes that determine these traits are highly conserved between Arabidopsis and crop plants. The germline-specific evolution system enables the generation of a diverse library of mutant lines for screening purposes. While some mutations can be easily screened in petri dishes by naked eye, certain phenotypic screenings may require more extensive work and equipment compared to screening for herbicide resistance. For instance, stresses such as salt, temperature fluctuations, pathogens, and various other stressors can lead to an increase in cytoplasmic Ca2+ levels and a decline in photosystem II activity [40,41,42]. In such cases, we can screen for mutants affecting plant responses to these stresses by utilizing an aequorin-expressing background [43, 44], which detects Ca2+ signals and exhibits fluorescence, or screen for mutants with altered photosystem II activity using Imaging-PAM [45]. Anyway, the workload for screening in Arabidopsis is typically smaller than that in crop plants. Furthermore, there is potential to leverage phenomics tools in combination with image recognition using machine learning to establish a series of efficient high-throughput phenotypic screening systems.
Conclusions
Base editing holds significant promise for artificial evolution, but so far, its full potential has been underestimated. To unlock this potential, the availability of expansive mutant pools with diverse genetic profiles and reliable high-throughput screening systems is crucial. The germline-specific evolution system we established in this study has demonstrated the capacity to generate a wide range of genetic variations by introducing unconventional base substitutions and expanding editing window and to perform large scale screening in a space/labor-saving way. This system has successfully addressed the challenge of premature termination, a significant hurdle in CBE-mediated protein evolution. Notably, the application of this system has enabled the efficient production of herbicide-resistant variants, showcasing its potential benefits for crop breeding.
Methods
Plant material and growth conditions
Arabidopsis ecotype Columbia (Col-0) was used for transformation and phenotypic analysis. All Arabidopsis plants were grown in a growth chamber with a 22 °C, 16/8 h day/night photoperiod, and 70% relative humidity. Rice variety Zhonghua11 was used for Agrobacterium-mediated transformation of rice callus cells, as previously described [46]. Hygromycin at 50 μg/mL was used to select transgenic calli, and transgenic plantlets were regenerated on selective medium 10 ~ 12 weeks later. The transgenic plantlets were grown at 28 °C and 16/8 h day/night photoperiod in bottles.
Plasmid construction
To generate plasmid pHEE901-A3A, pHEE901 was digested with XbaI and SacI to generate a 12.8 kb fragment, and pH-A3A-PBE plasmid was digested with AvrII and SacI, and a 5.1-kb fragment was recovered [13]. Then, the two fragments were ligated using T4 DNA ligase to produce pHEE901-A3A.
To generate the pHEE901-A3A-NG plasmid, we firstly digested the pHEE901-A3A using MluI and SacI. Next, we ligated this plasmid backbone with a 1.2-kb fragment that was digested from pHUE411-NG using the same enzymes to obtain an intermedia plasmid lacking uracil glycosylase inhibitor (UGI). Then, the UGI sequence was amplified from pHEE901-A3A using the primer pair UGI-F and UGI-R and inserted into the SacI-digested intermedia plasmid by Gibson assembly to generate the final pHEE901-A3A-NG plasmid.
To obtain Arabidopsis transgenic lines, 109 pairs of oligonucleotides encoding the designed sgRNAs were synthesized, annealed, and cloned into BsaI-digested pHEE901, pHEE901-A3A, or spCas9-NG by T4 DNA ligase. The sequences of the synthesized DNA oligos and all primers used in this study are listed in Additional file 2: Table S5.
For the generation of genome edited rice ALS mutants, the diagram and sequence of pegRNA and nicking sgRNA are listed in Additional file 1: Fig. S18. These sequences were synthesized and digested with BsaI and then inserted into the BsaI digested Prime editing vector enPPE2 [47].
Agrobacterium-mediated transformation of Arabidopsis
Constructs were transformed into Agrobacterium tumefaciens strain GV3101 and introduced into Arabidopsis by the floral dip method as described [48]. Seeds were harvested one-and-a-half months later and selected on 1/2 MS plates containing 25 μg/mL hygromycin. Resistant seedlings were transplanted to soil 2 weeks later.
Agrobacterium-mediated transformation of rice callus cells
Constructs were transformed into Agrobacterium tumefaciens strain EHA105 by electroporation. Agrobacterium-mediated transformation of Zhonghua11 callus cells was carried out as reported [46]. Three days after cocultivation, the calluses were moved to medium with 50 μg/mL hygromycin to select positive transformants.
Assessment of base editing efficiency
To measure the editing efficiency of T1 plants, samples of WT and T1 plants were collected. The editing efficiency of T2 plants was measured using samples collected from the progenies of WT T1 lines. DNA was extracted by heating samples at 95 °C for 10 min in DNA extraction buffer (50 mM Tris–HCl pH 7.5, 300 mM NaCl, 300 mM sucrose). After centrifugation at 12,000 g for 5 min, supernatants were collected and used as templates for PCR using SanTaq PCR Mix (Sangon Biotech). PCR products were sequenced by Sanger sequencing or Hi-TOM sequencing [18]. The sequencing results were analyzed using the BioEdit software.
The mutation type (homozygous, heterozygous, biallelic or chimeric) of each T1 and T2 plant was determined by Hi-TOM sequencing. Reads that represented less than 10% of the results of Hi-TOM sequencing were considered aerosol contaminants or technical errors. Plants that contained one allele were identified as WT or homozygote; plants had two alleles with within-twice frequency differences were considered to have two equal alleles, i.e., to be heterozygous or biallelic; otherwise, they were considered chimeric.
The frequencies of DNA alleles, protein variants, and single base or amino acid changes were calculated as number of plants with these changes versus the total number of plants examined.
Selection of herbicide-resistant plants
To screen the herbicide-resistant EPSPS mutants, seeds of T1 plants were sown on 1/2 MS plates supplemented with 100 mM glyphosate. To select the herbicide-resistant ALS mutants, seeds of T1 plants were planted on 1/2 MS plates containing 0.02 mg/L chlorsulfuron, 0.03 mg/L bispyribac, or 0.05 mg/L imazapic. To obtain the herbicide-resistant HPPD mutants, seeds of T1 plants were sown on 1/2 MS plates supplemented with 0.75 mg/L pyrasulfotole, 0.1 mg/L topramezone, 0.1 mg/L mesotrione, 0.1 mg/L tembotrione, or 0.04 mg/L isoxaflutole. Resistant seedlings were transplanted to soil two weeks later for DNA extraction and seed collection. Hi-TOM sequencing was used to determine the mutant alleles in herbicide-resistant plants.
BSR assay for ALS mutants
Seeds of the P197F, P197F/R198C, P197F/R198C/R199C, P197V, P197V/R198C, and P197V/R198C/R199C mutants were sown on 1/2 MS plates containing different concentrations of herbicides (low dosages: 0.04 mg/L imazapic, 0.04 mg/L pyroxsulam, 0.03 mg/L bispyribac, 0.3 mg/L chlorsulfuron or 0.3 mg/L flucarbazone-sodium; high dosages: 0.08 mg/L imazapic, 0.08 mg/L pyroxsulam, 0.06 mg/L bispyribac, 0.6 mg/L chlorsulfuron or 0.6 mg/L flucarbazone-sodium) for the BSR assay. After 2 weeks, photographs were taken, and the fresh weights of individual plants were recorded.
To investigate whether P171V/R172C/R173C, the equivalent OsALS mutant to P197V/R198C/R199C of AtALS, conferred BSR to ALS herbicides in rice, genome edited plants with modified gene sequences encoding the corresponding P171F and P171V/R172C/R173C variants of OsALS1-using prime editor in T1 rice plants were grown on 1/2 MS medium for 2 days and transferred to 1/2 MS solid medium containing 0.9 mg/L imazapic, 0.4 mg/L pyroxsulam, 1.68 mg/L bispyribac, 8 mg/L chlorsulfuron, or 3.15 mg/L flucarbazone-sodium for 16 days to check for growth.
BSR assay for HPPD mutants
Seeds of T2 of HPPD mutants were sown on 1/2 MS plates supplemented with 0.2 mg/L pyrasulfotole, 0.03 mg/L topramezone, 0.02 mg/L isoxaflutole, 0.03 mg/L mesotrione, or 0.03 mg/L tembotrione for 9 days.
Availability of data and materials
All study data are included in the article and supporting information.
References
Chen K, Wang Y, Zhang R, Zhang H, Gao C. CRISPR/Cas genome editing and precision plant breeding in agriculture. Annu Rev Plant Biol. 2019;70:667–97.
Kuang Y, Li S, Ren B, Yan F, Spetz C, Li X, Zhou X, Zhou H. Base-editing-mediated artificial evolution of OsALS1 in planta to develop novel herbicide-tolerant rice germplasms. Mol Plant. 2020;13:565–72.
Rodriguez-Leal D, Lemmon ZH, Man J, Bartlett ME, Lippman ZB. Engineering quantitative trait variation for crop improvement by genome editing. Cell. 2017;171(470–480): e478.
Song X, Meng X, Guo H, Cheng Q, Jing Y, Chen M, Liu G, Wang B, Wang Y, Li J, Yu H. Targeting a gene regulatory element enhances rice grain yield by decoupling panicle number and size. Nat Biotechnol. 2022;40:1403–11.
Xu R, Liu X, Li J, Qin R, Wei P. Identification of herbicide resistance OsACC1 mutations via in planta prime-editing-library screening in rice. Nat Plants. 2021;7:888–92.
Anzalone AV, Koblan LW, Liu DR. Genome editing with CRISPR-Cas nucleases, base editors, transposases and prime editors. Nat Biotechnol. 2020;38:824–44.
Anzalone AV, Randolph PB, Davis JR, Sousa AA, Koblan LW, Levy JM, Chen PJ, Wilson C, Newby GA, Raguram A, Liu DR. Search-and-replace genome editing without double-strand breaks or donor DNA. Nature. 2019;576:149–57.
Komor AC, Kim YB, Packer MS, Zuris JA, Liu DR. Programmable editing of a target base in genomic DNA without double-stranded DNA cleavage. Nature. 2016;533:420–4.
Rees HA, Liu DR. Base editing: precision chemistry on the genome and transcriptome of living cells. Nat Rev Genet. 2018;19:770–88.
Koornneef M, Meinke D. The development of Arabidopsis as a model plant. Plant J. 2010;61:909–21.
Somerville C, Koornneef M. A fortunate choice: the history of Arabidopsis as a model plant. Nat Rev Genet. 2002;3:883–9.
Chen Y, Wang Z, Ni H, Xu Y, Chen Q, Jiang L. CRISPR/Cas9-mediated base-editing system efficiently generates gain-of-function mutations in Arabidopsis. Sci China Life Sci. 2017;60:520–3.
Zong Y, Song Q, Li C, Jin S, Zhang D, Wang Y, Qiu JL, Gao C: Efficient C-to-T base editing in plants using a fusion of nCas9 and human APOBEC3A. Nat Biotechnol. 2018;36:950–3.
Wu HW, Deng S, Xu H, Mao HZ, Liu J, Niu QW, Wang H, Chua NH. A noncoding RNA transcribed from the AGAMOUS (AG) second intron binds to CURLY LEAF and represses AG expression in leaves. New Phytol. 2018;219:1480–91.
Schruff MC, Spielman M, Tiwari S, Adams S, Fenby N, Scott RJ. The AUXIN RESPONSE FACTOR 2 gene of Arabidopsis links auxin signalling, cell division, and the size of seeds and other organs. Development. 2006;133:251–61.
Perotti VE, Larran AS, Palmieri VE, Martinatto AK, Alvarez CE, Tuesca D, Permingeat HR. A novel triple amino acid substitution in the EPSPS found in a high-level glyphosate-resistant Amaranthus hybridus population from Argentina. Pest Manag Sci. 2019;75:1242–51.
Wang Z, Ji H, Yuan B, Wang S, Su C, Yao B, Zhao H, Li X. ABA signalling is fine-tuned by antagonistic HAB1 variants. Nat Commun. 2015;6:8138.
Liu Q, Wang C, Jiao X, Zhang H, Song L, Li Y, Gao C, Wang K. Hi-TOM: a platform for high-throughput tracking of mutations induced by CRISPR/Cas systems. Sci China Life Sci. 2019;62:1–7.
Jiang Y, Chai Y, Qiao D, Wang J, Xin C, Sun W, Cao Z, Zhang Y, Zhou Y, Wang XC, Chen QJ. Optimized prime editing efficiently generates glyphosate-resistant rice plants carrying homozygous TAP-IVS mutation in EPSPS. Mol Plant. 2022;15:1646–9.
Qiao D, Wang J, Lu MH, Xin C, Chai Y, Jiang Y, Sun W, Cao Z, Guo S, Wang XC, Chen QJ. Optimized prime editing efficiently generates heritable mutations in maize. J Integr Plant Biol. 2023;65:900–6.
Garcia MD, Nouwens A, Lonhienne TG, Guddat LW. Comprehensive understanding of acetohydroxyacid synthase inhibition by different herbicide families. Proc Natl Acad Sci U S A. 2017;114:E1091–100.
Lin HY, Dong J, Dong J, Yang WC, Yang GF. Insights into 4-hydroxyphenylpyruvate dioxygenase-inhibitor interactions from comparative structural biology. Trends Biochem Sci. 2023;48:568–84.
Nishimasu H, Shi X, Ishiguro S, Gao L, Hirano S, Okazaki S, Noda T, Abudayyeh OO, Gootenberg JS, Mori H, et al. Engineered CRISPR-Cas9 nuclease with expanded targeting space. Science. 2018;361:1259–62.
Schindler CEM, Hollenbach E, Mietzner T, Schleifer KJ, Zacharias M. Free energy calculations elucidate substrate binding, gating mechanism, and tolerance-promoting mutations in herbicide target 4-hydroxyphenylpyruvate dioxygenase. Protein Sci. 2019;28:1048–58.
Li C, Zhang R, Meng X, Chen S, Zong Y, Lu C, Qiu JL, Chen YH, Li J, Gao C. Targeted, random mutagenesis of plant genes with dual cytosine and adenine base editors. Nat Biotechnol. 2020;38:875–82.
Liu X, Qin R, Li J, Liao S, Shan T, Xu R, Wu D, Wei P. A CRISPR-Cas9-mediated domain-specific base-editing screen enables functional assessment of ACCase variants in rice. Plant Biotechnol J. 2020;18:1845–7.
Ren Q, Sretenovic S, Liu G, Zhong Z, Wang J, Huang L, Tang X, Guo Y, Liu L, Wu Y, et al. Improved plant cytosine base editors with high editing activity, purity, and specificity. Plant Biotechnol J. 2021;19:2052–68.
Wang H, He Y, Wang Y, Li Z, Hao J, Song Y, Wang M, Zhu JK. Base editing-mediated targeted evolution of ACCase for herbicide-resistant rice mutants. J Integr Plant Biol. 2022;64:2029–32.
Zhang A, Shan T, Sun Y, Chen Z, Hu J, Hu Z, Ming Z, Zhu Z, Li X, He J, et al. Directed evolution rice genes with randomly multiplexed sgRNAs assembly of base editors. Plant Biotechnol J. 2023;21:2597–610.
Ma Y, Zhang J, Yin W, Zhang Z, Song Y, Chang X. Targeted AID-mediated mutagenesis (TAM) enables efficient genomic diversification in mammalian cells. Nat Methods. 2016;13:1029–35.
Zhao D, Li J, Li S, Xin X, Hu M, Price MA, Rosser SJ, Bi C, Zhang X. Glycosylase base editors enable C-to-A and C-to-G base changes. Nat Biotechnol. 2021;39:35–40.
Tong H, Wang X, Liu Y, Liu N, Li Y, Luo J, Ma Q, Wu D, Li J, Xu C, Yang H. Programmable A-to-Y base editing by fusing an adenine base editor with an N-methylpurine DNA glycosylase. Nat Biotechnol. 2023;41:1080–4.
Ye L, Zhao D, Li J, Wang Y, Li B, Yang Y, Hou X, Wang H, Wei Z, Liu X, et al. Glycosylase-based base editors for efficient T-to-G and C-to-G editing in mammalian cells. Nat Biotechnol. 2024:1–10. https://doi.org/10.1038/s41587-023-02050-w.
Chen L, Hong M, Luan C, Gao H, Ru G, Guo X, Zhang D, Zhang S, Li C, Wu J, et al. Adenine transversion editors enable precise, efficient A*T-to-C*G base editing in mammalian cells and embryos. Nat Biotechnol. 2024;42:638–50.
Gaudelli NM, Komor AC, Rees HA, Packer MS, Badran AH, Bryson DI, Liu DR. Programmable base editing of A*T to G*C in genomic DNA without DNA cleavage. Nature. 2017;551:464–71.
Grunewald J, Zhou R, Lareau CA, Garcia SP, Iyer S, Miller BR, Langner LM, Hsu JY, Aryee MJ, Joung JK. A dual-deaminase CRISPR base editor enables concurrent adenine and cytosine editing. Nat Biotechnol. 2020;38:861–4.
Sakata RC, Ishiguro S, Mori H, Tanaka M, Tatsuno K, Ueda H, Yamamoto S, Seki M, Masuyama N, Nishida K, et al. Base editors for simultaneous introduction of C-to-T and A-to-G mutations. Nat Biotechnol. 2020;38:865–9.
Xie J, Huang X, Wang X, Gou S, Liang Y, Chen F, Li N, Ouyang Z, Zhang Q, Ge W, et al. ACBE, a new base editor for simultaneous C-to-T and A-to-G substitutions in mammalian systems. BMC Biol. 2020;18:131.
Zhang X, Zhu B, Chen L, Xie L, Yu W, Wang Y, Li L, Yin S, Yang L, Hu H, et al. Dual base editor catalyzes both cytosine and adenine base conversions in human cells. Nat Biotechnol. 2020;38:856–60.
Chauhan J, Prathibha MD, Singh P, Choyal P, Mishra UN, Saha D, Kumar R, Anuragi H, Pandey S, Bose B, et al: Plant photosynthesis under abiotic stresses: Damages, adaptive, and signaling mechanisms. Plant Stress. 2023;10:100296.
McAinsh MR, Pittman JK. Shaping the calcium signature. New Phytol. 2009;181:275–94.
Pirayesh N, Giridhar M, Ben Khedher A, Vothknecht UC, Chigri F. Organellar calcium signaling in plants: an update. Biochim Biophys Acta Mol Cell Res. 2021;1868: 118948.
Alonso MT, Rodriguez-Prados M, Navas-Navarro P, Rojo-Ruiz J, Garcia-Sancho J. Using aequorin probes to measure Ca(2+) in intracellular organelles. Cell Calcium. 2017;64:3–11.
Sun S, Zhang X, Chen K, Zhu X, Zhao Y. Screening for Arabidopsis mutants with altered Ca(2+) signal response using aequorin-based Ca(2+) reporter system. STAR Protoc. 2021;2: 100558.
Keller JM, Frieboes MJ, Jodecke L, Kappel S, Wulff N, Rindfleisch T, Sandoval-Ibanez O, Gerlach I, Thiele W, Bock R, et al. Eukaryote-specific assembly factor DEAP2 mediates an early step of photosystem II assembly in Arabidopsis. Plant Physiol. 2023;193:1970–86.
Zhang Y, Li J, Gao C. Generation of stable transgenic rice (Oryza sativa L.) by Agrobacterium-mediated transformation. Curr Protoc Plant Biol. 2016;1:235–46.
Li J, Chen L, Liang J, Xu R, Jiang Y, Li Y, Ding J, Li M, Qin R, Wei P. Development of a highly efficient prime editor 2 system in plants. Genome Biol. 2022;23:161.
Zhang X, Henriques R, Lin SS, Niu QW, Chua NH. Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method. Nat Protoc. 2006;1:641–6.
Acknowledgements
We thank Qijun Chen from China Agricultural University for sharing the pHEE401E plasmid.
Review history
The review history is available as Additional file 3.
Peer review information
Tim Sands was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
Funding
This work was supported by the Taishan Scholar Foundation of Shandong Province (TSQN202103160 to Huawei Zhang), the Excellent Youth Foundation of Shandong Scientific Committee (ZR202103010168 to Huawei Zhang), and the National Key R&D Program of China (2022YFD1201700 to Lijing Liu).
Author information
Authors and Affiliations
Contributions
H.Z. and L.L. designed the experiments. X.W., W.P., C.S., H.Y., Z.C. F.Y., G.M., and Y. S. performed the experiments. X.W., W.P., C.S., and R.Z. analyzed the results. H.Z., L.L., and C.G. supervised the project and wrote the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors’ patent based on the results reported in this paper has been granted (patent number ZL202211484789.3) and does not prevent other from using the described method to screen.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13059_2024_3358_MOESM1_ESM.docx
Additional file 1: Fig. S1. Schematic overview of the pipeline for generating large scale mutant pool. Fig. S2. Schematic overview of targets used for checking the editing efficiency of pHEE901-A3A. Fig. S3. The mutated alleles generated by pHEE901-A3A at AG-T2, ARF2 and HAB1 target sites in the T1 generation. Fig. S4. The mutated alleles generated by pHEE901-A3A at AG and EPSPS sites in the T1 generation. Fig. S5. The mutated alleles generated by pHEE901-A3A at AG and EPSPS sites in the T2 generation. Fig. S6. The amino acids sequences at AG and EPSPS targets in the T1 and T2 generation mutants. Fig. S7. The DNA sequences encoding the effective mutated proteins at AG target in the T2 generation. Fig. S8. The alignment of ALSs from Arabidopsis and rice. Fig. S9. The numbers of T1 plants for each sgRNA used for ALS evolution. Fig. S10. Mutation forms at the DNA and protein levels at the ALS targets detected in herbicide-resistant plants. Fig. S11. No obvious growth penalty was observed inALS mutants. Fig. S12. The herbicide resistant phenotypes of different ALS mutants. Fig. S13. No obvious growth penalty was observed in OsALS genome edited rice mutants. Fig. S14. The herbicide resistant phenotypes of OsALS genome edited rice mutants. Fig. S15. Sequence alignment of HPPDs from Arabidopsis and crops. Fig. S16. Mutation forms at the DNA and protein levels at the HPPD targets detected in herbicide-resistant plants. Fig. S17. Amino acid substitutions detected in herbicide resistant HPPD mutant forms. Fig. S18. The diagram and sequences used for the prime editing of ALS in rice.
13059_2024_3358_MOESM2_ESM.xlsx
Additional file 2: Table S1. The ratios of different mutation types at the five endogenous targets in the T1 generation. Table S2. The base editing efficiency at each nucleotide in the protospacers of five targets in T1 plants. Table S3. The ratios of different mutation types at AG and EPSPS sites in the T2 generation. Table S4. The genotypes and mutation forms of herbicide resistant plants screened out by the efficient germline-specific evolution system. Table S5. The primers used in this study.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Wang, X., Pan, W., Sun, C. et al. Creating large-scale genetic diversity in Arabidopsis via base editing-mediated deep artificial evolution. Genome Biol 25, 215 (2024). https://doi.org/10.1186/s13059-024-03358-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13059-024-03358-9