Next Article in Journal
Identification of the Actin-Binding Region and Binding to Host Plant Apple Actin of Immunodominant Transmembrane Protein of ‘Candidatus Phytoplasma mali’
Previous Article in Journal
Posttranscriptional Regulation of the Plasminogen Activation System by Non-Coding RNA in Cancer
Previous Article in Special Issue
Integration of QTL Mapping and Whole Genome Sequencing Identifies Candidate Genes for Alkalinity Tolerance in Rice (Oryza sativa)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fine Mapping and Cloning of a qRA2 Affect the Ratooning Ability in Rice (Oryza sativa L.)

Rice Research Institute, Fujian Academy of Agricultural Sciences, Fujian High Quality Rice Research and Development Center, Fuzhou 350018, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(2), 967; https://doi.org/10.3390/ijms24020967
Submission received: 15 November 2022 / Revised: 17 December 2022 / Accepted: 20 December 2022 / Published: 4 January 2023
(This article belongs to the Special Issue Research Advances in Rice Breeding and Genetics)

Abstract

:
Ratooning ability is a key factor that influences the ratoon rice yield in areas where light and temperature are not sufficient for second-season rice. Near-isogenic lines (NILs) are the most powerful tools for the detection and precise mapping of quantitative trait loci (QTLs). In this study, using 176 NILs, we identified a novel QTL for ratooning ability in NIL128. First, we mapped the QTL between the markers Indel12-29 and Indel12-31, which encompass a region of 233 kb. The rice genome annotation indicated the existence of three candidate genes in this region that may be related to ratooning ability. Through gene prediction and cDNA sequencing, we speculated that the target gene of ratooning ability is LOC_Os02g51930 which encodes cytokinin glucosyl transferases (CGTs), hereafter named qRA2. Further analysis showed that qra2 was a 1-bp substitution in the first exon in NIL128, which resulted in the premature termination of qRA2. The results of the knockdown experiment showed that the Jiafuzhan knockdown mutants exhibited the ratooning ability phenotype of NIL128. Interestingly, the qRA2 gene was found to improve ratooning ability without affecting major agronomic traits. These results will help us better understand the genetic basis of rice ratooning ability and provide a valuable gene resource for breeding strong ratoon rice varieties.

1. Introduction

Rice (Oryza sativa L.) is one of the major staple foods worldwide, feeding more than 65 percent of China’s population [1]. Maintaining a sustained increase in rice production is particularly important to meet the needs of a growing population, thereby ensuring national food security [2]. However, due to the rapid growth of the world’s population, rice production will have to increase significantly in the next half century to meet the needs of world population growth. Obviously, this increase in rice production requires higher yields and more frequent harvests on existing land [3]. This is mainly because it is difficult to increase rice yields further by expanding the area under cultivation [4].
Ratoon rice, related to ratooning ability, is a second crop formed by the regeneration of tillers from the remaining nodal buds of the stubble after the first harvest, thus increasing the total yield by increasing the frequency of the harvests [5,6]. Compared with other rice planting modes in China, the annual production yield of ratoon rice is significantly higher than that of single-season rice, while the net economic benefit is higher than that of double-season rice [7]. Therefore, the rapid development of regrown rice in China and around the world, both in applied research and basic research, has become a worldwide research hotspot in recent years [8,9,10].
To date, few studies have been reported on the ratooning ability of rice, and its genetic basis is still unclear. Six QTLs for ratooning ability were detected in a doubled haploid population [11]. Using advanced backcross lines, a QTL called qRAT5 was detected between InDel3 and RM249 on chromosome five [12]. One QTL for ratooning ability was identified using an SSR marker [13]. Two QTLs for ratooning ability were identified by a recombinant inbred line population [14]. Three QTLs for ratooning ability were detected by a recombinant inbred line population [15]. Recently, a QTL named qRA5 was located between 14.12 Mb and 14.65 Mb with a distance of 3.1 Mb using an introgression line population [16].
The genes encoding cytokinin glycosyl transferase (CGT) and their potential applications in crop yields have received little attention. Gene coding for a xylosyl transferase was first isolated from Phaseolus lunatus [17]. Li et al. [18] identified the first cytokinin glycosyltransferase Os6 from rice, and further research revealed that Os6 was expressed in both the vegetative and reproductive tissues of rice. Recent studies have shown that CGT is a key regulator of cytokinin homeostasis and has potential value for wheat genetic improvement [19].
Therefore, it is the most effective way to improving the yield of ratoon rice will require better understanding of the specific functions of ratooning ability-related genes. The strength of ratooning ability plays an important role in the stable and high-yield cultivation of ratoon rice, however, the identification and genetic mechanism of ratooning ability-related genes remain unclear. Consequently, further studies are needed to elucidate the genetic mechanism of ratoon rice formation, cultivate high-ratoon rice varieties, and promote the development of the ratoon rice industry by optimizing cultivation techniques. In this study, Hui1586, a japonica cultivar, was used as the donor parent, and an indica cultivar, Jiafuzhan, was used as the recipient parent. Ratooning ability was evaluated using 176 NILs. The gene qRA2, which affects the ratoon ability-related traits, was validated using NIL128, and its genetic background effect was evaluated. Furthermore, qRA2 was fine mapped within 233 kb by a secondary population, and LOC_Os02g51930 was identified as a candidate gene. Using clustered regularly interspaced short palindromic repeats/CRISPR-associated nuclease 9 (CRISPR/Cas9) genome editing technology, we knocked out the qRA2 gene in the background of Jiafuzhan and observed that the knockout mutants exhibited the ratooning ability phenotype of NIL128. Interestingly, the qRA2 gene was found to improve ratooning ability without affecting major agronomic traits. These results indicate that qRA2 has good application prospects in ratoon rice breeding in the future.

