search for




 

Identification of Heterosis QTLs for Yield and Yield-Related Traits in Indica-Japonica Recombinant Inbred Lines of Rice (Oryza sativa L.)
Plant Breeding and Biotechnology 2017;5:371-389
Published online December 1, 2017
© 2017 Korean Society of Breeding Science.

Chang-Kug Kim1,†, Sang-Ho Chu2,†, Han Yong Park3, Jeonghwan Seo4, Backki Kim4, Gileung Lee4, Hee-Jong Koh4,*, and Joong Hyoun Chin5,*

1Genomics Division, National Institute of Agricultural Sciences, Rural Development Administration, Jeonju 54874, Korea, 2Center for Crop Breeding on Omics and Artificial Intelligence, Kongju National University, Yesan 32439, Korea, 3Department of Bioresources Engineering, College of Life Sciences, Sejong University, Seoul 05006, Korea, 4Department of Plant Science and Research Institute for Agriculture and Life Sciences, Plant Genomics and Breeding Institute, Seoul National University, Seoul 08826, Korea, 5Graduate School of Integrated Bioindustry, Sejong University, Seoul 05006, Korea
Correspondence to: Joong Hyoun Chin, jhchin@sejong.ac.kr, Tel: +82-2-6935-3897
Received November 22, 2017; Revised November 24, 2017; Accepted November 24, 2017.
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract

Supplying sufficient rice to growing populations is a global challenge. Hybrid indica rice varieties exploiting heterosis have increased yields, but inter-subspecific crosses between indica and japonica varieties are hampered by sterility. Examination and genetic understanding of yield heterosis in indica/japonica crosses addressing yield barriers are basic requirements. In this study, QTLs for heterosis of yield traits were identified in indica-japonica recombinant inbred lines (RILs) using a total of 178 RILs originating from Dasanbyeo (indica) × TR22183 (japonica) (DT-RILs) and their backcrossed populations. Nine of sixty-six major quantitative trait loci (QTLs) identified in DT-RILs exhibited heterosis. Heterosis QTLs clustered with other traits on chromosomes 1, 4, and 8, and clusters were conserved between different RILs. The clusters contained several known yield enhancement genes/QTLs. Specific heterotic allele combinations contributed to four major heterosis QTLs, particularly for panicle and spikelet number traits. Heterosis for yield and yield-related traits was explained by the harmonized effects of overdominance, dominance, and epistatic interactions in inter-subspecific breeding populations.

Keywords : Rice, Heterosis, Yield, QTL
INTRODUCTION

