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.

Gene regulatory network inference methodology for genomic and transcriptomic data acquired in genetically related heterozygote individuals.

MOTIVATION: Inferring gene regulatory networks in non-independent genetically related panels is a methodological challenge. This hampers evolutionary and biological studies using heterozygote individuals such as in wild sunflower populations or cultivated hybrids. RESULTS: First, we simulated 100 datasets of gene expressions and polymorphisms, displaying the same gene expression distributions, heterozygosities and heritabilities as in our dataset including 173 genes and 353 genotypes measured in sunflower hybrids. Secondly, we performed a meta-analysis based on six inference methods [least absolute shrinkage and selection operator (Lasso), Random Forests, Bayesian Networks, Markov Random Fields, Ordinary Least Square and fast inference of networks from directed regulation (Findr)] and selected the minimal density networks for better accuracy with 64 edges connecting 79 genes and 0.35 area under precision and recall (AUPR) score on average. We identified that triangles and mutual edges are prone to errors in the inferred networks. Applied on classical datasets without heterozygotes, our strategy produced a 0.65 AUPR score for one dataset of the DREAM5 Systems Genetics Challenge. Finally, we applied our method to an experimental dataset from sunflower hybrids. We successfully inferred a network composed of 105 genes connected by 106 putative regulations with a major connected component. AVAILABILITY AND IMPLEMENTATION: Our inference methodology dedicated to genomic and transcriptomic data is available at https://forgemia.inra.fr/sunrise/inference_methods. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Optimized cultivar deployment improves the efficiency and stability of sunflower crop production at national scale.

KEY MESSAGE: Crop simulation helps to analyze environmental impacts on crops and provides year-independent context information. This information is of major importance when deciding which cultivar to choose at sowing time. Plant breeding programs design new crop cultivars which, while developed for distinct populations of environments, are nevertheless grown over large areas during their time in the market. Over its cultivation area, the crop is exposed to highly diverse stress patterns caused by climatic uncertainty and multiple management options, which often leads to decreased expected crop performance. In this study, we aim to assess how finer spatial management of genetic resources could reduce the yield variance explained by genotype × environment interactions in a set of cropping environments and ultimately improve the efficiency and stability of crop production. We used modeling and simulation to predict the crop performance resulting from the interaction between cultivar growth and development, climate and soil conditions, and management practices. We designed a computational experiment that evaluated the performance of a collection of commercial sunflower cultivars in a realistic population of cropping conditions in France, built from extensive agricultural surveys. Distinct farming locations sharing similar simulated abiotic stress patterns were clustered together to specify environment types. We then used optimization methods to search for cultivars × environments combinations leading to increased yield expectations. Results showed that a single cultivar choice adapted to the most frequent environment-type in the population is a robust strategy. However, the relevance of cultivar recommendations to specific locations was gradually increasing with the knowledge of pedo-climatic conditions. We argue that this approach while being operational on current genetic material could act synergistically with plant breeding as more diverse material could enable access to cultivars with distinctive traits, more adapted to specific conditions.

Genome-wide and comparative phylogenetic analysis of senescence-associated NAC transcription factors in sunflower (Helianthus annuus).

BACKGROUND: Leaf senescence delay impacts positively in grain yield by maintaining the photosynthetic area during the reproductive stage and during grain filling. Therefore a comprehensive understanding of the gene families associated with leaf senescence is essential. NAC transcription factors (TF) form a large plant-specific gene family involved in regulating development, senescence, and responses to biotic and abiotic stresses. The main goal of this work was to identify sunflower NAC TF (HaNAC) and their association with senescence, studying their orthologous to understand possible functional relationships between genes of different species. RESULTS: To clarify the orthologous relationships, we used an in-depth comparative study of four divergent taxa, in dicots and monocots, with completely sequenced genomes (Arabidopsis thaliana, Vitis vinifera, Musa acuminata and Oryza sativa). These orthologous groups provide a curated resource for large scale protein sequence annotation of NAC TF. From the 151 HaNAC genes detected in the latest version of the sunflower genome, 50 genes were associated with senescence traits. These genes showed significant differential expression in two contrasting lines according to an RNAseq assay. An assessment of overexpressing the Arabidopsis line for HaNAC001 (a gene of the same orthologous group of Arabidopsis thaliana ORE1) revealed that this line displayed a significantly higher number of senescent leaves and a pronounced change in development rate. CONCLUSIONS: This finding suggests HaNAC001 as an interesting candidate to explore the molecular regulation of senescence in sunflower.

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.

