Genetic control of abiotic stress-related specialized metabolites in sunflower.

BACKGROUND: Abiotic stresses in plants include all the environmental conditions that significantly reduce yields, like drought and heat. One of the most significant effects they exert at the cellular level is the accumulation of reactive oxygen species, which cause extensive damage. Plants possess two mechanisms to counter these molecules, i.e. detoxifying enzymes and non-enzymatic antioxidants, which include many classes of specialized metabolites. Sunflower, the fourth global oilseed, is considered moderately drought resistant. Abiotic stress tolerance in this crop has been studied using many approaches, but the control of specialized metabolites in this context remains poorly understood. Here, we performed the first genome-wide association study using abiotic stress-related specialized metabolites as molecular phenotypes in sunflower. After analyzing leaf specialized metabolites of 450 hybrids using liquid chromatography-mass spectrometry, we selected a subset of these compounds based on their association with previously known abiotic stress-related quantitative trait loci. Eventually, we characterized these molecules and their associated genes. RESULTS: We putatively annotated 30 compounds which co-localized with abiotic stress-related quantitative trait loci and which were associated to seven most likely candidate genes. A large proportion of these compounds were potential antioxidants, which was in agreement with the role of specialized metabolites in abiotic stresses. The seven associated most likely candidate genes, instead, mainly belonged to cytochromes P450 and glycosyltransferases, two large superfamilies which catalyze greatly diverse reactions and create a wide variety of chemical modifications. This was consistent with the high plasticity of specialized metabolism in plants. CONCLUSIONS: This is the first characterization of the genetic control of abiotic stress-related specialized metabolites in sunflower. By providing hints concerning the importance of antioxidant molecules in this biological context, and by highlighting some of the potential molecular mechanisms underlying their biosynthesis, it could pave the way for novel applications in breeding. Although further analyses will be required to better understand this topic, studying how antioxidants contribute to the tolerance to abiotic stresses in sunflower appears as a promising area of research.

Sunflower Hybrids and Inbred Lines Adopt Different Physiological Strategies and Proteome Responses to Cope with Water Deficit.

Sunflower is a hybrid crop that is considered moderately drought-tolerant and adapted to new cropping systems required for the agro-ecological transition. Here, we studied the impact of hybridity status (hybrids vs. inbred lines) on the responses to drought at the molecular and eco-physiological level exploiting publicly available datasets. Eco-physiological traits and leaf proteomes were measured in eight inbred lines and their sixteen hybrids grown in the high-throughput phenotyping platform Phenotoul-Heliaphen. Hybrids and parental lines showed different growth strategies: hybrids grew faster in the absence of water constraint and arrested their growth more abruptly than inbred lines when subjected to water deficit. We identified 471 differentially accumulated proteins, of which 256 were regulated by drought. The amplitude of up- and downregulations was greater in hybrids than in inbred lines. Our results show that hybrids respond more strongly to water deficit at the molecular and eco-physiological levels. Because of presence/absence polymorphism, hybrids potentially contain more genes than their parental inbred lines. We propose that detrimental homozygous mutations and the lower number of genes in inbred lines lead to a constitutive defense mechanism that may explain the lower growth of inbred lines under well-watered conditions and their lower reactivity to water deficit.

Similar Transcriptomic Responses to Early and Late Drought Stresses Produce Divergent Phenotypes in Sunflower (Helianthus annuus L.).