Rice is one of the most important staple foods worldwide. Yearly global rice consumption is estimated to reach 500 million tonnes by 2050 (Abdullah et al. 2008). Arable land area per capita has declined substantially in recent decades, from 0.37 ha in 1961 to 0.196 ha in 2014 ( http://data.worldbank.org/indicator/%20AG.LND.ARBL.HA.PC), and supplying sufficient food to a growing population is a major challenge for future human development. Utilization of high yield potential traits will be needed to facilitate the development of high yield rice production systems to overcome the global rice deficit. One of the most effective approaches was used in China, where hybrid rice programs utilizing heterosis contributed to increased yields of up to 13.5 tonnes/ha in 2012 (FAORAP and APSA 2014).

Heterosis, also known as hybrid vigor, is observed in plants, animals, and fungi and is characterized by the superior phenotypic performance of a hybrid compared with its parents. Hybrid vigor is observed in a range of environmental conditions, and phenotypic performance is stably observed. Despite substantial research on heterosis, the genetic underpinnings of the phenomenon are not yet fully understood. Dominance, overdominance, and epistasis (or digenic interactions) are three major genetic components involved in heterosis (Li et al. 2001; Luo et al. 2001; Li et al. 2008; Shen et al. 2014). Dominance and overdominance are interactions of different alleles at the same locus, and epistasis describes interactions between different loci.

Currently available commercial rice hybrids utilizing heterosis optimize the yield potential of crosses between indica varieties. Several recent studies investigated the potential of indica-japonica hybrids (Chin et al. 2011; Xin et al. 2011; Chu et al. 2012; Dan et al. 2014); however, severe hybrid sterility was observed, creating a major barrier to developing super high-yielding hybrid rice (Chu et al. 2012; Dan et al. 2014). The relationship between hybrid sterility and heterosis is unclear, but yield decreases in F1 plants of inter-subspecific crosses were attributed to severely low seed set. A method of yield adjustment, namely, presupposition yield, was suggested as a means of experimentally minimizing the effect of hybrid sterility (Chu et al. 2012).

Broad-sense heritability of yield-related traits was 21.47–79.85% for seven yield-related traits of testcross populations derived from indica-japonica chromosome segment substitution lines (Tao et al. 2016). The contribution of additive effects (narrow sense heritability) to heterosis was particularly large for grain yield and thousand grain weight (TGW). Dominant effects for spikelet and panicle number also had a major impact on heterosis.

In this study, two recombinant inbred lines (RILs) derived from indica (Korean Tongil type) / japonica crosses were used to investigate the genetic contributions of dominance, overdominance, and epistasis to yield and yield-related traits. Major quantitative trait loci (MQTLs) and EpQTLs (epistatic QTLs or digenic interactions) were identified for each RIL. Comparisons of QTLs/genes in two RILs allowed identification of putative heterosis QTL hotspots. MQTL analysis used average values for the three heterotic allele types (AiAi/AiAj, AiAj/AiAi, and AiAj/AiAj) present in a BC1F1 population. The genetic effects of each allele type were then analyzed separately to understand which heterotic allele types were the major contributors to yield and yield-related traits. Lastly, the direct effects of yield-related traits on yield heterosis QTLs were determined.

MATERIALS AND METHODS

Plant materials

A total of 179 F8 RILs derived using the single-seed descendant method from a cross between Dasanbyeo (Tongil type, indica) and TR22183 (temperate japonica) developed by Seoul National University in Korea (DT-RILs) were used in this study (Chin et al. 2011; Jiang et al. 2011). Each RIL was crossed with Dasanbyeo as the female parent, and 178 BC1F1 lines (except for line 40) were produced (BC1F1 line). F1 plants of Dasanbyeo and TR22183 were generated and cultivated alongside the BC1F1 line to allow identification of specific heterosis traits.

Experiment design and phenotyping

All plant materials were grown at Seoul National University Experimental Farm. Field preparation and crop cultivation were as described previously (Chu et al. 2012). Nine traits were measured and scored: culm length (CL, cm), days to heading (DTH, days), panicle number per plant (PN), panicle length (PL, cm), spikelet number per panicle (SN), spikelet fertility (SF, %), TGW (g), grain yield per plant (GYP), and presupposition yield per plant (PYP). Of these, PN, SN, SF, and TGW were regarded as yield-component traits. To remove the effect of hybrid sterility, PYP was estimated and applied where substantial hybrid sterility was observed in BC1F1 lines (Chu et al. 2012). PYP was calculated if the specific BC1F1 line exhibited fertility <90%. As SF had a proportional effect on GYP, an adjustment for grain yield value was used (Supplementary Fig. S1). The equation used was PYP = GYP × 90 / SF. Phenotyping results of RILs and BC1F1s and calculated mid-parent heterosis (MPH) values are described as the ‘RIL’, ‘BCF1’, and ‘MPH’ datasets, respectively.

Linkage map construction and QTL mapping analysis

Genomic DNA was extracted from 2 cm fresh leaves from seedlings at the three-leaf stage using cetyl-trimethylammonium-bromide (CTAB) as previously described (Rogers and Bendich 1988). Across the 12 chromosomes, 162 SSR and subspecies-specific STS markers (Chin et al. 2007) that showed polymorphisms between Dasanbyeo and TR22183 were used to construct a linkage map with JoinMap version 3.0 software (Stam 1993). Only markers with LOD scores >3.0 were used.

MPH values were calculated for each RIL using the following equation, as given in Chu et al. (2012): (value of MPH) = (average value of the corresponding BC1F1) – (average value of RIL and Dasanbyeo). All datasets were analyzed using QTLmapper version 1.6 for the identification of MQTLs and QTLmapper version 2.0 for EpQTLs (Wang et al. 1999). For MQTL analysis, interval mapping using a mixed model approach was applied. For EpQTLs (digenic interactions), only significant two-locus interactions below P < 0.001 were used for further analysis.

Raw genotyping and phenotyping data from another set of indica-japonica RILs, MT-RILs (Milyang 23 × Tong 88-7) (Chu et al. 2012), were re-analyzed using the same software and options to allow comparative analysis with DT-RILs.

Determination and detection of heterosis QTLs

The genetic components of heterosis were suggested to be overdominance, dominance, and epistasis (Li et al. 2008). Different types of dominance (overdominance, partial-to-complete dominance, and underdominance) can be determined by estimating additive effect (a) and dominance effect (d) in RIL and MPH datasets, respectively (Supplementary Fig. S1). Theoretically, sum of dominance and additive effects (a + d) estimated from hetero alleles of linked loci can act as another indicator of heterosis QTLs in the BCF1 dataset. The QTLs detected in both the BCF1 and MPH datasets were regarded as candidates for heterosis QTLs. QTLs detected only in RILs and/or the BC1F1 population were classified as additive QTLs (‘A’ type). Dominant QTLs (‘D’ type) and overdominant QTLs (‘OD’ type) were determined using the values of d/a and a + d. QTLs with 0 < d/a < 1 or (a + d) ≥ 2d were designated as ‘D’ type QTLs. QTLs with d/a > 1, or 2d > a + d, or those detected only in MPH datasets, were designated as ‘OD’ type QTLs (Mei et al. 2005; Chu et al. 2012). QTLs detected in MPH datasets but showing negative ‘d’ values were defined as underdominance QTLs (‘UD’ type). Of the detected QTLs, only ‘D’, ‘OD’, and ‘UD’ type QTLs were used for subsequent comparative studies and direct effect analysis.

The allele types that contributed to heterosis in the three heterotic combinations (homo/hetero, hetero/homo, and hetero/hetero) were identified for MQTLs for GYP and PYP. The additive values of yield and yield-related traits (PN, SN, SF, and TGW) of the three heterotic allele types in BC1F1 populations were calculated and compared. All statistical analyses were conducted using Minitab version 14.

Comparative mapping study and direct effect analysis

Re-analyzed QTL results and the map of MT-RILs reported by Chu et al. (2012) were compared with those of DT-RILs. Primer sequences of linked markers for all the identified heterosis QTLs were searched using GRAMENE Blast tools ( http://ensembl.gramene.org, December 15, 2015) to identify their physical positions in the Nipponbare genome. QTLs anchored in those regions were compared with the heterosis QTLs from this study using the QTARO database ( http://qtaro.abr.affrc.go.jp/).

SPSS-AMOS version 21 was used to estimate the standardized direct and indirect effects of yield-related traits (PN, SN, SF, and TGW). The contribution fractions of MQTLs and EpQTLs for each trait were calculated and analyzed together with the standardized direct effects of the yield-related trait components. For each GYP and PYP QTL in two RILs, the traits with the strongest positive effects were proposed if the standardized direct effect value was >0.5.

RESULTS

Trait performance with respect to yield and yield-related traits

Significant differences were observed between Dasanbyeo (P1) and TR22183 (P2) for four traits (DTH, CL, PN, and SN), with no significant differences seen for the remaining traits (PL, SF, TGW, GYP, and PYP) (Table 1). TR22183 exhibited earlier heading (88 days after transplanting) than Dasanbyeo but maintained the same productivity. Values of CL, PL, SN, TGW, GYP, and PYP were significantly elevated in F1 plants, but very low fertility was observed. In DT-RILs, transgressive segregation was observed in CL, PL, SN, SF, TGW, GYP, and PYP. With the exception of PN and SF, the trait value distributions differed between the BC1F1 and RIL populations (Fig. 1). The median panicle trait values were smaller for BC1F1 than for F1 plants in SN, PL, TGW, and PYP. After fertility recovery by backcrossing, median SF and GYP values were larger in BC1F1 than in F1 plants.

MQTL identification

A total of 66 QTLs (15 QTLs for yield traits and 51 QTLs for the other traits) were identified in the RIL, BCF1, or MPH datasets for the traits used in this study (Table 2). Of these, 19 and 14 QTLs were identified in the BCF1 and MPH datasets, respectively.

Culm length (CL)

Eight QTLs in RIL, four QTLs in BCF1, and one QTL in MPH explained 81.4%, 36.4%, and 21.6% of total variance, respectively. No common QTLs were identified across all datasets, but two QTLs of BCF1 (qCL1.5 and qCL4.3) were adjacent to RIL QTLs (qCL1.4 and qCL4.2). One overdominance (‘OD’ type) QTL, which had the effect of extending CL by 2.11 cm, was detected on chromosome 9 (qCL9.1). All other QTLs were classified as additive (‘A’ type) QTLs.

Days to heading (DTH)

Four QTLs in RIL and one QTL in BCF1 explained 42.3% and 12.4% of total variance, respectively. No QTLs showing the effects of heterotic alleles were identified.

Panicle number per plant (PN)

Three QTLs in RIL and one QTL in BCF1 explained 23.5% and 16.3% of total variance, respectively. No heterosis QTLs were identified for this trait.

Panicle length (PL)

Three QTLs in RIL, two QTLs in BCF1, and two QTLs in MPH were responsible for 59.1%, 16.4%, and 24.8% of total variance, respectively. Two heterosis QTLs were identified in MPH. One dominance (‘D’ type) QTL with relatively large additive (36.8% phenotypic variance explained; PVE) and dominance effects (15.4% PVE) was identified in RIL on chromosome 8 (qPL8.1). Another overdominance QTL (qPL4.1), which had a smaller dominance effect (9.4% PVE) than qPL8.1, was identified in MPH.

Spikelet number per panicle (SN)

Five QTLs in RIL, two QTLs in BCF1, and one QTL in MPH explained 43.9%, 17.9%, and 10.2% of total phenotypic variance, respectively. One overdominance heterosis QTL (qSN9.1), which resulted in a 6.2% increase in SN, was identified on chromosome 9 in MPH.

Spikelet fertility (SF)

Two QTLs in RIL and two QTLs in BCF1 were responsible for 13.7% and 27.6% of total phenotypic variance, respectively. For SF, examination of hybrid sterility in BCF1 is more meaningful than consideration of heterosis QTL. Two QTLs in BCF1 (qSF4.1 and qSF6.1) contributed to 7.86% and 19.69% reductions in SF, respectively.

Thousand grain weight (TGW)

Eight QTLs in RIL, two QTLs in BCF1, and six QTLs in MPH explained 50.4%, 22.4%, and 43.9% of total phenotypic variance, respectively. Three overdominance heterosis QTLs (qTGW3.3, qTGW4.1, and qTGW10.1) and two dominance heterosis QTLs (qTGW1.1 and qTGW7.1) were identified. One underdominance QTL (qTGW4.2) with a negative dominance effect (−0.41) was identified on chromosome 4, adjacent to qSF4.1. Excepting qTGW4.2, the total PVE of the five heterosis QTLs was 41.7%.

Grain yield per plant (GYP)

Three QTLs in RIL, three QTLs in BCF1, and one QTL in MPH were responsible for 30.7%, 34.4%, and 16.0% of total phenotypic variance, respectively. No heterosis QTLs were identified for GYP in DT-RIL. One identified QTL on chromosome 5 (qGYP5.1) showed a negative dominance effect (−2.72%).

Presupposition yield per plant (PYP)

Three QTLs in RIL, two QTLs in BCF1, and three QTLs in MPH explained 27.3%, 15.1%, and 23.4% of total phenotypic variance. Two QTLs (qPYP8.1 and qPYP9.1) were also QTLs for GYP (qGYP8.1 and qGYP9.1). One QTL for GYP on chromosome 11 (qGYP11.1), which showed a negative additive effect (−0.9), was not detected in the PYP dataset, as expected. One QTL on chromosome 2 (qPYP2.1) was not found for GYP. Three overdominance heterosis QTLs were identified on chromosomes 1, 7, and 8. The removal of hybrid sterility effects allowed possible QTLs for GYP to be identified.

The QTLs identified only in the BCF1 population were not classified as heterosis QTLs in this study. However, many of these QTLs showed positive values of combined additive and dominant effects: all four QTLs for CL; one QTL for DTH; one QTL for PN; one of two QTLs for PL (qPL9.1); one of two QTLs for SN (qSN4.2); one of two QTLs for TGW (qTGW1.2), and both QTLs for PYP.

Many of the QTLs were clustered within specific regions on chromosome 1 (RM297-S01143A), chromosome 4 (S04120-S04129B), and chromosome 8 (RM25-RM72) (Fig. 2, left chromosome shown).

EpQTL identification

When all traits were considered, 23 EpQTLs in RIL, 15 EpQTLs in BCF1, and two EpQTLs in MPH were identified in this study (Table 3). The total PVE of EpQTLs for each trait was 10.7–27.5% for RIL, 11.2–26.8% for BCF1, and 7.9–22.0% for MPH. No EpQTLs were identified for SN and PYP in RIL or for PN and TGW in BCF1.

For grain yield traits, different EpQTLs for GYP were identified in RIL and BCF1. Only one EpQTL, between chromosome 1 (S01143A-S01157B, qSF1.1 of RIL) and chromosome 8 (RM72-S08052B, qGYP8.1 of RIL), contributed positively to higher grain yield (1.75 g/plant). For PYP, one EpQTL between chromosome 3 (S03010B-S03020) and chromosome 10 (S10071-RM590, qSN10.1 of RIL) had positive additive effects.

One EpQTL for PL and one for SN were identified in the MPH dataset. No MQTLs were identified within the chromosomal regions linked with EpQTLs for PL. However, for SN, qPL12.1 in RIL interacted with the region on chromosome 4.

Trait contribution to heterosis QTLs for grain yield

First, to identify the yield-component traits contributing substantially to the three overdominance heterosis QTLs and one underdominance QTL for grain yield in DT-RIL, the additive effects of each yield-component trait on the QTLs were presented in radar charts (Fig. 3). Interval mapping through backcrossing was used to identify QTLs between markers. There were three types of heterotic allele combination for each marker (AiAi/AiAj, AiAj/AiAi, or AiAj/AiAj), and the additive effects of each allele type were calculated by dividing the value of the corresponding allele type with that of the AiAi/AiAi (homo) allele genotype. Additive effects identified in this study were in the −15–15% range (Table 4).

For qPYP1.1, the AiAi/AiAj (S01140/S01143A) genotypes did not contribute to heterosis. However, AiAj/AiAi genotypes had a significantly larger PN, and significantly or highly significantly smaller SN, SF, and GYP, than homoallele genotypes. Likewise, AiAj/AiAj genotypes had significantly higher PN and PYP values than homoallele genotypes. This suggested that the PN of AiAj/AiAj genotypes might be the major factor explaining the significant heterotic effect on PYP.

For qPYP7.1, AiAi/AiAj (RM234/S07101) genotypes showed negative additive effects on PYP. The AiAj/AiAi and AiAj/AiAj genotypes had significantly higher PYP values than homoallele genotypes. The additive effect of AiAj/AiAi on SN might explain the significantly positive additive effect of PYP. Although a significant positive additive effect of SF was found for AiAi/AiAj, the PYP values of the genotypes were negative, thus demonstrating that SF alone could not be a factor of heterosis.

For qPYP8.2, AiAi/AiAj (S08090/S08106) and AiAj/AiAj showed significant additive effects for PYP, which may have been caused by low hybrid fertility. AiAj/AiAj genotypes showed a positive effect of SF, but PYP did not show additive effects.

Unlike the PYP cases, no heterosis QTLs were identified in GYP. However, one underdominance GYP QTL (qGYP5.1) was analyzed to see which allele types were contributory. Only AiAj/AiAj (S05014B/RM267) showed a significant negative effect on SF and GYP. Two other heterotic allele genotypes of linked markers (AiAj/AiAi and AiAi/AiAj) contributed to significantly higher SF and GYP values. This result indicates the importance of studying different heterotic allele combinations of linked markers for underdominance QTLs.

Second, the direct effects of each trait on grain yield were analyzed using the path analysis module of SPSS-AMOS software (Table 5). Only DTH, PN, SN, SF, and TGW were considered as yield-related traits for this analysis, because CL was a vegetative trait and PL was highly correlated with SN. In each heterosis QTL for grain yield, the standardized direct effects of each trait >0.5 or <−0.5 were considered as positive effect QTLs. As expected from the radial chart analysis, PN and SN were the most important traits explaining heterosis of qPYP1.1. Likewise, PN, SN, and TGW variously contributed to qPYP7.1 for each allele combination. qPYP8.2 was related to PN and SN for the AiAj/AiAi type only. PN, SN, and SF showed direct effects on GYP via the contribution of the AiAj/AiAi and AiAi/AiAj genotypes at qGYP5.1. SF was a negatively effective trait in all heterosis QTLs for PYP.

DISCUSSION

Heterosis QTLs of intra- and inter-locus interactions

Heterosis is governed by various genetic effects including overdominance effects, dominance effects, additive effects, and their interactions. In this study, indica-japonica RILs were used to identify heterosis QTLs for grain yield and its related traits. Theoretically, heterotic effects (inter-locus heterosis) should be explained solely by the genetic interactions between different recombinant blocks originating from different parents. However, two-QTL interactions (or digenic interactions) are difficult for plant breeders to utilize due to their complexity and the difficulties in dealing with the large numbers of samples in breeding programs. Accessing heterosis through intralocus interactions between different alleles (dominance and overdominance) may thus be more amenable. Intra-and inter-locus heterosis can be assessed using BC1F1 populations, which were used to study overdominance and dominance effects in this study. Inter-locus heterogeneity is more prevalent in RIL rather than BC1F1 populations, and estimation of intra-locus heterosis may be improved using BCF1 and MPH datasets. The heterotic effects on inter-locus and intra-locus interactions governing grain yield may not be identified at the same loci in different populations. For example, two EpQTLs for GYP in RIL were unidentified in the BCF1 and MPH datasets (Tables 2 and 3). On the other hand, two EpQTLs for GYP and two EpQTLs for PYP in MT-RILs (RILs derived from Milyang 23 and Tong 88-7; Chu et al. 2012) were not identified in heterosis MQTL analysis using BCF1 and MPH datasets in DT-RILs (Supplementary Tables S1 and S2).

In DT-RIL, there were nine heterosis MQTLs for all traits and one underdominance heterosis MQTL for GYP (Table 2). EpQTLs for heterosis are rarely detected if the population size is small, because the effects of multi-locus allele combinations require a large number of samples and meaningful statistical analysis for detection. Thus, the subsequent heterosis QTL analyses in this study were conducted using only MQTL results. Nevertheless, there were clearly no EpQTLs for GYP or PYP found using MPH datasets in DT-RIL. Conversely, in MPH datasets of MT-RIL, three overdominance MQTLs for GYP and two for PYP were identified (Supplementary Table S1).

Common heterosis QTL clusters in two indica-japonica RILs indicating relatedness of heterosis for yield and for yield-related traits

MQTLs for different traits were co-located at some chromosomal regions (Fig. 2), and several QTL clusters were observed in MT-RILs, on chromosome 1 (S01143A-S01157B), chromosome 4 (RM348-S04128), chromosome 6 (S06015-S06019B, RM50-RM527, and S06087-S06100), and chromosome 7 (S07035-S07055B). In both RILs, QTL clusters included QTLs identified in RIL, BCF1, and MPH datasets for several yield-related traits. Thus, an additive effect of preferred alleles in RILs would be expected in marker-assisted breeding programs that utilized the heterosis QTLs in the QTL clusters. This may explain the plant growth vigor seen in recombinant plants generated through introgression of donor QTLs in some generations of breeding programs.

Heterosis QTLs for each trait were compared between the two RILs (Supplementary Table S3), and no common heterosis QTLs were found. However, the QTL clusters may represent common loci for heterosis (Fig. 2). QTL clusters in DT-RIL (RM297-S01143A) and MT-RIL (S01143A-S01157B) were at adjacent regions on chromosome 1, and might share common genetic elements for heterosis. This was also observed for S04120-S04129B (CL, SF, and TGW) in DT-RIL and RM348-S04128 (PL, GYP, and SN) in MT-RIL on chromosome 4. Three QTL regions linked to grain yield traits are close in comparative linkage maps and should be studied further to see whether common regions for heterosis are shared between the RILs (Fig. 2): RM584-RM204 in DT-RIL and S06015-S06018 in MT-RIL on chromosome 6, RM253-S06031 (SF in DT-RIL) and RM50-RM527 (PYP in MT-RIL) on chromosome 6, and RM21-S11117 (GYP in DT-RIL) and S11112-S11117B (PYP in MT-RIL) on chromosome 11.

Previously reported QTLs that are linked to the identified heterosis QTLs in DT and MT-RILs and are linked to a range of agronomic traits are listed in Supplementary Table S4. QTL locations were compared between the physical map locations of the markers in this study and the linked markers from various mapping populations. In DT- and MT-RILs, qTGW4.2, qGYP8.1, qSN8.1, and qCL9.1/qSN9.1 were not previously reported as yield or yield-related QTLs. With the exception of these, all the QTL regions in this study were previously linked to yield or yield-related traits. For example, qSN1.1 of MT-RIL was co-located with gp1a in F2 (grains per panicle, Shen et al. 2003), qPN1 in BC4F4 (panicles per plant, Tian et al. 2006), and others. Thus, some previously identified QTLs linked to different traits may be co-located with or the same as QTLs identified in this study. This may explain the significant positive correlations observed between some yield traits. For example, with respect to the linked QTLs given above, grain number (gp1a) is correlated with spikelet number (qSN1.1) and panicle number (qPN1.1).

Previous studies of RILs rarely discussed whether these genes/QTLs exhibited heterosis. Some previously reported QTLs/genes were co-located with the heterosis QTLs identified in this study, and further analysis is needed to determine whether heterosis is seen in heterotic alleles. For DT-RIL, these were as follows: qTGW3.3 (gw3a, grain weight, Yu et al. 1997); qGYP5.1 (S31, hybrid sterility, Li et al. 2005; qSF5.2, hybrid sterility, Reflinur et al. 2012); qPYP7.1 (qPLLN7, pollen fertility, Yan et al. 2003; hd7a, heading date, Yu et al. 2002); qTGW7.1 (S23, hybrid sterility, Sobrizal et al. 2000); qPL8.1 (RH8, heading date, Li et al. 2016; qPL8, Hd5, heading date, Lin et al. 2003; PL, Miyata et al. 2007); and qPYP8.2 (WFP, panicle branching and yield, Miura et al. 2010). For MT-RIL, these were qCL1.2 (sd1, plant height, Li et al. 2003); qGYP6.1 (DPL2, sterility, Mizuta et al. 2010; qSPTF6, SF, Yan et al. 2003); qPYP6.1 (spf6, SF, Song et al. 2005; S5, SF, Yan et al. 2000); qCL7.2 (Hd4, DTH, Lin et al. 2003); and qPYP11.1 (FON3, floral organ number, Jiang et al. 2007). These QTLs/genes were reported in various population types as well as in backcrossed populations, implying that QTLs showing additive effects in homozygous alleles may give heterotic effects.

Comparative analysis indicated that the negative additive values of some QTLs of DT- and MT-RILs could be partially explained by hybrid sterility QTLs/genes. qGYP5.1, qGYP6.1, and qGYP10.1 were co-located with known sterility QTLs (qSF5.2 and S31) and QTLs for pollen fertility (Doi et al. 1998; Tan et al. 1998). However, there were no known hybrid sterility QTLs/genes for qTGW4.2, qSN8.1, qTGW9.1, and qGYP10.1. Conversely, qPYP6.1 in MT-RIL did not show negative additive values, indicating a lack of allele-induced hybrid sterility between Milyang 23 and Tong 88-7.

Tao et al. (2016) also reported some heterosis QTLs from testcrossed CSSLs (Chromosome Segment Substitution Lines). The QTL hotspots on chromosomes 1, 6, and 8 identified in this study were also identified by Tao et al. (2016).

Three heterotic allele genotypes on heterosis QTLs of MT-RILs

For QTL determination, the effects of three heterotic allele combinations (AiAi/AiAj, AiAj/AiAi, and AiAj/AjAj) were combined, hampering the identification of the allele types that had the most substantial impact on phenotypic values (Fig. 3 and Supplementary Fig. S2).

The range of additive values was larger for MT-RILs than for DT-RILs, indicating that the effect of MQTLs was larger for MT-RILs than for DT-RILs. Five QTLs for GYP and two QTLs for PYP were investigated in this study (Supplementary Fig. S2 and Supplementary Table S5). For qGYP1.2 (RM259/RM600), the heterotic allele of RM259 contributed to the additive effect of SF and heterosis of GYP. The heterotic alleles linked to both the markers contributed to the increase of PN. A QTL for panicle number (ppp1.1) was reported on qGYP1.2 (Cho et al. 2007). This was not significant (P = 0.051), but a significant positive effect was seen for heterotic allele type for SF contribution to qGYP5.1. On qPYP6.1 (S06015/S06018), all three heterotic allele types exhibited similar contributions to heterosis of PYP. SN was the most important contributory trait for PYP increase, though SF was significantly low for all three types. A QTL for SN was identified on qPYP6.1 (Zhuang et al. 2001). At qGYP8.1 (S08059/S08066), heterotic allele effects of S08059 and S08066 gave contrasting effects for SF and GYP. S08059 heterotic allele type positively impacted the additive effect of GYP. On qPYP11.1 (S11112/S11117B), the AiAi/AiAj and AiAj/AiAj types showed significantly low SF but high PYP. The presupposition yield calculation was influenced by S11117B heterotic allele type. FON3 (Jiang et al. 2007) reported a QTL for floral organ numbers in the F2 population on the same region. However, it was unclear which trait contributed to the PYP increase of qPYP11.1 in this population. Low SF on qGYP6.1 (RM50/RM527) was impacted by RM527 heterotic allele type, giving significant positive additive effects of PYP and significant negative effects of GYP. Likewise, qGYP10.1 (S10019/S10026C) also showed an underdominance heterosis type.

Direct effect analysis for MQTLs and EpQTLs

It is desirable to identify the yield-related traits or yield-component traits that are most influential on yield. In this study, to simplify the complicated mediation networks, only the standardized direct effect of each yield-related trait on PYP and GYP was estimated by path analysis (Table 5).

The direct effect of each trait was estimated for each different allele type. One homoallele type (AiAi/AiAi) in RIL and three heterotic allele types (AiAi/AiAj, AiAj/AiAi, and AiAj/AiAj) in BCF1 populations were analyzed separately. PN and SN were the yield-related traits with most direct impact on heterosis QTLs in DT-RILs. For qPYP1.1, the strongest trait, with an estimated direct effect >0.5, was SN, particularly for AiAj/AiAi and AiAj/AiAj types (62 of 179 RILs; 34.6% (Table 4). However, for the same locus, PN was the most effective trait for qPYP1.1 with respect to the homoallele (AiAi/AiAi) and complete heterotic allele (AiAj/AiAj) types. The same approach could be applied to other heterosis yield QTLs. Different combinations of yield-related traits provided direct effects for each qPYP7.1 allele type. Only the AiAj/AiAi allele type of qPYP8.2 was linked with PN and SN in BCF1. Conversely, for qGYP5.1, only the AiAi/AiAj type was directly impacted by PN, without the negative effect of SF.

A similar approach was used for MT-RILs (Supplementary Table S6). Unlike with DT-RILs, the positive direct effect on SF was apparent across many QTLs. A direct negative effect of SF was found only for qPYP6.1 in all allele types. qGYP1.2, qGYP5.1, and qGYP10.1 all positively impacted SF. Conversely, SN had the most positive effect in AiAj/AiAj of qGYP1.2, AiAj/AiAj of qPYP6.1, AiAj/AiAi of qGYP6.1, and AiAj/AiAj of qGYP10.1. TGW was influenced by AiAi/AiAj of qPYP6.1 and AiAj/AiAi of qGYP6.1. PN and DTH directly affected qPYP6.1 in partial heterotic allele types (AiAi/AiAj or AiAj/AiAi).

In brief, the standardized direct effects of different heterotic allele types on yield-related traits differed between two indica-japonica RILs. Heterosis effects were identified using the accumulated additive effects of PN, SN, SF, and TGW (Dan et al. 2015). For each QTL, the yield-related trait with the most positive impact could be identified by this analysis. Followed by the development of near-isogenic lines (NILs) for each QTL, the most important yield-related traits could next be studied to estimate the consequent yield heterosis.

The contribution of MQTLs and EpQTLs to heterosis can be estimated by the sum of total PVEs for the identified QTLs. Fig. 4 shows the total effects of PVE of each yield-related trait (center), PYP (outer sides), and GYP (inner sides) connected by arrows representing the direct effects. The total PVEs of yield-related traits were larger in DT-RILs than in the BCF1 population. However, total PVEs of GYP and PYP were larger in the BCF1 population than in RILs. In the DT-BCF1 population, contributions of EpQTLs for GYP (61.2%) and PYP (36.4%) were large, and consequently, the heritability (total PVE) of yield heterosis was larger than GYP (45.8%) and PYP (27.3%) in RILs. This differed for MT-RILs (Supplementary Fig. S3), where heritability of yield heterosis was smaller than in RILs. This implies that the contribution of the dominance effect was larger than that of epistasis between Dasanbyeo (indica) and TR22183 (japonica).

Suggestions for plant breeding strategies utilizing heterosis QTLs

Heterosis QTLs identified by current software provide average values for different allele types of heterotic alleles. The additive values of different heterotic allele combinations could be calculated from the identified QTLs to identify the allele type that contributed to the preferable yield heterosis. Hybrid sterility, which is common in indica-japonica crosses in rice, hampers the interpretation of heterotic effects of yield traits in specific heterotic allele types at those loci. The breakdown of allele types may help resolve this.

The heterosis QTLs identified in the two RILs in this study mostly coincided with previously reported yield and yield-related traits. This similarity was also noted when specifically expressed genes of F1 plants were compared to the linkage map of Nipponbare/93-11 (Venu et al. 2014). This also showed that heterosis could be explained by certain specific yield-related traits associated with specific loci. These loci might be detected as MQTLs in conventional QTL studies.

The identified EpQTLs explained the presence of epistatic interactions between loci for heterosis. To utilize EpQTLs for plant breeding, marker-assisted recurrent selection or genomic selection strategies would be more effective rather than the simple marker-assisted backcrossing strategies used for MQTLs.

Common heterosis QTLs identified across different populations and studies would be highly valuable for heterosis breeding programs; however, further quantitation of the complex genetic networks for each yield-related trait is needed to design breeding programs with high predictability.

ACKNOWLEDGEMENTS

This study was supported by a grant from the Next-Generation BioGreen 21 Program (no. PJ01102401) of the Rural Development Administration, Korea. This research was also supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (no. 2017R1D1A1B04034862).

Supplementary Information
Figures
Fig. 1. Phenotypic distribution of yield and yield-related traits of Dasanbyeo (P1), TR22183 (P2), their respective F1 values, and their recombinant inbred lines (RILs). (a) Culm length, (b) days to heading, (c) panicle number, (d) panicle length, (e) spikelet number per panicle, (f) spikelet fertility, (g) thousand grain weight, (h) grain yield per plant, and (i) presupposition yield per plant.
Fig. 2. Comparative linkage maps of DT-RILs (left, DT) and MT-RILs (right, MT) for the 12 rice chromosomes. Vertical gray bars represent chromosomes. Black dots on each chromosome represent the estimated position of the centromere. QTLs for each trait identified in this study are located next to the linked markers on each chromosome. The legend box below the chromosomes indicates the symbols corresponding to each trait. Vertical dotted lines divide the QTL groups identified in different populations (RILs, BCF1, and MPH). QTLs with black symbols were identified in MPH, gray in BCF1, and blank in RILs. Upward or downward arrows next to QTLs designate the QTLs as overdominant or underdominant, respectively. Dotted lines connect the same markers used in the two RILs.
Fig. 3. Phenotypic additive performance (%) of three different heterotic allele types on yield heterosis QTLs identified in DT-RILs. The radial charts for each yield heterosis QTL show two yield (GYP and PYP) and four yield-related traits (PN, SN, SF, and TGW). Blue (AiAi/AiAj), orange (AiAj/AiAi), and gray (AiAj/AiAj) lines indicate the additive value in percentage over the value AiAi/AiAi (Dasanbyeo) type (dotted closed hexagonal line). Asterisks represent the statistical significance level (*** for P < 0.001; ** for 0.001 < P < 0.01; and * for 0.01 < P < 0.05). The upward/downward arrows indicate the significant additive values of the corresponding heterotic allele types.
Fig. 4. Contribution of MQTL and EpQTL for each trait to yield in DT-RILs. Pie chart size for each trait represents the sum of total PVE of MQTL (dark gray) and EpQTL (light gray) identified in this study. Five yield-related traits were connected by arrows with numbers (standardized direct effect) to GYP (inner side) and PYP (outer side). The direct effects of each trait were shown in two populations: RIL (left) and BCF1 (right).
Tables

Phenotypic performance and heterosis of traits in parents: F1s, RILs, and BCF1 populations of DT-RIL.

DTHz)CLPLPNSNSFTGWGYPPYP
P1 (Dasanbyeo)113.0 ± 0.7 a58.2 ± 2.6 a24.9 ± 1.1 a13.4 ± 1.8 a156.6 ± 10.4 a82.0 ± 7.5 a21.8 ± 0.7 a27.6 ± 2.0 a30.3 a
P2 (TR22183)88.0 ± 1.0 b71.4 ± 4.1 b25.4 ± 1.5 a6.8 ± 1.1 b179.3 ± 16.2 b89.0 ± 7.6 a22.2 ± 0.7 a28.7 ± 4.8 a29.0 a
MP = (P1 + P2) / 2100.5 a64.8 a25.2 a10.1 a168.0 a85.5 a22.0 a28.2 a29.7 a
F1100.0 ± 1.0 a81.7 ± 3.3 b27.4 ± 0.9 b8.9 ± 2.6 b236.3 ± 12.6 b47.3 ± 9.0 b30.5 ± 0.2 b33.0 ± 2.3 b54.1 b
MPHy) = F1 – MP−0.516.92.3−1.265.9−28.06.65.582.5
RILs101.3 ± 8.3 (88.0–130.0)69.1 ± 12.0 (44.2–113.8)24.2 ± 2.8 (18.2–32.3)10.3 ± 2.4 (6.0–16.6)175.3 ± 41.8 (68.0–313.8)80.3 ± 11.1 (37.4–96.8)24.2 ± 3.1 (16.8–31.8)27.9 ± 6.8 (15.1–48)31.8 ± 6.8 (19.2–56.1)
BCF1s (P1 × RILs)106.6 ± 6.1 (97.0–127.0)80.0 ± 7.2 (65.0–107.9)25.9 ± 1.6 (22.4–30.3)9.8 ± 1.6 (7.0–15.4)201.0 ± 29.3 (94.0–266.8)81.9 ± 14.8 (24.8–96.0)27.7 ± 1.6 (23.5–31.8)37.6 ± 6.2 (16.5–52.3)41.9 ± 6.0 (30.9–70.0)
BC MPHy) = BCF1 – (RIL + P1) / 2−0.2 ± 6.0 (−18.5–28.5)16.4 ± 5.2 (−3.3–38.5)1.3 ± 1.5 (−1.7–5.7)−2.0 ± 1.3 (−4.7–2.6)35.0 ± 24.9 (−53.6–91.3)−0.1 ± 15.1 (−84.7–24.0)4.7 ± 1.4 (0.2–9.1)9.5 ± 6.8 (−26.8–26.0)11.1 ± 5.7 (−1.2–34.9)

