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

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.

Leveraging Multiomics Insights and Exploiting Wild Relatives' Potential for Drought and Heat Tolerance in Maize

Climate change, particularly drought and heat stress, may slash agricultural productivity by 25.7% by 2080, with maize being the hardest hit. Therefore, unraveling the molecular nature of plant responses to these stressors is vital for the development of climate-smart maize. This manuscript's primary objective was to examine how maize plants respond to these stresses, both individually and in combination. Additionally, the paper delved into harnessing the potential of maize wild relatives as a valuable genetic resource and leveraging AI-based technologies to boost maize resilience. The role of multiomics approaches particularly genomics and transcriptomics in dissecting the genetic basis of stress tolerance was also highlighted. The way forward was proposed to utilize a bunch of information obtained through omics technologies by an interdisciplinary state-of-the-art forward-looking big-data, cyberagriculture system, and AI-based approach to orchestrate the development of climate resilient maize genotypes.

High-dimensional multi-omics measured in controlled conditions are useful for maize platform and field trait predictions

KEY MESSAGE

Transcriptomics and proteomics information collected on a platform can predict additive and non-additive effects for platform traits and additive effects for field traits. The effects of climate change in the form of drought, heat stress, and irregular seasonal changes threaten global crop production. The ability of multi-omics data, such as transcripts and proteins, to reflect a plant's response to such climatic factors can be capitalized in prediction models to maximize crop improvement. Implementing multi-omics characterization in field evaluations is challenging due to high costs. It is, however, possible to do it on reference genotypes in controlled conditions. Using omics measured on a platform, we tested different multi-omics-based prediction approaches, using a high dimensional linear mixed model (MegaLMM) to predict genotypes for platform traits and agronomic field traits in a panel of 244 maize hybrids. We considered two prediction scenarios: in the first one, new hybrids are predicted (CV-NH), and in the second one, partially observed hybrids are predicted (CV-POH). For both scenarios, all hybrids were characterized for omics on the platform. We observed that omics can predict both additive and non-additive genetic effects for the platform traits, resulting in much higher predictive abilities than GBLUP. It highlights their efficiency in capturing regulatory processes in relation to growth conditions. For the field traits, we observed that the additive components of omics only slightly improved predictive abilities for predicting new hybrids (CV-NH, model MegaGAO) and for predicting partially observed hybrids (CV-POH, model GAOxW-BLUP) in comparison to GBLUP. We conclude that measuring the omics in the fields would be of considerable interest in predicting productivity if the costs of omics drop significantly.

Envirotyping within a multi-environment trial allowed identifying genetic determinants of winter oilseed rape yield stability

KEY MESSAGE

A comprehensive environmental characterization allowed identifying stable and interactive QTL for seed yield: QA09 and QC09a were detected across environments; whereas QA07a was specifically detected on the most stressed environments. A main challenge for rapeseed consists in maintaining seed yield while adapting to climate changes and contributing to environmental-friendly cropping systems. Breeding for cultivar adaptation is one of the keys to meet this challenge. Therefore, we propose to identify the genetic determinant of seed yield stability for winter oilseed rape using GWAS coupled with a multi-environmental trial and to interpret them in the light of environmental characteristics. Due to a comprehensive characterization of a multi-environmental trial using 79 indicators, four contrasting envirotypes were defined and used to identify interactive and stable seed yield QTL. A total of four QTLs were detected, among which, QA09 and QC09a, were stable (detected at the multi-environmental trial scale or for different envirotypes and environments); and one, QA07a, was specifically detected into the most stressed envirotype. The analysis of the molecular diversity at QA07a showed a lack of genetic diversity within modern lines compared to older cultivars bred before the selection for low glucosinolate content. The results were discussed in comparison with other studies and methods as well as in the context of breeding programs.

Primary, seminal and lateral roots of maize show type-specific growth and hydraulic responses to water deficit

The water uptake capacity of a root system is determined by its architecture and hydraulic properties, which together shape the root hydraulic architecture. Here, we investigated root responses to water deficit (WD) in seedlings of a maize (Zea mays) hybrid line (B73H) grown in hydroponic conditions, taking into account the primary root (PR), the seminal roots (SR), and their respective lateral roots. WD was induced by various polyethylene glycol concentrations and resulted in dose-dependent inhibitions of axial and lateral root growth, lateral root formation, and hydraulic conductivity (Lpr), with slightly distinct sensitivities to WD between PR and SR. Inhibition of Lpr by WD showed a half-time of 5 to 6 min and was fully (SR) or partially (PR) reversible within 40 min. In the two root types, WD resulted in reduced aquaporin expression and activity, as monitored by mRNA abundance of 13 plasma membrane intrinsic protein (ZmPIP) isoforms and inhibition of Lpr by sodium azide, respectively. An enhanced suberization/lignification of the epi- and exodermis was observed under WD in axial roots and in lateral roots of the PR but not in those of SR. Inverse modeling revealed a steep increase in axial conductance in root tips of PR and SR grown under WD that may be due to the decreased growth rate of axial roots in these conditions. Overall, our work reveals that these root types show quantitative differences in their anatomical, architectural, and hydraulic responses to WD, in terms of sensitivity, amplitude and reversibility. This distinct functionalization may contribute to integrative acclimation responses of whole root systems to soil WD.

Maize green leaf area index dynamics: genetic basis of a new secondary trait for grain yield in optimal and drought conditions

KEY MESSAGE

Green Leaf Area Index dynamics is a promising secondary trait for grain yield and drought tolerance. Multivariate GWAS is particularly well suited to identify the genetic determinants of the green leaf area index dynamics. Improvement of maize grain yield is impeded by important genotype-environment interactions, especially under drought conditions. The use of secondary traits, that are correlated with yield, more heritable and less prone to genotype-environment interactions, can increase breeding efficiency. Here, we studied the genetic basis of a new secondary trait: the green leaf area index (GLAI) dynamics over the maize life cycle. For this, we used an unmanned aerial vehicle to characterize the GLAI dynamics of a diverse panel in well-watered and water-deficient trials in two years. From the dynamics, we derived 24 traits (slopes, durations, areas under the curve), and showed that six of them were heritable traits representative of the panel diversity. To identify the genetic determinants of GLAI, we compared two genome-wide association approaches: a univariate (single-trait) method and a multivariate (multi-trait) method combining GLAI traits, grain yield, and precocity. The explicit modeling of correlation structure between secondary traits and grain yield in the multivariate mixed model led to 2.5 times more associations detected. A total of 475 quantitative trait loci (QTLs) were detected. The genetic architecture of GLAI traits appears less complex than that of yield with stronger-effect QTLs that are more stable between environments. We also showed that a subset of GLAI QTLs explains nearly one fifth of yield variability across a larger environmental network of 11 water-deficient trials. GLAI dynamics is a promising grain yield secondary trait in optimal and drought conditions, and the detected QTLs could help to increase breeding efficiency through a marker-assisted approach.

A review of artificial intelligence-assisted omics techniques in plant defense: current trends and future directions

Plants intricately deploy defense systems to counter diverse biotic and abiotic stresses. Omics technologies, spanning genomics, transcriptomics, proteomics, and metabolomics, have revolutionized the exploration of plant defense mechanisms, unraveling molecular intricacies in response to various stressors. However, the complexity and scale of omics data necessitate sophisticated analytical tools for meaningful insights. This review delves into the application of artificial intelligence algorithms, particularly machine learning and deep learning, as promising approaches for deciphering complex omics data in plant defense research. The overview encompasses key omics techniques and addresses the challenges and limitations inherent in current AI-assisted omics approaches. Moreover, it contemplates potential future directions in this dynamic field. In summary, AI-assisted omics techniques present a robust toolkit, enabling a profound understanding of the molecular foundations of plant defense and paving the way for more effective crop protection strategies amidst climate change and emerging diseases.

