Prediction of single-cross hybrid performance for grain yield and grain dry matter content in maize using AFLP markers associated with QTL

Optimizing selection based on BLUPs or BLUEs in multiple sets of genotypes differing in their population parameters

KEY MESSAGE

Selection response in truncation selection across multiple sets of candidates hinges on their post-selection proportions, which can deviate grossly from their initial proportions. For BLUPs, using a uniform threshold for all candidates maximizes the selection response, irrespective of differences in population parameters. Plant breeding programs typically involve multiple families from either the same or different populations, varying in means, genetic variances and prediction accuracy of BLUPs or BLUEs for true genetic values (TGVs) of candidates. We extend the classical breeder's equation for truncation selection from single to multiple sets of genotypes, indicating that the expected overall selection response ( Δ G Tot ) for TGVs depends on the selection response within individual sets and their post-selection proportions. For BLUEs, we show that maximizingΔ G Tot requires thresholds optimally tailored for each set, contingent on their population parameters. For BLUPs, we prove thatΔ G Tot is maximized by applying a uniform threshold across all candidates from all sets. We provide explicit formulas for the origin of the selected candidates from different sets and show that their proportions before and after selection can differ substantially, especially for sets with inferior properties and low proportion. We discuss implications of these results for (a) optimum allocation of resources to training and prediction sets and (b) the need to counteract narrowing the genetic variation under genomic selection. For genomic selection of hybrids based on BLUPs of GCA of their parent lines, selecting distinct proportions in the two parent populations can be advantageous, if these differ substantially in the variance and/or prediction accuracy of GCA. Our study sheds light on the complex interplay of selection thresholds and population parameters for the selection response in plant breeding programs, offering insights into the effective resource management and prudent application of genomic selection for improved crop development.

Portability of genomic predictions trained on sparse factorial designs across two maize silage breeding cycles

KEY MESSAGE

We validated the efficiency of genomic predictions calibrated on sparse factorial training sets to predict the next generation of hybrids and tested different strategies for updating predictions along generations. Genomic selection offers new prospects for revisiting hybrid breeding schemes by replacing extensive phenotyping of individuals with genomic predictions. Finding the ideal design for training genomic prediction models is still an open question. Previous studies have shown promising predictive abilities using sparse factorial instead of tester-based training sets to predict single-cross hybrids from the same generation. This study aims to further investigate the use of factorials and their optimization to predict line general combining abilities (GCAs) and hybrid values across breeding cycles. It relies on two breeding cycles of a maize reciprocal genomic selection scheme involving multiparental connected reciprocal populations from flint and dent complementary heterotic groups selected for silage performances. Selection based on genomic predictions trained on a factorial design resulted in a significant genetic gain for dry matter yield in the new generation. Results confirmed the efficiency of sparse factorial training sets to predict candidate line GCAs and hybrid values across breeding cycles. Compared to a previous study based on the first generation, the advantage of factorial over tester training sets appeared lower across generations. Updating factorial training sets by adding single-cross hybrids between selected lines from the previous generation or a random subset of hybrids from the new generation both improved predictive abilities. The CDmean criterion helped determine the set of single-crosses to phenotype to update the training set efficiently. Our results validated the efficiency of sparse factorial designs for calibrating hybrid genomic prediction experimentally and showed the benefit of updating it along generations.

Genetic resources and breeding of maize for Striga resistance: a review

The potential yield of maize (Zea mays L.) and other major crops is curtailed by several biotic, abiotic, and socio-economic constraints. Parasitic weeds, Striga spp., are major constraints to cereal and legume crop production in sub-Saharan Africa (SSA). Yield losses reaching 100% are reported in maize under severe Striga infestation. Breeding for Striga resistance has been shown to be the most economical, feasible, and sustainable approach for resource-poor farmers and for being environmentally friendly. Knowledge of the genetic and genomic resources and components of Striga resistance is vital to guide genetic analysis and precision breeding of maize varieties with desirable product profiles under Striga infestation. This review aims to present the genetic and genomic resources, research progress, and opportunities in the genetic analysis of Striga resistance and yield components in maize for breeding. The paper outlines the vital genetic resources of maize for Striga resistance, including landraces, wild relatives, mutants, and synthetic varieties, followed by breeding technologies and genomic resources. Integrating conventional breeding, mutation breeding, and genomic-assisted breeding [i.e., marker-assisted selection, quantitative trait loci (QTL) analysis, next-generation sequencing, and genome editing] will enhance genetic gains in Striga resistance breeding programs. This review may guide new variety designs for Striga-resistance and desirable product profiles in maize.

Predicting heterosis via genetic distance and the number of SNPs in selected segments of chromosomes in maize

A reliable method is needed for predicting heterosis to help maize (Zea mays L.) breeders develop new hybrids more efficiently. The objectives of this study were to 1) investigate if the numbers of selected PEUS SNPs (the SNP in the Promoters (1 kb upstream of the start codon), Exons, Untranslated region (UTR), and Stop codons) could be used for predicting MPH or BPH of GY; 2) if the number of PEUS SNPs is a better predictor of MPH and/or BPH of GY than genetic distance (GD). A line × tester experiment was conducted with 19 elite maize inbreds from three heterotic groups, which were crossed with five testers. The multi-location trial data on GY were recorded. Whole-genome resequencing of the 24 inbreds was carried out. After filtration, a total of 58,986,791 SNPs were called with high confidence. Selected SNPs in the promoters, exons, untranslated region (UTRs), and stop codons (PEUS SNPs) were counted, and the GD was calculated. The correlation between heterozygous PEUS SNPs/GD and mean MPH, BPH of GY revealed that 1) both the number of heterozygous PEUS SNP and the GD were highly correlated to both MPH_GY and BPH_GY at p<0.01 with correlation coefficients for the number of heterozygous PEUS SNP being higher than that for GD; 2) the mean number of heterozygous PEUS SNPs was also highly correlated with mean BPH_GY or mean MPH_GY (p<0.05) in the 95 crosses grouped by either male or female parents, implying that inbreds can be selected before making the actual crosses in the field. We concluded that the number of heterozygous PEUS SNPs would be a better predictor of MPH_GY and BPH_GY than GD. Hence, maize breeders could use heterozygous PEUS SNPs to select inbreds with high heterosis potential before actually making the crosses, thus improving the breeding efficiency.

Genomic prediction of hybrid performance: comparison of the efficiency of factorial and tester designs used as training sets in a multiparental connected reciprocal design for maize silage

Calibrating a genomic selection model on a sparse factorial design rather than on tester designs is advantageous for some traits, and equivalent for others. In maize breeding, the selection of the candidate inbred lines is based on topcross evaluations using a limited number of testers. Then, a subset of single-crosses between these selected lines is evaluated to identify the best hybrid combinations. Genomic selection enables the prediction of all possible single-crosses between candidate lines but raises the question of defining the best training set design. Previous simulation results have shown the potential of using a sparse factorial design instead of tester designs as the training set. To validate this result, a 363 hybrid factorial design was obtained by crossing 90 dent and flint inbred lines from six segregating families. Two tester designs were also obtained by crossing the same inbred lines to two testers of the opposite group. These designs were evaluated for silage in eight environments and used to predict independent performances of a 951 hybrid factorial design. At a same number of hybrids and lines, the factorial design was as efficient as the tester designs, and, for some traits, outperformed them. All available designs were used as both training and validation set to evaluate their efficiency. When the objective was to predict single-crosses between untested lines, we showed an advantage of increasing the number of lines involved in the training set, by (1) allocating each of them to a different tester for the tester design, or (2) reducing the number of hybrids per line for the factorial design. Our results confirm the potential of sparse factorial designs for genomic hybrid breeding.

Genomic prediction of hybrid crops allows disentangling dominance and epistasis

We revisited, in a genomic context, the theory of hybrid genetic evaluation models of hybrid crosses of pure lines, as the current practice is largely based on infinitesimal model assumptions. Expressions for covariances between hybrids due to additive substitution effects and dominance and epistatic deviations were analytically derived. Using dense markers in a GBLUP analysis, it is possible to split specific combining ability into dominance and across-groups epistatic deviations, and to split general combining ability (GCA) into within-line additive effects and within-line additive by additive (and higher order) epistatic deviations. We analyzed a publicly available maize data set of Dent × Flint hybrids using our new model (called GCA-model) up to additive by additive epistasis. To model higher order interactions within GCAs, we also fitted "residual genetic" line effects. Our new GCA-model was compared with another genomic model which assumes a uniquely defined effect of genes across origins. Most variation in hybrids is accounted by GCA. Variances due to dominance and epistasis have similar magnitudes. Models based on defining effects either differently or identically across heterotic groups resulted in similar predictive abilities for hybrids. The currently used model inflates the estimated additive genetic variance. This is not important for hybrid predictions but has consequences for the breeding scheme-e.g. overestimation of the genetic gain within heterotic group. Therefore, we recommend using GCA-model, which is appropriate for genomic prediction and variance component estimation in hybrid crops using genomic data, and whose results can be practically interpreted and used for breeding purposes.

Identifying quantitative trait loci for the general combining ability of yield-relevant traits in maize

Maize is the most important staple crop worldwide. Many of its agronomic traits present with a high level of heterosis. Combining ability was proposed to exploit the rule of heterosis, and general combining ability (GCA) is a crucial measure of parental performance. In this study, a recombinant inbred line population was used to construct testcross populations by crossing with four testers based on North Carolina design II. Six yield-relevant traits were investigated as phenotypic data. GCA effects were estimated for three scenarios based on the heterotic group and the number of tester lines. These estimates were then used to identify quantitative trait loci (QTL) and dissect genetic basis of GCA. A higher heritability of GCA was obtained for each trait. Thus, testing in early generation of breeding may effectively select candidate lines with relatively superior GCA performance. The GCA QTL detected in each scenario was slightly different according to the linkage mapping. Most of the GCA-relevant loci were simultaneously detected in all three datasets. Therefore, the genetic basis of GCA was nearly constant although discrepant inbred lines were appointed as testers. In addition, favorable alleles corresponding to GCA could be pyramided via marker-assisted selection and made available for maize hybrid breeding.

Genetic Dissection of the General Combining Ability of Yield-Related Traits in Maize

Maize yield components including row number, kernel number per row, kernel thickness, kernel width, kernel length, 100-kernel weight, and volume weight affect grain yield directly. Previous studies mainly focused on dissecting the genetic basis of per se performances for yield-related traits, but the genetic basis of general combining ability (GCA) for these traits is still unclear. In the present study, 328 RILs were crossed as males to two testers according to the NCII mating design, resulting in a hybrid panel composed of 656 hybrids. Both the hybrids and parental lines were evaluated in four environments in 2015 and 2016. Correlation analysis showed the performances of GCA effects were significantly correlated to the per se performances of RILs for all yield-related traits (0.17 ≤ r ≤ 0.64, P > 0.01). Only 17 of 95 QTL could be detected for both per se performances of RILs and GCA effects for eight yield-related traits. The QTL qKN7-1 and qHKW1-3, which could explain more than 10% of the variation in the GCA effects of KN and HKW, were also detected for per se performances for the traits. The pleiotropic loci qRN3-1 and qRN6, which together explained 14.92% of the observed variation in GCA effects for RN, were associated with the GCA effects of KW and HKW, but not with per se performances for these traits. In contrast, Incw1, which was related to seed weight in maize, was mapped to the region surrounding MK2567 at the qHKW5-2 locus, but no GCA effect was detected. The QTL identified in present study for per se performances and corresponding GCA effects for yield-related traits might be useful for maize hybrid breeding.

Revisiting hybrid breeding designs using genomic predictions: simulations highlight the superiority of incomplete factorials between segregating families over topcross designs

Simulations showed that hybrid performances issued from an incomplete factorial between segregating families of two heterotic groups enable to calibrate genomic predictions of hybrid value more efficiently than tester-based designs. Genomic selection offers new opportunities to revisit hybrid breeding by replacing extensive phenotyping of hybrid combinations by genomic predictions. A key question remains to identify the best design to calibrate genomic prediction models. We proposed to use single-cross hybrids issued from an incomplete factorial design between segregating populations and compared this strategy with a conventional approach based on topcross evaluation. Two multiparental segregating populations of lines, each specific of one heterotic group, were simulated. Hybrids considered as training sets were generated using either (1) a parental line from the opposite group as tester or (2) following an incomplete factorial design. Different specific combining ability (SCA) proportions were simulated by considering different levels of group divergence and dominance effects for the simulated QTL. For the incomplete factorial design, for a same number of hybrids, we considered different numbers of parental lines and different contributions of lines (one to four) to calibration hybrids. We evaluated for different training set sizes prediction accuracies of new hybrids and genetic gains along three generations. At a given training set size, factorial design was as efficient (considering accuracy) as tester design in additive scenarios, but significantly outperformed tester design when SCA was present. The contribution number of each parental line to the incomplete factorial design had a small impact on accuracies. Our simulations confirmed experimental results and showed that calibrating models on hybrids between two multiparental populations is a cost-efficient way to perform genomic predictions in both groups, opening prospects for revisiting reciprocal recurrent selection schemes.

Hybrid Breeding for MLN Resistance: Heterosis, Combining Ability, and Hybrid Prediction

Prior knowledge on heterosis and quantitative genetic parameters on maize lethal necrosis (MLN) can help the breeders to develop numerous resistant or tolerant hybrids with optimum resources. Our objectives were to (1) estimate the quantitative genetic parameters for MLN disease severity, (2) investigate the efficiency of the prediction of hybrid performance based on parental per se and general combining ability (GCA) effects, and (3) examine the potential of hybrid prediction for MLN resistance or tolerance based on markers. Fifty elite maize inbred lines were selected based on their response to MLN under artificial inoculation. Crosses were made in a half diallel mating design to produce 307 F1 hybrids. All hybrids were evaluated in MLN quarantine facility in Naivasha, Kenya for two seasons under artificial inoculation. All 50 inbreds were genotyped with genotyping-by-sequencing (GBS) SNPs. The phenotypic variation was significant for all traits and the heritability was moderate to high. We observed that hybrids were superior to the mean performance of the parents for disease severity (−14.57%) and area under disease progress curve (AUDPC) (14.9%). Correlations were significant and moderate between line per se and GCA; and mean of parental value with hybrid performance for both disease severity and AUDPC value. Very low and negative correlation was observed between parental lines marker based genetic distance and heterosis. Nevertheless, the correlation of GCA effects was very high with hybrid performance which can suggests as a good predictor of MLN resistance. Genomic prediction of hybrid performance for MLN is high for both traits. We therefore conclude that there is potential for prediction of hybrid performance for MLN. Overall, the estimated quantitative genetic parameters suggest that through targeted approach, it is possible to develop outstanding lines and hybrids for MLN resistance.

Identification of heterotic groups in South-Asian-bred hybrid parents of pearl millet

Pearl millet breeding programs can use this heterotic group information on seed and restorer parents to generate new series of pearl millet hybrids having higher yields than the existing hybrids. Five hundred and eighty hybrid parents, 320 R- and 260 B-lines, derived from 6 pearl millet breeding programs in India, genotyped following RAD-GBS (about 0.9 million SNPs) clustered into 12 R- and 7 B-line groups. With few exceptions, hybrid parents of all the breeding programs were found distributed across all the marker-based groups suggesting good diversity in these programs. Three hundred and twenty hybrids generated using 37 (22 R and 15 B) representative parents, evaluated for grain yield at four locations in India, showed significant differences in yield, heterosis, and combining ability. Across all the hybrids, mean mid- and better-parent heterosis for grain yield was 84.0% and 60.5%, respectively. Groups G12 B × G12 R and G10 B × G12 R had highest heterosis of about 10% over best check hybrid Pioneer 86M86. The parents involved in heterotic hybrids were mainly from the groups G4R, G10B, G12B, G12R, and G13B. Based on the heterotic performance and combining ability of groups, 2 B-line (HGB-1 and HGB-2) and 2 R-line (HGR-1 and HGR-2) heterotic groups were identified. Hybrids from HGB-1 × HGR-1 and HGB-2 × HGR-1 showed grain yield heterosis of 10.6 and 9.3%, respectively, over best hybrid check. Results indicated that parental groups can be formed first by molecular markers, which may not predict the best hybrid combination, but it can reveal a practical value of assigning existing and new hybrid pearl millet parental lines into heterotic groups to develop high-yielding hybrids from the different heterotic groups.

