Cyril Bauland

Maximizing the Reliability of Genomic Selection by Optimizing the Calibration Set of Reference Individuals: Comparison of Methods in Two Diverse Groups of Maize Inbreds ( Zea mays L.)

Genomic selection refers to the use of genotypic information for predicting breeding values of selection candidates. A prediction formula is calibrated with the genotypes and phenotypes of reference individuals constituting the calibration set. The size and the composition of this set are essential parameters affecting the prediction reliabilities. The objective of this study was to maximize reliabilities by optimizing the calibration set. Different criteria based on the diversity or on the prediction error variance (PEV) derived from the realized additive relationship matrix–best linear unbiased predictions model (RA–BLUP) were used to select the reference individuals. For the latter, we considered the mean of the PEV of the contrasts between each selection candidate and the mean of the population (PEVmean) and the mean of the expected reliabilities of the same contrasts (CDmean). These criteria were tested with phenotypic data collected on two diversity panels of maize (Zea mays L.) genotyped with a 50k SNPs array. In the two panels, samples chosen based on CDmean gave higher reliabilities than random samples for various calibration set sizes. CDmean also appeared superior to PEVmean, which can be explained by the fact that it takes into account the reduction of variance due to the relatedness between individuals. Selected samples were close to optimality for a wide range of trait heritabilities, which suggests that the strategy presented here can efficiently sample subsets in panels of inbred lines. A script to optimize reference samples based on CDmean is available on request.

Intraspecific variation of recombination rate in maize.

BackgroundIn sexually reproducing organisms, meiotic crossovers ensure the proper segregation of chromosomes and contribute to genetic diversity by shuffling allelic combinations. Such genetic reassortment is exploited in breeding to combine favorable alleles, and in genetic research to identify genetic factors underlying traits of interest via linkage or association-based approaches. Crossover numbers and distributions along chromosomes vary between species, but little is known about their intraspecies variation.ResultsHere, we report on the variation of recombination rates between 22 European maize inbred lines that belong to the Dent and Flint gene pools. We genotype 23 doubled-haploid populations derived from crosses between these lines with a 50 k-SNP array and construct high-density genetic maps, showing good correspondence with the maize B73 genome sequence assembly. By aligning each genetic map to the B73 sequence, we obtain the recombination rates along chromosomes specific to each population. We identify significant differences in recombination rates at the genome-wide, chromosome, and intrachromosomal levels between populations, as well as significant variation for genome-wide recombination rates among maize lines. Crossover interference analysis using a two-pathway modeling framework reveals a negative association between recombination rate and interference strength.ConclusionsTo our knowledge, the present work provides the most comprehensive study on intraspecific variation of recombination rates and crossover interference strength in eukaryotes. Differences found in recombination rates will allow for selection of high or low recombining lines in crossing programs. Our methodology should pave the way for precise identification of genes controlling recombination rates in maize and other organisms.

Usefulness of Multiparental Populations of Maize (Zea mays L.) for Genome-Based Prediction

The efficiency of marker-assisted prediction of phenotypes has been studied intensively for different types of plant breeding populations. However, one remaining question is how to incorporate and counterbalance information from biparental and multiparental populations into model training for genome-wide prediction. To address this question, we evaluated testcross performance of 1652 doubled-haploid maize (Zea mays L.) lines that were genotyped with 56,110 single nucleotide polymorphism markers and phenotyped for five agronomic traits in four to six European environments. The lines are arranged in two diverse half-sib panels representing two major European heterotic germplasm pools. The data set contains 10 related biparental dent families and 11 related biparental flint families generated from crosses of maize lines important for European maize breeding. With this new data set we analyzed genome-based best linear unbiased prediction in different validation schemes and compositions of estimation and test sets. Further, we theoretically and empirically investigated marker linkage phases across multiparental populations. In general, predictive abilities similar to or higher than those within biparental families could be achieved by combining several half-sib families in the estimation set. For the majority of families, 375 half-sib lines in the estimation set were sufficient to reach the same predictive performance of biomass yield as an estimation set of 50 full-sib lines. In contrast, prediction across heterotic pools was not possible for most cases. Our findings are important for experimental design in genome-based prediction as they provide guidelines for the genetic structure and required sample size of data sets used for model training.

Linkage disequilibrium with linkage analysis of multiline crosses reveals different multiallelic QTL for hybrid performance in the flint and dent heterotic groups of maize.