Ameliorating the effects of multiple stresses on agronomic traits in crops: modern biotechnological and omics approaches

While global climate change poses a significant environmental threat to agriculture, the increasing population is another big challenge to food security. To address this, developing crop varieties with increased productivity and tolerance to biotic and abiotic stresses is crucial. Breeders must identify traits to ensure higher and consistent yields under inconsistent environmental challenges, possess resilience against emerging biotic and abiotic stresses and satisfy customer demands for safer and more nutritious meals. With the advent of omics-based technologies, molecular tools are now integrated with breeding to understand the molecular genetics of genotype-based traits and develop better climate-smart crops. The rapid development of omics technologies offers an opportunity to generate novel datasets for crop species. Identifying genes and pathways responsible for significant agronomic traits has been made possible by integrating omics data with genetic and phenotypic information. This paper discusses the importance and use of omics-based strategies, including genomics, transcriptomics, proteomics and phenomics, for agricultural and horticultural crop improvement, which aligns with developing better adaptability in these crop species to the changing climate conditions.

Potential application of elastic nets for shared polygenicity detection with adapted threshold selection

Current research suggests that hundreds to thousands of single nucleotide polymorphisms (SNPs) with small to modest effect sizes contribute to the genetic basis of many disorders, a phenomenon labeled as polygenicity. Additionally, many such disorders demonstrate polygenic overlap, in which risk alleles are shared at associated genetic loci. A simple strategy to detect polygenic overlap between two phenotypes is based on rank-ordering the univariate p-values from two genome-wide association studies (GWASs). Although high-dimensional variable selection strategies such as Lasso and elastic nets have been utilized in other GWAS analysis settings, they are yet to be utilized for detecting shared polygenicity. In this paper, we illustrate how elastic nets, with polygenic scores as the dependent variable and with appropriate adaptation in selecting the penalty parameter, may be utilized for detecting a subset of SNPs involved in shared polygenicity. We provide theory to better understand our approaches, and illustrate their utility using synthetic datasets. Results from extensive simulations are presented comparing the elastic net approaches with the rank ordering approach, in various scenarios. Results from simulations studies exhibit one of the elastic net approaches to be superior when the correlations among the SNPs are high. Finally, we apply the methods on two real datasets to illustrate further the capabilities, limitations and differences among the methods.

Plasticity QTLs specifically contribute to the genotype × water availability interaction in maize

KEY MESSAGE

Multi-trial genome wide association study of plasticity indices allow to detect QTLs specifically involved in the genotype x water availability interaction. Concerns regarding high maize yield losses due to increasing occurrences of drought events are growing, and breeders are still looking for molecular markers for drought tolerance. However, the genetic determinism of traits in response to drought is highly complex and identification of causal regions is a tremendous task. Here, we exploit the phenotypic data obtained from four trials carried out on a phenotyping platform, where a diversity panel of 254 maize hybrids was grown under well-watered and water deficit conditions, to investigate the genetic bases of the drought response in maize. To dissociate drought effect from other environmental factors, we performed multi-trial genome-wide association study on well-watered and water deficit phenotypic means, and on phenotypic plasticity indices computed from measurements made for six ecophysiological traits. We identify 102 QTLs and 40 plasticity QTLs. Most of them were new compared to those obtained from a previous study on the same dataset. Our results show that plasticity QTLs cover genetic regions not identified by QTLs. Furthermore, for all ecophysiological traits, except one, plasticity QTLs are specifically involved in the genotype by water availability interaction, for which they explain between 60 and 100% of the variance. Altogether, QTLs and plasticity QTLs captured more than 75% of the genotype by water availability interaction variance, and allowed to find new genetic regions. Overall, our results demonstrate the importance of considering phenotypic plasticity to decipher the genetic architecture of trait response to stress.

Natural plant growth and development achieved in the IPK PhenoSphere by dynamic environment simulation

In plant science, the suboptimal match of growing conditions hampers the transfer of knowledge from controlled environments in glasshouses or climate chambers to field environments. Here we present the PhenoSphere, a plant cultivation infrastructure designed to simulate field-like environments in a reproducible manner. To benchmark the PhenoSphere, the effects on plant growth of weather conditions of a single maize growing season and of an averaged season over three years are compared to those of a standard glasshouse and of four years of field trials. The single season simulation proves superior to the glasshouse and the averaged season in the PhenoSphere: The simulated weather regime of the single season triggers plant growth and development progression very similar to that observed in the field. Hence, the PhenoSphere enables detailed analyses of performance-related trait expression and causal biological mechanisms in plant populations exposed to weather conditions of current and anticipated future climate scenarios.

Genomic prediction and association mapping of maize grain yield in multi-environment trials based on reaction norm models

Genotype-by-environment interaction (GEI) is among the greatest challenges for maize breeding programs. Strong GEI limits both the prediction of genotype performance across variable environmental conditions and the identification of genomic regions associated with grain yield. Incorporating GEI into yield prediction models has been shown to improve prediction accuracy of yield; nevertheless, more work is needed to further understand this complex interaction across populations and environments. The main objectives of this study were to: 1) assess GEI in maize grain yield based on reaction norm models and predict hybrid performance across a gradient of environmental (EG) conditions and 2) perform a genome-wide association study (GWAS) and post-GWAS analyses for maize grain yield using data from 2014 to 2017 of the Genomes to Fields initiative hybrid trial. After quality control, 2,126 hybrids with genotypic and phenotypic data were assessed across 86 environments representing combinations of locations and years, although not all hybrids were evaluated in all environments. Heritability was greater in higher-yielding environments due to an increase in genetic variability in these environments in comparison to the low-yielding environments. GWAS was carried out for yield and five single nucleotide polymorphisms (SNPs) with the highest magnitude of effect were selected in each environment for follow-up analyses. Many candidate genes in proximity of selected SNPs have been previously reported with roles in stress response. Genomic prediction was performed to assess prediction accuracy of previously tested or untested hybrids in environments from a new growing season. Prediction accuracy was 0.34 for cross validation across years (CV0-Predicted EG) and 0.21 for cross validation across years with only untested hybrids (CV00-Predicted EG) when compared to Best Linear Unbiased Prediction (BLUPs) that did not utilize genotypic or environmental relationships. Prediction accuracy improved to 0.80 (CV0-Predicted EG) and 0.60 (CV00-Predicted EG) when compared to the whole-dataset model that used the genomic relationships and the environmental gradient of all environments in the study. These results identify regions of the genome for future selection to improve yield and a methodology to increase the number of hybrids evaluated across locations of a multi-environment trial through genomic prediction.

Genetic diversity and selection signatures in a gene bank panel of maize inbred lines from Southeast Europe compared with two West European panels

Southeast Europe (SEE) is a very important maize-growing region, comparable to the Corn belt region of the United States, with similar dent germplasm (dent by dent hybrids). Historically, this region has undergone several genetic material swaps, following the trends in the US, with one of the most significant swaps related to US aid programs after WWII. The imported accessions used to make double-cross hybrids were also mixed with previously adapted germplasm originating from several more distant OPVs, supporting the transition to single cross-breeding. Many of these materials were deposited at the Maize Gene Bank of the Maize Research Institute Zemun Polje (MRIZP) between the 1960s and 1980s. A part of this Gene Bank (572 inbreds) was genotyped with Affymetrix Axiom Maize Genotyping Array with 616,201 polymorphic variants. Data were merged with two other genotyping datasets with mostly European flint (TUM dataset) and dent (DROPS dataset) germplasm. The final pan-European dataset consisted of 974 inbreds and 460,243 markers. Admixture analysis showed seven ancestral populations representing European flint, B73/B14, Lancaster, B37, Wf9/Oh07, A374, and Iodent pools. Subpanel of inbreds with SEE origin showed a lack of Iodent germplasm, marking its historical context. Several signatures of selection were identified at chromosomes 1, 3, 6, 7, 8, 9, and 10. The regions under selection were mined for protein-coding genes and were used for gene ontology (GO) analysis, showing a highly significant overrepresentation of genes involved in response to stress. Our results suggest the accumulation of favorable allelic diversity, especially in the context of changing climate in the genetic resources of SEE.

