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

Genetic architecture of phosphorus use efficiency across diverse environmental conditions: insights from maize elite and landrace lines

Phosphorus is an essential nutrient for all crops. Thus, a better understanding of the genetic control of phosphorus use efficiency evident in physiological, developmental, and morphological traits and its environmental plasticity is required to establish the basis for maintaining or enhancing yield while making agriculture more sustainable. In this study, we utilized a diverse panel of maize (Zea mays L.), including 398 elite and landrace lines, phenotyped across three environments and two phosphorus fertilization treatments. We performed genome-wide association mapping for 13 traits, including phosphorus uptake and allocation, that showed a strong environment dependency in their expression. Our results highlight the complex genetic architecture of phosphorus use efficiency as well as the substantial differences between the evaluated genetic backgrounds. Despite harboring more of the identified quantitative trait loci, almost all of the favorable alleles from landraces were found to be present in at least one of the two elite heterotic groups. Notably, we also observed trait-specific genetic control even among biologically related characteristics, as well as a substantial plasticity of the genetic architecture of several traits in response to the environment and phosphorus fertilization. Collectively, our work illustrates the difficulties in improving phosphorus use efficiency, but also presents possible solutions for the future contribution of plant breeding to improve the phosphorus cycle.

metaGE: Investigating genotype x environment interactions through GWAS meta-analysis

Elucidating the genetic components of plant genotype-by-environment interactions is of key importance in the context of increasing climatic instability, diversification of agricultural practices and pest pressure due to phytosanitary treatment limitations. The genotypic response to environmental stresses can be investigated through multi-environment trials (METs). However, genome-wide association studies (GWAS) of MET data are significantly more complex than that of single environments. In this context, we introduce metaGE, a flexible and computationally efficient meta-analysis approach for jointly analyzing single-environment GWAS of any MET experiment. The metaGE procedure accounts for the heterogeneity of quantitative trait loci (QTL) effects across the environmental conditions and allows the detection of QTL whose allelic effect variations are strongly correlated to environmental cofactors. We evaluated the performance of the proposed methodology and compared it to two competing procedures through simulations. We also applied metaGE to two emblematic examples: the detection of flowering QTLs whose effects are modulated by competition in Arabidopsis and the detection of yield QTLs impacted by drought stresses in maize. The procedure identified known and new QTLs, providing valuable insights into the genetic architecture of complex traits and QTL effects dependent on environmental stress conditions. The whole statistical approach is available as an R package.

Combined Genome-Wide Association Study and Linkage Analysis for Mining Candidate Genes for the Kernel Row Number in Maize (Zea mays L.)

Kernel row number (KRN) is one of the key traits that significantly affect maize yield and productivity. Therefore, investigating the candidate genes and their functions in regulating KRN provides a theoretical basis and practical direction for genetic improvement in maize breeding, which is vital for increasing maize yield and understanding domestication. In this study, three recombinant inbred line (RIL) populations were developed using the parental lines AN20, YML1218, CM395, and Ye107, resulting in a multiparent population comprising a total of 490 F9 RILs. Phenotypic evaluation of the RILs for KRN was performed in three distinct environments. The heritability estimates of the RILs ranged from 81.40% to 84.16%. Genotyping-by-sequencing (GBS) of RILs identified 569,529 high-quality single nucleotide polymorphisms (SNPs). Combined genome-wide association study (GWAS) and linkage analyses revealed 120 SNPs and 22 quantitative trait loci (QTLs) which were significantly associated with KRN in maize. Furthermore, two novel candidate genes, Zm00001d042733 and Zm00001d042735, regulating KRN in maize were identified, which were located in close proximity to the significant SNP3-178,487,003 and overlapping the interval of QTL qKRN3-1. Zm00001d042733 encodes ubiquitin carboxyl-terminal hydrolase and Zm00001d042735 encodes the Arabidopsis Tóxicos en Levadura family of proteins. This study identified novel candidate loci and established a theoretical foundation for further functional validation of candidate genes. These findings deepen our comprehension of the genetic mechanisms that underpin KRN and offer potential applications of KRN-related strategies in developing maize varieties with higher yield.

Modeling QTL-by-environment interactions for multi-parent populations

Multi-parent populations (MPPs) are attractive for genetic and breeding studies because they combine genetic diversity with an easy-to-control population structure. Most methods for mapping QTLs in MPPs focus on the detection of QTLs in single environments. Little attention has been given to mapping QTLs in multienvironment trials (METs) and to detecting and modeling QTL-by-environment interactions (QEIs). We present mixed model approaches for the detection and modeling of consistent versus environment-dependent QTLs, i.e., QTL-by-environment interaction (QEI). QTL effects are assumed to be normally distributed with variances expressing consistency or dependence on environments and families. The entries of the corresponding design matrices are functions of identity-by-descent (IBD) probabilities between parents and offspring and follow from the parental origin of offspring DNA. A polygenic effect is added to the models to account for background genetic variation. We illustrate the wide applicability of our method by analyzing several public MPP datasets with observations from METs. The examples include diallel, nested association mapping (NAM), and multi-parent advanced inter-cross (MAGIC) populations. The results of our approach compare favorably with those of previous studies that used tailored methods.

Genetic mapping reveals new loci and alleles for flowering time and plant height using the double round-robin population of barley

Flowering time and plant height are two critical determinants of yield potential in barley (Hordeum vulgare). Despite their role in plant physiological regulation, a complete overview of the genetic complexity of flowering time and plant height regulation in barley is still lacking. Using a double round-robin population originated from the crossings of 23 diverse parental inbred lines, we aimed to determine the variance components in the regulation of flowering time and plant height in barley as well as to identify new genetic variants by single and multi-population QTL analyses and allele mining. Despite similar genotypic variance, we observed higher environmental variance components for plant height than flowering time. Furthermore, we detected new QTLs for flowering time and plant height. Finally, we identified a new functional allelic variant of the main regulatory gene Ppd-H1. Our results show that the genetic architecture of flowering time and plant height might be more complex than reported earlier and that a number of undetected, small effect, or low-frequency genetic variants underlie the control of these two traits.

Characterization of adaptation mechanisms in sorghum using a multireference back-cross nested association mapping design and envirotyping

Identifying the genetic factors impacting the adaptation of crops to environmental conditions is of key interest for conservation and selection purposes. It can be achieved using population genomics, and evolutionary or quantitative genetics. Here we present a sorghum multireference back-cross nested association mapping population composed of 3,901 lines produced by crossing 24 diverse parents to 3 elite parents from West and Central Africa-back-cross nested association mapping. The population was phenotyped in environments characterized by differences in photoperiod, rainfall pattern, temperature levels, and soil fertility. To integrate the multiparental and multi-environmental dimension of our data we proposed a new approach for quantitative trait loci (QTL) detection and parental effect estimation. We extended our model to estimate QTL effect sensitivity to environmental covariates, which facilitated the integration of envirotyping data. Our models allowed spatial projections of the QTL effects in agro-ecologies of interest. We utilized this strategy to analyze the genetic architecture of flowering time and plant height, which represents key adaptation mechanisms in environments like West Africa. Our results allowed a better characterization of well-known genomic regions influencing flowering time concerning their response to photoperiod with Ma6 and Ma1 being photoperiod-sensitive and the region of possible candidate gene Elf3 being photoperiod-insensitive. We also accessed a better understanding of plant height genetic determinism with the combined effects of phenology-dependent (Ma6) and independent (qHT7.1 and Dw3) genomic regions. Therefore, we argue that the West and Central Africa-back-cross nested association mapping and the presented analytical approach constitute unique resources to better understand adaptation in sorghum with direct application to develop climate-smart varieties.

