Identification of Heterosis-Associated Stable QTLs for Ear-Weight-Related Traits in an Elite Maize Hybrid Zhengdan 958 by Design III

Heterosis plays a decisive role in maize production worldwide, but its genetic basis remains unclear. In this study, we explored heterosis for ear-weight (EW)-related traits using a North Carolina Experiment III design (Design III) population derived from the elite maize hybrid Zhengdan 958. Quantitative trait loci (QTL) analysis was conducted based on phenotypic data collected from five environments using a high-density linkage map that consisted of 905 single nucleotide polymorphisms (SNP). A total of 38 environmentally stable QTLs were detected, and the numbers for the Z1 and Z2 populations were 18 and 20, respectively. All environmentally stable QTLs for Z2 were characterized by the overdominance effect (OD), which indicated that overdominance was one of the most important contributors to the heterosis of EW-related traits. Consistent with the significant positive correlations between EW-related traits, 9 genomic regions with overlapped QTLs for different traits were found and were located on chromosomes 1 (1), 3 (2), 4 (3), 7 (1), 8 (1), and 9 (1). Compared to previous reports, we found that the genomic regions for heterosis were not always congruent between different hybrids, which suggested that the combination of heterotic loci in different hybrids was genotype-dependent. Collectively, these data provided further evidence that the potential utilization of QTLs for heterosis may be feasible by pyramiding if we treat the QTLs as inherited units.


INTRODUCTION
Heterosis or hybrid vigor, which was denoted as the high-ranking performance of the F 1 hybrid relative to their parents, has played a decisive role in crop production for almost a century (East, 1908;Shull, 1908). However, constant efforts from different researchers have not led to consistent conclusions about the genetic mechanism of heterosis (Stuber, 1994;Hua et al., 2003;Radoev et al., 2008;Liang et al., 2015), which has hindered its use in crop improvement programs. To date, three main hypotheses have been proposed to explain the genetic basis of heterosis, including dominance and overdominance based on individual gene loci as well as epistasis based on the interactions among the genes. The dominance hypothesis emphasized the complementation of slightly deleterious recessive alleles that lie in the inbred parents (Jones, 1917). The overdominance hypothesis indicated that the performance of heterozygosity was superior to that of the homozygous condition (Shull, 1908;Hull, 1945). Epistasis attributed heterosis to interactions of genes from different loci (Richey, 1942;Powers, 1944).
Over the past two decades, quantitative genetics and the advance of molecular markers have provided powerful tools for the genetic analysis of heterosis in crops (Edwards et al., 1992). For example, Hua et al. (2003) found 33 single heterotic loci (HL) and that dominance by dominance interactions could adequately explain the genetic basis of heterosis in an elite rice hybrid using an "immortalized F 2. " In contrast, in a reanalysis of their dataset with an ultra-high density SNP bin map, Zhou et al. (2012) concluded that overdominance/pseudo-overdominance was very important to heterosis of yield and that dominance by dominance interactions is important for heterosis of grain weight and tillers per plant in rice. Recently, Huang et al. (2016) found support for the partial dominance of heterozygous locus for yieldrelated traits and better-parent heterosis for overall performance by conducting a GWAS analysis of 10,074 F 2 lines from 17 representative hybrid rice crosses. In rapeseed, epistasis along with all levels of dominance (partial dominance, dominance, and overdominance) were declared to be responsible for the expression of heterosis (Radoev et al., 2008;Li et al., 2012).
To date, several QTLs (genes) related to heterosis were fine mapped and/or cloned in various species. Using near isogenic lines (NILs), He et al. (2006) narrowed the yield-improving QTL qGY2-1 down to a 102.9-kb region on rice chromosome 2 and found allelic expression variation in this gene cluster. Krieger et al. (2010) reported a gene named SINGLE FLOWER TRUSS (SFT) in tomato, which was the first example of a single overdominant gene for yield. Heterozygosity for loss-of-function alleles of SFT increases the yield by up to 60%. Xue et al. (2008) isolated a quantitative trait locus Ghd7 from an elite rice hybrid, and this locus encoded a CCT domain protein. Analysis of backcross (BC) populations revealed that the dominance effect (D) was observed for the number of grains per panicle, plant height, and heading date because heterozygotes were close to the performance of the higher parent. In palm, Singh et al. (2013) cloned a SHELL gene that is a homolog of the MADS-box gene SEEDSTICK. This gene can improve the yield of mesocarp oil via heterodimerization, which provided a genetic explanation for single gene heterosis.
Maize has been applied as a model species for a long time in exploring heterosis for its high heterosis degree (Jones, 1922;Duvick, 2001). Grain yield and yield-related traits in maize were analyzed in various experimental designs. Adopting Abbreviations: QTL, Quantitative trait locus; RIL, Recombinant inbred line; TC, Testcross; ERN, Ear row number; ED, Ear diameter; RSN, Number of seeds per row; EL, Ear length; HSW, One hundred seed weight; ESN, Ear seed number, ESW, Ear seed weight; EW, Ear weight; MPH, Mid-parent heterosis; D * , Average degree of dominance; AA, Additive by additive interaction effects; DD, Dominance by dominance interaction effects. Design III, Stuber et al. (1992) did pioneering work to clarify the heterosis of grain yields and plant height with the aid of molecular markers using separate backcrosses. Their results showed that overdominance (or pseudo-overdominance) plays an especially important role in the phenomenon of heterosis. However, Cockerham and Zeng (1996) obtained a different conclusion for an alternative method of data analysis using the same data mentioned above, which considered dominance along with epistatic effects between linked QTLs to be the principal elements of heterosis in maize. The dominance hypothesis of grain yield and its related traits were also supported by the results of Tang et al. (2010), using an "IF 2 " design in an elite maize hybrid. Jiang et al. (2015) also concluded that dominance and epistatic QTLs are important for the heterosis of kernel-shaperelated traits using a triple testcross population in maize. The ever-changing clues for heterosis in previous studies indicated that the genetic basis of heterosis is intricate and that further studies should be pursued.
The hybrid Zhengdan 958 is one of the most widely grown commercial varieties in China currently (∼500 million hectares between 2001 and 2015) for its high planting density, good stresstolerance, and high but stable yield (Li et al., 2009;Lai et al., 2010). To date, only one study based on F 2:3 families derived from Zhengdan 958 was reported (Guo et al., 2011), but no research focusing on the genetic basis of heterosis underling Zhengdan 958 was found in the literature. The main objectives of this study were (1) to assess the level of heterosis for EW and its components; (2) to detect the QTLs and evaluate the effects related to heterosis; (3) to examine the role of epistasis in heterosis; and (4) to compare heterotic QTLs with previous studies.

