The Investigation of Perennial Sunflower Species (Helianthus L.) Mitochondrial Genomes.

The genus Helianthus is a diverse taxonomic group with approximately 50 species. Most sunflower genomic investigations are devoted to economically valuable species, e.g., H. annuus, while other Helianthus species, especially perennial, are predominantly a blind spot. In the current study, we have assembled the complete mitogenomes of two perennial species: H. grosseserratus (273,543 bp) and H. strumosus (281,055 bp). We analyzed their sequences and gene profiles in comparison to the available complete mitogenomes of H. annuus. Except for sdh4 and trnA-UGC, both perennial sunflower species had the same gene content and almost identical protein-coding sequences when compared with each other and with annual sunflowers (H. annuus). Common mitochondrial open reading frames (ORFs) (orf117, orf139, and orf334) in sunflowers and unique ORFs for H. grosseserratus (orf633) and H. strumosus (orf126, orf184, orf207) were identified. The maintenance of plastid-derived coding sequences in the mitogenomes of both annual and perennial sunflowers and the low frequency of nonsynonymous mutations point at an extremely low variability of mitochondrial DNA (mtDNA) coding sequences in the Helianthus genus.

Jerusalem artichoke (Helianthus tuberosus L.) as a medicinal plant and its natural products.

Identifying the nutritional and health properties of Helianthus tuberosus, and learning more about this valuable species. It is believed that increased consumption of Jerusalem artichoke (JA) products is related to low blood pressure. One of many questions to answer is whether supplementation of inulin and inulin derivatives obtained from Helianthus tuberosus tubers and aerial parts can be used as antidiabetic, anti-carcinogenic, anti-fungistatic, anti-constipation, body mass-reducing, metabolism-improving agents. We ran a search in Medline, Web of Science, Scopus, Agricola, EBSCO - Food Science Source, Europe PMC, PBL, Cochrane Central Register of Controlled Trials until March 2020. We also browsed reference lists of articles and previous reviews. No language limitations were applied. Jerusalem artichoke (JA) has multiple applications thanks to its rich chemical composition, resistance to biotic and abiotic factors, as: functional food, bioactive ingredient, raw material for the production of ethanol and butanol, succinic, citric and lactic acid. It can be used in medicine and the pharmaceutical industry, because it contains anti-fungistatic, anti-carcinogenic and antioxidant components, and the production of the raw material is easy and inexpensive. It also lowers high cholesterol, triglycerides and high glucose levels; facilitates weight loss; detoxes the organism (e.g. alcohol, heavy metals, radionuclides); lowers uric acid levels; has immunostimulating properties; protects the gastric mucosa, prevents constipation; prevents acne; improves metabolism in lipid disorders; reduces body mass; has cytotoxic properties in breast cancer. It also helps in cardiovascular diseases, chronic infectious diseases; chronic fatigue syndrome; gut flora disorders; immune system disorders. A number of Jerusalem artichoke-derived products were discussed.

Genetic Dissection of Phomopsis Stem Canker Resistance in Cultivated Sunflower Using High Density SNP Linkage Map.

Phomopsis stem canker (PSC) caused by Diaporthe helianthi is increasingly becoming a global threat for sunflower production. In this study, the genetic basis of PSC resistance was investigated in a recombinant inbred line (RIL) population developed from a cross between HA 89 (susceptible) and HA-R3 (resistant). The RIL population was evaluated for PSC disease incidence (DI) in seven screening trials at multiple locations during 2016-2018. The distribution of PSC DI in the RIL population was continuous, confirming a polygenic inheritance of the trait. A moderately high broad-sense heritability (H2, 0.76) was estimated for the trait across environments. In the combined analysis, both the genotype and the genotype × environment interactions were highly significant. A linkage map spanning 1505.33 cM was constructed using genotyping-by-sequencing derived markers. Marker-trait association analysis identified a total of 15 quantitative trait loci (QTL) associated with PSC resistance on 11 sunflower chromosomes, each explaining between 5.24 and 17.39% of the phenotypic variation. PSC resistance QTL were detected in two genomic regions each on chromosomes 3, 5, 13, and 17, while one QTL each was detected in the remaining seven chromosomes. Tightly linked single nucleotide polymorphism (SNP) markers flanking the PSC resistance QTL will facilitate marker-assisted selection in PSC resistance sunflower breeding.

Skim-Sequencing Reveals the Likely Origin of the Enigmatic Endangered Sunflower Helianthus schweinitzii.