2. Results

2.1. Identification and Analysis of QTLs for Ratooning Ability in the NILs

To evaluate the potential advantages of the NILs for QTL detection, phenotypic variations in ratooning ability were observed in 176 NILs, Jiafuzhan, and Hui1586 at the Fuzhou Experimental Station in July 2019 (Fuzhou, Fujian, China), and it was found that NIL128 exhibited a lower ratooning ability than Jiafuzhan (Figure 1). There was no significant difference between the ratooning ability of the remaining 175 NILs and that of Jiafuzhan.
We further analysed the reasons for the changes in the ratooning ability of NIL128. WinQTLCart 2.5 software was used to analyse the relevant QTLs. The results showed that the LOD score of the QTL for the ratooning ability of NIL128 was 28.26, and the explained phenotypic variance (R2) was 17.82%. Further analysis showed that the ratooning abilities of Jiafuzhan and Hui1586 were 137.5% and 52.3%, respectively, while that of NIL128 was 48.2%. According to the t-test, the differences in the ratooning ability of these strains reached a significant level (Table 1; Figure 1).
Phenotypic comparisons between NIL128 and Jiafuzhan were performed at the Fuzhou experimental station in June 2021 (Fuzhou, Fujian, China) (Table 1). Apart from the ratooning ability, the only two agronomic traits that showed differences were the seed set rate and yield per plant between the NIL128 line and the Jiafuzhan cultivar. However, there were no significant differences in the remaining agronomic traits (Table 1).

2.2. Genetic Analysis of NIL128 for Ratooning Ability

To conduct a genetic analysis of the ratooning ability of NIL128, NIL128 was hybridized with Jiafuzhan. The results showed that F1 hybridization exhibited the phenotype of Jiafuzhan of ratooning ability, and the F2 population phenotype was consistent with the Mendelian separation ratio (Table 2). Segregation between the Jiafuzhan and NIL128 phenotypes fit a 3:1 segregation ratio in the two F2 populations (χ2 = 0.134–0.456, p > 0.5). The results showed that the regeneration ability phenotype of NIL128 was inherited by invisible traits.

2.3. Analysis of the Genetic Background of NIL128

The genetic background of NIL128 homozygous material was analysed using 296 primers that showed polymorphism between Hui1586 and Jiafuzhan cultivars [20]. The results showed that five markers exhibited homozygous alleles of Hui 1586, including Ind1-3 on chromosome 1, RM7286 on chromosome 2, Ind3-18 on chromosome 3, RM566 on chromosome 9, and RM19 on chromosome 12. The remaining 291 markers showed the genetic backgrounds of the Jafuzhan variety (Supplementary Figure S1). Therefore, we speculate that the NIL128 phenotype for ratooning ability may originate from these five loci.

2.4. Linkage Analysis of QTLs for Ratooning Ability in NIL128

To determine which gene causes the NIL128 phenotype, we located the QTL of NIL128 for ratooning ability based on five potential candidate sites. A total of 325 recessive plants from the F2 population were used for linkage analysis between five markers and the QTL. One of five SSR markers, RM3795, was found to be linked to the trait among the 325 F2 individuals.

2.5. Initial Localization of QTLs for Ratooning Ability

To initially locate the QTL, published markers were used around RM3795. The QTL for ratooning ability was located between the molecular markers RM1386 and RM7286 at a distance of 3.1 Mb by linkage labeling (Figure 2a), and each molecular marker corresponding to the physical distance obtained according to each molecular marker sequence.

2.6. Fine Mapping of QTLs for Ratooning Ability

To fine map, the QTL, fourteen polymorphic markers were selected from 42 newly developed molecular markers (Table 3). Six polymorphic markers were used to genotype all recombinant strains. The results showed that the QTL was located between the molecular markers Indle2-20 and Indle2-24 on chromosome 2 at a physical distance of 233 kb (Figure 2b and Table 3).

2.7. Candidate Genes in the 233-kb Region

