Skip to main content
  • Research article
  • Open access
  • Published:

Genetic dissection of yield-related traits and mid-parent heterosis for those traits in maize (Zea mays L.)

Abstract

Background

Utilization of heterosis in maize could be critical in maize breeding for boosting grain yield. However, the genetic architecture of heterosis is not fully understood. To dissect the genetic basis of yield-related traits and heterosis in maize, 301 recombinant inbred lines derived from 08 to 641 × YE478 and 298 hybrids from the immortalized F2 (IF2) population were used to map quantitative trait loci (QTLs) for nine yield-related traits and mid-parent heterosis.

Results

We observed 156 QTLs, 28 pairs of loci with epistatic interaction, and 10 significant QTL × environment interactions in the inbred and hybrid mapping populations. The high heterosis in F1 and IF2 populations for kernel weight per ear (KWPE), ear weight per ear (EWPE), and kernel number per row (KNPR) matched the high percentages of QTLs (over 50%) for those traits exhibiting overdominance, whereas a notable predominance of loci with dominance effects (more than 70%) was observed for traits that show low heterosis such as cob weight per ear (CWPE), rate of kernel production (RKP), ear length (EL), ear diameter (ED), cob diameter, and row number (RN). The environmentally stable QTL qRKP3–2 was identified across two mapping populations, while qKWPE9, affecting the trait mean and the mid-parent heterosis (MPH) level, explained over 18% of phenotypic variations. Nine QTLs, qEWPE9–1, qEWPE10–1, qCWPE6, qEL8, qED2–2, qRN10–1, qKWPE9, qKWPE10–1, and qRKP4–3, accounted for over 10% of phenotypic variation. In addition, QTL mapping identified 95 QTLs that were gathered together and integrated into 33 QTL clusters on 10 chromosomes.

Conclusions

The results revealed that (1) the inheritance of yield-related traits and MPH in the heterotic pattern improved Reid (PA) × Tem-tropic I (PB) is trait-dependent; (2) a large proportion of loci showed dominance effects, whereas overdominance also contributed to MPH for KNPR, EWPE, and KWPE; (3) marker-assisted selection for markers at genomic regions 1.09–1.11, 2.04, 3.08–3.09, and 10.04–10.05 contributed to hybrid performance per se and heterosis and were repeatedly reported in previous studies using different heterotic patterns is recommended.

Background

Utilization of heterosis in maize is of great importance for boosting grain yield [1,2,3,4]. Hybrids in maize accounted for 65% of total maize cultivation by the late twentieth century and had contributed to a quadrupling of annual maize production [1]. Numerous breeders have been interested in heterosis for many years; however, comprehension of the genetics implicated in heterosis for grain yield still remains elusive. Therefore, investigating and assessing genetic mechanisms of maize heterosis for grain yield would provide a paved route for understanding those phenomena and help to optimize breeding for grain yield in different heterotic groups.

Through earlier inbred selection and experimental breeding [5, 6], maize breeders realized the importance of germplasm resources and heterosis. Previous studies [2, 7,8,9,10,11,12,13] used molecular markers, phenotypic identification, and pedigree information to divide different germplasm into different heterotic groups. The heterotic pattern Stiff Stalk Synthetic (SS) × Non-Stiff Stalk (NSS) has been widely used for maize breeding [12, 14,15,16]. Li and Wang [17], based on former studies, proposed five heterotic groups, namely Tangsipingtou (TSPT), Lancaster, Lancaster-like, the improved Reid (PA), and Tem-tropic I (PB). Extensive research has been focused on the TSPT × Reid heterotic pattern because it is one of the predominant heterotic patterns in northern China and has been used for the production of elite hybrids such as Zhengdan958 [18]. However, this heterotic pattern is rarely used in southwestern China, where alternative, scarcely investigated, heterotic patterns such as PB × PA are more suitable [17]. Therefore, the current study investigated the genetics of heterosis for yield-related traits in a PB × PA hybrid and the genetics of the yield-related traits per se.

Heterosis is defined as the superiority of F1 hybrids in their performance over their parents [19, 20] and has been attributed to three different genetic hypotheses: dominance [21, 22], overdominance [19, 23], and/or epistasis [24, 25]. Specific loci for heterosis reported in previous studies corroborated each one of these hypotheses [26,27,28,29,30,31,32,33,34,35,36,37,38]. Those loci were identified using different designs such as the North Carolina design III with F3 families or recombinant inbred lines (RILs), immortalized F2 (IF2) populations, or chromosome segment substitution lines (CSSLs). Those studies suggested different degrees of dominance for grain yield and inconsistency of heterotic loci across different heterotic patterns. For example, Frascaroli et al. [29] used the RILs from a cross between two elite inbred lines, B73 and H99, and three testcross populations obtained by crossing 142 RILs to each parent and F1 found a significant correlation between heterozygosity level and phenotypic performance. Schön et al. [32] demonstrated high congruency of QTLs for heterosis of grain yield in three populations derived from the crosses of the heterotic pattern Iowa Stiff Stalk × Lancaster Sure Crop and proposed that diverse alleles were fixed in each heterotic group. Tang et al. [31] identified 13 heterotic loci for yield-related traits, and Guo et al. [34] reported that the genetic basis of heterosis for grain yield of the elite maize hybrid Yuyu22 relies on the cumulative effects of dominance, overdominance, and epistasis. Li et al. [37] detected 38 heterotic loci for ear-related traits and suggested a genotype-dependent combination of heterotic loci in different hybrids. Notably, Wang et al. [35, 38] identified 169 heterotic loci associated with grain yield and yield components in two CSSL test populations. Therefore, the results pertaining to genetics of grain yield and heterosis using crosses from well-studied heterotic patterns, such as Reid × Lancaster and TSPT × Reid, cannot be extrapolated to other heterotic patterns such PA × PB.

In this study, in order to dissect genetic architecture of grain yield-related traits in maize and the contribution of different genetic effects to heterosis for grain yield-related traits, we evaluated 301 RILs from 08 to 641 × YE478 in four environments and an IF2 population produced from those RILs under three environments. The elite line YE478, belonging to the PA heterotic group [17], is one of the hybrid foundation parents in China. The other parent, 08–641 (SAU08–641), which is derived from PB (or Tem-tropic I heterotic group) germplasm, is a typical elite inbred line widely used in southwest China since 1998. PA × PB is one of the predominant heterotic patterns in southwestern China. The main objectives of this research were (1) to detect and estimate QTLs with additive and additive-additive effects on yield-related traits using an RIL population (2) to identify QTLs with additive, dominant, and/or epistatic effect on yield-related traits using a hybrid population derived from the RIL population (IF2 population) (3) to identify QTLs for mid-parent heterosis (MPH) using both mapping populations and (4) to identify the most suitable genomic regions associated with grain yield and MPH to be managed by marker-assisted selection.

Results

The performance of grain yield-related traits and heterosis across environments

The parental line Ye478 had significantly higher mean RKP and kernel number per row (KNPR) compared with those of 08–641 (Additional file 1: Figure S1). The RIL and IF2 populations showed considerable variations for yield-related traits (Fig. 1). The means of the F1 were significantly higher than the means of any of the parental lines for ear weight per ear (EWPE), cob weight (CWPE), ear length (EL), ear diameter (ED), kernel weight per ear (KWPE), KNPR, and RKP. MPH of the hybrid Ye478 × 08–641 was < 16% for ED and RKP, 43.5% for EL, and > 70% for EWPE, KWPE, CWPE, and KNPR. Average MPH values across hybrids of the IF2 population were similar to MPH values for F1 (Additional file 4: Table S2). These results suggested a strong heterosis level for EWPE, KWPE, KNPR, and CWPE. The RIL and IF2 populations thus seem suitable for dissection of genetic effects involved in yield-related traits and heterosis for these traits.