Resolving the origin of endangered taxa is an essential component of conservation. This information can be used to guide efforts of bolstering genetic diversity, and also enables species recovery and future evolutionary studies. Here, we used low-coverage whole genome sequencing to clarify the origin of Helianthus schweinitzii, an endangered tetraploid sunflower that is endemic to the Piedmont Plateau in the eastern United States. We surveyed four accessions representing four populations of H. schweinitzii and 38 accessions of six purported parental species. Using de novo approaches, we assembled 87,004 bp of the chloroplast genome and 6770 bp of the nuclear 35S rDNA. Phylogenetic reconstructions based on the chloroplast genome revealed no reciprocal monophyly of taxa. In contrast, nuclear rDNA data strongly supported the currently accepted sections of the genus Helianthus. Information from combined cpDNA and rDNA provided evidence that H. schweinitzii is likely an allo-tetraploid that formed as a result of hybridization between the diploids Helianthus giganteus and Helianthus microcephalus.

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.

Different histories of two highly variable LTR retrotransposons in sunflower species.

In the Helianthus genus, very large intra- and interspecific variability related to two specific retrotransposons of Helianthus annuus (Helicopia and SURE) exists. When comparing these two sequences to sunflower sequence databases recently produced by our lab, the Helicopia family was shown to belong to the Maximus/SIRE lineage of the Sirevirus genus of the Copia superfamily, whereas the SURE element (whose superfamily was not even previously identified) was classified as a Gypsy element of the Ogre/Tat lineage of the Metavirus genus. Bioinformatic analysis of the two retrotransposon families revealed their genomic abundance and relative proliferation timing. The genomic abundance of these families differed significantly among 12 Helianthus species. The ratio between the abundance of long terminal repeats and their reverse transcriptases suggested that the SURE family has relatively more solo long terminal repeats than does Helicopia. Pairwise comparisons of Illumina reads encoding the reverse transcriptase domain indicated that SURE amplification may have occurred more recently than that of Helicopia. Finally, the analysis of population structure based on the SURE and Helicopia polymorphisms of 32 Helianthus species evidenced two subpopulations, which roughly corresponded to species of the Helianthus and Divaricati/Ciliares sections. However, a number of species showed an admixed structure, confirming the importance of interspecific hybridisation in the evolution of this genus. In general, these two retrotransposon families differentially contributed to interspecific variability, emphasising the need to refer to specific families when studying genome evolution.

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.

Genotyping-by-sequencing targeting of a novel downy mildew resistance gene Pl 20 from wild Helianthus argophyllus for sunflower (Helianthus annuus L.).

KEY MESSAGE: Genotyping-by-sequencing revealed a new downy mildew resistance gene, Pl 20 , from wild Helianthus argophyllus located on linkage group 8 of the sunflower genome and closely linked to SNP markers that facilitate the marker-assisted selection of resistance genes. Downy mildew (DM), caused by Plasmopara halstedii, is one of the most devastating and yield-limiting diseases of sunflower. Downy mildew resistance identified in wild Helianthus argophyllus accession PI 494578 was determined to be effective against the predominant and virulent races of P. halstedii occurring in the United States. The evaluation of 114 BC1F2:3 families derived from the cross between HA 89 and PI 494578 against P. halstedii race 734 revealed that single dominant gene controls downy mildew resistance in the population. Genotyping-by-sequencing analysis conducted in the BC1F2 population indicated that the DM resistance gene derived from wild H. argophyllus PI 494578 is located on the upper end of the linkage group (LG) 8 of the sunflower genome, as was determined single nucleotide polymorphism (SNP) markers associated with DM resistance. Analysis of 11 additional SNP markers previously mapped to this region revealed that the resistance gene, named Pl 20 , co-segregated with four markers, SFW02745, SFW09076, S8_11272025, and S8_11272046, and is flanked by SFW04358 and S8_100385559 at an interval of 1.8 cM. The newly discovered P. halstedii resistance gene has been introgressed from wild species into cultivated sunflower to provide a novel gene with DM resistance. The homozygous resistant individuals were selected from BC2F2 progenies with the use of markers linked to the Pl 20 gene, and these lines should benefit the sunflower community for Helianthus improvement.

Molecular cloning, phylogenetic analysis, and expression patterns of LATERAL SUPPRESSOR-LIKE and REGULATOR OF AXILLARY MERISTEM FORMATION-LIKE genes in sunflower (Helianthus annuus L.).

