Molecular demographic history of the annual sunflowers Helianthus annuus and H. petiolaris--large effective population sizes and rates of long-term gene flow

Repeatability of adaptation in sunflowers reveals that genomic regions harbouring inversions also drive adaptation in species lacking an inversion

Local adaptation commonly involves alleles of large effect, which experience fitness advantages when in positive linkage disequilibrium (LD). Because segregating inversions suppress recombination and facilitate the maintenance of LD between locally adapted loci, they are also commonly found to be associated with adaptive divergence. However, it is unclear what fraction of an adaptive response can be attributed to inversions and alleles of large effect, and whether the loci within an inversion could still drive adaptation in the absence of its recombination-suppressing effect. Here, we use genome-wide association studies to explore patterns of local adaptation in three species of sunflower: Helianthus annuus, Helianthus argophyllus, and Helianthus petiolaris, which each harbour a large number of species-specific inversions. We find evidence of significant genome-wide repeatability in signatures of association to phenotypes and environments, which are particularly enriched within regions of the genome harbouring an inversion in one species. This shows that while inversions may facilitate local adaptation, at least some of the loci can still harbour mutations that make substantial contributions without the benefit of recombination suppression in species lacking a segregating inversion. While a large number of genomic regions show evidence of repeated adaptation, most of the strongest signatures of association still tend to be species-specific, indicating substantial genotypic redundancy for local adaptation in these species.

Mitogenomic Research of Silverleaf Sunflower (Helianthus argophyllus) and Its Interspecific Hybrids

Interspecific hybridization is widespread for sunflowers, both in wild populations and commercial breeding. One of the most common species that can efficiently cross with Helianthus annuus is the Silverleaf sunflower-Helianthus argophyllus. The current study carried out structural and functional organization analyses of mitochondrial DNA in H. argophyllus and the interspecific hybrid, H. annuus (VIR114A line) × H. argophyllus. The complete mitogenome of H. argophyllus counts 300,843 bp, has a similar organization to the mitogenome of cultivated sunflower, and holds SNPs typical for wild sunflowers. RNA editing analysis predicted 484 sites in H. argophyllus mitochondrial CDS. The mitochondrial genome of the H. annuus × H. argophyllus hybrid is identical to the maternal line (VIR114A). We expected that significant rearrangements in the mitochondrial DNA of the hybrid would occur, due to the frequent recombination. However, the hybrid mitogenome lacks rearrangements, presumably due to the preservation of nuclear-cytoplasmic interaction paths.

Re-evaluating Homoploid Reticulate Evolution in Helianthus Sunflowers

Sunflowers of the genus Helianthus are models for hybridization research and contain three of the best-studied examples of homoploid hybrid speciation. To understand a broader picture of hybridization within the annual sunflowers, we used whole-genome resequencing to conduct a phylogenomic analysis and test for gene flow between lineages. We find that all annual sunflower species tested have evidence of admixture, suggesting hybridization was common during the radiation of the genus. Support for the major species tree decreases with increasing recombination rate, consistent with hybridization and introgression contributing to discordant topologies. Admixture graphs found hybridization to be associated with the origins of the three putative hybrid species (Helianthus anomalus, Helianthus deserticola, and Helianthus paradoxus). However, the hybridization events are more ancient than suggested by previous work. Furthermore, H. anomalus and H. deserticola appear to have arisen from a single hybridization event involving an unexpected donor, rather than through multiple independent events as previously proposed. This means our results are consistent with, but not definitive proof of, two ancient independent homoploid hybrid speciation events in the genus. Using a broader data set that covers the whole Helianthus genus, including perennial species, we find that signals of introgression span the genus and beyond, suggesting highly divergent introgression and/or the sorting of ancient haplotypes. Thus, Helianthus can be viewed as a syngameon in which largely reproductively isolated species are linked together by occasional or frequent gene flow.

Uncovering the genomic basis of an extraordinary plant invasion

Invasive species are a key driver of the global biodiversity crisis, but the drivers of invasiveness, including the role of pathogens, remain debated. We investigated the genomic basis of invasiveness in Ambrosia artemisiifolia (common ragweed), introduced to Europe in the late 19th century, by resequencing 655 ragweed genomes, including 308 herbarium specimens collected up to 190 years ago. In invasive European populations, we found selection signatures in defense genes and lower prevalence of disease-inducing plant pathogens. Together with temporal changes in population structure associated with introgression from closely related Ambrosia species, escape from specific microbial enemies likely favored the plant's remarkable success as an invasive species.

The genome sequencing and comparative analysis of a wild kiwifruit Actinidia eriantha

The current kiwifruit industry is mainly based on the cultivars derived from the species Actinidia chinensis (Ac) which may bring risks such as canker disease. Introgression of desired traits from wild relatives is an important method for improving kiwifruit cultivars. Actinidia eriantha (Ae) is a particularly important taxon used for hybridization or introgressive breeding of new kiwifruit cultivars because of its valued species-specific traits. Here, we assembled a chromosome-scale high-quality genome of a Ae sample which was directly collected from its wild populations. Our analysis revealed that 41.3% of the genome consists of repetitive elements, comparable to the percentage in Ac and Ae cultivar "White" genomes. The genomic structural variation, including the presence/absence-variation (PAV) of genes, is distinct between Ae and Ac, despite both sharing the same two kiwifruit-specific whole genome duplication (WGD) events. This suggests that a post-WGD divergence mechanism occurred during their evolution. We further investigated genes involved in ascorbic acid biosynthesis and disease-resistance of Ae, and we found introgressive genome could contribute to the complex relationship between Ae and other representative kiwifruit taxa. Collectively, the Ae genome offers valuable genetic resource to accelerate kiwifruit breeding applications.