Fig. 1
figure 1

The histogram for yield-related traits and mid-parent heterosis (MPH) for those traits in the RILs and the IF2 population. EWPE, ear weight per ear; CWPE, cob weight per ear; EL, ear length; ED, ear diameter; CD, cob diameter; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per ear; RKP, rate of kernel production

In the combined analyses of variance of the RIL and IF2 populations across all environments, the sources of variation for genotype, environment, and genotype × environment interaction were highly significant (p < 0.01) for most yield-related traits (Table 1). The broad-sense heritability (h2) for these traits ranged from 60.8% for KWPE to 91.1% for RN. Heritability for EWPE and KWPE was greater in the IF2 population than in the RIL population. The correlations among most yield-related traits and among MPH values for those traits were significant (Fig. 2a; Additional file 5: Table S3). EWPE in the RIL and IF2 population was moderately correlated with CWPE, KWPE, EL, ED, and KNPR (Pearson’s correlation coefficients > 0.5). MPH for EWPE was moderately and positively correlated with MPH for CWPE, KWPE, EL, ED, and KNPR. The hierarchical cluster analysis based on the standardized data from the RILs classified the yield-related traits into two unrooted groups, with traits within each group being, in general, moderately to highly correlated based on Pearson’s correlation coefficients (Fig. 2b).

Table 1 Mean squares from the combined analyses of variance in the RILs under four environments and the IF2 population across three environments
Fig. 2
figure 2

The correlation relationship among yield-related traits, and mid-parent heterosis (MPH) for those traits, and QTL distribution in the different datasets. a Pearson’s correlation coefficients among yield-related traits and MPH (|r2| > = 0.5). The yellow lines indicate significant correlation found in the MPH dataset; Blue line indicates significant correlation in the RILs; green line indicates significant correlation in both the IF2 and the MPH dataset; black lines indicate significant correlations found across the RILs, the IF2 and the MPH dataset. The number below each trait indicates the number of the correlated traits. b Hierarchical cluster analysis (“hclust”) of phenotypic traits based on the standardized data in the RILs. The short-scale indicates Euclidean distance. c Frequency distribution of QTLs with PVE (phenotypic variance explained by each QTL) detected in the RILs, IF2 population, and the MPH dataset. The vertical and dashed lines indicate the median values. d The QTL numbers for nine yield-related traits identified in the RILs, IF2 population, and the MPH dataset. e Venn diagram for the 156 QTLs detected using averaged means for yield-related traits (EWPE, ear weight per ear; CWPE, cob weight per ear; EL, ear length; ED, ear diameter; CD, cob diameter; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per ear; RKP, rate of kernel production) across environments in the RILs and IF2 populations, and the MPH dataset

Identification of QTLs

In the current study, 156 QTLs were found for averaged yield-related traits and MPH for those traits across all environments (Table 2; Fig. 3). Ninety-six QTLs detected across environments were also identified in at least one environment via single environment analysis (Additional file 6: Table S4). These QTLs were distributed on all chromosomes, although two QTL–concentrated regions were identified at bin 3.08–3.09 and 10.04–10.05. Each QTL explained 0.38 to 33.99% of phenotypic variance (PVE), and over half of QTLs accounted each for less than 5% of PVE (Fig. 2c). The QTL number and the total PVE greatly differed among traits and populations (Table 2; Fig. 2d, e). Nine QTLs were detected simultaneously in the RIL and IF2 populations, while 13 QTLs detected in the IF2 population and in the MPH dataset co-localized. Only two QTLs found in the RIL population were also positioned in the same regions when QTL analyses were performed for MPH. These results suggested that genetics involved in hybrid performance per se could be, in part, related to that implicated in MPH and it is in agreement with the moderate to high correlation coefficients between hybrid performance per se in the IF2 population and the corresponding MPH levels reported for all traits, except for RKP (Additional file 7: Table S5).

Table 2 Main features of QTLs for yield-related traits detected in the RILs, the IF2 population, and for the mid-parental heterosis (MPH)
Fig. 3
figure 3

Distribution of QTLs for yield-related traits detected in the RIL and IF2 populations and for heterosis. The red rectangles represent QTLs with negative effects (additive in the RIL population, dominant and additive in the IF2 population and dominance effects for mid-parent heterosis (MPH)), while the blue rectangles represent positive effects. The width of the rectangles indicates the confidence interval of QTLs. The heatmap of QTL density across ten chromosomes was indicated at the bottom of the figure and below the X-axis. EWPE, ear weight per ear; CWPE, cob weight per ear; EL, ear length; ED, ear diameter; CD, cob diameter; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per ear; RKP, rate of kernel production

Dominance and overdominance effects play crucial roles in heterosis for yield-related traits

The contribution of dominance and overdominance to the heterosis phenomena was assessed by the (|D/A|) ratio for each locus [34, 39] using single-marker analysis (SMA) and inclusive composite interval mapping (ICIM) (Fig. 4). Overdominance was detected for over 60% of genomic markers for EWPE and KWPE in both the IF2 and MPH/RIL datasets, more than 50% for KNPR and EL, over 40% for CWPE and ED, over 30% for CD and RKP, and over 15% for RN. Based on SMA, the number of significant markers in the IF2 population showing overdominance accounted for 66.6% for KNPR, more than 40% for EWPE and KWPE, more than 15% for EL, RKP, and ED, and less than 2% for CWPE, CD, and RN, respectively. Additionally, approximately 78% of QTLs for KWPE exhibited overdominance, as did 72.7% of QTLs for EWPE, 50% of QTLs for KNPR, and less than 30% of QTLs for the remaining traits. Among traits, we also found a strong association between the proportion of QTLs with overdominance and the ratio of MPH/mid-parent for F1 (r2 = 0.951***) or the MPH level across hybrids from the IF2 population (r2 = 0.966***). We found, in general, significant but low correlations between performance per se or MPH of individuals in the IF2 population and their percentages of heterozygous marker loci (Additional file 8: Table S6).

Fig. 4
figure 4

Distribution of markers in the whole genome and significant loci for yield-related traits with dominance and over-dominance effect. The ratio D/A was calculated in two ways: both additive and dominance genetic effects (A and D) were estimated from the IF2 population (IF2) or additive effects (A) were estimated from RILs and dominance effects (D) from the MPH dataset (MPH/RIL). SMA, single-marker analysis; ICIM, inclusive composite interval mapping. PO, PD, ND, and NO, indicate positive overdominance, positive dominance, negative dominance, and negative overdominance effects, respectively. EWPE, ear weight per ear; CWPE, cob weight per ear; EL, ear length; ED, ear diameter; CD, cob diameter; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per ear; RKP, rate of kernel production

Analysis of EPI and QEI