Plant Materials
The high-heterosis hybrid Zhengdan 958 has been planted all over China. Its parental inbreeds belong to two different heterotic groups: Zheng 58 belongs to the PA heterotic group, a subgroup of SS, and Chang 7-2 belongs to TSPT heterotic group, a subgroup of NSS (Bai et al., 2015;Ma et al., 2016). Following the testcross (TC) progeny production scheme (Comstock and Robinson, 1952), 174 RILs (F 7 ) derived from the hybrid Zhengdan 958 were used as pollen parents to cross the parental lines Zheng 58 [TC (Zheng 58)] and Chang 7-2 [TC (Chang 7-2)]. Due to a seed shortage in a few crosses, 162 RILs and their corresponding testcrosses (TCs) were employed.

Field Experiments
The two populations of TC progeny along with other materials (i.e., the parental lines, Zhengdan 958 and RILs) were fieldtested in 2012 and 2013 on the experimental farm of Jilin Academy of Agricultural Sciences (Jilin province, China) and the experimental farm of the Xinjiang Academy of Agricultural Sciences (Xinjiang province, China) and, in 2012, on the Agronomy Farm of Jinghai (Tianjin, China) ( Figure S1 and Table  S1). The TCs were planted in a randomized complete block design with three replications at each location. Each plot included rows that were 4 m long with 0.67 m of space between rows. The population density was 45,000 plants per hectare. The 162 RILs, hybrid and the two parental lines of Zhengdan 958 were also planted using the same experimental design, which was near the TC experiment. All fields were well-watered with broad irrigation and rainfall. Other field management policies followed local standard practices.
Ten ears from consecutive plants in each plot were harvested and air-dried after maturity. Data for the following traits were collected: ear row number (ERN), ear diameter (ED), number of seeds per row (RSN), ear length (EL), one hundred seed weight (HSW), ear seed number (ESN), ear seed weight (ESW), and ear weight (EW).