Highly Replicated Evolution of Parapatric Ecotypes

Parallel evolution of ecotypes occurs when selection independently drives the evolution of similar traits across similar environments. The multiple origins of ecotypes are often inferred based on a phylogeny that clusters populations according to geographic location and not by the environment they occupy. However, the use of phylogenies to infer parallel evolution in closely related populations is problematic because gene flow and incomplete lineage sorting can uncouple the genetic structure at neutral markers from the colonization history of populations. Here, we demonstrate multiple origins within ecotypes of an Australian wildflower, Senecio lautus. We observed strong genetic structure as well as phylogenetic clustering by geography and show that this is unlikely due to gene flow between parapatric ecotypes, which was surprisingly low. We further confirm this analytically by demonstrating that phylogenetic distortion due to gene flow often requires higher levels of migration than those observed in S. lautus. Our results imply that selection can repeatedly create similar phenotypes despite the perceived homogenizing effects of gene flow.

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights

BACKGROUND AND AIMS

Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species.

METHODS

We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples.

KEY RESULTS

We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small.

CONCLUSIONS

We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

Evolutionary history of Hemerocallis in Japan inferred from chloroplast and nuclear phylogenies and levels of interspecific gene flow

The perennial herb genus Hemerocallis (Asphodelaceae) shows four flowering types: diurnal half-day, diurnal one-day, nocturnal half-day, and nocturnal one-day flowering. These flowering types are corresponding to their main pollinators, and probably act as a primary mechanism of reproductive isolation. To examine how the four flowering types diverged, we reconstructed the phylogeny of the Japanese species of Hemerocallis using 1615 loci of nuclear genome-wide SNPs and 2078 bp sequences of four cpDNA regions. We also examined interspecific gene flows among taxa by an Isolation-with-Migration model and a population structure analysis. Our study revealed an inconsistency between chloroplast and nuclear genome phylogenies, which may have resulted from chloroplast capture. Each of the following five clusters is monophyletic and clearly separated on the nuclear genome-wide phylogenetic tree: (I) two nocturnal flowering species with lemon-yellow flowers, H. citrina (half-day flowering) and H. lilioasphodelus (one-day flowering); (II) a diurnal one-day flowering species with yellow-orange flowers, H. middendorffii; (III) a variety of a diurnal half-day flowering species with reddish orange flowers, H. fulva var. disticha; (IV) another variety of a diurnal half-day flowering species with reddish orange flowers, H. fulva var. aurantiaca, and a diurnal one-day flowering species with yellow-orange flowers, H. major; (V) a diurnal half-day flowering species with yellow-orange flowers, H. hakuunensis. The five clusters are consistent with traditional phenotype-based taxonomy (cluster I, cluster II, and clusters III-V correspond to Hemerocallis sect. Hemerocallis, Capitatae, and Fulvae, respectively). These findings could indicate that three flowering types (nocturnal flowering, diurnal one-day flowering, and diurnal half-day flowering) diverged in early evolutionary stages of Hemerocallis and subsequently a change from diurnal half-day flowering to diurnal one-day flowering occurred in a lineage of H. major. While genetic differentiation among the five clusters was well maintained, significant gene flow was detected between most pairs of taxa, suggesting that repeated hybridization played a role in the evolution of those taxa.

Standing variation rather than recent adaptive introgression probably underlies differentiation of the texanus subspecies of Helianthus annuus

