Association mapping and the genomic consequences of selection in sunflower

The genomic basis of nitrogen utilization efficiency and trait plasticity to improve nutrient stress tolerance in cultivated sunflower

Maintaining crop productivity is challenging as population growth, climate change, and increasing fertilizer costs necessitate expanding crop production to poorer lands whilst reducing inputs. Enhancing crops' nutrient use efficiency is thus an important goal, but requires a better understanding of related traits and their genetic basis. We investigated variation in low nutrient stress tolerance in a diverse panel of cultivated sunflower genotypes grown under high and low nutrient conditions, assessing relative growth rate (RGR) as performance. We assessed variation in traits related to nitrogen utilization efficiency (NUtE), mass allocation, and leaf elemental content. Across genotypes, nutrient limitation generally reduced RGR. Moreover, there was a negative correlation between vigor (RGR in control) and decline in RGR in response to stress. Given this trade-off, we focused on nutrient stress tolerance independent of vigor. This tolerance metric correlated with the change in NUtE, plasticity for a suite of morphological traits, and leaf element content. Genome-wide associations revealed regions associated with variation and plasticity in multiple traits, including two regions with seemingly additive effects on NUtE change. Our results demonstrate potential avenues for improving sunflower nutrient stress tolerance independent of vigor, and highlight specific traits and genomic regions that could play a role in enhancing tolerance.

Cultivating potential: Harnessing plant stem cells for agricultural crop improvement

Meristems are stem cell-containing structures that produce all plant organs and are therefore important targets for crop improvement. Developmental regulators control the balance and rate of cell divisions within the meristem. Altering these regulators impacts meristem architecture and, as a consequence, plant form. In this review, we discuss genes involved in regulating the shoot apical meristem, inflorescence meristem, axillary meristem, root apical meristem, and vascular cambium in plants. We highlight several examples showing how crop breeders have manipulated developmental regulators to modify meristem growth and alter crop traits such as inflorescence size and branching patterns. Plant transformation techniques are another innovation related to plant meristem research because they make crop genome engineering possible. We discuss recent advances on plant transformation made possible by studying genes controlling meristem development. Finally, we conclude with discussions about how meristem research can contribute to crop improvement in the coming decades.

Candidate pathway association and genome-wide association approaches reveal alternative genetic architectures of carotenoid content in cultivated sunflower (Helianthus annuus)

Premise

The explosion of available genomic data poses significant opportunities and challenges for genome-wide association studies. Current approaches via linear mixed models (LMM) are straightforward but prevent flexible assumptions of an a priori genomic architecture, while Bayesian sparse LMMs (BSLMMs) allow this flexibility. Complex traits, such as specialized metabolites, are subject to various hierarchical effects, including gene regulation, enzyme efficiency, and the availability of reactants.

Methods

To identify alternative genetic architectures, we examined the genetic architecture underlying the carotenoid content of an association mapping panel of Helianthus annuus individuals using multiple BSLMM and LMM frameworks.

Results

The LMMs of genome-wide single-nucleotide polymorphisms (SNPs) identified a single transcription factor responsible for the observed variations in the carotenoid content; however, a BSLMM of the SNPs with the bottom 1% of effect sizes from the results of the LMM identified multiple biologically relevant quantitative trait loci (QTLs) for carotenoid content external to the known (annotated) carotenoid pathway. A candidate pathway analysis (CPA) suggested a β-carotene isomerase to be the enzyme with the highest impact on the observed carotenoid content within the carotenoid pathway.

Discussion

While traditional LMM approaches suggested a single unknown transcription factor associated with carotenoid content variation in sunflower petals, BSLMM proposed several QTLs with interpretable biological relevance to this trait. In addition, the CPA allowed for the dissection of the regulatory vs. biosynthetic genetic architectures underlying this metabolic trait.

Genetic diversity in Sickleweed (Falcaria vulgaris) and using stepwise regression to identify marker associated with traits

One of the well-known medicinal plants in the Falcaria genus is Sickleweed. Falcaria species exhibit a high degree of genetic variability, posing challenges in the examination of genetic diversity due to the significant potential for hybridization and introgression among them. Utilizing morphological traits and molecular markers may prove to be a valuable approach in evaluating and harnessing germplasm, considering the current obstacles faced in breeding this medicinal herb. In 2021, fifteen Sickleweed populations were cultivated in pots under field conditions, employing a randomized complete block design with three replications. This aimed to assess genetic diversity and conduct marker-trait association analyses utilizing morpho-physiological characteristics and SSR markers. The Sickleweed populations displayed considerable genetic diversity across all traits. Through cluster analysis of traits and the utilization of the UPGMA method based on the Gower distance matrix, the population was classified into three distinct clusters. Upon examining all genotypes, 52 polymorphic bands were detected, with an average of 8.68 bands per primer. The average expected heterozygosity across all loci was 0.864, while the average PIC was 0.855. Molecular data analysis employing the Jaccard similarity index and UPGMA method revealed the division of Sickleweed populations into two major groups. Furthermore, the results of molecular variance analysis indicated that variation within the population exceeded that between populations. Thirty-two SSR fragments were found to be significantly associated with genomic regions controlling the studied traits, determined through the application of stepwise regression. Selection based on molecular markers offers a rapid method for breeding programs, with the genetic information obtained from these markers playing a crucial role. Therefore, alongside traits, selecting superior genotypes and populations of high value in breeding programs becomes feasible. The findings highlight that certain markers are linked to multiple traits, emphasizing the critical importance of this characteristic in plant breeding for the simultaneous improvement of numerous traits. The study's insights regarding markers hold potential for application in Sickleweed breeding programs.

Genome-Wide Association Mapping in Sunflower (Helianthus annuus) Reveals Common Loci and Putative Candidate Genes for Resistance to Diaporthe gulyae and D. helianthi Causing Phomopsis Stem Canker

Diaporthe gulyae and D. helianthi cause Phomopsis stem canker of sunflower (Helianthus annuus L.) in the United States. Because Phomopsis stem canker did not gain importance until the disease epidemic in 2010, limited studies were conducted to understand the genetic basis of sunflower resistance to D. gulyae and D. helianthi. The objectives of this study were to evaluate the United States Department of Agriculture cultivated accessions for resistance to D. gulyae and D. helianthi as well as to utilize genome-wide association studies (GWAS) to identify quantitative trait loci (QTLs) and putative candidate genes underlying those loci common to both organisms. For each fungus, 213 accessions were screened in a complete randomized design in the greenhouse and the experiment was repeated once. Six plants per accession were inoculated with a single isolate of D. gulyae or D. helianthi at four to six true leaves using the mycelium-contact inoculation method. At 15 days (D. gulyae) and 30 days (D. helianthi) postinoculation, accessions were evaluated for disease severity and compared with the susceptible confection inbred PI 552934. GWAS identified 28 QTLs common to the two fungi, and 24 genes overlapped close to these QTLs. Additionally, it was observed that the resistance QTLs derived mainly from landraces rather than from wild species. Seventeen putative candidate genes associated with resistance to D. gulyae or D. helianthi were identified that may be related to plant-pathogen interactions. These findings advanced our understanding of the genetic basis of resistance to D. gulyae and D. helianthi and will help develop resources for genomics-assisted breeding.

Association mapping for broomrape resistance in sunflower

INTRODUCTION

Sunflower breeding for resistance to the parasitic plant sunflower broomrape (Orobanche cumana Wallr.) requires the identification of novel resistance genes. In this research, we conducted a genome-wide association study (GWAS) to identify QTLs associated with broomrape resistance.

METHODS

The marker-trait associations were examined across a germplasm set composed of 104 sunflower accessions. They were genotyped with a 600k AXIOM® genome-wide array and evaluated for resistance to three populations of the parasite with varying levels of virulence (races EFR, FGV, and GTK) in two environments.

RESULTS AND DISCUSSION

The analysis of the genetic structure of the germplasm set revealed the presence of two main groups. The application of optimized treatments based on the general linear model (GLM) and the mixed linear model (MLM) allowed the detection of 14 SNP markers significantly associated with broomrape resistance. The highest number of marker-trait associations were identified on chromosome 3, clustered in two different genomic regions of this chromosome. Other associations were identified on chromosomes 5, 10, 13, and 16. Candidate genes for the main genomic regions associated with broomrape resistance were studied and discussed. Particularly, two significant SNPs on chromosome 3 associated with races EFR and FGV were found at two tightly linked SWEET sugar transporter genes. The results of this study have confirmed the role of some QTL on resistance to sunflower broomrape and have revealed new ones that may play an important role in the development of durable resistance to this parasitic weed in sunflower.

