Using crop growth model stress covariates and AMMI decomposition to better predict genotype-by-environment interactions.

KEY MESSAGE: We propose new methods to predict genotype × environment interaction by selecting relevant environmental covariates and using an AMMI decomposition of the interaction. Farmers are asked to produce more efficiently and to reduce their inputs in the context of climate change. They have to face more and more limiting factors that can combine in numerous stress scenarios. One solution to this challenge is to develop varieties adapted to specific environmental stress scenarios. For this, plant breeders can use genomic predictions coupled with environmental characterization to identify promising combinations of genes in relation to stress covariates. One way to do it is to take into account the genetic similarity between varieties and the similarity between environments within a mixed model framework. Molecular markers and environmental covariates (EC) can be used to estimate relevant covariance matrices. In the present study, based on a multi-environment trial of 220 European elite winter bread wheat (Triticum aestivum L.) varieties phenotyped in 42 environments, we compared reference regression models potentially including ECs, and proposed alternative models to increase prediction accuracy. We showed that selecting a subset of ECs, and estimating covariance matrices using an AMMI decomposition to benefit from the information brought by the phenotypic records of the training set are promising approaches to better predict genotype-by-environment interactions (G × E). We found that using a different kinship for the main genetic effect and the G × E effect increased prediction accuracy. Our study also demonstrates that integrative stress indexes simulated by crop growth models are more efficient to capture G × E than climatic covariates.

Studying the genetic basis of drought tolerance in sorghum by managed stress trials and adjustments for phenological and plant height differences.

Managed environments in the form of well watered and water stressed trials were performed to study the genetic basis of grain yield and stay green in sorghum with the objective of validating previously detected QTL. As variations in phenology and plant height may influence QTL detection for the target traits, QTL for flowering time and plant height were introduced as cofactors in QTL analyses for yield and stay green. All but one of the flowering time QTL were detected near yield and stay green QTL. Similar co-localization was observed for two plant height QTL. QTL analysis for yield, using flowering time/plant height cofactors, led to yield QTL on chromosomes 2, 3, 6, 8 and 10. For stay green, QTL on chromosomes 3, 4, 8 and 10 were not related to differences in flowering time/plant height. The physical positions for markers in QTL regions projected on the sorghum genome suggest that the previously detected plant height QTL, Sb-HT9-1, and Dw2, in addition to the maturity gene, Ma5, had a major confounding impact on the expression of yield and stay green QTL. Co-localization between an apparently novel stay green QTL and a yield QTL on chromosome 3 suggests there is potential for indirect selection based on stay green to improve drought tolerance in sorghum. Our QTL study was carried out with a moderately sized population and spanned a limited geographic range, but still the results strongly emphasize the necessity of corrections for phenology in QTL mapping for drought tolerance traits in sorghum.

A mixed model QTL analysis for sugarcane multiple-harvest-location trial data.

Sugarcane-breeding programs take at least 12 years to develop new commercial cultivars. Molecular markers offer a possibility to study the genetic architecture of quantitative traits in sugarcane, and they may be used in marker-assisted selection to speed up artificial selection. Although the performance of sugarcane progenies in breeding programs are commonly evaluated across a range of locations and harvest years, many of the QTL detection methods ignore two- and three-way interactions between QTL, harvest, and location. In this work, a strategy for QTL detection in multi-harvest-location trial data, based on interval mapping and mixed models, is proposed and applied to map QTL effects on a segregating progeny from a biparental cross of pre-commercial Brazilian cultivars, evaluated at two locations and three consecutive harvest years for cane yield (tonnes per hectare), sugar yield (tonnes per hectare), fiber percent, and sucrose content. In the mixed model, we have included appropriate (co)variance structures for modeling heterogeneity and correlation of genetic effects and non-genetic residual effects. Forty-six QTLs were found: 13 QTLs for cane yield, 14 for sugar yield, 11 for fiber percent, and 8 for sucrose content. In addition, QTL by harvest, QTL by location, and QTL by harvest by location interaction effects were significant for all evaluated traits (30 QTLs showed some interaction, and 16 none). Our results contribute to a better understanding of the genetic architecture of complex traits related to biomass production and sucrose content in sugarcane.

A mixed-model approach to association mapping using pedigree information with an illustration of resistance to Phytophthora infestans in potato.

Association or linkage disequilibrium (LD)-based mapping strategies are receiving increased attention for the identification of quantitative trait loci (QTL) in plants as an alternative to more traditional, purely linkage-based approaches. An attractive property of association approaches is that they do not require specially designed crosses between inbred parents, but can be applied to collections of genotypes with arbitrary and often unknown relationships between the genotypes. A less obvious additional attractive property is that association approaches offer possibilities for QTL identification in crops with hard to model segregation patterns. The availability of candidate genes and targeted marker systems facilitates association approaches, as will appropriate methods of analysis. We propose an association mapping approach based on mixed models with attention to the incorporation of the relationships between genotypes, whether induced by pedigree, population substructure, or otherwise. Furthermore, we emphasize the need to pay attention to the environmental features of the data as well, i.e., adequate representation of the relations among multiple observations on the same genotypes. We illustrate our modeling approach using 25 years of Dutch national variety list data on late blight resistance in the genetically complex crop of potato. As markers, we used nucleotide binding-site markers, a specific type of marker that targets resistance or resistance-analog genes. To assess the consistency of QTL identified by our mixed-model approach, a second independent data set was analyzed. Two markers were identified that are potentially useful in selection for late blight resistance in potato.

QTL methodology for response curves on the basis of non-linear mixed models, with an illustration to senescence in potato.

The improvement of quantitative traits in plant breeding will in general benefit from a better understanding of the genetic basis underlying their development. In this paper, a QTL mapping strategy is presented for modelling the development of phenotypic traits over time. Traditionally, crop growth models are used to study development. We propose an integration of crop growth models and QTL models within the framework of non-linear mixed models. We illustrate our approach with a QTL model for leaf senescence in a diploid potato cross. Assuming a logistic progression of senescence in time, two curve parameters are modelled, slope and inflection point, as a function of QTLs. The final QTL model for our example data contained four QTLs, of which two affected the position of the inflection point, one the senescence progression-rate, and a final one both inflection point and rate.