z)DTH (days to heading, days), CL (culm length, cm), PL (panicle length, cm), PN (panicle number per plant), SN (spikelet number per panicle), SF (fertility, %) TGW (thousand grain weight), GYP (grain yield per plant), PYP (presupposition yield per plant); See the text for the detailed information.

y)MPH: mid-parent heterosis; see the main text for the detailed information.


Major QTLs identified in DT-RIL, BCF1, and MPH datasets.

TraitNameChr.z)IntervalRILBCF1MPHQTL typeu)|d / a| |2d > a + d|



LODy)ax)R2 (%)PVEw) (%)LODa + dR2 (%)PVE (%)LODdv)R2 (%)PVE (%)
CLqCL1.11RM1-RM12981.43.621.926.5836.421.6A
qCL1.31RM297-S011408.424.2411.5A
qCL1.41S01143A-S01157B3.65−3.106.1A
qCL1.51S01157B-S011605.132.047.46A
qCL3.13S03115-S031206.68−3.568.1A
qCL3.23S03145-RM1436.58−3.156.4A
qCL4.14S04023-S040487.55−3.658.5A
qCL4.24S04120-S041287.93−3.357.2A
qCL4.34S04128-RM12710.12.7313.34A
qCL8.18RM25-RM7227.166.8429.9A
qCL9.19S09000A-S090065.032.228.814.712.1115.1OD1.902
qCL11.111RM21-S111176.68−3.688.7A
DTHqDTH3.13S03145-RM1434.42.015.642.312.4No analysisA
qDTH5.15RM334-RM313.091.754.2A
qDTH8.18RM25-RM7217.34.5528.5A
qDTH11.111S11045-S11055A2.66−1.714.0A
qDTH11.211RM229-RM214.792.2112.37A
PNqPN1.11RM297-S0114023.56.280.6016.2516.3A
qPN1.21S01140-S01143A2.24−0.514.4A
qPN9.19S09073-S09075A6.070.8010.6A
qPN11.111S11004A-S11006A4.23−0.718.5A
PLqPL1.11RM297-S0114059.12.35−0.436.8716.424.8A
qPL4.14S04077A-RM4515.380.459.4OD
qPL8.18RM25-RM7219.861.7236.87.830.5815.4D0.337
qPL9.19S09049-S090583.160.509.48A
qPL11.111RM332-S110287.55−0.878.5A
qPL12.112S12091-RM27012.561.0613.8A
SNqSN1.11RM297-S0114043.92.8−7.396.2417.910.2A
qSN1.21S01140-S01143A7.54−15.5113.1A
qSN1.31S01167-S011814.419.845.3A
qSN4.14S04077A-RM4514.32−10.245.7A
qSN4.24RM451-S04097B6.2810.1011.66A
qSN8.18RM25-RM729.2316.4114.7A
qSN9.19S09000A-S090062.566.386.2OD
qSN10.110S10071-RM5906.6713.139.4A
SFqSF1.11S01143A-S01157B4.392.324.813.727.6No analysisA
qSF4.14S04128-RM1274.98−4.177.86A
qSF6.16RM253-S060319.14−6.6019.69A
qSF10.110S10058A-S100717.253.168.9A
TGWqTGW1.11RM297-S0114013.60.928.850.422.46.120.5313.443.9D0.576
qTGW1.21S01167-S011814.72−0.563.35.460.448.39A
qTGW3.13S03002-S03010B5.420.593.6A
qTGW3.23S03046-S0304811.11−1.0010.5A
qTGW3.33S03099-S031152.80.4710.4OD
qTGW4.14S04077A-RM4513.020.397.2OD
qTGW4.24RM127-S04129B3.440.417.9UD
qTGW5.15S05004A-S050094.66−0.593.7A
qTGW5.25S05080-S0510113.630.9910.37.89−0.5713.96A
qTGW7.17RM248-RM1727.130.694.92.410.304.3D0.436
qTGW10.110RM258-S10058A3.020.366.4OD
qTGW12.112RM4A-S12011B11.140.918.7A
GYPqGYP5.15S05014B-RM26730.734.47.742.7216.016UD
qGYP6.16RM584-RM2042.56−1.636.67A
qGYP8.18RM25-RM7220.83.1822.4A
qGYP9.19S09073-S09075A8.211.716.5A
qGYP10.110S10019-S10026C6.71−2.8820.67A
qGYP11.111S11078-RM2292.17−0.901.8A
qGYP11.211RM21-S111172.29−1.687.06A
PYPqPYP1.11S01140-S01143A27.315.13.061.506.523.4OD
qPYP2.12S02026-RM712.921.474.5A
qPYP7.17RM234-S071013.971.698.3OD
qPYP8.18RM25-RM729.442.8416.7A
qPYP8.28S08090-S081062.971.718.5OD
qPYP8.38S08106-S081074.181.697.96A
qPYP9.19S09075A-S09093A5.271.726.1A
qPYP11.111S11055A-S110642.481.607.14A