Quantitative Trait Locus Mapping of the Combining Ability for Yield-Related Traits in Wild Rice Oryza longistaminata

In decades of hybrid rice breeding, the combining ability has been successfully used to evaluate excellent parental lines and predict heterosis. However, previous studies for the combining ability mainly focused on cultivated rice and rarely involved wild rice. In this study, for the first time, we identified 20 new quantitative trait loci (QTLs) for the combining ability in wild rice using a North Carolina II mating design. Among them, qGCA1, one of the major QTLs that can significantly improve the general combining ability of the plant height, spikelet number, and yield per plant, was delimited to an interval of about 72 kb on chromosome 1. qSCA8, another major QTL, which can significantly improve the specific combining ability of the seed-setting rate and yield per plant, was located in an interval of about 90 kb on chromosome 8. These QTLs discovered from wild rice will provide new ideas to explain the genetic mechanism of the combining ability and establish the basis for breeding of high-combining-ability rice.

Across-years prediction of hybrid performance in maize using genomics

Inclusion of historical training data improved the genomics-based prediction of performance of maize hybrids, the extent depending on the phenotypic trait and genotype-by-year interaction. Prediction of hybrid performance using existing phenotypic data on previous hybrids combined with molecular data collected on the parent lines allows to identify the most promising candidates from a huge number of possible hybrids at an early stage. Phenotypic data on yield and dry matter of 1970 grain maize hybrids from 19 years of a public breeding program were aggregated considering the underlying structure of factorial sets of hybrids. Pedigree records and 50 K SNP data were collected on their 170 Dent and 127 Flint parent lines. The performance of untested hybrids was predicted by best linear unbiased predictors (BLUP) on basis of pedigree or genomic data. For composition of training sets (TRN) and test sets (TST), three schemes for collecting factorials from specific years were employed which resulted in 490 scenarios. For each scenario, the predictive ability and genomic relationship between TRN and TST hybrids were determined. For extended TRNs, where earlier years were successively added to the TRN, the maximum relationship increased and the predictive ability improved, with the extent of the latter depending on the phenotypic trait and its genotype-by-year interaction. Genomic BLUP outperformed pedigree BLUP and better utilized the early years' data, especially for prediction of hybrids from factorials in a more distant future. This study on hybrid prediction in grain maize illustrated that including historical phenotypic data for training, although consisting of less related genotypes, can improve genomic prediction and enables optimization of hybrid variety development.

Dissecting the Genetic Basis Underlying Combining Ability of Plant Height Related Traits in Maize

Maize plant height related traits including plant height, ear height, and internode number are tightly linked with biomass, planting density, and grain yield in the field. Previous studies have focused on understanding the genetic basis of plant architecture traits per se, but the genetic basis of combining ability remains poorly understood. In this study, 328 recombinant inbred lines were inter-group crossed with two testers to produce 656 hybrids using the North Carolina II mating design. Both of the parental lines and hybrids were evaluated in two summer maize-growing regions of China in 2015 and 2016. QTL mapping highlighted that 7 out of 16 QTL detected for RILs per se could be simultaneously detected for general combining ability (GCA) effects, suggesting that GCA effects and the traits were genetically controlled by different sets of loci. Among the 35 QTL identified for hybrid performance, 57.1% and 28.5% QTL overlapped with additive/GCA and non-additive/SCA effects, suggesting that the small percentage of hybrid variance due to SCA effects in our design. Two QTL hotspots, located on chromosomes 5 and 10 and including the qPH5-1 and qPH10 loci, were validated for plant height related traits by Ye478 derivatives. Notably, the qPH5-1 locus could simultaneously affect the RILs per se and GCA effects while the qPH10, a major QTL (PVE > 10%) with pleiotropic effects, only affected the GCA effects. These results provide evidence that more attention should be focused on loci that influence combining ability directly in maize hybrid breeding.

Small RNA-based prediction of hybrid performance in maize

BACKGROUND

Small RNA (sRNA) sequences are known to have a broad impact on gene regulation by various mechanisms. Their performance for the prediction of hybrid traits has not yet been analyzed. Our objective was to analyze the relation of parental sRNA expression with the performance of their hybrids, to develop a sRNA-based prediction approach, and to compare it to more common SNP and mRNA transcript based predictions using a factorial mating scheme of a maize hybrid breeding program.

RESULTS

Correlation of genomic differences and messenger RNA (mRNA) or sRNA expression differences between parental lines with hybrid performance of their hybrids revealed that sRNAs showed an inverse relationship in contrast to the other two data types. We associated differences for SNPs, mRNA and sRNA expression between parental inbred lines with the performance of their hybrid combinations and developed two prediction approaches using distance measures based on associated markers. Cross-validations revealed parental differences in sRNA expression to be strong predictors for hybrid performance for grain yield in maize, comparable to genomic and mRNA data. The integration of both positively and negatively associated markers in the prediction approaches enhanced the prediction accurary. The associated sRNAs belong predominantly to the canonical size classes of 22- and 24-nt that show specific genomic mapping characteristics.

CONCLUSION

Expression profiles of sRNA are a promising alternative to SNPs or mRNA expression profiles for hybrid prediction, especially for plant species without reference genome or transcriptome information. The characteristics of the sRNAs we identified suggest that association studies based on breeding populations facilitate the identification of sRNAs involved in hybrid performance.

Parental Expression Variation of Small RNAs Is Negatively Correlated with Grain Yield Heterosis in a Maize Breeding Population

Heterosis refers to a quantitative phenomenon in which F1 hybrid trait values exceed the mean of the parental values in a positive direction. Generally, it is dependent on a high degree of heterozygosity, which is maintained in hybrid breeding by developing parental lines in separate, genetically distinct heterotic groups. The mobility of small RNAs (sRNAs) that mediate epigenetic regulation of gene expression renders them promising candidates for modulating the action of combined diverse genomes in trans-and evidence already indicates their contribution to transgressive phenotypes. By sequencing small RNA libraries of a panel of 21 maize parental inbred lines we found a low overlap of 35% between the sRNA populations from both distinct heterotic groups. Surprisingly, in contrast to genetic or gene expression variation, parental sRNA expression variation is negatively correlated with grain yield (GY) heterosis. Among 0.595 million expressed sRNAs, we identified 9,767, predominantly 22- and 24-nt long sRNAs, which showed an association of their differential expression between parental lines and GY heterosis of the respective hybrids. Of these sRNAs, 3,485 or 6,282 showed an association with high or low GY heterosis, respectively, thus the low heterosis associated group prevailing at 64%. The heterosis associated sRNAs map more frequently to genes that show differential expression between parental lines than reference sets. Together these findings suggest that trans-chromosomal actions of sRNAs in hybrids might add up to a negative contribution in heterosis formation, mediated by unfavorable gene expression regulation. We further revealed an exclusive accumulation of 22-nt sRNAs that are associated with low GY heterosis in pericentromeric genomic regions. That recombinational suppression led to this enrichment is indicated by its close correlation with low recombination rates. The existence of this enrichment, which we hypothesize resulted from the separated breeding of inbred lines within heterotic groups, may have implications for hybrid breeding strategies addressing the recombinational constraints characteristic of complex crop genomes.

Reciprocal Genetics: Identifying QTL for General and Specific Combining Abilities in Hybrids Between Multiparental Populations from Two Maize (Zea mays L.) Heterotic Groups