Identifying QTLs involved in hybrid performance and heterotic group complementarity: new GWAS models applied to factorial and admixed diallel maize hybrid panels

KEY MESSAGE

An original GWAS model integrating the ancestry of alleles was proposed and allowed the detection of background specific additive and dominance QTLs involved in heterotic group complementarity and hybrid performance. Maize genetic diversity is structured into genetic groups selected and improved relative to each other. This process increases group complementarity and differentiation over time and ensures that the hybrids produced from inter-group crosses exhibit high performances and heterosis. To identify loci involved in hybrid performance and heterotic group complementarity, we introduced an original association study model that disentangles allelic effects from the heterotic group origin of the alleles and compared it with a conventional additive/dominance model. This new model was applied on a factorial between Dent and Flint lines and a diallel between Dent-Flint admixed lines with two different layers of analysis: within each environment and in a multiple-environment context. We identified several strong additive QTLs for all traits, including some well-known additive QTLs for flowering time (in the region of Vgt1/2 on chromosome 8). Yield trait displayed significant non-additive effects in the diallel panel. Most of the detected Yield QTLs exhibited overdominance or, more likely, pseudo-overdominance effects. Apparent overdominance at these QTLs contributed to a part of the genetic group complementarity. The comparison between environments revealed a higher stability of additive QTL effects than non-additive ones. Several QTLs showed variations of effects according to the local heterotic group origin. We also revealed large chromosomic regions that display genetic group origin effects. Altogether, our results illustrate how admixed panels combined with dedicated GWAS modeling allow the identification of new QTLs that could not be revealed by a classical hybrid panel analyzed with traditional modeling.

Genomic prediction for the Germplasm Enhancement of Maize project

The Germplasm Enhancement of Maize (GEM) project was initiated in 1993 as a cooperative effort of public- and private-sector maize (Zea mays L.) breeders to enhance the genetic diversity of the U.S. maize crop. The GEM project selects progeny lines with high topcross yield potential from crosses between elite temperate lines and exotic parents. The GEM project has released hundreds of useful breeding lines based on phenotypic selection within selfing generations and multienvironment yield evaluations of GEM line topcrosses to elite adapted testers. Developing genomic selection (GS) models for the GEM project may contribute to increases in the rate of genetic gain. Here we evaluated the prediction ability of GS models trained on 6 yr of topcross evaluations from the two GEM programs in Raleigh, NC, and Ames, IA, documenting prediction abilities ranging from 0.36 to 0.75 for grain yield and from 0.78 to 0.96 for grain moisture when models were cross-validated within program and heterotic group. Predicted genetic gain from GS ranged from 0.95 to 2.58 times the gain from phenotypic selection. Prediction ability across program and heterotic group was generally poorer than within groups. Based on observed genomic relationships between GEM breeding lines and their tropical ancestors, GS for either yield or moisture would reduce recovery of exotic germplasm only slightly. Using GS models trained within program, the GEM programs should be able to more effectively deliver on its mission to broaden the genetic base of U.S. germplasm.

An Integration of Linkage Mapping and GWAS Reveals the Key Genes for Ear Shank Length in Maize

Ear shank length (ESL) has significant effects on grain yield and kernel dehydration rate in maize. Herein, linkage mapping and genome-wide association study were combined to reveal the genetic architecture of maize ESL. Sixteen quantitative trait loci (QTL) were identified in the segregation population, among which five were repeatedly detected across multiple environments. Meanwhile, 23 single nucleotide polymorphisms were associated with the ESL in the association panel, of which four were located in the QTL identified by linkage mapping and were designated as the population-common loci. A total of 42 genes residing in the linkage disequilibrium regions of these common variants and 12 of them were responsive to ear shank elongation. Of the 12 genes, five encode leucine-rich repeat receptor-like protein kinases, proline-rich proteins, and cyclin11, respectively, which were previously shown to regulate cell division, expansion, and elongation. Gene-based association analyses revealed that the variant located in Cyclin11 promoter affected the ESL among different lines. Cyclin11 showed the highest expression in the ear shank 15 days after silking among diverse tissues of maize, suggesting its role in modulating ESL. Our study contributes to the understanding of the genetic mechanism underlying maize ESL and genetic modification of maize dehydration rate and kernel yield.

statgenMPP: an R package implementing an IBD-based mixed model approach for QTL mapping in a wide range of multi-parent populations

MOTIVATION

Multi-parent populations (MPPs) are popular for QTL mapping because they combine wide genetic diversity in parents with easy control of population structure, but a limited number of software tools for QTL mapping are specifically developed for general MPP designs.

RESULTS

We developed an R package called statgenMPP, adopting a unified identity-by-descent (IBD)-based mixed model approach for QTL analysis in MPPs. The package offers easy-to-use functionalities of IBD calculations, mixed model solutions and visualizations for QTL mapping in a wide range of MPP designs, including diallele, nested-association mapping populations, multi-parent advanced genetic inter-cross populations and other complicated MPPs with known crossing schemes.

AVAILABILITY AND IMPLEMENTATION

The R package statgenMPP is open-source and freely available on CRAN at https://CRAN.R-project.org/package=statgenMPP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Modeling allelic diversity of multiparent mapping populations affects detection of quantitative trait loci

The search for quantitative trait loci that explain complex traits such as yield and drought tolerance has been ongoing in all crops. Methods such as biparental quantitative trait loci mapping and genome-wide association studies each have their own advantages and limitations. Multiparent advanced generation intercross populations contain more recombination events and genetic diversity than biparental mapping populations and are better able to estimate effect sizes of rare alleles than association mapping populations. Here, we discuss the results of using a multiparent advanced generation intercross population of doubled haploid maize lines created from 16 diverse founders to perform quantitative trait loci mapping. We compare 3 models that assume bi-allelic, founder, and ancestral haplotype allelic states for quantitative trait loci. The 3 methods have differing power to detect quantitative trait loci for a variety of agronomic traits. Although the founder approach finds the most quantitative trait loci, all methods are able to find unique quantitative trait loci, suggesting that each model has advantages for traits with different genetic architectures. A closer look at a well-characterized flowering time quantitative trait loci, qDTA8, which contains vgt1, highlights the strengths and weaknesses of each method and suggests a potential epistatic interaction. Overall, our results reinforce the importance of considering different approaches to analyzing genotypic datasets, and shows the limitations of binary SNP data for identifying multiallelic quantitative trait loci.

Multiallelic models for QTL mapping in diverse polyploid populations

Abstract