According to the gene annotation database (http://rice.plantbiology.msu.edu/; http://www.tigr.org/ and accessed on 13 November 2020), there were fifty-one annotated genes in the 233-kb region (Figure 2c), and almost every gene had a corresponding full-length cDNA. After further analysis, we found that three genes were speculated to be related to ratooning ability: LOC_Os02g51900, LOC_Os02g51910, and LOC_Os02g51930. All three genes encode cytokinin-O-glucosyltransferase 2.

2.8. Sequence Analyses of the QTL for Ratooning Ability

To investigate which gene was responsible for the phenotype, we then sequenced three candidate genes of Jafuchian and NIL128. The results showed that a 2-bp substitution (493: G⇨T and 1297: T⇨G) was found in LOC_Os02g51930 (Figure 2c), and there was no difference in the sequence of the remaining two genes. Therefore, we speculate that qra2 was identified as having a weak ratooning ability in NIL128, while LOC_Os02g51930 corresponds to qRA2.
Analysis of the structure of open reading fragments (ORFs) showed that qRA2 had two exons; compared with the reference sequence Nipponbare (Nip), there was a base substitution (1297: T⇨G) at the second exon that caused the amino acid to change from C to G (Supplementary Figures S2 and S3). Further analysis revealed that qra2 was a 1-bp substitution (493: G⇨T) in the 493 rd base of the first exon, which resulted in the premature termination of qRA2 (Supplementary Figures S2 and S3).

2.9. qra2 Is Responsible for the Ratooning Ability Phenotype of NIL128

To confirm whether qra2 is associated with ratooning ability, we investigated whether knocking out qRA2 in Jiafuzhan results in the ratooning ability phenotype of NIL128. A sequence-specific guided RNA (sgRNA) knockout of the qRA2 gene was designed using the CRISPR/Cas9 gene editing system. Two homozygous lines (qPA2KO-line1 and qPA2KO-line2) were obtained from two independent events and confirmed by sequencing the insertion or deletion of their main carrying targets in December 2020 (Table 4).
In April 2021, qRA2KO-line1, qRA2KO-line1, Jiafuzhan, Hui1586, and NIL128 were planted at the Fuzhou Experimental Station. We then investigated the ratooning ability phenotypes of these two homozygous lines (qPA2KO-line1 and qPA2KO-line2) in the first season and ratooning season and found that both lines showed the phenotype of NIL128 in the ratooning season. In addition, there was no significant difference between the two knockout lines (qPA2KO-line1 and qPA2KO-line2) and Jiafuzhan in terms of the main agronomic characteristics, including plant height, panicle length, number of effective panicles, number of spikelets per panicle, seed setting rate, 1000-grain weight, grain length, grain width, and yield per plant (Figure 3 and Table 1). Therefore, we speculated that the loss of qRA2 function was the main reason leading to the weak ratooning ability phenotype of NIL128, and qRA2 may not affect the main agronomic traits of rice.

2.10. Tissue Expression Analysis of qRA2

To further understand the function of qRA2, quantitative reverse transcription-PCR (qRT−PCR) was used to detect the expression patterns of qRA2 at different developmental stages of rice. The results showed that qRA2 was expressed in all tissues tested here, including roots, shoots, and leaves of two- and four-week-old panicles of 0.5–1 cm, 1–3 cm, 3–5 cm, and 5–10 cm length, along with germinating and mature seeds and callus, however, it was predominantly expressed in the seed (germination and mature) and four-week-old shoots (Figure 4).

3. Discussion

3.1. qRA2, a New Gene Associated with Ratooning Ability

In areas where light and temperature are not suitable for the growth of second-season rice, it is important to expand the planting area of reproducing rice. However, the ratooning ability is a complex quantitative trait that is controlled by multiple genes and greatly influenced by the environment. NILs have been used in plant genetic studies, QTL isolation, and fine mapping of genome-wide target traits [21]. The main feature of the NIL population is that each NIL contains only one or more donor chromosome fragments, which has a unique advantage for QTL identification, and by masking background genetic noise, QTLs can be visualized as a single Mendelian factor. Using the population of NILs, several QTLs have been identified or cloned [22,23,24,25,26,27,28]. In this study, we constructed 176 NILs using Jiafuzhan as the genetic background. Using these 176 NILs, we mapped the ratooning ability QTL and designated it qRA2. Map-based cloning and knockout experiments confirmed that the weak ratooning ability of NIL128 was caused by the loss of CGT function.
CGT is significantly related to plant responses to internal and external environments. From the perspective of the application, progress has been made in controlling related agronomic traits by regulating the level of cytokinin within agronomic species [19]. It was found that in rice when the expression of glucose transferase was reduced by RNAi, both grain number/panicle number and grain size were increased [29]. In addition, in a study of wheat, it was found that among the five wheat varieties with different yields, the varieties with higher yields had a lower glucose transferase expression [30]. Recent studies have shown that CGT is a key regulator of cytokinin homeostasis and has potential value for wheat genetic improvement [19]. Our study revealed that NIL128, compared with Jiafuzhan, had a lower seed setting rate, resulting in a decrease in yield. However, we did not find any effect on the yield in the two knockout lines of Jiafuzhan (Table 1). We speculated that this was probably due to the different genetic backgrounds between NIL128 and the knockout lines, and there were five differences in genetic background analysis between NIL128 and Jiafuzhan (Supplementary Figure S1). Although our study revealed that knockout of CGT did not affect rice yield traits, it remains to be further verified whether changes in qRA2 expression levels and overexpression affect rice yield-related traits.

3.2. The Application Prospect of the qRA2 Gene in Rice Breeding

Breeding is a combination of science and art. In the breeding process, we need to accumulate good genes together. However, we have some traits that are difficult to evaluate directly in the field and require molecular techniques, such as quality, resistance, and ratooning ability. At the same time, some genes are also associated with undesirable agronomic traits, so it is difficult to use them directly. Studies have shown that some resistance R genes are associated with undesirable agronomic traits, often at the cost of “sacrificing” yield, making it impossible to achieve both yield and resistance [31,32]. In this study, our data indicated that the qRA2 gene not only enhanced the ratooning ability of rice but also did not affect the important agronomic traits of rice (Table 1 and Figure 3). Meanwhile, the identification of QTLs related to ratooning ability and the introduction of QTLs into excellent rice varieties are effective ways to improve rice viability. Together, the results indicate that qRA2 has good application prospects for breeding strong ratooning ability varieties in the future. For example, to improve the ratooning ability of rice varieties, qRA2 can be transferred into different rice materials.

4. Materials and Methods

4.1. Plant Material

The indica rice cultivar Jiafuzhan and the japonica rice cultivar Hui1586 were preserved at the Rice Research Institute, Fujian Academy of Agricultural Sciences, China. In the early stage, 176 NILs were constructed with the indica cultivar Jiafuzhan as the acceptor and the japonica cultivar Hui1586 as the donor [20].

4.2. Investigation of Ratooning Ability

In this study, the final ratooning rate was used to express the ratooning ability, and the final ratooning rate (%) = the number of effective ears of reproducible rice/the number of effective ears of full-season rice × 100. Five effective panicles were randomly sampled in the middle of each plot, and the average of three replicates was used as the investigation result [33].

4.3. QTL Analysis of Ratooning Ability

In this study, WinQTLCart 2.5 software (http://statgen.ncsu.edu/qtlcart/WQTLCart.htm) [34] was used to achieve composite interval mapping (CIM) for QTL detection [35]. The confidence interval was defined as the 1-LOD reduction region near the peak LOD site of the identified QTL, and LOD ≥ 2.5 was used as the threshold for QTL identification.

4.4. Identification of QTLs for Ratooning Ability

In April 2019, Jiafuzhan, Hui 1586, and 176 NILs were planted in paddy fields under natural conditions at the experimental farm of the Fujian Academy of Agricultural Sciences (Fuzhou, Fujian, China). For parents and each NIL, a total of 64 plants were planted in 6 rows. All plants were grown using standard commercial methods, with row spacings of 13.3 cm and 26.4 cm, and field management techniques followed conventional planting methods.
Five plants were selected from the middle of each plot to study the characteristics of ratooning ability, and QTLs were determined according to the significance of differences between parents. In addition, the related agronomic characteristics, mainly plant height, panicle length, number of effective panicles, number of spikelets per panicle, seed setting rate, 1000-grain weight, grain length, grain width, and yield per plant, were measured at the maturity stage.

4.5. Construction of the Mapping Population

NIL128 was hybridized with Jafuzhan to construct a localization population at the Fuzhou Experimental Station in July 2019 (Fuzhou, Fujian, China), and the F1 plants were grown at the Sanya Experimental Station in November 2019 (Sanya, Hainan, China). The F2 localization population was formed by F1 self-crossing, and all F2 plants were planted at the Fuzhou Experimental Station in April 2020 (Fuzhou, Fujian, China). The primary linkage of ratooning ability QTLs was obtained from 45 recessive plants in the F2 population, and 325 recessive plants in the F2 population were selected for fine localization.

4.6. Molecular Marker Amplification and Detection

DNA was extracted from rice leaves by the CTAB method and modified appropriately [36]. For PCR amplification, every 20-μL reaction mixture contained 30 ng DNA (2-μL), 0.4 μM of each primer (1-μL), 2× Es Tag MasterMix (Dye) (10-μL), and 6-μL sterile deionized water [20]. The amplified PCR products were stained by 3% agarose gel electrophoresis with ethidium bromide staining as the marker [37].

4.7. Molecular Mapping of QTLs for Ratooning Ability

The bands of Hui 1586 and Jiafuzhan were denoted as 1 and 3, respectively, and the heterozygotes were denoted as 2. MAPMAKER version 3.0 software was used to analyse the association between QTL loci and molecular markers. MapDraw V2.1 was used to estimate genetic distance [38], and the linkage map construction method was the same as that previously reported [39]. The common primers used near RM3795 are references to this database (http://archive.gramene.org/markers/ and accessed on 12 August 2020), and Indel markers were designed on the basis of published rice genome sequences, and the polymorphisms between NIL128 and Jiafuzhan were predicted by comparing the sequences of Nipponbare (http://rgp.dna.affrc.go.jp/ and accessed on 1 September 2020) and 93-11 (http://rice.genomics.org.cn/ and accessed on 11 September 2020). Possible polymorphic primers were designed according to the different sequences of the target region sequences, and polymorphisms between NIL128 and Jiafuzhan were further detected.
In order to obtain each molecular marker in rice chromosome physical distance, we first opened the database (http://plants.ensembl.org and accessed on 3 September 2021), clicked the “BLAST” program, pasted the preprimer or postprimer sequence of primers on both sides of a marker into the box of “Enter the Query Sequence”, selected the species “Rice”, and clicked “Run” for retrieval. The specific position of the molecular marker on the physical map of rice can be obtained.

4.8. Physical Map Construction and Bioinformatics Analysis of QTLs for Ratooning Ability

The physical map of QTLs for ratooning ability was constructed based on the information of the International Rice Genome Sequencing Project (IRGSP, http://rgp.dna.affrc.go.jp/IRGSP/index.html and accessed on 9 October 2020). The clones were anchored with the target gene-linked markers, and sequence alignment was then carried out using the pairwise Basic Local Alignment Search Tool (http://www.ncbi.nlm.nih.gov/blast/bl2seq/b12.html and accessed on 9 October 2020). The candidate genes were predicted according to the existing sequence annotation database (http://rice.plantbiology.msu.edu/; http://www.tigr.org/ and accessed on 13 November 2020).

4.9. Targeted Mutagenesis of qRA2 in Rice with CRISPR/Cas9

We used the CRISPR plant database and website to design highly specific gRNA sequences at 267 bp on the first exon of the qRA2 gene (Supplementary Table S1) [40]. Deletion and insertion of target genes were detected by PCR technology, and PCR amplification products were sequenced from the transgenic CRISPR-edited lines to determine the location and type of target gene mutations [41]. Primers for the CRISPR/Cas9 study are listed in Supplementary Table S1. Sequence primers for qRA2 gene cDS amplification are shown in Supplementary Table S1.

5. Conclusions

In this study, a novel QTL for ratooning ability was mapped to a 233-kb locus region. Then, CRISPR/Cas9 knockout experiments confirmed that qra2 is responsible for the ratooning ability phenotype of NIL128. Importantly, the qRA2 gene, encoding CGT, was found to improve ratooning ability without affecting major agronomic traits, indicative of potential applications of qRA2 in rice breeding programs. Taken together, our results will help us better understand the genetic basis of rice ratooning ability and provide a valuable gene resource for breeding strong ratoon rice varieties.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24020967/s1.

Author Contributions

All authors contributed to the study’s conception and design. D.Y. planned and carried out the experiments, data collection, and wrote the manuscript with input from all authors. N.H. and F.H. were involved in conducting experiments, data collection, and analyses. All authors discussed the results and contributed to the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the Special Fund for Agro-scientific Research in the Public Interest of Fujian Province (No. 2023R11010016-3), the Fujian Provincial Natural Science Foundation (No. 2022J01450 and 2021J01471), the Upgrading and Construction of the Platform of Fuzhou Branch Center of National Rice Improvement Center (No. CXPT20230003), Major Science and Technology Projects of Fujian Province (No. 2020NZ08016), Science and Technology Innovation Team (No. CXTD2021001 and CXTD2021005-3), 5511 Collaborative Engineering Project (No. KXXYJBG0021), and the 100 Talent Plans of Fujian Province.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated during this study are included in this published article and its Supplementary Materials files, and the raw data used or analyzed during the current study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, M.; Zou, Y. Integrating mechanization with agronomy and breeding for food security in China. Field Crops Res. 2018, 224, 22–27. [Google Scholar] [CrossRef]
  2. Peng, S.; Tang, Q.; Zou, Y. Current status and challenges of rice production in China. Plant Prod. Sci. 2009, 12, 3–8. [Google Scholar] [CrossRef] [Green Version]
  3. Ray, D.K.; Foley, J.A. Increasing global crop harvest frequency: Recent trends and future directions. Environ. Res. Lett. 2013, 8, 044041. [Google Scholar] [CrossRef]
  4. Deng, N.; Grassini, P.; Yang, H.; Huang, J.; Cassman, K.G.; Peng, S. Closing yield gaps for rice self-sufficiency in China. Nat. Commun. 2019, 10, 1725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Harrell, D.L.; Bond, J.A.; Blanche, S. Evaluation of main-crop stubble height on ratoon rice growth and development. Field Crops Res. 2009, 114, 396–403. [Google Scholar] [CrossRef]
  6. Wang, W.; He, A.; Jiang, G.; Sun, H.; Jiang, M.; Man, J.; Cui, K.; Nie, L. Ratoon rice technology: A green and resource-efficient way for rice production. Adv. Agron. 2020, 159, 135–167. [Google Scholar]
  7. Yuan, S.; Cassman, K.G.; Huang, J.; Peng, S.; Grassini, P. Can ratoon cropping improve resource use efficiencies and profitability of rice in central China? Field Crops Res. 2019, 234, 66–72. [Google Scholar] [CrossRef]
  8. Wang, F.; Huang, J.; Peng, S. Research and development of mechanized rice ratooning technology in China. China Rice 2021, 27, 1–6. [Google Scholar]
  9. Yu, X.; Tao, X.; Liao, J.; Liu, S.; Xu, L.; Yuan, S.; Zhang, Z.; Wang, F.; Deng, N.; Huang, J.; et al. Predicting potential cultivation region and paddy area for ratoon rice production in China using Maxent model. Field Crops Res. 2022, 275, 108372. [Google Scholar] [CrossRef]
  10. Zheng, C.; Wang, Y.; Yuan, S.; Xiao, S.; Sun, Y.; Huang, J.; Peng, S. Heavy soil drying during mid-to-late grain filling stage of the main crop to reduce yield loss of the ratoon crop in a mechanized rice ratooning system. Crop J. 2022, 10, 280–285. [Google Scholar] [CrossRef]
  11. Tian, Z.; Shen, L.; Yuan, Z.; Lu, C.; Chen, Y.; Zhou, K.; Zhu, L. Identification of QTLs for ratooning ability and grain yield traits of rice and analysis of thair genetic effectes. Acta Agron. Sin. 1997, 23, 289–295. [Google Scholar]
  12. Ji, S.D.; Luo, X.; Ahn, S.N. Mapping QTL for ratooning ability in advanced backcross lines from an Oryza sativa O. rufipogon cross. Korean J. Agric. Sci. 2014, 41, 1–7. [Google Scholar] [CrossRef]
  13. Zheng, J.; Li, Y.; Lin, W. Identification of QTLs for ratooning ability and grain yield traits in ratoon rice based on SSR marker. Mol. Plant Breed. 2004, 2, 342–347. [Google Scholar]
  14. Yang, C.; Wang, Y.; Tu, B.; Li, T.; Hu, L.; Li, S. QTL analysis of rice ratooning ability and related agronomic traits by using RIL populations. Acta Agron. Sin. 2012, 38, 1240–1246. [Google Scholar] [CrossRef]
  15. Yang, L.; Jiang, K.; Zhang, T.; Zhao, J.; Guo, X.; Yang, Q.; Wan, X.; You, S.; Cao, Y.; Luo, J.; et al. QTL analysis of rice ratooning ability and grain yield traits in a RIL population of rice. Southwest China J. Agri. Sci. 2013, 26, 2179–2184. [Google Scholar]
  16. Hu, H.; Gao, R.; He, L.; Liang, F.; Li, Z.; Xu, J.; Yang, L.; Wang, C.; Liu, Z.; Xu, J.; et al. Genetic dissection of rice ratooning ability using an introgression line population and substitution mapping of a pleiotropic quantitative trait locus qRA5. Plants 2022, 11, 1134. [Google Scholar] [CrossRef]
  17. Martin, R.C.; Mok, M.C.; Mok, D.W. A gene encoding the cytokinin enzyme zeatin O-xylosyltransferase of Phaseolus vulgaris. Plant Physiol. 1999, 120, 553–557. [Google Scholar] [CrossRef] [Green Version]
  18. Li, P.; Lei, K.; Li, Y.; He, X.; Wang, S.; Liu, R.; Ji, L.; Hou, B. Identification and characterization of the first cytokinin glycosyltransferase from rice. Rice 2019, 12, 19. [Google Scholar] [CrossRef]
  19. Chen, L.; Zhao, J.; Song, J.; Jameson, P.E. Cytokinin glucosyl transferases, key regulators of cytokinin homeostasis, have potential value for wheat improvement. Plant Biotechnol. J. 2021, 19, 878–896. [Google Scholar] [CrossRef]
  20. He, N.; Zhan, G.; Huang, F.; Abou-Elwafa, S.F.; Yang, D. Fine mapping and cloning of a major QTL qph12, which simultaneously affects the plant height, panicle length, spikelet number and yield in rice (Oryza sativa L.). Front. Plant Sci. 2022, 13, 878558. [Google Scholar] [CrossRef]
  21. Yang, D.; Ye, X.; Zheng, X.; Cheng, C.; Ye, N.; Huang, F. Development and evaluation of chromosome segment substitution lines carrying overlapping chromosome segments of the whole wild rice genome. Front. Plant Sci. 2016, 7, 01737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ashikari, M.; Sakakibara, H.; Lin, S.; Yamamoto, T.; Takashi, T.; Nishimura, A.; Angeles, E.R.; Qian., Q.; Kitano., H.; Matsuoka., M. Cytokinin oxidase regulates rice grain production. Science 2005, 309, 741–745. [Google Scholar] [CrossRef] [PubMed]
  23. Fan, C.; Xing, Y.; Mao, H.; Lu, T.; Han, B.; Xu, C.; Li, X.; Zhang, Q. GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor. Appl. Genet. 2006, 112, 1164–1171. [Google Scholar] [CrossRef] [PubMed]
  24. Song, X.; Huang, W.; Shi, M.; Zhu, M.; Lin, H. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat. Genet. 2007, 39, 623–630. [Google Scholar] [CrossRef] [PubMed]
  25. Xue, W.; Xing, Y.; Weng, X.; Zhao, Y.; Tang, W.; Wang, L.; Zhou, H.; Yu, S.; Xu, C.; Li, X.; et al. Natural variation in Ghd7 is an important regulator of heading date and yield potential in rice. Nat. Genet. 2008, 40, 761–767. [Google Scholar] [CrossRef]
  26. Furuta, T.; Uehara, K.; Rosalyn, B.; Angeles-Shim, R.B.; Shim, J.; Ashikari, M.; Takashi, T. Development and evaluation of chromosome segment substitution lines (CSSLs) carrying chromosome segments derived from Oryza rufipogon in the genetic background of Oryza sativa L. Breed. Sci. 2014, 63, 468–475. [Google Scholar] [CrossRef] [Green Version]
  27. Henry, A.; Swamy, B.P.M.; Dixit, S.; Torres, R.D.; Batoto, T.C.; Manalili, M.; Anantha, M.S.; Mandal, N.P.; Kumar, A. Physiological mechanisms contributing to the QTL-combination effects on improved performance of IR64 rice NILs under drought. J. Exp. Bot. 2015, 66, 1787–1799. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, D.; Cheng, C.; Zheng, X.; Ye, X.; Ye, N.; Huang, F. Identification and fine mapping of a major QTL, qHD19, that plays pleiotropic roles in regulating the heading date in rice. Mol. Breed. 2020, 40, 234. [Google Scholar] [CrossRef]
  29. Shang, X.; Xie, R.; Tian, H.; Wang, Q.; Guo, F. Putative zeatin O-glucosyltransferase OscZOG1 regulates root and shoot development and formation of agronomic traits in rice. J. Integr. Plant Biol. 2016, 58, 627–641. [Google Scholar] [CrossRef]
  30. Nguyen, H.N.; Perry, L.; Kisiala, A.; Olechowski, H.; Emery, R.J.N. Cytokinin activity during early kernel development corresponds positively with yield potential and later stage ABA accumulation in field-grown wheat (Triticum aestivum L.). Planta 2020, 252, 76. [Google Scholar] [CrossRef]
  31. Nelson, R.; Wiesner-Hanks, T.; Wisser, R.; Balint-Kurti., P. Navigating complexity to breed disease-resistant crops. Nat. Rev. Genet. 2018, 19, 21–33. [Google Scholar] [CrossRef]
  32. Wang, J.; Zhou, L.; Shi, H.; Chern, M.; Yu, H.; Yi, H.; He, M.; Yin, J.; Zhu, X.; Li, Y.; et al. A single transcription factor promotes both yield and immunity in rice. Science 2018, 361, 1026–1028. [Google Scholar] [CrossRef] [Green Version]
  33. Yang, C.; Wang, Y.; Tu, B.; Zeng, Y.; Li, S. QTL analysis of cold resistance and ratooning ability by using indica-japonica RIL populations in rice. Chin. J. Rice Sci. 2012, 26, 741–745. [Google Scholar]
  34. Wang, S.; Basten, C.J.; Gaffney, P.; Zeng, Z. Windows QTL Cartographer Version 2.5; North Carolina State University Bioinformatics Research Center: Raleigh, NC, USA, 2012. [Google Scholar]
  35. Zeng, Z. Theoretical basis for separation of multiple linked gene effects in mapping quantitative trait loci. PNAS 1993, 90, 10972–10976. [Google Scholar] [CrossRef] [Green Version]
  36. Murray, M.G.; Thompson, W.F. Rapid isolation of high molecular weight plant DNA. Nucleic. Acids Res. 1980, 8, 4321–4325. [Google Scholar] [CrossRef] [Green Version]
  37. Panaud, O.; Chen, X.; Mccouch, S.R. Development of microsatellite and characterization of simple sequence length polymorphism (SSLP) in rice (Oryza sativa L.). Mol. Gen. Genet. 1996, 252, 597–607. [Google Scholar] [CrossRef]
  38. Liu, R.H.; Meng, J.L. MapDraw: A microsoft excel macro for drawing genetic linkage maps based on given genetic linkage data. Hereditas 2003, 25, 317–321. [Google Scholar]
  39. Rahman, M.L.; Chu, S.; Choi, M.; Qiao, Y.; Jiang, W.; Piao, R.; Khanam, S.; Cho, Y.; Jeung, J.; Jena, K.; et al. Identification of QTLs for some agronomic traits in rice using an introgression line from Oryaza minuta. Mo. Cell. 2007, 24, 16–26. [Google Scholar]
  40. Xie, K.; Zhang, J.; Yang, Y. Genome-wide prediction of highly specific guide RNA spacers for CRISPR-Cas9-mediated genome editing in model plants and major crops. Mol. Plant 2014, 7, 923–926. [Google Scholar] [CrossRef] [Green Version]
  41. Ma, X.; Chen, L.; Zhu, Q.; Chen, Y.; Liu, Y. Rapid decoding of sequence-specific nuclease-induced heterozygous and biallelic mutations by direct sequencing of PCR products. Mol. Plant 2015, 8, 1285–1287. [Google Scholar] [CrossRef]
Figure 1. Ratooning ability comparison of Jiafuzhan, Hui1586, and NIL128. The phenotypes of Jiafuzhan, Hui1586, and NIL128 during the ratooning period are shown.
Figure 1. Ratooning ability comparison of Jiafuzhan, Hui1586, and NIL128. The phenotypes of Jiafuzhan, Hui1586, and NIL128 during the ratooning period are shown.
Ijms 24 00967 g001
Figure 2. Physical maps and structural comparison of qRA2. (a): Primary mapping of qRA2. The QTL was mapped to the region between markers RM1386 and RM7286. (b): Further mapping of qRA2, which was localized at a 233-kb region between the markers Indel2-20 and Indel2-24, and the recombinant number between the markers. In addition, the GenBack was marked below each primer marker. (c): qRA2 has two exons; compared with Nip, qRA2 showed a 1-bp substitution in the second exon, and qra2 showed a 1-bp substitution in the first and second exons.
Figure 2. Physical maps and structural comparison of qRA2. (a): Primary mapping of qRA2. The QTL was mapped to the region between markers RM1386 and RM7286. (b): Further mapping of qRA2, which was localized at a 233-kb region between the markers Indel2-20 and Indel2-24, and the recombinant number between the markers. In addition, the GenBack was marked below each primer marker. (c): qRA2 has two exons; compared with Nip, qRA2 showed a 1-bp substitution in the second exon, and qra2 showed a 1-bp substitution in the first and second exons.
Ijms 24 00967 g002
Figure 3. qRA2knockout lines generated using CRISPR/Cas9 technology and exhibiting the NIL128 phenotype. (a): Two independent events (designated qPA2KOline1 and qPA2KOline2) were generated using the CRISPR/Cas9 system and verified by sequencing. (b): Plant phenotype of Jiafuzhan and two knockout lines. (c): Ratooning ability comparison of Jiafuzhan and two knockout lines.
Figure 3. qRA2knockout lines generated using CRISPR/Cas9 technology and exhibiting the NIL128 phenotype. (a): Two independent events (designated qPA2KOline1 and qPA2KOline2) were generated using the CRISPR/Cas9 system and verified by sequencing. (b): Plant phenotype of Jiafuzhan and two knockout lines. (c): Ratooning ability comparison of Jiafuzhan and two knockout lines.
Ijms 24 00967 g003
Figure 4. Determination of the expression patterns of qRA2 via qRT−PCR. RNA samples were extracted from different tissues of Jiafuzhan, including roots, shoots, and leaves of two-, four- and six-week-old seedlings, spikelets of 0.5–1 cm, 1–3 cm, 3–5 cm, and 5–10 cm in length, germinating and mature seeds, and calli. Data represent the mean and standard deviation of three biological replicates. Three technical replicates for each biological sample were used. The error bar represents the standard deviation (SD) of the value from three independent biological samples.
Figure 4. Determination of the expression patterns of qRA2 via qRT−PCR. RNA samples were extracted from different tissues of Jiafuzhan, including roots, shoots, and leaves of two-, four- and six-week-old seedlings, spikelets of 0.5–1 cm, 1–3 cm, 3–5 cm, and 5–10 cm in length, germinating and mature seeds, and calli. Data represent the mean and standard deviation of three biological replicates. Three technical replicates for each biological sample were used. The error bar represents the standard deviation (SD) of the value from three independent biological samples.
Ijms 24 00967 g004
Table 1. Comparison of the main agronomical traits of Jiafuzhan, NIL128, and the qRA2KO knockout mutant lines.
Table 1. Comparison of the main agronomical traits of Jiafuzhan, NIL128, and the qRA2KO knockout mutant lines.
TraitsJiafuzhanHui1586NIL128qRA2KO-Line1qRA2KO-Line2
Plant height (cm)118.32 ± 3.2186.46 ± 1.94 **119.43 ± 3.32118.36 ± 3.21119.28 ± 3.42
Panicle length (cm)28.12 ± 1.2329.35 ± 1.16 *29.12 ± 1.3127.82 ± 1.2628.22 ± 1.41
Number of effective panicle9.84 ± 1.0410.06 ± 1.189.75 ± 1.0210.01 ± 1.029.92 ± 1.09
Spikelets per panicle163.26 ± 4.56145.46 ± 3.88 *164.12 ± 4.21163.41 ± 4.32164.46 ± 4.16
Seed setting rate (%)92.52 ± 1.1693.74 ± 1.6686.22 ± 1.19 *91.52 ± 1.2693.12 ± 1.13
1000-grain weight (g)23.12 ± 0.5626.12 ± 0.54 *23.22 ± 0.3823.42 ± 0.6123.36 ± 0.52
Grain length (mm)10.84 ± 0.128.38 ± 0.22 *10.72 ± 0.1410.88 ± 0.1510.88 ± 0.14
Grain width (mm)2.62 ± 0.073.74 ± 0.08 **2.64 ± 0.062.62 ± 0.092.60 ± 0.08
Yield per plant (g)34.36 ± 0.9235.82 ± 1.01 *32.04 ± 0.98 *35.06 ± 1.0235.49 ± 0.94
Ratooning ability (%)137.5 ± 5.8552.3 ± 4.92 **49.2 ± 4.83 **46.8 ± 5.22 **47.2 ± 4.96 **
* and ** indicate the significance levels of the differences between Jiafuzhan, NIL128 and knockout lines at p < 0.05 and p < 0.01, respectively. The data was derived from the trial that was performed at the Fuzhou experimental station in June 2021.
Table 2. Segregations of F2 populations crossed by NIL128.
Table 2. Segregations of F2 populations crossed by NIL128.
CrossesF1 Phenotype F2 Populationχ2 (3:1)p
The Phenotype of JiafuzhanThe Phenotype of NIL128Total
Plants
NIL128/JiafuzhanThe phenotype of Jiafuzhan98832513130.336 *>0.9
* Denote the segregation ratio of normal plants to mutant plants complied with 3:1 at 0.05 significant probability level. The data was derived from the trial that was performed at the Sanya Experimental Station in April 2020 and at the Fuzhou Experimental Station in August 2020.
Table 3. Indel and SSR molecular markers used for fine mapping of qra2.
Table 3. Indel and SSR molecular markers used for fine mapping of qra2.
MarkerSequence of Forward PrimerSequence of Reverse Primer
RM7286CAGAACAATTCGACCGCTTCGGCTTGAGAGCGTTTGTAGG
RM1386CTACTCCCTAGTTGGCAGCGTCTCCTGCAGGTACGTGCC
RM3795CATTTGCATGGAGAGGATAGTCATCTTCATTTCATTTCACC
Indle2-1GTCTTGGAATTAAATGCTGCGACCAAGATAAATGACAGGC
Indle2-3TACAGGTTAAGAGGAGGCAACAGTTGCAGGATTTATCTGT
Indle2-7AAGATATATAAACGCGTCGGATTGTCCTAAACGTACTGCC
Indle2-11TCCATTGTCTTTCATCCTCTGTGGTACCGATCAATTTCAG
Indle2-12TCTTTTGACAATTTCCCATTCGCAATGACCTTATCTGATT
Indle2-15AGACGATGCGTAGACAAGATTGTATCGTCGTCTTTTGTTG
Indle2-17TGAGACTAAGCGAAGGTAGGCTCAAGCTTACAATTGACCC
Indle2-20TTTCTCCTGAATGAATTGCTCCATCCTTTCTTCTGGTACA
Indle2-24TCACTACACAAACATACAACGCATGTCAACTAAATGGGTT
Indle2-28TATTCGGCCTCTATATCCAACTTTCTCTAATAAAAAGATA
Indle2-33CTTCTTCGACGACCTGGGCAATTTCCAAGCTCTCTCC
Indle2-35AGATATCGCAAGTTAGGCTGTAGGGCAAAAGTTAAAACGA
Indle2-38TTCGATCCTTTTTAATCACCCACATGTTGTGTTGTAGATG
Indle2-40TTTGAAGCTATATGCGTCTGATGTGTGCCATGATTTACAA
Table 4. Mutation site of three target mutant lines.
Table 4. Mutation site of three target mutant lines.
LineTarget TypeMutation Site
qRA2KO-line1gRNAs GGCAAGGCAGGTCTCCTTGGGGAGCGGCACAGCGAC
(Insert 2 bp)
qRA2KO-line2gRNAsGGCAAGGCAGGTCTCCTTGG----GCGGCACAGCGACGG
(Deletion 1 bp)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

He, N.; Huang, F.; Yang, D. Fine Mapping and Cloning of a qRA2 Affect the Ratooning Ability in Rice (Oryza sativa L.). Int. J. Mol. Sci. 2023, 24, 967. https://doi.org/10.3390/ijms24020967

AMA Style

He N, Huang F, Yang D. Fine Mapping and Cloning of a qRA2 Affect the Ratooning Ability in Rice (Oryza sativa L.). International Journal of Molecular Sciences. 2023; 24(2):967. https://doi.org/10.3390/ijms24020967

Chicago/Turabian Style

He, Niqing, Fenghuang Huang, and Dewei Yang. 2023. "Fine Mapping and Cloning of a qRA2 Affect the Ratooning Ability in Rice (Oryza sativa L.)" International Journal of Molecular Sciences 24, no. 2: 967. https://doi.org/10.3390/ijms24020967

APA Style

He, N., Huang, F., & Yang, D. (2023). Fine Mapping and Cloning of a qRA2 Affect the Ratooning Ability in Rice (Oryza sativa L.). International Journal of Molecular Sciences, 24(2), 967. https://doi.org/10.3390/ijms24020967

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop