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.

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.

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.

Spatial Genetic Structure and Pathogenic Race Composition at the Field Scale in the Sunflower Downy Mildew Pathogen, Plasmopara halstedii.

Yield losses in sunflower crops caused by Plasmopara halstedii can be up to 100%, depending on the cultivar susceptibility, environmental conditions, and virulence of the pathogen population. The aim of this study was to investigate the genetic and phenotypic structure of a sunflower downy mildew agent at the field scale. The genetic diversity of 250 P. halstedii isolates collected from one field in southern France was assessed using single-nucleotide polymorphisms (SNPs) and single sequence repeats (SSR). A total of 109 multilocus genotypes (MLG) were identified among the 250 isolates collected in the field. Four genotypes were repeated more than 20 times and spatially spread over the field. Estimates of genetic relationships among P. halstedii isolates using principal component analysis and a Bayesian clustering approach demonstrated that the isolates are grouped into two main genetic clusters. A high level of genetic differentiation among clusters was detected (FST = 0.35), indicating overall limited exchange between them, but our results also suggest that recombination between individuals of these groups is not rare. Genetic clusters were highly related to pathotypes, as previously described for this pathogen species. Eight different races were identified (100, 300, 304, 307, 703, 704, 707, and 714), with race 304 being predominant and present at most of the sites. The co-existence of multiple races at the field level is a new finding that could have important implications for the management of sunflower downy mildew. These data provide the first population-wide picture of the genetic structure of P. halstedii at a fine spatial scale.

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.