Quantitative trait locus (QTL) analysis allows to identify regions responsible for a trait and to associate alleles with their effect on phenotypes. When using biallelic markers to find these QTL regions, two alleles per QTL are modelled. This assumption might be close to reality in specific biparental crosses but is unrealistic in situations where broader genetic diversity is studied. Diversity panels used in genome-wide association studies or multi-parental populations can easily harbour multiple QTL alleles at each locus, more so in the case of polyploids that carry more than two alleles per individual. In such situations a multiallelic model would be closer to reality, allowing for different genetic effects for each potential allele in the population. To obtain such multiallelic markers we propose the usage of haplotypes, concatenations of nearby SNPs. We developed “mpQTL” an R package that can perform a QTL analysis at any ploidy level under biallelic and multiallelic models, depending on the marker type given. We tested the effect of genetic diversity on the power and accuracy difference between bi-allelic and multiallelic models using a set of simulated multiparental autotetraploid, outbreeding populations. Multiallelic models had higher detection power and were more precise than biallelic, SNP-based models, particularly when genetic diversity was higher. This confirms that moving to multi-allelic QTL models can lead to improved detection and characterization of QTLs.

Key message

QTL detection in populations with more than two functional QTL alleles (which is likely in multiparental and/or polyploid populations) is more powerful when using multiallelic models, rather than biallelic models.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04607-z.

A new Bayesian approach for QTL mapping of family data

In this paper, we propose a new Bayesian approach for QTL mapping of family data. The main purpose is to model a phenotype as a function of QTLs' effects. The model considers the detailed familiar dependence and it does not rely on random effects. It combines the probability for Mendelian inheritance of parents' genotype and the correlation between flanking markers and QTLs. This is an advance when compared with models which use only Mendelian segregation or only the correlation between markers and QTLs to estimate transmission probabilities. We use the Bayesian approach to estimate the number of QTLs, their location and the additive and dominance effects. We compare the performance of the proposed method with variance component and LASSO models using simulated and GAW17 data sets. Under tested conditions, the proposed method outperforms other methods in aspects such as estimating the number of QTLs, the accuracy of the QTLs' position and the estimate of their effects. The results of the application of the proposed method to data sets exceeded all of our expectations.

An IBD-based mixed model approach for QTL mapping in multiparental populations

The identity-by-descent (IBD)-based mixed model approach introduced in this study can detect quantitative trait loci (QTLs) referring to the parental origin and simultaneously account for multilevel relatedness of individuals within and across families. This unified approach is proved to be a powerful approach for all kinds of multiparental population (MPP) designs. Multiparental populations (MPPs) have become popular for quantitative trait loci (QTL) detection. Tools for QTL mapping in MPPs are mostly developed for specific MPPs and do not generalize well to other MPPs. We present an IBD-based mixed model approach for QTL mapping in all kinds of MPP designs, e.g., diallel, Nested Association Mapping (NAM), and Multiparental Advanced Generation Intercross (MAGIC) designs. The first step is to compute identity-by-descent (IBD) probabilities using a general Hidden Markov model framework, called reconstructing ancestry blocks bit by bit (RABBIT). Next, functions of IBD information are used as design matrices, or genetic predictors, in a mixed model approach to estimate variance components for multiallelic genetic effects associated with parents. Family-specific residual genetic effects are added, and a polygenic effect is structured by kinship relations between individuals. Case studies of simulated diallel, NAM, and MAGIC designs proved that the advanced IBD-based multi-QTL mixed model approach incorporating both kinship relations and family-specific residual variances (IBD.MQMkin_F) is robust across a variety of MPP designs and allele segregation patterns in comparison to a widely used benchmark association mapping method, and in most cases, outperformed or behaved at least as well as other tools developed for specific MPP designs in terms of mapping power and resolution. Successful analyses of real data cases confirmed the wide applicability of our IBD-based mixed model methodology.

Joint analysis of days to flowering reveals independent temperate adaptations in maize

Domesticates are an excellent model for understanding biological consequences of rapid climate change. Maize (Zea mays ssp. mays) was domesticated from a tropical grass yet is widespread across temperate regions today. We investigate the biological basis of temperate adaptation in diverse structured nested association mapping (NAM) populations from China, Europe (Dent and Flint) and the United States as well as in the Ames inbred diversity panel, using days to flowering as a proxy. Using cross-population prediction, where high prediction accuracy derives from overall genomic relatedness, shared genetic architecture, and sufficient diversity in the training population, we identify patterns in predictive ability across the five populations. To identify the source of temperate adapted alleles in these populations, we predict top associated genome-wide association study (GWAS) identified loci in a Random Forest Classifier using independent temperate-tropical North American populations based on lines selected from Hapmap3 as predictors. We find that North American populations are well predicted (AUC equals 0.89 and 0.85 for Ames and USNAM, respectively), European populations somewhat well predicted (AUC equals 0.59 and 0.67 for the Dent and Flint panels, respectively) and that the Chinese population is not predicted well at all (AUC is 0.47), suggesting an independent adaptation process for early flowering in China. Multiple adaptations for the complex trait days to flowering in maize provide hope for similar natural systems under climate change.

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.

The influence of QTL allelic diversity on QTL detection in multi-parent populations: a simulation study in sugar beet

BACKGROUND

Multi-parent populations (MPPs) are important resources for studying plant genetic architecture and detecting quantitative trait loci (QTLs). In MPPs, the QTL effects can show various levels of allelic diversity, which can be an important factor influencing the detection of QTLs. In MPPs, the allelic effects can be more or less specific. They can depend on an ancestor, a parent or the combination of parents in a cross. In this paper, we evaluated the effect of QTL allelic diversity on the QTL detection power in MPPs.

RESULTS

We simulated: a) cross-specific QTLs; b) parental and ancestral QTLs; and c) bi-allelic QTLs. Inspired by a real application in sugar beet, we tested different MPP designs (diallel, chessboard, factorial, and NAM) derived from five or nine parents to explore the ability to sample genetic diversity and detect QTLs. Using a fixed total population size, the QTL detection power was larger in MPPs with fewer but larger crosses derived from a reduced number of parents. The use of a larger set of parents was useful to detect rare alleles with a large phenotypic effect. The benefit of using a larger set of parents was however conditioned on an increase of the total population size. We also determined empirical confidence intervals for QTL location to compare the resolution of different designs. For QTLs representing 6% of the phenotypic variation, using 1600 F₂ offspring individuals, we found average 95% confidence intervals over different designs of 49 and 25 cM for cross-specific and bi-allelic QTLs, respectively.

CONCLUSIONS

MPPs derived from less parents with few but large crosses generally increased the QTL detection power. Using a larger set of parents to cover a wider genetic diversity can be useful to detect QTLs with a reduced minor allele frequency when the QTL effect is large and when the total population size is increased.

A combination of linkage mapping and GWAS brings new elements on the genetic basis of yield-related traits in maize across multiple environments