We estimated the total PVE by the epistatic loci in the whole-genome search using ICIM-EPI (Additional file 2: Figure S2). In the IF2 population, epistatic loci had higher PVE than single-locus QTLs for CWPE (50% vs 41%), ED (24% vs 20%), CD (56% vs 32%), and KNPR (53% vs 37%). In the MPH dataset, the total PVE by epistatic loci was greater than by single-locus QTLs for all traits except for EWPE (20% vs 49%), CWPE (23% vs 28%), KWPE (43% vs 52%), and RKP (16% vs 43%), indicating that EPI had a greater impact on these traits than single-locus QTLs, regardless of the statistical significance. The current findings also suggest that the cumulative effects resulting from combining the epistatic loci and single-locus QTLs explained a large proportion of PVE. Meanwhile, based on mixed-model-based composite interval mapping (MCIM) with the full QTL model, 28 pairs of marker intervals showed significant epistatic interaction effects; 13 additive by additive (AA) interactions were detected in the RIL population, 10 epistatic interactions in the IF2 population, and five digenic interactions in the MPH dataset (Additional file 9: Table S7). We observed 23 pairs of significant AA interactions, three for additive by dominance (AD) interactions, seven for dominance by additive (DA) interactions, and six for dominance by dominance (DD) interactions. Each marker interval interacted with one to five loci. In addition, a total of nine loci related to four traits exhibited significant additive × environment interactions in the RIL and IF2 populations, seven in the RILs, and two in the IF2 population; only one locus, qKWPE2–2, for the MPH level for KWPE showed a significant dominance × environment interaction in the IF2 population (Additional file 10: Table S8). Many of these marker intervals involved in epistasis and QEI also contained QTLs with main effects on yield-related traits or MPH. For instance, the interval PZE-101184757/SYN25826 was involved in three significant AA interactions and contains qRKP1–1.

QTL co-localization across traits

Overlapping QTLs for different traits were assigned to the same QTL cluster. In total, 95 QTLs were grouped into 33 QTL clusters distributed across all chromosomes (Table 3). Each QTL cluster contained QTLs for two to six traits. The cluster QC3–3 harboring six QTLs affected the means of EL, RN, ED, KNPR, RKP, and KWPE, as well as the MPH for RKP. QC10–3, within a 5 cM region on chromosome 10, had a pleiotropic effect on EWPE, KNPR, CWPE, and the MPH level of ED and CD. As expected from the moderate to high correlation between hybrid performance per se and MPH, some loci associated with yield-related traits and MPH gathered together. Nearly 60% of QTL clusters were associated with two traits. Some of these clusters were specific for a determined population or for MPH. For example, clusters QC1–5 (bin 1.09–1.11) affected mid-parent heterosis for EWPE and KWPE. Contrarily, many other clusters comprised QTLs for yield-related traits and heterosis for those traits. Cluster QC1–4 (bin 1.08) affected KWPE and MPH for RKP. Genomic region 3.09 affected KNPR, RKP, and the MPH level of KWPE and EWPE. Interestingly, some of these clusters co-located in the interval of loci showing epistasis and QTL × environment interactions. For instance, QC7–1 (bin 7.02) had a pleiotropic effect on RKP, EWPE, KNPR, KWPE, and the MPH level for KWPE. This locus interacted with other loci exhibiting epistatic effects for RN in the RILs and for KWPE in the IF2 population. This hotspot also interacted with E1 (2014JH) for RKP in the RILs. These results further confirmed the complex genetic architecture of yield-related traits.

Table 3 QTL clusters for yield-related traits and for heterosis for those traits detected in the RILand IF2 populations

Discussion

Heterotic loci with dominance and epistasis and the complex genetic network of “crosstalk” contribute to grain yield

Previous studies [29, 31, 34, 40, 41] proposed that better-parent heterosis was yield-related loci displaying dominance, overdominance, and/or epistasis. In the present study, over half of QTLs exhibited overdominance for KWPE, EWPE, and KNPR, which is similar to previous results reported by Frascaroli et al. [29] for overdominant loci for grain yield and kernel number per plant and by Guo et al. [34] for single markers and QTLs for grain yield. Such a high average degree of overdominance for grain yield and kernel number per plant could be, in part, explained by linkage in the repulsion of loci of the dominance type in mapping populations [29]. In accordance with this, greater number of genome markers with overdominance for EL and ED was previously observed in the F2:3 population, from which the current RIL population has been derived [42], than in the present mapping populations. The degree of dominance for the traits was well in accordance with the level of F1 superiority over mid-parent performance and the type of inheritance of different traits varied among those traits, which was consistent with other studies [29, 33, 34, 37, 43]. High correlation coefficients were observed between the proportion of QTLs displaying overdominance and the MPH levels across different traits. For example, the F1 did not differ from the inbred parent lines 08–641 and Ye478 for RN and CD, and accordingly QTLs, such as qRN1–2 and qRN7–2 displaying similar additive effects, detected for those traits showed low or no significant dominance effects and were mainly detected in RIL and IF2 populations [34]. The QTLs for EWPE, KNPR, and KWPE in IF2 population and the MPH dataset displayed high dominance effects. In general, similar dominance effects were found in the IF2 population and in the MPH dataset, but epistatic effects could interfere in the estimation of dominance effects and cause discrepancies between results obtained in different populations [44]. Discrepancies between additive effects estimated in the RIL and IF2 populations were even more conspicuous. For example, additive effects estimations for qCWPE6 were different when obtained from RIL or IF2 evaluations. These findings demonstrated that the loci displaying dominance contributed to yield-related traits and heterosis for those traits, although overdominance had an outstanding contribution to heterosis for EWPE, KWPE, and KNPR, which presented the highest values for MPH.

The level of whole marker loci with heterozygosity and phenotypic variations (hybrid performance per se or MPH) were weakly correlated, whereas Frascaroli et al. [29] and Larièpe et al. [33] demonstrated moderate-to-high correlations between them. Such low correlations indicated that the overall heterozygosity contributed little to hybrid performance per se or heterosis [45, 46]. Huang et al. [46] demonstrated that phenotypic variations for yield-related traits are highly correlated with the number of accumulated superior gene alleles. In contrast, numerous epistatic effects and QTL × environment interactions for yield-related traits have been addressed in the literature [29, 31, 34, 37, 47,48,49,50]. The total PVE values by single-locus QTLs and epistasis together approach the high heritability for these yield-related traits [34], indicating the cumulative effects of different genetic effects [34, 41]. The QTL model, which was used to estimate the loci showing additive, dominance, and EPI effects, identified 28 pairs of significant loci with epistatic effects and 10 loci with QTL × environment interaction in the two mapping populations. Of those, only five pairs of epistatic loci were observed for MPH. In addition, many of these loci showing interactions did not co-map with putative QTLs. Similar features were also reported by Yang et al. [51]. EPI and QEI are also involved in the complex genetic web for grain yield via “cross-talk” among loci [45, 47, 48, 52]. Collectively, the cumulative effects of various genetic effects discussed above seem to explain the genetic mechanism of grain yield and better-parent performance [34, 41].

QTL hotspots regulate the architecture of yield-related traits and heterosis