VegAnn, Vegetation Annotation of multi-crop RGB images acquired under diverse conditions for segmentation

Applying deep learning to images of cropping systems provides new knowledge and insights in research and commercial applications. Semantic segmentation or pixel-wise classification, of RGB images acquired at the ground level, into vegetation and background is a critical step in the estimation of several canopy traits. Current state of the art methodologies based on convolutional neural networks (CNNs) are trained on datasets acquired under controlled or indoor environments. These models are unable to generalize to real-world images and hence need to be fine-tuned using new labelled datasets. This motivated the creation of the VegAnn - Vegetation Annotation - dataset, a collection of 3775 multi-crop RGB images acquired for different phenological stages using different systems and platforms in diverse illumination conditions. We anticipate that VegAnn will help improving segmentation algorithm performances, facilitate benchmarking and promote large-scale crop vegetation segmentation research.

A root functional-structural model allows assessment of the effects of water deficit on water and solute transport parameters

Root water uptake is driven by a combination of hydrostatic and osmotic forces. Water transport was characterized in primary roots of maize seedlings grown hydroponically under standard and water deficit (WD) conditions, as induced by addition of 150 g l-1 polyethylene glycol 8000 (water potential= -0.336 MPa). Flow measurements were performed using the pressure chamber technique in intact roots or on progressively cut root system architectures. To account for the concomitant transport of water and solutes in roots under WD, we developed within realistic root system architectures a hydraulic tree model integrating both solute pumping and leak. This model explains the high spontaneous sap exudation of roots grown in standard conditions, the non-linearity of pressure-flow relationships, and negative fluxes observed under WD conditions at low external hydrostatic pressure. The model also reveals the heterogeneity of driving forces and elementary radial flows throughout the root system architecture, and how this heterogeneity depends on both plant treatment and water transport mode. The full set of flow measurement data obtained from individual roots grown under standard or WD conditions was used in an inverse modeling approach to determine their respective radial and axial hydraulic conductivities. This approach allows resolution of the dramatic effects of WD on these two components.

Integrated omics approaches for flax improvement under abiotic and biotic stress: Current status and future prospects

Flax (Linum usitatissimum L.) or linseed is one of the important industrial crops grown all over the world for seed oil and fiber. Besides oil and fiber, flax offers a wide range of nutritional and therapeutic applications as a feed and food source owing to high amount of α-linolenic acid (omega-3 fatty acid), lignans, protein, minerals, and vitamins. Periodic losses caused by unpredictable environmental stresses such as drought, heat, salinity-alkalinity, and diseases pose a threat to meet the rising market demand. Furthermore, these abiotic and biotic stressors have a negative impact on biological diversity and quality of oil/fiber. Therefore, understanding the interaction of genetic and environmental factors in stress tolerance mechanism and identification of underlying genes for economically important traits is critical for flax improvement and sustainability. In recent technological era, numerous omics techniques such as genomics, transcriptomics, metabolomics, proteomics, phenomics, and ionomics have evolved. The advancements in sequencing technologies accelerated development of genomic resources which facilitated finer genetic mapping, quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection in major cereal and oilseed crops including flax. Extensive studies in the area of genomics and transcriptomics have been conducted post flax genome sequencing. Interestingly, research has been focused more for abiotic stresses tolerance compared to disease resistance in flax through transcriptomics, while the other areas of omics such as metabolomics, proteomics, ionomics, and phenomics are in the initial stages in flax and several key questions remain unanswered. Little has been explored in the integration of omic-scale data to explain complex genetic, physiological and biochemical basis of stress tolerance in flax. In this review, the current status of various omics approaches for elucidation of molecular pathways underlying abiotic and biotic stress tolerance in flax have been presented and the importance of integrated omics technologies in future research and breeding have been emphasized to ensure sustainable yield in challenging environments.

Physiological adaptive traits are a potential allele reservoir for maize genetic progress under challenging conditions

Combined phenomic and genomic approaches are required to evaluate the margin of progress of breeding strategies. Here, we analyze 65 years of genetic progress in maize yield, which was similar (101 kg ha-1 year-1) across most frequent environmental scenarios in the European growing area. Yield gains were linked to physiologically simple traits (plant phenology and architecture) which indirectly affected reproductive development and light interception in all studied environments, marked by significant genomic signatures of selection. Conversely, studied physiological processes involved in stress adaptation remained phenotypically unchanged (e.g. stomatal conductance and growth sensitivity to drought) and showed no signatures of selection. By selecting for yield, breeders indirectly selected traits with stable effects on yield, but not physiological traits whose effects on yield can be positive or negative depending on environmental conditions. Because yield stability under climate change is desirable, novel breeding strategies may be needed for exploiting alleles governing physiological adaptive traits.

Identification of environment types and adaptation zones with self-organizing maps; applications to sunflower multi-environment data in Europe

We evaluate self-organizing maps (SOM) to identify adaptation zones and visualize multi-environment genotypic responses. We apply SOM to multiple traits and crop growth model output of large-scale European sunflower data. Genotype-by-environment interactions (G × E) complicate the selection of well-adapted varieties. A possible solution is to group trial locations into adaptation zones with G × E occurring mainly between zones. By selecting for good performance inside those zones, response to selection is increased. In this paper, we present a two-step procedure to identify adaptation zones that starts from a self-organizing map (SOM). In the SOM, trials across locations and years are assigned to groups, called units, that are organized on a two-dimensional grid. Units that are further apart contain more distinct trials. In an iterative process of reweighting trial contributions to units, the grid configuration is learnt simultaneously with the trial assignment to units. An aggregation of the units in the SOM by hierarchical clustering then produces environment types, i.e. trials with similar growing conditions. Adaptation zones can subsequently be identified by grouping trial locations with similar distributions of environment types across years. For the construction of SOMs, multiple data types can be combined. We compared environment types and adaptation zones obtained for European sunflower from quantitative traits like yield, oil content, phenology and disease scores with those obtained from environmental indices calculated with the crop growth model Sunflo. We also show how results are affected by input data organization and user-defined weights for genotypes and traits. Adaptation zones for European sunflower as identified by our SOM-based strategy captured substantial genotype-by-location interaction and pointed to trials in Spain, Turkey and South Bulgaria as inducing different genotypic responses.

Genome-Wide Association Study Reveals Genetic Architecture and Candidate Genes for Yield and Related Traits under Terminal Drought, Combined Heat and Drought in Tropical Maize Germplasm

