Introduction

Genetic erosion was proposed by Harlan1 in 1975 to describe the genetic resources after the Green revolution, which involved directional selection of semi-dwarf genes in crop plants. The semi-dwarf genes improved plant architecture in response to heavy use of nitrogen fertilizer2,3. With the modernization of agriculture and rapid urbanization, only a few modern varieties are now cultivated, in contrast to the large number of landraces used in rice production prior to the 1980s4. This has not only led to a genetic bottleneck, but has also intensified the permanent loss of many landraces5,6.

In China, the genetic erosion in rice is even more serious due to three main reasons: (1) from 1930–1960, several rice landraces (such as NanTe and Shenglixian) were dominant in the main rice-producing region of the Yangtze rice valley. (2) From 1950–1970, spontaneous Semi-Dwarf1 (SD1) mutant alleles (sd1) were derived from the main cultivars7, such as NanTe and Zhaiyeqing, or from local Taiwan varieties Dijiaowujian and Aizaizhan, and then were backcrossed into a few main cultivars. The resulting plants underwent heavy selection at the sd1 locus as well as other loci to give decreased plant height for better lodging resistance and adaptation to high Nitrogen fertilizer condition and gaining high harvest index. The main cultivars selected in the 1930s–1950s contained the sd1 mutant alleles to confer shorter plants and were used as the main breeding materials for new varieties8,9. (3) From 1970 to today, hybrid rice varieties are used to take advantage of the heterosis between two genetically distant rice varieties with semi-dwarf plant architecture and to enhance nitrogen utilization. However, increases in yield have been limited since the late 1990s. In addition, pedigree analysis showed that the parents of the hybrids had similar allelic variations and many of the same traceable linkage blocks, which when combined, may not improve yield- and resistance-related traits10,11. Together, these activities resulted in a very narrow genetic pool in the main modern Chinese rice varieties.

Relieving genetic erosion and improving the yield of modern rice varieties to satisfy food supply demands and foster sustainable development is an urgent issue. In our study, we used a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) genome editing technique to edit Semi-Dwarf1 (SD1) and Photosensitivity5 (SE5) in the elite landraces Kasalath and TTP (TeTePu), which contain many desirable agronomic traits such as tolerance to low phosphorous12 and broad-spectrum resistance to diseases and insects13. Our results showed that precise targeting of SD1 for gene editing in Kasalath or TTP resulted in new lines with a semi-dwarf plant architecture, which is desired in modern rice varieties, and maintained most of the desired agronomic traits of their progenitors. We show that using gene editing on available landraces can rapidly increase genetic diversity and produce new varieties that satisfy current production requirements.

Results

Kasalth and TTP sd1 mutants are resistant to lodging and show good field performance