z)Chr: chromosome.

y)LOD: logarithm of odds numbers.

x)a: additive effect.

w)PVE: phenotypic variance explained.

v)d: dominant effect.

U)A: additive QTL, D: dominant QTL, OD: overdominant QTL, UD: underdominant QTL. See the main text for the details.


EpQTLs of DT-RIL and BCF1 populations and MPH datasets.

Traitz)Chr.Interval iChr.Interval jLODAiy)AjAAijH^2(Ai)H^2(Aj)H^2(AAij)PVEx) (%)
RILCL3S03002-S03010B9S09000A-S090069.921.33**2.480.0070.02615.8
4S04048-S040584S04097B-RM3484.792.340.023
5S05036-RM28912RM235-RM177.122.960.037
7RM336-S0707611S11028-S1104510.961.9***2.140.0150.019
7RM248-RM1729S09065-S0907311.511.03*−3.570.0040.054
DTH8S08107-S08012010S10071-RM5906.422.950.10710.7
PN6S06035-RM59638S08075-S080907.21−0.780.09617.1
8S08052B-S080668S08090-S081064.8−0.720.075
PL1S01167-S011819S09065-S090735.610.670.03316.8
3S03099-S0311510S10058A-S100714.660.620.028
5S05009-S05014B6RM469-RM2255.080.630.029
6S06074-RM307RM336-S070768.630.78***−0.760.0450.043
7RM346-S070848S08107-S0801206.020.87*−0.690.0560.035
SF1RM488-RM2462S02085-RM3185.213.380.07627.5
6RM469-RM2259S09075A-S09093A9.754.580.140
8RM25-RM728S08090-S081067.371.91**2.990.0240.059
TGW2RM71-RM4382S02085-RM3187.07−0.5***−0.560.0210.02325.8
3S03120-S031308S08090-S0810611.51.020.085
6RM3183-S060748S08075-S0809010.250.890.065
8RM25-RM7210S10026C-RM25813.430.37**0.99***0.610.0110.0800.031
11S11028-S1104512S12091-RM2707.42−0.810.054
GYP5S05054-S050648RM25-RM7216.632.96***1.520.1550.04115.1
6RM527-RM318311RM229-RM218.152.490.110
BCF1CL1S01157B-S0116012RM277-S120918.45−1.93***−1.880.0480.04611.2
5S05009-S05014B9S09000A-S090068.27−2.47***2.270.0790.067
DTH2RM475-S0208511S11018-RM3324.77−1*2.010.0170.06813.7
6S06053-RM52711RM229-RM216.98−1.99***2.040.0660.070
PL1RM128-RM2978RM25-RM725.460.44***0.530.0550.08022.8
3S03099-S031154RM348-S041206.53−0.720.148
SN1S01143A-S01157B8S08052B-S080665.4310.930.14814.8
SF1RM128-RM29711RM21-S111174.623.32**3.04*−5.620.0210.0170.05914.7
1RM297-S0114011RM21-S111175.924.47***3.58**−6.010.0370.0240.068
4S04128-RM12711S11045-S11055A5.954.05***−3.250.0310.020
GYP1S01143A-S01157B8RM72-S08052B4.75−0.84*1.750.0130.05626.8
2RM475-S020853S03120-S031305.461.22**−2.220.0270.089
5S05036-RM2896RM469-RM2254.77−2.60.123
PYP3S03010B-S0302010S10071-RM5904.7920.11421.3
4RM348-S041209S09062B-S090654.87−1.860.099
MPHPL1RM128-RM2977S07024-S070384.80.390.0797.9
SN4S04048-S0405812S12091-RM2705.6212.940.22022.0