Using GWAS and QTL mapping identified 100 QTL and 138 SNPs, which control yield-related traits in maize. The candidate gene GRMZM2G098557 was further validated to regulate ear row number by using a segregation population. Understanding the genetic basis of yield-related traits contributes to the improvement of grain yield in maize. This study used an inter-mated B73 × Mo17 (IBM) Syn10 doubled-haploid (DH) population and an association panel to identify the genetic loci responsible for nine yield-related traits in maize. Using quantitative trait loci (QTL) mapping, 100 QTL influencing these traits were detected across different environments in the IBM Syn10 DH population, with 25 co-detected in multiple environments. Using a genome-wide association study (GWAS), 138 single-nucleotide polymorphisms (SNPs) were identified as correlated with these traits (P < 2.04E-06) in the association panel. Twenty-one pleiotropic QTL/SNPs were identified to control different traits in both populations. A combination of QTL mapping and GWAS uncovered eight significant SNPs (PZE-101097575, PZE-103169263, ZM011204-0763, PZE-104044017, PZE-104123110, SYN8062, PZE-108060911, and PZE-102043341) that were co-located within seven QTL confidence intervals. According to the eight co-localized SNPs by the two populations, 52 candidate genes were identified, among which the ear row number (ERN)-associated SNP SYN8062 was closely linked to SBP-transcription factor 7 (GRMZM2G098557). Several SBP-transcription factors were previously demonstrated to modulate maize ERN. We then validated the phenotypic effects of SYN8062 in the IBM Syn10 DH population, indicating that the ERN of the lines with the A-allele in SYN8062 was significantly (P < 0.05) larger than that of the lines with the G-allele in SYN8062 in each environment. These findings provide valuable information for understanding the genetic mechanisms of maize grain yield formation and for improving molecular marker-assisted selection for the high-yield breeding of maize.

Accounting for Group-Specific Allele Effects and Admixture in Genomic Predictions: Theory and Experimental Evaluation in Maize

Populations structured into genetic groups may display group-specific linkage disequilibrium, mutations, and/or interactions between quantitative trait loci and the genetic background. These factors lead to heterogeneous marker effects affecting the efficiency of genomic prediction, especially for admixed individuals. Such individuals have a genome that is a mosaic of chromosome blocks from different origins, and may be of interest to combine favorable group-specific characteristics. We developed two genomic prediction models adapted to the prediction of admixed individuals in presence of heterogeneous marker effects: multigroup admixed genomic best linear unbiased prediction random individual (MAGBLUP-RI), modeling the ancestry of alleles; and multigroup admixed genomic best linear unbiased prediction random allele effect (MAGBLUP-RAE), modeling group-specific distributions of allele effects. MAGBLUP-RI can estimate the segregation variance generated by admixture while MAGBLUP-RAE can disentangle the variability that is due to main allele effects from the variability that is due to group-specific deviation allele effects. Both models were evaluated for their genomic prediction accuracy using a maize panel including lines from the Dent and Flint groups, along with admixed individuals. Based on simulated traits, both models proved their efficiency to improve genomic prediction accuracy compared to standard GBLUP models. For real traits, a clear gain was observed at low marker densities whereas it became limited at high marker densities. The interest of including admixed individuals in multigroup training sets was confirmed using simulated traits, but was variable using real traits. Both MAGBLUP models and admixed individuals are of interest whenever group-specific SNP allele effects exist.

Multi-parent multi-environment QTL analysis: an illustration with the EU-NAM Flint population

Multi-parent populations multi-environment QTL experiments data should be analysed jointly to estimate the QTL effect variation within the population and between environments. Commonly, QTL detection in multi-parent populations (MPPs) data measured in multiple environments (ME) is done by analyzing genotypic values 'averaged' across environments. This method ignores the environment-specific QTL (QTLxE) effects. Running separate single environment analyses is a possibility to measure QTLxE effects, but those analyses do not model the genetic covariance due to the use of the same genotype in different environments. In this paper, we propose methods to analyse MPP-ME QTL experiments using simultaneously the data from several environments and modelling the genotypic covariance. Using data from the EU-NAM Flint population, we show that these methods estimate the QTLxE effects and that they can improve the quality of the QTL detection. Those methods also have a larger inference power. For example, they can be extended to integrate environmental indices like temperature or precipitation to better understand the mechanisms behind the QTLxE effects. Therefore, our methodology allows the exploitation of the full MPP-ME data potential: to estimate QTL effect variation (a) within the MPP between sub-populations due to different genetic backgrounds and (b) between environments.

Maize adaptation across temperate climates was obtained via expression of two florigen genes

Expansion of the maize growing area was central for food security in temperate regions. In addition to the suppression of the short-day requirement for floral induction, it required breeding for a large range of flowering time that compensates the effect of South-North gradients of temperatures. Here we show the role of a novel florigen gene, ZCN12, in the latter adaptation in cooperation with ZCN8. Strong eQTLs of ZCN8 and ZCN12, measured in 327 maize lines, accounted for most of the genetic variance of flowering time in platform and field experiments. ZCN12 had a strong effect on flowering time of transgenic Arabidopsis thaliana plants; a path analysis showed that it directly affected maize flowering time together with ZCN8. The allelic composition at ZCN QTLs showed clear signs of selection by breeders. This suggests that florigens played a central role in ensuring a large range of flowering time, necessary for adaptation to temperate areas.

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.

Optimized breeding strategies to harness genetic resources with different performance levels

Background

The narrow genetic base of elite germplasm compromises long-term genetic gain and increases the vulnerability to biotic and abiotic stresses in unpredictable environmental conditions. Therefore, an efficient strategy is required to broaden the genetic base of commercial breeding programs while not compromising short-term variety release. Optimal cross selection aims at identifying the optimal set of crosses that balances the expected genetic value and diversity. We propose to consider genomic selection and optimal cross selection to recurrently improve genetic resources (i.e. pre-breeding), to bridge the improved genetic resources with elites (i.e. bridging), and to manage introductions into the elite breeding population. Optimal cross selection is particularly adapted to jointly identify bridging, introduction and elite crosses to ensure an overall consistency of the genetic base broadening strategy.

Results

We compared simulated breeding programs introducing donors with different performance levels, directly or indirectly after bridging. We also evaluated the effect of the training set composition on the success of introductions. We observed that with recurrent introductions of improved donors, it is possible to maintain the genetic diversity and increase mid- and long-term performances with only limited penalty at short-term. Considering a bridging step yielded significantly higher mid- and long-term genetic gain when introducing low performing donors. The results also suggested to consider marker effects estimated with a broad training population including donor by elite and elite by elite progeny to identify bridging, introduction and elite crosses.

Conclusion

Results of this study provide guidelines on how to harness polygenic variation present in genetic resources to broaden elite germplasm.

Disentangling group specific QTL allele effects from genetic background epistasis using admixed individuals in GWAS: An application to maize flowering

When handling a structured population in association mapping, group-specific allele effects may be observed at quantitative trait loci (QTLs) for several reasons: (i) a different linkage disequilibrium (LD) between SNPs and QTLs across groups, (ii) group-specific genetic mutations in QTL regions, and/or (iii) epistatic interactions between QTLs and other loci that have differentiated allele frequencies between groups. We present here a new genome-wide association (GWAS) approach to identify QTLs exhibiting such group-specific allele effects. We developed genetic materials including admixed progeny from different genetic groups with known genome-wide ancestries (local admixture). A dedicated statistical methodology was developed to analyze pure and admixed individuals jointly, allowing one to disentangle the factors causing the heterogeneity of allele effects across groups. This approach was applied to maize by developing an inbred "Flint-Dent" panel including admixed individuals that was evaluated for flowering time. Several associations were detected revealing a wide range of configurations of allele effects, both at known flowering QTLs (Vgt1, Vgt2 and Vgt3) and new loci. We found several QTLs whose effect depended on the group ancestry of alleles while others interacted with the genetic background. Our GWAS approach provides useful information on the stability of QTL effects across genetic groups and can be applied to a wide range of species.