Data Analysis
The percentage of heterosis was analyzed in the basic generations as mid-parent heterosis (MPH), which was computed as MPH = (F 1 -MP)/MP × 100, where MP represented the mid-parent value. Following the methods reported by Comstock and Robinson (1952) and Melchinger et al. (2007), the crosses of RILs to their parental lines Zheng 58 and Chang 7-2 were indicated as L 1i and L 2i (i = 1∼162), respectively. The linear transformations were Z 1i = (L 1i +L 2i )/2 and Z 2i = L 2i -L 1i . A combined ANOVA over five environments was calculated to estimate variance components. Additive Variances (V A ) within Z 1 and dominance variances (V D ) within Z 2 were used to estimate the average degree of dominance D * as (V D /2V A ) 0.5 , which stood for the degree of dominance over all separating loci (Cockerham and Zeng, 1996;Frascaroli et al., 2007;Melchinger et al., 2007).
The adjusted mean (Best Linear Unbiased Prediction, BLUP) values across five environments were calculated with the PROC MIXED procedure in SAS (SAS Institute Inc., North Carolina, USA). Broad-sense heritability (h 2 B ) were calculated as h 2 B = σ 2 g /(σ 2 g + σ 2 ge /n+σ 2 /nr), where σ 2 g is the genetic variance, σ 2 ge is the genotype by environment interaction variance, σ 2 is the error variance, n is the number of environments, and r is the number of replications of each experiment (Knapp et al., 1985;Churchill and Doerge, 1994). Correlation coefficients among traits were estimated using adjusted mean values for both Z 1 and Z 2 .

Genotyping and Linkage Analyses
The Zhengdan 958 RIL population along with the two parents were genotyped using a Maize SNP50 BeadChip (Ganal et al., 2011). Ten seeds from each genotype were germinated, and then young leaves were used for DNA extraction. The DNA quality was checked before genotyping. SNP genotyping was performed using GoldenGate assays (Illumina, San Diego, CA, USA) according to the manufacturer's protocol. SNP allele clustering and genotype calling was performed using Genome Studio v. 2011.1 software (Illumina). A genetic linkage map was constructed using MSTMap software (Wu et al., 2008).

QTL Analysis
For each Z s (s = 1, 2) population, the trait averaged values of three replicates for each environment were used for QTL analysis. The adjusted mean (BLUP) values for each trait across five environments were used for a combined analysis.
Composite interval mapping (CIM) implemented in Windows QTL Cartographer version 2.5 (Zeng, 1994;Wang et al., 2006) was employed. Model 6 from Zmapqtl was used, and 5 markers were set up as the number of cofactors. The QTLs were scanned with a 0.5-cM interval between markers. Permutation tests with a minimum of 1,000 replicates were adopted to determine the thresholds for the logarithm of odds (LOD) scores of putative QTL (Churchill and Doerge, 1994). The mapped QTL in Z 1 and Z 2 reflect the augmented additive effects a * i and augmented dominance effects d * i , separately (Melchinger et al., 2007). According to the scale of dominance degree commonly used in previous studies (Stuber et al., 1987;Jiang et al., 2015), the dominance degree ratios were estimated as |d * i /a * i | = augmented dominance effects/augmented additive effects: A, . QTL confidence intervals were determined based on positions ±2 LOD away from the peaks of the likelihood ratios (LRs) (Zhai et al., 2016). QTLs with overlapping confidence intervals were treated as congruent.
Based on the mixed model approach described by Wang et al. (1999), digenic epistasis QTL was analyzed using QTLMapper. The epistatic effects observed in Z 1 and Z 2 stood for additive by additive (AA) and dominance by dominance (DD) interaction effects, respectively.