z)See Table 1 for the trait abbreviations.

y)Ai and Aj are the main effects of locus i and locus j. AAij is the epistatic effect between loci i and j, as defined by Mei et al. (2005). Significance level at

*P ≤ 0.05

**P ≤ 0.01

***P ≤ 0.001.

x)Percentage of the total variation explained by AAij. Significant at P ≤ 0.001.


Percent additive effect of each yield-related trait to its parent haplotype (AiAi) and yield traits in BCF1 of DT-RIL.

QTLMarkerHaplotypenPNSNSFTGWGYPPYP






Az) (%)Py) (T < t)A (%)P (T < t)A (%)P (T < t)A (%)P (T < t)A (%)P (T < t)A (%)P (T < t)
qPYP1.1S01140AiAiS01143AAiAj9−2.340.170−1.240.358−1.200.394−0.890.368−1.520.391−1.320.463
AiAjAiAi1212.130.021−11.750.046−8.750.0480.920.180−8.440.0040.410.459
AiAjAiAj509.690.001−2.130.188−8.870.0080.780.220−0.630.4308.390.002
qGYP5.1S05014BAiAiRM267AiAj224.240.170−5.950.0226.650.0001.490.1456.270.043−0.970.386
AiAjAiAi95.080.1932.970.2775.760.032−0.040.4939.320.0053.820.207
AiAjAiAj462.590.195−0.200.471−11.180.0042.160.013−9.550.0042.690.177
qPYP7.1RM234AiAiS07101AiAj115.040.197−4.870.0926.170.013−1.010.3523.510.180−2.900.180
AiAjAiAi104.620.1998.820.042−6.380.1652.120.2332.590.33810.100.029
AiAjAiAj720.160.4763.240.098−4.570.069−0.370.3240.650.4136.170.006
qPYP8.2S08090AiAiS08106AiAj205.290.0561.560.349−7.720.0680.310.405−0.170.48311.920.011
AiAjAiAi352.950.202−1.990.2701.910.1630.180.4401.530.281−0.150.475
AiAjAiAj574.490.0780.030.496−9.380.004−1.120.147−3.650.1457.700.002