Multiparental designs combined with dense genotyping of parents have been proposed as a way to increase the diversity and resolution of quantitative trait loci (QTL) mapping studies, using methods combining linkage disequilibrium information with linkage analysis (LDLA). Two new nested association mapping designs adapted to European conditions were derived from the complementary dent and flint heterotic groups of maize (Zea mays L.). Ten biparental dent families (N = 841) and 11 biparental flint families (N = 811) were genotyped with 56,110 single nucleotide polymorphism markers and evaluated as test crosses with the central line of the reciprocal design for biomass yield, plant height, and precocity. Alleles at candidate QTL were defined as (i) parental alleles, (ii) haplotypic identity by descent, and (iii) single-marker groupings. Between five and 16 QTL were detected depending on the model, trait, and genetic group considered. In the flint design, a major QTL (R2 = 27%) with pleiotropic effects was detected on chromosome 10, whereas other QTL displayed milder effects (R2 < 10%). On average, the LDLA models detected more QTL but generally explained lower percentages of variance, consistent with the fact that most QTL display complex allelic series. Only 15% of the QTL were common to the two designs. A joint analysis of the two designs detected between 15 and 21 QTL for the five traits. Of these, between 27 for silking date and 41% for tasseling date were significant in both groups. Favorable allelic effects detected in both groups open perspectives for improving biomass production.

Genome-wide analysis of yield in Europe: allelic effects as functions of drought and heat scenarios

A genome-wide analysis of maize yield in identify genomic regions associated with adaptation to scenarios with drought or heat stresses. Assessing the genetic variability of plant performance under heat and drought scenarios can contribute to reduce the negative effects of climate change. We propose here an approach that consisted of (1) clustering time courses of environmental variables simulated by a crop model in current (35 years × 55 sites) and future conditions into six scenarios of temperature and water deficit as experienced by maize (Zea mays L.) plants; (2) performing 29 field experiments in contrasting conditions across Europe with 244 maize hybrids; (3) assigning individual experiments to scenarios based on environmental conditions as measured in each field experiment; frequencies of temperature scenarios in our experiments corresponded to future heat scenarios (+5°C); (4) analyzing the genetic variation of plant performance for each environmental scenario. Forty-eight quantitative trait loci (QTLs) of yield were identified by association genetics using a multi-environment multi-locus model. Eight and twelve QTLs were associated to tolerances to heat and drought stresses because they were specific to hot and dry scenarios, respectively, with low or even negative allelic effects in favorable scenarios. Twenty-four QTLs improved yield in favorable conditions but showed nonsignificant effects under stress; they were therefore associated with higher sensitivity. Our approach showed a pattern of QTL effects expressed as functions of environmental variables and scenarios, allowing us to suggest hypotheses for mechanisms and candidate genes underlying each QTL. It can be used for assessing the performance of genotypes and the contribution of genomic regions under current and future stress situations and to accelerate breeding for drought-prone environments.

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

Understanding genetic architecture of hybrid performances is important for species showing heterosis. Giraud et al. evaluated an... 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.

Linkage Analysis and Association Mapping QTL Detection Models for Hybrids Between Multiparental Populations from Two Heterotic Groups: Application to Biomass Production in Maize (Zea mays L.)

Identification of quantitative trait loci (QTL) involved in the variation of hybrid value is of key importance for cross-pollinated species such as maize (Zea mays L.). In a companion paper, we illustrated a new QTL mapping population design involving a factorial mating between two multiparental segregating populations. Six biparental line populations were developed from four founder lines in the Dent and Flint heterotic groups. They were crossed to produce 951 hybrids and evaluated for silage performances. Previously, a linkage analysis (LA) model that assumes each founder line carries a different allele was used to detect QTL involved in General and Specific Combining Abilities (GCA and SCA, respectively) of hybrid value. This previously introduced model requires the estimation of numerous effects per locus, potentially affecting QTL detection power. Using the same design, we compared this “Founder alleles” model to two more parsimonious models, which assume that (i) identity in state at SNP alleles from the same heterotic group implies identity by descent (IBD) at linked QTL (“SNP within-group” model) or (ii) identity in state implies IBD, regardless of population origin of the alleles (“Hybrid genotype” model). This last model assumes biallelic QTL with equal effects in each group. It detected more QTL on average than the two other models but explained lower percentages of variance. The “SNP within-group” model appeared to be a good compromise between the two other models. These results confirm the divergence between the Dent and Flint groups. They also illustrate the need to adapt the QTL detection model to the complexity of the allelic variation, which depends on the trait, the QTL, and the divergence between the heterotic groups.