Genomic prediction with a maize collaborative panel: identification of genetic resources to enrich elite breeding programs

Collaborative diversity panels and genomic prediction seem relevant to identify and harness genetic resources for polygenic trait-specific enrichment of elite germplasms. In plant breeding, genetic diversity is important to maintain the pace of genetic gain and the ability to respond to new challenges in a context of climatic and social expectation changes. Many genetic resources are accessible to breeders but cannot all be considered for broadening the genetic diversity of elite germplasm. This study presents the use of genomic predictions trained on a collaborative diversity panel, which assembles genetic resources and elite lines, to identify resources to enrich an elite germplasm. A maize collaborative panel (386 lines) was considered to estimate genome-wide marker effects. Relevant predictive abilities (0.40-0.55) were observed on a large population of private elite materials, which supported the interest of such a collaborative panel for diversity management perspectives. Grain-yield estimated marker effects were used to select a donor that best complements an elite recipient at individual loci or haplotype segments, or that is expected to give the best-performing progeny with the elite. Among existing and new criteria that were compared, some gave more weight to the donor-elite complementarity than to the donor value, and appeared more adapted to long-term objective. We extended this approach to the selection of a set of donors complementing an elite population. We defined a crossing plan between identified donors and elite recipients. Our results illustrated how collaborative projects based on diversity panels including both public resources and elite germplasm can contribute to a better characterization of genetic resources in view of their use to enrich elite germplasm.

Improving Short- and Long-Term Genetic Gain by Accounting for Within-Family Variance in Optimal Cross-Selection

The implementation of genomic selection in recurrent breeding programs raises the concern that a higher inbreeding rate could compromise the long-term genetic gain. An optimized mating strategy that maximizes the performance in progeny and maintains diversity for long-term genetic gain is therefore essential. The optimal cross-selection approach aims at identifying the optimal set of crosses that maximizes the expected genetic value in the progeny under a constraint on genetic diversity in the progeny. Optimal cross-selection usually does not account for within-family selection, i.e., the fact that only a selected fraction of each family is used as parents of the next generation. In this study, we consider within-family variance accounting for linkage disequilibrium between quantitative trait loci to predict the expected mean performance and the expected genetic diversity in the selected progeny of a set of crosses. These predictions rely on the usefulness criterion parental contribution (UCPC) method. We compared UCPC-based optimal cross-selection and the optimal cross-selection approach in a long-term simulated recurrent genomic selection breeding program considering overlapping generations. UCPC-based optimal cross-selection proved to be more efficient to convert the genetic diversity into short- and long-term genetic gains than optimal cross-selection. We also showed that, using the UCPC-based optimal cross-selection, the long-term genetic gain can be increased with only a limited reduction of the short-term commercial genetic gain.

Multiple genetic loci affect place learning and memory performance in Drosophila melanogaster

Learning and memory are critical functions for all animals, giving individuals the ability to respond to changes in their environment. Within populations, individuals vary, however the mechanisms underlying this variation in performance are largely unknown. Thus, it remains to be determined what genetic factors cause an individual to have high learning ability and what factors determine how well an individual will remember what they have learned. To genetically dissect learning and memory performance, we used the Drosophila synthetic population resource (DSPR), a multiparent mapping resource in the model system Drosophila melanogaster, consisting of a large set of recombinant inbred lines (RILs) that naturally vary in these and other traits. Fruit flies can be trained in a "heat box" to learn to remain on one side of a chamber (place learning) and can remember this (place memory) over short timescales. Using this paradigm, we measured place learning and memory for ~49 000 individual flies from over 700 DSPR RILs. We identified 16 different loci across the genome that significantly affect place learning and/or memory performance, with 5 of these loci affecting both traits. To identify transcriptomic differences associated with performance, we performed RNA-Seq on pooled samples of seven high performing and seven low performing RILs for both learning and memory and identified hundreds of genes with differences in expression in the two sets. Integrating our transcriptomic results with the mapping results allowed us to identify nine promising candidate genes, advancing our understanding of the genetic basis underlying natural variation in learning and memory performance.

Construction of Genetic Linkage Maps in Multiparental Populations

Construction of genetic linkage maps has become a routine step for mapping quantitative trait loci (QTL), particularly in animal and plant breeding populations. Many multiparental populations have recently been produced to increase genetic diversity and QTL mapping resolution. However, few software packages are available for map construction in these populations. In this paper, we build a general framework for the construction of genetic linkage maps from genotypic data in diploid populations, including bi- and multiparental populations, cross-pollinated (CP) populations, and breeding pedigrees. The framework is implemented as an automatic pipeline called magicMap, where the maximum multilocus likelihood approach utilizes genotypic information efficiently. We evaluate magicMap by extensive simulations and eight real datasets: one biparental, one CP, four multiparent advanced generation intercross (MAGIC), and two nested association mapping (NAM) populations, the number of markers ranging from a few hundred to tens of thousands. Not only is magicMap the only software capable of accommodating all of these designs, it is more accurate and robust to missing genotypes and genotyping errors than commonly used packages.

Usefulness Criterion and Post-selection Parental Contributions in Multi-parental Crosses: Application to Polygenic Trait Introgression

Predicting the usefulness of crosses in terms of expected genetic gain and genetic diversity is of interest to secure performance in the progeny and to maintain long-term genetic gain in plant breeding. A wide range of crossing schemes are possible including large biparental crosses, backcrosses, four-way crosses, and synthetic populations. In silico progeny simulations together with genome-based prediction of quantitative traits can be used to guide mating decisions. However, the large number of multi-parental combinations can hinder the use of simulations in practice. Analytical solutions have been proposed recently to predict the distribution of a quantitative trait in the progeny of biparental crosses using information of recombination frequency and linkage disequilibrium between loci. Here, we extend this approach to obtain the progeny distribution of more complex crosses including two to four parents. Considering agronomic traits and parental genome contribution as jointly multivariate normally distributed traits, the usefulness criterion parental contribution (UCPC) enables to (i) evaluate the expected genetic gain for agronomic traits, and at the same time (ii) evaluate parental genome contributions to the selected fraction of progeny. We validate and illustrate UCPC in the context of multiple allele introgression from a donor into one or several elite recipients in maize (Zea mays L.). Recommendations regarding the interest of two-way, three-way, and backcrosses were derived depending on the donor performance. We believe that the computationally efficient UCPC approach can be useful for mate selection and allocation in many plant and animal breeding contexts.

Unravelling the genetic basis of Fusarium seedling rot resistance in the MAGIC maize population: novel targets for breeding

Fungal infection by Fusarium verticillioides is cause of prevalent maize disease leading to substantial reductions in yield and grain quality worldwide. Maize resistance to the fungus may occur at different developmental stages, from seedling to maturity. The breeding of resistant maize genotypes may take advantage of the identification of quantitative trait loci (QTL) responsible for disease resistance already commenced at seedling level. The Multi-parent Advance Generation Intercross (MAGIC) population was used to conduct high-definition QTL mapping for Fusarium seedling rot (FSR) resistance using rolled towel assay. Infection severity level, seedling weight and length were measured on 401 MAGIC maize recombinant inbred lines (RILs). QTL mapping was performed on reconstructed RIL haplotypes. One-fifth of the MAGIC RILs were resistant to FSR and 10 QTL were identified. For FSR, two QTL were detected at 2.8 Mb and 241.8 Mb on chromosome 4, and one QTL at 169.6 Mb on chromosome 5. Transcriptomic and sequencing information generated on the MAGIC founder lines was used to guide the identification of eight candidate genes within the identified FSR QTL. We conclude that the rolled towel assay applied to the MAGIC maize population provides a fast and cost-effective method to identify QTL and candidate genes for early resistance to F. verticillioides in maize.