z)Percent additive effect of the corresponding haplotype to its parent haplotype (AiAi).

y)Probability of significance assuming H0 = 0 by Student’s t-test.


Standardized direct effect of yield-related traits to yield in DT-RILs.

QTLAllele typeDTHz)PNSNSFTGWMost positively effective trait (>0.5 or <−0.5)
qPYP1.1AiAi/AiAi0.0830.5120.413−0.4360.108PN
AiAi/AiAj−0.2870.2410.242−0.6860.149(SF)a
AiAj/AiAi−0.3450.0480.548−0.4990.124SN
AiAj/AiAj0.1140.5230.505−0.5420.191PN and SN
qPYP7.1AiAi/AiAi−0.0590.5060.681−0.2750.14PN and SN
AiAi/AiAj0.3241.6331.0510.2580.538PN, SN and TGW
AiAj/AiAi−0.0460.3570.631−0.2720.047SN
AiAj/AiAj−0.0050.5510.281−0.680.191PN, (SF)
qPYP8.2AiAi/AiAi−0.2250.7190.574−0.0580.212PN and SN
AiAi/AiAj−0.3330.2190.246−0.8280.131(SF)
AiAj/AiAi0.1260.8520.831−0.0680.146PN and SN
AiAj/AiAj0.2890.4860.317−0.50.168(SF)
qGYP5.1AiAi/AiAi0.030.5920.510.6820.161PN, SN and SF
AiAi/AiAj0.240.7590.2660.1990.256PN
AiAj/AiAi−0.2640.8421.0420.6230.252PN, SN and SF
AiAj/AiAj−0.0740.3310.3440.9610.095SF

