Original Research ARTICLE
QTL for Maize Midparent Heterosis in the Heterotic Pattern American Dent × European Flint under Corn Borer Pressure
- 1Misión Biológica de Galicia, CSIC, Pontevedra, Spain
Despite the importance of heterosis and the efforts to comprehend this phenomenon, its molecular bases are still unknown. In this study, we intended to detect Quantitative trait loci (QTL) for mid-parent heterosis under infestation with the Mediterranean corn borer (MCB, Sesamia nonagrioides Lef.) using a North Carolina design III approach with a RIL population derived from a European flint inbred (EP42) × American dent inbred (A637) cross. QTL for heterosis of kernel yield have been positioned in regions corresponding to previously identified QTL for the same trait in different backgrounds. These results reinforce the high congruency of genes controlling heterosis across populations, even when populations have been developed from different heterotic patterns. A high percentage of genetic variation for mid-parent heterosis (Z2) for kernel yield could not be explained. Furthermore, genomic regions involved in heterosis for yield and plant height were not found despite the high genetic correlation between Z2 transformations for kernel yield and plant height. The moderate power in detecting QTL for mid-parent heterosis suggests that many genes with low augmented dominance effects contribute to the genetic architecture of mid-parent heterosis; dominance and additive-additive epistatic effects could also contribute to heterosis. However, results from this and previous studies suggest that the region 8.03–8.05 deserves special attention in future works in order to fine map loci involved in mid-parent heterosis for yield.
Heterosis is defined as the greater biomass, speed of development or fertility in heterozygotes, polyploids or hybrids compared to their parents (Schnable and Springer, 2013). Despite the importance of this phenomenom, its molecular basis is still unknown. In maize, mid-parent heterosis for yield, defined as the superiority of the hybrid over the mean of the parental inbreds, is particularly important because maize hybrids are extensively grown around the world. Most studies on the basis of maize mid-parent heterosis have been specifically focused in a unique heterotic pattern, Stiff Stalk Synthetic × Lancaster Sure Crop (Stuber et al., 1992; Frascaroli et al., 2007; Franco Garcia et al., 2008; Springer et al., 2009; Schon et al., 2010; Barber et al., 2012; Ding et al., 2014; Paschold et al., 2014), and conclusions from those studies cannot be generalized to other hybrids. Therefore, identifying quantitative trait loci (QTL) involved in heterosis in different genetic backgrounds is recommended as it has been rarely done (Tang et al., 2010; Lariepe et al., 2012; Guo et al., 2014). The heterotic pattern American dent × European Flint is particularly interesting because is the preferred pattern for hybrid development in Northern and Central Europe and Japan (Messmer et al., 1992; Moreno-Gonzalez et al., 2000; Technow et al., 2014).
The North Carolina design III and the immortilized F2 population design have been used to study the QTL involved in maize heterosis (Stuber et al., 1992; Cockerham and Zeng, 1996; Frascaroli et al., 2007; Tang et al., 2010; Guo et al., 2014). To detect QTL for mid-parent heterosis with an immortalized F2 population, a population of random crosses among RILs along with the RILs per se should be tested. Therefore, although the North Carolina design III with RILs does not allow to separate the contribution of dominance and epistasis to heterosis (Melchinger et al., 2007); it has been more extensively used because it does not require to test the RILs per se, but the backcrosses of those RILs to both parents (Frascaroli et al., 2007; Schon et al., 2010; He et al., 2012; Lariepe et al., 2012). The original North Carolina III mating design was developed by Comstock and Robinson (1948) to estimate the average level of dominance. This mating design has primarily been used in F2 populations to determine the effect of linkages on the estimates of additive and dominance variances, and the average level of dominance using progenies developed by backcrossing individual S0 plants of the F2 population to both parents (Hallauer and Miranda, 1981). In the present study, we attempted to detect QTL for mid-parent heterosis using a North Carolina design III approach with a RIL population derived from a cross between a European flint (EP42) and an American dent (A637) inbreds. The evaluations were conducted under infestation with the Mediterranean corn borer (MCB, Sesamia nonagrioides Lef.) to induce stalk damage, a major limiting factor for maize yield in the Mediterranean region. The cross EP42 × A637 presented high specific combining ability for yield under MCB infestation (Butron et al., 1999) and the backcrosses of RILs derived from this cross are ideal materials to study quantitative genetics of mid-parent heterosis for yield under infestation with MCB. Genes involved in heterosis under high pest pressure could be clearly different from those implicated in heterosis at optimal conditions because experimental data suggests that a suppression of defense response gene activities are important for generating the hybrid vigor phenotype due to the possible competition for resource distribution between plant defense and plant growth (Groszmann et al., 2015).
Materials and Methods
Plant Material and Genotypic Data
A population of 146 RILs derived from the cross of the European flint inbred line EP42 and the American dent inbred line A637 was developed (Samayoa et al., 2014) and 136 RILs were successfully backcrossed to both parents for mapping QTL for heterosis under infestation with MCB eggs. The mechanisms of resistance can be classified into three groups: antixenosis, antibiosis and tolerance (Painter, 1951). Inbreds EP42 and A637 were both classified as susceptible under MCB attack attending to their antibiosis levels, but A637 was tolerant (Butron et al., 1998) and showed good general combining ability for yield (Butron et al., 1999); while the cross EP42 × A637 presented significant and favorable specific combining abilities for yield and yield loss under MCB infestation (Butron et al., 1999). The linkage map used has been previously published by Samayoa et al. (2014).
The backcrosses of the 136 RILs to both parents were evaluated at Pontevedra (42°24′ N, 8°38′ W, and 20 m above sea level), Spain, in 2010 and 2011. Backcrosses of RILs were arranged in 17 sets comprised each one of 8 RILs backcrossed to both parents, backcrosses of the same RIL were kept in adjacent plots. Sets and RILs within sets were randomized in each replication (two per trial). The trials were hand planted and each experimental plot consisted of one row spaced 0.8 m apart from the other row with 13 two-kernel hills spaced 0.18 m apart. Plots were overplanted and thinned, obtaining a final density of ~70,000 plant ha−1. The evaluations were performed under artificial infestation with eggs of MCB obtained at the Misión Biologica de Galicia by rearing MCB as described by Eizaguirre and Albajes (1992) and Khan and Saxena (1997). Five plants of each plot were infested with ~40 MCB eggs placed between the stem and the sheath of a basal leaf. Data collected were: days to silking computed as the number of days from planting to 50% of plants in the plot showing silks; ear and plant heights on five representative plants as the length (in cm) from the ground to the main ear and from the ground to the top of the plant, respectively; stem lodging defined as the percentage of plants in the plot with the stem broken below the main ear; kernel resistance to MCB attack on the ears of the infested plants collected at harvest according to a subjective visual resistance scale from 1 to 9 (1 indicates completely damaged and 9 indicates no damage); tunnel length as the average stem tunnel length made by borers on the five infested plants; kernel yield of the plot expressed in Mg ha−1 at 140 g H2O kg−1 (infested and non-infested plants were considered); and kernel moisture at harvest as g of water in 100 g of kernels.
The individual analysis of each trial was made with the SAS mixed model procedure (PROC MIXED) (SAS Institute Inc, 2011) considering replications as random effects and sets and RILs within sets as fixed effects. The combined analysis across years was made considering years and replications as random effects. A Best linear unbiased estimate (BLUE) was obtained to estimate mean phenotypic values of each RIL backcrossed to P1 (EP42) and P2 (A637) in each year and across years.
To express the mid-parent heterosis as the sum of individual QTL effects, Melchinger et al. (2007) defined a new type of heterotic genetic effect that includes the dominance effect of QTLi minus half the sum of its additive × additive epistatic interactions and was named augmented dominance effect. The same authors defined another genetic effect (named augmented additive effect) that includes the additive effect for QTLi minus half the sum of dominance × additive epistatic interactions with all other QTL and corresponds exactly to the net contribution of QTLi to parental differences. They also demonstrate that contrasts Z1 and Z2 in a modification of design III of Comstock and Robinson (Comstock and Robinson, 1952) using RILs instead of a F2 account for augmented additive and dominance effects, respectively. The linear transformations Z1 and Z2 are defined as follow:
Z1 = (H1 + H2)/2 (trait mean across each pair of RIL backcrosses),
Z2 = (H1 – H2)/2 (half the trait difference between each pair of RIL backcrosses),
H1 and H2 being the phenotypic observations on backcrosses of each RIL with the parental inbreds P1 and P2, respectively.
As a genome scan with Z2 localizes genomic regions affecting mid-parent heterosis and with Z1 identifies QTL that contribute to parental differences, QTL analyses were performed with values of Z1 and Z2 for each RIL and trait across years. In addition, a joint QTL analysis of Z1 and Z2 was implemented. Genetic variances and heritabilities (ĥ2) across environments were estimated for Z1 and Z2 for each trait on a family-mean basis as described previously by Holland et al. (2003). The augmented degree of dominance (D*) was estimated as the square root of the division between genetic variances for Z2 and Z1. The genetic correlation coefficients among Z1 tranformations for the different trais, among Z2 tranformations and between Z1 and Z2 for each specific trait were computed following Holland (2006).
The QTL analyses were performed with the software package PlabMQTL (Utz, 2012). Composite interval mapping approach (CIM) was conducted for the QTL detection and to estimate QTL effects. According to a previously executed permutation test with 1000 random reshuffles, a LOD threshold of 2.5 was chosen to declare significant a putative QTL assuming an experiment-wise error <30%. The phenotypic variance explained by the QTL model was estimated by the adjusted coefficient of determination () which accounts for the total proportion of the phenotypic variance explained by all detected QTL in the final fit. Then, the phenotypic variance explained by an individual putative QTLi was calculated as:
Where n = total of QTL in the final fit and = partial coefficient of determination, estimated for the ith QTL detected (Zhu et al., 2004). The proportion of the genotypic variance () explained by all detected QTL in final fit was estimated from the ratio = /ĥ2 and the proportion of the genetic variance explained by each individual QTLi was estimated as = /ĥ2.
A 5-fold cross validation (CV) approach was employed for evaluating QTL mapping results (Utz et al., 2000). The data set (DS) was split at random into 5 data subsets; four subsets were combined to form de estimation set (ES) and the remaining subset formed the test set (TS). 1000 runs of CV were performed in order to determine the QTL frequency and bias of QTL effect estimate at the position of a QTL detected in the original data set (Melchinger et al., 2004). The proportion of the bias in the effect estimation of each individual QTL due to genotypic samplings was calculated as the difference between the average estimates in ES and TS divided by the estimate in ES.
Results and Discussion
Genetic variances for Z1 were significantly different from zero for all traits except for stalk lodging, tunnel length and kernel resistance; while genetic variances for Z2 were significantly different from zero for yield, yield-related traits, such as plant and ear heights, days to silking and stalk lodging (Table 1). The augmented degree of dominance was significantly higher than 1 for kernel yield, and significantly lower than 1 for days to silking, kernel moisture, and plant an ear heights. These results emphasize the great importance of heterosis for kernel yield and suggest that overdominance could play an important role as it has been reported in similar studies, although it should be taken into account that the augmented degree of dominance is biased by QTL linkage and epistasis and, consequently, does not exactly reflect gene action at individual loci (Melchinger et al., 2007; Schon et al., 2010; Lariepe et al., 2012). Heterosis was also important for yield-related traits, such as plant and ear heights and days to silking, but, as the augmented degree of dominance was significantly <1, dominance would be more important than over-dominance effects for these traits. In accordance to Lariepe et al. (2012), the higher importance of heterosis for days to silking in this cross compared to that in crosses between Reid and Lancaster inbreds could be due to the earliness of materials used and could emphasize the importance of studying heterosis in different backgrounds to get a better picture of the genetic architecture of the trait.
Table 1. Genetic variance ± standard error and heritability estimations for the transformations Z1 and Za2.
Heritabilities for Z1 of agronomical traits, such as days to silking, kernel moisture, plant and ear heights and kernel yield were low to moderate; while heritabilities for Z2 were higher for plant and ear heights and kernel yield, slightly lower for days to silking, and considerably lower for kernel moisture than heritability estimates for the corresponding Z1 transformations (Table 1). Heritabilities did not differ from zero for Z1 and Z2 for related-resistance traits, although in a previous work heritability for tunnel length of the RILs per se was low but significantly different from zero (Samayoa et al., 2014). The discrepancy between the heritabilities computed for tunnel length of the RILs per se and for Z1 in this study could be consequence of opposite effects for additive × dominance epistasis and additivity. It could also be due to a higher experimental error in the current experiment based on the number of genotypes tested (backcrosses of RILs to both parents) that doubled the number of RILs per se. As additive × additive epistasis effects were not significant for tunnel length in a previous study (Butron et al., 2009), the absence of heritability for Z2 for tunnel length would indicate that dominance effects are not important for resistance to stem tunneling by MCB as it has been already suggested (Butron et al., 1999; Cartea et al., 1999).
As the environment × genotype interaction was not significant for Z1 and Z2 for yield and resistance-related traits, QTL analysis was not independently done for each environment (Data not shown). No QTL for Z1 for kernel yield were found, while two were found for Z1 for days to silking that explained the 32% of the genetic variation, and one for Z1 for kernel moisture, tunnel length, plant and ear heights and kernel resistance that explained the 17, 72, 20, 21, and 18% of genetic variation, respectively (Table 2). The QTL identified for Z1 for tunnel length at 9.03–9.04 co-localized with QTL for Z1 for both plant and ear heights. There was not a complete agreement between QTL detected for Z1 and QTL detected for the per se traits in a previous study (Samayoa et al., 2014); although evaluations of the RILs per se and backcrossed to the parental inbreds were performed in adjacent trials. Higher experimental errors expected in the current experiment compared to those in which the RILs per se were evaluated could be responsible for detecting fewer QTL for Z1. In addition, additive and additive × additive epistasis effects were confounded in the study by Samayoa et al. (2014); while additive and additive × dominance epistasis effects are undistinguished in the current study and both contribute to the augmented additive (a*) effect. Therefore, epistasis effects in this particular cross could be important enough to counteract additive effects. Previous studies have (Yan et al., 2006; Guo et al., 2014) already pointed out that digenic interactions at the two-locus level might play an important role in the genetic basis of maize heterosis. In any case, QTL for Z1 for tunnel length, and plant and ear heights positioned at 9.03–9.04, QTL for Z1 for kernel resistance positioned at 5.05–5.06, and QTL for Z1 for days to silking positioned at 8.05–8.06 were placed in the same bin or in adjacent bins to those where QTL for these traits per se were previously found (Samayoa et al., 2014) corroborating the presence of genes with additive effects for these traits in those regions.
Table 2. Summary of QTL for transformations Z1 and Za2 of several traits evaluated under infestation with MCB in the backcrosses of a RIL population derived from EP42 xA637.
Quantitative trait loci (QTL) for mid-parent heterosis (Z2 transformation) were found for the agronomical traits stem lodging, kernel moisture, kernel yield, plant height, and ear height and explained the 42, 100, 39, 11, and 48% of genetic variance, respectively (Table 2). QTLs for mid-parent heterosis for kernel yield were located at bins 1.02–1.03, 304–3.05, 6.06, and 8.05, were stable across years, and explained the 8.5, 6.2, 8.9, and 15.3 of the genetic variance, respectively (Tables 2, 3). The accuracy of effect estimation was low for QTL at 1.02–1.03 and 3.04–3.05, moderate for the QTL at 6.06 and high for the QTL at 8.05 because the estimations of bias proportions derived from cross-validation runs were 0.56, 0.58, 0.30, and 0.08, respectively. In addition, the level of occurrence of these QTLs in cross-validations was moderate and high for QTL at 6.06 and 8.05, respectively. When the Akaike information criterion (AIC), a less conservative criterion, is used instead of the modified Bayesian information criterion (BIC) for fitting a model that includes all QTLs for heterosis for kernel yield, the percentage of genetic variance explained increases to 52% and additional QTL were found at 2.04–2.06, 3.03, and 5.03–5.04 (Data not shown). Schon et al. (2010) also detected QTLs for the Z2 transformation for kernel yield in all these regions (same or adjacent bins) in at least one of three populations derived from Reid × Lancaster inbred crosses. However, Lariepe et al. (2012) detected nine QTL for Z2 for kernel yield that explained the 53% of phenotypic variance in crosses among an European flint and two dent inbreds, but only three of the genomic regions detected (the centromeric regions of chromosomes 3 and 5 and the long arm of chromosome 8) bear QTL for Z2 in the current study. Lastly, Tang et al. (2010) located QTL for mid-parent heterosis for yield in bins 1.01, 1.10, and 8.03 using an immortalized F2 population derived from the cross between a Chinese dent inbred and an exotic flint inbred (Tang et al., 2010; Guo et al., 2014). Then, QTL for kernel yield heterosis under high insect pressure were located in genomic regions known for including loci involved in heterosis at optimal conditions, but these results cannot serve to declare that, in general, competition between hybrid vigor and defense mechanisms does not exist because, in the current study, genetic variability for resistance-related traits was low. The number of QTLs detected for yield heterosis were lower than those detected by Schon et al. (2010) in dent crosses, similar to those detected by Lariepe et al. (2012) in early populations derived from European flint × dent and dent × dent crosses and larger than those detected by Tang et al. (2010) in a dent × flint cross. Assuming similar population sizes, the higher power to detect QTL for mid-parent heterosis for yield in the Reid × Lancaster heterotic pattern than in other heterotic patterns evaluated would relate to the larger magnitude of heterosis (absolute value) in the Reid × Lancaster pattern compared to the others.
Table 3. Stability of QTLs for Z1 and Z2 transformationa across years: effects and percentage of phenotypic variability explained by each QTL in each year.
Results from this and previous studies suggest that the region 8.03–8.05 deserves special attention in future works in order to fine map loci involved in mid-parent heterosis for yield. This genomic region harbors QTL with important additive augmented (a*) effects for days to silking and dominance augmented (d*) effects for ear height and kernel yield in this study. Previously, in the same region, significant QTL for Z2 for yield were identified in several populations (Schon et al., 2010; Tang et al., 2010; Lariepe et al., 2012) and is one of the genomic regions with high density of QTL peaks for yield based on a MetaQTL analysis involving results from 44 manuscripts (Martinez et al., 2016). In addition, Samayoa et al. (2014) also detected important QTL for days to silking and plant and ear heights in the same region. This region presents high homology with the region of chromosome 6 [these regions correspond to duplicated regions of an ancient tetraploid (Schnable et al., 2011)] where a QTL for yield heterosis has been located in the current study and in that performed by Schon et al. (2010).
No overlapping QTL for Z1 and Z2 for kernel yield were found although genetic correlation between both parameters was moderate (0.50) suggesting that genomic regions containing QTL with augmented additive and dominance effects for yield went undetected due to the low magnitude of genetic effects (Table 4). In a similar way, although the correlation coefficients between Z1 for yield and Z1 for ear height and between Z2 for yield and Z2 for plant and ear heights were high and significant (Table 4), only two QTL for Z2 for kernel yield (at 8.05 and 5.03–5.04) co-localized with QTL for Z2 for ear height. These facts along with the lack of power for detecting QTL for Z1 for kernel yield and Z2 for days to silking, although genetic variability for these transformations significantly differed from zero (Table 1), suggest that many genes with low augmented additive and/or dominance effects can contribute to kernel yield and mid-parent heterosis for yield, at least in this particular cross. In addition, the augmented dominance effects of some loci that contribute to heterosis could be consequence of additive-additive epistatic effects. Lariepe et al. (2012) showed that when a QTL for Z2 for yield overlap with a QTL for yield with no significant dominance effect, the QTL probably presents significant additive-by-additive interaction effects. Therefore, the QTL detected for Z2 for yield at 5.03–5.04 that was already reported as important for yield per se (Samayoa et al., 2014) could rather contribute to Z2 due to additive-by-additive interactions than to dominance effects. The significant contribution of additive-by-additive interactions to heterosis in this particular cross was already reported in a previous study (Butron et al., 2009).
Table 4. Genetic correlation coefficients ± standard errors between Z1 and Z2 for yield and the corresponding Z1 and Z2 transformationsa for other traits.
As conclusions, the low to moderate power of the experiment to detect QTL for mid-parent heterosis suggests that many genes with low augmented dominance effects can contribute to the genetic architecture of mid-parent heterosis. Furthermore, dominance and additive-additive epistatic effects could also contribute to heterosis. The genomic regions detected for mid-parent heterosis for kernel yield in this particular cross (European flint × American dent inbreds) had been detected as significant for heterosis in other heterotic patterns reinforcing the idea of high congruency of heterosis genes across populations. Results from this and previous studies suggest that the region 8.03–8.05 deserves special attention in future works in order to fine map loci involved in mid-parent heterosis for yield.
RM and AB designed the work; LS did the collection and statistical analysis of data; and AB drafted the manuscript that was critically revised by LS and RM.
This work was supported by the National Plan for Research and Development of Spain (projects AGL2012-33415 and AGL2015-67313-C2-1-R, both funded in part by the European Regional Development Fund).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
LS acknowledges a contract JAE-Predoc from the Spanish Council for Scientific Research (CSIC). We thank Jeffrey Dunne for his careful language review and comments.
MCB, Mediterranean corn borer; QTL, Quantitative trait loci; RIL, Recombinant inbred line; CV, Cross validation; DS, Data set; ES, Estimation set; TS, Test set.
Barber, W. T., Zhang, W., Win, H., Varala, K. K., Dorweiler, J. E., Hudson, M. E., et al. (2012). Repeat associated small RNAs vary among parents and following hybridization in maize. Proc. Natl. Acad. Sci. U.S.A. 109, 10444–10449. doi: 10.1073/pnas.1202073109
Butron, A., Malvar, R. A., Velasco, P., Vales, M. I., and Ordas, A. (1999). Combining abilities for maize stem antibiosis, yield loss, and yield under infestation and non-infestation with pink stem borer. Crop Sci. 39, 691–696. doi: 10.2135/cropsci1999.0011183X003900020015x
Butron, A., Sandoya, G., Revilla, P., and Malvar, R. A. (2009). Genetics of resistance to the pink stem borer (Sesamia nonagrioides) in maize (Zea mays). Ann. Appl. Biol. 154, 205–217. doi: 10.1111/j.1744-7348.2008.00284.x
Cartea, M. E., Malvar, R. A., Butron, A., Vales, M. I., and Ordas, A. (1999). Inheritance of antibiosis to Sesamia nonagrioides (Lepidoptera: Noctuidae) in maize. J. Econ. Entomol. 92, 994–998. doi: 10.1093/jee/92.4.994
Comstock, R. E., and Robinson, H. F. (1948). The components of genetic variance in populations of biparental progenies and their use in estimating the average degree of dominance. Biometrics 4, 254–266. doi: 10.2307/3001412
Comstock, R. E., and Robinson, H. F. (1952). “Estimation of average dominance of genes,” in Heterosis: A Reord of Researches Directed Toward Explaining and Utilizing the Vigor of Hybrids, ed J. W. Gowen (Ames, IA: Iowa State College Press), 494–516.
Ding, H., Qin, C., Luo, X., Li, L., Chen, Z., Liu, H., et al. (2014). Heterosis in early maize ear inflorescence development: a genome-wide transcription analysis for two maize inbred lines and their hybrid. Int. J. Mol. Sci. 15, 13892–13915. doi: 10.3390/ijms150813892
Franco Garcia, A. A., Wang, S., Melchinger, A. E., and Zeng, Z.-B. (2008). Quantitative trait loci mapping and the genetic basis of heterosis in maize and rice. Genetics 180, 1707–1724. doi: 10.1534/genetics.107.082867
Frascaroli, E., Cane, M. A., Landi, P., Pea, G., Gianfranceschi, L., Villa, M., et al. (2007). Classical genetic and quantitative trait loci analyses of heterosis in a maize hybrid between two elite inbred lines. Genetics 176, 625–644. doi: 10.1534/genetics.106.064493
Groszmann, M., Gonzalez-Bayon, R., Lyons, R. L., Greaves, I. K., Kazan, K., Peacock, W. J., et al. (2015). Hormone-regulated defense and stress response networks contribute to heterosis in Arabidopsis F1 hybrids. Proc. Natl. Acad. Sci. U.S.A. 112, E6397–E6406. doi: 10.1073/pnas.1519926112
Guo, T., Yang, N., Tong, H., Pan, Q., Yang, X., Tang, J., et al. (2014). Genetic basis of grain yield heterosis in an immortalized F-2 maize population. Theor. Appl. Genet. 127, 2149–2158. doi: 10.1007/s00122-014-2368-x
Holland, J. B. (2006). Estimating genotypic correlations and their standard errors using multivariate restricted maximum likelihood estimation with SAS Proc MIXED. Crop Sci. 46, 642–654. doi: 10.2135/cropsci2005.0191
Khan, Z. R., and Saxena, R. C. (1997). Use of surrogate stem for eliciting ovipositional response of Busseola fusca (Lepidoptera: Noctuidae). J. Econ. Entomol. 90, 1426–1429. doi: 10.1093/jee/90.5.1426
Lariepe, A., Mangin, B., Jasson, S., Combes, V., Dumas, F., Jamin, P., et al. (2012). 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 190, U795–U835. doi: 10.1534/genetics.111.133447
Martinez, A. K., Soriano, J. M., Tuberosa, R., Koumproglou, R., Jahrmann, T., and Salvi, S. (2016). Yield QTLome distribution correlates with gene density in maize. Plant Sci. 242, 300–309. doi: 10.1016/j.plantsci.2015.09.022
Melchinger, A. E., Utz, H. F., Piepho, H. P., Zeng, Z. B., and Schon, C. C. (2007). The role of epistasis in the manifestation of heterosis: a systems-oriented approach. Genetics 177, 1815–1825. doi: 10.1534/genetics.107.077537
Melchinger, A. E., Utz, H. F., and Schön, C. C. (2004). QTL analyses of complex traits with cross validation, bootstrapping and other biometric methods. Euphytica 137, 1–11. doi: 10.1023/B:EUPH.0000040498.48379.68
Messmer, M. M., Melchinger, A. E., Boppenmaier, J., Brunklausjung, E., and Herrmann, R. G. (1992). Relationships among early European maize inbreds: I. Genetic diversity among flint and dent lines revealed by RFLPS. Crop Sci. 32, 1301–1309. doi: 10.2135/cropsci1992.0011183X003200060001x
Moreno-Gonzalez, J., Martinez, I., Brichette, I., Lopez, A., and Castro, P. (2000). Breeding potential of European flint and US corn belt dent maize populations for forage use. Crop Sci. 40, 1588–1595. doi: 10.2135/cropsci2000.4061588x
Paschold, A., Larson, N. B., Marcon, C., Schnable, J. C., Yeh, C.-T., Lanz, C., et al. (2014). Nonsyntenic genes drive highly dynamic complementation of gene expression in maize hybrids. Plant Cell 26, 3939–3948. doi: 10.1105/tpc.114.130948
Schnable, J. C., Springer, N. M., and Freeling, M. (2011). Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss. Proc. Natl. Acad. Sci. U.S.A. 108, 4069–4074. doi: 10.1073/pnas.1101368108
Schon, C. C., Dhillon, B. S., Utz, H. F., and Melchinger, A. E. (2010). High congruency of QTL positions for heterosis of grain yield in three crosses of maize. Theor. Appl. Genet. 120, 321–332. doi: 10.1007/s00122-009-1209-9
Springer, N. M., Ying, K., Fu, Y., Ji, T., Yeh, C.-T., Jia, Y., et al. (2009). Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content. PLoS Genet. 5:734. doi: 10.1371/journal.pgen.1000734
Stuber, C. W., Lincoln, S. E., Wolff, D. W., Helentjaris, T., and Lander, E. S. (1992). Identification of genetic-factors contributing to heterosis in a hybrid from 2 elite maize inbred lines using molecular markers. Genetics 132, 823–839.
Tang, J., Yan, J., Ma, X., Teng, W., Wu, W., Dai, J., et al. (2010). Dissection of the genetic basis of heterosis in an elite maize hybrid by QTL mapping in an immortalized F-2 population. Theor. Appl. Genet. 120, 333–340. doi: 10.1007/s00122-009-1213-0
Technow, F., Schrag, T. A., Schipprack, W., Bauer, E., Simianer, H., and Melchinger, A. E. (2014). Genome properties and prospects of genomic prediction of hybrid performance in a breeding program of maize. Genetics 197, U1343–U1469. doi: 10.1534/genetics.114.165860
Utz, H. F. (2012). PlabMQTL- Software for Meta-QTL Analysis with Composite Interval Mapping. version 0.9. Stuttgart: Institute of plant breeding, seed science, and population genetics, University of Hohenheim.
Utz, H. F., Melchinger, A. E., and Schön, C. C. (2000). Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics 154, 1839–1849.
Yan, J. B., Tang, H., Huang, Y. Q., Zheng, Y. L., and Li, J. S. (2006). Quantitative trait loci mapping and epistatic analysis for grain yield and yield components using molecular markers with an elite maize hybrid. Euphytica 149, 121–131. doi: 10.1007/s10681-005-9060-9
Keywords: Zea mays, yield, heterosis, quantitative trait loci, QTL
Citation: Samayoa LF, Malvar RA and Butrón A (2017) QTL for Maize Midparent Heterosis in the Heterotic Pattern American Dent × European Flint under Corn Borer Pressure. Front. Plant Sci. 8:573. doi: 10.3389/fpls.2017.00573
Received: 07 November 2016; Accepted: 29 March 2017;
Published: 19 April 2017.
Edited by:Ann E. Stapleton, University of North Carolina at Wilmington, USA
Reviewed by:Abhishek Bohra, Indian Institute of Pulses Research, India
Antoni Rafalski, DuPont (United States), USA
Copyright © 2017 Samayoa, Malvar and Butrón. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana Butrón, email@example.com
†Present Address: Luis F. Samayoa, Department of Crop and Soil Sciences, North Carolina State University, Raleigh, NC, USA