Heterosis and Population Performance
The average field performances of eight EW-related traits for the hybrid Zhengdan 958 and its parental lines (Zheng 58 and Chang 7-2) for five environments are listed in Table 1. Zheng 58 had higher RSN and HSW compared to Chang 7-2, whereas ERN, ESN, and ED of Chang 7-2 were significantly higher than those of Zheng 58 (P < 0.01). Compared to parental lines, hybrid Zhengdan 958 exhibited overwhelming superiority in all eight traits. Notably, EW and ESW showed high MPH (138.61 and 152.87%), followed by ESN, RSN,and EL (67.21,56.25,and 39.60%), and ED and ERN were relatively smaller (22.09 and 16.80%).
Mean values and broad sense heritability (h 2 B ) of Z 1 and Z 2 for each trait are listed in Table 2. The h 2 B in the Z 1 population ranged from 0.80 to 0.96, among which ERN, RSN, EL, HSW, and ESN had high heritability (h 2 B > 0.90). For the Z 2 population, h 2 B varied from 0.71 to 0.91. Remarkably, the h 2 B of ESW and EW in Z 2 was higher than in Z 1, whereas the rest were lower compared to Z 1 .
Variance analysis of Z 1 and Z 2 showed that V A and V D for all traits were significant (P < 0.01) ( Table 2). Furthermore, we calculated the average degree of dominance (D * ) for each trait. The results revealed that the D * was > 1 for RSN, ESN, ESW, and EW and was < 1 for ERN, EL, and HSW.
Correlation coefficients among the eight traits within Z 1 and Z 2 are listed in Table 3. Significantly positive correlations were observed between EW and six traits for both Z 1 and Z 2 , including HSW, RSN, ED, EL, ESN, and ESW. Notably, HSW had positive correlations with ED and EL and ESW but negative correlations with ERN and ESN in Z 1 . In addition, EL was positively  correlated with ED and ERN in Z 2 but was negatively correlated in Z 1 .

Construction of the SNP-Based Genetic Linkage Map
In total, 905 SNP markers exhibiting polymorphisms between Zheng 58 and Chang 7-2 were adopted to construct the genetic linkage map (Appendix A). Of the 905 markers, 214 (23.6%) showed distortion segregation at P < 0.05, and 127 (14.0%) showed distortion segregation at P < 0.01. However, a previous study concluded that distorted markers will not just have a great effect on QTL detection (Zhang et al., 2010). As a result, these SNP markers were assigned to 10 chromosomes, spanning 2402.0 cM, with an average of 2.65 cM between adjacent markers ( Figure  S2 and Table S2).

Mapping Environmentally Stable QTLs for Z 1 and Z 2
Based on a genetic linkage map of 905 SNP markers, 483 QTLs were detected for eight traits for the Z 1 and Z 2 populations, which were distributed on all 10 chromosomes (Appendix B). We defined a QTL detected within two or more environments and in the combined analysis as an "environmentally stable QTL." According to this criterion, 38 environmentally stable QTLs were detected on chromosomes 1, 3, 4, 5, 7, 8, 9, and 10 ( Figure 1 and Table 4). Forty-five QTLs associated with ERN were detected. Five environmentally stable QTLs for ERN were identified on chromosomes 3, 7, and 10, and they were designated qERN.A-3.2, qERN.A-3.3, qERN.A-7.1, qERN.A-7.2, and qERN.A-10.1. Parental line Zheng 58 contributed an increased additive effect (A) for qERN.A-3.2 and qERN.A-3.3, which explained the 9.8 and 9.0% variation observed in the combined analysis. In contrast, parental line Chang 7-2 contributed increased additive effects (A) for qERN.A-7.1, qERN.A-7.2, and qERN.A-10.1, which explained 9.2, 7.0, and 10.7% of the variation in the combined analysis, respectively. No environmentally stable QTL was detected for Z 2 .
Sixty-four QTLs associated with ED were detected. Two environmentally stable QTLs were identified on chromosome 4 (qED.A-4.4, qED.B-4.5). Parental line Chang 7-2 contributed additive effects (A) for the increased ED of qED.A-4.4, which explained 8.2% of the ED variation for the combined analysis. qED.B-4.5 showed an overdominance effect (OD) for Z 2 , which explained 14.4% of the variation for the combined analysis.

Analysis of Digenic Interaction across the Entire Genome
Previous studies reported that epistasis also played a certain role in heterosis (Frascaroli et al., 2007;Jiang et al., 2015). In this study, we thus also analyzed digenic interaction. In total, 206 pairs of digenic epistatic QTLs were detected, and the numbers for Z 1 and Z 2 were 122 and 84, respectively (Appendix C). For each trait, both additive by additive (AA) and dominance by dominance digenic interactions were observed. The total variation explained by all digenic interactions either for Z 1 or Z 2 , was < 25% for most traits. The highest values of accumulated R 2 were found for HSW for Z 1 (41.13%, E3) and for ERN for Z 2 (35.97%, E1), respectively.
Interestingly, several digenic epistatic regions were mapped to the confidence intervals of QTLs (Appendix C). For example, the genetic region flanked by SNP18775 and SNP18742 on chromosome 4 was found to share an additive QTL for RSN via an AA interaction with the genetic region on chromosome 5. In FIGURE 1 | Continued black vertical bar with a red triangle represents additive QTLs with positive alleles from parent Zheng 58; a red vertical bar with a red triangle represents additive QTLs with positive alleles from parent Chang 7-2. Double-headed arrows represent the genomic regions characterized by QTL clusters. Red shadows on the physical map indicate the corresponding positions of each QTL. The verticals in different colors alongside the physical map indicate known heterotic-related QTLs from different studies Frascaroli et al., 2007;Tang et al., 2010;Guo et al., 2011Guo et al., , 2014.
addition, the genetic region flanked by SNP7302 and SNP7373 on chromosome 1 was found to interact with a genetic region flanked by SNP54321 and SNP25914 on chromosome 5 via a DD interaction, and both genetic regions shared QTL for HSW. However, it is necessary to note that most of digenic epistatic regions did not co-localize with the regions of main-effect QTLs.