Maize (Zea mays L.) production is constrained by drought and heat stresses. The combination of these two stresses is likely to be more detrimental. To breed for maize cultivars tolerant of these stresses, 162 tropical maize inbred lines were evaluated under combined heat and drought (CHD) and terminal drought (TD) conditions. The mixed linear model was employed for the genome-wide association study using 7834 SNP markers and several phenotypic data including, days to 50% anthesis (AD) and silking (SD), husk cover (HUSKC), and grain yield (GY). In total, 66, 27, and 24 SNPs were associated with the traits evaluated under CHD, TD, and their combined effects, respectively. Of these, four single nucleotide polymorphism (SNP) markers (SNP_161703060 on Chr01, SNP_196800695 on Chr02, SNP_195454836 on Chr05, and SNP_51772182 on Chr07) had pleiotropic effects on both AD and SD under CHD conditions. Four SNPs (SNP_138825271 (Chr03), SNP_244895453 (Chr04), SNP_168561609 (Chr05), and SNP_62970998 (Chr06)) were associated with AD, SD, and HUSKC under TD. Twelve candidate genes containing phytohormone cis-acting regulating elements were implicated in the regulation of plant responses to multiple stress conditions including heat and drought. The SNPs and candidate genes identified in the study will provide invaluable information for breeding climate smart maize varieties under tropical conditions following validation of the SNP markers.

Preparation and Curation of Phenotypic Datasets

Based on case studies, in this chapter we discuss the extent to which the number and identity of quantitative trait loci (QTL) identified from genome-wide association studies (GWAS) are affected by curation and analysis of phenotypic data. The chapter demonstrates through examples the impact of (1) cleaning of outliers, and of (2) the choice of statistical method for estimating genotypic mean values of phenotypic inputs in GWAS. No cleaning of outliers resulted in the highest number of dubious QTL, especially at loci with highly unbalanced allelic frequencies. A trade-off was identified between the risk of false positives and the risk of missing interesting, yet rare alleles. The choice of the statistical method to estimate genotypic mean values also affected the output of GWAS analysis, with reduced QTL overlap between methods. Using mixed models that capture spatial trends, among other features, increased the narrow-sense heritability of traits, the number of identified QTL and the overall power of GWAS analysis. Cleaning and choosing robust statistical models for estimating genotypic mean values should be included in GWAS pipelines to decrease both false positive and false negative rates of QTL detection.

Wavelet Screening: a novel approach to analyzing GWAS data

Background

Traditional methods for single-variant genome-wide association study (GWAS) incur a substantial multiple-testing burden because of the need to test for associations with a vast number of single-nucleotide polymorphisms (SNPs) simultaneously. Further, by ignoring more complex joint effects of nearby SNPs within a given region, these methods fail to consider the genomic context of an association with the outcome.

Results

To address these shortcomings, we present a more powerful method for GWAS, coined ‘Wavelet Screening’ (WS), that greatly reduces the number of tests to be performed. This is achieved through the use of a sliding-window approach based on wavelets to sequentially screen the entire genome for associations. Wavelets are oscillatory functions that are useful for analyzing the local frequency and time behavior of signals. The signals can then be divided into different scale components and analyzed separately. In the current setting, we consider a sequence of SNPs as a genetic signal, and for each screened region, we transform the genetic signal into the wavelet space. The null and alternative hypotheses are modeled using the posterior distribution of the wavelet coefficients. WS is enhanced by using additional information from the regression coefficients and by taking advantage of the pyramidal structure of wavelets. When faced with more complex genetic signals than single-SNP associations, we show via simulations that WS provides a substantial gain in power compared to both the traditional GWAS modeling and another popular regional association test called SNP-set (Sequence) Kernel Association Test (SKAT). To demonstrate feasibility, we applied WS to a large Norwegian cohort (N=8006) with genotypes and information available on gestational duration.

Conclusions

WS is a powerful and versatile approach to analyzing whole-genome data and lends itself easily to investigating various omics data types. Given its broader focus on the genomic context of an association, WS may provide additional insight into trait etiology by revealing genes and loci that might have been missed by previous efforts.

Impacts of environmental conditions, and allelic variation of cytosolic glutamine synthetase on maize hybrid kernel production

Cytosolic glutamine synthetase (GS1) is the enzyme mainly responsible of ammonium assimilation and reassimilation in maize leaves. The agronomic potential of GS1 in maize kernel production was investigated by examining the impact of an overexpression of the enzyme in the leaf cells. Transgenic hybrids exhibiting a three-fold increase in leaf GS activity were produced and characterized using plants grown in the field. Several independent hybrids overexpressing Gln1-3, a gene encoding cytosolic (GS1), in the leaf and bundle sheath mesophyll cells were grown over five years in different locations. On average, a 3.8% increase in kernel yield was obtained in the transgenic hybrids compared to controls. However, we observed that such an increase was simultaneously dependent upon both the environmental conditions and the transgenic event for a given field trial. Although variable from one environment to another, significant associations were also found between two GS1 genes (Gln1-3 and Gln1-4) polymorphic regions and kernel yield in different locations. We propose that the GS1 enzyme is a potential lead for producing high yielding maize hybrids using either genetic engineering or marker-assisted selection. However, for these hybrids, yield increases will be largely dependent upon the environmental conditions used to grow the plants.

Amiour et al. use a multi-year field trial evaluation and association mapping to determine if increased enzyme activity and native allelic variations at the GS1 loci in maize contribute to differences in grain yield. Overexpression of GS1 and polymorphisms in the corresponding loci were associated with kernel yield, indicating that GS1 expression can directly control kernel production and that GS1 has a potential lead in the production of high yielding maize hybrids depending on environmental conditions.

Optimizing Genomic-Enabled Prediction in Small-Scale Maize Hybrid Breeding Programs: A Roadmap Review

The usefulness of genomic prediction (GP) for many animal and plant breeding programs has been highlighted for many studies in the last 20 years. In maize breeding programs, mostly dedicated to delivering more highly adapted and productive hybrids, this approach has been proved successful for both large- and small-scale breeding programs worldwide. Here, we present some of the strategies developed to improve the accuracy of GP in tropical maize, focusing on its use under low budget and small-scale conditions achieved for most of the hybrid breeding programs in developing countries. We highlight the most important outcomes obtained by the University of São Paulo (USP, Brazil) and how they can improve the accuracy of prediction in tropical maize hybrids. Our roadmap starts with the efforts for germplasm characterization, moving on to the practices for mating design, and the selection of the genotypes that are used to compose the training population in field phenotyping trials. Factors including population structure and the importance of non-additive effects (dominance and epistasis) controlling the desired trait are also outlined. Finally, we explain how the source of the molecular markers, environmental, and the modeling of genotype-environment interaction can affect the accuracy of GP. Results of 7 years of research in a public maize hybrid breeding program under tropical conditions are discussed, and with the great advances that have been made, we find that what is yet to come is exciting. The use of open-source software for the quality control of molecular markers, implementing GP, and envirotyping pipelines may reduce costs in an efficient computational manner. We conclude that exploring new models/tools using high-throughput phenotyping data along with large-scale envirotyping may bring more resolution and realism when predicting genotype performances. Despite the initial costs, mostly for genotyping, the GP platforms in combination with these other data sources can be a cost-effective approach for predicting the performance of maize hybrids for a large set of growing conditions.

Applications of Multi-Omics Technologies for Crop Improvement