The wild sunflower (Helianthus annuus) plants develop a highly branched form with numerous small flowering heads. The origin of a no branched sunflower, producing a single large head, has been a key event in the domestication process of this species. The interaction between hormonal factors and several genes organizes the initiation and outgrowth of axillary meristems (AMs). From sunflower, we have isolated two genes putatively involved in this process, LATERAL SUPPRESSOR (LS)-LIKE (Ha-LSL) and REGULATOR OF AXILLARY MERISTEM FORMATION (ROX)-LIKE (Ha-ROXL), encoding for a GRAS and a bHLH transcription factor (TF), respectively. Typical amino acid residues and phylogenetic analyses suggest that Ha-LSL and Ha-ROXL are the orthologs of the branching regulator LS and ROX/LAX1, involved in the growth habit of both dicot and monocot species. qRT-PCR analyses revealed a high accumulation of Ha-LSL transcripts in roots, vegetative shoots, and inflorescence shoots. By contrast, in internodal stems and young leaves, a lower amount of Ha-LSL transcripts was observed. A comparison of transcription patterns between Ha-LSL and Ha-ROXL revealed some analogies but also remarkable differences; in fact, the gene Ha-ROXL displayed a low expression level in all organs analyzed. In situ hybridization (ISH) analysis showed that Ha-ROXL transcription was strongly restricted to a small domain within the boundary zone separating the shoot apical meristem (SAM) and the leaf primordia and in restricted regions of the inflorescence meristem, beforehand the separation of floral bracts from disc flower primordia. These results suggested that Ha-ROXL may be involved to establish a cell niche for the initiation of AMs as well as flower primordia. The accumulation of Ha-LSL transcripts was not restricted to the boundary zones in vegetative and inflorescence shoots, but the mRNA activity was expanded in other cellular domains of primary shoot apical meristem as well as AMs. In addition, Ha-LSL transcript accumulation was also detected in leaves and floral primordia at early stages of development. These results were corroborated by qRT-PCR analyses that evidenced high levels of Ha-LSL transcripts in very young leaves and disc flowers, suggesting a role of Ha-LSL for the early outgrowth of lateral primordia.

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.

The Evolution of Haploid Chromosome Numbers in the Sunflower Family.

Chromosome number changes during the evolution of angiosperms are likely to have played a major role in speciation. Their study is of utmost importance, especially now, as a probabilistic model is available to study chromosome evolution within a phylogenetic framework. In the present study, likelihood models of chromosome number evolution were fitted to the largest family of flowering plants, the Asteraceae. Specifically, a phylogenetic supertree of this family was used to reconstruct the ancestral chromosome number and infer genomic events. Our approach inferred that the ancestral chromosome number of the family is n = 9. Also, according to the model that best explained our data, the evolution of haploid chromosome numbers in Asteraceae was a very dynamic process, with genome duplications and descending dysploidy being the most frequent genomic events in the evolution of this family. This model inferred more than one hundred whole genome duplication events; however, it did not find evidence for a paleopolyploidization at the base of this family, which has previously been hypothesized on the basis of sequence data from a limited number of species. The obtained results and potential causes of these discrepancies are discussed.

Revisiting a classic case of introgression: hybridization and gene flow in Californian sunflowers.

During invasion, colonizing species can hybridize with native species, potentially swamping out native genomes. However, theory predicts that introgression will often be biased into the invading species. Thus, empirical estimates of gene flow between native and invasive species are important to quantify the actual threat of hybridization with invasive species. One classic example of introgression occurs in California, where Helianthus bolanderi was thought to be a hybrid between the serpentine endemic Helianthus exilis and the congeneric invader Helianthus annuus. We used genotyping by sequencing to look for signals of introgression and population structure. We find that H. bolanderi and H. exilis form one genetic clade, with weak population structure that is associated with geographic location rather than soil composition and likely represent a single species, not two. Additionally, while our results confirmed early molecular analysis and failed to support the hybrid origin of H. bolanderi, we did find evidence for introgression mainly into the invader H. annuus, as predicted by theory.

Evolutionary Divergences in Root Exudate Composition among Ecologically-Contrasting Helianthus Species.