z)DTH: days to heading, PN: panicle number per plant, SN: spikelet number per panicle, SF: fertility, TGW: thousand grain weight.


References
  1. Abdullah, A, Kobayashi, H, Matsumura, I, and Itoh, S (2008). World rice demand towards 2050: Impact of decreasing demand of per capita rice consumption for China and India. Japan and East Asian Regionalism, Hassan, A, Nasrudin, M, and Akhir, M, ed. Kuala Lumpur: East Asian Studies Department, pp. 1-17
  2. Chin, JH, Chu, SH, Jiang, W, Cho, Y-I, Basyirin, R, and Brar, DS (2011). identification of QTLs for hybrid fertility in inter-subspecific crosses of rice (Oryza sativa L.). Genes Genom. 33, 39-48.
    CrossRef
  3. Chin, JH, Kim, JH, Jiang, W, Chu, SH, Woo, MO, and Han, L (2007). Identification of subspecies-specific STS markers and their association with segregation distortion in rice (Oryza sativa L.). J Crop Sci Biotech. 10, 175-184.
  4. Cho, YG, Kang, HJ, Lee, JS, Lee, YT, Lim, SJ, and Gauch, H (2007). Identification of quantitative trait loci in rice for yield, yield components, and agronomic traits across years and locations. Crop Sci. 47, 2403-2417.
    CrossRef
  5. Chu, SH, Jiang, W, Lee, J, Chin, JH, and Koh, HJ (2012). QTL analyses of heterosis for grain and grain yield-related traits in indica-japonica crosses of rice (Oryza sativa L). Genes Genom. 34, 367-377.
    CrossRef
  6. Dan, Z, Hu, J, Zhou, W, Yao, G, Zhu, R, and Huang, W (2015). Hierarchial additive effects on heterosis in rice (Oryza sativa L.). Front Plant Sci. 6, 738.
    CrossRef
  7. Dan, Z, Liu, P, Huang, W, Zhou, W, Yao, G, and Hu, J (2014). Balance between a higher degree of heterosis and increased reproductive Isolation: a strategic design for breeding inter-subspecific hybrid rice. PLoS ONE. 9, e9312.
    CrossRef
  8. Doi, K, Yoshimura, A, and Iwata, N (1998). RFLP mapping and QTL analysis of heading date and pollen sterility using backcross populations between Oryza sativa L. and Oryza glaberrima Steud. Breed Sci. 48, 395-399.
  9. FAORAP and APSA (FAO Regional Office for Asia and the Pacific and Asia-Pacific Seed Association) 2014. Hybrid rice development in Asia. Assessment of limitations and potential, Proceedings of the Regional Expert Consultation, 2–3 July 2014, Bangkok, pp.35-46.
  10. Jiang, L, Zhang, W, Xia, Z, Jiang, G, Qian, Q, and Li, A (2007). A paracentric inversion suppresses genetic recombination at the FON3 locus with breakpoints corresponding to sequence gaps on rice chromosome 11L. Mol Genet Genomics. 277, 263-272.
    CrossRef
  11. Jiang, W, Jin, Y-M, Lee, J, Lee, K-I, Piao, R, and Han, L (2011). Quantitative trait loci for cold tolerance of rice recombinant inbred lines in low temperature environments. Mol Cells. 32, 579-587.
    Pubmed KoreaMed CrossRef
  12. Li, DT, Zhu, SS, Jiang, L, and Wan, JM (2005). Mapping for a new locus causing hybrid sterility in a China landrace (Oryza sativa L.). Rice Genetics Newsletters. 22, 20.
  13. Li, D, Huang, Z, Song, S, Xin, Y, Mao, D, and Lv, Q (2016). Integrated analysis of phenome, genome transcriptome of hybrid rice uncovered multiple heterosis-related loci for yield increase. Proc Natl Acad Sci USA. 113, E6026-E6035.
    CrossRef
  14. Li, L, Lu, K, Chen, Z, Mu, T, Hu, Z, and Li, X (2008). Dominance, overdominance and epistasis condition the heterosis in two heterotic rice hybrids. Genetics. 180, 1725-1742.
    Pubmed KoreaMed CrossRef
  15. Li, ZK, Luo, LJ, Mei, HW, Shu, QY, Tabien, R, and Zhong, DB (2001). Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. I. Biomass and grain yield. Genetics. 158, 1737-1753.
    Pubmed KoreaMed
  16. Li, ZK, Yu, SB, Lafitte, HR, Huang, N, Courtois, B, and Hittalmani, S (2003). QTL × environment interactions in rice. I. heading date and plant height. Theor Appl Genet. 108, 141-153.
    Pubmed CrossRef
  17. Lin, H, Liang, ZW, Sasaki, T, and Yano, M (2003). Fine mapping and characterization of quantitative trait loci Hd4 and Hd5 controlling heading date in rice. Breed Sci. 53, 51-59.
    CrossRef
  18. Luo, LJ, Li, ZK, Mei, HW, Shu, QY, Tabien, R, and Zhong, DB (2001). Overdominant epistatic loci are the primary genetic basis of inbreeding depression and heterosis in rice. II. grain yield components. Genetics. 158, 1755-1771.
    Pubmed KoreaMed
  19. Mei, HW, Li, ZK, Shu, QY, Guo, LB, Wang, YP, and Yu, XQ (2005). Genetic actions of QTLs affecting several agronomic traits resolved in a recombinant inbred rice population and two backcross populations. Theor Appl Genet. 110, 649-659.
    Pubmed CrossRef
  20. Miura, K, Ikeda, M, Matsubara, A, Song, XJ, Ito, M, and Asano, K (2010). OsSPL14 promotes panicle branching and higher grain productivity in rice. Nat Genet. 42, 545-549.
    Pubmed CrossRef
  21. Miyata, M, Yamamoto, T, Komori, T, and Nitta, N (2007). Marker-assisted selection and evaluation of the QTL for stigma exsertion under japonica rice genetic background. Theor Appl Genet. 114, 539-548.
    CrossRef
  22. Mizuta, Y, Harushima, Y, and Kurata, N (2010). Rice pollen hybrid incompatibility caused by reciprocal gene loss of duplicated genes. Proc Natl Acad Sci USA. 107, 20417-20422.
    Pubmed KoreaMed CrossRef
  23. Reflinur, , Chin, JH, Jang, SM, Kim, B, Lee, J, and Koh, HJ (2012). QTLs for hybrid fertility and their association with female and male sterility in rice. Genes Genom. 34, 355-365.
    CrossRef
  24. Rogers, SO, and Bendich, AJ (1988). Extraction of DNA from plant tissues. Plant Mol Biol Man. 6, 1-10.
  25. Shen, B, Zhunag, JY, Zhang, KQ, Xia, QQ, Sheng, CX, and Zheng, KL (2003). QTLs mapping of leaf traits and root vitality in a recombinant inbred line population of rice. Yi Chuan Xue Bao. 30, 1133-1139.
  26. Shen, G, Zhan, W, Chen, H, and Xing, Y (2014). Dominance and epistasis are the main contributors to heterosis for plant height in rice. Plant Sci. 215, 11-18.
    Pubmed CrossRef
  27. Sobrizal, , Matsuzaki, Y, Sanchez, PL, Ikeda, K, and Yoshimura, A (2000). Mapping of F1 pollen semi-sterility gene found in backcross progeny of Oryza sativa L. and Oryza glumaepatula Steud. Rice Genetics Newsletters. 17, 61-62.
  28. Song, X, Qiu, SQ, Xu, CG, Li, XH, and Zhang, Q (2005). Genetic dissection of embryo sac fertility, pollen fertility, and their contributions to spikelet fertility of intersubspecific hybrids in rice. Theor Appl Genet. 110, 205-211.
    Pubmed CrossRef
  29. Stam, P (1993). Construction of integrated genetic linkage maps by means of a new computer package: Join Map. Plant J. 3, 739-744.
    CrossRef
  30. Tan, XL, Vanavichit, A, Amornsilpa, S, and Trangoonrung, S (1998). Genetic analysis of rice CMS-WA fertility restoration based on QTL mapping. Theor Appl Genet. 97, 994-999.
    CrossRef
  31. Tao, Y, Zhu, J, Xu, J, Wang, L, Gu, H, and Zhou, R (2016). Exploitation of heterosis loci for yield and yield components in rice using chromosome segment substitution lines. Sci Rep. 6, 36802.
    Pubmed KoreaMed CrossRef
  32. Tian, F, Li, DJ, Fu, Q, Zhu, ZF, Fu, YC, and Wang, XK (2006). Construction of introgression lines carrying wild rice (Oryza rufipogon Griff.) segments in cultivated rice (Oryza sativa L.) background and characterization of introgressed segments associated with yield-related traits. Theor Appl Genet. 112, 570-580.
    CrossRef
  33. Venu, RC, Jia, JMY, Liu, G, Jia, MH, Nobuta, K, and Sreerekha, MV (2014). Identification of candidate genes associated with positive and negative heterosis in rice. PLoS ONE. 9, e95178.
    Pubmed KoreaMed CrossRef
  34. Wang, DL, Zhu, J, Li, ZK, and Paterson, AH (1999). Mapping QTLs with epistatic effects and QTL x environment interactions by mixed linear model approaches. Theor Appl Genet. 99, 1255-1264.
    CrossRef
  35. Xin, XY, Wang, WX, Yang, JS, and Luo, XJ (2011). Genetic analysis of heterotic loci detected in a cross between indica and japonica rice (Oryza sativa L.). Breed Sci. 61, 380-388.
    CrossRef
  36. Yan, CJ, Liang, GH, Gu, SL, Yi, CD, Lu, JF, and Li, X (2003). Molecular marker analysis and genetic basis for sterility of typical indica/japonica hybrids. Yi Chuan Xue Bao. 30, 267-276.
    Pubmed
  37. Yan, CJ, Liang, GH, Zhu, LH, and Gu, MH (2000). RFLP analysis on wide compatibility genes in rice variety dular of ecotype aus. Yi Chuan Xue Bao. 27, 409-417.
  38. Yu, SB, Li, JX, Xu, CG, Tan, YF, Gao, YJ, and Li, XH (1997). Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci USA. 94, 9226-9231.
    Pubmed KoreaMed CrossRef
  39. Yu, SB, Li, JX, Xu, CG, Tan, YF, Li, XH, and Zhang, Q (2002). Identification of quantitative trait loci and epistatic interactions for plant height and heading date in rice. Theor Appl Genet. 104, 619-625.
    CrossRef
  40. Zhuang, JY, Fan, YY, Wu, JL, Xia, YW, and Zheng, KL (2001). Comparison of the detection of QTL for yield traits in different generations of a rice cross using two mapping approaches. Yi Chuan Xue Bao. 28, 458-464.
    Pubmed