Heterotic patterns of primary and secondary metabolites in the oilseed crop Brassica juncea

Heterosis refers to the superior performance of F1 hybrids over their respective parental inbred lines. Although the genetic and expression basis of heterosis have been previously investigated, the metabolic basis for this phenomenon is poorly understood. In a preliminary morphological study in Brassica juncea, we observed significant heterosis at the 50% flowering stage, wherein both the growth and reproduction of F1 reciprocal hybrids were greater than that of their parents. To identify the possible metabolic causes or consequences of this heterosis, we carried out targeted LC-MS analysis of 48 primary (amino acids and sugars) and secondary metabolites (phytohormones, glucosinolates, flavonoids, and phenolic esters) in five developmental tissues at 50% flowering in hybrids and inbred parents. Principal component analysis (PCA) of metabolites clearly separated inbred lines from their hybrids, particularly in the bud tissues. In general, secondary metabolites displayed more negative heterosis values in comparison to primary metabolites. The tested primary and secondary metabolites displayed both additive and non-additive modes of inheritance in F1 hybrids, wherein the number of metabolites showing an additive mode of inheritance were higher in buds and siliques (52.77-97.14%) compared to leaf tissues (47.37-80%). Partial least regression (PLS) analysis further showed that primary metabolites, in general, displayed higher association with morphological parameters in F1 hybrids. Overall, our results are consistent with a resource-cost model for heterosis in B. juncea, where metabolite allocation in hybrids appears to favor growth, at the expense of secondary metabolism.

Technological advances in maize breeding: past, present and future

Maize has for many decades been both one of the most important crops worldwide and one of the primary genetic model organisms. More recently, maize breeding has been impacted by rapid technological advances in sequencing and genotyping technology, transformation including genome editing, doubled haploid technology, parallelled by progress in data sciences and the development of novel breeding approaches utilizing genomic information. Herein, we report on past, current and future developments relevant for maize breeding with regard to (1) genome analysis, (2) germplasm diversity characterization and utilization, (3) manipulation of genetic diversity by transformation and genome editing, (4) inbred line development and hybrid seed production, (5) understanding and prediction of hybrid performance, (6) breeding methodology and (7) synthesis of opportunities and challenges for future maize breeding.

Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism

Plants have served as sources providing humans with metabolites for food and nutrition, biomaterials for living, and treatment for pain and disease. Plants produce a huge array of metabolites, with an immense diversity at both the population and individual levels. Dissection of the genetic bases for metabolic diversity has attracted increasing research attention. The concept of genome-wide association study (GWAS) was extended to studies on the diversity of plant metabolome that benefitted from the development of mass-spectrometry-based analytical systems and genome sequencing technologies. Metabolic genome-wide association study (mGWAS) is one of the most powerful tools for global identification of genetic determinants for diversity of plant metabolism. Recently, mGWAS has been performed for various species with continuous improvements, providing deeper insights into the genetic bases of metabolic diversity. In this review, we discuss fully the achievements to date and remaining challenges that are associated with both mGWAS and mGWAS-based multi-dimensional analysis. We begin with a summary of GWAS and its development based on statistical methods and populations. As variation in targeted traits is essential for GWAS, we review metabolic diversity and its rise at both the population and individual levels. Subsequently, the application of mGWAS for plants and its corresponding achievements are fully discussed. We address the current knowledge on mGWAS-based multi-dimensional analysis and emerging insights into the diversity of metabolism.

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.

Genetic Loci Governing Androgenic Capacity in Perennial Ryegrass (Lolium perenne L.)

Immature pollen can be induced to switch developmental pathways from gametogenesis to embryogenesis and subsequently regenerate into homozygous, diploid plants. Such androgenic production of doubled haploids is particularly useful for species where inbreeding is hampered by effective self-incompatibility systems. Therefore, increasing the generally low androgenic capacity of perennial ryegrass (Lolium perenne L.) germplasm would enable the efficient production of homozygous plant material, so that a more effective exploitation of heterosis through hybrid breeding schemes can be realized. Here, we present the results of a genome-wide association study in a heterozygous, multiparental population of perennial ryegrass (n = 391) segregating for androgenic capacity. Genotyping-by-sequencing was used to interrogate gene- dense genomic regions and revealed over 1,100 polymorphic sites. Between one and 10 quantitative trait loci (QTL) were identified for anther response, embryo and total plant production, green and albino plant production and regeneration. Most traits were under polygenic control, although a major QTL on linkage group 5 was associated with green plant regeneration. Distinct genetic factors seem to affect green and albino plant recovery. Two intriguing candidate genes, encoding chromatin binding domains of the developmental phase transition regulator, Polycomb Repressive Complex 2, were identified. Our results shed the first light on the molecular mechanisms behind perennial ryegrass microspore embryogenesis and enable marker-assisted introgression of androgenic capacity into recalcitrant germplasm of this forage crop of global significance.

Phenomics allows identification of genomic regions affecting maize stomatal conductance with conditional effects of water deficit and evaporative demand

Stomatal conductance is central for the trades-off between hydraulics and photosynthesis. We aimed at deciphering its genetic control and that of its responses to evaporative demand and water deficit, a nearly impossible task with gas exchanges measurements. Whole-plant stomatal conductance was estimated via inversion of the Penman-Monteith equation from data of transpiration and plant architecture collected in a phenotyping platform. We have analysed jointly 4 experiments with contrasting environmental conditions imposed to a panel of 254 maize hybrids. Estimated whole-plant stomatal conductance closely correlated with gas-exchange measurements and biomass accumulation rate. Sixteen robust quantitative trait loci (QTLs) were identified by genome wide association studies and co-located with QTLs of transpiration and biomass. Light, vapour pressure deficit, or soil water potential largely accounted for the differences in allelic effects between experiments, thereby providing strong hypotheses for mechanisms of stomatal control and a way to select relevant candidate genes among the 1-19 genes harboured by QTLs. The combination of allelic effects, as affected by environmental conditions, accounted for the variability of stomatal conductance across a range of hybrids and environmental conditions. This approach may therefore contribute to genetic analysis and prediction of stomatal control in diverse environments.

Genomic Prediction Within and Across Biparental Families: Means and Variances of Prediction Accuracy and Usefulness of Deterministic Equations