Genome-Wide Association Studies in Sunflower: Towards Sclerotinia sclerotiorum and Diaporthe/Phomopsis Resistance Breeding

Diseases caused by necrotrophic fungi, such as the cosmopolitan Sclerotinia sclerotiorum and the Diaporthe/Phomopsis complex, are among the most destructive diseases of sunflower worldwide. The lack of complete resistance combined with the inefficiency of chemical control makes assisted breeding the best strategy for disease control. In this work, we present an integrated genome-wide association (GWA) study investigating the response of a diverse panel of sunflower inbred lines to both pathogens. Phenotypic data for Sclerotinia head rot (SHR) consisted of five disease descriptors (disease incidence, DI; disease severity, DS; area under the disease progress curve for DI, AUDPCI, and DS, AUDPCS; and incubation period, IP). Two disease descriptors (DI and DS) were evaluated for two manifestations of Diaporthe/Phomopsis: Phomopsis stem canker (PSC) and Phomopsis head rot (PHR). In addition, a principal component (PC) analysis was used to derive transformed phenotypes as inputs to a univariate GWA (PC-GWA). Genotypic data comprised a panel of 4269 single nucleotide polymorphisms (SNP), generated via genotyping-by-sequencing. The GWA analysis revealed 24 unique marker-trait associations for SHR, 19 unique marker-trait associations for Diaporthe/Phomopsis diseases, and 7 markers associated with PC1 and PC2. No common markers were found for the response to the two pathogens. Nevertheless, epistatic interactions were identified between markers significantly associated with the response to S. sclerotiorum and Diaporthe/Phomopsis. This suggests that, while the main determinants of resistance may differ for the two pathogens, there could be an underlying common genetic basis. The exploration of regions physically close to the associated markers yielded 364 genes, of which 19 were predicted as putative disease resistance genes. This work presents the first simultaneous evaluation of two manifestations of Diaporthe/Phomopsis in sunflower, and undertakes a comprehensive GWA study by integrating PSC, PHR, and SHR data. The multiple regions identified, and their exploration to identify candidate genes, contribute not only to the understanding of the genetic basis of resistance, but also to the development of tools for assisted breeding.

HeliantHOME, a public and centralized database of phenotypic sunflower data