The Role of Allelic Interactions in Heterosis
Heterosis has a revolutionary influence on the maize breeding program. There has been considerable interest in the genetic basis of heterosis. Previous studies have detected multiple QTLs related to maize heterosis based on various experimental designs. However, the number and QTL positions were not always congruent among different crosses Lu et al., 2003;Frascaroli et al., 2007;Tang et al., 2010;Larièpe et al., 2012;Guo et al., 2014;Jiang et al., 2015). Our present study was conducted based on a genetic linkage map of 905 SNP markers in five environments, which can precisely define the positions and stability of QTLs. Of the 483 QTLs detected, thirty-eight environmentally stable QTLs for eight traits (ERN, ED, RSN, EL, HSW, ESN, ESW, and EW) were selected. To be specific, 18 of them were detected for Z 1 , and 20 were for Z 2 . Notably, all 20 environmentally stable QTLs detected for Z 2 were characterized by positive OD. Interestingly, the heterosis level (MPH) and D * of different traits were globally met with the proportion of "overdominant" QTL, which revealed a good consistency of classical genetic analysis and QTL analysis. For example, the traits with the highest heterosis level and average degree of dominance, particularly EW and ESW, were the ones that showed the highest proportion of the OD. In contrast, traits with the lowest heterosis level and average degree of dominance, such as ERN, were the ones that showed the lowest proportion of the OD. These results may suggest the distinct genetic architectures of studied traits. The phenomenon of the overwhelming superiority of OD was reported by several research groups (Frascaroli et al., 2007;Larièpe et al., 2012). However, pseudo-overdominance cannot be excluded. For example, Graham et al. (1997) dissected an overdominant QTL on chromosome 5 associated with grain yield into two linked dominant QTLs by fine mapping. Li et al. (2015) showed that two separate loci with a repulsion linkage could appear as a single locus with an overdominance mode of inheritance. Therefore, to uncover whether the overdominance effects found are true, fine-mapping strategies should be used in future work.  Stuber et al. (1992); 2 Frascaroli et al. (2007); 3 Tang et al. (2010); 4 Guo et al. (2011);5 Guo et al. (2014).
Comparison analysis revealed that the heterosis related genomic regions in this study were also reported for yield and/or yield-related traits in different studies (Figure 1 and Table 4). For instance, genomic region 1.2 was detected to have pleiotropic effects for EW, ESW, ESN, EL, and RSN in the present study and was reported to affect grain yield and the number of kernels per plant by Stuber et al. (1992) and Frascaroli et al. (2007Frascaroli et al. ( , 2009 as well as by Larièpe et al. (2012) for grain yield. The genomic region 3.2 that we highlighted for ESN, ESW, and EW was determined to be a heterotic locus for plant height, number of kernels per plant, and grain yield by Frascaroli et al. (2007Frascaroli et al. ( , 2009. Region 8 appeared to be involved in EW, ESW, RSN, and EL, and Frascaroli et al. (2007) reported that it was a QTL cluster for grain yield, 100-seed weight, number of kernels per plant and plant height, in addition to a report by Tang et al. (2010) for grain yield. Remarkably, we found that genomic regions for heterosis of yield or yield-related traits were not always congruent between any two studies, which suggested that the combination of heterotic loci in different hybrids was genotype-dependent (Figure 1 and Table 4). Collectively, it can be seen that QTL analysis of different hybrids would broaden our understanding of the mechanism of heterosis. Moreover, these data provided further evidence for the notion that potential utilization of QTLs for heterosis is feasible by pyramiding if we treat them as inherited units, which deserves further investigation (Lu et al., 2003).

Epistasis Effect on Heterosis
Epistasis, which is an important genetic phenomenon of interactions between non-allelic genes, has been reported to exert certain roles in heterosis in maize (Frascaroli et al., 2007;Ma et al., 2007;Tang et al., 2010). In this study, the number of epistatic interactions in a single environment ranged from 0 to 5 for a given trait. The degree of epistatic interactions varied in different environments but at a relatively low level with R 2 < 25% for total epistatic interactions, either for AA or DD. The low level of epistasis effects was documented in the literature. For example, Ma et al. (2007) reported that 44 pairs of interactions for yield and its components were found and that the total contribution of the digenic interactions was 7%. Notably, only a few cases involved loci in epistasis interactions that were co-localized with main effect QTLs, as confirmed in previous studies (Frascaroli et al., 2007;Ma et al., 2007;Zhang et al., 2014). One of the causal reasons may result from the limitation of Design III in separating the QTL effect and their epistatic interactions with other QTLs in the analysis of heterosis (Melchinger et al., 2008). Additionally, no environmentally stable epistasis interaction was detected for any trait, which suggested that epistatic interactions were susceptible to environmental influence.

Pleiotropic of Environmentally Stable QTLs
Because multiplication effects are a main cause of heterosis for grain yield in maize, it was difficult to dissect the genetic components underlying this complex trait. In this study, we chose EW and its secondary component traits for detailed analysis because these traits were less complex and because the results may be more credible and intuitive (Williams, 1959;Hua et al., 2003). Our data revealed that significant positive correlations were observed between EW and EW-related traits (Table 3).
Theoretically, QTLs for EW and its components would be overlapped to some degree. As expected, a total of 10 genomic regions contained 33 co-localized environmentally stable QTLs, among which region 7.1 contained two tightly linked QTLs for same trait (ERN), whereas the other 9 genomic regions harbored QTLs for two or more different traits (Figure 1 and Table 4). For example, genomic region 1.2 affected five traits (EW, ESW, ESN, EL, and RSN) simultaneously, and all of them showed an overdominance effect OD. Region 8 harbored four QTLs for RSN, EL, ESW, and EW, with OD, and regions 3.2, 4.2, and 9 harbored overdominant QTLs for 3, 2, and 2 traits, respectively. Interestingly, region 4.1 contained two tightly linked overdominant QTLs for HSW and two additive QTLs for RSN and ED. In addition, additive QTLs harbored in genomic regions 3.1, 4.3, and 7.3 contributed to two EW-related traits, which agreed with the significant correlation between them. However, similar to the case of previous studies (Kusterer et al., 2007;Larièpe et al., 2012), we cannot determine that QTLs in the same genomic region were pleiotropic and/or closely linked QTLs for the reason of limited samples and relatively large confidence intervals for QTL positions.

The Usefulness of Design III in Dissecting the Genetic Basis of Heterosis
In this study, we adopted a variant of Design III from the single cross between maize inbred lines Zheng 58 and Chang 7-2 to analyze QTLs contributing to heterosis. The detection by the modified Design III can not only determine the precise locations of heterotic QTL but also estimate the augmented QTL effect a * i and d * i (Melchinger et al., 2007). The Design III population also possessed several advantages. First, the whole population could be recreated, which would allow experiments with replications under many different environments to be conducted, and this population could even be used for alternate experimental schemes. Second, similar to immortalized F 2 populations, the phenotypic values used to evaluate heterosis come from hybrids instead of progenies, which may lead to inbreeding depression (Hua et al., 2003). Third, as shown in previous studies, this approach offers opportunities for analyzing heterosis (Hua et al., 2003;Frascaroli et al., 2007;Jiang et al., 2015). Recently, Guo et al. (2011) evaluated a set of 231 F 2:3 families derived from the same hybrid Zhengdan 958 at two different plant densities to analyze the genetic basis of 12 yield-related traits. A comparison analysis revealed that two heterotic QTL regions were detected in both studies, including region 1.2 and 4.2 (Figure 1 and Table 4). For example, the ones on chromosome 1.08 for Ear weight, Grain weight per ear, Ear length, 100-kernel weight, and Kernel number per row was related to EW, ESW, ESN, EL, and RSN in our study. Nevertheless, the heterotic QTLs with an OD in regions 1.1, 3.2, 4.1, and 9 in our study were not reported by Guo et al. (2011). Collectively, our data showed that Design III was more powerful for analyzing heterotic QTLs than F 2:3 families.