Metabolomic characterization of sunflower leaf allows discriminating genotype groups or stress levels with a minimal set of metabolic markers.

INTRODUCTION: Plant and crop metabolomic analyses may be used to study metabolism across genetic and environmental diversity. Complementary analytical strategies are useful for investigating metabolic changes and searching for biomarkers of response or performance. METHODS AND OBJECTIVES: The experimental material consisted in eight sunflower lines with two line status, four restorers (R, used as males) and four maintainers (B, corresponding to females) routinely used for sunflower hybrid varietal production, respectively to complement or maintain the cytoplasmic male sterility PET1. These lines were either irrigated at full soil capacity (WW) or submitted to drought stress (DS). Our aim was to combine targeted and non-targeted metabolomics to characterize sunflower leaf composition in order to investigate the effect of line status genotypes and environmental conditions and to find the best and smallest set of biomarkers for line status and stress response using a custom-made process of variables selection. RESULTS: Five hundred and eighty-eight metabolic variables were measured by using complementary analytical methods such as 1H-NMR, MS-based profiles and targeted analyses of major metabolites. Based on statistical analyses, a limited number of markers were able to separate WW and DS samples in a more discriminant manner than previously published physiological data. Another metabolic marker set was able to discriminate line status. CONCLUSION: This study underlines the potential of metabolic markers for discriminating genotype groups and environmental conditions. Their potential use for prediction is discussed.

Different genetic architectures underlie crop responses to the same pathogen: the {Helianthus annuus * Phoma macdonaldii} interaction case for black stem disease and premature ripening.

BACKGROUND: Phoma macdonaldii has been reported as the causal agent of black stem disease (BS) and premature ripening (PR) on sunflower. PR is considered as the most widespread and detrimental disease on sunflower in France. While genetic variability and QTL mapping for partial resistance of sunflower to stem, collar and roots attacks have been reported on plantlets in controlled conditions, this work aims to describe the genetic variability in a subset of a sunflower lines, and for the first time to map QTL involved in PR resistance evaluated in field conditions using controlled inoculation. RESULTS: An efficient and reliable method for inoculation used in field experiments induced stem base necrosis on up to 98% of all plants. A significant genetic variability for PR resistance in the field was detected among the 20 inbred lines of the core collection tested across the two years. For QTL mapping, the PR resistance evaluation was performed on two recombinant inbred lines (RIL) populations derived from the crosses XRQxPSC8 and FUxPAZ2 in two different years. QTL analyses were based on a newly developed consensus genetic map comprising 1007 non-redundant molecular markers. In each of the two RIL populations, different QTL involved in PR partial sunflower resistance were detected. The most significant QTL were detected 49 days post infection (DPI) on LG10 (LOD 7.7) and on LG7 (LOD 12.1) in the XRQxPSC8 and FUxPAZ2 RIL population, respectively. In addition, different QTL were detected on both populations for PR resistance measured between 14 and 35 DPI. In parallel, the incidence of natural attack of P. macdonaldii resulting in BS disease was recorded, showing that in these populations, the genetic of resistance to both diseases is not governed by the same factors. CONCLUSION: This work provides the first insights on the genetic architecture of sunflower PR resistance in the field. Moreover, the separate studies of symptoms on different organs and in time series allowed the identification of a succession of genetic components involved in the sunflower resistance to PR and BS diseases caused by Phoma macdonaldii along the development of the {plant * pathogen} interaction.

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.