Several plant and animal species of agricultural importance are commercialized as hybrids to take advantage of the heterosis phenomenon. Understanding the genetic architecture of hybrid performances is therefore of key importance. We developed two multiparental maize (Zea mays L.) populations, each corresponding to an important heterotic group (dent or flint) and comprised of six connected biparental segregating populations of inbred lines (802 and 822 lines for each group, respectively) issued from four founder lines. Instead of using "testers" to evaluate their hybrid values, segregating lines were crossed according to an incomplete factorial design to produce 951 dent-flint hybrids, evaluated for four biomass production traits in eight environments. QTL detection was carried out for the general-combining-ability (GCA) and specific-combining-ability (SCA) components of hybrid value, considering allelic effects transmitted from each founder line. In total, 42 QTL were detected across traits. We detected mostly QTL affecting GCA, 31% (41% for dry matter yield) of which also had mild effects on SCA. The small impact of dominant effects is consistent with the known differentiation between the dent and flint heterotic groups and the small percentage of hybrid variance due to SCA observed in our design (∼20% for the different traits). Furthermore, most (80%) of GCA QTL were segregating in only one of the two heterotic groups. Relative to tester-based designs, use of hybrids between two multiparental populations appears highly cost efficient to detect QTL in two heterotic groups simultaneously. This presents new prospects for selecting superior hybrid combinations with markers.

Genomic Prediction of Sunflower Hybrids Oil Content

Prediction of hybrid performance using incomplete factorial mating designs is widely used in breeding programs including different heterotic groups. Based on the general combining ability (GCA) of the parents, predictions are accurate only if the genetic variance resulting from the specific combining ability is small and both parents have phenotyped descendants. Genomic selection (GS) can predict performance using a model trained on both phenotyped and genotyped hybrids that do not necessarily include all hybrid parents. Therefore, GS could overcome the issue of unknown parent GCA. Here, we compared the accuracy of classical GCA-based and genomic predictions for oil content of sunflower seeds using several GS models. Our study involved 452 sunflower hybrids from an incomplete factorial design of 36 female and 36 male lines. Re-sequencing of parental lines allowed to identify 468,194 non-redundant SNPs and to infer the hybrid genotypes. Oil content was observed in a multi-environment trial (MET) over 3 years, leading to nine different environments. We compared GCA-based model to different GS models including female and male genomic kinships with the addition of the female-by-male interaction genomic kinship, the use of functional knowledge as SNPs in genes of oil metabolic pathways, and with epistasis modeling. When both parents have descendants in the training set, the predictive ability was high even for GCA-based prediction, with an average MET value of 0.782. GS performed slightly better (+0.2%). Neither the inclusion of the female-by-male interaction, nor functional knowledge of oil metabolism, nor epistasis modeling improved the GS accuracy. GS greatly improved predictive ability when one or both parents were untested in the training set, increasing GCA-based predictive ability by 10.4% from 0.575 to 0.635 in the MET. In this scenario, performing GS only considering SNPs in oil metabolic pathways did not improve whole genome GS prediction but increased GCA-based prediction ability by 6.4%. Our results show that GS is a major improvement to breeding efficiency compared to the classical GCA modeling when either one or both parents are not well-characterized. This finding could therefore accelerate breeding through reducing phenotyping efforts and more effectively targeting for the most promising crosses.

General and specific combining abilities in a maize (Zea mays L.) test-cross hybrid panel: relative importance of population structure and genetic divergence between parents

General and specific combining abilities of maize hybrids between 288 inbred lines and three tester lines were highly related to population structure and genetic distance inferred from SNP data. Many studies have attempted to provide reliable and quick methods to identify promising parental lines and combinations in hybrid breeding programs. Since the 1950s, maize germplasm has been organized into heterotic groups to facilitate the exploitation of heterosis. Molecular markers have proven efficient tools to address the organization of genetic diversity and the relationship between lines or populations. The aim of the present work was to investigate to what extent marker-based evaluations of population structure and genetic distance may account for general (GCA) and specific (SCA) combining ability components in a population composed of 800 inter and intra-heterotic group hybrids obtained by crossing 288 inbred lines and three testers. Our results illustrate a strong effect of groups identified by population structure analysis on both GCA and SCA components. Including genetic distance between parental lines of hybrids in the model leads to a significant decrease of SCA variance component and an increase in GCA variance component for all the traits. The latter suggests that this approach can be efficient to better estimate the potential combining ability of inbred lines when crossed with unrelated lines, and limits the consequences of tester choice. Significant residual GCA and SCA variance components of models taking into account structure and/or genetic distance highlight the variation available for breeding programs within structure groups.

Genomic Prediction of Sunflower Hybrids Oil Content

Prediction of hybrid performance using incomplete factorial mating designs is widely used in breeding programs including different heterotic groups. Based on the general combining ability (GCA) of the parents, predictions are accurate only if the genetic variance resulting from the specific combining ability is small and both parents have phenotyped descendants. Genomic selection (GS) can predict performance using a model trained on both phenotyped and genotyped hybrids that do not necessarily include all hybrid parents. Therefore, GS could overcome the issue of unknown parent GCA. Here, we compared the accuracy of classical GCA-based and genomic predictions for oil content of sunflower seeds using several GS models. Our study involved 452 sunflower hybrids from an incomplete factorial design of 36 female and 36 male lines. Re-sequencing of parental lines allowed to identify 468,194 non-redundant SNPs and to infer the hybrid genotypes. Oil content was observed in a multi-environment trial (MET) over 3 years, leading to nine different environments. We compared GCA-based model to different GS models including female and male genomic kinships with the addition of the female-by-male interaction genomic kinship, the use of functional knowledge as SNPs in genes of oil metabolic pathways, and with epistasis modeling. When both parents have descendants in the training set, the predictive ability was high even for GCA-based prediction, with an average MET value of 0.782. GS performed slightly better (+0.2%). Neither the inclusion of the female-by-male interaction, nor functional knowledge of oil metabolism, nor epistasis modeling improved the GS accuracy. GS greatly improved predictive ability when one or both parents were untested in the training set, increasing GCA-based predictive ability by 10.4% from 0.575 to 0.635 in the MET. In this scenario, performing GS only considering SNPs in oil metabolic pathways did not improve whole genome GS prediction but increased GCA-based prediction ability by 6.4%. Our results show that GS is a major improvement to breeding efficiency compared to the classical GCA modeling when either one or both parents are not well-characterized. This finding could therefore accelerate breeding through reducing phenotyping efforts and more effectively targeting for the most promising crosses.

Genomic Prediction of Single Crosses in the Early Stages of a Maize Hybrid Breeding Pipeline

Prediction of single-cross performance has been a major goal of plant breeders since the beginning of hybrid breeding. Recently, genomic prediction has shown to be a promising approach, but only limited studies have examined the accuracy of predicting single-cross performance. Moreover, no studies have examined the potential of predicting single crosses among random inbreds derived from a series of biparental families, which resembles the structure of germplasm comprising the initial stages of a hybrid maize breeding pipeline. The main objectives of this study were to evaluate the potential of genomic prediction for identifying superior single crosses early in the hybrid breeding pipeline and optimize its application. To accomplish these objectives, we designed and analyzed a novel population of single crosses representing the Iowa Stiff Stalk synthetic/non-Stiff Stalk heterotic pattern commonly used in the development of North American commercial maize hybrids. The performance of single crosses was predicted using parental combining ability and covariance among single crosses. Prediction accuracies were estimated using cross-validation and ranged from 0.28 to 0.77 for grain yield, 0.53 to 0.91 for plant height, and 0.49 to 0.94 for staygreen, depending on the number of tested parents of the single cross and genomic prediction method used. The genomic estimated general and specific combining abilities showed an advantage over genomic covariances among single crosses when one or both parents of the single cross were untested. Overall, our results suggest that genomic prediction of single crosses in the early stages of a hybrid breeding pipeline holds great potential to redesign hybrid breeding and increase its efficiency.