Over 60% of QTLs gathered together and were integrated into 33 clusters distributed on 10 chromosomes. We found a good consistency for QTL clusters between our study and former studies (Table 3). Thirteen pleiotropic regions in the present study were consistently mapped to the clusters in all chromosomes, except in chromosome 4, that were defined by Chen et al. [53] in their meta-analysis of yield-related QTLs. Notably, many of these clusters related to MPH co-located with QTLs found in previous studies [29, 31, 34, 35] for better parent performance (Table 3). Cluster QC1–2 associated with EL, CD, and CWPE overlapped with the hotspot region for ear weight, grain weight, and KNPR between 38 and 42 Mb on chromosome 1, which was highlighted by Guo et al. [49], and with a genomic region containing QTLs for heterosis for grain yield, ear weight, kernel weight, and ear length reported by Wang et al. [35]. The heterosis locus QC1–5 (1.09–1.11) for EWPE and KWPE co-located with QTLs for heterosis for EWPE and KRP presented by Wang et al. [35]; with a heterosis locus for grain yield at approximately 271 Mbp found by Tang et al. [31] using the heterotic pattern Mixed × Flint and by Larièpe et al. [33] using the pattern European flint × Dent; and with QTLs for 100-kernel weight and grain yield in an IF2 population reported by Guo et al. [34]. The genomic region 2.04 (QC2–2) controlling EWPE and the MPH level for EWPE and KWPE co-located with a QTL involved in heterosis for EL, RN, KWPE, KRP, EWPE, and grain yield by Wang et al. [35], and it was nearby the region implicated in heterosis for grain yield in the Reid × Lancaster pattern reported by Frascaroli et al. [29] and the meta-QTL for grain yield, ear-related traits, and kernel-related traits highlighted by Chen et al. [53]. Cluster QC3–3, being particularly interesting in the current study because of its association with heterosis for several traits, overlapped with a QTL hotspot for heterosis for EL, EWPE, RN, KRP, KWPE, and grain yield located between 209 and 213 Mbp on chromosome 3 reported by Wang et al. [35]. Genomic region 3.09 (QC3–4), related to KNPR and RKP and the MPH level for KWPE and EWPE overlapped with the environmentally stable region for heterosis for KNPR, ear weight, and kernel weight between 218 and 228 Mbp on chromosome 3 reported by Li et al. [37] in the heterotic pattern Reid × TSPT. Another hotspot region, QC10–3, related to five traits co-mapped with a heterotic locus that affects grain yield and kernel number presented by Frascaroli et al. [29]. In addition, genomic regions 2.04, 7.02, and 9.03, associated with heterosis for KWPE, contain QTLs for grain yield heterosis in the Reid × Lancaster heterotic pattern published by Schön et al. [32]. The cluster QC6–3 for RKP, KWPE, and the MPH level for EL overlapped with the heterotic loci for grain yield found by Samayoa et al. [36] in the population derived from an American dent × European flint cross. These results suggested the consistency of hotspot regions across different heterotic patterns, which may be used for genetic improvement of grain yield in the future.

In contrast, hybrid performance per se for yield-related traits and the MPH level for those traits were moderately to highly correlated and 13 QTLs found in the IF2 population co-localized with QTLs for MPH. These findings indicated the likely common genetic basis between grain yield and its MPH level. The results also partly supported that heterosis loci for grain yield are not independent from loci for grain yield, thus agreeing with Guo et al. [34] and refuting the conclusions by Tang et al. [31] who proposed independence between genomics involved in performance per se and heterosis.

Genetic basis for grain yield and its components in different populations from the same parents

Our study revealed that QTLs found in different populations derived from the same cross are not always consistent with previous studies [34, 54,55,56]. Among the 156 detected QTLs for yield-related traits and heterosis, only one QTL, qRN10–2, was previously found across F2:3 families derived from the cross 08–641 × YE478 [60]. In addition, only nine QTLs with |D/A| < 1 were simultaneously detected in the RILs and the IF2 population. Lack of congruency for QTLs found in populations derived from the same cross but with different homozygous levels could be due to the different genetic effects acting in each population, but also to the different environments used for evaluations and biases caused by population sampling [29].

Conclusions

Heterosis for yield-related traits in maize is one of the main issues for maize breeders, and investigation and assessment of heterotic loci involved in specific heterotic patterns are of particular interest. The inheritance of yield-related traits and MPH varied among different traits. A large proportion of the loci with dominance effect had a greater effect on most traits compared with loci with other genetic effects, whereas overdominance also contributed greatly to MPH for KNPR, EWPE, and KWPE in the current heterotic pattern. As QTL hotspots at 1.09–1.11, 2.04, 3.08–3.09, and 10.04–10.05 encompass genomic regions where other authors have found QTLs for heterosis across different heterotic patterns, markers in those regions could be used in marker-assisted selection programs for increased yield.

Methods

Plant material and field experiments

The mapping populations used in this study consisted of 301 RILs derived from a cross between the maize inbred lines 08–641 and Ye478. The elite inbred 08–641 comes from the Maize Research Institute of Sichuan Agricultural University, and Ye478 was provided by the Project of National Major Basic Dairy Research “973” Plan, which is developed by the Laizhou Academy of Agricultural Sciences [17]. The parental line 08–641 was crossed with Ye478 in the winter of 2010 to produce a set of RILs (F8) by using the single seed descent method. The 301 RILs and the two parental lines were evaluated in four environments. Three environments were located at the Xishuangbanna maize breeding base of the Maize Research Institute of Sichuan Agricultural University, Jinghong (JH 21°95 N, 100°76E; average daytime temperature 18.6–21.9 °C; average rainfall per year 1200–1700 mm; tropical humid monsoon climate), in Yunnan province, in March 2014 (JH2014), March 2015 (JH2015), and March 2016 (JH2016), and one environment was located at the Modern Agriculture Research and Development Center of Sichuan Agricultural University, Chongzhou (CZ 30°33′N, 103°38′E; average daytime temperature 15.9 °C; average rainfall per year 1012.4 mm; sub-tropical humid monsoon climate), in Sichuan province, April 2016 (CZ2016). An IF2 population was developed with 320 RILs following the procedure of Hua et al. [57] and Tang et al. [31] The 320 RILs were randomly divided into two groups of 160 RILs. Then, the lines of both groups each RIL was only involved in a single cross, yielding 160 F1 crosses. We repeated the procedure mentioned above twice and constructed 320 crosses that constituted the IF2 population. After genotyping the RILs, those with heterozygosity above 20% and their respective crosses were removed. After filtering the RILs, the IF2 population consisted of 298 crosses that were evaluated under three environments (2015JH, 2016JH, and 2016CZ). In each environment, each trial was performed following a randomized complete block design with two repetitions. Plant density was 57,000 plants per ha. Each plot consisted of a single row. Each row with 14 plants was 3 m long, with 0.80 m space between rows. The trials with two mapping populations were adjacent and occupied a uniform and square parcel. Field management was in accordance with local practices.

Phenotypic measurements and analysis

In each plot, 10 well-pollinated ears were randomly chosen from all ears harvested. The nine yield-related traits measured (Additional file 3: Table S1) were: EL (cm), ED (mm), CD (mm), RN (count), KNPR, EWPE (g), CWPE (g), KWPE (g), and RKP. Means for each replication were used for further analyses. The distribution of the phenotypic traits and the Pearson’s correlation, using PROC CORR from the statistical software package SPSS 17.0 (SPSS, Inc., Chicago, IL, USA), were conducted with the means of the phenotypic traits in the 301 RILs and 298 hybrids of the IF2 population across all environments. A hierarchical cluster analysis (“hclust”) was carried out for phenotypic traits, based on the standardized data from the RILs [58, 59]. Combined analyses of variance of the RILs and the IF2 population were computed for each trait using the GLM procedure in SPSS 17.0 (SPSS, Inc., Chicago, IL, USA) with genotype as a fixed effect and replication and environment as random effects. Broad-sense heritability (h2) for the RILs and IF2 population was computed on an entry mean basis as described by Hallauer and Miranda [60]:

$$ {h}^2={\sigma}_g^2/\left({\sigma}_g^2+\frac{\sigma_{ge}^2}{n}+\frac{\sigma^2}{nb}\right), $$
(1)

where \( {\sigma}_g^2 \) represents the genetic variance, \( {\sigma}_{ge}^2 \) represents the genotype × environment interaction variance, σ2 represents the error variance, b represents the number of replications, and n is the number of environments. The 90% confidence interval of heritability (h2) was determined according to Knapp et al. [61]. Each phenotypic trait of the 301 RILs and each trait of the IF2 population were analyzed following two mixed models fitted by restricted maximum likelihood:

$$ {Y}_{mk}=\mu +{G}_m+{R}_k+{\varepsilon}_{mk}, $$
(2)
$$ {\mathrm{Y}}_{\mathrm{mik}}=\upmu +{G}_m+{GE}_{mi}+{E}_i+{R}_k+{\varepsilon}_{mi k}, $$
(3)

where Ymk is the phenotypic value of genotype m in replication k; Ymik indicates the phenotypic value of genotype m in environment i and replication k; μ is the overall mean of the RIL population or the IF2 population; Gm is the random effect of genotype m; GEmi is the random effect of the interaction between genotype m and environment i; Ei is the random effect of environment i; Rk is the random effect of replication k; and εmk and εmik denote the random experimental error. The best linear unbiased predictor values (BLUPs) with Eqs. (2) and (3) were used for single and combined environment analysis, respectively. The adjusted mean BLUP values of each RIL across four environments and each hybrid across three environments were used for QTL mapping. The analyses were conducted with the R program for statistical computing [59, 62]. Besides, MPH was estimated as [40]:

$$ {\mathrm{MPH}}_{12}={\mathrm{F}}_{12}-\left({\mathrm{P}}_1+{\mathrm{P}}_2\right)/2 $$
(4)

where F12 is the genotypic value of the F12 individual from the IF2 population, and (P1 + P2)/2 is the average of BLUP values of the corresponding parents (P1 and P2 from the RIL population) estimated from the RIL evaluation. The genotypic dataset for MPH depends on dominance and epistatic effects, but epistatic effects were not included and an additive-dominance model was assumed for the IF2 population [31, 34].

Molecular linkage map construction

Following the modified CTAB protocol [63], DNA was isolated from 7-day-old seedling leaves of the 301 RILs and the parents grown in shade (10 plants per RIL as a bulk). The oligonucleotide pool assay consisted of 3072 well-distributed, high-quality single nucleotide polymorphisms (SNPs) from all 10 maize chromosomes that were selected from 56,110 SNPs in 513 maize inbred lines developed by the National Maize Improvement Center of China using Illumina GoldenGate technology. The protocol for genotyping SNPs using an Illumina BeadStation 500 G (Illumina, San Diego, CA, USA) was conducted as described by Fan et al. [64]. Marker data from the RILs were filtered for heterozygous data points (< 20%), missing data points (< 20%), and segregation distortion (in accordance with the expected Mendelian segregation ratio of 1:1). A total of 683 SNPs selected for their uniform distribution throughout all 10 maize chromosomes were used to construct the linkage map with a total genetic length of 1786.1 cM and an average interval distance of 2.61 cM. The genotypes of each cross from the IF2 population were deduced from the marker genotypes of their RIL parents, with QTL mapping for the IF2 population and the MPH dataset performed using the molecular linkage map of the RIL population [31, 57]. The genetic map for the RILs and IF2 population was developed using MapDisto 1.7.5 [65]. The Kosambi mapping function was used for converting recombination frequencies to genetic distances [66]. The average proportion of homozygous loci for YE478 variants, homozygous loci for 08–641 variants, and heterozygous loci in each individual of the IF2 population was 25.2, 24.7, and 50.1%, respectively, according to the expected ratio in an F2 population.

QTL mapping

QTL mapping for each trait in the RIL and IF2 populations and for MPH using both populations was performed using QTL ICIMapping software [67] with ICIM [68, 69]. For QTL detection on ICIM-ADD, the largest p-value for entering variables in stepwise regression of residuals on marker variables, PIN was set to 0.001. The largest p-value for removing variables is assumed to be two times the PIN. The threshold logarithm of the odds (LOD) score to declare significant QTL was calculated using the 1000 permutation test at Type I error of 0.05 with the step size of 1 cM was 2.5 for the RILs, IF2 population, and the MPH dataset since very minor differences were observed in LOD score among them. Under the additive-dominance model assumed, only QTLs with dominant effect were detected for MPH [34]. We also performed single marker analysis (SMA) to detect QTLs with the same empirical threshold LOD mentioned previously, and the degree of dominance was calculated as a ratio of dominance to additive effects (|D/A|) for each QTL [34]. Loci with |D/A| values greater than 1.26 displayed overdominance; otherwise, the QTL was treated as dominance type [39]. The QTLs were nominated as follows: a lowercase letter “q” standing for QTL, the trait abbreviation, one number standing for the chromosome and another for the physical order or the QTL within the chromosome. The QTLs for different traits detected in overlapping or adjacent intervals were regarded as a QTL cluster. The graphical presentation of linkage maps and QTLs were generated with the R program [59]. To scan for significant marker interval interactions, epistasis interaction (EPI, digenic interactions; AA, AD, DA, and DD) was identified with the step in scanning of 5 cM and PIN = 0.0001 on ICIM-EPI. An empirical threshold LOD for digenic epistasis interaction was set at 5. Meanwhile, the QTL Network program 2.1 based on MCIM [70,71,72] was also used to identify QTL with AA effects in the RILs; AA, AD, DA, and DD effects in the IF2 population using combined data across all environments and the MPH dataset using the average data across environments. QTL × environment interaction effects (QEI) in the RIL and IF2 population were estimated for yield-related traits via combined analysis across all environments. The testing window and filtration window size was set at 10 cM with a walk speed set to 2 cM for one-dimensional genome scans. The F-statistic and the critical F-value were estimated with the help of 1000-permutations test [73] with p < 0.05 as experiment-wise significance level for candidate interval selection. The final full QTL model for each trait incorporated significant additive, dominance, and epistatic effects, as well as their interactions with environments. Bonferroni correction was used to compute the comparison-wise significance threshold assuming an experiment-wise error < 0.05 [74]. We also calculated the percentage of heterozygous loci for each RIL analyzed the correlation relationships between phenotypic performance and heterozygosity, and conducted regression analyses of phenotypic performance on heterozygosity [59].

Availability of data and materials

Supporting data are available in Additional file 11, and materials are available from the authors upon request.

Abbreviations

AA:

Additive × additive interaction

AD:

Additive × dominance interaction

BLUP:

Best linear unbiased predictor

CD:

Cob diameter

CWPE:

Cob weight per ear

DA:

Dominance × additive interaction

DD:

Dominance × dominance interaction

ED:

Ear diameter

EL:

Ear length

EPI:

Epistatic interaction

EWPE:

Ear weight per ear

ICIM:

Inclusive composite interval mapping

IF2 :

The immortalized F2

KNPR:

Kernel number per row

KWPE:

Kernel weight per ear

MCIM:

Mixed-model-based composite interval mapping

MPH:

Mid-parent heterosis

QTL:

Quantitative trait loci

RIL:

Recombinant inbred lines

RKP:

Rate of kernel production

RN:

Row number

SMA:

Single marker analysis

SNP:

Single nucleotide polymorphism