A major application of genomic prediction (GP) in plant breeding is the identification of superior inbred lines within families derived from biparental crosses. When models for various traits were trained within related or unrelated biparental families (BPFs), experimental studies found substantial variation in prediction accuracy (PA), but little is known about the underlying factors. We used SNP marker genotypes of inbred lines from either elite germplasm or landraces of maize (Zeamays L.) as parents to generate in silico 300 BPFs of doubled-haploid lines. We analyzed PA within each BPF for 50 simulated polygenic traits, using genomic best linear unbiased prediction (GBLUP) models trained with individuals from either full-sib (FSF), half-sib (HSF), or unrelated families (URF) for various sizes (Ntrain) of the training set and different heritabilities (h2). In addition, we modified two deterministic equations for forecasting PA to account for inbreeding and genetic variance unexplained by the training set. Averaged across traits, PA was high within FSF (0.41–0.97) with large variation only for Ntrain<50 and h2<0.6. For HSF and URF, PA was on average ∼40–60% lower and varied substantially among different combinations of BPFs used for model training and prediction as well as different traits. As exemplified by HSF results, PA of across-family GP can be very low if causal variants not segregating in the training set account for a sizeable proportion of the genetic variance among predicted individuals. Deterministic equations accurately forecast the PA expected over many traits, yet cannot capture trait-specific deviations. We conclude that model training within BPFs generally yields stable PA, whereas a high level of uncertainty is encountered in across-family GP. Our study shows the extent of variation in PA that must be at least reckoned with in practice and offers a starting point for the design of training sets composed of multiple BPFs.

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.

Predicting genomic selection efficiency to optimize calibration set and to assess prediction accuracy in highly structured populations

KEY MESSAGE

We propose a criterion to predict genomic selection efficiency for structured populations. This criterion is useful to define optimal calibration set and to estimate prediction reliability for multiparental populations. Genomic selection refers to the use of genotypic information for predicting the performance of selection candidates. It has been shown that prediction accuracy depends on various parameters including the composition of the calibration set (CS). Assessing the level of accuracy of a given prediction scenario is of highest importance because it can be used to optimize CS sampling before collecting phenotypes, and once the breeding values are predicted it informs the breeders about the reliability of these predictions. Different criteria were proposed to optimize CS sampling in highly diverse panels, which can be useful to screen collections of genotypes. But plant breeders often work on structured material such as biparental or multiparental populations, for which these criteria are less adapted. We derived from the generalized coefficient of determination (CD) theory different criteria to optimize CS sampling and to assess the reliability associated to predictions in structured populations. These criteria were evaluated on two nested association mapping (NAM) populations and two highly diverse panels of maize. They were efficient to sample optimized CS in most situations. They could also estimate at least partly the reliability associated to predictions between NAM families, but they could not estimate differences in the reliability associated to the predictions of NAM families using the highly diverse panels as calibration sets. We illustrated that the CD criteria could be adapted to various prediction scenarios including inter and intra-family predictions, resulting in higher prediction accuracies.

Genetic Gain Increases by Applying the Usefulness Criterion with Improved Variance Prediction in Selection of Crosses

A crucial step in plant breeding is the selection and combination of parents to form new crosses. Genome-based prediction guides the selection of high-performing parental lines in many crop breeding programs which ensures a high mean performance of progeny. To warrant maximum selection progress, a new cross should also provide a large progeny variance. The usefulness concept as measure of the gain that can be obtained from a specific cross accounts for variation in progeny variance. Here, it is shown that genetic gain can be considerably increased when crosses are selected based on their genomic usefulness criterion compared to selection based on mean genomic estimated breeding values. An efficient and improved method to predict the genetic variance of a cross based on Markov chain Monte Carlo samples of marker effects from a whole-genome regression model is suggested. In simulations representing selection procedures in crop breeding programs, the performance of this novel approach is compared with existing methods, like selection based on mean genomic estimated breeding values and optimal haploid values. In all cases, higher genetic gain was obtained compared with previously suggested methods. When 1% of progenies per cross were selected, the genetic gain based on the estimated usefulness criterion increased by 0.14 genetic standard deviation compared to a selection based on mean genomic estimated breeding values. Analytical derivations of the progeny genotypic variance-covariance matrix based on parental genotypes and genetic map information make simulations of progeny dispensable, and allow fast implementation in large-scale breeding programs.

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-based-breeding tools for tropical maize improvement

Maize has traditionally been the main staple diet in the Southern Asia and Sub-Saharan Africa and widely grown by millions of resource poor small scale farmers. Approximately, 35.4 million hectares are sown to tropical maize, constituting around 59% of the developing worlds. Tropical maize encounters tremendous challenges besides poor agro-climatic situations with average yields recorded <3 tones/hectare that is far less than the average of developed countries. On the contrary to poor yields, the demand for maize as food, feed, and fuel is continuously increasing in these regions. Heterosis breeding introduced in early 90 s improved maize yields significantly, but genetic gains is still a mirage, particularly for crop growing under marginal environments. Application of molecular markers has accelerated the pace of maize breeding to some extent. The availability of array of sequencing and genotyping technologies offers unrivalled service to improve precision in maize-breeding programs through modern approaches such as genomic selection, genome-wide association studies, bulk segregant analysis-based sequencing approaches, etc. Superior alleles underlying complex traits can easily be identified and introgressed efficiently using these sequence-based approaches. Integration of genomic tools and techniques with advanced genetic resources such as nested association mapping and backcross nested association mapping could certainly address the genetic issues in maize improvement programs in developing countries. Huge diversity in tropical maize and its inherent capacity for doubled haploid technology offers advantage to apply the next generation genomic tools for accelerating production in marginal environments of tropical and subtropical world. Precision in phenotyping is the key for success of any molecular-breeding approach. This article reviews genomic technologies and their application to improve agronomic traits in tropical maize breeding has been reviewed in detail.

How do the type of QTL effect and the form of the residual term influence QTL detection in multi-parent populations? A case study in the maize EU-NAM population

In the QTL analysis of multi-parent populations, the inclusion of QTLs with various types of effects can lead to a better description of the phenotypic variation and increased power. For the type of QTL effect in QTL models for multi-parent populations (MPPs), various options exist to define them with respect to their origin. They can be modelled as referring to close parental lines or to further away ancestral founder lines. QTL models for MPPs can also be characterized by the homo- or heterogeneity of variance for polygenic effects. The most suitable model for the origin of the QTL effect and the homo- or heterogeneity of polygenic effects may be a function of the genetic distance distribution between the parents of MPPs. We investigated the statistical properties of various QTL detection models for MPPs taking into account the genetic distances between the parents of the MPP. We evaluated models with different assumptions about the QTL effect and the form of the residual term using cross validation. For the EU-NAM data, we showed that it can be useful to mix in the same model QTLs with different types of effects (parental, ancestral, or bi-allelic). The benefit of using cross-specific residual terms to handle the heterogeneity of variance was less obvious for this particular data set.

Increased Power To Dissect Adaptive Traits in Global Sorghum Diversity Using a Nested Association Mapping Population