Prediction of hybrid performance in maize with a ridge regression model employed to DNA markers and mRNA transcription profiles

Background

Ridge regression models can be used for predicting heterosis and hybrid performance. Their application to mRNA transcription profiles has not yet been investigated. Our objective was to compare the prediction accuracy of models employing mRNA transcription profiles with that of models employing genome-wide markers using a data set of 98 maize hybrids from a breeding program.

Results

We predicted hybrid performance and mid-parent heterosis for grain yield and grain dry matter content and employed cross validation to assess the prediction accuracy. Prediction with a ridge regression model using random effects for mRNA transcription profiles resulted in similar prediction accuracies than employing the model to DNA markers. For hybrids, of which none of the parental inbred lines was part of the training set, the ridge regression model did not reach the prediction accuracy that was obtained with a model using transcriptome-based distances.

Conclusion

We conclude that mRNA transcription profiles are a promising alternative to DNA markers for hybrid prediction, but further studies with larger data sets are required to investigate the superiority of alternative prediction models.

Genomic Prediction of Testcross Performance in Canola (Brassica napus)

Genomic selection (GS) is a modern breeding approach where genome-wide single-nucleotide polymorphism (SNP) marker profiles are simultaneously used to estimate performance of untested genotypes. In this study, the potential of genomic selection methods to predict testcross performance for hybrid canola breeding was applied for various agronomic traits based on genome-wide marker profiles. A total of 475 genetically diverse spring-type canola pollinator lines were genotyped at 24,403 single-copy, genome-wide SNP loci. In parallel, the 950 F1 testcross combinations between the pollinators and two representative testers were evaluated for a number of important agronomic traits including seedling emergence, days to flowering, lodging, oil yield and seed yield along with essential seed quality characters including seed oil content and seed glucosinolate content. A ridge-regression best linear unbiased prediction (RR-BLUP) model was applied in combination with 500 cross-validations for each trait to predict testcross performance, both across the whole population as well as within individual subpopulations or clusters, based solely on SNP profiles. Subpopulations were determined using multidimensional scaling and K-means clustering. Genomic prediction accuracy across the whole population was highest for seed oil content (0.81) followed by oil yield (0.75) and lowest for seedling emergence (0.29). For seed yieId, seed glucosinolate, lodging resistance and days to onset of flowering (DTF), prediction accuracies were 0.45, 0.61, 0.39 and 0.56, respectively. Prediction accuracies could be increased for some traits by treating subpopulations separately; a strategy which only led to moderate improvements for some traits with low heritability, like seedling emergence. No useful or consistent increase in accuracy was obtained by inclusion of a population substructure covariate in the model. Testcross performance prediction using genome-wide SNP markers shows considerable potential for pre-selection of promising hybrid combinations prior to resource-intensive field testing over multiple locations and years.

OsPRR37 and Ghd7 are the major genes for general combining ability of DTH, PH and SPP in rice

Artificial selection of high yield crops and better livestock is paramount importance in breeding programs. Selection of elite parents with preferred traits from a phalanx of inbred lines is extremely laborious, time-consuming and highly random. General combining ability (GCA) was proposed and has been widely used for the evaluation of parents in hybrid breeding for more than half a century. However, the genetic and molecular basis of GCA has been largely overlooked. Here, we present two pleotropic QTLs are accounting for GCA of days to heading (DTH), plant height (PH) and spikelet per panicle (SPP) using an F₂-based NCII design, the BC₃F₂ population as well as a set of nearly isogenic lines (NILs) with five testers. Both GCA1 and GCA2 were loss-of-function gene in low-GCA parent and gain-of-function gene in high-GCA parent, encoding the putative Pseudo-Response Regulators, OsPRR37 and Ghd7, respectively. Overexpression of GCA1 in low-GCA parent significantly increases GCA effects in three traits. Our results demonstrate that two GCA loci associate with OsPRR37 and Ghd7 and reveal that the genes responsible for important agronomic traits could simultaneously account for GCA effects.

Genetic Architecture of Ear Fasciation in Maize (Zea mays) under QTL Scrutiny

MAIZE EAR FASCIATION

Knowledge of the genes affecting maize ear inflorescence may lead to better grain yield modeling. Maize ear fasciation, defined as abnormal flattened ears with high kernel row number, is a quantitative trait widely present in Portuguese maize landraces.

MATERIAL AND METHODS

Using a segregating population derived from an ear fasciation contrasting cross (consisting of 149 F2:3 families) we established a two location field trial using a complete randomized block design. Correlations and heritabilities for several ear fasciation-related traits and yield were determined. Quantitative Trait Loci (QTL) involved in the inheritance of those traits were identified and candidate genes for these QTL proposed.

RESULTS AND DISCUSSION

Ear fasciation broad-sense heritability was 0.73. Highly significant correlations were found between ear fasciation and some ear and cob diameters and row number traits. For the 23 yield and ear fasciation-related traits, 65 QTL were identified, out of which 11 were detected in both environments, while for the three principal components, five to six QTL were detected per environment. Detected QTL were distributed across 17 genomic regions and explained individually, 8.7% to 22.4% of the individual traits or principal components phenotypic variance. Several candidate genes for these QTL regions were proposed, such as bearded-ear1, branched silkless1, compact plant1, ramosa2, ramosa3, tasselseed4 and terminal ear1. However, many QTL mapped to regions without known candidate genes, indicating potential chromosomal regions not yet targeted for maize ear traits selection.

CONCLUSIONS

Portuguese maize germplasm represents a valuable source of genes or allelic variants for yield improvement and elucidation of the genetic basis of ear fasciation traits. Future studies should focus on fine mapping of the identified genomic regions with the aim of map-based cloning.

Genome properties and prospects of genomic prediction of hybrid performance in a breeding program of maize

Maize (Zea mays L.) serves as model plant for heterosis research and is the crop where hybrid breeding was pioneered. We analyzed genomic and phenotypic data of 1254 hybrids of a typical maize hybrid breeding program based on the important Dent × Flint heterotic pattern. Our main objectives were to investigate genome properties of the parental lines (e.g., allele frequencies, linkage disequilibrium, and phases) and examine the prospects of genomic prediction of hybrid performance. We found high consistency of linkage phases and large differences in allele frequencies between the Dent and Flint heterotic groups in pericentromeric regions. These results can be explained by the Hill-Robertson effect and support the hypothesis of differential fixation of alleles due to pseudo-overdominance in these regions. In pericentromeric regions we also found indications for consistent marker-QTL linkage between heterotic groups. With prediction methods GBLUP and BayesB, the cross-validation prediction accuracy ranged from 0.75 to 0.92 for grain yield and from 0.59 to 0.95 for grain moisture. The prediction accuracy of untested hybrids was highest, if both parents were parents of other hybrids in the training set, and lowest, if none of them were involved in any training set hybrid. Optimizing the composition of the training set in terms of number of lines and hybrids per line could further increase prediction accuracy. We conclude that genomic prediction facilitates a paradigm shift in hybrid breeding by focusing on the performance of experimental hybrids rather than the performance of parental lines in test crosses.

Genome-wide meta-analysis of maize heterosis reveals the potential role of additive gene expression at pericentromeric loci

BACKGROUND

The identification of QTL involved in heterosis formation is one approach to unravel the not yet fully understood genetic basis of heterosis - the improved agronomic performance of hybrid F1 plants compared to their inbred parents. The identification of candidate genes underlying a QTL is important both for developing markers and determining the molecular genetic basis of a trait, but remains difficult owing to the large number of genes often contained within individual QTL. To address this problem in heterosis analysis, we applied a meta-analysis strategy for grain yield (GY) of Zea mays L. as example, incorporating QTL-, hybrid field-, and parental gene expression data.

RESULTS