Multiple "omics" approaches have emerged as successful technologies for plant systems over the last few decades. Advances in next-generation sequencing (NGS) have paved a way for a new generation of different omics, such as genomics, transcriptomics, and proteomics. However, metabolomics, ionomics, and phenomics have also been well-documented in crop science. Multi-omics approaches with high throughput techniques have played an important role in elucidating growth, senescence, yield, and the responses to biotic and abiotic stress in numerous crops. These omics approaches have been implemented in some important crops including wheat (Triticum aestivum L.), soybean (Glycine max), tomato (Solanum lycopersicum), barley (Hordeum vulgare L.), maize (Zea mays L.), millet (Setaria italica L.), cotton (Gossypium hirsutum L.), Medicago truncatula, and rice (Oryza sativa L.). The integration of functional genomics with other omics highlights the relationships between crop genomes and phenotypes under specific physiological and environmental conditions. The purpose of this review is to dissect the role and integration of multi-omics technologies for crop breeding science. We highlight the applications of various omics approaches, such as genomics, transcriptomics, proteomics, metabolomics, phenomics, and ionomics, and the implementation of robust methods to improve crop genetics and breeding science. Potential challenges that confront the integration of multi-omics with regard to the functional analysis of genes and their networks as well as the development of potential traits for crop improvement are discussed. The panomics platform allows for the integration of complex omics to construct models that can be used to predict complex traits. Systems biology integration with multi-omics datasets can enhance our understanding of molecular regulator networks for crop improvement. In this context, we suggest the integration of entire omics by employing the "phenotype to genotype" and "genotype to phenotype" concept. Hence, top-down (phenotype to genotype) and bottom-up (genotype to phenotype) model through integration of multi-omics with systems biology may be beneficial for crop breeding improvement under conditions of environmental stresses.

Genetic Dissection of Grain Yield of Maize and Yield-Related Traits Through Association Mapping and Genomic Prediction

High yield is the primary objective of maize breeding. Genomic dissection of grain yield and yield-related traits contribute to understanding the yield formation and improving the yield of maize. In this study, two genome-wide association study (GWAS) methods and genomic prediction were made on an association panel of 309 inbred lines. GWAS analyses revealed 22 significant trait-marker associations for grain yield per plant (GYP) and yield-related traits. Genomic prediction analyses showed that reproducing kernel Hilbert space (RKHS) outperformed the other four models based on GWAS-derived markers for GYP, ear weight, kernel number per ear and row, ear length, and ear diameter, whereas genomic best linear unbiased prediction (GBLUP) showed a slight superiority over other modes in most subsets of the trait-associated marker (TAM) for thousand kernel weight and kernel row number. The prediction accuracy could be improved when significant single-nucleotide polymorphisms were fitted as the fixed effects. Integrating information on population structure into the fixed model did not improve the prediction performance. For GYP, the prediction accuracy of TAMs derived from fixed and random model Circulating Probability Unification (FarmCPU) was comparable to that of the compressed mixed linear model (CMLM). For yield-related traits, CMLM-derived markers provided better accuracies than FarmCPU-derived markers in most scenarios. Compared with all markers, TAMs could effectively improve the prediction accuracies for GYP and yield-related traits. For eight traits, moderate- and high-prediction accuracies were achieved using TAMs. Taken together, genomic prediction incorporating prior information detected by GWAS could be a promising strategy to improve the grain yield of maize.

Discovery of beneficial haplotypes for complex traits in maize landraces

Genetic variation is of crucial importance for crop improvement. Landraces are valuable sources of diversity, but for quantitative traits efficient strategies for their targeted utilization are lacking. Here, we map haplotype-trait associations at high resolution in ~1000 doubled-haploid lines derived from three maize landraces to make their native diversity for early development traits accessible for elite germplasm improvement. A comparative genomic analysis of the discovered haplotypes in the landrace-derived lines and a panel of 65 breeding lines, both genotyped with 600k SNPs, points to untapped beneficial variation for target traits in the landraces. The superior phenotypic performance of lines carrying favorable landrace haplotypes as compared to breeding lines with alternative haplotypes confirms these findings. Stability of haplotype effects across populations and environments as well as their limited effects on undesired traits indicate that our strategy has high potential for harnessing beneficial haplotype variation for quantitative traits from genetic resources.

Genetic variations present in landraces are critical for crop genetic improvement. Here, the authors map haplotype-trait associations in ~1000 doubled haploid lines derived from three European maize landraces and identify beneficial haplotypes for quantitative traits that are not present in breeding lines.

Simulating the effect of flowering time on maize individual leaf area in contrasting environmental scenarios

The quality of yield prediction is linked to that of leaf area. We first analysed the consequences of flowering time and environmental conditions on the area of individual leaves in 127 genotypes presenting contrasting flowering times in fields of Europe, Mexico, and Kenya. Flowering time was the strongest determinant of leaf area. Combined with a detailed field experiment, this experiment showed a large effect of flowering time on the final leaf number and on the distribution of leaf growth rate and growth duration along leaf ranks, in terms of both length and width. Equations with a limited number of genetic parameters predicted the beginning, end, and maximum growth rate (length and width) for each leaf rank. The genotype-specific environmental effects were analysed with datasets in phenotyping platforms that assessed the effects (i) of the amount of intercepted light on leaf width, and (ii) of temperature, evaporative demand, and soil water potential on leaf elongation rate. The resulting model was successfully tested for 31 hybrids in 15 European and Mexican fields. It potentially allows prediction of the vertical distribution of leaf area of a large number of genotypes in contrasting field conditions, based on phenomics and on sensor networks.

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.

Root architecture and hydraulics converge for acclimation to changing water availability

Because of intense transpiration and growth, the needs of plants for water can be immense. Yet water in the soil is most often heterogeneous if not scarce due to more and more frequent and intense drought episodes. The converse context, flooding, is often associated with marked oxygen deficiency and can also challenge the plant water status. Under our feet, roots achieve an incredible challenge to meet the water demand of the plant's aerial parts under such dramatically different environmental conditions. For this, they continuously explore the soil, building a highly complex, branched architecture. On shorter time scales, roots keep adjusting their water transport capacity (their so-called hydraulics) locally or globally. While the mechanisms that directly underlie root growth and development as well as tissue hydraulics are being uncovered, the signalling mechanisms that govern their local and systemic adjustments as a function of water availability remain largely unknown. A comprehensive understanding of root architecture and hydraulics as a whole (in other terms, root hydraulic architecture) is needed to apprehend the strategies used by plants to optimize water uptake and possibly improve crops regarding this crucial trait.

Reconstruction of Networks with Direct and Indirect Genetic Effects

Genetic variance of a phenotypic trait can originate from direct genetic effects, or from indirect effects, i.e., through genetic effects on other traits, affecting the trait of interest. This distinction is often of great importance, for example, when trying to improve crop yield and simultaneously control plant height. As suggested by Sewall Wright, assessing contributions of direct and indirect effects requires knowledge of (1) the presence or absence of direct genetic effects on each trait, and (2) the functional relationships between the traits. Because experimental validation of such relationships is often unfeasible, it is increasingly common to reconstruct them using causal inference methods. However, most current methods require all genetic variance to be explained by a small number of quantitative trait loci (QTL) with fixed effects. Only a few authors have considered the "missing heritability" case, where contributions of many undetectable QTL are modeled with random effects. Usually, these are treated as nuisance terms that need to be eliminated by taking residuals from a multi-trait mixed model (MTM). But fitting such an MTM is challenging, and it is impossible to infer the presence of direct genetic effects. Here, we propose an alternative strategy, where genetic effects are formally included in the graph. This has important advantages: (1) genetic effects can be directly incorporated in causal inference, implemented via our PCgen algorithm, which can analyze many more traits; and (2) we can test the existence of direct genetic effects, and improve the orientation of edges between traits. Finally, we show that reconstruction is much more accurate if individual plant or plot data are used, instead of genotypic means. We have implemented the PCgen-algorithm in the R-package pcgen.

The Genetic Architecture for Phenotypic Plasticity of the Rice Grain Ionome

