Quanya Tan1,2†
Suhong Bu1,2†
Guodong Chen1,2Zhenguang Yan1,2
Zengyuan Chang1,2Haitao Zhu1,2
Weifeng Yang1,2Penglin Zhan1,2Shaojun Lin1,2
Liang Xiong1,2Songliang Chen1,2
Guifu Liu1,2
Zupei Liu1,2Shaokui Wang1,2*
Guiquan Zhang1,2*- 1Guangdong Provincial Key Laboratory of Plant Molecular Breeding, State Key Laboratory for Conservation and Utilization of Subtropical Agro-Bioresources, South China Agricultural University, Guangzhou, China
- 2Guangdong Laboratory for Lingnan Modern Agriculture, South China Agricultural University, Guangzhou, China
Asian cultivated rice is a self-pollinating crop, which has already lost some traits of natural outcrossing in the process of domestication. However, male sterility lines (MSLs) need to have a strong outcrossing ability to produce hybrid seeds by outcrossing with restorer lines of male parents in hybrid rice seed production. Stigma exsertion rate (SER) is a trait related to outcrossing ability. Reconstruction of the high-SER trait is essential in the MSL breeding of rice. In previous studies, we detected eighteen quantitative trait loci (QTLs) for SER from Oryza sativa, Oryza glaberrima, and Oryza glumaepatula using single-segment substitution lines (SSSLs) in the genetic background of Huajingxian 74 (HJX74). In this study, eleven of the QTLs were used to develop pyramiding lines. A total of 29 pyramiding lines with 2–6 QTLs were developed from 10 SSSLs carrying QTLs for SER in the HJX74 genetic background. The results showed that the SER increased with increasing QTLs in the pyramiding lines. The SER in the lines with 5–6 QTLs was as high as wild rice with strong outcrossing ability. The epistasis of additive by additive interaction between QTLs in the pyramiding lines was less-than-additive or negative effect. One QTL, qSER3a-sat, showed minor-effect epistasis and increased higher SER than other QTLs in pyramiding lines. The detection of epistasis of QTLs on SER uncovered the genetic architecture of SER, which provides a basis for using these QTLs to improve SER levels in MSL breeding. The reconstruction of the high-SER trait will help to develop the MSLs with strong outcrossing ability in rice.
Introduction
Every cultivated crop was once wild. In the process of domestication, crop species increased their productivity and narrowed their genetic base. The result is that many crops contain only a small fraction of the genetic variation of their wild relatives (Zamir, 2001). Asian cultivated rice (Oryza sativa) is a self-pollinating crop with <1% natural cross-pollination (Virmani and Athwal, 1973). During domestication, the homozygosity of varieties led to a steady decline in the outcrossing ability of cultivated rice. In contrast, African cultivated rice (Oryza glaberrima), which has been domesticated for a shorter period of time, has more wild characteristics in floral traits and outcrossing ability (Marathi and Jena, 2015; Marathi et al., 2015). Wild Oryza species were more outcrossed due to their larger stigmas, longer styles, higher stigma exsertion, and longer spikelet opening period, with perennial species being more outcrossed than annual species (Oka and Morishima, 1967; Parmar et al., 1979; Marathi and Jena, 2015; Marathi et al., 2015). Therefore, the cultivated rice has lost some traits of natural outcrossing during domestication (Parmar et al., 1979).
Hybrid rice was successfully applied in China in the 1970s. Since then, hybrid rice has been widely planted in China and around the world (Yuan and Virmani, 1988). The hybrid rice exhibits strong heterosis, which plays an important role in improving rice yield (Cheng et al., 2007; Zhang et al., 2021; Zhang, 2022). Since male sterility lines (MSLs) are used to produce hybrid seeds by outcrossing with restorer lines of male parents, the outcrossing ability of MSLs is a key factor in improving hybrid seed yield. Therefore, the improvement of outcrossing ability is an important goal of MSL breeding. Stigma exsertion rate (SER) is an evaluating indicator of outcrossing ability and is controlled by quantitative trait loci (QTLs). Over the past decades, many QTLs controlling SER and related traits from different genetic resources have been located in the rice genome (Marathi and Jena, 2015; Zhou et al., 2017a; Liu et al., 2019; Tan et al., 2022). Some of the QTLs were detected in relative species of the cultivated rice, such as O. glaberrima (Tan et al., 2022), Oryza rufipogon (Xiong et al., 1999; Li et al., 2001; Uga et al., 2003a; Huang et al., 2012; Bakti and Tanaka, 2019; Zou et al., 2020), Oryza barthii and Oryza meridionalis (Zou et al., 2020), Oryza longistaminata (Li et al., 2010), and Oryza glumaepatula (Tan et al., 2020). The detection of QTLs for SER laid a foundation for the improvement of outcrossing ability of MSLs.
Understanding the genetic architecture of traits is a prerequisite for trait reconstruction. Genetic architecture of quantitative traits includes the numbers and genome locations of genes affecting a trait, the magnitude of their effects, and the relative contributions of additive, dominant, and epistatic gene effects (Kroymann and Mitchell-Olds, 2005; Holland, 2007). Epistasis is defined as interactions between genome-wide loci, impacting target traits either additively increasing each other's phenotypic effects or non-additively counteracting the effects of the major loci, via other loci that interact epistatically (Mackay, 2014; Misra et al., 2021). Epistasis complicates the genotype–phenotype relationship because it causes hidden quantitative genetic variation in natural populations and may lead to the small additive effects (Holland, 2007; Mackay, 2014). The role of epistasis in the genetic architecture of quantitative traits is controversial due to the fact that most genetic variation for quantitative traits is additive (Mackay, 2014). In QTL mapping approaches, first-order effects are commonly fitted before second- and higher-order (epistatic) effects, which makes it statistically difficult to detect epistatic effects (Gaertner et al., 2012; Doust et al., 2014). The studies in Arabidopsis and rice suggested that epistatic QTL effects were more important than additive QTL for fitness traits (Malmberg et al., 2005; Mei et al., 2005). By contrast, the studies designed to explicitly model epistatic interactions in maize revealed that epistasis was of less or only moderate importance for quantitative traits (Schon et al., 2004; Mihaljevic et al., 2005; Blanc et al., 2006; Holland, 2007).
Like near isogenic lines (NILs), single-segment substitution lines (SSSLs) carry only one chromosome substitution segment from donors in a recipient genetic background (Zhang et al., 2004; Keurentjes et al., 2007; Zhang, 2021). We have developed an SSSL library, which includes 2,360 SSSLs derived from 43 donors of 7 species of rice AA genome in the Huajingxian 74 (HJX74) genetic background (Zhang et al., 2004; Xi et al., 2006; Zhang, 2021). The HJX74-SSSLs were widely used to detect QTLs for traits of agronomic importance (Zhang et al., 2012; Zhu et al., 2014, 2018; Yang et al., 2016, 2021a,b; Zhou et al., 2017b; Pan et al., 2021; Zhan et al., 2021; Fu et al., 2022), to clone genes of functional importance, and to mine alleles of functional variants (Zeng et al., 2006; Teng et al., 2012; Wang et al., 2012; Fang et al., 2019; Sui et al., 2019; Gao et al., 2021; Zhan et al., 2022). Using the HJX74-SSSL library as a platform for rice design, a series of cytoplasmic male sterility (CMS), maintainer and restorer lines, and wide-compatible indica lines (WCILs) were developed (Dai et al., 2015, 2016; Luan et al., 2019; Guo et al., 2022). Recently, the SSSLs were used to detect QTLs controlling SER in rice. Eighteen QTLs for SER from O. sativa, O. glaberrima, and O. glumaepatula were detected in the HJX74-SSSLs (Tan et al., 2020, 2021, 2022). In this study, eleven of the QTLs were used to develop 2- to 6-QTL pyramiding lines in the HJX74 genetic background. The results showed that the SER increased with increasing QTLs in the pyramiding lines. The epistasis of additive by additive interaction between QTLs in 2- to 6-QTL lines was less-than-additive or negative effect. The pyramiding lines carrying 5–6 QTLs for SER showed as high SER as wild rice. The reconstruction of the high SER trait will help to develop MSLs with strong outcrossing ability in rice.
Materials and Methods
Plant Materials
In total, eleven of the QTLs for SER were used to develop QTL-pyramiding lines. Among the 11 QTLs, four were from O. sativa (Tan et al., 2021), four were from O. glaberrima (Tan et al., 2022), and three were from O. glumaepatula (Tan et al., 2020). The QTLs were located in the substitution segments of SSSLs in the HJX74 genetic background. In total, twenty-one SSSLs carrying the 11 QTLs were used in this study, 10 of which were used as parents to develop pyramiding lines (Supplementary Table 1, Supplementary Figure 1).
Field Experiments
All plant materials were planted in the farm of South China Agricultural University, Guangzhou (23°07′N, 113°15′E). The materials were planted in 2015–2020, two cropping seasons per year. The first cropping season (FCS) was from late February to mid-July, and the second cropping season (SCS) was from late July to mid-November. Germinated seeds were sown in a seedling bed, and seedlings were transplanted to the paddy field in a single seedling. Field cultivation and controlling of diseases and insect pests were followed by conventional methods in the local area.
Genotyping of the QTLs for SER
All markers in the substitution segments carrying QTLs for SER were selected from the previous QTL mapping studies. The substitution segments and the QTLs for SER in SSSLs and pyramiding lines were detected following the previous studies (Tan et al., 2020, 2021, 2022). To distinguish different loci with the same name and to show the allelic origin of the QTLs, the first three letters of the donor species name were added to the suffix of the QTL names (Supplementary Table 1, Supplementary Figure 1).
Phenotyping and Statistical Analysis
The SER and agronomic traits were investigated following the previous studies (Tan et al., 2020, 2021, 2022). For the statistical analysis, data of percentages were converted to the arcsine square root. The student's t-test was used to compare the two sets of data. Dunnett's t-test was used for multiple group comparison with the control group. The least significance range (LSR) was used for the multiple range test among multiple groups (Duncan, 1955). SPSS statistics 23.0 and OriginPro 9.0 were used for the data analysis and figure making (https://www.originlab.com).
Estimation of Additive Effects and Epistatic Effects of QTLs
The additive effect of a QTL for SER was defined as the genotypic value of the homozygous genotype of a QTL in the HJX74 genetic background. Therefore, the additive effect of each locus was the difference in mean values of SER between HJX74 and an SSSL carrying a target QTL, and the additive effects of two- or multiple-locus genotypes were the difference in mean values of SER between HJX74 and pyramiding lines carrying two or multiple QTLs in the HJX74 genetic background. Epistasis of additive by additive interactions among QTLs for SER in 2- to 6-QTL lines was detected using homozygous genotypes. Epistatic effects among QTLs were estimated by the formula, , where i is an epistasis among the pyramided QTLs, Pn is a phenotype of a pyramiding line harboring n of QTLs, P0 is a phenotype of HJX74, ai (1 ≤ i ≤ n) is an additive effect of a single QTL at the ith QTL. Epistatic effects among QTLs in pyramiding lines were tested in Student's t-test under null hypothesis (H0) i = 0.
Results
Genotypes in the Lines With 1- to 6-QTLs for SER
In total, eleven QTLs for SER were used to develop QTL-pyramiding lines. The QTLs were located in the substitution segments of SSSLs with the HJX74 genetic background. Three SSSLs, namely, A35, A42, and A88, each carried two SER-QTLs in a substitution segment, while the other 18 SSSLs all carried one SER-QTL in their substitution segments (Supplementary Table 1, Supplementary Figure 1). The 1-QTL lines (1QLs), which carried 1 QTL for SER in a substitution segment, were mutually crossed to develop 2-QTL lines (2QLs). The 2QLs were crossed with 1QLs to develop 3-QTL lines (3QLs) or crossed with other 2QLs to develop 4-QTL lines (4QLs). In this way, 5-QTL lines (5QLs) and 6-QTL lines (6QLs) were developed (Supplementary Figure 2).
A total of 29 QTL-pyramiding lines were developed from the SSSLs, which carried different combinations of QTLs for SER in the HJX74 genetic background (Table 1, Supplementary Table 2), including five 2QLs (Supplementary Figure 3), ten 3QLs (Supplementary Figure 4), seven 4QLs (Supplementary Figure 5), five 5QLs (Supplementary Figure 6), and two 6QLs (Supplementary Figure 7).
Table 1. SER-QTL combinations in pyramiding lines.
Phenotypes in the Lines With 1- to 6-QTLs for SER
The SER in the lines with 1- to 6-QTLs were investigated in three cropping seasons (Figures 1A–D). The average SER for 1-QTL lines (1QLs) was 48.1%, ranging from 42.5 to 56.5%, which was significantly higher than the 29.2% of HJX74 (Figures 1C,D, Supplementary Table 3). In seven 2QLs including two SSSLs each carrying 2 QTLs in a single substitution segment and five QTL-pyramiding lines, the average SER was 54.6%, ranging from 46.9 to 63.2%, which was 6.5% higher than the SER in 1QLs. In ten 3QLs, the average SER was 67.4%, ranging from 60.0 to 78.1%, which was 12.8% higher than the SER in 2QLs. In seven 4QLs, the average SER was 70.6%, ranging from 63.4 to 76.5%, which was 3.2% higher than the SER in 3QLs. In five 5QLs, the average SER was 77.1%, ranging from 73.3 to 80.9%, which was 6.5% higher than the SER in 4QLs. In two 6QLs, the average SER was 88.8%, ranging from 86.0 to 91.7%, which was 11.7% higher than the SER in 5QLs (Figures 1C,D, Supplementary Table 4). These results showed that SERs increased with the increase of QTLs in pyramiding lines.
Figure 1. Stigma exsertion rate (SER) of 1- to 6-QTL lines. (A) Plant types in the 2QL-2, 4QL-2, and 6QL-2 lines and HJX74. Scale bar, 15 cm. (B) The phenotypes of stigma exsertion in the 2QL-2, 4QL-2, and 6QL-2 lines and HJX74. Scale bar, 1 cm. The red arrows point to the exposed stigma. (C) SER of 1- to 6-QTL lines. HJX74 as a control. Data are shown as the mean ± S.E. of three cropping seasons. Capital letters indicate the significance of differences at the 1% levels. (D) Frequency distributions of SER in 1- to 6-QTL lines.
Phenotypes of eight agronomic traits in 2- to 6-QTL lines were investigated. Most of the traits in these lines, such as heading date, plant height, grain number per panicle, panicle length, and seed setting rate, were not significantly different from those of HJX74. In grain length, 13 lines were significantly different from HJX74, while the other 18 lines were not significantly different. In grain width, 14 lines were significantly different from HJX74, while the other 17 lines were not significantly different (Supplementary Table 5). These results showed that the pyramiding lines had a similar genetic background to HJX74, except for the genotype of SER. The grain size, including grain length and grain width, had no significant effect on SER.
Additive Effects of QTLs for SER in 1- to 6-QTL Lines
The additive effects of 11 QTLs on SER varied greatly. The average additive effect of 11 QTLs on SER was 18.9%, ranging from 13.3 to 27.3%. Based on the additive effects, 11 QTLs could be divided into three levels. There were 7 QTLs at the low-effect level with the average additive effect of 16.2%, ranging from 13.3 to 17.5%. There were two QTLs at the moderate-effect level with the average additive effect of 20.8%, ranging from 20.2 to 21.4%. At the high-effect level, the additive effects of two QTLs were 25.5 and 27.3%, respectively, and 26.4% on average (Figure 2A, Supplementary Table 6). The 11 QTLs were from three donors. Four QTLs from O. sativa showed the low additive effects with the average additive effect of 16.6%, ranging from 14.8 to 17.5%. The average additive effect of four QTLs from O. glaberrima was 17.1%, ranging from 13.3 to 21.4% with three low-effect QTLs and one moderate-effect QTL. The average additive effect of three QTLs from O. glumaepatula was 24.3%, ranging from 20.2 to 27.3% with two high-effect QTLs and one moderate-effect QTL (Figure 2A). The magnitude of additive effects of QTLs for SER was consistent with the evolutionary levels of species.
Figure 2. Additive effects of QTLs on stigma exsertion rate (SER). (A) Additive effects of 11 QTLs on SER. (B) Additive effects of QTL combinations on SER in 1- to 6-QTL lines. Data are shown as the mean ± S.E. of three cropping seasons. Capital letters indicate the significance of differences at the 1% levels.
Additive effects of 2- to 6-QTLs were calculated from the SER in 2- to 6-QTL lines. The average additive effect of 2-QTL combinations on SER was 25.4% in seven 2QLs, ranging from 17.7 to 34.1%. The average additive effect of 3-QTL combinations on SER was 38.3% in ten 3QLs, ranging from 30.9 to 48.9%. The average additive effect of 4-QTL combinations on SER was 41.4% in seven 4QLs, ranging from 34.3 to 47.3%. The average additive effect of 5-QTL combinations on SER was 47.9% in five 5QLs, ranging from 44.2 to 51.8%. The average additive effect of 6-QTL combinations on SER was 59.7% in two 6QLs, ranging from 56.8 to 62.6% (Figure 2B, Supplementary Table 7). These results showed that the total additive effects of QTLs increased with increasing QTLs in 1- to 6-QTL lines.
Epistasis of Additive by Additive Interaction Among QTLs on SER in 2- to 6-QTL Lines
Although the additive effects of QTLs on SER were increasing in 2- to 6-QTL lines, they were generally below the expected values (Figure 2B, Supplementary Table 4). It suggested that there were epistatic effects of QTLs in the 2- to 6-QTL lines. Among the seven 2QLs, the average epistatic effect per QTL combination was −7.8%, ranging from 0.9 to −14.1%, where the average epistatic effect per QTL was −3.9%, but the epistatic effect was not significant. Among the ten 3QLs, the average epistatic effect per QTL combination was −14.6%, ranging from −3.9 to −21.4%, where the average epistatic effect per QTL was −4.9%, and the epistatic effect was significant at the 0.05 level. Among the seven 4QLs, the average epistatic effect per QTL combination was −25.2%, ranging from −17.0 to −31.2%, where the average epistatic effect per QTL was −6.3%, and the epistatic effect was significant at the 0.01 level. Among the five 5QLs, the average epistatic effect per QTL combination was −40.4%, ranging from −32.4 to −49.0%, where the average epistatic effect per QTL was −8.1%, and the epistatic effect was significant at the 0.001 level. Among the two 6QLs, the average epistatic effect per QTL combination was −47.7%, ranging from −41.7 to −53.7%, where the average epistatic effect per QTL was −8.0%, and the epistatic effect was significant at the 0.001 level (Figures 3A,B, Supplementary Table 8). These results showed that the additive effect of multiple QTLs was smaller than the sum of the additive effects of the corresponding single QTL. The epistatic interactions between QTLs in pyramiding lines reduced the total additive effect of QTL combinations. The epistatic effects of QTLs increased with increasing QTLs in 2- to 6-QTL lines. Therefore, the epistasis was less-than-additive or negative effect.
Figure 3. Epistatic effects of QTLs on stigma exsertion rate (SER). (A) The total epistatic effects of QTL combinations on SER in 2- to 6-QTL lines. (B) The average epistatic effects per QTL on SER in 2- to 6-QTL lines. *P ≤ 0.05; **P ≤ 0.01; ***P ≤ 0.001. (C) Percentage of epistatic effects in genetic effects for each QTL in 2- to 6-QTL lines. Data are shown as the mean ± S.E. of three cropping seasons. Capital letters indicate the significance of differences at the 1% levels.
In 2- to 6-QTL lines, the percentage of epistatic effects in genetic effects for each QTL were 20.7, 26.0, 33.4, 42.9, and 42.4%, respectively (Figure 3C). The results showed that although there was epistasis, additivity was the major component of QTL interaction.
QTL qSER3a-sat Showed a Minor Epistasis Effect
In the 2- to 4-QTL lines, there were two genotypes with and without qSER3a-sat. In 2QLs, the average epistatic effect of three lines with qSER3a-sat was −2.1%, while that of four lines without qSER3a-sat was −12.0%. In 3QLs, the average epistatic effect of six lines with qSER3a-sat was −11.9%, while that of four lines without qSER3a-sat was −18.8%. In 4QLs, the average epistatic effect of four lines with qSER3a-sat was −21.4%, while that of three lines without qSER3a-sat was −30.4%. In the three sets of QTL lines, the average epistatic effects in the lines with qSER3a-sat were significantly less than those without qSER3a-sat (Figure 4A, Supplementary Tables 8, 9). In 2- to 4-QTL lines, the epistatic effects of qSER3a-sat in the qSER3a-sat groups were 3.9, 0.7, and 1.4%, respectively, while the average epistatic effects per QTL in the qSER3a-sat-free groups were −6.0, −6.3, and −7.6%, respectively (Figure 4B, Supplementary Table 9). Differing from other QTLs, the epistatic effects of qSER3a-sat on SER were minor and positive. As a result, the lines with qSER3a-sat had greater SER than those without qSER3a-sat. In 2QLs, the average SER of the lines with qSER3a-sat was 58.6%, while that of the lines without qSER3a-sat was 51.6%. In 3QLs, the average SER of the lines with qSER3a-sat was 71.4%, while that of the lines without qSER3a-sat was 61.5%. In 4QLs, the average SER of the lines with qSER3a-sat was 72.5%, while that of the lines without qSER3a-sat was 68.0% (Figure 4C, Supplementary Table 4). These results showed that qSER3a-sat significantly increased SER than other QTLs in the pyramiding lines due to its minor epistasis effect.
Figure 4. Epistatic effect of qSER3a-sat on stigma exsertion rate (SER). (A) Epistatic effects of the QTLs in the 2- to 4-QTL lines with and without qSER3a-sat. (B) Epistatic effects of a single QTL of qSER3a-sat and other QTLs in 2- to 4-QTL lines. (C) SER in the 2- to 4-QTL lines with and without qSER3a-sat. Data are shown as the mean ± S.E. of three cropping seasons. *P ≤ 0.05; **P ≤ 0.01; ***P ≤ 0.001.
Discussion
The Reduced Additive Effect of QTLs on SER May Be Responsible for the Degradation of Outcrossing Ability in Cultivated Rice
The emergence of cultivated crops is the result of natural and artificial selection. During the domestication of cultivated rice, the selection of varieties led to homozygosity of genotypes. At the same time, genotypic homozygosity reinforced the inbreeding of varieties, leading to a gradual decline in their outcrossing ability. Therefore, the degradation of outcrossing ability is an inevitable outcome of the domestication of cultivated rice. Carter et al. (2005) studied the effects of different forms of epistasis on the response to directional selection using a multilinear epistatic model and showed that the main factor determining the evolution of the additive variance, and thus the evolvability, is directional epistasis. Doust et al. (2014) argued that domestication is a multifaceted evolutionary process, involving changes in individual genes, gene interactions, and emergent phenotypes, and found that alleles favored during domestication tend to have larger phenotypic effects and are relatively insensitive to genetic background and environmental influences compared with wild progenitor alleles. In this study, the additive effects of 11 QTLs could be divided into three levels. The four QTLs from O. sativa, an Asian cultivated rice, were the low additive effects. O. glaberrima, the recently domesticated African cultivated rice, had three low-effect QTLs and one moderate-effect QTL. Among the three QTLs from O. glumaepatula, a wild species, two were high-effect QTLs and one was moderate-effect QTL (Figure 2A). The magnitude of additive effects of QTLs for SER was consistent with the evolutionary levels of species. In the process of domestication, the additive effects of QTLs on SER were reduced in cultivated rice. Therefore, the reduction in the additive effect of QTLs on SER during domestication may be an important factor in the degradation of the outcrossing ability of cultivated rice.
The Epistasis of Additive by Additive Interaction Between QTLs on SER Was Less-Than-Additive
Mapping epistatic interactions is experimentally, statistically, and computationally challenging. The classical populations used to map QTLs are poorly efficient to detect epistasis. The use of secondary mapping populations, such as chromosome segment substitution lines (CSSLs), introgression lines (ILs), and NILs, in which a region containing the QTL is introgressed into the isogenic background of one of the parental lines, and the QTLs are narrowed down to a small genomic interval by recombination in successive generations, facilitates the analysis of epistasis between naturally occurring variants (Causse et al., 2007; Mackay, 2014). Eshed and Zamir (1996) analyzed interactions between individual Lycopersicon pennellii chromosome segments introgressed into an otherwise homogeneous genetic background of L. esculentum in a half diallele scheme. Of the 180 tested interactions, 28% were epistatic on both linear and geometric scales, and the detected epistasis was predominately less-than-additive. Although the frequency of epistasis was high, additivity was the major component in the interaction of pairs of QTLs. Kroymann and Mitchell-Olds (2005) reported that a small chromosome interval with no effect on the growth rate of A. thaliana NILs contained two epistatically interacting QTLs affecting growth, for one of which the effect on growth was in opposite directions in the different genetic backgrounds. Secondary mapping populations were also used to detect epistasis of QTLs in tomato (Causse et al., 2007) and rice (Uwatoko et al., 2008; Chen et al., 2014; Qin et al., 2015; Yang et al., 2018; Misra et al., 2021). In this study, we used pyramiding lines with 2- to 6-QTLs for SER developed by SSSLs to detect epistasis at two or more loci. Due to the SSSLs and pyramiding lines being homozygous genotypes in the HJX74 genetic background, only the epistasis of additive by additive interaction was detected. The results showed that while the epistasis of QTLs increased with increasing QTLs in 2- to 6-QTL lines, the total additive effects of QTLs increased (Figures 2B, 3). These results indicated that the epistatic interactions of the QTLs only changed the magnitude of effects, in which the phenotype of one locus is suppressed by genotypes at the other locus. Therefore, the epistasis of additive by additive interaction among QTLs in 2- to 6-QTL lines was less-than-additive or negative effect. The results also suggest that detecting epistasis of additive-additive interactions in pyramiding lines developed from SSSLs of the same genetic background is a simpler and more effective scheme.
The High-SER Trait Lost in the Domestication of Cultivated Rice Can Be Reconstructed
In the past decades, great progress has been made in the heterosis utilization of hybrid rice (Yuan, 2017; Zhang, 2022). Improving the outcrossing ability of MSLs is very important for hybrid rice seed production. In the past two decades, a larger number of QTLs controlling SER and related traits have been identified on 12 chromosomes of the rice genome from various genetic resources (Marathi and Jena, 2015; Zhou et al., 2017a; Liu et al., 2019). Recently, eighteen QTLs for SER from O. sativa, O. glaberrima, and O. glumaepatula were detected in the SSSLs in the HJX74 genetic background (Tan et al., 2020, 2021, 2022). It laid a foundation for improving the outcrossing ability of MSLs of cultivated rice using the QTLs for SER. Most MSLs and wild rice were found to have the SER of 70% or more (Virmani and Athwal, 1973; Ying and Zhang, 1989; Uga et al., 2003b; Li et al., 2014; Lou et al., 2014; Guo et al., 2017; Zou et al., 2020). In this study, we found that despite the presence of widespread epistasis, additivity remained a major component of QTL interactions, which resulted in an increase in the SER with increasing QTLs in the pyramiding lines. When carrying 5–6 QTLs, the pyramiding lines had the SER level of most MSLs and wild rice (Figure 1). In addition, it was found that qSER3a-sat had a minor epistasis, which was more effective for improving the SER (Figure 4). Our study suggests that it is possible to reconstruct the high-SER trait, which was lost during the domestication of cultivated rice. The understanding of the genetic architecture of SER in rice lays the foundation for reconstructing the high-SER trait in MSL breeding.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
GZ and SW designed and supervised the work. QT and SB performed most of the experiments and prepared the experimental data. GC, ZY, ZC, WY, PZ, SL, LX, and SC conducted a part of the experiments. HZ, GL, and ZL prepared the experimental materials and assisted in the experiments. GZ analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by grants from the Major Program of Transgenic New Variety Breeding of China (2014ZX08009-037B) and from the National Natural Science Foundation of China (91435207 and 91735304).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.921700/full#supplementary-material
References
Bakti, C., and Tanaka, J. (2019). Detection of dominant QTLs for stigma exsertion ratio in rice derived from Oryza rufipogon accession ‘W0120’. Breed. Sci. 69, 143–150. doi: 10.1270/jsbbs.18139
Blanc, G., Charcosset, A., Mangin, B., Gallais, A., and Moreau, L. (2006). Connected populations for detecting quantitative trait loci and testing for epistasis: an application in maize. Theor. Appl. Genet. 113, 206–224. doi: 10.1007/s00122-006-0287-1
Carter, A. J., Hermisson, J., and Hansen, T. F. (2005). The role of epistatic gene interactions in the response to selection and the evolution of evolvability. Theor. Popul. Biol. 68, 179–196. doi: 10.1016/j.tpb.2005.05.002
Causse, M., Chaib, J., Lecomte, L., Buret, M., and Hospital, F. (2007). Both additivity and epistasis control the genetic variation for fruit quality traits in tomato. Theor. Appl. Genet. 115, 429–442. doi: 10.1007/s00122-007-0578-1
Chen, J., Li, X., Cheng, C., Wang, Y., Qin, M., Zhu, H., et al. (2014). Characterization of epistatic interaction of QTLs LH8 and EH3 controlling heading date in rice. Sci. Rep. 4, 4263. doi: 10.1038/srep04263
Cheng, S., Cao, L., Zhuang, J., Chen, S., Zhan, X., Fan, Y., et al. (2007). Super hybrid rice breeding in China: achievements and prospects. J. Integr. Plant Biol. 49, 805–810. doi: 10.1111/j.1744-7909.2007.00514.x
Dai, Z., Lu, Q., Luan, X., Cai, J., Zhu, H., Liu, Z., et al. (2015). Development of a platform for breeding by design of CMS lines based on an SSSL library in rice (Oryza sativa L.). Euphytica 205, 63–72. doi: 10.1007/s10681-015-1384-5
Dai, Z., Lu, Q., Luan, X., Ouyang, L., Guo, J., Liang, J., et al. (2016). Development of a platform for breeding by design of CMS restorer lines based on an SSSL library in rice (Oryza sativa L.). Breed. Sci. 66, 768–775. doi: 10.1270/jsbbs.16044
Doust, A. N., Lukens, L., Olsen, K. M., Mauro-Herrera, M., Meyer, A., and Rogers, K. (2014). Beyond the single gene: how epistasis and gene-by-environment effects influence crop domestication. Proc. Natl. Acad. Sci. U.S.A. 111, 6178–6183. doi: 10.1073/pnas.1308940110
Duncan, D. B. (1955). Multiple range and multiple F tests. Biometrics 11, 1–42. doi: 10.2307/3001478
Eshed, Y., and Zamir, D. (1996). Less-than-additive epistatic interactions of quantitative trait loci in tomato. Genetics 143, 1807–1817. doi: 10.1093/genetics/143.4.1807
Fang, C., Li, L., He, R., Wang, D., Wang, M., Hu, Q., et al. (2019). Identification of S23 causing both interspecific hybrid male sterility and environment-conditioned male sterility in rice. Rice 12, 10. doi: 10.1186/s12284-019-0271-4
Fu, Y., Zhao, H., Huang, J., Zhu, H., Luan, X., Bu, S., et al. (2022). Dynamic analysis of QTLs on plant height with single segment substitution lines in rice. Sci. Rep. 12, 5465. doi: 10.1038/s41598-022-09536-8
Gaertner, B. E., Parmenter, M. D., Rockman, M.V., Kruglyak, L., and Phillips, P. C. (2012). More than the sum of its parts: a complex epistatic network underlies natural variation in thermal preference behavior in Caenorhabditis elegans. Genetics 192, 1533. doi: 10.1534/genetics.112.142877
Gao, M., He, Y., Yin, X., Zhong, X., Yan, B., Wu, Y., et al. (2021). Ca2+ sensor-mediated ROS scavenging suppresses rice immunity and is exploited by a fungal effector. Cell 184, 5391. doi: 10.1016/j.cell.2021.09.009
Guo, J., Li, Y., Xiong, L., Yan, T., Zou, J., Dai, Z., et al. (2022). Development of wide-compatible indica lines (WCILs) by pyramiding multiple neutral alleles of indica-japonica hybrid sterility loci. Front. Plant Sci. 13, 890568. doi: 10.3389/fpls.2022.890568
Guo, L., Qiu, F., Gandhi, H., Kadaru, S., De Asis, E. J., Zhuang, J., et al. (2017). Genome-wide association study of outcrossing in cytoplasmic male sterile lines of rice. Sci. Rep. 7, 3223. doi: 10.1038/s41598-017-03358-9
Holland, J. B. (2007). Genetic architecture of complex traits in plants. Curr. Opin. Plant Biol. 10, 156–161. doi: 10.1016/j.pbi.2007.01.003
Huang, X., Kurata, N., Wei, X., Wang, Z., Wang, A., Zhao, Q., et al. (2012). A map of rice genome variation reveals the origin of cultivated rice. Nature 490, 497–501. doi: 10.1038/nature11532
Keurentjes, J. J. B., Bentsink, L., Alonso-Blanco, C., Hanhart, C. J., Vries, H. B., Effgen, S., et al. (2007). Development of a near-isogenic line population of Arabidopsis thaliana and comparison of mapping power with a recombinant inbred line population. Genetics 175, 891–905. doi: 10.1534/genetics.106.066423
Kroymann, J., and Mitchell-Olds, T. (2005). Epistasis and balanced polymorphism influencing complex trait variation. Nature 435, 95–98. doi: 10.1038/nature03480
Li, C., Sun, C., Mu, P., Chen, L., and Wang, X. (2001). QTL analysis of anther length and ratio of stigma exsertion, two key traits of classification for cultivated rice (Oryza sativa L.) and common wild rice (O. rufipogon Griff.). Acta Genet. Sin. 28, 746–751.
Li, H., Gao, F., Zeng, L., Li, Q., Lu, X., Li, Z., et al. (2010). QTL analysis of rice stigma morphology using an introgression line from Oryza longistaminata L. Mol. Plant Breed. 8, 1082–1089. doi: 10.3969/mpb.008.001082
Li, P., Feng, F., Zhang, Q., Chao, Y., Gao, G., and He, Y. (2014). Genetic mapping and validation of quantitative trait loci for stigma exsertion rate in rice. Mol. Breed. 34, 2131–2138. doi: 10.1007/s11032-014-0168-2
Liu, Y., Zhang, A., Wang, F., Kong, D., Li, M., Bi, J., et al. (2019). Fine mapping a quantitative trait locus, qSER-7, that controls stigma exsertion rate in rice (Oryza sativa L.). Rice 12, 46. doi: 10.1186/s12284-019-0304-z
Lou, J., Yue, G. H., Yang, W. Q., Mei, H. W., Luo, L. J., and Lu, H. J. (2014). Mapping QTLs influencing stigma exertion in rice. Bulg. J. Agric. Sci. 20, 1450–1456.
Luan, X., Dai, Z., Yang, W., Tan, Q., Lu, Q., Guo, J., et al. (2019). Breeding by design of CMS lines on the platform of SSSL library in rice. Mol. Breed. 39, 126. doi: 10.1007/s11032-019-1028-x
Mackay, T. F. C. (2014). Epistasis and quantitative traits: using model organisms to study gene-gene interactions. Nat. Rev. Genet. 15, 22–33. doi: 10.1038/nrg3627
Malmberg, R. L., Held, S., Waits, A., and Mauricio, R. (2005). Epistasis for fitness-related quantitative traits in Arabidopsis thaliana grown in the field and in the greenhouse. Genetics 171, 2013–2027. doi: 10.1534/genetics.105.046078
Marathi, B., and Jena, K. K. (2015). Floral traits to enhance outcrossing for higher hybrid seed production in rice: present status and future prospects. Euphytica 201, 1–14. doi: 10.1007/s10681-014-1251-9
Marathi, B., Ramos, J., Hechanova, S. L., Oane, R. H., and Jena, K. K. (2015). SNP genotyping and characterization of pistil traits revealing a distinct phylogenetic relationship among the species of Oryza. Euphytica 201, 131–148. doi: 10.1007/s10681-014-1213-2
Mei, H., Li, Z., Shu, Q., Guo, L., Wang, Y., Yu, X., et al. (2005). Gene actions of QTLs affecting several agronomic traits resolved in a recombinant inbred rice population and two backcross populations. Theor. Appl. Genet. 110, 649–659. doi: 10.1007/s00122-004-1890-7
Mihaljevic, R., Utz, H. F., and Melchinger, A. E. (2005). No evidence for epistasis in hybrid and per se performance of elite European flint maize inbreds from generation means and QTL analyses. Crop Sci. 45, 2605–2613. doi: 10.2135/cropsci2004.0760
Misra, G., Badoni, S., Parween, S., Singh, R. K., Leung, H., Ladejobi, O., et al. (2021). Genome-wide association coupled gene to gene interaction studies unveil novel epistatic targets among major effect loci impacting rice grain chalkiness. Plant Biotechnol. J. 19, 910–925. doi: 10.1111/pbi.13516
Oka, H., and Morishima, H. (1967). Variations in the breeding systems of a wild rice, Oryza perennis. Evolution 21, 249–258. doi: 10.1111/j.1558-5646.1967.tb00153.x
Pan, Z., Tan, B., Cao, G., Zheng, R., Liu, M., Zeng, R., et al. (2021). Integrative QTL identification, fine mapping and candidate gene analysis of a major locus qLTG3a for seed low-temperature germinability in rice. Rice 14, 103. doi: 10.1186/s12284-021-00544-2
Parmar, K. S., Siddiq, E. A., and Swaminathan, M. S. (1979). Variation in components of flowering behaviour of rice. Indian J. Genet. Plant Breed. 39, 542–550.
Qin, M., Zhao, X., Ru, J., Zhang, G., and Ye, G. (2015). Bigenic epistasis between QTLs for heading date in rice analyzed using single segment substitution lines. Field Crop. Res. 178, 16–25. doi: 10.1016/j.fcr.2015.03.020
Schon, C. C., Utz, H. F., Groh, S., Truberg, B., Openshaw, S., and Melchinger, A. E. (2004). Quantitative trait locus mapping based on resampling in a vast maize testcross experiment and its relevance to quantitative genetics for complex traits. Genetics 167, 485–498. doi: 10.1534/genetics.167.1.485
Sui, F., Zhao, D., Zhu, H., Gong, Y., Tang, Z., Huang, X., et al. (2019). Map-based cloning of a new total loss-of-function allele of OsHMA3 causes high cadmium accumulation in rice grain. J. Exp. Bot. 70, 2857–2871. doi: 10.1093/jxb/erz093
Tan, Q., Wang, C., Luan, X., Zheng, L., Ni, Y., Yang, W., et al. (2021). Dissection of closely linked QTLs controlling stigma exsertion rate in rice by substitution mapping. Theor. Appl. Genet. 134, 1253–1262. doi: 10.1007/s00122-021-03771-9
Tan, Q., Zhu, H., Liu, H., Ni, Y., Wu, S., Luan, X., et al. (2022). Fine mapping of QTLs for stigma exsertion rate from Oryza glaberrima by chromosome segment substitution. Rice Sci. 29, 55–66. doi: 10.1016/j.rsci.2021.12.005
Tan, Q., Zou, T., Zheng, M., Ni, Y., Luan, X., Li, X., et al. (2020). Substitution mapping of the major quantitative trait loci controlling stigma exsertion rate from Oryza glumaepatula. Rice 13, 37. doi: 10.1186/s12284-020-00397-1
Teng, B., Zeng, R., Wang, Y., Liu, Z., Zhang, Z., Zhu, H., et al. (2012). Detection of allelic variation at the Wx locus with single-segment substitution lines in rice (Oryza sativa L.). Mol. Breed. 30, 583–595. doi: 10.1007/s11032-011-9647-x
Uga, Y., Fukuta, Y., Cai, H.W., Iwata, H., Ohsawa, R., Morishima, H., et al. (2003a). Mapping QTLs influencing rice floral morphology using recombinant inbred lines derived from a cross between Oryza sativa L. and Oryza rufipogon Griff. Theor. Appl. Genet. 107, 218–226. doi: 10.1007/s00122-003-1227-y
Uga, Y., Fukuta, Y., Ohsawa, R., and Fujimura, T. (2003b). Variations of floral traits in Asian cultivated rice (Oryza sativa L.) and its wild relatives (O. rufipogon Griff.). Breed. Sci. 53, 345–352. doi: 10.1270/jsbbs.53.345
Uwatoko, N., Onishi, A., Ikeda, Y., Kontani, M., Sasaki, A., Matsubara, K., et al. (2008). Epistasis among the three major flowering time genes in rice: coordinate changes of photoperiod sensitivity, basic vegetative growth and optimum photoperiod. Euphytica 163, 167–175. doi: 10.1007/s10681-007-9584-2
Virmani, S. S., and Athwal, D. S. (1973). Genetic variability in floral characteristics influencing outcrossing in Oryza sativa L. Crop Sci. 13, 66–67. doi: 10.2135/cropsci1973.0011183X001300010019x
Wang, S., Wu, K., Yuan, Q., Liu, X., Liu, Z., Lin, X., et al. (2012). Control of grain size, shape and quality by OsSPL16 in rice. Nat. Genet. 44, 950–954. doi: 10.1038/ng.2327
Xi, Z., He, F., Zeng, R., Zhang, Z., Ding, X., Li, W., et al. (2006). Development of a wide population of chromosome single-segment substitution lines in the genetic background of an elite cultivar of rice (Oryza sativa L.). Genome 49, 476–484. doi: 10.1139/g06-005
Xiong, L., Liu, K., Dai, X., Xu, C., and Zhang, Q. (1999). Identification of genetic factors controlling domestication-related traits of rice using an F2 population of a cross between Oryza sativa and O. rufipogon. Theor. Appl. Genet. 98, 243–251. doi: 10.1007/s001220051064
Yang, T., Zhang, S., Zhao, J., Liu, Q., Huang, Z., Mao, X., et al. (2016). Identification and pyramiding of QTLs for cold tolerance at the bud bursting and the seedling stages by use of single segment substitution lines in rice (Oryza sativa L.). Mol. Breed. 36, 96. doi: 10.1007/s11032-016-0520-9
Yang, W., Liang, J., Hao, Q., Luan, X., Tan, Q., Lin, S., et al. (2021a). Fine mapping of two grain chalkiness QTLs sensitive to high temperature in rice. Rice 14, 33. doi: 10.1186/s12284-021-00476-x
Yang, W., Xiong, L., Liang, J., Hao, Q., Luan, X., Tan, Q., et al. (2021b). Substitution mapping of two closely linked QTLs on chromosome 8 controlling grain chalkiness in rice. Rice 14, 85. doi: 10.1186/s12284-021-00526-4
Yang, Z., Jin, L., Zhu, H., Wang, S., Zhang, G., and Liu, G. (2018). Analysis of epistasis among QTLs on heading date based on single segment substitution lines in rice. Sci. Rep. 8, 3059. doi: 10.1038/s41598-018-20690-w
Ying, C., and Zhang, Q. (1989) Studies on the character of stigma exsertion among some of Oryza species. Chin. J. Rice Sci. 3, 62–66.
Yuan, L. (2017). Progress in super-hybrid rice breeding. Crop J. 5, 100–102. doi: 10.1016/j.cj.2017.02.001
Yuan, L., and Virmani, S. (1988). “Status of hybrid rice research and development,” in Hybrid Rice. eds Smith, W. H., Bostian, L. R., and Cervantes, E. (Manila: Internationa Rice Research Institute), 7–24.
Zamir, D. (2001). Improving plant breeding with exotic genetic libraries. Nat. Rev. Genet. 2, 983–989. doi: 10.1038/35103590
Zeng, R., Zhang, Z., He, F., Xi, Z., Akshay, T., Shi, J., et al. (2006). Identification of multiple alleles at the Wx locus and development of single segment substitution lines for the alleles in rice. Rice Sci. 13, 9–14. doi: 10.1360/aps040074
Zhan, P., Ma, S., Xiao, Z., Li, F., Wei, X., Lin, S., et al. (2022). Natural variations in grain length 10 (GL10) regulate rice grain size. J. Genet. Genomics. 49, 405–413. doi: 10.1016/j.jgg.2022.01.008
Zhan, P., Wei, X., Xiao, Z., Wang, X., Ma, S., Lin, S., et al. (2021). GW10, a member of P450 subfamily regulates grain size and grain number in rice. Theor. Appl. Genet. 134, 3941–3950. doi: 10.1007/s00122-021-03939-3
Zhang, G. (2021). Target chromosome-segment substitution: a way to breeding by design in rice. Crop J. 9, 658–668. doi: 10.1016/j.cj.2021.03.001
Zhang, G. (2022). The next generation of rice: inter-subspecific indica-japonica hybrid rice. Front. Plant Sci. 13, 857896. doi: 10.3389/fpls.2022.857896
Zhang, G., Zeng, R., Zhang, Z., Ding, X., Li, W., Liu, G., et al. (2004). The construction of a library of single segment substitution lines in rice (Oryza sativa L.). Rice Genet. Newsl. 21, 85–87.
Zhang, S., Huang, X., and Han, B. (2021). Understanding the genetic basis of rice heterosis: advances and prospects. Crop J. 9, 688–692. doi: 10.1016/j.cj.2021.03.01
Zhang, Y., Yang, J., Shan, Z., Chen, S., Qiao, W., Zhu, X., et al. (2012). Substitution mapping of QTLs for blast resistance with SSSLs in rice (Oryza sativa L.). Euphytica 184, 141–150. doi: 10.1007/s10681-011-0601-0
Zhou, H., Li, P., Xie, W., Hussain, S., Li, Y., Xia, D., et al. (2017a). Genome-wide association analyses reveal the genetic basis of stigma exsertion in rice. Mol. Plant 10, 634–644. doi: 10.1016/j.molp.2017.01.001
Zhou, Y., Xie, Y., Cai, J., Liu, C., Zhu, H., Jiang, R., et al. (2017b). Substitution mapping of QTLs controlling seed dormancy using single segment substitution lines derived from multiple cultivated rice donors in seven cropping seasons. Theor. Appl. Genet. 130, 1191–1205. doi: 10.1007/s00122-017-2881-9
Zhu, H., Li, Y., Liang, J., Luan, X., Xu, P., Wang, S., et al. (2018). Analysis of QTLs on heading date based on single segment substitution lines in rice (Oryza sativa L.). Sci. Rep. 8, 13232. doi: 10.1038/s41598-018-31377-7
Zhu, Y., Zuo, S., Chen, Z., Chen, X., Li, G., Zhang, Y., et al. (2014). Identification of two major rice sheath blight resistance QTLs, qSB1-1HJX74 and qSB11HJX74, in field trials using chromosome segment substitution lines. Plant Dis. 98, 1112–1121. doi: 10.1094/PDIS-10-13-1095-RE
Keywords: outcrossing, stigma exsertion, QTL-pyramiding, epistasis, trait reconstruction, rice
Citation: Tan Q, Bu S, Chen G, Yan Z, Chang Z, Zhu H, Yang W, Zhan P, Lin S, Xiong L, Chen S, Liu G, Liu Z, Wang S and Zhang G (2022) Reconstruction of the High Stigma Exsertion Rate Trait in Rice by Pyramiding Multiple QTLs. Front. Plant Sci. 13:921700. doi: 10.3389/fpls.2022.921700
Received: 16 April 2022; Accepted: 05 May 2022;
Published: 07 June 2022.
Edited by:
Dayun Tao, Yunnan Academy of Agricultural Sciences, ChinaReviewed by:
Ryo Ishikawa, Kobe University, JapanJiawu Zhou, Yunnan Academy of Agricultural Sciences, China
Copyright © 2022 Tan, Bu, Chen, Yan, Chang, Zhu, Yang, Zhan, Lin, Xiong, Chen, Liu, Liu, Wang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shaokui Wang, c2hhb2t1aXdhbmcmI3gwMDA0MDtzY2F1LmVkdS5jbg==; Guiquan Zhang, Z3F6aGFuZyYjeDAwMDQwO3NjYXUuZWR1LmNu
†These authors have contributed equally to this work