For the identification of genes underlying known heterotic QTL, we made use of tight associations between gene expression pattern and the trait of interest, identified by correlation analyses. Using this approach genes strongly associated with heterosis for GY were discovered to be clustered in pericentromeric regions of the complex maize genome. This suggests that expression differences of sequences in recombination-suppressed regions are important in the establishment of heterosis for GY in F1 hybrids and also in the conservation of heterosis for GY across genotypes. Importantly functional analysis of heterosis-associated genes from these genomic regions revealed over-representation of a number of functional classes, identifying key processes contributing to heterosis for GY. Based on the finding that the majority of the analyzed heterosis-associated genes were addtitively expressed, we propose a model referring to the influence of cis-regulatory variation on heterosis for GY by the compensation of fixed detrimental expression levels in parents.

CONCLUSIONS

The study highlights the utility of a meta-analysis approach that integrates phenotypic and multi-level molecular data to unravel complex traits in plants. It provides prospects for the identification of genes relevant for QTL, and also suggests a model for the potential role of additive expression in the formation and conservation of heterosis for GY via dominant, multigenic quantitative trait loci. Our findings contribute to a deeper understanding of the multifactorial phenomenon of heterosis, and thus to the breeding of new high yielding varieties.

Next generation characterisation of cereal genomes for marker discovery

Cereal crops form the bulk of the world’s food sources, and thus their importance cannot be understated. Crop breeding programs increasingly rely on high-resolution molecular genetic markers to accelerate the breeding process. The development of these markers is hampered by the complexity of some of the major cereal crop genomes, as well as the time and cost required. In this review, we address current and future methods available for the characterisation of cereal genomes, with an emphasis on faster and more cost effective approaches for genome sequencing and the development of markers for trait association and marker assisted selection (MAS) in crop breeding programs.

Identification of combining ability loci for five yield-related traits in maize using a set of testcrosses with introgression lines

Combining ability is essential for hybrid breeding in crops. However, the genetic basis of combining ability remains unclear and has been seldom investigated. Identifying molecular markers associated with this complex trait would help to understand its genetic basis and provide useful information for hybrid breeding in maize. In this study, we identified genetic loci of general combining ability (GCA) and specific combining ability (SCA) for five yield-related traits under three environments using a set of testcrosses with introgression lines (ILs). GCA or SCA of the five yield-related traits of the ILs was estimated by the performance of testcrosses with four testers from different heterotic groups. Genetic correlations between GCA of the traits and the corresponding traits per se were not significant or not strong, suggesting that the genetic basis between them is different. A total of 56 significant loci for GCA and 21 loci for SCA were commonly identified in at least two environments, and only 5 loci were simultaneously controlling GCA and SCA, indicating that the genetic basis of GCA and SCA is different. For all of the traits investigated, positive and significant correlations between the number of GCA loci in the ILs and the performance of the corresponding GCA of the ILs were detected, implying that pyramiding GCA loci would have positive effect on the performance of GCA. Results in this study would be useful for maize hybrid breeding.

Performance prediction of F1 hybrids between recombinant inbred lines derived from two elite maize inbred lines

Selection of recombinant inbred lines (RILs) from elite hybrids is a key method in maize breeding especially in developing countries. The RILs are normally derived by repeated self-pollination and selection. In this study, we first investigated the accuracy of different models in predicting the performance of F(1) hybrids between RILs derived from two elite maize inbred lines Zong3 and 87-1, and then compared these models through simulation using a wider range of genetic models. Results indicated that appropriate prediction models depended on genetic architecture, e.g., combined model using breeding value and genome-wide prediction (BV+GWP) has the highest prediction accuracy for high V(D)/V(A) ratio (>0.5) traits. Theoretical studies demonstrated that different components of genetic variance were captured by different prediction models, which in turn explained the accuracy of these models in predicting the F(1) hybrid performance. Based on genome-wide prediction model (GWP), 114 untested F(1) hybrids possibly having higher grain yield than the original F(1) hybrid Yuyu22 (the single cross between Zong3 and 87-1) have been identified and recommended for further field test.

Genetic diversity analysis of elite European maize (Zea mays L.) inbred lines using AFLP, SSR, and SNP markers reveals ascertainment bias for a subset of SNPs

Recent advances in high-throughput sequencing technologies have triggered a shift toward single-nucleotide polymorphism (SNP) markers. A systematic bias can be introduced if SNPs are ascertained in a small panel of genotypes and then used for characterizing a larger population (ascertainment bias). With the objective of evaluating a potential ascertainment bias of the Illumina MaizeSNP50 array with respect to elite European maize dent and flint inbred lines, we compared the genetic diversity among these materials based on 731 amplified fragment length polymorphisms (AFLPs), 186 simple sequence repeats (SSRs), 41,434 SNPs of the MaizeSNP50 array (SNP-A), and two subsets of it, i.e., 30,068 Panzea (SNP-P) and 11,366 Syngenta markers (SNP-S). We evaluated the bias effects on major allele frequency, allele number, gene diversity, modified Roger's distance (MRD), and on molecular variance (AMOVA). We revealed ascertainment bias in SNP-A, compared to AFLPs and SSRs. It affected especially European flint lines analyzed with markers (SNP-S) specifically developed to maximize differences among North American dent germplasm. The bias affected all genetic parameters, but did not substantially alter the relative distances between inbred lines within groups. For these reasons, we conclude that the SNP markers of the MaizeSNP50 array can be employed for breeding purposes in the investigated material. However, attention should be paid in case of comparisons between genotypes belonging to different heterotic groups. In this case, it is advisable to prefer a marker subset with potentially low ascertainment bias, like in our case the SNP-P marker set.

Genomic prediction of hybrid performance in maize with models incorporating dominance and population specific marker effects

Identifying high performing hybrids is an essential part of every maize breeding program. Genomic prediction of maize hybrid performance allows to identify promising hybrids, when they themselves or other hybrids produced from their parents were not tested in field trials. Using simulations, we investigated the effects of marker density (10, 1, 0.3 marker per mega base pair, Mbp(-1)), convergent or divergent parental populations, number of parents tested in other combinations (2, 1, 0), genetic model (including population-specific and/or dominance marker effects or not), and estimation method (GBLUP or BayesB) on the prediction accuracy. We based our simulations on marker genotypes of Central European flint and dent inbred lines from an ongoing maize breeding program. To simulate convergent or divergent parent populations, we generated phenotypes by assigning QTL to markers with similar or very different allele frequencies in both pools, respectively. Prediction accuracies increased with marker density and number of parents tested and were higher under divergent compared with convergent parental populations. Modeling marker effects as population-specific slightly improved prediction accuracy under lower marker densities (1 and 0.3 Mbp(-1)). This indicated that modeling marker effects as population-specific will be most beneficial under low linkage disequilibrium. Incorporating dominance effects improved prediction accuracies considerably for convergent parent populations, where dominance results in major contributions of SCA effects to the genetic variance among inter-population hybrids. While the general trends regarding the effects of the aforementioned influence factors on prediction accuracy were similar for GBLUP and BayesB, the latter method produced significantly higher accuracies for models incorporating dominance.

Genomic prediction of hybrid performance in maize with models incorporating dominance and population specific marker effects

Identifying high performing hybrids is an essential part of every maize breeding program. Genomic prediction of maize hybrid performance allows to identify promising hybrids, when they themselves or other hybrids produced from their parents were not tested in field trials. Using simulations, we investigated the effects of marker density (10, 1, 0.3 marker per mega base pair, Mbp(-1)), convergent or divergent parental populations, number of parents tested in other combinations (2, 1, 0), genetic model (including population-specific and/or dominance marker effects or not), and estimation method (GBLUP or BayesB) on the prediction accuracy. We based our simulations on marker genotypes of Central European flint and dent inbred lines from an ongoing maize breeding program. To simulate convergent or divergent parent populations, we generated phenotypes by assigning QTL to markers with similar or very different allele frequencies in both pools, respectively. Prediction accuracies increased with marker density and number of parents tested and were higher under divergent compared with convergent parental populations. Modeling marker effects as population-specific slightly improved prediction accuracy under lower marker densities (1 and 0.3 Mbp(-1)). This indicated that modeling marker effects as population-specific will be most beneficial under low linkage disequilibrium. Incorporating dominance effects improved prediction accuracies considerably for convergent parent populations, where dominance results in major contributions of SCA effects to the genetic variance among inter-population hybrids. While the general trends regarding the effects of the aforementioned influence factors on prediction accuracy were similar for GBLUP and BayesB, the latter method produced significantly higher accuracies for models incorporating dominance.