Cultivated sunflower (Helianthus annuus L.) exhibits numerous phenotypic and transcriptomic responses to drought. However, the ways in which these responses vary with differences in drought timing and severity are insufficiently understood. We used phenotypic and transcriptomic data to evaluate the response of sunflower to drought scenarios of different timing and severity in a common garden experiment. Using a semi-automated outdoor high-throughput phenotyping platform, we grew six oilseed sunflower lines under control and drought conditions. Our results reveal that similar transcriptomic responses can have disparate phenotypic effects when triggered at different developmental time points. Leaf transcriptomic responses, however, share similarities despite timing and severity differences (e.g., 523 differentially expressed genes (DEGs) were shared across all treatments), though increased severity elicited greater differences in expression, particularly during vegetative growth. Across treatments, DEGs were highly enriched for genes related to photosynthesis and plastid maintenance. A co-expression analysis identified a single module (M8) enriched in all drought stress treatments. Genes related to drought, temperature, proline biosynthesis, and other stress responses were overrepresented in this module. In contrast to transcriptomic responses, phenotypic responses were largely divergent between early and late drought. Early-stressed sunflowers responded to drought with reduced overall growth, but became highly water-acquisitive during recovery irrigation, resulting in overcompensation (higher aboveground biomass and leaf area) and a greater overall shift in phenotypic correlations, whereas late-stressed sunflowers were smaller and more water use-efficient. Taken together, these results suggest that drought stress at an earlier growth stage elicits a change in development that enables greater uptake and transpiration of water during recovery, resulting in higher growth rates despite similar initial transcriptomic responses.

Exploring gene networks in two sunflower lines with contrasting leaf senescence phenotype using a system biology approach.

BACKGROUND: Leaf senescence is a complex process, controlled by multiple genetic and environmental variables. In sunflower, leaf senescence is triggered abruptly following anthesis thereby limiting the capacity of plants to keep their green leaf area during grain filling, which subsequently has a strong impact on crop yield. Recently, we performed a selection of contrasting sunflower inbred lines for the progress of leaf senescence through a physiological, cytological and molecular approach. Here we present a large scale transcriptomic analysis using RNA-seq and its integration with metabolic profiles for two contrasting sunflower inbred lines, R453 and B481-6 (early and delayed senescence respectively), with the aim of identifying metabolic pathways associated to leaf senescence. RESULTS: Gene expression profiles revealed a higher number of differentially expressed genes, as well as, higher expression levels in R453, providing evidence for early activation of the senescence program in this line. Metabolic pathways associated with sugars and nutrient recycling were differentially regulated between the lines. Additionally, we identified transcription factors acting as hubs in the co-expression networks; some previously reported as senescence-associated genes in model species but many are novel candidate genes. CONCLUSIONS: Understanding the onset and the progress of the senescence process in crops and the identification of these new candidate genes will likely prove highly useful for different management strategies to mitigate the impact of senescence on crop yield. Functional characterization of candidate genes will help to develop molecular tools for biotechnological applications in breeding crop yield.

Data describing the eco-physiological responses of twenty-four sunflower genotypes to water deficit.

This article presents experimental data describing the physiology and morphology of sunflower plants subjected to water deficit. Twenty-four sunflower genotypes were selected to represent genetic diversity within cultivated sunflower and included both inbred lines and their hybrids. Drought stress was applied to plants in pots at the vegetative stage using the high-throughput phenotyping platform Heliaphen at INRA Toulouse (France). Here, we provide data including specific leaf area, osmotic potential and adjustment, carbon isotope discrimination, leaf transpiration, plant architecture: plant height, leaf number, stem diameter. We also provide leaf areas of individual organs through time and growth rate during the stress period, environmental data such as temperatures, wind and radiation during the experiment. These data differentiate both treatment and the different genotypes and constitute a valuable resource to the community to study adaptation of crops to drought and the physiological basis of heterosis. It is available on the following repository: https://doi.org/10.25794/phenotype/er6lPW7V.

Comparison of GWAS models to identify non-additive genetic control of flowering time in sunflower hybrids.