Plant roots exude numerous metabolites into the soil that influence nutrient availability. Although root exudate composition is hypothesized to be under selection in low fertility soils, few studies have tested this hypothesis in a phylogenetic framework. In this study, we examined root exudates of three pairs of Helianthus species chosen as phylogenetically-independent contrasts with respect to native soil nutrient availability. Under controlled environmental conditions, seedlings were grown to the three-leaf-pair stage, then transferred to either high or low nutrient treatments. After five days of nutrient treatments, we used gas chromatography-mass spectrometry for analysis of root exudates, and detected 37 metabolites across species. When compared in the high nutrient treatment, species native to low nutrient soils exhibited overall higher exudation than their sister species native to high nutrient soils in all three species pairs, providing support for repeated evolutionary shifts in response to native soil fertility. Species native to low nutrient soils and those native to high nutrient soils responded similarly to low nutrient treatments with increased exudation of organic acids (fumaric, citric, malic acids) and glucose, potentially as a mechanism to enhance nutrition acquisition. However, species native to low nutrient soils also responded to low nutrient treatments with a larger decrease in exudation of amino acids than species native to high nutrient soils in all three species pairs. This indicates that species native to low nutrient soils have evolved a unique sensitivity to changes in nutrient availability for some, but not all, root exudates. Overall, these repeated evolutionary divergences between species native to low nutrient soils and those native to high nutrient soils provide evidence for the adaptive value of root exudation, and its plasticity, in contrasting soil environments.

Phylogenetic structural equation modelling reveals no need for an 'origin' of the leaf economics spectrum.

The leaf economics spectrum (LES) is a prominent ecophysiological paradigm that describes global variation in leaf physiology across plant ecological strategies using a handful of key traits. Nearly a decade ago, Shipley et al. (2006) used structural equation modelling to explore the causal functional relationships among LES traits that give rise to their strong global covariation. They concluded that an unmeasured trait drives LES covariation, sparking efforts to identify the latent physiological trait underlying the 'origin' of the LES. Here, we use newly developed phylogenetic structural equation modelling approaches to reassess these conclusions using both global LES data as well as data collected across scales in the genus Helianthus. For global LES data, accounting for phylogenetic non-independence indicates that no additional unmeasured traits are required to explain LES covariation. Across datasets in Helianthus, trait relationships are highly variable, indicating that global-scale models may poorly describe LES covariation at non-global scales.

Species tree estimation of diploid Helianthus (Asteraceae) using target enrichment.

PREMISE OF THE STUDY: The sunflower genus Helianthus has long been recognized as economically significant, containing species of agricultural and horticultural importance. Additionally, this genus displays a large range of phenotypic and genetic variation, making Helianthus a useful system for studying evolutionary and ecological processes. Here we present the most robust Helianthus phylogeny to date, laying the foundation for future studies of this genus. METHODS: We used a target enrichment approach across 37 diploid Helianthus species/subspecies with a total of 103 accessions. This technique garnered 170 genes used for both coalescent and concatenation analyses. The resulting phylogeny was additionally used to examine the evolution of life history and growth form across the genus. KEY RESULTS: Coalescent and concatenation approaches were largely congruent, resolving a large annual clade and two large perennial clades. However, several relationships deeper within the phylogeny were more weakly supported and incongruent among analyses including the placement of H. agrestis, H. cusickii, H. gracilentus, H. mollis, and H. occidentalis. CONCLUSIONS: The current phylogeny supports three major clades including a large annual clade, a southeastern perennial clade, and another clade of primarily large-statured perennials. Relationships among taxa are more consistent with early phylogenies of the genus using morphological and crossing data than recent efforts using single genes, which highlight the difficulties of phylogenetic estimation in genera known for reticulate evolution. Additionally, conflict and low support at the base of the perennial clades may suggest a rapid radiation and/or ancient introgression within the genus.

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.

RNA-seq analysis and de novo transcriptome assembly of Jerusalem artichoke (Helianthus tuberosus Linne).

Jerusalem artichoke (Helianthus tuberosus L.) has long been cultivated as a vegetable and as a source of fructans (inulin) for pharmaceutical applications in diabetes and obesity prevention. However, transcriptomic and genomic data for Jerusalem artichoke remain scarce. In this study, Illumina RNA sequencing (RNA-Seq) was performed on samples from Jerusalem artichoke leaves, roots, stems and two different tuber tissues (early and late tuber development). Data were used for de novo assembly and characterization of the transcriptome. In total 206,215,632 paired-end reads were generated. These were assembled into 66,322 loci with 272,548 transcripts. Loci were annotated by querying against the NCBI non-redundant, Phytozome and UniProt databases, and 40,215 loci were homologous to existing database sequences. Gene Ontology terms were assigned to 19,848 loci, 15,434 loci were matched to 25 Clusters of Eukaryotic Orthologous Groups classifications, and 11,844 loci were classified into 142 Kyoto Encyclopedia of Genes and Genomes pathways. The assembled loci also contained 10,778 potential simple sequence repeats. The newly assembled transcriptome was used to identify loci with tissue-specific differential expression patterns. In total, 670 loci exhibited tissue-specific expression, and a subset of these were confirmed using RT-PCR and qRT-PCR. Gene expression related to inulin biosynthesis in tuber tissue was also investigated. Exsiting genetic and genomic data for H. tuberosus are scarce. The sequence resources developed in this study will enable the analysis of thousands of transcripts and will thus accelerate marker-assisted breeding studies and studies of inulin biosynthesis in Jerusalem artichoke.

