1. Introduction
Rice growth and yield are influenced by light-temperature interaction, affecting flowering, fruiting, and yield formation. The flowering period, a complex trait governed by multiple genes [
1], plays a crucial role in rice adaptation and yield optimization [
2] by maximizing light and temperature resources [
3]. Modifying the flowering period expands the cultivation area for high-quality rice varieties, addressing production challenges. Rice blast, caused by
Magnaporthe oryzae, poses a significant threat to rice yield stability [
4]. Developing disease-resistant rice varieties is pivotal for effective pathogen control.
Traditional rice breeding is time-consuming due to limited heritable variation and extensive phenotypic screening. Genetic engineering techniques, particularly CRISPR/Cas9, which is a cutting-edge gene editing tool and widely applied in crop genetic research for revolutionize trait improvement by precise DNA-level editing, accelerating breeding with enhanced efficiency and targeted outcomes. The approach offers rapid enhancements in flowering date, grain quality, stress resistance, and crop yield, minimizing time costs [
5,
6,
7,
8,
9]. Noteworthy studies include editing the
DTH8 gene to alter tasselling periods of
Geng rice 99-25 and obtaining
DTH8 mutant material with an earlier tasselling period [
10]. In addition, this technology has been used to enhance resistant starch content through SS3a and SS3b editing, achieving a significant increase of 4.7%-5.0% in single knockouts, and the double-knocked out mutant exhibiting significant differences, with the resistant starch increasing 9.5%-9.7%, creating new nutritious rice germplasm [
11]. Wang Kejian's team at the China National Rice Research Institute employed CRISPR/Cas9 gene editing technology to knock out four reproduction-related genes,
PAIR1,
REC8,
OSD1, and
MTL, in the hybrid rice Chunyou 84, achieving fusion-free reproduction of hybrid rice from 0 to 1, which is of great theoretical and practical significance [
12]. Zhou Wenjia used CRISPR/Cas9 to edit the flowering period gene
Hd2 and the aroma gene
Badh2 of Sui-Geng 14, to shorten the flowering period [
13]. Zhou employed CRISPR/Cas9 to edit broad-spectrum rice blast resistance genes
Bsr-d1,
Pi21, and
ERF922. Single- or triple-mutants displayed high resistance to rice blast, with
Pi21 or
erf922 single mutants exhibiting enhanced resistance to leaf blight [
14] Plant gene editing technology offers avenues for creating high-quality crop varieties, demonstrating the feasibility of directly editing genes in superior crop varieties for improved commercial crops. Cite a case, rice flowering is controlled by multiple genes, with allele variants influencing regional adaptability under varying day lengths.
Hd2, a core gene in rice photoperiodic flowering regulation system, functions as a flowering repressor under long daylight conditions [
15]. Its expression level correlates significantly with flowering time, regulating rice flowering by enhancing photoreceptivity [
16]
Ghd7, a major QTL, controls grain number per spike, plant height, and flowering period in rice [
17,
18]. In the photoperiodic pathway under prolonged sunlight,
Ghd7 represses
Ehd1 expression, inhibiting flowering through the Ghd7-Ehd1-Hd3a/RFT1 pathway, and upregulates
OsCOL10 and
OsMFT1 as downstream flowering repressors [
19,
20]. Similarly,
DTH8, a member of the rice HAP family (containing the HAP3D subunit), affects flowering, plant height, and grain number per spike [
21] under extended sunlight, repressing flowering through protein interactions with
Hd1,
Ghd7,
DTH7, and
HAPL1 [
22].
Previous studies have demonstrated the
Ehd1-dependent photoperiodic flowering pathway, three flowering repressor genes,
Hd2,
Ghd7, and
DTH8 negatively regulated during the flowering period in rice varieties. Delaying rice flowering under prolonged sunlight conditions produced increased plant height and number of grains per spike, leading to higher yields. Therefore, it is important to shorten the rice flowering period by editing
Hd2,
Ghd7, and
DTH8, to genetically improve the flowering period of
Geng rice variety.
Pi21, a recessive rice blast resistance gene [
23], contributes to persistent rice blast resistance in Japanese rice varieties. The recessive version of the
Pi21 gene has a deletion of the first and second proline domains corresponding to 21 bp and 48 bp, respectively, associated with its resistance to rice blast [
24].
Pi21 negatively regulates disease resistance, and its loss-of-function shows resistance to rice blast, and the use of this resistance is expected to overcome the problem of rice varieties prone to lose rice blast resistance.
Jiahe212 (JH212) is a traditional Geng rice variety renowned for its superior quality and serves as the maintainer line for Jiahe212A (JH212A), a high-quality hybrid line widely used in southern China. However, hybrid combinations involving JH212A often exhibit prolonged growth periods and low heat tolerance, hindering their widespread adoption. To address this, shortening the flowering period can mitigate heat stress during sowing, allowing for more flexible planting arrangements and facilitating the expansion of cultivation areas, especially northward. In this study, we employed the CRISPR/Cas9 system to conduct multi-gene editing targeting key genes in rice flowering regulation (Hd2, Ghd7, and DTH8), the rice blast resistance gene Pi21, and the flavor gene Badh2 in JH212. The objective was to develop homozygous and stable mutant lines that enhance disease resistance and expedite the flowering period while preserving rice quality. These modifications aim to enhance breeding flexibility and promote wider adoption of high-quality Xian/Geng interspecific hybrid rice varieties.
3. Discussion
Rice yield and adaptability are influenced by factors like spiking duration, disease resistance of specific physiological minor varieties, and cropping schedules. While extending the spiking period is commonly seen as a method to boost yields, there is growing interest in short-duration high-quality rice varieties due to their flexible sowing nature, which allows for expansion without significant yield loss. The spiking period, a critical physiological phase determining rice variety adaptability [
25], has garnered attention from breeders. Improving the flowering period of existing
Geng varieties can expand high-quality rice cultivation and address flowering period-related challenges in popularized varieties. Maximizing rice yield relies on cultivating varieties with optimal flowering periods for the planting region. Gene editing offers a rapid and stable breeding pathway favored by experts [
26], boasting high efficiency, stability, and speed compared to traditional breeding methods like hybridization and backcrossing. CRISPR/Cas9 gene editing technology has rapidly advanced in recent years, extensively applied in creating germplasm resources [
27], genetic enhancement [
28], fertility improvement [
29], disease resistance [
14], and quality improvement of crops [
29]. Its adoption in rice research has become instrumental in unraveling gene functions and molecular mechanisms [
26].
In this study, we employed CRISPR/Cas9 gene editing technology to modify
Hd2,
Ghd7, and
DTH8 in the rice flowering pathway and rice blast resistance gene
Pi21. We successfully created new allelic variants in the JH212 background, establishing novel breeding materials for the flowering variation. In the T
0 generation, 88.23% efficiency was achieved, with 30 out of 34 transgenic seedlings showing positive results. Sequenced analyses revealed high mutation frequencies:
Hd2 (86.6%),
Ghd7 (83.3%),
DTH8 (83.3%), and
Pi21 (86.6%), predominantly involving base insertions and deletions. Our results align with previous studies: Li et al. achieved 77.8% simultaneous editing of
Hd2,
Ghd7, and
DTH8 across seven cultivars [
30]; Yang et al. obtained an 86.7% mutation rate for
Pi21; Xu et al. reported mutation frequencies of 75% for
Pita, 85% for
Pi21, and 65% for
ERF922 [
31]. These data fully proved the effectiveness of our new-constructed vector.
Important agronomic traits in rice are also jointly regulated by multiple genes or gene families. Using the CRISPR-Cas9 multiple knockout system, multiple genes can be targeted simultaneously to obtain plants with variable traits and to investigate the interactions among genes [
32], offering new materials for crop breeding. In this study, we found that the probability of single gene mutation was very limited, and many of the mutations occurred at the same time, primarily as polygenic mutations. Notably, we identified a total of eight homozygous mutation types for
Ghd7, six homozygous mutation types for
Hd2, nine homozygous mutation types for
Pi21, and four homozygous mutation types for
DTH8. These results are basically consistent with previous researches. Yang et al. utilized CRISPR/Cas9 to simultaneously knockdown eight genes highly expressed in the glutenin gene family and obtained seven mutants with varying combinations of mutations without affecting the content of storage substances and rice appearance, while the rice protein content was down-regulated to varying degrees [
33]. Shen et al successfully developed co-knockout vectors for eight genes associated with agronomic traits in rice using CRISPR-Cas9, and obtained homozygous six-mutant, seven-mutant, and eight-mutant strains [
34]. This approach provides a strategy for the rapid introduction of genetic diversity in the development of crop breeding. In constructing knockout lines, we have obtained materials to further study the dosage effect of alleles on phenotypes and genetic interactions between non-alleles, but also expanded the possibilities for the use of CRISPR/Cas9 gene editing technology in targeted editing of the genomes of rice varieties, to obtain polymerization of favorable alleles. In future work, we can better utilize these haplotypes to enhance the diversity of flowering times. Moreover, by integrating several genes with flowering-promoting impacts, it is possible to breed varieties with very early flowering times.
In the 17 homozygous mutant lines obtained, the four-gene homozygous mutant line JH-C18, the three-gene homozygous mutant line JH-C15, and the two-gene homozygous mutant line JH-C31 were chosen for the following experiments by observing phenotypes, and investigation of agronomic traits. In analyzing the expression of flowering-related genes in the gene-edited homozygous lines, the expression of
Hd3a in the homozygous mutant lines JH-C15, JH-C18, and JH-C31 was significantly higher than that of JH212, and the transcript levels of
Hd2,
Ghd7,
DTH8, and
Pi21, might be impacted after being target-edited by the Cas9 system, partially degrading the RNAs in vivo, producing a significant down-regulation of their gene expression in JH-C15, JH-C18, and JH-C31. Previous studies reported that in the
Ehd1-dependent photoperiodic flowering pathway,
Hd2,
Ghd7, and
DTH8 are the primary genes negatively regulated by flowering in rice varieties, delaying rice flowering under long sunlight conditions. Li et al. designed specific targets for three flowering genes, including
Hd2,
Ghd7, and
DTH8, showing that the knockout progeny could reduce the flowering period to varying degrees [
30]. Zhou et al. employed gene editing technology to target the flavor gene
Badh2 and the flowering period gene
Hd2 to obtain improved rice varieties with improved flavor and early maturity [
35]. Multi-target editing of
Pi21 and the flavor suppressor gene
OsBadh2 resulted in highly significant reductions in expression relative to the wild-type, and improved rice blast resistance and accumulation of 2-AP. In this study, the knockdown of
Hd2,
Ghd7, and
DTH8 resulted in significant advancement in the flowering period of mutant lines compared to the wild-type, indicating that the CRISPR/Cas9-mediated gene editing at the flowering period can achieve targeted fertility improvement, facilitate the introduction of high-quality rice resources, and accelerate local breeding programs. The analysis of genetic and molecular mechanisms of flowering regulation in rice and the development of flowering regulatory pathways will offer breeders more targeted and effective options for selecting the most suitable ecotypes for local growing areas and provide optimal forms and combinations of flowering gene alleles for future molecular design breeding. The use of CRISPR/Cas9 can regulate gene expression, creating new alleles that do not exist in nature and offering broad application prospects in future breeding of new breakthrough and diversified rice varieties.
Recently, the CRISPR/Cas9 gene editing system has become increasingly advanced. Compared to previous generations of gene editing technology, the CRISPR/Cas9 system offers advantages such as easier construction, lower cost, and simpler operation. However, CRISPR/Cas9 is not without flaws, and its high off-target rate has been a subject of criticism. In this study, the
Badh2 gene target experienced off-target effects. To address this issue, Fu et al. developed a series of progressively shorter guide RNAs (gRNAs) for the
EGFP reporter gene, including 15, 17, 19, or 20 complementary nucleotides. gRNAs with 17 or 18 complementary nucleotides could function at the intended target sites, reducing the off-target effects of paired Cas9 endonucleases without compromising the efficiency of targeted genome editing [37] Research has shown that gRNAs can direct Cas9 for on-target editing via base pairing with "seed sequences" located 8-12 nucleotides upstream of the protospacer adjacent motif (PAM). Therefore, appropriately lowering the length of the gRNA by 1-3 nucleotides can effectively reduce the off-target probability without impacting its targeting ability [38]. In addition, we selected physiological minor seeds with significant resistance against rice blast in varieties such as ZH11 and ZH8015 to inoculate JH212 with rice blast in this study [
31,
36]. The final plants did not get the relevant phenotypes, leading to the failure of the experiment. The reason for this failure may be due to the fact that JH212 is a variety with good blast resistance, or the fact that the selected physiological seed is not suitable for the inoculation, or the genetic effect of mutation type were not strong enough. In conclusion, despite certain limitations, CRISPR/Cas9 has attracted increasing attention from researchers due to its unique ability to precisely modify genes. This has prompted in-depth studies to overcome its challenges and examine its full potential. The three homozygous mutant lines obtained in this study, namely, JH-C15, JH-C18 and JH-C31, especially the last two, could serve as potential resources for expanding the promotion and application areas of JH212.
Figure 1.
Target site location data. a-e correspond to Hd2, Ghd7, DTH8, Pi21, and Badh2, respectively. The letters in the red box indicate the PAM sequence.
Figure 1.
Target site location data. a-e correspond to Hd2, Ghd7, DTH8, Pi21, and Badh2, respectively. The letters in the red box indicate the PAM sequence.
Figure 2.
Expression vector construction procedure.
Figure 2.
Expression vector construction procedure.
Figure 3.
PCR screening T0 generation positive strains. M, 2K DNA marker; 1-34 is the test line number.
Figure 3.
PCR screening T0 generation positive strains. M, 2K DNA marker; 1-34 is the test line number.
Figure 4.
Mutant genotype of the positive T1 plants tested.
Figure 4.
Mutant genotype of the positive T1 plants tested.
Figure 5.
PCR screening of T2 generation mutant lines without transgenic components and the phenotype of the three selected homozygotes lines. a and b, PCR screening of the Cas9 vector region and Hygromycin region in the single plants sampled from the three selected homozygotes lines. M, 2K DNA marker; 1-3, are negative control, distilled water, positive control, 4-16 are strains of JH-C15, 17- 30 are strains of JH-C18, a31-a45 and 31-45 are strains of JH-C31. c, e and g, the phenotype of JH-C15, JH-C18 and JH-C31 at flowering stage; d, f and h, the phenotype of JH-C15, JH-C18 and JH-C31 at harvest stage.
Figure 5.
PCR screening of T2 generation mutant lines without transgenic components and the phenotype of the three selected homozygotes lines. a and b, PCR screening of the Cas9 vector region and Hygromycin region in the single plants sampled from the three selected homozygotes lines. M, 2K DNA marker; 1-3, are negative control, distilled water, positive control, 4-16 are strains of JH-C15, 17- 30 are strains of JH-C18, a31-a45 and 31-45 are strains of JH-C31. c, e and g, the phenotype of JH-C15, JH-C18 and JH-C31 at flowering stage; d, f and h, the phenotype of JH-C15, JH-C18 and JH-C31 at harvest stage.
Figure 6.
qPCR analysis of flowering related genes in the selected homozygotes lines.
Figure 6.
qPCR analysis of flowering related genes in the selected homozygotes lines.
Figure 8.
Analysis of grain shape and agronomic traits of the selected homozygotes lines.
Figure 8.
Analysis of grain shape and agronomic traits of the selected homozygotes lines.
Table 1.
Off-target frequency assessment.
Table 1.
Off-target frequency assessment.
Target |
Position |
Strand |
GC(%) |
Region |
Potential off-target sites (Max score) |
Pairing with sg RNA (>= 8 nt) |
Hd2 |
738 - 757 |
+ |
55 |
CDS1 |
0.074 |
None |
Ghd7 |
366 - 385 |
+ |
60 |
CDS2 |
0.25 |
None |
DTH8 |
986 - 1005 |
+ |
55 |
CDS1 |
0.238 |
None |
Pi21 |
503 - 522 |
+ |
55 |
CDS1 |
0.133 |
None |
Badh2 |
499 - 518 |
+ |
60 |
CDS2 |
0.083 |
None |
Table 2.
Mutation frequency of positive plants in T0 generation.
Table 2.
Mutation frequency of positive plants in T0 generation.
Gene |
No. of plants |
Number of mutant strains |
Mutation rate |
Hd2 |
30 |
26/30 |
86.66% |
Ghd7 |
30 |
25/30 |
83.33% |
DTH8 |
30 |
25/30 |
83.33% |
Pi21 |
30 |
26/30 |
86.66% |
Table 3.
Genotype combination analysis of T1 homozygous lines.
Table 3.
Genotype combination analysis of T1 homozygous lines.
Line No.\Gene |
Hd2 |
Ghd7 |
DTH8 |
Pi21 |
Badh2 |
Flowering after sowing (Days) |
JH212 |
WT |
WT |
WT |
WT |
WT |
88.8±3.5 |
JH-C15 |
Allele1: +1bp |
WT |
Allele1: -4bp |
Allele3: +1bp |
WT |
57.38±2.75** |
JH-C17 |
Allele1: +1bp |
Allele1: +1bp |
WT |
Allele6: -2bp |
WT |
53.22±2.75** |
JH-C18 |
Allele6: -46bp |
Allele1: +1bp |
Allele1: -1bp |
Allele1: +1bp |
WT |
61.4±2.25** |
JH-C19 |
Allele2: +1bp |
Allele7: -28bp |
WT |
Allele2: +1bp |
WT |
56.55±3.25** |
JH-C22 |
Allele2: +1bp |
Allele1: +1bp |
WT |
Allele9: -19bp |
WT |
56.55±4.25** |
JH-C23 |
Allele3: +1bp |
Allele5: -5bp |
WT |
Allele2: +1bp |
WT |
50.5±2.75** |
JH-C24 |
Allele3: +1bp |
Allele4: -4bp |
WT |
Allele9: -19bp |
WT |
55.5±3.25** |
JH-C26 |
Allele1: +1bp |
Allele5: +1bp |
WT |
Allele6: -2bp |
WT |
50.5±2.25** |
JH-C27-5 |
Allele2: +1bp |
WT |
WT |
Allele7: -4bp |
WT |
51.55±4.25** |
JH-C27-2 |
WT |
Allele7: -28bp |
WT |
Allele7: -4bp |
WT |
52.55±3.25** |
JH-C28-5 |
Allele4: +1bp |
Allele2: +1bp |
WT |
Allele8: -7bp |
WT |
60.55±2.75** |
JH-C31 |
Allele5: -5bp |
WT |
WT |
Allele4: -1bp |
WT |
59.07±2.75** |
JH-C32 |
Allele2: +1bp |
WT |
WT |
Allele4: -1bp |
WT |
60.55±3.25** |
JH-C34 |
Allele5: -5bp |
Allele5: -5bp |
WT |
Allele5: -1bp |
WT |
63.55±2.75** |
JH-C37 |
Allele3: +1bp |
Allele3: -3bp |
WT |
Allele2: +1bp |
WT |
59.55±3.25** |
JH-C38 |
Allele2: +1bp |
Allele8: -45bp |
WT |
Allele9: -19bp |
WT |
55.5±3.75** |
JH-C40 |
Allele4: +1bp |
Allele3: -3bp |
WT |
Allele4: -1bp |
WT |
60.55±2.75** |