Adaptation of domesticated species to diverse agroclimatic regions has led to abundant trait diversity. However, the resulting population structure and genetic heterogeneity confounds association mapping of adaptive traits. To address this challenge in sorghum [Sorghum bicolor (L.) Moench]-a widely adapted cereal crop-we developed a nested association mapping (NAM) population using 10 diverse global lines crossed with an elite reference line RTx430. We characterized the population of 2214 recombinant inbred lines at 90,000 SNPs using genotyping-by-sequencing. The population captures ∼70% of known global SNP variation in sorghum, and 57,411 recombination events. Notably, recombination events were four- to fivefold enriched in coding sequences and 5' untranslated regions of genes. To test the power of the NAM population for trait dissection, we conducted joint linkage mapping for two major adaptive traits, flowering time and plant height. We precisely mapped several known genes for these two traits, and identified several additional QTL. Considering all SNPs simultaneously, genetic variation accounted for 65% of flowering time variance and 75% of plant height variance. Further, we directly compared NAM to genome-wide association mapping (using panels of the same size) and found that flowering time and plant height QTL were more consistently identified with the NAM population. Finally, for simulated QTL under strong selection in diversity panels, the power of QTL detection was up to three times greater for NAM vs. association mapping with a diverse panel. These findings validate the NAM resource for trait mapping in sorghum, and demonstrate the value of NAM for dissection of adaptive traits.

Strong epistatic and additive effects of linked candidate SNPs for Drosophila pigmentation have implications for analysis of genome-wide association studies results

Background

The mapping resolution of genome-wide association studies (GWAS) is limited by historic recombination events and effects are often assigned to haplotype blocks rather than individual SNPs. It is not clear how many of the SNPs in the block, and which ones, are causative. Drosophila pigmentation is a powerful model to dissect the genetic basis of intra-specific and inter-specific phenotypic variation. Three tightly linked SNPs in the t-MSE enhancer have been identified in three D. melanogaster populations as major contributors to female abdominal pigmentation. This enhancer controls the expression of the pigmentation gene tan (t) in the abdominal epidermis. Two of the three SNPs were confirmed in an independent study using the D. melanogaster Genetic Reference Panel established from a North American population.

Results

We determined the functional impact of SNP1, SNP2, and SNP3 using transgenic lines to test all possible haplotypes in vivo. We show that all three candidate SNPs contribute to female Drosophila abdominal pigmentation. Interestingly, only two SNPs agree with the effect predicted by GWAS; the third one goes in the opposite direction because of linkage disequilibrium between multiple functional SNPs. Our experimental design uncovered strong additive effects for the three SNPs, but we also found significant epistatic effects explaining up to 11% of the total variation.

Conclusions

Our results suggest that linked causal variants are important for the interpretation of GWAS and functional validation is needed to understand the genetic architecture of traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1262-7) contains supplementary material, which is available to authorized users.

Loci Contributing to Boric Acid Toxicity in Two Reference Populations of Drosophila melanogaster

Populations maintain considerable segregating variation in the response to toxic, xenobiotic compounds. To identify variants associated with resistance to boric acid, a commonly-used household insecticide with a poorly understood mechanism of action, we assayed thousands of individuals from hundreds of strains. Using the Drosophila Synthetic Population Resource (DSPR), a multi-parental population (MPP) of inbred genotypes, we mapped six QTL to short genomic regions containing few protein-coding genes (3–188), allowing us to identify plausible candidate genes underlying resistance to boric acid toxicity. One interval contains multiple genes from the cytochrome P450 family, and we show that ubiquitous RNAi of one of these genes, Cyp9b2, markedly reduces resistance to the toxin. Resistance to boric acid is positively correlated with caffeine resistance. The two phenotypes additionally share a pair of QTL, potentially suggesting a degree of pleiotropy in the genetic control of resistance to these two distinct xenobiotics. Finally, we screened the Drosophila Genetic Reference Panel (DGRP) in an attempt to identify sequence variants within mapped QTL that are associated with boric acid resistance. The approach was largely unsuccessful, with only one QTL showing any associations at QTL-specific 20% False Discovery Rate (FDR) thresholds. Nonetheless, these associations point to a potential candidate gene that can be targeted in future validation efforts. Although the mapping data resulting from the two reference populations do not clearly overlap, our work provides a starting point for further genetic dissection of the processes underlying boric acid toxicity in insects.

High-throughput SNP genotyping of modern and wild emmer wheat for yield and root morphology using a combined association and linkage analysis

Durum wheat (Triticum turgidum var. durum Desf.) is a major world crop that is grown primarily in areas of the world that experience periodic drought, and therefore, breeding climate-resilient durum wheat is a priority. High-throughput single nucleotide polymorphism (SNP) genotyping techniques have greatly increased the power of linkage and association mapping analyses for bread wheat, but as yet there is no durum wheat-specific platform available. In this study, we evaluate the new 384HT Wheat Breeders Array for its usefulness in tetraploid wheat breeding by genotyping a breeding population of F₆ hybrids, derived from multiple crosses between T. durum cultivars and wild and cultivated emmer wheat accessions. Using a combined linkage and association mapping approach, we generated a genetic map including 1345 SNP markers, and identified markers linked to 6 QTLs for coleoptile length (2), heading date (1), anthocyanin accumulation (1) and osmotic stress tolerance (2). We also developed a straightforward approach for combining genetic data from multiple families of limited size that will be useful for evaluating and mapping pre-existing breeding material.

Genome-wide Association Studies in Maize: Praise and Stargaze

Genome-wide association study (GWAS) has become a widely accepted strategy for decoding genotype-phenotype associations in many species thanks to advances in next-generation sequencing (NGS) technologies. Maize is an ideal crop for GWAS and significant progress has been made in the last decade. This review summarizes current GWAS efforts in maize functional genomics research and discusses future prospects in the omics era. The general goal of GWAS is to link genotypic variations to corresponding differences in phenotype using the most appropriate statistical model in a given population. The current review also presents perspectives for optimizing GWAS design and analysis. GWAS analysis of data from RNA, protein, and metabolite-based omics studies is discussed, along with new models and new population designs that will identify causes of phenotypic variation that have been hidden to date. The joint and continuous efforts of the whole community will enhance our understanding of maize quantitative traits and boost crop molecular breeding designs.

Estimating parent-specific QTL effects through cumulating linked identity-by-state SNP effects in multiparental populations

The emergence of multiparental mapping populations enabled plant geneticists to gain deeper insights into the genetic architecture of major agronomic traits and to map quantitative trait loci (QTLs) controlling the expression of these traits. Although the investigated mapping populations are similar, one open question is whether genotype data should be modelled as identical by state (IBS) or identical by descent (IBD). Whereas IBS simply makes use of raw genotype scores to distinguish alleles, IBD data are derived from parental offspring information. We report on comparing IBS and IBD by applying two multiple regression models on four traits studied in the barley nested association mapping (NAM) population HEB-25. We observed that modelling parent-specific IBD genotypes produced a lower number of significant QTLs with increased prediction abilities compared with modelling IBS genotypes. However, at lower trait heritabilities the IBS model produced higher prediction abilities. We developed a method to estimate multiallelic QTL effects in multiparental populations from simple biallelic IBS data. This method is based on cumulating IBS-derived single-nucleotide polymorphism (SNP) effect estimates in a defined genetic region surrounding a QTL. Comparing the resulting parent-specific QTL effects with those obtained from IBD approaches revealed high accordance that could be confirmed through simulations. The method turned out to be also applicable to a barley multiparent advanced generation inter-cross (MAGIC) population. The 'cumulation method' represents a universal approach to differentiate parent-specific QTL effects in multiparental populations, even if no IBD information is available. In future, the method could further benefit from the availability of much denser SNP maps.

Genome-Wide Analysis of Yield in Europe: Allelic Effects Vary with Drought and Heat Scenarios

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.