Stéphane D. Nicolas

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.

Homeologous Recombination Plays a Major Role in Chromosome Rearrangements That Occur During Meiosis of Brassica napus Haploids

Chromosomal rearrangements can be triggered by recombination between distinct but related regions. Brassica napus (AACC; 2n = 38) is a recent allopolyploid species whose progenitor genomes are widely replicated. In this article, we analyze the extent to which chromosomal rearrangements originate from homeologous recombination during meiosis of haploid B. napus (n = 19) by genotyping progenies of haploid × euploid B. napus with molecular markers. Our study focuses on three pairs of homeologous regions selected for their differing levels of divergence (N1/N11, N3/N13, and N9/N18). We show that a high number of chromosomal rearrangements occur during meiosis of B. napus haploid and are transmitted by first division restitution (FDR)-like unreduced gametes to their progeny; half of the progeny of Darmor-bzh haploids display duplications and/or losses in the chromosomal regions being studied. We demonstrate that half of these rearrangements are due to recombination between regions of primary homeology, which represents a 10- to 100-fold increase compared to the frequency of homeologous recombination measured in euploid lines. Some of the other rearrangements certainly result from recombination between paralogous regions because we observed an average of one to two autosyndetic A–A and/or C–C bivalents at metaphase I of the B. napus haploid. These results are discussed in the context of genome evolution of B. napus.

Genetic Regulation of Meiotic Cross-Overs between Related Genomes in Brassica napus Haploids and Hybrids

Although the genetic regulation of recombination in allopolyploid species plays a pivotal role in evolution and plant breeding, it has received little recent attention, except in wheat (Triticum aestivum). PrBn is the main locus that determines the number of nonhomologous associations during meiosis of microspore cultured Brassica napus haploids (AC; 19 chromosomes). In this study, we examined the role played by PrBn in recombination. We generated two haploid × euploid populations using two B. napus haploids with differing PrBn (and interacting genes) activity. We analyzed molecular marker transmission in these two populations to compare genetic changes, which have arisen during meiosis. We found that cross-over number in these two genotypes was significantly different but that cross-overs between nonhomologous chromosomes showed roughly the same distribution pattern. We then examined genetic recombination along a pair of A chromosomes during meiosis of B. rapa × B. napus AAC and AACC hybrids that were produced with the same two B. napus genotypes. We observed significant genotypic variation in cross-over rates between the two AAC hybrids but no difference between the two AACC hybrids. Overall, our results show that PrBn changes the rate of recombination between nonhomologous chromosomes during meiosis of B. napus haploids and also affects homologous recombination with an effect that depends on plant karyotype.

Recovering Power in Association Mapping Panels with Variable Levels of Linkage Disequilibrium

Association mapping has permitted the discovery of major QTL in many species. It can be applied to existing populations and, as a consequence, it is generally necessary to take into account structure and relatedness among individuals in the statistical model to control false positives. We analytically studied power in association studies by computing noncentrality parameter of the tests and its relationship with parameters characterizing diversity (genetic differentiation between groups and allele frequencies) and kinship between individuals. Investigation of three different maize diversity panels genotyped with the 50k SNPs array highlighted contrasted average power among panels and revealed gaps of power of classical mixed models in regions with high linkage disequilibrium (LD). These gaps could be related to the fact that markers are used for both testing association and estimating relatedness. We thus considered two alternative approaches to estimating the kinship matrix to recover power in regions of high LD. In the first one, we estimated the kinship with all the markers that are not located on the same chromosome than the tested SNP. In the second one, correlation between markers was taken into account to weight the contribution of each marker to the kinship. Simulations revealed that these two approaches were efficient to control false positives and were more powerful than classical models.

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.

Exploiting the Genetic Diversity of Maize using a Combined Metabolomic, Enzyme Activity Profiling, and Metabolic Modelling Approach to Link Leaf Physiology to Kernel Yield

A metabolomic, biochemical, fluxomic, and metabolic modeling approach developed using 19 genetically distant maize lines links leaf physiology to kernel yield. A combined metabolomic, biochemical, fluxomic, and metabolic modeling approach was developed using 19 genetically distant maize (Zea mays) lines from Europe and America. Considerable differences were detected between the lines when leaf metabolic profiles and activities of the main enzymes involved in primary metabolism were compared. During grain filling, the leaf metabolic composition appeared to be a reliable marker, allowing a classification matching the genetic diversity of the lines. During the same period, there was a significant correlation between the genetic distance of the lines and the activities of enzymes involved in carbon metabolism, notably glycolysis. Although large differences were observed in terms of leaf metabolic fluxes, these variations were not tightly linked to the genome structure of the lines. Both correlation studies and metabolic network analyses allowed the description of a maize ideotype with a high grain yield potential. Such an ideotype is characterized by low accumulation of soluble amino acids and carbohydrates in the leaves and high activity of enzymes involved in the C4 photosynthetic pathway and in the biosynthesis of amino acids derived from glutamate. Chlorogenates appear to be important markers that can be used to select for maize lines that produce larger kernels.

Independent introductions and admixtures have contributed to adaptation of European maize and its American counterparts

Through the local selection of landraces, humans have guided the adaptation of crops to a vast range of climatic and ecological conditions. This is particularly true of maize, which was domesticated in a restricted area of Mexico but now displays one of the broadest cultivated ranges worldwide. Here, we sequenced 67 genomes with an average sequencing depth of 18x to document routes of introduction, admixture and selective history of European maize and its American counterparts. To avoid the confounding effects of recent breeding, we targeted germplasm (lines) directly derived from landraces. Among our lines, we discovered 22,294,769 SNPs and between 0.9% to 4.1% residual heterozygosity. Using a segmentation method, we identified 6,978 segments of unexpectedly high rate of heterozygosity. These segments point to genes potentially involved in inbreeding depression, and to a lesser extent to the presence of structural variants. Genetic structuring and inferences of historical splits revealed 5 genetic groups and two independent European introductions, with modest bottleneck signatures. Our results further revealed admixtures between distinct sources that have contributed to the establishment of 3 groups at intermediate latitudes in North America and Europe. We combined differentiation- and diversity-based statistics to identify both genes and gene networks displaying strong signals of selection. These include genes/gene networks involved in flowering time, drought and cold tolerance, plant defense and starch properties. Overall, our results provide novel insights into the evolutionary history of European maize and highlight a major role of admixture in environmental adaptation, paralleling recent findings in humans.