To rapidly create new germplasm with decreased plant height in the traditional landraces Kasalath and TTP, we constructed the CRISPR/Cas9 vector for simultaneously targeting SD1 and Photoperiod-sensitivity-5 (SE5) based on a previous study14. The guide RNA (gRNA1 and gRNA2 targeting sites were designed in the first and second exon of SD1, and the gRNA3 targeting site was designed in the first exon of SE5. The three target sequences (Fig. 1a) were cloned into SK-gRNA, then two gRNAs were assembled into one intermediate vector (SK-gRNA1SD1 and SK-gRNA2SD1 were assembled into SK-gRNA1SD1-gRNA2SD1; SK-gRNA1SD1 and SK-gRNA3SE5 were assembled into SK-gRNA1SD1- gRNA3SE5). The single guide RNAs (sgRNAs) were cloned into the pC1300-Cas9 expression vector (Fig. 1b).

Figure 1
figure 1

Phenotypes of sd1 and se5/ sd1 mutants created by CRISPR-Cas9 in the Kasalath and TTP backgrounds. (a) SD1 and SE5 target loci. The gRNA1 and gRNA2 targeting sites were designed in the first and second exons of SD1, and the gRNA3 targeting site was in the first exon of SE5. The target sites are labeled in black lowercase letters. The protospacer adjacent motif (PAM) sequences are underlined and in red type. (b) The two vectors used to tranform rice. gRNA1 and gRNA2 were assembled into the pC1300-Cas9 expression vector for Kasalath tranformation; gRNA1 and gRNA3 were assembled into the pC1300-Cas9 expression vector for TTP tranformation. (c) Phenotypes of different sd1 mutants in the Kasalath background. (Bar: 15 cm). (d) The internode length proportion in kasalth, sd1-1Kas, sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas mutants. The bar in the image of whole plants represents 4 cm. (e) Comparison of plant height among wild type (Kasalath), sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas. Values in plant height are means±standard deviation (±sd.), n = 10. Bars with different letters are significantly different. Statistical differences among the agronomic traits were detected by Duncan’s multiple range test (p < 0.05). (f) Comparison of plot yield (kg of grain per 30 plants) among wild type (Kasalath), sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas. Bars with different letters are significantly different. Data are means ± standard deviation, n = 3. Statistical differences among the agronomic traits were detected by Duncan’s multiple range test (p < 0.05). (g) Phenotypes of sd1 and sd1/ se5 mutants in the TTP background. Two sd1 mutants and one se5/ sd1 double mutant in TTP. The bar represents 20 cm. (h) Diseased leaves of 50- to 60-day-old plants from the field evaluation for resistance to blast: NPB (upper), TTP (upper the second), sd1 or sd1/se5 mutants (below). The lesion area of the TTP and TTP mutants were near 0 for all races tested.

Full size image

The binary vectors targeting SD1 or SD1 and SE5 were used for genetic transformation. A total of 25 and 29 positive transgenic plants were obtained in the Kasalath and TTP T0 generation, respectively (Supplementary Table 1). We then sequenced the target regions in all positive transgenic plants to analyze the mutations of the target sites. In the gRNA1 and gRNA2 targeting site in Kasalath, the mutation rate was 44 and 68%, respectively. In the gRNA1 and gRNA3 targeting site in TTP, the mutation rate was 28 and 21%, respectively (Supplementary Table 1). This may reveal different editing efficiencies in different genetic backgrounds or accessions.

To eliminate the selective maker gene and T–DNA element, and to evaluate the phenotype of the new sd1 mutants, T0 progeny plants were analyzed by DNA sequencing and self-pollinated to generate T1 progeny plants. T2 progeny plants that were either homozygous sd1/sd1 or T-DNA free plants were analyzed by phenotypic screening and genotyping (Fig. S1). Sequencing analysis of Kasalath lines revealed that sd1-1Kas, sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas had a 1-nucleotide (nt), 2-nt, 1-nt, 3-nt, and 61-nt deletions at the predicted editing site, respectively. The sd1-1Kas line had a 1-nt deletion next to the gRNA1 protospacer-adjacent motif (PAM), and the sd1-4Kas line had a 1-nt deletion next to the gRNA1 PAM and a 2-nt deletion next to the gRNA2 PAM, creating the same stop codon, leading to premature termination. The sd1-2Kas line had a 2-nt deletion next to the gRNA2 PAM motif leading to premature termination. The sd1-3Kas line had a 1-nt deletion next to the gRNA2 PAM motif leading to premature termination. The sd1-5Kas line had a 61-nt deletion next to the gRNA2 PAM motif leading to premature termination (Fig. S2a). All of these mutations would lead to a functionally altered proteins, and have different mutation sites compared with previously reported sd1 alleles. The previously reported sd1 spontaneous mutation sites were mostly detected in the gibberellin (GA) 20-oxidase (76-371) motif, varying from a deletion to a single base substitution, indicating that this motif is important for SD1 function. In our study, the new CRISPR/Cas9 mutation sites were also affected in this region, and the lines showed a dwarf phenotype (Fig. 1c), which supported this hypothesis.

We next evaluated the phenotype of the sd1 mutants. Although the plant height of the sd1-1Kas, sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas lines was shorter than the WT (Fig. 1c,e), the reduction in internode length was proportional in the mutants compared to the corresponding internode length of the WT (Fig. 1d). As the sd1-1Kas and sd1-4Kas lines contain the same base deletion at the first editing target, leading to premature transcription termination and a similar mutant phenotype, we chose the sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas lines for further research.

To test whether these four allelic variations impact yield, we measured yield-related agronomic traits including plant height, grain number per panicle, etc. The average plant height of the sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas lines were 82.5, 87.62, 87.16, and 91.5 cm, respectively (Fig. 1e). Additionally, the grain number per panicle in all the new sd1 mutants was decreased by various degrees when compared with the WT, and the sd1-2Kas line had a significant reduction in the number of spikelets. Conversely, tiller number was notably increased in the mutant lines compared with the WT (Table 1). The heading date was delayed from 3–5 days in the sd1 mutants (Table 1). For example, the heading date of the sd1-2Kas line was delayed by 4.3 days, and the heading date of the sd1-5Kas line was delayed by 3 days (Table 1) which showed different effects in the different sd1 alleles. We calculated the average yield from 30 plants from each of the sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas lines. We found that the sd1-2Kas line had a strong reduction in yield when compared to the WT and may not be suitable for breeding. The sd1-3Kas and sd1-5Kas lines showed similar yields to the WT and may be potential breeding materials.

Table 1 Agronomic traits comparison between Kasalath and its sd1 mutant lines.
Full size table

To study the lines in response to nitrogen fertility, we treated Kasalath and the sd1-3Kas and sd1-5Kas lines with different concentrations of nitrogen fertilizer in the field. Under high nitrogen, Kasalath is more sensitive than Kasalath sd1 to nitrogen (Fig. S3a, Table 2). Additionally, Kasalath was more susceptible to lodging than the sd1-5Kas line (Figs. S3b, S4). These results suggest that the sd1 mutation confers better lodging resistance and better nitrogen utilization.

Table 2 Yield-related traits comparion between Kasalath and its sd1 mutant lines under different Nitrogen treatment condition.
Full size table

Simultaneously, we obtained new sd1 single mutants and sd1 se5 double mutants by gene editing of SD1 and SE5 in the TTP background carrying no selection marker (Fig. S1b). PCR check SD1 gene of these stable inherited mutant lines confirmed there are different mutation of SD1 or SE5 in these corresponding mutant lines (Figs. 1g, S2b). As TTP has broad-spectrum resistance to the rice blast fungus, we inoculated the leaf sheaths of TTP and its mutants with mix of Magnaporthe oryzae strain CH102 and CH184. and found that the sd1 mutants had slightly increased resistance to rice blast compared with TTP (Fig. 1h). In addition, TTP is a landrace of Vietnam which is a thermo-sensitive and late-flowering in main China rice production area. Knock-out Se5 gene (LOC_Os06g40080) in TTP lead to 2 weeks early flowering than TTP in Hangzhou, China. Overall, using the CRISPR-Cas9 gene editing technique, we obtained new sd1 mutant alleles in the elite landraces Kasalath and TTP.

The desirable agronomic traits of the progenitors were maintained in the new sd1 mutant rice lines

To investigate whether the sd1 mutation affects the GA response, we tested the GA response of the sd1-3Kas, sd1-4Kas, and sd1-5Kas lines. Treatment with 1, 10, and 100 mM GA3 had similar effects in Kasalath and its mutants (Fig. S5a,c). The effect of different concentrations of GA3 on the length of the second internode was also similar between the WT and the mutants (Fig. S5b), indicating that the sd1-3Kas, sd1-4Kas, and sd1-5Kas lines have a normal sensitivity to GA. Consistently, in GA assay of TTP, sd1-1-TTP,sd1-2-TTP, we got similar physiological results with sd1-3Kas, sd1-4Kas, and sd1-5Kas lines (Fig. S6). All these results showed that these sd1mutants attained by Crispr-Cas9 are also GA deficiency mutants as reported those spontaneous sd1 mutants.

Seed dormancy has been reported to be associated with red grain color. Previous studies have shown that Kasalath, which has red grain, contains a cluster of quantitative trait loci at Os07g11020/Rc for seed dormancy compared to the white rice cultivar Nipponbare which is 14 bp deletion at Rc coding region15. An appropriate level of seed dormancy can decrease pre-harvest sprouting, and thus improve grain yield and quality. Therefore, we analyzed whether the sd1-2Kas, sd1-4Kas, and sd1-5Kas lines have increased seed dormancy like Kasalath (Fig. 2a–c). Seed dormancy was evaluated by measuring the germination of seeds after imbibing for 6 days.

Figure 2
figure 2

Comparison of germination among Kasalath, sd1-2Kas, sd1-4Kas, sd1-5Kas. (a) Images of caryopsis morphologies of Kasalath and the sd1 mutants. (b) Germination rates of Kasalath, sd1-2Kas, sd1-4Kas, and sd1-5Kas seeds at different concentrations of GA. Germination was evaluated at 6 d after imbibition and with three samples (100 seeds/sample). (c) Germination phenotypes of Kasalath, sd1-2Kas, sd1-4Kas, sd1-5Kas. The germination phenotypes are shown at 6 d after imbibition. Values means ± standard deviation (±sd.), n > = 3. Statistical differences were detected by Duncan’s multiple range test (p < 0.05).

Full size image

We first sequenced Os07g11020/Rc in Nipponbare, Kasalath, and the sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas mutants and found a 14-bp deletion in the rc from rice variety 93-11 and Nipponbare. Meanwhile, we tested two dormancy-breaking treatments to evaluate the difference in the germination response among Kasalath and its mutants and found that the mutation of sd1 can, to some degree, increase seed dormancy (Fig. 2b,c), which is similar to findings in a previous study16. Together, our results suggest that the sd1 mutation in Kasalath may decrease the endogenous GA concentration and enhance seed dormancy, which may be beneficial to rice production during harvest under high temperatures and wet conditions.

Phosphorus deficiency has a detrimental impact on plant growth. Phosphatic fertilizers could relieve phosphorus deficiency, but low use efficiency of available phosphorus in rice varieties is a bottleneck and could lead to environmental consequences. Therefore, it is important to use landraces containing genes for high phosphorus use efficiency when breeding new rice varieties. A major quantitative trait locus for tolerance to phosphorus deficiency, Pup1 (also named Pstol1), was identified in Kasalath12. This gene is absent from Nipponbare. In order to explore whether the sd1 mutation in Kasalath impacts phosphorus use efficiency, we conducted phenotypic analyses of Kasalath, sd1-3Kas, sd1-5Kas, in nutrient solution with low P (0.5 mg/L), CK (10 mg/L) and high P (25 mg/L) hydroponics solution for 18 days (Fig. 3a,b). Under low P, the root lengths and surface area of sd1-3Kas, sd1-5Kas were increased compare with Kasalath but decline in the CK and high-P concentration (Fig. 3a,b). These results demonstrated that the sd1 mutation does not disrupt the function of PSTOL1 in the Kasalath mutants under low P.

Figure 3
figure 3

30-day-old of root growth response under different concentration of P. (a) Total root length and surface area of Kasalath, sd1-3Kas, and sd1-5Kas in low -P (0.5 mg/L), CK (10 mg/L) and high-P (25 mg/L) hydroponics solution for 20 days. Error bars indicate standard error. (b) The root architecture of Kasalath, sd1-3Kas, and sd1-5Kas plants under different concentration of P. Values means±standard deviation (±sd), n > = 3. Statistical differences were detected by Duncan’s multiple range test (p < 0.05).

Full size image

Flowering date (also known as heading date) is an important agronomic trait in rice. Photoperiod regulation is an important factor for controlling heading date. Ghd8(Grain Yield, Plant Height and Heading Date8) is a major quantitative trait locus associated with pleiotropic effects on grain yield, heading date, and plant height17. Mutation of Ghd8 leads to earlier flowering. Heading date 1 (Hd1)18 represses flowering under long-day (LD) conditions and induces it under short-day (SD) conditions. Mutants of Hd1 in Kasalath exhibited no photoperiod response under LD conditions. Knockout mutants of OsGSK1 which is an orthologue of Arabidopsis BIN2 showed enhanced tolerance to cold, heat, salt, and drought stresses19. OsGSK1 in Kasalath contains many single nucleotide polymorphisms (SNPs) and insertions/deletions compared with Nipponbare and 93-11. S5n can restore the sterility in indica and japonica hybrids20, which is an important resource when utilizing heterosis between these rice subspecies. S5n- wide compatible gene exits in Kasalath. In order to test whether the CRISPR-Cas9 editing technique affected Ghd8, hd1, OsGSK1, and S5n and the phenotypic consequences, we sequenced Rc-HLH, Ghd8, hd1, and OsGSK in Kasalath and the sd1-2Kas, sd1-3Kas, sd1-4Kas, and sd1-5Kas mutants. In line with previous results, Rc-HLH contains a 14-bp insertion in Kasalath and its mutants. Ghd8 contained the same SNPs and deletion in Kasalath and its mutants. hd1 contained the same insertion mutation in Kasalath and its mutants. OsGSK1 contained the same insertion in Kasalath and its mutants. In addition, the loss of function of S5n was observed in Kasalath and its mutants (Fig. S7).

In the TTP mutant rice lines, we sequenced the Pi54 gene21 which confer high degree of resistance to diverse isolates of M. oryzae and found no changes in these rice lines, which is consistent with the resistance test in the field (Fig. S8).All these results show that the sd1 mutation in Kasalath does not have negative effects on other genes associated with many desirable agronomic traits.

In summary, The CRISPR-Cas9 editing technique can be used to more rapidly create the sd1 mutation in desirable germplasm. Two sd1 alleles in Kasalath and one double sd1 se5 mutant created in this study can serve as potential materials for breeding. Creating additional sd1 alleles in other desirable landraces could help to improve genetic diversity in rice and benefit rice production.

Materials and Methods

Plant materials and Measurements of agronomic traits

The background of transgenic plants is Kasalath or TTP and all the rice plants were grown in the paddy fields under natural conditions in Hangzhou or LingShui (China National Rice Research Institute, China). The agronomic traits were analyze after rice harvested, a total of 10 or 30 randomly chosen rice plant were used to measure plant height, grain number per panicle and tiller number, heading date.

Plasmid construction and Plant transformation

The three target sites were designed for knock out of SD1 or SE5 genes using the CRISPR/Cas9 system. The gRNA1SD1 (digested with Kpn I/BamH I) and gRNA2SD1 (digested with Kpn I/Bgl II) were assembled into one intermediate vector. Similar methods were used to assemble gRNA1SD1and gRNA3SE5 into one intermediate vector. The two intermediate vector (digested with Kpn I/Bgl II) was assembled to the pC1300-Cas9 binary vector (digested with Kpn I/BamH I), respectively. The target sequences are provided in Supplementary Table S2 in Supporting Information. The pC1300-Cas9 binary vector loading two sgRNAs was used for genetic transformation via the Agrobacterium-mediated transformation (strain EHA105) method for generating transgenic rice, according to Japonica rice and Indica rice transformation methods22,23. Detection of mutations Genomic DNA of transgenic plants was extracted from approximately 100 mg leaf tissue of rice via the cetyltrimethylammoniumbromide (CTAB) method. PCR was conducted with KOD FX DNA polymerase (Toyobo, Japan) to amplify the fragments surrounding the three target sites. The DNA fragments were sequenced by the Sanger method and analyzed by the degenerate sequence decoding method24.

Seed germination rate measurement

Seed dormancy was evaluated by germination of seed samples after-ripened for 5 days. Harvested rice was dried at 42 °C for 2 days and stored at room temperature before testing. Seeds of Kasalath and T2 of the sd1-2Kas, sd1-4Kas, sd1-5Kas were used in this study. Seed dormancy was evaluated by germination testing. A sample of approximately 100 well-developed seeds was distributed on filter paper in a in 23 cm × 23 cm aseptic culture dishes with 150 ml distilled water at 30 °C and 100% relative humidity. An artificial climate chamber (POX-330B-22H, Life Apparatus Co., Ningbo, China) was used for the seed dormancy treatments. Germination was evaluated visually by protrusion of the radicle from the hull by more than 3 mm and counted daily from day 2 to day 6 or at day 6. A test was replicated three times, and germination percentages were averaged for genetic analysis.

GA assay on 93-11, Nipponbare, Kasalath and new sd1 mutant lines

Seeds of Kasalath, T2 of the sd1-4Kas, sd1-5Kas or TTP, sd1-TTP were pre-germinated in the dark at room temperature. After 3 days, germinated seeds pre-germinated in the dark at room temperature were transferred to Yoshida culture solution with 0 mM, 1 mM,10 mM,100 mM GA3 respectively. The solution was replaced every 3 days. The plant height and second internode length of seedlings (10 DAG) were analyzed.

Root scan of Kasalath, sd1-3Kas, sd1-5Kas grown in hydroponics

Seeds of the Kasalath and T2 of the sd1-3Kas, sd1-5Kas were pre-germinated in Petri dishes in the dark at room temperature. After 3 days, germinated seeds were transferred to Yoshida culture solution grow for 7 days. Then plants were transferred to Yoshida culture solution with low -P (0.5 mg/L), CK (10 mg/L) and high-P (25 mg/L) respectively12. The solution was replaced every 3 days. The root length of seedlings (30 DAG) were checked.

N fertilizer treatments

Field experiments were conducted in the same field, HangZhou ZheJiang Province, China. The four N treatments were 0 (No N), 8 (Low N), 14 (moderate N) and 20 (high N) kg N ha−1. In each treatment, N was applied at the basal, tillering and panicle initiation stages at a ratio of 4: 3: 3. Phosphorus was applied as a basal fertilizer at a rate of 5.6 kg P ha−1 and K was applied equally between the basal and panicle initiation stages at 9.6 kg ha−1. Kasalath, sd1-3Kas, sd1-5Kas were used in the N treatments25.

Assessment of disease phenotypes with Magnaportheoryzae

Fifty-day-old of NPB and TTP, sd1-1TTP, sd1-2TTP, sd1-3/se5TTP plants were inoculated with Magnaporthe oryzae spore suspension (1 × 105 spores/ml). After 7 days post inoculation, disease reaction of each rice line was photoed26.