The ionome of the rice grain is crucial for the health of populations that consume rice as a staple food. However, the contribution of phenotypic plasticity to the variation of rice grain ionome and the genetic architecture of phenotypic plasticity are poorly understood. In this study, we investigated the rice grain ionome of a rice diversity panel in up to eight environments. A considerable proportion of phenotypic variance can be attributed to phenotypic plasticity. Then, phenotypic plasticity and mean phenotype were quantified using Bayesian Finlay-Wilkinson regression, and a significant correlation between them was observed. However, the genetic architecture of mean phenotype was distinct from that of phenotypic plasticity. Also, the correlation between them was mainly attributed to the phenotypic divergence between rice subspecies. Furthermore, the results of whole-genome regression analysis showed that the genetic loci related to phenotypic plasticity can explain a considerable proportion of the phenotypic variance in some environments, especially for Cd, Cu, Mn, and Zn. Our study not only sheds light on the genetic architecture of phenotypic plasticity of the rice grain ionome but also suggests that the genetic loci which related to phenotypic plasticity are valuable in rice grain ionome improvement breeding.

Seed Weight as a Covariate in Association and Prediction Studies for Biomass Traits in Maize Seedlings

BACKGROUND

The seedling stage has received little attention in maize breeding to identify genotypes tolerant to water deficit. The aim of this study is to evaluate incorporation of seed weight (expressed as hundred kernel weight, HKW) as a covariate into genomic association and prediction studies for three biomass traits in a panel of elite inbred lines challenged by water withholding at seedling stage.

METHODS

109 genotyped-by-sequencing (GBS) elite maize inbreds were phenotyped for HKW and planted in controlled conditions (16/8 day/night, 25 °C, 50% RH, 200 µMol/m2/s) in trays filled with soil. Plants in control (C) were watered every two days, while watering was stopped for 10 days in water withholding (WW). Fresh weight (FW), dry weight (DW), and dry matter content (DMC) were measured.

RESULTS

Adding HKW as a covariate increased the power of detection of associations in FW and DW by 44% and increased genomic prediction accuracy in C and decreased in WW.

CONCLUSIONS

Seed weight was effectively incorporated into association studies for biomass traits in maize seedlings, whereas the incorporation into genomic predictions, particularly in water-stressed plants, was not worthwhile.

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.

Can Metabolite- and Transcript-Based Selection for Drought Tolerance in Solanum tuberosum Replace Selection on Yield in Arid Environments?

Climate models predict an increased likelihood of drought, demanding efficient selection for drought tolerance to maintain yield stability. Classic tolerance breeding relies on selection for yield in arid environments, which depends on yield trials and takes decades. Breeding could be accelerated by marker-assisted selection (MAS). As an alternative to genomic markers, transcript and metabolite markers have been suggested for important crops but also for orphan corps. For potato, we suggested a random-forest-based model that predicts tolerance from leaf metabolite and transcript levels with a precision of more than 90% independent of the agro-environment. To find out how the model based selection compares to yield-based selection in arid environments, we applied this approach to a population of 200 tetraploid Solanum tuberosum ssp. tuberosum lines segregating for drought tolerance. Twenty-four lines were selected into a phenotypic subpopulation (PP_t) for superior tolerance based on relative tuber starch yield data from three drought stress trials. Two subpopulations with superior (MP_t) and inferior (MP_s) tolerance were selected based on drought tolerance predictions based on leaf metabolite and transcript levels from two sites. The 60 selected lines were phenotyped for yield and drought tolerance in 10 multi-environment drought stress trials representing typical Central European drought scenarios. Neither selection affected development or yield potential. Lines with superior drought tolerance and high yields under stress were over-represented in both populations selected for superior tolerance, with a higher number in PP_t compared to MP_t. However, selection based on leaf metabolites may still be an alternative to yield-based selection in arid environments as it works on leaves sampled in breeder's fields independent of drought trials. As the selection against low tolerance was ineffective, the method is best used in combination with tools that select against sensitive genotypes. Thus, metabolic and transcript marker-based selection for drought tolerance is a viable alternative to the selection on yield in arid environments.

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

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.

Using environmental clustering to identify specific drought tolerance QTLs in bread wheat (T. aestivum L.)

Environmental clustering helps to identify QTLs associated with grain yield in different water stress scenarios. These QTLs could be useful for breeders to improve grain yields and increase genetic resilience in marginal environments. Drought is one of the main abiotic stresses limiting winter bread wheat growth and productivity around the world. The acquisition of new high-yielding and stress-tolerant varieties is therefore necessary and requires improved understanding of the physiological and genetic bases of drought resistance. A panel of 210 elite European varieties was evaluated in 35 field trials. Grain yield and its components were scored in each trial. A crop model was then run with detailed climatic data and soil water status to assess the dynamics of water stress in each environment. Varieties were registered from 1992 to 2011, allowing us to test timewise genetic progress. Finally, a genome-wide association study (GWAS) was carried out using genotyping data from a 280 K SNP chip. The crop model simulation allowed us to group the environments into four water stress scenarios: an optimal condition with no water stress, a post-anthesis water stress, a moderate-anthesis water stress and a high pre-anthesis water stress. Compared to the optimal water condition, grain yield losses in the stressed conditions were 3.3%, 12.4% and 31.2%, respectively. This environmental clustering improved understanding of the effect of drought on grain yields and explained 20% of the G × E interaction. The greatest genetic progress was obtained in the optimal condition, mostly represented in France. The GWAS identified several QTLs, some of which were specific of the different water stress patterns. Our results make breeding for improved drought resistance to specific environmental scenarios easier and will facilitate genetic progress in future environments, i.e., water stress environments.

Exome sequences and multi-environment field trials elucidate the genetic basis of adaptation in barley

Broadening the genetic base of crops is crucial for developing varieties to respond to global agricultural challenges such as climate change. Here, we analysed a diverse panel of 371 domesticated lines of the model crop barley to explore the genetics of crop adaptation. We first collected exome sequence data and phenotypes of key life history traits from contrasting multi-environment common garden trials. Then we applied refined statistical methods, including some based on exomic haplotype states, for genotype-by-environment (G×E) modelling. Sub-populations defined from exomic profiles were coincident with barley's biology, geography and history, and explained a high proportion of trial phenotypic variance. Clear G×E interactions indicated adaptation profiles that varied for landraces and cultivars. Exploration of circadian clock-related genes, associated with the environmentally adaptive days to heading trait (crucial for the crop's spread from the Fertile Crescent), illustrated complexities in G×E effect directions, and the importance of latitudinally based genic context in the expression of large-effect alleles. Our analysis supports a gene-level scientific understanding of crop adaption and leads to practical opportunities for crop improvement, allowing the prioritisation of genomic regions and particular sets of lines for breeding efforts seeking to cope with climate change and other stresses.

Genotyping-by-sequencing and SNP-arrays are complementary for detecting quantitative trait loci by tagging different haplotypes in association studies

BACKGROUND

Single Nucleotide Polymorphism (SNP) array and re-sequencing technologies have different properties (e.g. calling rate, minor allele frequency profile) and drawbacks (e.g. ascertainment bias). This lead us to study their complementarity and the consequences of using them separately or combined in diversity analyses and Genome-Wide Association Studies (GWAS). We performed GWAS on three traits (grain yield, plant height and male flowering time) measured in 22 environments on a panel of 247 F1 hybrids obtained by crossing 247 diverse dent maize inbred lines with a same flint line. The 247 lines were genotyped using three genotyping technologies (Genotyping-By-Sequencing, Illumina Infinium 50 K and Affymetrix Axiom 600 K arrays).

RESULTS

The effects of ascertainment bias of the 50 K and 600 K arrays were negligible for deciphering global genetic trends of diversity and for estimating relatedness in this panel. We developed an original approach based on linkage disequilibrium (LD) extent in order to determine whether SNPs significantly associated with a trait and that are physically linked should be considered as a single Quantitative Trait Locus (QTL) or several independent QTLs. Using this approach, we showed that the combination of the three technologies, which have different SNP distributions and densities, allowed us to detect more QTLs (gain in power) and potentially refine the localization of the causal polymorphisms (gain in resolution).