The origins of geographic races in wide-ranging species are poorly understood. In Texas, the texanus subspecies of Helianthus annuus has long been thought to have acquired its defining phenotypic traits via introgression from a local congener, H. debilis, but previous tests of this hypothesis were inconclusive. Here, we explore the origins of H. a. texanus using whole genome sequencing data from across the entire range of H. annuus and possible donor species, as well as phenotypic data from a common garden study. We found that although it is morphologically convergent with H. debilis, H. a. texanus has conflicting signals of introgression. Genome wide tests (Patterson's D and TreeMix) only found evidence of introgression from H. argophyllus (sister species to H. annuus and also sympatric), but not H. debilis, with the exception of one individual of 109 analysed. We further scanned the genome for localized signals of introgression using PCAdmix and found minimal but nonzero introgression from H. debilis and significant introgression from H. argophyllus in some populations. Given the paucity of introgression from H. debilis, we argue that the morphological convergence observed in Texas is probably from standing genetic variation. We also found that genomic differentiation in H. a. texanus is mostly driven by large segregating inversions, several of which have signatures of natural selection based on haplotype frequencies.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.

Comparative analysis of chloroplast genomes of seven perennial Helianthus species

Sequencing and a comparative analysis of the complete chloroplast genomes of seven perennial Helianthus species were carried out. The chloroplast genomes have a typical quadripartite structure, including large and small single regions and a pair of inverted repeats. Genome sizes were between 151,152 bp and 151,289 bp. The genome of H. grosseserratus was the smallest, while that of H. microcephalus was the largest. The size variation of the chloroplast genomes is substantially due to the change in the length of simple sequence repeats (SSRs) in non-coding regions. An analysis of these SSRs revealed 35 polymorphic loci (average PIC value > 0.5) that can be used to examine ecological and evolutionary processes in wild Helianthus species. Eight divergence hotspots, including five intergenic regions (petN-psbM, clpP intron, rps3-rpl16, ndhD-ccsA, and ndhF-rpl32) and three gene regions (rbcL, ycf1, and ndhF) were also identified in Helianthus chloroplast genomes. The evolutionary selection pressure analysis revealed a strong purifying selection. Only the rbcL gene experienced efficiency of positive selection at the annual/perennial transitions. The inverted repeat (IR)/single copy (SC) boundaries were identical in all of these (Helianthus) species. In general, the comparison of the genomes revealed low levels of sequence variability (Pi = 0.00051). This indicates that the chloroplast genomes of the studied perennial species of Helianthus, in addition to purifying selection, are closely related and have a recent divergence time.

Implications for evolutionary trends from the pairing frequencies among golden-winged and blue-winged warblers and their hybrids

Extensive range loss for the Golden-winged Warbler (Vermivora chrysoptera) has occurred in areas of intrusion by the Blue-winged Warbler (V. cyanoptera) potentially related to their close genetic relationship. We compiled data on social pairing from nine studies for 2,679 resident Vermivora to assess evolutionary divergence. Hybridization between pure phenotypes occurred with 1.2% of resident males for sympatric populations. Pairing success rates for Golden-winged Warblers was 83% and for Blue-winged Warblers was 77%. Pairing success for the hybrid Brewster's Warbler was significantly lower from both species at 54%, showing sexual selection against hybrids. Backcross frequencies for Golden-winged Warblers at 4.9% were significantly higher than for Blue-winged Warblers at 1.7%. More frequent backcrossing by Golden-winged Warblers, which produces hybrid phenotypes, may contribute to the replacement of Golden-winged by Blue-winged Warblers. Reproductive isolation due to behavioral isolation plus sexual selection against hybrids was 0.960. Our analyses suggest that plumage differences are the main driving force for this strong isolation with reduced hybrid fitness contributing to a lesser degree. The major impact of plumage differences to reproductive isolation is compatible with genomic analyses (Current Biology, 2016, 26, 2313), which showed the largest genetic difference between these phenotypes occurred with plumage genes. These phenotypes have maintained morphological, behavioral, and ecological differences during two centuries of hybridization. Our estimate of reproductive isolation supports recognition of these phenotypes as two species. The decline and extirpation of the Golden-winged Warbler in almost all areas of recent sympatry suggest that continued coexistence of both species will require eco-geographic isolation.

Massive haplotypes underlie ecotypic differentiation in sunflowers

Species often include multiple ecotypes that are adapted to different environments¹. However, it is unclear how ecotypes arise and how their distinctive combinations of adaptive alleles are maintained despite hybridization with non-adapted populations^2-4. Here, by resequencing 1,506 wild sunflowers from 3 species (Helianthus annuus, Helianthus petiolaris and Helianthus argophyllus), we identify 37 large (1-100 Mbp in size), non-recombining haplotype blocks that are associated with numerous ecologically relevant traits, as well as soil and climate characteristics. Limited recombination in these haplotype blocks keeps adaptive alleles together, and these regions differentiate sunflower ecotypes. For example, haplotype blocks control a 77-day difference in flowering between ecotypes of the silverleaf sunflower H. argophyllus (probably through deletion of a homologue of FLOWERING LOCUS T (FT)), and are associated with seed size, flowering time and soil fertility in dune-adapted sunflowers. These haplotypes are highly divergent, frequently associated with structural variants and often appear to represent introgressions from other-possibly now-extinct-congeners. These results highlight a pervasive role of structural variation in ecotypic adaptation.

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.

Character displacement of a learned behaviour and its implications for ecological speciation

Cultural evolution may accelerate population divergence and speciation, though most support for this hypothesis is restricted to scenarios of allopatric speciation driven by random cultural drift. By contrast, the role of cultural evolution in non-allopatric speciation (i.e. speciation with gene flow) has received much less attention. One clade in which cultural evolution may have figured prominently in speciation with gene flow includes the conifer-seed-eating finches in the red crossbill (Loxia curvirostra) complex. Here we focus on Cassia crossbills (Loxia sinesciuris; an ecotype recently split taxonomically from red crossbills) that learn social contact calls from their parents. Previous work found that individuals modify their calls throughout life such that they become increasingly divergent from a closely related, sympatric red crossbill ecotype. This open-ended modification of calls could lead to character displacement if it causes population-level divergence in call structure that, in turn, reduces (maladaptive) heterospecific flocking. Heterospecific flocking is maladaptive because crossbills use public information from flockmates to assess resource quality, and feeding rates are depressed when flockmates differ in their ability to exploit a shared resource (i.e. when flockmates are heterospecifics). We confirm the predictions of character displacement by documenting substantial population-level divergence in Cassia crossbill call structure over just two decades and by using field experiments to demonstrate that Cassia and red crossbills differentially respond to these evolved differences in call structure, reducing heterospecific flock formation. Moreover, because crossbills choose mates from within flocks, a reduction in heterospecific flocking should increase assortative mating and may have been critical for speciation of Cassia crossbills in the face of ongoing gene flow in as few as 5000 years. Our results provide evidence for a largely neglected yet potentially widespread mechanism by which reproductive isolation can evolve between sympatric lineages as a byproduct of adaptive cultural evolution.

Ecological divergence plays an important role in strong but complex reproductive isolation in campions (Silene)

New species arise through the evolution of reproductive barriers between formerly interbreeding lineages. Yet, comprehensive assessments of potential reproductive barriers, which are needed to make inferences on processes driving speciation, are only available for a limited number of systems. In this study, we estimated individual and cumulative strengths of seven prezygotic and six postzygotic reproductive barriers between the recently diverged taxa Silene dioica (L.) Clairv. and S. latifolia Poiret using both published and new data. A combination of multiple partial reproductive barriers resulted in near-complete reproductive isolation between S. dioica and S. latifolia, consistent with earlier estimates of gene flow between the taxa. Extrinsic barriers associated with adaptive ecological divergence were most important, while intrinsic postzygotic barriers had moderate individual strength but contributed only little to total reproductive isolation. These findings are in line with ecological divergence as driver of speciation. We further found extensive variation in extrinsic reproductive isolation, ranging from sites with very strong selection against migrants and hybrids to intermediate sites where substantial hybridization is possible. This situation may allow for, or even promote, heterogeneous genetic divergence.

The Chrysanthemum nankingense Genome Provides Insights into the Evolution and Diversification of Chrysanthemum Flowers and Medicinal Traits

The Asteraceae (Compositae), a large plant family of approximately 24 000-35 000 species, accounts for ∼10% of all angiosperm species and contributes a lot to plant diversity. The most representative members of the Asteraceae are the economically important chrysanthemums (Chrysanthemum L.) that diversified through reticulate evolution. Biodiversity is typically created by multiple evolutionary mechanisms such as whole-genome duplication (WGD) or polyploidization and locally repetitive genome expansion. However, the lack of genomic data from chrysanthemum species has prevented an in-depth analysis of the evolutionary mechanisms involved in their diversification. Here, we used Oxford Nanopore long-read technology to sequence the diploid Chrysanthemum nankingense genome, which represents one of the progenitor genomes of domesticated chrysanthemums. Our analysis revealed that the evolution of the C. nankingense genome was driven by bursts of repetitive element expansion and WGD events including a recent WGD that distinguishes chrysanthemum from sunflower, which diverged from chrysanthemum approximately 38.8 million years ago. Variations of ornamental and medicinal traits in chrysanthemums are linked to the expansion of candidate gene families by duplication events including paralogous gene duplication. Collectively, our study of the assembled reference genome offers new knowledge and resources to dissect the history and pattern of evolution and diversification of chrysanthemum plants, and also to accelerate their breeding and improvement.

Using Genomic Location and Coalescent Simulation to Investigate Gene Tree Discordance in Medicago L

Several well-documented evolutionary processes are known to cause conflict between species-level phylogenies and gene-level phylogenies. Three of the most challenging processes for species tree inference are incomplete lineage sorting, hybridization and gene duplication, which may result in unwarranted comparisons of paralogous genes. Several existing methods have dealt with these processes but none has yet been able to untangle all three at once. Here, we propose a stepwise method by which these processes can be discerned using information on genomic location coupled with coalescent simulations. In the first step, highly discordant genes within genomic blocks (putative paralogs) are identified and excluded from the data set and, in the second step, blocks of linked genes are grouped according to their hybrid history. Existing multispecies coalescent software can then be applied to recover the principal tree(s) that make up the species tree/network without violating the underlying model. The potential of the approach is evaluated on simulated data derived from a species network composed of nine species, of which one is of hybrid origin, and displaying a single-gene duplication that leads to paralogous comparisons. We apply our method to an empirical set of 12 genes from 7 species sampled in the plant genus Medicago that display phylogenetic discordance. We identify the causes of the discordance and demonstrate that the Medicago orbicularis lineage experienced an episode of ancient hybridization. Our results show promise as a new way to explore phylogenetic sequence data that can significantly improve species tree inference in presence of hybridization and undetected paralogy or other causes leading to extremely discordant gene trees. [Coalescent simulation; gene tree; genomic location; hybridization; incomplete lineage sorting; paralogy; phylogenetic incongruence; principal tree; species tree.].

The evolution of recombination rates in finite populations during ecological speciation

Recombination can impede ecological speciation with gene flow by mixing locally adapted genotypes with maladapted migrant genotypes from a divergent population. In such a scenario, suppression of recombination can be selectively favoured. However, in finite populations evolving under the influence of random genetic drift, recombination can also facilitate adaptation by reducing Hill-Robertson interference between loci under selection. In this case, increased recombination rates can be favoured. Although these two major effects on recombination have been studied individually, their joint effect on ecological speciation with gene flow remains unexplored. Using a mathematical model, we investigated the evolution of recombination rates in two finite populations that exchange migrants while adapting to contrasting environments. Our results indicate a two-step dynamic where increased recombination is first favoured (in response to the Hill-Robertson effect), and then disfavoured, as the cost of recombining locally with maladapted migrant genotypes increases over time (the maladaptive gene flow effect). In larger populations, a stronger initial benefit for recombination was observed, whereas high migration rates intensify the long-term cost of recombination. These dynamics may have important implications for our understanding of the conditions that facilitate incipient speciation with gene flow and the evolution of recombination in finite populations.

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.

Lineages of Silene nutans developed rapid, strong, asymmetric postzygotic reproductive isolation in allopatry

Reproductive isolation can rise either as a consequence of genomic divergence in allopatry or as a byproduct of divergent selection in parapatry. To determine whether reproductive isolation in gynodioecious Silene nutans results from allopatric divergence or from ecological adaptation following secondary contact, we investigated the pattern of postzygotic reproductive isolation and hybridization in natural populations using two phylogeographic lineages, western (W1) and eastern (E1). Experimental crosses between the lineages identified strong, asymmetric postzygotic isolation between the W1 and the E1 lineages, independent of geographic overlap. The proportion of ovules fertilized, seeds aborted, and seeds germinated revealed relatively little effect on the fitness of hybrids. In contrast, hybrid mortality was high and asymmetric: while half of the hybrid seedlings with western lineage mothers died, nearly all hybrid seedlings with E1 mothers died. This asymmetric mortality mirrored the proportion of chlorotic seedlings, and is congruent with cytonuclear incompatibility. We found no evidence of hybridization between the lineages in regions of co-occurrence using nuclear and plastid markers. Together, our results are consistent with the hypothesis that strong postzygotic reproductive isolation involving cytonuclear incompatibilities arose in allopatry. We argue that the dynamics of cytonuclear gynodioecy could facilitate the evolution of reproductive isolation.

The emergence of the hyperinvasive vine, Mikania micrantha (Asteraceae), via admixture and founder events inferred from population transcriptomics

Biological invasions that involve well-documented rapid adaptations to new environments provide unequalled opportunities for testing evolutionary hypotheses. Mikania micrantha Kunth (Asteraceae), a perennial herbaceous vine native to tropical Central and South America, successfully invaded tropical Asia in the early 20th century. It is regarded as one of the most aggressive weeds in the world. To elucidate the molecular and evolutionary processes underlying this invasion, we extensively sampled this weed throughout its invaded range in South-East and South Asia and surveyed its genetic structure using variants detected from population transcriptomics. Clustering results suggest that more than one source population contributed to this invasion. Computer simulations using genomewide genetic variation support a scenario of admixture and founder events during invasion. The genes differentially expressed between native and invasive populations were found to be involved in oxidative and high light intensity stress responses, pointing to a possible ecological mechanism of adaptation. Our results provide a foundation for further detailed mechanistic and population studies of this ecologically and economically important invasion. This line of research promises to provide new mitigation strategies for invasive species as well as insights into mechanisms of adaptation.

Escape to Ferality: The Endoferal Origin of Weedy Rice from Crop Rice through De-Domestication

Domestication is the hallmark of evolution and civilization and harnesses biodiversity through selection for specific traits. In regions where domesticated lines are grown near wild relatives, congeneric sources of aggressive weedy genotypes cause major economic losses. Thus, the origins of weedy genotypes where no congeneric species occur raise questions regarding management effectiveness and evolutionary mechanisms responsible for weedy population success. Since eradication in the 1970s, California growers avoided weedy rice through continuous flood culture and zero-tolerance guidelines, preventing the import, presence, and movement of weedy seeds. In 2003, after decades of no reported presence in California, a weedy rice population was confirmed in dry-seeded fields. Our objectives were to identify the origins and establishment of this population and pinpoint possible phenotypes involved. We show that California weedy rice is derived from a different genetic source among a broad range of AA genome Oryzas and is most recently diverged from O. sativa temperate japonica cultivated in California. In contrast, other weedy rice ecotypes in North America (Southern US) originate from weedy genotypes from China near wild Oryza, and are derived through existing crop-wild relative crosses. Analyses of morphological data show that California weedy rice subgroups have phenotypes like medium-grain or gourmet cultivars, but have colored pericarp, seed shattering, and awns like wild relatives, suggesting that reversion to non-domestic or wild-like traits can occur following domestication, despite apparent fixation of domestication alleles. Additionally, these results indicate that preventive methods focused on incoming weed sources through contamination may miss burgeoning weedy genotypes that rapidly adapt, establish, and proliferate. Investigating the common and unique evolutionary mechanisms underlying global weed origins and subsequent interactions with crop relatives sheds light on how weeds evolve and addresses broader questions regarding the stability of selection during domestication and crop improvement.

Demographic history of speciation in a Senecio altitudinal hybrid zone on Mt. Etna

Hybrid zones typically form as a result of species coming into secondary contact, but can also be established in situ as an ecotonal hybrid zone, a situation which has been reported far less frequently. An altitudinal hybrid zone on Mount Etna between two ragwort species (the low elevation Senecio chrysanthemifolius and high elevation S. aethnensis) could potentially represent either of these possibilities. However, a scenario of secondary contact vs. speciation with gene flow has not been explicitly tested. Here, we test these alternatives and demonstrate that the data do not support secondary contact. Furthermore, we report that the previous analyses of speciation history of these species were based on admixed populations, which has led to inflated estimates of ongoing, interspecific gene flow. Our new analyses, based on 'pure' S. aethnensis and S. chrysanthemifolius populations, reveal gene exchange of less than one effective migrant per generation, a level low enough to allow the species to accumulate neutral, genomewide differences. Overall, our results are consistent with a scenario of speciation with gene flow and a divergence time which coincides with the rise of Mt. Etna to altitudes above 2000 m (~150 KY). Further work to quantify the role of adaptation to contrasting environments of high and low altitudes will be needed to support the scenario of recent ecological speciation in this system.

Maintenance of Species Boundaries Despite Ongoing Gene Flow in Ragworts

The role of hybridization between diversifying species has been the focus of a huge amount of recent evolutionary research. While gene flow can prevent speciation or initiate species collapse, it can also generate new hybrid species. Similarly, while adaptive divergence can be wiped out by gene flow, new adaptive variation can be introduced via introgression. The relative frequency of these outcomes, and indeed the frequency of hybridization and introgression in general are largely unknown. One group of closely-related species with several documented cases of hybridization is the Mediterranean ragwort (genus: Senecio) species-complex. Examples of both polyploid and homoploid hybrid speciation are known in the clade, although their evolutionary relationships and the general frequency of introgressive hybridization among them remain unknown. Using a whole genome gene-space dataset comprising eight Senecio species we fully resolve the phylogeny of these species for the first time despite phylogenetic incongruence across the genome. Using a D-statistic approach, we demonstrate previously unknown cases of introgressive hybridization between multiple pairs of taxa across the species tree. This is an important step in establishing these species as a study system for diversification with gene flow, and suggests that introgressive hybridization may be a widespread and important process in plant evolution.

Reconstructing the demographic history of divergence between European river and brook lampreys using approximate Bayesian computations

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey (L. planeri) show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.

Heterogeneous genome divergence, differential introgression, and the origin and structure of hybrid zones

Hybrid zones have been promoted as windows on the evolutionary process and as laboratories for studying divergence and speciation. Patterns of divergence between hybridizing species can now be characterized on a genomewide scale, and recent genome scans have focused on the presence of 'islands' of divergence. Patterns of heterogeneous genomic divergence may reflect differential introgression following secondary contact and provide insights into which genome regions contribute to local adaptation, hybrid unfitness and positive assortative mating. However, heterogeneous genome divergence can also arise in the absence of any gene flow, as a result of variation in selection and recombination across the genome. We suggest that to understand hybrid zone origins and dynamics, it is essential to distinguish between genome regions that are divergent between pure parental populations and regions that show restricted introgression where these populations interact in hybrid zones. The latter, more so than the former, reveal the likely genetic architecture of reproductive isolation. Mosaic hybrid zones, because of their complex structure and multiple contacts, are particularly good subjects for distinguishing primary intergradation from secondary contact. Comparisons among independent hybrid zones or transects that involve the 'same' species pair can also help to distinguish between divergence with gene flow and secondary contact. However, data from replicate hybrid zones or replicate transects do not reveal consistent patterns; in a few cases, patterns of introgression are similar across independent transects, but for many taxa, there is distinct lack of concordance, presumably due to variation in environmental context and/or variation in the genetics of the interacting populations.

The genome sequence of the outbreeding globe artichoke constructed de novo incorporating a phase-aware low-pass sequencing strategy of F1 progeny

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F₁ individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.

Genome Sequence and Architecture of the Tobacco Downy Mildew Pathogen Peronospora tabacina

Peronospora tabacina is an obligate biotrophic oomycete that causes blue mold or downy mildew on tobacco (Nicotiana tabacum). It is an economically important disease occurring frequently in tobacco-growing regions worldwide. We sequenced and characterized the genomes of two P. tabacina isolates and mined them for pathogenicity-related proteins and effector-encoding genes. De novo assembly of the genomes using Illumina reads resulted in 4,016 (63.1 Mb, N50 = 79 kb) and 3,245 (55.3 Mb, N50 = 61 kb) scaffolds for isolates 968-J2 and 968-S26, respectively, with an estimated genome size of 68 Mb. The mitochondrial genome has a similar size (approximately 43 kb) and structure to those of other oomycetes, plus several minor unique features. Repetitive elements, primarily retrotransposons, make up approximately 24% of the nuclear genome. Approximately 18,000 protein-coding gene models were predicted. Mining the secretome revealed approximately 120 candidate RxLR, six CRN (candidate effectors that elicit crinkling and necrosis), and 61 WY domain-containing proteins. Candidate RxLR effectors were shown to be predominantly undergoing diversifying selection, with approximately 57% located in variable gene-sparse regions of the genome. Aligning the P. tabacina genome to Hyaloperonospora arabidopsidis and Phytophthora spp. revealed a high level of synteny. Blocks of synteny show gene inversions and instances of expansion in intergenic regions. Extensive rearrangements of the gene-rich genomic regions do not appear to have occurred during the evolution of these highly variable pathogens. These assemblies provide the basis for studies of virulence in this and other downy mildew pathogens.

Dynamics of asymmetrical hybridization in North American wood ferns: reconciling patterns of inheritance with gametophyte reproductive biology

Hybridization is an important evolutionary force in plants, but the mechanisms underlying it have not been well studied for many groups. In particular, the drivers of non-random patterns of interspecific gene flow (asymmetrical hybridization) remain poorly understood, especially in the seed-free vascular plants. Here, we examine patterns of asymmetrical hybridization in two widespread fern hybrids from eastern North America and study the role of gametophyte ecology in the determination of hybridization bias. We characterized the maternal parentage of > 140 hybrid sporophytes by sequencing a c. 350-bp region of chloroplast DNA (cpDNA). To identify factors contributing to patterns of asymmetrical hybridization, we cultured gametophytes of the parental species and evaluated critical aspects of their reproductive biology. We found that asymmetrical hybridization was prevalent across the populations of both hybrids. Reproductive traits varied across species and suggest that selfing potential, antheridiogen responsiveness, sperm dispersal capacity and gamete size all contribute to the mediation of the direction of hybridization in this group. Our findings suggest that asymmetrical hybridization in ferns is driven by an array of reproductive traits. This study helps to sharpen and define a mechanistic understanding of patterns of hybridization in this group and demonstrates the importance of considering gametophyte biology when studying evolutionary processes in ferns.

The genomic architecture of population divergence between subspecies of the European rabbit

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

Rapid speciation with gene flow following the formation of Mt. Etna

Environmental or geological changes can create new niches that drive ecological species divergence without the immediate cessation of gene flow. However, few such cases have been characterized. On a recently formed volcano, Mt. Etna, Senecio aethnensis and S. chrysanthemifolius inhabit contrasting environments of high and low altitude, respectively. They have very distinct phenotypes, despite hybridizing promiscuously, and thus may represent an important example of ecological speciation “in action,” possibly as a response to the rapid geological changes that Mt. Etna has recently undergone. To elucidate the species’ evolutionary history, and help establish the species as a study system for speciation genomics, we sequenced the transcriptomes of the two Etnean species, and the outgroup, S. vernalis, using Illumina sequencing. Despite the species’ substantial phenotypic divergence, synonymous divergence between the high- and low-altitude species was low (dS = 0.016 ± 0.017 [SD]). A comparison of species divergence models with and without gene flow provided unequivocal support in favor of the former and demonstrated a recent time of species divergence (153,080 ya ± 11,470 [SE]) that coincides with the growth of Mt. Etna to the altitudes that separate the species today. Analysis of dN/dS revealed wide variation in selective constraint between genes, and evidence that highly expressed genes, more “multifunctional” genes, and those with more paralogs were under elevated purifying selection. Taken together, these results are consistent with a model of ecological speciation, potentially as a response to the emergence of a new, high-altitude niche as the volcano grew.

The bright side of transposons in crop evolution

The past decades have revealed an unexpected yet prominent role of so-called 'junk DNA' in the regulation of gene expression, thereby challenging our view of the mechanisms underlying phenotypic evolution. In particular, several mechanisms through which transposable elements (TEs) participate in functional genome diversity have been depicted, bringing to light the 'TEs bright side'. However, the relative contribution of those mechanisms and, more generally, the importance of TE-based polymorphisms on past and present phenotypic variation in crops species remain poorly understood. Here, we review current knowledge on both issues, and discuss how analyses of massively parallel sequencing data combined with statistical methodologies and functional validations will help unravelling the impact of TEs on crop evolution in a near future.

Genomic islands of divergence are not affected by geography of speciation in sunflowers

Genomic studies of speciation often report the presence of highly differentiated genomic regions interspersed within a milieu of weakly diverged loci. The formation of these speciation islands is generally attributed to reduced inter-population gene flow near loci under divergent selection, but few studies have critically evaluated this hypothesis. Here, we report on transcriptome scans among four recently diverged pairs of sunflower (Helianthus) species that vary in the geographical context of speciation. We find that genetic divergence is lower in sympatric and parapatric comparisons, consistent with a role for gene flow in eroding neutral differences. However, genomic islands of divergence are numerous and small in all comparisons, and contrary to expectations, island number and size are not significantly affected by levels of interspecific gene flow. Rather, island formation is strongly associated with reduced recombination rates. Overall, our results indicate that the functional architecture of genomes plays a larger role in shaping genomic divergence than does the geography of speciation.

Parallel evolution of local adaptation and reproductive isolation in the face of gene flow

Parallel evolution of similar phenotypes provides strong evidence for the operation of natural selection. Where these phenotypes contribute to reproductive isolation, they further support a role for divergent, habitat-associated selection in speciation. However, the observation of pairs of divergent ecotypes currently occupying contrasting habitats in distinct geographical regions is not sufficient to infer parallel origins. Here we show striking parallel phenotypic divergence between populations of the rocky-shore gastropod, Littorina saxatilis, occupying contrasting habitats exposed to either wave action or crab predation. This divergence is associated with barriers to gene exchange but, nevertheless, genetic variation is more strongly structured by geography than by ecotype. Using approximate Bayesian analysis of sequence data and amplified fragment length polymorphism markers, we show that the ecotypes are likely to have arisen in the face of continuous gene flow and that the demographic separation of ecotypes has occurred in parallel at both regional and local scales. Parameter estimates suggest a long delay between colonization of a locality and ecotype formation, perhaps because the postglacial spread of crab populations was slower than the spread of snails. Adaptive differentiation may not be fully genetically independent despite being demographically parallel. These results provide new insight into a major model of ecologically driven speciation.

The timetable for allopolyploidy in flowering plants

BACKGROUND

Our understanding of the processes and dynamics of allopolyploid speciation, the long-term consequences of ploidal change, and the genetic and chromosomal changes in new emerged allopolyploids has substantially increased during the past few decades. Yet we remain uncertain about the time since lineage divergence when two taxa are capable of spawning such entities. Indeed, the matter has seemed intractable. Knowledge of the window of opportunity for allopolyploid production is very important because it provides temporal insight into a key evolutionary process, and a temporal reference against which other modes of speciation may be measured.

SCOPE

This Viewpoint paper reviews and integrates published information on the crossability of herbaceous species and the fertility of their hybrids in relation to species' divergence times. Despite limitations in methodology and sampling, the estimated times to hybrid sterility are somewhat congruent across disparate lineages. Whereas the waiting time for hybrid sterility is roughly 4-5 million years, the waiting time for cross-incompatibility is roughly 8-10 million years, sometimes considerably more. Strict allopolyploids may be formed in the intervening time window. The progenitors of several allopolyploids diverged between 4 and 6 million years before allopolyploid synthesis, as expected. This is the first study to propose a general temporal framework for strict allopolyploidy. This Viewpoint paper hopefully will stimulate interest in studying the tempo of speciation and the tempo of reproductive isolation in general.

Genomic divergence during speciation driven by adaptation to altitude

Even though Darwin's "On the Origin of Species" implied selection being the main driver of species formation, the role of natural selection in speciation remains poorly understood. In particular, it remains unclear how selection at a few genes can lead to genomewide divergence and the formation of distinct species. We used a particularly attractive clear-cut case of recent plant ecological speciation to investigate the demography and genomic bases of species formation driven by adaptation to contrasting conditions. High-altitude Senecio aethnensis and low-altitude S. chrysanthemifolius live at the extremes of a mountain slope on Mt. Etna, Sicily, and form a hybrid zone at intermediate altitudes but remain morphologically distinct. Genetic differentiation of these species was analyzed at the DNA polymorphism and gene expression levels by high-throughput sequencing of transcriptomes from multiple individuals. Out of ≈ 18,000 genes analyzed, only a small number (90) displayed differential expression between the two species. These genes showed significantly elevated species differentiation (FST and Dxy), consistent with diversifying selection acting on these genes. Genomewide genetic differentiation of the species is surprisingly low (FST = 0.19), while ≈ 200 genes showed significantly higher (false discovery rate < 1%; mean outlier FST > 0.6) interspecific differentiation and evidence for local adaptation. Diversifying selection at only a handful of loci may be enough for the formation and maintenance of taxonomically well-defined species, despite ongoing gene flow. This provides an explanation of why many closely related species (in plants, in particular) remain phenotypically and ecologically distinct despite ongoing hybridization, a question that has long puzzled naturalists and geneticists alike.

Signatures of demography and recombination at coding genes in naturally-distributed populations of Arabidopsis lyrata subsp. petraea

Demography impacts the observed standing level of genetic diversity present in populations. Distinguishing the relative impacts of demography from selection requires a baseline of expressed gene variation in naturally occurring populations. Six nuclear genes were sequenced to estimate the patterns and levels of genetic diversity in natural Arabidopsis lyrata subsp. petraea populations that differ in demographic histories since the Pleistocene. As expected, northern European populations have genetic signatures of a strong population bottleneck likely due to glaciation during the Pleistocene. Levels of diversity in the northern populations are about half of that in central European populations. Bayesian estimates of historical population size changes indicate that central European populations also have signatures of population size change since the last glacial maxima, suggesting that these populations are not as stable as previously thought. Time since divergence amongst northern European populations is higher than amongst central European populations, suggesting that the northern European populations were established before the Pleistocene and survived glaciation in small separated refugia. Estimates of demography based on expressed genes are complementary to estimates based on microsatellites and transposable elements, elucidating temporal shifts in population dynamics and confirming the importance of marker selection for tests of demography.

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation-migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π(sil)=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F(st)=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.

The Population Genomics of Sunflowers and Genomic Determinants of Protein Evolution Revealed by RNAseq

Few studies have investigated the causes of evolutionary rate variation among plant nuclear genes, especially in recently diverged species still capable of hybridizing in the wild. The recent advent of Next Generation Sequencing (NGS) permits investigation of genome wide rates of protein evolution and the role of selection in generating and maintaining divergence. Here, we use individual whole-transcriptome sequencing (RNAseq) to refine our understanding of the population genomics of wild species of sunflowers (Helianthus spp.) and the factors that affect rates of protein evolution. We aligned 35 GB of transcriptome sequencing data and identified 433,257 polymorphic sites (SNPs) in a reference transcriptome comprising 16,312 genes. Using SNP markers, we identified strong population clustering largely corresponding to the three species analyzed here (Helianthus annuus, H. petiolaris, H. debilis), with one distinct early generation hybrid. Then, we calculated the proportions of adaptive substitution fixed by selection (alpha) and identified gene ontology categories with elevated values of alpha. The “response to biotic stimulus” category had the highest mean alpha across the three interspecific comparisons, implying that natural selection imposed by other organisms plays an important role in driving protein evolution in wild sunflowers. Finally, we examined the relationship between protein evolution (d_N/d_S ratio) and several genomic factors predicted to co-vary with protein evolution (gene expression level, divergence and specificity, genetic divergence [F_ST], and nucleotide diversity pi). We find that variation in rates of protein divergence was correlated with gene expression level and specificity, consistent with results from a broad range of taxa and timescales. This would in turn imply that these factors govern protein evolution both at a microevolutionary and macroevolutionary timescale. Our results contribute to a general understanding of the determinants of rates of protein evolution and the impact of selection on patterns of polymorphism and divergence.

Sunflower genetic, genomic and ecological resources

Long a major focus of genetic research and breeding, sunflowers (Helianthus) are emerging as an increasingly important experimental system for ecological and evolutionary studies. Here, we review the various attributes of wild and domesticated sunflowers that make them valuable for ecological experimentation and describe the numerous publicly available resources that have enabled rapid advances in ecological and evolutionary genetics. Resources include seed collections available from germplasm centres at the USDA and INRA, genomic and EST sequences, mapping populations, genetic markers, genetic and physical maps and other forward- and reverse-genetic tools. We also discuss some of the key evolutionary, genetic and ecological questions being addressed in sunflowers, as well as gaps in our knowledge and promising areas for future research.

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome-scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well-defined categories of locus types (e.g., neutral vs. highly selected). Our multi-locus analysis reveals cross-genomic variation in selective constraints on gene flow and locus-specific flexibility in the permeability of the interspecies membrane.

The sunflower (Helianthus annuus L.) genome reflects a recent history of biased accumulation of transposable elements

Aside from polyploidy, transposable elements are the major drivers of genome size increases in plants. Thus, understanding the diversity and evolutionary dynamics of transposable elements in sunflower (Helianthus annuus L.), especially given its large genome size (∼3.5 Gb) and the well-documented cases of amplification of certain transposons within the genus, is of considerable importance for understanding the evolutionary history of this emerging model species. By analyzing approximately 25% of the sunflower genome from random sequence reads and assembled bacterial artificial chromosome (BAC) clones, we show that it is composed of over 81% transposable elements, 77% of which are long terminal repeat (LTR) retrotransposons. Moreover, the LTR retrotransposon fraction in BAC clones harboring genes is disproportionately composed of chromodomain-containing Gypsy LTR retrotransposons ('chromoviruses'), and the majority of the intact chromoviruses contain tandem chromodomain duplications. We show that there is a bias in the efficacy of homologous recombination in removing LTR retrotransposon DNA, thereby providing insight into the mechanisms associated with transposable element (TE) composition in the sunflower genome. We also show that the vast majority of observed LTR retrotransposon insertions have likely occurred since the origin of this species, providing further evidence that biased LTR retrotransposon activity has played a major role in shaping the chromatin and DNA landscape of the sunflower genome. Although our findings on LTR retrotransposon age and structure could be influenced by the selection of the BAC clones analyzed, a global analysis of random sequence reads indicates that the evolutionary patterns described herein apply to the sunflower genome as a whole.