References

  1. Duvick DN. Heterosis: feeding people and protecting natural resources. In: Coors JG, Pandey S, editors. Proceedings of the International Symposium on the Genetics and Exploitation of Heterosis in Crops, CIMMYT, Mexico City, 17–22 Aug 1997. Madison: ASA, CSSA, and SSSA; 1999. p. 19–29.

    Google Scholar 

  2. Duvick DN. Long-term selection in a commercial hybrid maize breeding program. In: Janick J, editor. Plant breeding reviews 24, part 2. Hoboken: Wiley; 2004. p. 109–51.

    Google Scholar 

  3. Duvick DN. The contribution of breeding to yield advances in maize. Adv Agron. 2005;86:83–145.

    Article  Google Scholar 

  4. Ci X, Li M, Liang X, Xie Z, Zhang D, Li X, et al. Genetic contribution to advanced yield for maize hybrids released from 1970 to 2000 in China. Crop Sci. 2011;51:13–20.

    Article  Google Scholar 

  5. Moll RH, Lonnquist JH, Velez FJ, Johnson C. The relationship of heterosis and genetic divergence in maize. Genetics. 1965;52:139–44.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Moll RH, Salhuana WS, Robinson HF. Heterosis and genetic diversity in variety crosses of maize. Crop Sci. 1962;2:197–8.

    Article  Google Scholar 

  7. Smith OS, Smith JS, Bowen SL, Tenborg RA, Wall SJ. Similarities among a group of elite maize inbreds as measured by pedigree, F1grain yield, grain yield, heterosis, and RFLPs. Theor Appl Genet. 1990;80:833–40.

    Article  CAS  Google Scholar 

  8. Beck DL, Vasal SK, Crossa J. Heterosis and combining abilityamong subtropical and temperate intermediate-maturity maize germplasm. Crop Sci. 1991;31:68–73.

    Article  Google Scholar 

  9. Messmer MM, Melchinger AE, Boppenmaier J, Herrmann RG, Brunklaus-Jung E. RFLP analysis of early maturing European maize germplasm. I. Genetic diversity among flint and dent inbreds. Theor Appl Genet. 1992;83:1003–12.

    Article  CAS  Google Scholar 

  10. Messmer MM, Melchinger AE, Herrmann R, Boppenmaier J. Relationship among early maturing European maize hybrids. II Crop Sci. 1993;33:944–50.

    Article  Google Scholar 

  11. Troyer AF. Background of U.S. hybrid corn. Crop Sci. 1999;39:601–26.

    Article  Google Scholar 

  12. Troyer AF. Background of U.S. hybrid corn. II. Crop Sci. 2004;44:370–80.

    Article  Google Scholar 

  13. Teng WT, Cao JS, Chen YH, Liu XH, Jing XQ, Zhang FJ, et al. Analysis of maize heterotic groups and patterns during past decade in China. Sci Agric Sin. 2004;37:1804–11.

    Google Scholar 

  14. Lu H, Bernardo R. Molecular marker diversity among current and historical maize inbreeds. Theor Appl Genet. 2001;103:613–7.

    Article  CAS  Google Scholar 

  15. Smith CW, Betrán J, Runge E. Corn: origin, history, technology, and production. Hoboken: Wiley; 2004.

    Google Scholar 

  16. Mikel MA, Dudley JW. Evolution of north American dent corn from public to proprietary germplasm. Crop Sci. 2006;46:1193–205.

    Article  Google Scholar 

  17. Li Y, Wang TY. Germplasm base of maize breeding in China and formation of foundation parents. J Maize Sci. 2010;18:1–8.

    Google Scholar 

  18. Zhang S, Li X, Yuan L, Li M, Peng Z. Heterotic groups and exploitation of heterosis–methodology, strategy, and use in hybrid maize breeding in China. In: Srinivasan G, Zaidi PH, Prasanna BM, Gonzalez F, Lesnick K, editors. Proceedings of the Eight Asian Regional Maize Workshop: New technologies for the new millennium. Bangkok; 2002. p. 64–9.

  19. Shull GH. The composition of a field of maize. American Breeders Assoc Rep. 1908;4:296–301.

    Google Scholar 

  20. Shull GH. What is heterosis? Genetics. 1948;33:439–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Davenport CB. Degeneration, albinism and inbreeding. Science. 1908;28:454–5.

    Article  CAS  Google Scholar 

  22. Bruce AB. The Mendelian theory of heredity and the augmentation of vigor. Science. 1910;32:627–8.

    Article  CAS  Google Scholar 

  23. East EM. Inbreeding in corn. Rep Conn Agric Exp Stn. 1907;1908:419–28.

    Google Scholar 

  24. Powers L. An expansion of Jones’s theory for the explanation of heterosis. Am Nat. 1944;78:275–80.

    Article  Google Scholar 

  25. Williams W. Heterosis and the genetics of complex characters. Nature. 1959;184:527–30.

    Article  CAS  Google Scholar 

  26. Stuber CW, Lincoln SE, Wolff DW, Helentjaris T, Lander ES. Identification of genetic factors contributing to heterosis in a hybrid from two elite maize inbred lines using molecular markers. Genetics. 1992;132:823–39.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Cockerham CC, Zeng ZB. Design III with marker loci. Genetics. 1996;143:1437–56.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Lu H, Romero-Severson J, Bernardo R. Genetic basis of heterosis explored by simple sequence repeat markers in a random-mated maize population. Theor Appl Genet. 2003;107:494–502.

    Article  CAS  Google Scholar 

  29. Frascaroli E, MA C’, Landi P, Pea G, Gianfranceschi L, Villa M, et al. Classical genetic and quantitative trait loci analyses of heterosis in a maize hybrid between two elite inbred lines. Genetics. 2007;176:625–44.

    Article  CAS  Google Scholar 

  30. Garcia AAF, Wang SC, Melchinger AE, Zeng ZB. Quantitative trait loci mapping and the genetic basis of heterosis in maize and rice. Genetics. 2008;180:1707–24.

    Article  Google Scholar 

  31. Tang JH, Yan JB, Ma XQ, Teng WT, Wu WR, Dai JR, et al. Dissection of the genetic basis of heterosis in an elite maize hybrid by QTL mapping in an immortalized F2 population. Theor Appl Genet. 2010;120:333–40.

    Article  Google Scholar 

  32. Schön CC, Dhillon BS, Utz HF, Melchinger AE. High congruency of QTL positions for heterosis of grain yield in three crosses of maize. Send to Theor Appl Genet. 2010;120(2):321–32.

    Article  Google Scholar 

  33. Larièpe A, Mangin B, Jasson S, Combes V, Dumas F, Jamin P, et al. The genetic basis of heterosis: multiparental quantitative trait loci mapping reveals contrasted levels of apparent overdominance among traits of agronomical interest in maize (Zea mays L.). Genetics. 2012;190:795–835.

    Article  Google Scholar 

  34. Guo T, Yang N, Tong H, Pan Q, Yang X, Tang J, et al. Genetic basis of grain yield heterosis in an “immortalized F2” maize population. Theor Appl Genet. 2014;127:2149.

    Article  Google Scholar 

  35. Wang H, Zhang X, Yang H, Liu X, Li H, Yuan L, et al. Identification of heterotic loci associated with grain yield and its components using two CSSL test populations in maize. Sci Rep. 2016;6:38205.

    Article  CAS  Google Scholar 

  36. Samayoa LF, Malvar RA, Butrón A. QTL for maize Midparent Heterosis in the heterotic pattern American dent × European Flint under corn borer pressure. Front Plant Sci. 2017;8:573.

    Article  Google Scholar 

  37. Li H, Yang Q, Gao L, Zhang M, Ni Z, Zhang Y. Identification of heterosis-associated stable QTLs for ear-weight-related traits in an elite maize hybrid Zhengdan 958 by design III. Front Plant Sci. 2017;8:561.

    PubMed  PubMed Central  Google Scholar 

  38. Wang Y, Zhang X, Shi X, Sun C, Jin J, Tian R, et al. Heterotic loci identified for maize kernel traits in two chromosome segment substitution line test populations. Sci Rep. 2018;8:11101.

    Article  Google Scholar 

  39. Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. Essex: Longman; 1996.

    Google Scholar 

  40. Hua JP, Xing YZ, Wei WR, Xu CG, Sun XL, Yu SB, et al. Single-locus heterotic effects and dominance by dominance interactions can adequately explain the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci U S A. 2003;100:2574–9.

    Article  CAS  Google Scholar 

  41. Zhou G, Chen Y, Yao W, Zhang C, Xie W, Hua J, et al. Genetic composition of yield heterosis in an elite rice hybrid. Proc Natl Acad Sci U S A. 2012;109:15847–52.

    Article  CAS  Google Scholar 

  42. Liu Y, Hou X, Xiao Q, Yi Q, Bian S, Hu Y, et al. Genetic analysis in maize foundation parents with mapping population and testcross population: Ye478 carried more favorable alleles and using QTL information could improve foundation parents. Front Plant Sci. 2016;7:1417.

    PubMed  PubMed Central  Google Scholar 

  43. Flint-Garcia SA, Buckler ES, Tiffin P, Ersoz E, Springer NM. Heterosis is prevalent for multiple traits in diverse maize germplasm. PLoS One. 2009;4(10):e7433.

    Article  Google Scholar 

  44. Melchinger AE, Utz HF, Piepho HP, Zeng ZB, Schön CC. The role of epistasis in the manifestation of heterosis: a systems-oriented approach. Genetics. 2007;177(3):1815–25.

    Article  CAS  Google Scholar 

  45. Yu S, Li J, Xu C, Tan Y, Gao Y, Li X, et al. Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci U S A. 1997;94:9226–31.

    Article  CAS  Google Scholar 

  46. Huang X, Yang S, Gong J, Zhao Y, Feng Q, Gong H, et al. Genomic analysis of hybrid rice varieties reveals numerous superior alleles that contribute to heterosis. Nat Commun. 2015;6:6258.

    Article  CAS  Google Scholar 

  47. Gao YM, Zhu J. Mapping QTLs with digenic epistasis under multiple environments and predicting heterosis based on QTL effects. Theor Appl Genet. 2007;115:325.

    Article  CAS  Google Scholar 

  48. Ma X, Tang JH, Teng WT, Yan JB, Meng YJ, Li JS. Epistatic interaction is an important genetic basis of grain yield and its components in maize. Mol Breed. 2007;20:41–51.

    Article  Google Scholar 

  49. Guo JJ, Chen ZL, Liu ZP, Wang BB, Song WB, Li W, et al. Identification of genetic factors affecting plant density response through QTL mapping of yield component traits in maize (Zea mays L.). Euphytica. 2011;182:409.

    Article  Google Scholar 

  50. Peng B, Li Y, Wang Y, Liu C, Liu Z, Tan W, et al. QTL analysis for yield components and kernel-related traits in maize across multi-environments. Theor Appl Genet. 2011;122:1305–20.

    Article  Google Scholar 

  51. Yang C, Tang D, Qu J, Zhang L, Zhang L, Chen Z, et al. Genetic mapping of QTL for the sizes of eight consecutive leaves below the tassel in maize (Zea mays L.). Theor Appl Genet. 2016;129(11):2191–209.

    Article  CAS  Google Scholar 

  52. Jannink JL. Identifying quantitative trait locus by genetic background interactions in association studies. Genetics. 2007;176:553–61.

    Article  CAS  Google Scholar 

  53. Chen L, An Y, Li Y, Li C, Shi Y, Song Y, et al. Candidate loci for yield-related traits in maize revealed by a combination of MetaQTL analysis and regional association mapping. Front Plant Sci. 2017;8:2190.

    Article  Google Scholar 

  54. Austin DF, Lee M. Comparative mapping in F2:3 and F6:7 generations of quantitative trait loci for grain yield and yield components in maize. Theor Appl Genet. 1996;92:817–26.

    Article  CAS  Google Scholar 

  55. Li YL, Niu SZ, Dong YB, Cui DQ, Wang YZ, Liu YY, et al. Identification of trait-improving quantitative trait loci for grain yield components from a dent corn inbred line in an advanced backcross BC2F2 population and comparison with its F2:3 population in popcorn. Theor Appl Genet. 2007;115:129–40.

    Article  CAS  Google Scholar 

  56. Li JZ, Zhang ZW, Li YL, Wang QL, Zhou YG. QTL consistency and meta-analysis for grain yield components in three generations in maize. Theor Appl Genet. 2011;122:771–82.

    Article  CAS  Google Scholar 

  57. Hua JP, Xing YZ, Xu CG, Sun XL, Yu SB, Zhang QF. Genetic dissection of an elite rice hybrid revealed that heterozygotes are not always advantageous for performance. Genetics. 2002;162:1885–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Pan Q, Xu Y, Li K, Peng Y, Zhan W, Li W, et al. The genetic basis of plant architecture in 10 maize recombinant inbred line populations. Plant Physiol. 2017;175:858–73.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. R Development Core Team. R: a language and environment for statistical computing. 2010.

    Google Scholar 

  60. Hallauer AR, Mirando FJB. Quantitative genetics in maize breeding. 2nd ed. Ames: Iowa State University; 1988.

    Google Scholar 

  61. Knapp S, Stroup W, Ross W. Exact confidence intervals for heritability on a progeny mean basis. Crop Sci. 1985;25:192–4.

    Article  Google Scholar 

  62. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.

    Article  Google Scholar 

  63. Chen DH, Ronald PC. A rapid DNA minipreparation method suitable for AFLP and other PCR applications. Plant Mol Biol Report. 1999;17:53–7.

    Article  CAS  Google Scholar 

  64. Fan JB, Gunderson KL, Bibikova M, Yeakley JM, Chen J, Wickham-Garcia E, et al. Illumina universal bead arrays. Methods Enzymol. 2006;410:57–73.

    Article  CAS  Google Scholar 

  65. Mathias L. MapDisto: fast and efficient computation of genetic linkage maps. Mol Breed. 2012;30:1231–5.

    Article  Google Scholar 

  66. Kosambi DD. The estimation of map distances from recombination values. Ann Eugenics. 1943;12:172–5.

    Article  Google Scholar 

  67. Meng L, Li H, Zhang L, Wang J. QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015;3(3):269–83.

    Article  Google Scholar 

  68. Li H, Ye G, Wang J. A modified algorithm for the improvement of composite interval mapping. Genetics. 2007;175:361–74.

    Article  Google Scholar 

  69. Li H, Ribaut JM, Li Z, Wang J. Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations. Theor Appl Genet. 2008;116:243–60.

    Article  Google Scholar 

  70. Wang D, Zhu J, Li Z, Paterson A. Mapping QTLs with epistatic effects and QTL × environment interactions by mixed linear model approaches. Theor Appl Genet. 1999;99:1255–64.

    Article  Google Scholar 

  71. Yang J, Hu C, Hu H, Yu R, Xia Z, Ye X, et al. QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations. Bioinformatics. 2008;24:721–3.

    Article  Google Scholar 

  72. Yang J, Zhu J, Williams RW. Mapping the genetic architecture of complex traits in experimental populations. Bioinformatics. 2007;23:1527–36.

    Article  CAS  Google Scholar 

  73. Doerge RW, Churchill GA. Permutation tests for multiple loci affecting a quantitative character. Genetics. 1996;142:285e294.

    Google Scholar 

  74. Rice WR. Analyzing tables of statistical tests. Evolution. 1989;43:223–5.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Dr. Ana Butrón and Dr. Pedro Revilla of Misión Biológica de Galicia (CSIC) for valuable suggestions and careful corrections, and also thank the help of Dr. Zhengjie Chen of Maize Research Institute for the construction of the IF2 population, and all students for participating in the fieldwork.