KEY MESSAGE: This study compares five models of GWAS, to show the added value of non-additive modeling of allelic effects to identify genomic regions controlling flowering time of sunflower hybrids. Genome-wide association studies are a powerful and widely used tool to decipher the genetic control of complex traits. One of the main challenges for hybrid crops, such as maize or sunflower, is to model the hybrid vigor in the linear mixed models, considering the relatedness between individuals. Here, we compared two additive and three non-additive association models for their ability to identify genomic regions associated with flowering time in sunflower hybrids. A panel of 452 sunflower hybrids, corresponding to incomplete crossing between 36 male lines and 36 female lines, was phenotyped in five environments and genotyped for 2,204,423 SNPs. Intra-locus effects were estimated in multi-locus models to detect genomic regions associated with flowering time using the different models. Thirteen quantitative trait loci were identified in total, two with both model categories and one with only non-additive models. A quantitative trait loci on LG09, detected by both the additive and non-additive models, is located near a GAI homolog and is presented in detail. Overall, this study shows the added value of non-additive modeling of allelic effects for identifying genomic regions that control traits of interest and that could participate in the heterosis observed in hybrids.

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.

Genomic Prediction of Sunflower Hybrids Oil Content.

Prediction of hybrid performance using incomplete factorial mating designs is widely used in breeding programs including different heterotic groups. Based on the general combining ability (GCA) of the parents, predictions are accurate only if the genetic variance resulting from the specific combining ability is small and both parents have phenotyped descendants. Genomic selection (GS) can predict performance using a model trained on both phenotyped and genotyped hybrids that do not necessarily include all hybrid parents. Therefore, GS could overcome the issue of unknown parent GCA. Here, we compared the accuracy of classical GCA-based and genomic predictions for oil content of sunflower seeds using several GS models. Our study involved 452 sunflower hybrids from an incomplete factorial design of 36 female and 36 male lines. Re-sequencing of parental lines allowed to identify 468,194 non-redundant SNPs and to infer the hybrid genotypes. Oil content was observed in a multi-environment trial (MET) over 3 years, leading to nine different environments. We compared GCA-based model to different GS models including female and male genomic kinships with the addition of the female-by-male interaction genomic kinship, the use of functional knowledge as SNPs in genes of oil metabolic pathways, and with epistasis modeling. When both parents have descendants in the training set, the predictive ability was high even for GCA-based prediction, with an average MET value of 0.782. GS performed slightly better (+0.2%). Neither the inclusion of the female-by-male interaction, nor functional knowledge of oil metabolism, nor epistasis modeling improved the GS accuracy. GS greatly improved predictive ability when one or both parents were untested in the training set, increasing GCA-based predictive ability by 10.4% from 0.575 to 0.635 in the MET. In this scenario, performing GS only considering SNPs in oil metabolic pathways did not improve whole genome GS prediction but increased GCA-based prediction ability by 6.4%. Our results show that GS is a major improvement to breeding efficiency compared to the classical GCA modeling when either one or both parents are not well-characterized. This finding could therefore accelerate breeding through reducing phenotyping efforts and more effectively targeting for the most promising crosses.

The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution.

The domesticated sunflower, Helianthus annuus L., is a global oil crop that has promise for climate change adaptation, because it can maintain stable yields across a wide variety of environmental conditions, including drought. Even greater resilience is achievable through the mining of resistance alleles from compatible wild sunflower relatives, including numerous extremophile species. Here we report a high-quality reference for the sunflower genome (3.6 gigabases), together with extensive transcriptomic data from vegetative and floral organs. The genome mostly consists of highly similar, related sequences and required single-molecule real-time sequencing technologies for successful assembly. Genome analyses enabled the reconstruction of the evolutionary history of the Asterids, further establishing the existence of a whole-genome triplication at the base of the Asterids II clade and a sunflower-specific whole-genome duplication around 29 million years ago. An integrative approach combining quantitative genetics, expression and diversity data permitted development of comprehensive gene networks for two major breeding traits, flowering time and oil metabolism, and revealed new candidate genes in these networks. We found that the genomic architecture of flowering time has been shaped by the most recent whole-genome duplication, which suggests that ancient paralogues can remain in the same regulatory networks for dozens of millions of years. This genome represents a cornerstone for future research programs aiming to exploit genetic diversity to improve biotic and abiotic stress resistance and oil production, while also considering agricultural constraints and human nutritional needs.

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach of crop modelling and genome-wide association.