Genomic studies often attempt to link natural genetic variation with important phenotypic variation. To succeed, robust and reliable phenotypic data, as well as curated genomic assemblies, are required. Wild sunflowers, originally from North America, are adapted to diverse and often extreme environments and have historically been a widely used model plant system for the study of population genomics, adaptation, and speciation. Moreover, cultivated sunflower, domesticated from a wild relative (Helianthus annuus) is a global oil crop, ranking fourth in production of vegetable oils worldwide. Public availability of data resources both for the plant research community and for the associated agricultural sector, are extremely valuable. We have created HeliantHOME ( http://www.helianthome.org ), a curated, public, and interactive database of phenotypes including developmental, structural and environmental ones, obtained from a large collection of both wild and cultivated sunflower individuals. Additionally, the database is enriched with external genomic data and results of genome-wide association studies. Finally, being a community open-source platform, HeliantHOME is expected to expand as new knowledge and resources become available.

Genomic regions associate with major axes of variation driven by gas exchange and leaf construction traits in cultivated sunflower (Helianthus annuus L.)

Stomata and leaf veins play an essential role in transpiration and the movement of water throughout leaves. These traits are thus thought to play a key role in the adaptation of plants to drought and a better understanding of the genetic basis of their variation and coordination could inform efforts to improve drought tolerance. Here, we explore patterns of variation and covariation in leaf anatomical traits and analyze their genetic architecture via genome-wide association (GWA) analyses in cultivated sunflower (Helianthus annuus L.). Traits related to stomatal density and morphology as well as lower-order veins were manually measured from digital images while the density of minor veins was estimated using a novel deep learning approach. Leaf, stomatal, and vein traits exhibited numerous significant correlations that generally followed expectations based on functional relationships. Correlated suites of traits could further be separated along three major principal component (PC) axes that were heavily influenced by variation in traits related to gas exchange, leaf hydraulics, and leaf construction. While there was limited evidence of colocalization when individual traits were subjected to GWA analyses, major multivariate PC axes that were most strongly influenced by several traits related to gas exchange or leaf construction did exhibit significant genomic associations. These results provide insight into the genetic basis of leaf trait covariation and showcase potential targets for future efforts aimed at modifying leaf anatomical traits in sunflower.

Combined GWAS and Transcriptome Analyses Provide New Insights Into the Response Mechanisms of Sunflower Against Drought Stress

Sunflower is one of the most important oil crops in the world, and drought stress can severely limit its production and quality. To understand the underlying mechanism of drought tolerance, and identify candidate genes for drought tolerance breeding, we conducted a combined genome-wide association studies (GWAS) and RNA-seq analysis. A total of 226 sunflower inbred lines were collected from different regions of China and other countries. Eight phenotypic traits were evaluated under control and drought stress conditions. Genotyping was performed using a Specific-Locus Amplified Fragment Sequencing (SLAF-seq) approach. A total of 934.08 M paired-end reads were generated, with an average Q30 of 91.97%. Based on the 243,291 polymorphic SLAF tags, a total of 94,162 high-quality SNPs were identified. Subsequent analysis of linkage disequilibrium (LD) and population structure in the 226 accessions was carried out based on the 94,162 high-quality SNPs. The average LD decay across the genome was 20 kb. Admixture analysis indicated that the entire population most likely originated from 11 ancestors. GWAS was performed using three methods (MLM, FarmCPU, and BLINK) simultaneously. A total of 80 SNPs showed significant associations with the 8 traits (p < 1.062 × 10-6). Next, a total of 118 candidate genes were found. To obtain more reliable candidate genes, RNA-seq analysis was subsequently performed. An inbred line with the highest drought tolerance was selected according to phenotypic traits. RNA was extracted from leaves at 0, 7, and 14 days of drought treatment. A total of 18,922 differentially expressed genes were obtained. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis showed up-regulated genes were mainly enriched in the branched-chain amino acid catabolic process, while the down-regulated genes were mainly enriched in the photosynthesis-related process. Six DEGs were randomly selected from all DEGs for validation; these genes showed similar patterns in RNA-seq and RT-qPCR analysis, with a correlation coefficient of 0.8167. Through the integration of the genome-wide association study and the RNA-sequencing, 14 candidate genes were identified. Four of them (LOC110885273, LOC110872899, LOC110891369, LOC110920644) were abscisic acid related protein kinases and transcription factors. These genes may play an important role in sunflower drought response and will be used for further study. Our findings provide new insights into the response mechanisms of sunflowers against drought stress and contribute to further genetic breeding.

Expression complementation of gene presence/absence polymorphisms in hybrids contributes importantly to heterosis in sunflower

INTRODUCTION

Numerous crops have transitioned to hybrid seed production to increase yields and yield stability through heterosis. However, the molecular mechanisms underlying heterosis and its stability across environments are not yet fully understood.

OBJECTIVES

This study aimed to (1) elucidate the genetic and molecular mechanisms underlying heterosis in sunflower, and (2) determine how heterosis is maintained under different environments.

METHODS

Genome-wide association (GWA) analyses were employed to assess the effects of presence/absence variants (PAVs) and stop codons on 16 traits phenotyped in the sunflower association mapping population at three locations. To link the GWA results to transcriptomic variation, we sequenced the transcriptomes of two sunflower cultivars and their F₁ hybrid (INEDI) under both control and drought conditions and analyzed patterns of gene expression and alternative splicing.

RESULTS

Thousands of PAVs were found to affect phenotypic variation using a relaxed significance threshold, and at most such loci the "absence" allele reduced values of heterotic traits, but not those of non-heterotic traits. This pattern was strengthened for PAVs that showed expression complementation in INEDI. Stop codons were much rarer than PAVs and less likely to reduce heterotic trait values. Hybrid expression patterns were enriched for the GO category, sensitivity to stimulus, but all genotypes responded to drought similarily - by up-regulating water stress response pathways and down-regulating metabolic pathways. Changes in alternative splicing were strongly negatively correlated with expression variation, implying that alternative splicing in this system largely acts to reinforce expression responses.

CONCLUSION

Our results imply that complementation of expression of PAVs in hybrids is a major contributor to heterosis in sunflower, consistent with the dominance model of heterosis. This mechanism can account for yield stability across different environments. Moreover, given the much larger numbers of PAVs in plant vs. animal genomes, it also offers an explanation for the stronger heterotic responses seen in the former.

Linkage disequilibrium and population structure in a core collection of Brassica napus (L.)

Estimation of genetic diversity in rapeseed is important for sustainable breeding program to provide an option for the development of new breeding lines. The objective of this study was to elucidate the patterns of genetic diversity within and among different structural groups, and measure the extent of linkage disequilibrium (LD) of 383 globally distributed rapeseed germplasm using 8,502 single nucleotide polymorphism (SNP) markers. We divided the germplasm collection into five subpopulations (P1 to P5) according to geographic and growth habit-related patterns. All subpopulations showed moderate genetic diversity (average H = 0.22 and I = 0.34). The pairwise Fst comparison revealed a great degree of divergence (Fst > 0.24) between most of the combinations. The rutabaga type showed highest divergence with spring and winter types. Higher divergence was also found between winter and spring types. Admixture model based structure analysis, principal component and neighbor-joining tree analysis placed all subpopulations into three distinct clusters. Admixed genotype constituted 29.24% of total genotypes, while remaining 70.76% belongs to identified clusters. Overall, mean linkage disequilibrium was 0.03 and it decayed to its half maximum within < 45 kb distance for whole genome. The LD decay was slower in C genome (< 93 kb); relative to the A genome (< 21 kb) which was confirmed by availability of larger haplotype blocks in C genome than A genome. The findings regarding LD pattern and population structure will help to utilize the collection as an important resource for association mapping efforts to identify genes useful in crop improvement as well as for selection of parents for hybrid breeding.

Genetic basis and dual adaptive role of floral pigmentation in sunflowers

Variation in floral displays, both between and within species, has been long known to be shaped by the mutualistic interactions that plants establish with their pollinators. However, increasing evidence suggests that abiotic selection pressures influence floral diversity as well. Here, we analyse the genetic and environmental factors that underlie patterns of floral pigmentation in wild sunflowers. While sunflower inflorescences appear invariably yellow to the human eye, they display extreme diversity for patterns of ultraviolet pigmentation, which are visible to most pollinators. We show that this diversity is largely controlled by cis-regulatory variation affecting a single MYB transcription factor, HaMYB111, through accumulation of ultraviolet (UV)-absorbing flavonol glycosides in ligules (the ‘petals’ of sunflower inflorescences). Different patterns of ultraviolet pigments in flowers are strongly correlated with pollinator preferences. Furthermore, variation for floral ultraviolet patterns is associated with environmental variables, especially relative humidity, across populations of wild sunflowers. Ligules with larger ultraviolet patterns, which are found in drier environments, show increased resistance to desiccation, suggesting a role in reducing water loss. The dual role of floral UV patterns in pollinator attraction and abiotic response reveals the complex adaptive balance underlying the evolution of floral traits.

Flowers are an important part of how many plants reproduce. Their distinctive colours, shapes and patterns attract specific pollinators, but they can also help to protect the plant from predators and environmental stresses.

Many flowers contain pigments that absorb ultraviolet (UV) light to display distinct UV patterns – although invisible to the human eye, most pollinators are able to see them. For example, when seen in UV, sunflowers feature a ‘bullseye’ with a dark centre surrounded by a reflective outer ring. The sizes and thicknesses of these rings vary a lot within and between flower species, and so far, it has been unclear what causes this variation and how it affects the plants.

To find out more, Todesco et al. studied the UV patterns in various wild sunflowers across North America by considering the ecology and molecular biology of different plants. This revealed great variation between the UV patterns of the different sunflower populations. Moreover, Todesco et al. found that a gene called HaMYB111 is responsible for the diverse UV patterns in the sunflowers. This gene controls how plants make chemicals called flavonols that absorb UV light.

Flavonols also help to protect plants from damage caused by droughts and extreme temperatures. Todesco et al. showed that plants with larger bullseyes had more flavonols, attracted more pollinators, and were better at conserving water. Accordingly, these plants were found in drier locations.

This study suggests that, at least in sunflowers, UV patterns help both to attract pollinators and to control water loss. These insights could help to improve pollination – and consequently yield – in cultivated plants, and to develop plants with better resistance to extreme weather. This work also highlights the importance of combining biology on small and large scales to understand complex processes, such as adaptation and evolution.

A Correlation Study on In Vitro Physiological Activities of Soybean Cultivars, 19 Individual Isoflavone Derivatives, and Genetic Characteristics

The functionality of soybeans is an important factor in the selection and utilization of excellent soybean cultivars, and isoflavones are representative functional substances in soybeans, which exhibit effects on antioxidants, estrogen activity, and cancer, and prevent cardiovascular diseases. This study analyzed ABTS, DPPH, estrogen, ER (ER) alpha, UCP-1, and NO inhibition activities in 48 types of soybean cultivars, as well as the relationship with 19 isolated types of individual isoflavone derivatives. Statistical analysis was conducted to find individual isoflavone derivatives affecting physiological activities, revealing the high correlation of three types of derivatives: genistein 7-O-(6″-O-acetyl)glucoside (6″-O-acetylgenistin), genistein 7-O-(2″-O-apiosyl)glucoside, and glycitein. Based on these results, 15 types of soybean cultivars were selected (one control type, seven yellow types, six black types, and one green type), which have both high physiological activities and a high content of individual isoflavone derivatives. In addition, these high correlations were further verified through a genome-wide association study (GWAS) to determine the association between activities, substances, and genetic characteristics. This study comprehensively describes the relationship between the specific physiological activities of soybean resources, individual isoflavone derivative substances, and SNPs, which will be utilized for in-depth research, such as selection of excellent soybean resources with specific physiological activities.

Genetic control of arbuscular mycorrhizal colonization by Rhizophagus intraradices in Helianthus annuus (L.)

Plant symbiosis with arbuscular mycorrhizal (AM) fungi provides many benefits, including increased nutrient uptake, drought tolerance, and belowground pathogen resistance. To develop a better understanding of the genetic architecture of mycorrhizal symbiosis, we conducted a genome-wide association study (GWAS) of this plant-fungal interaction in cultivated sunflower. A diversity panel of cultivated sunflower (Helianthus annuus L.) was phenotyped for root colonization under inoculation with the AM fungus Rhizophagus intraradices. Using a mixed linear model approach with a high-density genetic map, we identified genomic regions that are likely associated with R. intraradices colonization in sunflower. Additionally, we used a set of twelve diverse lines to assess the effect that inoculation with R. intraradices has on dried shoot biomass and macronutrient uptake. Colonization among lines in the mapping panel ranged from 0-70% and was not correlated with mycorrhizal growth response, shoot phosphorus response, or shoot potassium response among the Core 12 lines. Association mapping yielded three single-nucleotide polymorphisms (SNPs) that were significantly associated with R. intraradices colonization. This is the first study to use GWAS to identify genomic regions associated with AM colonization in an Asterid eudicot species. Three genes of interest identified from the regions containing these SNPs are likely related to plant defense.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

Genotyping and lipid profiling of 601 cultivated sunflower lines reveals novel genetic determinants of oil fatty acid content

BACKGROUND

Sunflower is an important oilseed crop domesticated in North America approximately 4000 years ago. During the last century, oil content in sunflower was under strong selection. Further improvement of oil properties achieved by modulating its fatty acid composition is one of the main directions in modern oilseed crop breeding.

RESULTS

We searched for the genetic basis of fatty acid content variation by genotyping 601 inbred sunflower lines and assessing their lipid and fatty acid composition. Our genome-wide association analysis based on the genotypes for 15,483 SNPs and the concentrations of 23 fatty acids, including minor fatty acids, revealed significant genetic associations for eleven of them. Identified genomic regions included the loci involved in rare fatty acids variation on chromosomes 3 and 14, explaining up to 34.5% of the total variation of docosanoic acid (22:0) in sunflower oil.

CONCLUSIONS

This is the first large scale implementation of high-throughput lipidomic profiling to sunflower germplasm characterization. This study contributes to the genetic characterization of Russian sunflower collections, which made a substantial contribution to the development of sunflower as the oilseed crop worldwide, and provides new insights into the genetic control of oil composition that can be implemented in future studies.

Detection of breeding signatures in wheat using a linkage disequilibrium-corrected mapping approach

Marker assisted breeding, facilitated by reference genome assemblies, can help to produce cultivars adapted to changing environmental conditions. However, anomalous linkage disequilibrium (LD), where single markers show high LD with markers on other chromosomes but low LD with adjacent markers, is a serious impediment for genetic studies. We used a LD-correction approach to overcome these drawbacks, correcting the physical position of markers derived from 15 and 135 K arrays in a diversity panel of bread wheat representing 50 years of breeding history. We detected putative mismapping of 11.7% markers and improved the physical alignment of 5.4% markers. Population analysis indicated reduced genetic diversity over time as a result of breeding efforts. By analysis of outlier loci and allele frequency change over time we traced back the 2NS/2AS translocation of Aegilops ventricosa to one cultivar, "Cardos" (registered in 1998) which was the first among the panel to contain this translocation. A "selective sweep" for this important translocation region on chromosome 2AS was found, putatively linked to plant response to biotic stress factors. Our approach helps in overcoming the drawbacks of incorrectly anchored markers on the wheat reference assembly and facilitates detection of selective sweeps for important agronomic traits.

Sunflower inflorescences absorb maximum light energy if they face east and afternoons are cloudier than mornings

The mature inflorescence of sunflowers (Helianthus annuus) orients eastward after its anthesis (the flowering period, especially the maturing of the stamens), from which point it no longer tracks the Sun. Although several hypothetical explanations have been proposed for the ecological functions of this east facing, none have been tested. Here we propose an atmospheric-optical explanation. Using (i) astronomical data of the celestial motion of the Sun, (ii) meteorological data of diurnal cloudiness for Boone County located in the region from which domesticated sunflowers originate, (iii) time-dependent elevation angle of mature sunflower heads, and (iv) absorption spectra of the inflorescence and the back of heads, we computed the light energy absorbed separately by the inflorescence and the back between anthesis and senescence. We found that the inflorescences facing east absorb the maximum radiation, being advantageous for seed production and maturation, furthermore west facing would be more advantageous than south facing. The reason for these is that afternoons are cloudier than mornings in the cultivation areas of sunflowers. Since the photosynthesizing green back of mature heads absorbs maximal energy when the inflorescence faces west, maximizing the energy absorbed by the back cannot explain the east facing of inflorescences. The same results were obtained for central Italy and Hungary, where mornings are also less cloudy than afternoons. In contrast, in south Sweden, where mornings are cloudier than afternoons, west-facing mature inflorescences would absorb the maximum light energy. We suggest that the domesticated Helianthus annuus developed an easterly final orientation of its mature inflorescence, because it evolved in a region with cloudier afternoons.

Key Traits and Genes Associate with Salinity Tolerance Independent from Vigor in Cultivated Sunflower

With rising food demands, crop production on salinized lands is increasingly necessary. Sunflower (Helianthus annuus), a moderately salt-tolerant crop, exhibits a tradeoff where more vigorous, high-performing genotypes have a greater proportional decline in biomass under salinity stress. Prior research has found deviations from this relationship across genotypes. Here, we identified the traits and genomic regions underlying variation in this expectation-deviation tolerance (the magnitude and direction of deviations from the expected effect of salinity). We grew a sunflower diversity panel under control and salt-stressed conditions and measured a suite of morphological (growth, mass allocation, plant and leaf morphology) and leaf ionomic traits. The genetic basis of variation and plasticity in these traits was investigated via genome-wide association, which also enabled the identification of genomic regions (i.e. haplotypic blocks) influencing multiple traits. We found that the magnitude and direction of plasticity in whole-root mass fraction, fine root mass fraction, and chlorophyll content, as well as leaf sodium and potassium content under saline conditions, were most strongly correlated with expectation-deviation tolerance. We identified multiple genomic regions underlying these traits as well as a single alpha-mannosidase gene directly associated with this tolerance metric. Our results show that, by taking the vigor-salinity effect tradeoff into account, we can identify unique traits and genes associated with salinity tolerance. Since these traits and genomic regions are distinct from those associated with high vigor (i.e. growth in benign conditions), they provide an avenue for increasing salinity tolerance in high-performing sunflower genotypes without compromising vigor.

Unveiling the genetic basis of Sclerotinia head rot resistance in sunflower

BACKGROUND

Sclerotinia sclerotiorum is a necrotrophic fungus that causes Sclerotinia head rot (SHR) in sunflower, with epidemics leading to severe yield losses. In this work, we present an association mapping (AM) approach to investigate the genetic basis of natural resistance to SHR in cultivated sunflower, the fourth most widely grown oilseed crop in the world.

RESULTS

Our association mapping population (AMP), which comprises 135 inbred breeding lines (ILs), was genotyped using 27 candidate genes, a panel of 9 Simple Sequence Repeat (SSR) markers previously associated with SHR resistance via bi-parental mapping, and a set of 384 SNPs located in genes with molecular functions related to stress responses. Moreover, given the complexity of the trait, we evaluated four disease descriptors (i.e, disease incidence, disease severity, area under the disease progress curve for disease incidence, and incubation period). As a result, this work constitutes the most exhaustive AM study of disease resistance in sunflower performed to date. Mixed linear models accounting for population structure and kinship relatedness were used for the statistical analysis of phenotype-genotype associations, allowing the identification of 13 markers associated with disease reduction. The number of favourable alleles was negatively correlated to disease incidence, disease severity and area under the disease progress curve for disease incidence, whereas it was positevily correlated to the incubation period.

CONCLUSIONS

Four of the markers identified here as associated with SHR resistance (HA1848, HaCOI_1, G33 and G34) validate previous research, while other four novel markers (SNP117, SNP136, SNP44, SNP128) were consistently associated with SHR resistance, emerging as promising candidates for marker-assisted breeding. From the germplasm point of view, the five ILs carrying the largest combination of resistance alleles provide a valuable resource for sunflower breeding programs worldwide.

Genetic Analysis of Domestication Parallels in Annual and Perennial Sunflowers (Helianthus spp.): Routes to Crop Development

Parallels exist between the domestication of new species and the improvement of various crops through selection on traits which favor the sowing, harvest and retention of yield potential and the directed efforts to improve their agronomics, disease resistance and quality characteristics. Common selection pressures may result in the parallel selection of orthologs underlying these traits and homologies between crop species can be exploited by plant breeders to improve germplasm. Perennial grains and oilseeds are a class of proposed crops for improving the diversity and sustainability of agricultural systems. Maximilian sunflower (Helianthus maximiliani Schrad.) is a perennial crop wild relative of sunflower (Helianthus annuus L.) and a candidate perennial oilseed species. Understanding parallels between cultivated H. annuus and H. maximiliani may provide new tools for the development of Maximilian sunflower and other wild relatives of sunflower as crops to enhance functional diversity in cropping systems. F₂ populations of Maximilian sunflower segregating for traits associated with the domestication ideotype of cultivated sunflower including branching architecture, capitulum morphology and flowering time were developed to investigate parallels between H. maximiliani and H. annuus. Genotype-by-sequencing (GBS) was employed to genotype novel Maximilian sunflower populations and perform quantitative-trait-loci (QTL) analysis. A total of 11 QTL in five regions were identified across 21 linkage groups using 4142 GBS derived single nucleotide polymorphism markers called using the sunflower reference genome as a guide. A major QTL on linkage group 17b, associated with aspects of floral development and apical dominance, was discovered and corresponds with a known domestication QTL hotspot in H. annuus and candidate genes were identified. This suggests the potential to exploit orthologs for neo-domestication of H. maximiliani for traits such as branching architecture, timing of anthesis, and capitulum size and morphology for the development of a perennial oilseed crop from wild relatives of cultivated sunflower.

Genome-Wide Association Mapping of Floral Traits in Cultivated Sunflower (Helianthus annuus)

Floral morphology and pigmentation are both charismatic and economically relevant traits associated with cultivated sunflower (Helianthus annuus L.). Recent work has linked floral morphology and pigmentation to pollinator efficiency and seed yield. Understanding the genetic architecture of such traits is essential for crop improvement, and gives insight into the role of genetic constraints in shaping floral diversity. A diversity panel of 288 sunflower genotypes was phenotyped for a variety of morphological, phenological, and color traits in both a greenhouse and a field setting. Association mapping was performed using 5788 SNP markers using a mixed linear model approach. Several dozen markers across 10 linkage groups were significantly associated with variation in morphological and color trait variation. Substantial trait plasticity was observed between greenhouse and field phenotyping, and associations differed between environments. Color traits mapped more strongly than morphology in both settings, with markers together explaining 16% of petal carotenoid content in the greenhouse, and 17% and 24% of variation in disc anthocyanin presence in the field and greenhouse, respectively. Morphological traits like disc size mapped more strongly in the field, with markers together explaining up to 19% of disc size variation. Loci identified here through association mapping within cultivated germplasm differ from those identified through biparental crosses between modern cultivated sunflower and either its wild progenitor or domesticated landraces. Several loci lie within genomic regions involved in domestication. Differences between phenotype expression under greenhouse and field conditions highlight the importance of plasticity in determining floral morphology and pigmentation.

Inheritance and molecular mapping of powdery mildew (Golovinomyces orontii) resistance gene(s) in sunflower (Helianthus annuus L.)

Sources of resistance to powdery mildew incited by Golovinomyces orontii have been identified in wild sunflowers and few exotic lines. The present investigation has been undertaken to study the inheritance of powdery mildew resistance and to map the quantitative trait loci (QTLs) governing resistance to powdery mildew in a multiple disease resistance line, TX16R (PI 642072). The inheritance was observed as a continuous distribution in a set of 264 F₂ population and 93 recombinant inbred lines (RILs) of a cross between a highly susceptible accession PS 2023 and TX16R. Screening of the two population sets was done with 484 sunflower-specific SSR primers of which 175 primers showed polymorphism between the parents. Based on the phenotyping and genotyping data, the linkage map was constructed with 93 RILs. The map spanned 1200 cM and included 64 markers distributed along the 17 sunflower chromosomes in the haploid set. Quantitative trait loci (QTL) analysis identified three genomic regions for resistance to powdery mildew, two of which mapped on chromosome 10 and one on chromosome 5. This is the first report on mapping of powdery mildew resistance in sunflower and paves the way in fine mapping and introgression of resistance for powdery mildew in sunflower through marker-assisted breeding.

Ancestral Reconstruction of Karyotypes Reveals an Exceptional Rate of Nonrandom Chromosomal Evolution in Sunflower

Mapping the chromosomal rearrangements between species can inform our understanding of genome evolution, reproductive isolation, and speciation. Here, we present a novel algorithm for identifying regions of synteny in pairs of genetic maps, which is implemented in the accompanying R package syntR. The syntR algorithm performs as well as previous ad hoc methods while being systematic, repeatable, and applicable to mapping chromosomal rearrangements in any group of species. In addition, we present a systematic survey of chromosomal rearrangements in the annual sunflowers, which is a group known for extreme karyotypic diversity. We build high-density genetic maps for two subspecies of the prairie sunflower, Helianthus petiolaris ssp. petiolaris and H. petiolaris ssp. fallax Using syntR, we identify blocks of synteny between these two subspecies and previously published high-density genetic maps. We reconstruct ancestral karyotypes for annual sunflowers using those synteny blocks and conservatively estimate that there have been 7.9 chromosomal rearrangements per million years, a high rate of chromosomal evolution. Although the rate of inversion is even higher than the rate of translocation in this group, we further find that every extant karyotype is distinguished by between one and three translocations involving only 8 of the 17 chromosomes. This nonrandom exchange suggests that specific chromosomes are prone to translocation and may thus contribute disproportionately to widespread hybrid sterility in sunflowers. These data deepen our understanding of chromosome evolution and confirm that Helianthus has an exceptional rate of chromosomal rearrangement that may facilitate similarly rapid diversification.

Genetic Diversity, Population Structure and Linkage Disequilibrium Assessment among International Sunflower Breeding Collections

Sunflower germplasm collections are valuable resources for broadening the genetic base of commercial hybrids and ameliorate the risk of climate events. Nowadays, the most studied worldwide sunflower pre-breeding collections belong to INTA (Argentina), INRA (France), and USDA-UBC (United States of America-Canada). In this work, we assess the amount and distribution of genetic diversity (GD) available within and between these collections to estimate the distribution pattern of global diversity. A mixed genotyping strategy was implemented, by combining proprietary genotyping-by-sequencing data with public whole-genome-sequencing data, to generate an integrative 11,834-common single nucleotide polymorphism matrix including the three breeding collections. In general, the GD estimates obtained were moderate. An analysis of molecular variance provided evidence of population structure between breeding collections. However, the optimal number of subpopulations, studied via discriminant analysis of principal components (K = 12), the bayesian STRUCTURE algorithm (K = 6) and distance-based methods (K = 9) remains unclear, since no single unifying characteristic is apparent for any of the inferred groups. Different overall patterns of linkage disequilibrium (LD) were observed across chromosomes, with Chr10, Chr17, Chr5, and Chr2 showing the highest LD. This work represents the largest and most comprehensive inter-breeding collection analysis of genomic diversity for cultivated sunflower conducted to date.

Phylogeography and the Evolutionary History of Sunflower (Helianthus annuus L.): Wild Diversity and the Dynamics of Domestication

Patterns of genetic variation in crops are the result of selection and demographic changes that occurred during their domestication and improvement. In many cases, we have an incomplete picture of the origin of crops in the context of their wild progenitors, particularly with regard to the processes producing observed levels of standing genetic variation. Here, we analyzed sequence diversity in cultivated sunflower (Helianthus annuus L.) and its wild progenitor (common sunflower, also H. annuus) to reconstruct phylogeographic relationships and population genetic/demographic patterns across sunflower. In common sunflower, south-north patterns in the distribution of nucleotide diversity and lineage splitting indicate a history of rapid postglacial range expansion from southern refugia. Cultivated sunflower accessions formed a clade, nested among wild populations from the Great Plains, confirming a single domestication event in central North America. Furthermore, cultivated accessions sorted by market type (i.e., oilseed vs. confectionery) rather than breeding pool, recapitulating the secondary development of oil-rich cultivars during its breeding history. Across sunflower, estimates of nucleotide diversity and effective population sizes suggest that cultivated sunflower underwent significant population bottlenecks following its establishment ~5000 years ago. The patterns inferred here corroborate those from previous studies of sunflower domestication, and provide a comprehensive overview of its evolutionary history.

Seed and floret size parameters of sunflower are determined by partially overlapping sets of quantitative trait loci with epistatic interactions

KEY MESSAGE

Floret and seed traits are moderately correlated phenotypically in modern sunflower cultivars, but the underlying genetics are mostly independent. Seed traits in particular are governed in part by epistatic effects among quantitative trait loci. Seed size is an important quality component in marketing commercial sunflower (Helianthus annuus L.), particularly for the in-shell confectionery market, where long and broad seed types are preferred as a directly consumed snack food globally. Floret size is also important because corolla tube length was previously shown to be inversely correlated with pollinator visitation, impacting bee foraging potential and pollinator services to the plant. Commercial sunflower production benefits from pollinator visits, despite being self-compatible, and bees are required in hybrid seed production, where "female" and "male" inbred lines are crossed at field scale. Issues with pollination of long-seed confectionery sunflower suggest that there may be an unfavorable correlation between seed and floret traits; thus, our objective was to determine the strength of the correlation between seed and floret traits, and confirm any co-localization of seed and floret trait loci using genome-wide association analysis in the SAM diversity panel of sunflower. Our results indicate that phenotypic correlations between seed and floret traits are generally low to moderate, regardless of market class, a component of population substructure. Association mapping results mirror the correlations: while a few loci overlap, many loci for the two traits are not overlapping or even adjacent. The genetics of these traits, while modestly quantitative and influenced by epistatic effects, are not a barrier to simultaneous improvement of seed length and pollinator-friendly floret traits. We conclude that breeding for large seed size, which is required for the confectionery seed market, is possible without producing florets too long for efficient use by pollinators, which promotes bee foraging and associated pollination services.

Genetic and phenotypic analyses indicate that resistance to flooding stress is uncoupled from performance in cultivated sunflower

Given the rising risk of extreme weather caused by climate change, enhancement of abiotic stress resistance in crops is increasingly urgent. But will the development of stress-resistant cultivars come at the cost of yield under ideal conditions? We hypothesize that this need not be inevitable, because resistance alleles with minimal pleiotropic costs may evade artificial selection and be retained in crop germplasm. Genome-wide association (GWA) analyses for variation in plant performance and flooding response were conducted in cultivated sunflower, a globally important oilseed. We observed broad variation in flooding responses among genotypes. Flooding resistance was not strongly correlated with performance in control conditions, suggesting no inherent trade-offs. Consistent with this finding, we identified a subset of loci conferring flooding resistance, but lacking antagonistic effects on growth. Genetic diversity loss at candidate genes underlying these loci was significantly less than for other resistance genes during cultivated sunflower evolution. Despite bottlenecks associated with domestication and improvement, low-cost resistance alleles remain within the cultivated sunflower gene pool. Thus, development of cultivars that are both flooding-tolerant and highly productive should be straightforward. Results further indicate that estimates of pleiotropic costs from GWA analyses explain, in part, patterns of diversity loss in crop genomes.

Development and Validation of Markers for the Fertility Restorer Gene Rf1 in Sunflower

Hybrid breeding in sunflowers based on CMS PET1 requires development of restorer lines carrying, in most cases, the restorer gene Rf1. Markers for marker-assisted selection have been developed, but there is still need for closer, more versatile, and co-dominant markers linked to Rf1. Homology searches against the reference sunflower genome using sequences of cloned markers, as well as Bacterial Artificial Chromosome (BAC)-end sequences of clones hybridizing to them, allowed the identification of two genomic regions of 30 and 3.9 Mb, respectively, as possible physical locations of the restorer gene Rf1 on linkage group 13. Nine potential candidate genes, encoding six pentatricopeptide repeat proteins, one tetratricopeptide-like helical domain, a probable aldehyde dehydrogenase 22A1, and a probable poly(A) polymerase 3 (PAPS3), were identified in these two genomic regions. Amplicon targeted next generation sequencing of these nine candidate genes for Rf1 was performed in an association panel consisting of 27 maintainer and 32 restorer lines and revealed the presence of 210 Single Nucleotide Polymorphisms (SNPs) and 67 Insertions/Deletions (INDELs). Association studies showed significant associations of 10 SNPs with fertility restoration (p-value < 10⁻⁴), narrowing Rf1 down to three candidate genes. Three new markers, one co-dominant marker 67N04_P and two dominant markers, PPR621.5R for restorer, and PPR621.5M for maintainer lines were developed and verified in the association panel of 59 sunflower lines. The versatility of the three newly developed markers, as well as of three existing markers for the restorer gene Rf1 (HRG01 and HRG02, Cleaved Amplified Polymorphic Sequence (CAPS)-marker H13), was analyzed in a large association panel consisting of 557 accessions.

Sunflower pan-genome analysis shows that hybridization altered gene content and disease resistance

Domesticated plants and animals often display dramatic responses to selection, but the origins of the genetic diversity underlying these responses remain poorly understood. Despite domestication and improvement bottlenecks, the cultivated sunflower remains highly variable genetically, possibly due to hybridization with wild relatives. To characterize genetic diversity in the sunflower and to quantify contributions from wild relatives, we sequenced 287 cultivated lines, 17 Native American landraces and 189 wild accessions representing 11 compatible wild species. Cultivar sequences failing to map to the sunflower reference were assembled de novo for each genotype to determine the gene repertoire, or 'pan-genome', of the cultivated sunflower. Assembled genes were then compared to the wild species to estimate origins. Results indicate that the cultivated sunflower pan-genome comprises 61,205 genes, of which 27% vary across genotypes. Approximately 10% of the cultivated sunflower pan-genome is derived through introgression from wild sunflower species, and 1.5% of genes originated solely through introgression. Gene ontology functional analyses further indicate that genes associated with biotic resistance are over-represented among introgressed regions, an observation consistent with breeding records. Analyses of allelic variation associated with downy mildew resistance provide an example in which such introgressions have contributed to resistance to a globally challenging disease.

Genomic sequence and copy number evolution during hybrid crop development in sunflowers

Hybrid crops, an important part of modern agriculture, rely on the development of male and female heterotic gene pools. In sunflowers, heterotic gene pools were developed through the use of crop-wild relatives to produce cytoplasmic male sterile female and branching, fertility restoring male lines. Here, we use genomic data from a diversity panel of male, female, and open-pollinated lines to explore the genetic changes brought during modern improvement. We find the male lines have diverged most from their open-pollinated progenitors and that genetic differentiation is concentrated in chromosomes, 8, 10 and 13, due to introgressions from wild relatives. Ancestral variation from open-pollinated varieties almost universally evolved in parallel for both male and female lines suggesting little or no selection for heterotic overdominance. Furthermore, we show that gene content differs between the male and female lines and that differentiation in gene content is concentrated in high FST regions. This means that the introgressions that brought branching and fertility restoration to the male lines, brought with them different gene content from the ancestral haplotypes, including the removal of some genes. Although we find no evidence that gene complementation genomewide is responsible for heterosis between male and female lines, several of the genes that are largely absent in either the male or female lines are associated with pathogen defense, suggesting complementation may be functionally relevant for crop breeders.

Genomic sequence and copy number evolution during hybrid crop development in sunflowers

Hybrid crops, an important part of modern agriculture, rely on the development of male and female heterotic gene pools. In sunflowers, heterotic gene pools were developed through the use of crop-wild relatives to produce cytoplasmic male sterile female and branching, fertility restoring male lines. Here, we use genomic data from a diversity panel of male, female, and open-pollinated lines to explore the genetic changes brought during modern improvement. We find the male lines have diverged most from their open-pollinated progenitors and that genetic differentiation is concentrated in chromosomes, 8, 10 and 13, due to introgressions from wild relatives. Ancestral variation from open-pollinated varieties almost universally evolved in parallel for both male and female lines suggesting little or no selection for heterotic overdominance. Furthermore, we show that gene content differs between the male and female lines and that differentiation in gene content is concentrated in high F_ST regions. This means that the introgressions that brought branching and fertility restoration to the male lines, brought with them different gene content from the ancestral haplotypes, including the removal of some genes. Although we find no evidence that gene complementation genomewide is responsible for heterosis between male and female lines, several of the genes that are largely absent in either the male or female lines are associated with pathogen defense, suggesting complementation may be functionally relevant for crop breeders.

Sunflower Genetics from Ancestors to Modern Hybrids-A Review

Domestication and the first steps of sunflower breeding date back more than 4000 years. As an interesting crop to humans, sunflower underwent significant changes in the past to finally find its place as one of the most significant oil crops today. Substantial progress has already been made in understanding how sunflower was domesticated. Recent advances in molecular techniques with improved experimental designs contributed to further understanding of the genetic and molecular basis underlying the architectural and phenotypic changes that occurred during domestication and improvements in sunflower breeding. Understanding the domestication process and assessing the current situation concerning available genotypic variations are essential in order for breeders to face future challenges. A review of the tools that are used for exploring the genetic and genome changes associated with sunflower domestication is given in the paper, along with a discussion of their possible implications on classical sunflower breeding techniques and goals.

Multiple genomic regions influence root morphology and seedling growth in cultivated sunflower (Helianthus annuus L.) under well-watered and water-limited conditions

With climate change and an ever-increasing human population threatening food security, developing a better understanding of the genetic basis of crop performance under stressful conditions has become increasingly important. Here, we used genome-wide association studies to genetically dissect variation in seedling growth traits in cultivated sunflower (Helianthus annuus L.) under well-watered and water-limited (i.e., osmotic stress) conditions, with a particular focus on root morphology. Water limitation reduced seedling size and produced a shift toward deeper rooting. These effects varied across genotypes, and we identified 13 genomic regions that were associated with traits of interest across the two environments. These regions varied in size from a single marker to 186.2 Mbp and harbored numerous genes, some of which are known to be involved in the plant growth/development as well as the response to osmotic stress. In many cases, these associations corresponded to growth traits where the common allele outperformed the rare variant, suggesting that selection for increased vigor during the evolution of cultivated sunflower might be responsible for the relatively high frequency of these alleles. We also found evidence of pleiotropy across multiple traits, as well as numerous environmentally independent genetic effects. Overall, our results indicate the existence of genetic variation in root morphology and allocation and further suggest that the majority of alleles associated with these traits have consistent effects across environments.

Specific LTR-Retrotransposons Show Copy Number Variations between Wild and Cultivated Sunflowers

The relationship between variation of the repetitive component of the genome and domestication in plant species is not fully understood. In previous work, variations in the abundance and proximity to genes of long terminal repeats (LTR)-retrotransposons of sunflower (Helianthus annuus L.) were investigated by Illumina DNA sequencingtocompare cultivars and wild accessions. In this study, we annotated and characterized 22 specific retrotransposon families whose abundance varies between domesticated and wild genotypes. These families mostly belonged to the Chromovirus lineage of the Gypsy superfamily and were distributed overall chromosomes. They were also analyzed in respect to their proximity to genes. Genes close to retrotransposon were classified according to biochemical pathways, and differences between domesticated and wild genotypes are shown. These data suggest that structural variations related to retrotransposons might have occurred to produce phenotypic variation between wild and domesticated genotypes, possibly by affecting the expression of genes that lie close to inserted or deleted retrotransposons and belong to specific biochemical pathways as those involved in plant stress responses.

Sunflower Hybrid Breeding: From Markers to Genomic Selection

In sunflower, molecular markers for simple traits as, e.g., fertility restoration, high oleic acid content, herbicide tolerance or resistances to Plasmopara halstedii, Puccinia helianthi, or Orobanche cumana have been successfully used in marker-assisted breeding programs for years. However, agronomically important complex quantitative traits like yield, heterosis, drought tolerance, oil content or selection for disease resistance, e.g., against Sclerotinia sclerotiorum have been challenging and will require genome-wide approaches. Plant genetic resources for sunflower are being collected and conserved worldwide that represent valuable resources to study complex traits. Sunflower association panels provide the basis for genome-wide association studies, overcoming disadvantages of biparental populations. Advances in technologies and the availability of the sunflower genome sequence made novel approaches on the whole genome level possible. Genotype-by-sequencing, and whole genome sequencing based on next generation sequencing technologies facilitated the production of large amounts of SNP markers for high density maps as well as SNP arrays and allowed genome-wide association studies and genomic selection in sunflower. Genome wide or candidate gene based association studies have been performed for traits like branching, flowering time, resistance to Sclerotinia head and stalk rot. First steps in genomic selection with regard to hybrid performance and hybrid oil content have shown that genomic selection can successfully address complex quantitative traits in sunflower and will help to speed up sunflower breeding programs in the future. To make sunflower more competitive toward other oil crops higher levels of resistance against pathogens and better yield performance are required. In addition, optimizing plant architecture toward a more complex growth type for higher plant densities has the potential to considerably increase yields per hectare. Integrative approaches combining omic technologies (genomics, transcriptomics, proteomics, metabolomics and phenomics) using bioinformatic tools will facilitate the identification of target genes and markers for complex traits and will give a better insight into the mechanisms behind the traits.

Genetic Architecture of Capitate Glandular Trichome Density in Florets of Domesticated Sunflower (Helianthus annuus L.)

Capitate glandular trichomes (CGT), one type of glandular trichomes, are most common in Asteraceae species. CGT can produce various secondary metabolites such as sesquiterpene lactones (STLs) and provide durable resistance to insect pests. In sunflower, CGT-based host resistance is effective to combat the specialist pest, sunflower moth. However, the genetic basis of CGT density is not well understood in sunflower. In this study, we identified two major QTL controlling CGT density in sunflower florets by using a F₄ mapping population derived from the cross HA 300 × RHA 464 with a genetic linkage map constructed from genotyping-by-sequencing data and composed of 2121 SNP markers. One major QTL is located on chromosome 5, which explained 11.61% of the observed phenotypic variation, and the second QTL is located on chromosome 6, which explained 14.06% of the observed phenotypic variation. The QTL effects and the association between CGT density and QTL support interval were confirmed in a validation population which included 39 sunflower inbred lines with diverse genetic backgrounds. We also identified two strong candidate genes in the QTL support intervals, and the functions of their orthologs in other plant species suggested their potential roles in regulating capitate glandular trichome density in sunflower. Our results provide valuable information to sunflower breeding community for developing host resistance to sunflower insect pests.

Main and epistatic QTL analyses for Sclerotinia Head Rot resistance in sunflower

Sclerotinia Head Rot (SHR), a disease caused by Sclerotinia sclerotiorum, is one of the most limiting factors in sunflower production. In this study, we identified genomic loci associated with resistance to SHR to support the development of assisted breeding strategies. We genotyped 114 Recombinant Inbred Lines (RILs) along with their parental lines (PAC2 -partially resistant-and RHA266 -susceptible-) by using a 384 single nucleotide polymorphism (SNP) Illumina Oligo Pool Assay to saturate a sunflower genetic map. Subsequently, we tested these lines for SHR resistance using assisted inoculations with S. sclerotiorum ascospores. We also conducted a randomized complete-block assays with three replicates to visually score disease incidence (DI), disease severity (DS), disease intensity (DInt) and incubation period (IP) through four field trials (2010-2014). We finally assessed main effect quantitative trait loci (M-QTLs) and epistatic QTLs (E-QTLs) by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively. As a result of this study, the improved map incorporates 61 new SNPs over candidate genes. We detected a broad range of narrow sense heritability (h2) values (1.86-59.9%) as well as 36 M-QTLs and 13 E-QTLs along 14 linkage groups (LGs). On LG1, LG10, and LG15, we repeatedly detected QTLs across field trials; which emphasizes their putative effectiveness against SHR. In all selected variables, most of the identified QTLs showed high determination coefficients, associated with moderate to high heritability values. Using markers shared with previous Sclerotinia resistance studies, we compared the QTL locations in LG1, LG2, LG8, LG10, LG11, LG15 and LG16. This study constitutes the largest report of QTLs for SHR resistance in sunflower. Further studies focusing on the regions in LG1, LG10, and LG15 harboring the detected QTLs are necessary to identify causal alleles and contribute to unraveling the complex genetic basis governing the resistance.

Association of functional markers with flowering time in lentil

In the present study, a diverse panel of 96 accessions of lentil germplasm was used to study flowering time over environments and to identify simple sequence repeat markers associated with flowering time through association mapping. The study showed high broad sense heritability estimate (h ² bs=0.93) for flowering time in lentil. Screening of 534 SSR markers resulted in an identification of 75 SSR polymorphic markers (13.9%) across studied genotypes. These markers amplified 266 loci and generated 697 alleles ranging from two to 16 alleles per locus. Model-based cluster analysis used for the determination of population structure resulted in the identification of two distinct subpopulations. Distribution of flowering time was ranged from 40 to 70 days in subpopulation I and from 54 to 69 days in subpopulation II and did not skew either late or early flowering time within a subpopulation. No admixture was observed within the subpopulations. Use of the most accepted maximum likelihood model (P3D mixed linear model with optimum compression) of MTA analysis showed significant association of 26 SSR markers with flowering time at <0.05 probability. The percent of phenotypic explained by each associated marker with flowering time ranged from 2.1 to 21.8% and identified QTLs for flowering time explaining high phenotypic variation across the environments (10.7-21.8%) or in a particular environment (10.2-21.4%). In the present study, 13 EST-SSR showed significant association with flowering time and explained large phenotypic variation (2.3-21.8%) compared to genomic SSR markers (2.1-10.2%). Hence, these markers can be used as functional markers in the lentil breeding program to develop short duration cultivars.

Structural variants in genes associated with human Williams-Beuren syndrome underlie stereotypical hypersociability in domestic dogs

Although considerable progress has been made in understanding the genetic basis of morphologic traits (for example, body size and coat color) in dogs and wolves, the genetic basis of their behavioral divergence is poorly understood. An integrative approach using both behavioral and genetic data is required to understand the molecular underpinnings of the various behavioral characteristics associated with domestication. We analyze a 5-Mb genomic region on chromosome 6 previously found to be under positive selection in domestic dog breeds. Deletion of this region in humans is linked to Williams-Beuren syndrome (WBS), a multisystem congenital disorder characterized by hypersocial behavior. We associate quantitative data on behavioral phenotypes symptomatic of WBS in humans with structural changes in the WBS locus in dogs. We find that hypersociability, a central feature of WBS, is also a core element of domestication that distinguishes dogs from wolves. We provide evidence that structural variants in GTF2I and GTF2IRD1, genes previously implicated in the behavioral phenotype of patients with WBS and contained within the WBS locus, contribute to extreme sociability in dogs. This finding suggests that there are commonalities in the genetic architecture of WBS and canine tameness and that directional selection may have targeted a unique set of linked behavioral genes of large phenotypic effect, allowing for rapid behavioral divergence of dogs and wolves, facilitating coexistence with humans.

Genome-wide divergence, haplotype distribution and population demographic histories for Gossypium hirsutum and Gossypium barbadense as revealed by genome-anchored SNPs

Use of 10,129 singleton SNPs of known genomic location in tetraploid cotton provided unique opportunities to characterize genome-wide diversity among 440 Gossypium hirsutum and 219 G. barbadense cultivars and landrace accessions of widespread origin. Using the SNPs distributed genome-wide, we examined genetic diversity, haplotype distribution and linkage disequilibrium patterns in the G. hirsutum and G. barbadense genomes to clarify population demographic history. Diversity and identity-by-state analyses have revealed little sharing of alleles between the two cultivated allotetraploid genomes, with a few exceptions that indicated sporadic gene flow. We found a high number of new alleles, representing increased nucleotide diversity, on chromosomes 1 and 2 in cultivated G. hirsutum as compared with low nucleotide diversity on these chromosomes in landrace G. hirsutum. In contrast, G. barbadense chromosomes showed negative Tajima's D on several chromosomes for both cultivated and landrace types, which indicate that speciation of G. barbadense itself, might have occurred with relatively narrow genetic diversity. The presence of conserved linkage disequilibrium (LD) blocks and haplotypes between G. hirsutum and G. barbadense provides strong evidence for comparable patterns of evolution in their domestication processes. Our study illustrates the potential use of population genetic techniques to identify genomic regions for domestication.

Range-wide phenotypic and genetic differentiation in wild sunflower

BACKGROUND

Divergent phenotypes and genotypes are key signals for identifying the targets of natural selection in locally adapted populations. Here, we used a combination of common garden phenotyping for a variety of growth, plant architecture, and seed traits, along with single-nucleotide polymorphism (SNP) genotyping to characterize range-wide patterns of diversity in 15 populations of wild sunflower (Helianthus annuus L.) sampled along a latitudinal gradient in central North America. We analyzed geographic patterns of phenotypic diversity, quantified levels of within-population SNP diversity, and also determined the extent of population structure across the range of this species. We then used these data to identify significantly over-differentiated loci as indicators of genomic regions that likely contribute to local adaptation.

RESULTS

Traits including flowering time, plant height, and seed oil composition (i.e., percentage of saturated fatty acids) were significantly correlated with latitude, and thus differentiated northern vs. southern populations. Average pairwise FST was found to be 0.21, and a STRUCTURE analysis identified two significant clusters that largely separated northern and southern individuals. The significant FST outliers included a SNP in HaFT2, a flowering time gene that has been previously shown to co-localize with flowering time QTL, and which exhibits a known cline in gene expression.

CONCLUSIONS

Latitudinal differentiation in both phenotypic traits and SNP allele frequencies is observed across wild sunflower populations in central North America. Such differentiation may play an important adaptive role across the range of this species, and could facilitate adaptation to a changing climate.

Diversity analysis and genomic prediction of Sclerotinia resistance in sunflower using a new 25 K SNP genotyping array

KEY MESSAGE

We have developed a SNP array for sunflower containing more than 25 K markers, representing single loci mostly in or near transcribed regions of the genome. The array was successfully applied to genotype a diversity panel of lines, hybrids, and mapping populations and represented well the genetic diversity of cultivated sunflower. Results of PCoA and population substructure analysis underlined the complexity of the genetic composition of current elite breeding material. The performance of this genotyping platform for genome-based prediction of phenotypes and detection of QTL with improved resolution could be demonstrated based on the re-evaluation of a population segregating for resistance to Sclerotinia midstalk rot. Given our results, the newly developed 25 K SNP array is expected to be of great utility for the most important applications in genome-based sunflower breeding and research.

ABSTRACT

Genotyping with a large number of molecular markers is a prerequisite to conduct genome-based genetic analyses with high precision. Here, we report the design and performance of a 25 K SNP genotyping array for sunflower (Helianthus annuus L.). SNPs were discovered based on variant calling in de novo assembled, UniGene-based contigs of sunflower derived from whole genome sequencing and amplicon sequences originating from four and 48 inbred lines, respectively. After inclusion of publically available transcriptome-derived SNPs, in silico design of the Illumina(®) Infinium iSelect HD BeadChip yielded successful assays for 22,299 predominantly haplotype-specific SNPs. The array was validated in a sunflower diversity panel including inbred lines, open-pollinated varieties, introgression lines, landraces, recombinant inbred lines, and F2 populations. Validation provided 20,502 high-quality bi-allelic SNPs with stable cluster performance whereby each SNP marker represents a single locus mostly in or near transcribed regions of the sunflower genome. Analyses of population structure and quantitative resistance to Sclerotinia midstalk rot demonstrate that this array represents a significant improvement over currently available genomic tools for genetic diversity analyses, genome-wide marker-trait association studies, and genetic mapping in sunflower.

Repetitive DNA and Plant Domestication: Variation in Copy Number and Proximity to Genes of LTR-Retrotransposons among Wild and Cultivated Sunflower (Helianthus annuus) Genotypes

The sunflower (Helianthus annuus) genome contains a very large proportion of transposable elements, especially long terminal repeat retrotransposons. However, knowledge on the retrotransposon-related variability within this species is still limited. We used next-generation sequencing (NGS) technologies to perform a quantitative and qualitative survey of intraspecific variation of the retrotransposon fraction of the genome across 15 genotypes—7 wild accessions and 8 cultivars—of H. annuus. By mapping the Illumina reads of the 15 genotypes onto a library of sunflower long terminal repeat retrotransposons, we observed considerable variability in redundancy among genotypes, at both superfamily and family levels. In another analysis, we mapped Illumina paired reads to two sets of sequences, that is, long terminal repeat retrotransposons and protein-encoding sequences, and evaluated the extent of retrotransposon proximity to genes in the sunflower genome by counting the number of paired reads in which one read mapped to a retrotransposon and the other to a gene. Large variability among genotypes was also ascertained for retrotransposon proximity to genes. Both long terminal repeat retrotransposon redundancy and proximity to genes varied among retrotransposon families and also between cultivated and wild genotypes. Such differences are discussed in relation to the possible role of long terminal repeat retrotransposons in the domestication of sunflower.

Population Genomics for Understanding Adaptation in Wild Plant Species

Darwin's theory of evolution by natural selection is the foundation of modern biology. However, it has proven remarkably difficult to demonstrate at the genetic, genomic, and population level exactly how wild species adapt to their natural environments. We discuss how one can use large sets of multiple genome sequences from wild populations to understand adaptation, with an emphasis on the small herbaceous plant Arabidopsis thaliana. We present motivation for such studies; summarize progress in describing whole-genome, species-wide sequence variation; and then discuss what insights have emerged from these resources, either based on sequence information alone or in combination with phenotypic data. We conclude with thoughts on opportunities with other plant species and the impact of expected progress in sequencing technology and genome engineering for studying adaptation in nature.

RNA-Seq Uncovers SNPs and Alternative Splicing Events in Asian Lotus (Nelumbo nucifera)

RNA-Seq is an efficient way to comprehensively identify single nucleotide polymorphisms (SNPs) and alternative splicing (AS) events from the expressed genes. In this study, we conducted transcriptome sequencing of four Asian lotus (Nelumbo nucifera) cultivars using Illumina HiSeq2000 platform to identify SNPs and AS events in lotus. A total of 505 million pair-end RNA-Seq reads were generated from four cultivars, of which 86% were mapped to the lotus reference genome. Using the four sets of data together, a total of 357,689 putative SNPs were identified with an average density of one SNP per 2.2 kb. These SNPs were located in 1,253 scaffolds and 15,016 expressed genes. A/G and C/T were the two major types of SNPs in the Asian lotus transcriptome. In parallel, a total of 177,540 AS events were detected in the four cultivars and were distributed in 64% of the expressed genes of lotus. The predominant type of AS events was alternative 5’ first exon, which accounted for 41.2% of all the observed AS events, and exon skipping only accounted for 4.3% of all AS. Gene Ontology analysis was conducted to analyze the function of the genes containing SNPs and AS events. Validation of selected SNPs and AS events revealed that 74% of SNPs and 80% of AS events were reliable, which indicates that RNA-Seq is an efficient approach to uncover gene-associated SNPs and AS events. A large number of SNPs and AS events identified in our study will facilitate further genetic and functional genomics research in lotus.

Sesquiterpene Lactone Composition of Wild and Cultivated Sunflowers and Biological Activity against an Insect Pest

Sesquiterpene lactones in sunflowers, Helianthus spp., are important to interactions with pathogens, weeds, and insects. Across a broad range of Helianthus annuus, differences in composition of sesquiterpene lactones extracted from disc florets were found between wild and cultivated sunflowers and also between distinct groups of inbreds used to produce sunflower hybrids. Discriminant function analysis showed the presence and relative abundance of argophyllone B, niveusin B, and 15-hydroxy-3-dehydrodesoxyfruticin were usually (75%) effective at classifying wild sunflowers, cultivated inbreds, and hybrids. Argophyllone B reduced the larval mass of the sunflower moth, Homeosoma electellum, by >30%, but only at a dose greater than that found in florets. Low doses of mixed extracts from cultivated florets produced a similar (≈40%) reduction in larval mass, suggesting combinations of sesquiterpene lactones act additively. Although the results support a role for sesquiterpene lactones in herbivore defense of cultivated sunflowers, additional information is needed to use these compounds purposefully in breeding.

Genome scans reveal candidate domestication and improvement genes in cultivated sunflower, as well as post-domestication introgression with wild relatives

The development of modern crops typically involves both selection and hybridization, but to date most studies have focused on the former. In the present study, we explore how both processes, and their interactions, have molded the genome of the cultivated sunflower (Helianthus annuus), a globally important oilseed. To identify genes targeted by selection during the domestication and improvement of sunflower, and to detect post-domestication hybridization with wild species, we analyzed transcriptome sequences of 80 genotypes, including wild, landrace, and modern lines of H. annuus, as well as two cross-compatible wild relatives, Helianthus argophyllus and Helianthus petiolaris. Outlier analyses identified 122 and 15 candidate genes associated with domestication and improvement, respectively. As in several previous studies, genes putatively involved in oil biosynthesis were the most extreme outliers. Additionally, several promising associations were observed with previously mapped quantitative trait loci (QTLs), such as branching. Admixture analyses revealed that all the modern cultivar genomes we examined contained one or more introgressions from wild populations, with every chromosome having evidence of introgression in at least one modern line. Cumulatively, introgressions cover c. 10% of the cultivated sunflower genome. Surprisingly, introgressions do not avoid candidate domestication genes, probably because of the reintroduction of branching.