Shared selective pressure and local genomic landscape lead to repeatable patterns of genomic divergence in sunflowers.

The repeated evolution of traits in organisms facing similar environmental conditions is considered to be fundamental evidence for the role of natural selection in moulding phenotypes. Yet, aside from case studies of parallel evolution and its genetic basis, the repeatability of evolution at the level of the whole genome remains poorly characterized. Here, through the use of transcriptome sequencing, we examined genomic divergence for three pairs of sister species of sunflowers. Two of the pairs (Helianthus petiolaris - H. debilis and H. annuus - H. argophyllus) have diverged along a similar latitudinal gradient and presumably experienced similar selective pressure. In contrast, a third species pair (H. exilis - H. bolanderi) diverged along a longitudinal gradient. Analyses of divergence, as measured in terms of FST, indicated little repeatability across the three pairs of species for individual genetic markers (SNPs), modest repeatability at the level of individual genes and the highest repeatability when large regions of the genome were compared. As expected, higher repeatability was observed for the two species pairs that have diverged along a similar latitudinal gradient, with genes involved in flowering time among the most divergent genes. Genes showing extreme low or high differentiation were more similar than genes showing medium levels of divergence, implying that both purifying and divergent selection contributed to repeatable patterns of divergence. The location of a gene along the chromosome also predicted divergence levels, presumably because of shared heterogeneity in both recombination and mutation rates. In conclusion, repeated genome evolution appeared to result from both similar selective pressures and shared local genomic landscapes.

Sequence validation of candidates for selectively important genes in sunflower.

Analyses aimed at identifying genes that have been targeted by past selection provide a powerful means for investigating the molecular basis of adaptive differentiation. In the case of crop plants, such studies have the potential to not only shed light on important evolutionary processes, but also to identify genes of agronomic interest. In this study, we test for evidence of positive selection at the DNA sequence level in a set of candidate genes previously identified in a genome-wide scan for genotypic evidence of selection during the evolution of cultivated sunflower. In the majority of cases, we were able to confirm the effects of selection in shaping diversity at these loci. Notably, the genes that were found to be under selection via our sequence-based analyses were devoid of variation in the cultivated sunflower gene pool. This result confirms a possible strategy for streamlining the search for adaptively-important loci process by pre-screening the derived population to identify the strongest candidates before sequencing them in the ancestral population.

Recent nonhybrid origin of sunflower ecotypes in a novel habitat.

The genomics of local adaptation is an increasingly active field, providing insights into the forces driving ecological speciation and the repeatability of evolution. Demography and gene flow play an important role in determining the paths by which parallel evolution occurs and the genomic signatures of adaptation. In the annual sunflowers, hybridization between species has repeatedly led to the colonization of extreme habitats, such as sand dunes. In a new case of adaptation to sand dunes that occurs in populations of H. petiolaris growing at Great Sand Dunes National Park and Preserve (Colorado), we wished to determine the age and long-term migration patterns of the system, as well as its ancestry. We addressed these questions with restriction-associated DNA (RAD) sequence data, aligned to a reference transcriptome. In an isolation with migration model using RAD sequences, coalescent analysis showed that the dune ecotype originated since the last ice age, which is very recent compared with the hybrid dune species, H. anomalus. Large effective population sizes and substantial numbers of gene migrants per generation between dune and nondune ecotypes explained the highly heterogeneous divergence observed among loci. Analysis of RAD-derived SNPs identified heterogeneous divergence between the dune and nondune ecotypes, as well as identifying its nearest relative. Our results did not support the hypothesis that the dune ecotype has hybrid ancestry, suggesting that adaptation of sunflowers to dunes has occurred by multiple mechanisms. The ancestry and long-term history of gene flow between incipient sunflower species provides valuable context for our understanding of ecological speciation and parallel adaptation.