CONCLUSIONS

Conceptually different technologies are complementary for detecting QTLs by tagging different haplotypes in association studies. Considering LD, marker density and the combination of different technologies (SNP-arrays and re-sequencing), the genotypic data available were most likely enough to well represent polymorphisms in the centromeric regions, whereas using more markers would be beneficial for telomeric regions.

Genetic Correlations Between Photosynthetic and Yield Performance in Maize Are Different Under Two Heat Scenarios During Flowering

Chlorophyll fluorescence (ChlF) parameters are reliable early stress indicators in crops, but their relations with yield are still not clear. The aims of this study are to examine genetic correlations between photosynthetic performance of JIP-test during flowering and grain yield (GY) in maize grown under two heat scenarios in the field environments applying quantitative genetic analysis, and to compare efficiencies of indirect selection for GY through ChlF parameters and genomic selection for GY. The testcrosses of 221 intermated recombinant inbred lines (IRILs) of the IBM Syn4 population were evaluated in six environments at two geographically distinctive locations in 3 years. According to day/night temperatures and vapor pressure deficit (VPD), the two locations in Croatia and Turkey may be categorized to the mild heat and moderate heat scenarios, respectively. Mild heat scenario is characterized by daytime temperatures often exceeding 33°C and night temperatures lower than 20°C while in moderate heat scenario the daytime temperatures often exceeded 33°C and night temperatures were above 20°C. The most discernible differences among the scenarios were obtained for efficiency of electron transport beyond quinone A (Q_A) [ET/(TR-ET)], performance index on absorption basis (PI_ABS) and GY. Under the moderate heat scenario, there were tight positive genetic correlations between ET/(TR-ET) and GY (0.73), as well as between PI_ABS and GY (0.59). Associations between the traits were noticeably weaker under the mild heat scenario. Analysis of quantitative trait loci (QTL) revealed several common QTLs for photosynthetic and yield performance under the moderate heat scenario corroborating pleiotropy. Although the indirect selection with ChlF parameters is less efficient than direct selection, ET/(TR-ET) and PI_ABS could be efficient secondary breeding traits for selection under moderate heat stress since they seem to be genetically correlated with GY in the stressed environments and not associated with yield performance under non-stressed conditions predicting GY during flowering. Indirect selection through PI_ABS was also shown to be more efficient than genomic selection in moderate heat scenario.

Genetic and environmental dissection of biomass accumulation in multi-genotype maize canopies

Multi-genotype canopies are frequent in phenotyping experiments and are of increasing interest in agriculture. Radiation interception efficiency (RIE) and radiation use efficiency (RUE) have low heritabilities in such canopies. We propose a revised Monteith equation that identifies environmental and genetic components of RIE and RUE. An environmental term, a component of RIE, characterizes the effect of the presence or absence of neighbours on light interception. The ability of a given plant to compete with its neighbours is then identified, which accounts for the genetic variability of RIE of plants having similar leaf areas. This method was used in three experiments in a phenotyping platform with 765 plants of 255 maize hybrids. As expected, the heritability of the environmental term was near zero, whereas that of the competitiveness term increased with phenological stage, resulting in the identification of quantitative trait loci. In the same way, RUE was dissected as an effect of intercepted light and a genetic term. This approach was used for predicting the behaviour of individual genotypes in virtual multi-genotype canopies. A large effect of competitiveness was observed in multi-genotype but not in single-genotype canopies, resulting in a bias for genotype comparisons in breeding fields.

Water Deficit-Responsive QTLs for Cell Wall Degradability and Composition in Maize at Silage Stage

The use of lignocellulosic biomass for animal feed or biorefinery requires the optimization of its degradability. Moreover, biomass crops need to be better adapted to the changing climate and in particular to periods of drought. Although the negative impact of water deficit on biomass yield has often been mentioned, its impact on biomass quality has only been recently reported in a few species. In the present study, we combined the mapping power of a maize recombinant inbred line population with robust near infrared spectroscopy predictive equations to track the response to water deficit of traits associated with biomass quality. The population was cultivated under two contrasted water regimes over 3 consecutive years in the south of France and harvested at silage stage. We showed that cell wall degradability and β-O-4-linked H lignin subunits were increased in response to water deficit, while lignin and p-coumaric acid contents were reduced. A mixed linear model was fitted to map quantitative trait loci (QTLs) for agronomical and cell wall-related traits. These QTLs were categorized as "constitutive" (QTL with an effect whatever the irrigation condition) or "responsive" (QTL involved in the response to water deficit) QTLs. Fifteen clusters of QTLs encompassed more than two third of the 213 constitutive QTLs and 13 clusters encompassed more than 60% of the 149 responsive QTLs. Interestingly, we showed that only half of the responsive QTLs co-localized with constitutive and yield QTLs, suggesting that specific genetic factors support biomass quality response to water deficit. Overall, our results demonstrate that water deficit favors cell wall degradability and that breeding of varieties that reconcile improved drought-tolerance and biomass degradability is possible.

Genome-wide association mapping of grain yield in a diverse collection of spring wheat (Triticum aestivum L.) evaluated in southern Australia

Wheat landraces, wild relatives and other 'exotic' accessions are important sources of new favorable alleles. The use of those exotic alleles is facilitated by having access to information on the association of specific genomic regions with desirable traits. Here, we conducted a genome-wide association study (GWAS) using a wheat panel that includes landraces, synthetic hexaploids and other exotic wheat accessions to identify loci that contribute to increases in grain yield in southern Australia. The 568 accessions were grown in the field during the 2014 and 2015 seasons and measured for plant height, maturity, spike length, spike number, grain yield, plant biomass, HI and TGW. We used the 90K SNP array and two GWAS approaches (GAPIT and QTCAT) to identify loci associated with the different traits. We identified 17 loci with GAPIT and 25 with QTCAT. Ten of these loci were associated with known genes that are routinely employed in marker assisted selection such as Ppd-D1 for maturity and Rht-D1 for plant height and seven of those were detected with both methods. We identified one locus for yield per se in 2014 on chromosome 6B with QTCAT and three in 2015, on chromosomes 4B and 5A with GAPIT and 6B with QTCAT. The 6B loci corresponded to the same region in both years. The favorable haplotypes for yield at the 5A and 6B loci are widespread in Australian accessions with 112 out of 153 carrying the favorable haplotype at the 5A locus and 136 out of 146 carrying the favorable haplotype at the 6A locus, while the favorable haplotype at 4B is only present in 65 out of 149 Australian accessions. The low number of yield QTL in our study corroborate with other GWAS for yield in wheat, where most of the identified loci have very small effects.

Heliaphen, an Outdoor High-Throughput Phenotyping Platform for Genetic Studies and Crop Modeling

Heliaphen is an outdoor platform designed for high-throughput phenotyping. It allows the automated management of drought scenarios and monitoring of plants throughout their lifecycles. A robot moving between plants growing in 15-L pots monitors the plant water status and phenotypes the leaf or whole-plant morphology. From these measurements, we can compute more complex traits, such as leaf expansion (LE) or transpiration rate (TR) in response to water deficit. Here, we illustrate the capabilities of the platform with two practical cases in sunflower (Helianthus annuus): a genetic and genomic study of the response of yield-related traits to drought, and a modeling study using measured parameters as inputs for a crop simulation. For the genetic study, classical measurements of thousand-kernel weight (TKW) were performed on a biparental population under automatically managed drought stress and control conditions. These data were used for an association study, which identified five genetic markers of the TKW drought response. A complementary transcriptomic analysis identified candidate genes associated with these markers that were differentially expressed in the parental backgrounds in drought conditions. For the simulation study, we used a crop simulation model to predict the impact on crop yield of two traits measured on the platform (LE and TR) for a large number of environments. We conducted simulations in 42 contrasting locations across Europe using 21 years of climate data. We defined the pattern of abiotic stresses occurring at the continental scale and identified ideotypes (i.e., genotypes with specific trait values) that are more adapted to specific environment types. This study exemplifies how phenotyping platforms can assist the identification of the genetic architecture controlling complex response traits and facilitate the estimation of ecophysiological model parameters to define ideotypes adapted to different environmental conditions.