Funding

The current research was supported by the Project of National Major Basic Dairy Research “973” Plan (#2014CB138202 and #2011CB100106) and Science and Technology Plan Projects in Sichuan Province (#2016JY0065). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

YBH, YHL, and XBH designed the study. YHL and XBH developed the RILs, and QY and HL constructed the IF2 population. All coauthors conducted the experiments of phenotypic identification in the field. YHL, XBH, and QY analyzed the data. YHL and YQ prepared the manuscript, and all authors read and approved the final manuscript.

Corresponding author

Correspondence to Yubi Huang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. The performance of parental inbreds of the RIL population (YE478 and 08–641) and the hybrid F1 across four environments. The Duncan multiple range test was used for the comparison of means. Genotypes with lower-case letters were significantly different at the 0.05 probability level. Means with same letters are not significantly different. The data are shown as means ± SD (standard deviation). EWPE, ear weight per ear; CWPE, cob weight per ear; EL, ear length; ED, ear diameter; CD, cob diameter; RN, row number; KNPR, kernel number per ear; KWPE, kernel weight per ear; RKP, rate of kernel production. (JPG 513 kb)

Additional file 2:

Figure S2. Phenotypic variance explained (PVE) by QTL and epistasis interaction (EPI, digenic interaction) effects for nine yield-related traits in the RILs, IF2 population and the MPH dataset. The total PVE for QTL and epistasis was estimated on ICIM-ADD and ICIM-EPI, respectively. ICIM, inclusive composite interval mapping. RIL, recombinant inbred line. IF2, immortalized F2. MPH, mid-parent heterosis. EWPE, ear weight per ear; CWPE, cob weight per ear; EL, ear length; ED, ear diameter; CD, cob diameter; RN, row number; KNPR, kernel number per ear; KWPE, kernel weight per ear; RKP, rate of kernel production. (PPTX 37 kb)

Additional file 3:

Table S1. List of abbreviations and definitions for the yield-related traits recorded at harvest in a sample of ten ears per plot. (DOCX 16 kb)

Additional file 4:

Table S2. The performance of yield-related traits in the RILs, the IF2 population, and the MPH dataset across all environments. EWPE, ear weight per ear; CWPE, cob weight per ear; ED, ear diameter; CD, cob diameter; EL, ear length; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per row; RKP, rate of kernel production. a F1, The cross of 08–641 × YE478. RIL, the recombinant inbred lines. b IF2 and MPH indicate the immortalized F2 population and mid-parent heterosis, respectively. c Mean ± SE represents mean across all environments ± standard error. d MPH (%) = (F1-(P1 + P2)/2)/((P1 + P2)/2)*100%. (DOCX 18 kb)

Additional file 5:

Table S3. Pearson’s correlation coefficients among yield-related traits in the RILsa across four environments, the IF2c populations across three environments, and the mid-parent heterosis (MPH) dataset. * Significant at the 0.05 probability level. ** Significant at the 0.01 probability level. The absolute values of Pearson’s correlation coefficients (r2) with bold are equal to or larger than 0.5. EWPE, ear weight per ear; CWPE, cob weight per ear; ED, ear diameter; CD, cob diameter; EL, ear length; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per row; RKP, rate of kernel production. a Pearson’s correlation coefficients (r2) in the first row for each trait are for the RILs. b Pearson’s correlation coefficients (r2) in the second row for each trait are for the IF2 population. c Pearson’s correlation coefficients (r2) in the third row for each trait are for mid-parent heterosis. (DOCX 1433 kb)

Additional file 6:

Table S4. Main feature of QTLs for yield-related traits and for mid-parental heterosis for those traits from the RILs and IF2 population via single environment analysis. a The nomination of QTL is made as follows: a “q” standing for the abbreviation of QTL, the abbreviation of the trait, an “S” standing for single environment analyses, a number standing for chromosome, and another for physical order. EWPE, ear weight per ear; CWPE, cob weight per ear; ED, ear diameter; CD, cob diameter; EL, ear length; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per row; RKP, rate of kernel production. b RIL, recombinant inbred lines; IF2, the immortalized F2; MPH, mid-parental heterosis derived from the RIL and IF2 populations. c Env, environment; E1, 2014Jinghong; E2, 2015Jinghong; E3, 2016Chongzhou; E4, 2016Jinghong. d Chr, chromosome. e A, additive effect; Negative additive values indicate that the allele for increasing trait value is contributed by the parent 08–641; Positive additive values indicate that the allele for increasing trait value is contributed by another parent Ye478. f D, dominance effect. g PVE, phenotypic variance explained by QTL. (DOCX 70 kb)

Additional file 7:

Table S5. Pearson’s correlation coefficients between phenotypic performance per se in IF2 population and mid-parent heterosis for those traits. *** Significant at p < 0.0001. (DOCX 15 kb)

Additional file 8:

Table S6. Pearson’s correlation coefficients between phenotypic variation and genome marker heterozygosity. a Correlation between performance per se in IF2 population and percentage of heterozygous loci. b Correlation between mid-parent heterosis and the percentage of heterozygous loci. ***, significant at p < 0.0001 level; **, significant at p < 0.01 level; *, significant at p < 0.05 level; NS, not significant. EWPE, ear weight per ear; CWPE, cob weight per ear; ED, ear diameter; CD, cob diameter; EL, ear length; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per row; RKP, rate of kernel production. (DOCX 16 kb)

Additional file 9:

Table S7. Epistatic interactions for yield-related traits and mid-parent heterosis detected in the RILs and IF2 population across all environments. a EWPE, ear weight per ear; CWPE, cob weight per ear; CD, cob diameter; EL, ear length; RN, row number; KNPR, kernel number per row; KWPE, kernel weight per row; RKP, rate of kernel production. b Interval of the first QTL i. c Chromosome ID at the first scanning position. d A_i Additive effect of the QTL i. e D_i Dominance effect of the QTL i. f Interval of the second QTL j. g Chromosome ID at the second scanning position. h A_j Additive effect of the QTL j. i D_j Dominance effect of the QTL j. j AA: Estimated additive by the additive effect of QTLs at the two scanning positions. k AD: Estimated additive by dominance effect of QTLs at the two scanning positions. m DA: Estimated dominance by the additive effect of QTLs at the two scanning positions. n DD: Estimated dominance by dominance effect of QTLs at the two scanning positions. l h2(%): Heritability for epistatic QTL effects. The interval in bold was co-located in QTLs associated with the same trait via combined analysis across all environments. *, **, *** indicate significance at p < 0.05, p < 0.01, and p < 0.0001, respectively. (DOCX 26 kb)

Additional file 10:

Table S8. QTL × environment interactions for yield-related traits detected in the RILs under four environments and the IF2 population under three environments. a RIL, recombinant inbred lines; IF2, the immortalized F2. b AE indicates the additive by designed environment interaction effect. E1, E2, E3, and E4, represent 2014JH, 2015JH, 2016CZ, and 2016JH, respectively. c h2(AE, %) is the contribution rate of additive by environment interaction effect. d DE indicates the dominance by designed environment interaction effect. E1, E2, E3, and E4, represent 2014JH, 2015JH, 2016CZ, and 2016JH, respectively. e h2(DE, %) is the contribution rate of dominance by environment interaction effect. *, **, *** indicate significance at p < 0.05, p < 0.01, and p < 0.0001, respectively. The interval in bold was co-located in QTLs associated with the same trait via combined analysis across all environments. (DOCX 21 kb)

Additional file 11:

Genotypic and phenotypic data of the RIL and IF populations for the QTL analyses. (XLSX 1583 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yi, Q., Liu, Y., Hou, X. et al. Genetic dissection of yield-related traits and mid-parent heterosis for those traits in maize (Zea mays L.). BMC Plant Biol 19, 392 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-2009-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-2009-2

Keywords