Partial least squares regression, support vector machine regression, and transcriptome-based distances for prediction of maize hybrid performance with gene expression data

The performance of hybrids can be predicted with gene expression data from their parental inbred lines. Implementing such prediction approaches in breeding programs promises to increase the efficiency of hybrid breeding. The objectives of our study were to compare the accuracy of prediction models employing multiple linear regression (MLR), partial least squares regression (PLS), support vector machine regression (SVM), and transcriptome-based distances (D(B)). For a factorial of 7 flint and 14 dent maize lines, the grain yield of the hybrids was assessed and the gene expression of the parental lines was profiled with a 56k microarray. The accuracy of the prediction models was measured by the correlation between predicted and observed yield employing two cross-validation schemes. The first modeled the prediction of hybrids when testcross data are available for both parental lines (type 2 hybrids), and the second modeled the prediction of hybrids when no testcross data for the parental lines were available (type 0 hybrids). MLR, SVM, and PLS resulted in a high correlation between predicted and observed yield for type 2 hybrids, whereas for type 0 hybrids D(B) had greater prediction accuracy. The regression methods were robust to the choice of the set of profiled genes and required only a few hundred genes. In contrast, for an accurate hybrid prediction with D(B), 1,000-1,500 genes were required, and the prediction accuracy depended strongly on the set of profiled genes. We conclude that for prediction within one set of genetic material MLR is a promising approach, and for transferring prediction models from one set of genetic material to a related one, the transcriptome-based distance D(B) is most promising.

QTL mapping of combining ability and heterosis of agronomic traits in rice backcross recombinant inbred lines and hybrid crosses

Background

Combining ability effects are very effective genetic parameters in deciding the next phase of breeding programs. Although some breeding strategies on the basis of evaluating combining ability have been utilized extensively in hybrid breeding, little is known about the genetic basis of combining ability. Combining ability is a complex trait that is controlled by polygenes. With the advent and development of molecular markers, it is feasible to evaluate the genetic bases of combining ability and heterosis of elite rice hybrids through QTL analysis.

Methodology/Principal Findings

In the present study, we first developed a QTL-mapping method for dissecting combining ability and heterosis of agronomic traits. With three testcross populations and a BCRIL population in rice, biometric and QTL analyses were conducted for ten agronomic traits. The significance of general combining ability and special combining ability for most of the traits indicated the importance of both additive and non-additive effects on expression levels. A large number of additive effect QTLs associated with performance per se of BCRIL and general combining ability, and dominant effect QTLs associated with special combining ability and heterosis were identified for the ten traits.

Conclusions/Significance

The combining ability of agronomic traits could be analyzed by the QTL mapping method. The characteristics revealed by the QTLs for combining ability of agronomic traits were similar with those by multitudinous QTLs for agronomic traits with performance per se of BCRIL. Several QTLs (1–6 in this study) were identified for each trait for combining ability. It demonstrated that some of the QTLs were pleiotropic or linked tightly with each other. The identification of QTLs responsible for combining ability and heterosis in the present study provides valuable information for dissecting genetic basis of combining ability.

Dissecting grain yield pathways and their interactions with grain dry matter content by a two-step correlation approach with maize seedling transcriptome

BACKGROUND

The importance of maize for human and animal nutrition, but also as a source for bio-energy is rapidly increasing. Maize yield is a quantitative trait controlled by many genes with small effects, spread throughout the genome. The precise location of the genes and the identity of the gene networks underlying maize grain yield is unknown. The objective of our study was to contribute to the knowledge of these genes and gene networks by transcription profiling with microarrays.

RESULTS

We assessed the grain yield and grain dry matter content (an indicator for early maturity) of 98 maize hybrids in multi-environment field trials. The gene expression in seedlings of the parental inbred lines, which have four different genetic backgrounds, was assessed with genome-scale oligonucleotide arrays. We identified genes associated with grain yield and grain dry matter content using a newly developed two-step correlation approach and found overlapping gene networks for both traits. The underlying metabolic pathways and biological processes were elucidated. Genes involved in sucrose degradation and glycolysis, as well as genes involved in cell expansion and endocycle were found to be associated with grain yield.

CONCLUSIONS

Our results indicate that the capability of providing energy and substrates, as well as expanding the cell at the seedling stage, highly influences the grain yield of hybrids. Knowledge of these genes underlying grain yield in maize can contribute to the development of new high yielding varieties.

Mixed model approaches for the identification of QTLs within a maize hybrid breeding program

Two outlines for mixed model based approaches to quantitative trait locus (QTL) mapping in existing maize hybrid selection programs are presented: a restricted maximum likelihood (REML) and a Bayesian Markov Chain Monte Carlo (MCMC) approach. The methods use the in-silico-mapping procedure developed by Parisseaux and Bernardo (2004) as a starting point. The original single-point approach is extended to a multi-point approach that facilitates interval mapping procedures. For computational and conceptual reasons, we partition the full set of relationships from founders to parents of hybrids into two types of relations by defining so-called intermediate founders. QTL effects are defined in terms of those intermediate founders. Marker based identity by descent relationships between intermediate founders define structuring matrices for the QTL effects that change along the genome. The dimension of the vector of QTL effects is reduced by the fact that there are fewer intermediate founders than parents. Furthermore, additional reduction in the number of QTL effects follows from the identification of founder groups by various algorithms. As a result, we obtain a powerful mixed model based statistical framework to identify QTLs in genetic backgrounds relevant to the elite germplasm of a commercial breeding program. The identification of such QTLs will provide the foundation for effective marker assisted and genome wide selection strategies. Analyses of an example data set show that QTLs are primarily identified in different heterotic groups and point to complementation of additive QTL effects as an important factor in hybrid performance.

Correlation between parental transcriptome and field data for the characterization of heterosis in Zea mays L

Heterosis is widely exploited in plant breeding, although its molecular basis is still not fully understood. For the characterization of this phenomenon and the development of transcriptome-based methods to predict hybrid performance (HP), we applied a microarray (46k) analysis of 21 European maize (Zea mays L.), 14 dent and 7 flint parental inbred lines. Expression profiles of the parental inbreds at the seedling stage were correlated with grain yield (GY) and grain dry matter content (GDMC) of 98 flint x dent factorial crosses at six locations. We observed highly significant correlations of the parental expression levels of certain differentially expressed genes with heterosis and HP for GY and also with HP for GDMC. This strong correlation provided first evidence toward a prediction potential of the genes and their expression levels. The identified gene set based on the parental transcriptome data revealed functional characteristics of HP and heterosis. Gene ontology (GO) analyses were performed to compare genes correlated with their expression pattern to HP for GY and GDMC, respectively. Between these gene groups, mostly different functional classes of genes were found to be enriched or underrepresented. The phenomenon of heterosis was characterized by the over- and underrepresentation of specific GO terms among heterosis-correlated genes.

Prediction of hybrid performance in maize using molecular markers and joint analyses of hybrids and parental inbreds

The identification of superior hybrids is important for the success of a hybrid breeding program. However, field evaluation of all possible crosses among inbred lines requires extremely large resources. Therefore, efforts have been made to predict hybrid performance (HP) by using field data of related genotypes and molecular markers. In the present study, the main objective was to assess the usefulness of pedigree information in combination with the covariance between general combining ability (GCA) and per se performance of parental lines for HP prediction. In addition, we compared the prediction efficiency of AFLP and SSR marker data, estimated marker effects separately for reciprocal allelic configurations (among heterotic groups) of heterozygous marker loci in hybrids, and imputed missing AFLP marker data for marker-based HP prediction. Unbalanced field data of 400 maize dent x flint hybrids from 9 factorials and of 79 inbred parents were subjected to joint analyses with mixed linear models. The inbreds were genotyped with 910 AFLP and 256 SSR markers. Efficiency of prediction (R (2)) was estimated by cross-validation for hybrids having no or one parent evaluated in testcrosses. Best linear unbiased prediction of GCA and specific combining ability resulted in the highest efficiencies for HP prediction for both traits (R (2) = 0.6-0.9), if pedigree and line per se data were used. However, without such data, HP for grain yield was more efficiently predicted using molecular markers. The additional modifications of the marker-based approaches had no clear effect. Our study showed the high potential of joint analyses of hybrids and parental inbred lines for the prediction of performance of untested hybrids.