Computational aspects underlying genome to phenome analysis in plants

Recent advances in genomics technologies have greatly accelerated the progress in both fundamental plant science and applied breeding research. Concurrently, high-throughput plant phenotyping is becoming widely adopted in the plant community, promising to alleviate the phenotypic bottleneck. While these technological breakthroughs are significantly accelerating quantitative trait locus (QTL) and causal gene identification, challenges to enable even more sophisticated analyses remain. In particular, care needs to be taken to standardize, describe and conduct experiments robustly while relying on plant physiology expertise. In this article, we review the state of the art regarding genome assembly and the future potential of pangenomics in plant research. We also describe the necessity of standardizing and describing phenotypic studies using the Minimum Information About a Plant Phenotyping Experiment (MIAPPE) standard to enable the reuse and integration of phenotypic data. In addition, we show how deep phenotypic data might yield novel trait-trait correlations and review how to link phenotypic data to genomic data. Finally, we provide perspectives on the golden future of machine learning and their potential in linking phenotypes to genomic features.

Genomic selection efficiency and a priori estimation of accuracy in a structured dent maize panel

Population structure affects genomic selection efficiency as well as the ability to forecast accuracy using standard GBLUP. Genomic prediction models usually assume that the individuals used for calibration belong to the same population as those to be predicted. Most of the a priori indicators of precision, such as the coefficient of determination (CD), were derived from those same models. But genetic structure is a common feature in plant species, and it may impact genomic selection efficiency and the ability to forecast prediction accuracy. We investigated the impact of genetic structure in a dent maize panel ("Amaizing Dent") using different scenarios including within- or across-group predictions. For a given training set size, the best accuracies were achieved when predicting individuals using a model calibrated on the same genetic group. Nevertheless, a diverse training set representing all the groups had a certain predictive efficiency for all the validation sets, and adding extra-group individuals was almost always beneficial. It underlines the potential of such a generic training set for dent maize genomic selection applications. Alternative prediction models, taking genetic structure explicitly into account, did not improve the prediction accuracy compared to GBLUP. We also investigated the ability of different indicators of precision to forecast accuracy in the within- or across-group scenarios. There was a global encouraging trend of the CD to differentiate scenarios, although there were specific combinations of target populations and traits where the efficiency of this indicator proved to be null. One hypothesis to explain such erratic performances is the impact of genetic structure through group-specific allele diversity at QTLs rather than group-specific allele effects.

Harnessing Phenotypic Plasticity to Improve Maize Yields

Plants can produce different phenotypes when exposed to different environments. Understanding the genetic basis of these plastic responses is crucial for crop breeding efforts. We discuss two recent studies that suggest that yield plasticity in maize has been under selection but is controlled by different genes than yield.

Using Biotechnology-Led Approaches to Uplift Cereal and Food Legume Yields in Dryland Environments

Drought and heat in dryland agriculture challenge the enhancement of crop productivity and threaten global food security. This review is centered on harnessing genetic variation through biotechnology-led approaches to select for increased productivity and stress tolerance that will enhance crop adaptation in dryland environments. Peer-reviewed literature, mostly from the last decade and involving experiments with at least two seasons' data, form the basis of this review. It begins by highlighting the adverse impact of the increasing intensity and duration of drought and heat stress due to global warming on crop productivity and its impact on food and nutritional security in dryland environments. This is followed by (1) an overview of the physiological and molecular basis of plant adaptation to elevated CO2 (eCO2), drought, and heat stress; (2) the critical role of high-throughput phenotyping platforms to study phenomes and genomes to increase breeding efficiency; (3) opportunities to enhance stress tolerance and productivity in food crops (cereals and grain legumes) by deploying biotechnology-led approaches [pyramiding quantitative trait loci (QTL), genomic selection, marker-assisted recurrent selection, epigenetic variation, genome editing, and transgene) and inducing flowering independent of environmental clues to match the length of growing season; (4) opportunities to increase productivity in C3 crops by harnessing novel variations (genes and network) in crops' (C3, C4) germplasm pools associated with increased photosynthesis; and (5) the adoption, impact, risk assessment, and enabling policy environments to scale up the adoption of seed-technology to enhance food and nutritional security. This synthesis of technological innovations and insights in seed-based technology offers crop genetic enhancers further opportunities to increase crop productivity in dryland environments.

Genome-wide association reveals novel genomic loci controlling rice grain yield and its component traits under water-deficit stress during the reproductive stage

A diversity panel comprising of 296 indica rice genotypes was phenotyped under non-stress and water-deficit stress conditions during the reproductive stage in the 2013 and 2014 dry seasons (DSs) at IRRI, Philippines. We investigated the genotypic variability for grain yield, yield components, and related traits, and conducted genome-wide association studies (GWAS) using high-density 45K single nucleotide polymorphisms. We detected 38 loci in 2013 and 64 loci in 2014 for non-stress conditions and 69 loci in 2013 and 55 loci in 2014 for water-deficit stress. Desynchronized flowering time confounded grain yield and its components under water-deficit stress in the 2013 experiment. Statistically corrected grain yield and yield component values using days to flowering helped to detect 31 additional genetic loci for grain yield, its components, and the harvest index in 2013. There were few overlaps in the detected loci between years and treatments, and when compared with previous studies using the same panel, indicating the complexity of yield formation under stress. Nevertheless, our analyses provided important insights into the potential links between grain yield with seed set and assimilate partitioning. Our findings demonstrate the complex genetic architecture of yield formation and we propose exploring the genetic basis of less complex component traits as an alternative route for further yield enhancement.

Landraces of snake melon, an ancient Middle Eastern crop, reveal extensive morphological and DNA diversity for potential genetic improvement

Background

Snake melon (Cucumis melo var. flexuosus, “Faqqous”) is a traditional and ancient vegetable in the Mediterranean area. A collection of landraces from 42 grower fields in Israel and Palestinian territories was grown and characterized in a “Common Garden” rain-fed experiment, at the morphological-horticultural and molecular level using seq-DArT markers.

Results

The different landraces (“populations”) showed extensive variation in morphology and quantitative traits such as yield and femaleness, and clustered into four horticultural varieties. Yield was assessed by five harvests along the season, with middle harvests producing the highest yields. Yield correlated with early vigor, and with femaleness, but not with late vigor. At the molecular level, 2784 SNP were produced and > 90% were mapped to the melon genome. Populations were very polymorphic (46–72% of the markers biallelic in a 4 individuals sample), and observed heterozygosity was higher than the expected, suggesting gene flow among populations and extensive cross pollination among individuals in the field. Genetic distances between landraces were significantly correlated with the geographical distance between collecting sites, and with long term March precipitation average; variation in yield correlated with April temperature maxima.

Conclusions

The extensive variation suggests that selection of local snake melon could result in yield improvement. Correlations between traits and climatic variables could suggest local adaptation of landraces to the diverse environment in which they evolved. This study stresses the importance of preserving this germplasm, and its potential for breeding better snake melons as an heirloom crop in our region.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0619-6) contains supplementary material, which is available to authorized users.