Understanding the genetic basis of phenotypic plasticity is crucial for predicting and managing climate change effects on wild plants and crops. Here, we combined crop modelling and quantitative genetics to study the genetic control of oil yield plasticity for multiple abiotic stresses in sunflower. First, we developed stress indicators to characterize 14 environments for three abiotic stresses (cold, drought and nitrogen) using the SUNFLO crop model and phenotypic variations of three commercial varieties. The computed plant stress indicators better explain yield variation than descriptors at the climatic or crop levels. In those environments, we observed oil yield of 317 sunflower hybrids and regressed it with three selected stress indicators. The slopes of cold stress norm reaction were used as plasticity phenotypes in the following genome-wide association study. Among the 65 534 tested Single Nucleotide Polymorphisms (SNPs), we identified nine quantitative trait loci controlling oil yield plasticity to cold stress. Associated single nucleotide polymorphisms are localized in genes previously shown to be involved in cold stress responses: oligopeptide transporters, lipid transfer protein, cystatin, alternative oxidase or root development. This novel approach opens new perspectives to identify genomic regions involved in genotype-by-environment interaction of a complex traits to multiple stresses in realistic natural or agronomical conditions.

RXLR and CRN Effectors from the Sunflower Downy Mildew Pathogen Plasmopara halstedii Induce Hypersensitive-Like Responses in Resistant Sunflower Lines.

Plasmopara halstedii is an obligate biotrophic oomycete causing downy mildew disease on sunflower, Helianthus annuus, an economically important oil crop. Severe symptoms of the disease (e.g., plant dwarfism, leaf bleaching, sporulation and production of infertile flower) strongly impair seed yield. Pl resistance genes conferring resistance to specific P. halstedii pathotypes were located on sunflower genetic map but yet not cloned. They are present in cultivated lines to protect them against downy mildew disease. Among the 16 different P. halstedii pathotypes recorded in France, pathotype 710 is frequently found, and therefore continuously controlled in sunflower by different Pl genes. High-throughput sequencing of cDNA from P. halstedii led us to identify potential effectors with the characteristic RXLR or CRN motifs described in other oomycetes. Expression of six P. halstedii putative effectors, five RXLR and one CRN, was analyzed by qRT-PCR in pathogen spores and in the pathogen infecting sunflower leaves and selected for functional analyses. We developed a new method for transient expression in sunflower plant leaves and showed for the first time subcellular localization of P. halstedii effectors fused to a fluorescent protein in sunflower leaf cells. Overexpression of the CRN and of 3 RXLR effectors induced hypersensitive-like cell death reactions in some sunflower near-isogenic lines resistant to pathotype 710 and not in susceptible corresponding lines, suggesting they could be involved in Pl loci-mediated resistances.

A biomarker based on gene expression indicates plant water status in controlled and natural environments.

Plant or soil water status is required in many scientific fields to understand plant responses to drought. Because the transcriptomic response to abiotic conditions, such as water deficit, reflects plant water status, genomic tools could be used to develop a new type of molecular biomarker. Using the sunflower (Helianthus annuus L.) as a model species to study the transcriptomic response to water deficit both in greenhouse and field conditions, we specifically identified three genes that showed an expression pattern highly correlated to plant water status as estimated by the pre-dawn leaf water potential, fraction of transpirable soil water, soil water content or fraction of total soil water in controlled conditions. We developed a generalized linear model to estimate these classical water status indicators from the expression levels of the three selected genes under controlled conditions. This estimation was independent of the four tested genotypes and the stage (pre- or post-flowering) of the plant. We further validated this gene expression biomarker under field conditions for four genotypes in three different trials, over a large range of water status, and we were able to correct their expression values for a large diurnal sampling period.