Transcriptome-based distance measures for grouping of germplasm and prediction of hybrid performance in maize

Grouping of germplasm and prediction of hybrid performance and heterosis are important applications in hybrid breeding programs. Gene expression analysis is a promising tool to achieve both tasks efficiently. Our objectives were to (1) investigate distance measures based on transcription profiles, (2) compare these with genetic distances based on AFLP markers, and (3) assess the suitability of transcriptome-based distances for grouping of germplasm and prediction of hybrid performance and heterosis in maize. We analyzed transcription profiles from seedlings of the 21 parental maize lines of a 7 x 14 factorial with a 46-k oligonucleotide array. The hybrid performance and heterosis of the 98 hybrids were assessed in field trials. In cluster and principal coordinate analyses for germplasm grouping, the transcriptome-based distances were as powerful as the genetic distances for separating flint from dent inbreds. Cross validation showed that prediction of hybrid performance with transcriptome-based distances using selected markers was more precise than earlier prediction models using DNA markers or general combining ability estimates using field data. Our results suggest that transcriptome-based prediction of hybrid performance and heterosis has a great potential to improve the efficiency of maize hybrid breeding programs.

Improved heterosis prediction by combining information on DNA- and metabolic markers

BACKGROUND

Hybrids represent a cornerstone in the success story of breeding programs. The fundamental principle underlying this success is the phenomenon of hybrid vigour, or heterosis. It describes an advantage of the offspring as compared to the two parental lines with respect to parameters such as growth and resistance against abiotic or biotic stress. Dominance, overdominance or epistasis based models are commonly used explanations.

CONCLUSION/SIGNIFICANCE

The heterosis level is clearly a function of the combination of the parents used for offspring production. This results in a major challenge for plant breeders, as usually several thousand combinations of parents have to be tested for identifying the best combinations. Thus, any approach to reliably predict heterosis levels based on properties of the parental lines would be highly beneficial for plant breeding.

METHODOLOGY/PRINCIPAL FINDINGS

Recently, genetic data have been used to predict heterosis. Here we show that a combination of parental genetic and metabolic markers, identified via feature selection and minimum-description-length based regression methods, significantly improves the prediction of biomass heterosis in resulting offspring. These findings will help furthering our understanding of the molecular basis of heterosis, revealing, for instance, the presence of nonlinear genotype-phenotype relationships. In addition, we describe a possible approach for accelerated selection in plant breeding.

Molecular marker-based prediction of hybrid performance in maize using unbalanced data from multiple experiments with factorial crosses

In hybrid breeding, the prediction of hybrid performance (HP) is extremely important as it is difficult to evaluate inbred lines in numerous cross combinations. Recent developments such as doubled haploid production and molecular marker technologies have enhanced the prospects of marker-based HP prediction to accelerate the breeding process. Our objectives were to (1) predict HP using a combined analysis of hybrids and parental lines from a breeding program, (2) evaluate the use of molecular markers in addition to phenotypic and pedigree data, (3) evaluate the combination of line per se data with marker-based estimates, (4) study the effect of the number of tested parents, and (5) assess the advantage of haplotype blocks. An unbalanced dataset of 400 hybrids from 9 factorial crosses tested in different experiments and data of 79 inbred parents were subjected to combined analyses with a mixed linear model. Marker data of the inbreds were obtained with 20 AFLP primer-enzyme combinations. Cross-validation was used to assess the performance prediction of hybrids of which no or only one parental line was testcross evaluated. For HP prediction, the highest proportion of explained variance (R (2)), 46% for grain yield (GY) and 70% for grain dry matter content (GDMC), was obtained from line per se best linear unbiased prediction (BLUP) estimates plus marker effects associated with mid-parent heterosis (TEAM-LM). Our study demonstrated that HP was efficiently predicted using molecular markers even for GY when testcross data of both parents are not available. This can help in improving greatly the efficiency of commercial hybrid breeding programs.

Impact of interpopulation divergence on additive and dominance variance in hybrid populations

We present a theoretical proof that the ratio of the dominance vs. the additive variance decreases with increasing genetic divergence between two populations. While the dominance variance is the major component of the variance due to specific combining ability (sigma(SCA)(2)), the additive variance is the major component of the variance due to general combining ability (sigma(GCA)(2)). Therefore, we conclude that interpopulation improvement becomes more efficient with divergent than with genetically similar heterotic groups, because performance of superior hybrids can be predicted on the basis of general combining ability effects.

Prediction of single-cross hybrid performance in maize using haplotype blocks associated with QTL for grain yield

Marker-based prediction of hybrid performance facilitates the identification of untested single-cross hybrids with superior yield performance. Our objectives were to (1) determine the haplotype block structure of experimental germplasm from a hybrid maize breeding program, (2) develop models for hybrid performance prediction based on haplotype blocks, and (3) compare hybrid performance prediction based on haplotype blocks with other approaches, based on single AFLP markers or general combining ability (GCA), under a validation scenario relevant for practical breeding. In total, 270 hybrids were evaluated for grain yield in four Dent x Flint factorial mating experiments. Their parental inbred lines were genotyped with 20 AFLP primer-enzyme combinations. Adjacent marker loci were combined into haplotype blocks. Hybrid performance was predicted on basis of single marker loci and haplotype blocks. Prediction based on variable haplotype block length resulted in an improved prediction of hybrid performance compared with the use of single AFLP markers. Estimates of prediction efficiency (R(2)) ranged from 0.305 to 0.889 for marker-based prediction and from 0.465 to 0.898 for GCA-based prediction. For inter-group hybrids with predominance of general over specific combining ability, the hybrid prediction from GCA effects was efficient in identifying promising hybrids. Considering the advantage of haplotype block approaches over single marker approaches for the prediction of inter-group hybrids, we see a high potential to substantially improve the efficiency of hybrid breeding programs.

Hybrid maize breeding with doubled haploids: II. Optimum type and number of testers in two-stage selection for general combining ability

Optimum allocation of test resources is of crucial importance for the efficiency of breeding programs. Our objectives were to (1) determine the optimum allocation of the number of lines, test locations, as well as number and type of testers in hybrid maize breeding using doubled haploids with two breeding strategies for improvement of general combining ability (GCA), (2) compare the maximum selection gain (DeltaG) achievable under both strategies, and (3) give recommendations for the optimum implementation of doubled haploids in commercial hybrid maize breeding. We calculated DeltaG by numerical integration for two two-stage selection strategies with evaluation of (1) testcross performance in both stages (BS1) or (2) line per se performance in the first stage followed by testcross performance in the second stage (BS2). Different assumptions were made regarding the budget, variance components (VCs), and the correlation between line per se performance and GCA. Selection gain for GCA increased with a broader genetic base of the tester. Hence, testers combining a large number of divergent lines are advantageous. However, in applied breeding programs, the use of single- or double-cross testers in the first and inbred testers in the second selection stage may be a good compromise between theoretical and practical requirements. With a correlation between line per se performance and GCA of 0.50, DeltaG for BS1 is about 5% higher than for BS2, if an economic weight of line per se performance is neglected. With increasing economic weight of line per se performance, relative efficiency of BS2 increased rapidly resulting in a superiority of BS2 over BS1 already for an economic weight for line per se performance larger than 0.1. Considering the importance of an economic seed production, an economic weight larger than 0.1 seems realistic indicating the necessity of separate breeding strategies for seed and pollen parent heterotic groups.