Hybridisation and chloroplast capture between distinct Themeda triandra lineages in Australia

Ecotypes are distinct populations within a species that are adapted to specific environmental conditions. Understanding how these ecotypes become established, and how they interact when reunited, is fundamental to elucidating how ecological adaptations are maintained. This study focuses on Themeda triandra, a dominant grassland species across Asia, Africa and Australia. It is the most widespread plant in Australia, where it has distinct ecotypes that are usually restricted to either wetter and cooler coastal regions or the drier and hotter interior. We generate a reference genome for T. triandra and use whole genome sequencing for over 80 Themeda accessions to reconstruct the evolutionary history of T. triandra and related taxa. Organelle phylogenies confirm that Australia was colonized by T. triandra twice, with the division between ecotypes predating their arrival in Australia. The nuclear genome provides evidence of differences in the dominant ploidal level and gene-flow among the ecotypes. In northern Queensland there appears to be a hybrid zone between ecotypes with admixed nuclear genomes and shared chloroplast haplotypes. Conversely, in the cracking claypans of Western Australia, there is cytonuclear discordance with individuals possessing the coastal chloroplast and interior clade nuclear genome. This chloroplast capture is potentially a result of adaptive introgression, with selection detected in the rpoC2 gene which is associated with water use efficiency. The reason that T. triandra is the most widespread plant in Australia appears to be a result of distinct ecotypic genetic variation and genome duplication, with the importance of each depending on the geographic scale considered.

Hybrid evolution repeats itself across environmental contexts in Texas sunflowers (Helianthus)

To what extent is evolution repeatable? Little is known about whether the evolution of hybrids is more (or less) repeatable than that of nonhybrids. We used field experimental evolution in annual sunflowers (Helianthus) in Texas to ask the extent to which hybrid evolution is repeatable across environments compared to nonhybrid controls. We created hybrids between Helianthus annuus (L.) and H. debilis (Nutt.) and grew plots of both hybrids and nonhybrid controls through eight generations at three sites in Texas. We collected seeds from each generation and grew each generation × treatment × home site combination at two final common gardens. We estimated the strength and direction of evolution in terms of fitness and 24 traits, tested for repeated versus nonrepeated evolution, and assessed overall phenotypic evolution across lineages and in relation to a locally adapted phenotype. Hybrids consistently evolved higher fitness over time, while controls did not, although trait evolution varied in strength across home sites. Repeated evolution was more evident in hybrids versus nonhybrid controls, and hybrid evolution was often in the direction of the locally adapted phenotype. Our findings have implications for both the nature of repeatability in evolution and the contribution of hybridization to evolution across environmental contexts.

Molecular signatures of parallel adaptive divergence causing reproductive isolation and speciation across two genera

Parallel evolution of reproductive isolation (PERI) provides strong evidence for natural selection playing a fundamental role in the origin of species. However, PERI has been rarely demonstrated for well established species drawn from different genera. In particular, parallel molecular signatures for the same genes in response to similar habitat divergence in such different lineages is lacking. Here, based on whole-genome sequencing data, we first explore the speciation process in two sister species of Carpinus (Betulaceae) in response to divergence for temperature and soil-iron concentration in habitats they occupy in northern and southwestern China, respectively. We then determine whether parallel molecular mutations occur during speciation in this pair of species and also in another sister-species pair of the related genus, Ostryopsis, which occupy similarly divergent habitats in China. We show that gene flow occurred during the origin of both pairs of sister species since approximately 9.8 or approximately 2 million years ago, implying strong natural selection during divergence. Also, in both species pairs we detected concurrent positive selection in a gene (LHY) for flowering time and in two paralogous genes (FRO4 and FRO7) of a gene family known to be important for iron tolerance. These changes were in addition to changes in other major genes related to these two traits. The different alleles of these particular candidate genes possessed by the sister species of Carpinus were functionally tested and indicated likely to alter flowering time and iron tolerance as previously demonstrated in the pair of Ostryopsis sister species. Allelic changes in these genes may have effectively resulted in high levels of prezygotic reproductive isolation to evolve between sister species of each pair. Our results show that PERI can occur in different genera at different timescales and involve similar signatures of molecular evolution at genes or paralogues of the same gene family, causing reproductive isolation as a consequence of adaptation to similarly divergent habitats.

•

PERI provides strong evidence for natural selection playing a fundamental role in the origin of species

•

PERI is rarely demonstrated for well-established species drawn from different genera

•

We detected PERI across two genera (Carpinus and Ostryopsis) in the family Betulaceae

•

PERI can occur in different genera at different timescales and involve molecular signatures at similar pathways

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.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

The Evolutionary Genomics of Serpentine Adaptation

Serpentine barrens are among the most challenging settings for plant life. Representing a perfect storm of hazards, serpentines consist of broadly skewed elemental profiles, including abundant toxic metals and low nutrient contents on drought-prone, patchily distributed substrates. Accordingly, plants that can tolerate the challenges of serpentine have fascinated biologists for decades, yielding important insights into adaptation to novel ecologies through physiological change. Here we highlight recent progress from studies which demonstrate the power of serpentine as a model for the genomics of adaptation. Given the moderate - but still tractable - complexity presented by the mix of hazards on serpentine, these venues are well-suited for the experimental inquiry of adaptation both in natural and manipulated conditions. Moreover, the island-like distribution of serpentines across landscapes provides abundant natural replicates, offering power to evolutionary genomic inference. Exciting recent insights into the genomic basis of serpentine adaptation point to a partly shared basis that involves sampling from common allele pools available from retained ancestral polymorphism or via gene flow. However, a lack of integrated studies deconstructing complex adaptations and linking candidate alleles with fitness consequences leaves room for much deeper exploration. Thus, we still seek the crucial direct link between the phenotypic effect of candidate alleles and their measured adaptive value - a prize that is exceedingly rare to achieve in any study of adaptation. We expect that closing this gap is not far off using the promising model systems described here.

Linking micro and macroevolution in the presence of migration

Understanding macroevolutionary patterns is central to evolutionary biology. This involves the process of divergence within a species, which starts at the microevolutionary level, for instance, when two subpopulations evolve towards different phenotypic optima. The speed at which these optima are reached is controlled by the degree of stabilising selection, which pushes the mean trait towards different optima in the different subpopulations, and ongoing migration that pulls the mean phenotype away from that optimum. Traditionally, macro phenotypic evolution is modelled by directional selection processes, but these models usually ignore the role of migration within species. Here, our goal is to reconcile the processes of micro and macroevolution by modelling migration as part of the speciation process. More precisely, we introduce an Ornstein-Uhlenbeck (OU) model where migration happens between two subpopulations within a branch of a phylogeny and this migration decreases over time as it happens during speciation. We then use this model to study the evolution of trait means along a phylogeny, as well as the way phenotypic disparity between species changes with successive epochs. We show that ignoring the effect of migration in sampled time-series data biases significantly the estimation of the selective forces acting upon it. We also show that migration decreases the expected phenotypic disparity between species and we analyse the effect of migration in the particular case of niche filling. We further introduce a method to jointly estimate selection and migration from time-series data. Our model extends traditional quantitative genetics results of selection and migration from a microevolutionary time frame to multiple speciation events at a macroevolutionary scale. Our results further support that not accounting for gene flow has important consequences in inferences at both the micro and macroevolutionary scale.

Parallel Speciation of Wild Rice Associated with Habitat Shifts

The occurrence of parallel speciation strongly implies the action of natural selection. However, it is unclear how general a phenomena parallel speciation is since it was only shown in a small number of animal species. In particular, the adaptive process and mechanisms underlying the process of parallel speciation remain elusive. Here, we used an integrative approach incorporating population genomics, common garden, and crossing experiments to investigate parallel speciation of the wild rice species Oryza nivara from O. rufipogon. We demonstrated that O. nivara originated multiple times from different O. rufipogon populations and revealed that different O. nivara populations have evolved similar phenotypes under divergent selection, a reflection of recurrent local adaptation of ancient O. rufipogon populations to dry habitats. Almost completed premating isolation was detected between O. nivara and O. rufipogon in the absence of any postmating barriers between and within these species. These results suggest that flowering time is a “magic” trait that contributes to both local adaptation and reproductive isolation in the origin of wild rice species. Our study thus demonstrates a convincing case of parallel ecological speciation as a consequence of adaptation to new environments.

Phylogeographic analysis suggests two origins for the riparian azalea Rhododendron indicum (L.) Sweet

Ecological speciation is an important factor in the diversification of plants. The distribution of the woody species Rhododendron indicum, which grows along rivers and is able to withstand water flow when rivers flood (i.e. it is a rheophyte), is disjunct, in contrast to the widespread distribution of its relative, Rhododendron kaempferi. This study aimed to elucidate the phylogenetic relationships between R. indicum and R. kaempferi and the evolutionary processes that gave rise to them. The sequences of three non-coding chloroplast DNA regions (total length 1977 bp) were obtained from 21 populations covering the ranges of the two species. In addition, genome-wide SNPs were genotyped from 20 populations using a genotyping by sequencing method. Leaf morphologies were measured for eight representative populations. Two chloroplast DNA haplotypes, which were detected in R. indicum, were shared between the two species. Genome-wide SNPs identified two lineages in R. indicum and these lineages did not constitute a monophyletic group. Each of these two lineages was related to geographically close populations of R. kaempferi. Leaf morphology, which is a characteristic feature in rheophytes, was not differentiated between the two lineages in R. indicum. The morphological similarity between the two heterogeneous lineages may be a result of parallel evolution from R. kaempferi or of introgressive hybridization between the species due to strong selective pressure imposed by flooding.

Can parallel ecological speciation be detected with phylogenetic analyses?

Distinguishing parallel divergence from vicariance scenarios is crucial to establish the determinism of natural selection in the formation of new species. It has been proposed that phylogenetic methods can be used to disentangle a single origin in allopatry and multiple origins in sympatry for ecological speciation. However, a key issue is to what extent introgression in a patchy environment may complicate the distinction between both origins through the analysis of variation at neutral markers. While divergence at environmentally-based selected loci retains the initial correlation with ecological variables, such association may be lost at neutral loci unlinked to any selected locus. Thus, neutral divergence might reflect in the long-term the molecular fingerprint of isolation by distance regardless of the model of speciation considered, and a question arises as to whether phylogenetic analyses of neutral markers are able or not to retrieve the signals acquired in the ancestral populations. Here, we use computer simulations to show that the detection of the original signal using a phylogenetic method strongly depends on the migration rates among populations. Recombination accelerates the loss of the initial phylogenetic signal, but this effect is rather small compared with the effect of migration, and only substantial when recombination is very large. For model species with reduced gene flow between distant populations and between populations adapted to different habitats, the phylogenetic approach is able to distinguish a single origin in allopatry from multiple origins in sympatry.

RADseq provides evidence for parallel ecotypic divergence in the autotetraploid Cochlearia officinalis in Northern Norway

Speciation encompasses a continuum over time from freely interbreeding populations to reproductively isolated species. Along this process, ecotypes - the result of local adaptation - may be on the road to new species. We investigated whether three autotetraploid Cochlearia officinalis ecotypes, adapted to different habitats (beach, estuary, spring), are genetically differentiated and result from parallel ecotypic divergence in two distinct geographical regions. We obtained genetic data from thousands of single nucleotide polymorphisms (SNPs) from restriction-site associated DNA sequencing (RADseq) and from six microsatellite markers for 12 populations to assess genetic divergence at ecotypic, geographic and population level. The genetic patterns support differentiation among ecotypes as suggested by morphology and ecology. The data fit a scenario where the ancestral beach ecotype has recurrently and polytopically given rise to the estuary and spring ecotypes. Several ecologically-relevant loci with consistent non-random segregating patterns are identified across the recurrent origins, in particular around genes related to salt stress. Despite being ecologically distinct, the Cochlearia ecotypes still represent an early stage in the process of speciation, as reproductive isolation has not (yet) developed. A sequenced annotated genome is needed to specifically target candidate genes underlying local adaptation.

The origins of reproductive isolation in plants

Reproductive isolation in plants occurs through multiple barriers that restrict gene flow between populations, but their origins remain uncertain. Work in the past decade has shown that postpollination barriers, such as the failure to form hybrid seeds or sterility of hybrid offspring, are often less strong than prepollination barriers. Evidence implicates multiple evolutionary forces in the origins of reproductive barriers, including mutation, stochastic processes and natural selection. Although adaptation to different environments is a common element of reproductive isolation, genomic conflicts also play a role, including female meiotic drive. The genetic basis of some reproductive barriers, particularly flower colour influencing pollinator behaviour, is well understood in some species, but the genetic changes underlying many other barriers, especially pollen-stylar interactions, are largely unknown. Postpollination barriers appear to accumulate at a faster rate in annuals compared with perennials, due in part to chromosomal rearrangements. Chromosomal changes can be important isolating barriers in themselves but may also reduce the recombination of genes contributing to isolation. Important questions for the next decade include identifying the evolutionary forces responsible for chromosomal rearrangements, determining how often prezygotic barriers arise due to selection against hybrids, and establishing the relative importance of genomic conflicts in speciation.

Multiple shifts to different pollinators fuelled rapid diversification in sexually deceptive Ophrys orchids

Episodes of rapid speciation provide unique insights into evolutionary processes underlying species radiations and patterns of biodiversity. Here we investigated the radiation of sexually deceptive bee orchids (Ophrys). Based on a time-calibrated phylogeny and by means of ancestral character reconstruction and divergence time estimation, we estimated the tempo and mode of this radiation within a state-dependent evolutionary framework. It appears that, in the Pleistocene, the evolution of Ophrys was marked by episodes of rapid diversification coinciding with shifts to different pollinator types: from wasps to Eucera bees to Andrena and other bees. An abrupt increase in net diversification rate was detected in three clades. Among these, two phylogenetically distant lineages switched from Eucera to Andrena and other bees in a parallel fashion and at about the same time in their evolutionary history. Lack of early radiation associated with the evolution of the key innovation of sexual deception suggests that Ophrys diversification was mainly driven by subsequent ecological opportunities provided by the exploitation of novel pollinator groups, encompassing many bee species slightly differing in their sex pheromone communication systems, and by spatiotemporal fluctuations in the pollinator mosaic.

Modeling the two-locus architecture of divergent pollinator adaptation: how variation in SAD paralogs affects fitness and evolutionary divergence in sexually deceptive orchids

Divergent selection by pollinators can bring about strong reproductive isolation via changes at few genes of large effect. This has recently been demonstrated in sexually deceptive orchids, where studies (1) quantified the strength of reproductive isolation in the field; (2) identified genes that appear to be causal for reproductive isolation; and (3) demonstrated selection by analysis of natural variation in gene sequence and expression. In a group of closely related Ophrys orchids, specific floral scent components, namely n-alkenes, are the key floral traits that control specific pollinator attraction by chemical mimicry of insect sex pheromones. The genetic basis of species-specific differences in alkene production mainly lies in two biosynthetic genes encoding stearoyl–acyl carrier protein desaturases (SAD) that are associated with floral scent variation and reproductive isolation between closely related species, and evolve under pollinator-mediated selection. However, the implications of this genetic architecture of key floral traits on the evolutionary processes of pollinator adaptation and speciation in this plant group remain unclear. Here, we expand on these recent findings to model scenarios of adaptive evolutionary change at SAD2 and SAD5, their effects on plant fitness (i.e., offspring number), and the dynamics of speciation. Our model suggests that the two-locus architecture of reproductive isolation allows for rapid sympatric speciation by pollinator shift; however, the likelihood of such pollinator-mediated speciation is asymmetric between the two orchid species O. sphegodes and O. exaltata due to different fitness effects of their predominant SAD2 and SAD5 alleles. Our study not only provides insight into pollinator adaptation and speciation mechanisms of sexually deceptive orchids but also demonstrates the power of applying a modeling approach to the study of pollinator-driven ecological speciation.

Incipient radiation within the dominant Hawaiian tree Metrosideros polymorpha

Although trees comprise a primary component of terrestrial species richness, the drivers and temporal scale of divergence in trees remain poorly understood. We examined the landscape-dominant tree, Metrosideros polymorpha, for variation at nine microsatellite loci across 23 populations on young Hawai'i Island, sampling each of the island's five varieties throughout its full geographic range. For four varieties, principal coordinate analysis revealed strong clustering of populations by variety across the 10 430 km(2) island, indicating partitioning of the species into multiple evolutionarily significant units. The single island-endemic form, riparian var. newellii, showed especially strong differentiation from other varieties despite occurring in sympatry with other varieties and likely evolved from a bog form on the oldest volcano, Kohala, within the past 500 000 years. Along with comparable riparian forms on other Pacific Islands, var. newellii appears to represent parallel incipient ecological speciation within Metrosideros. Greater genetic distance among the more common varieties on the oldest volcano and an inverse relationship between allelic diversity and substrate age appear consistent with colonization of Hawai'i Island by older, partially diverged varieties followed by increased hybridization among varieties on younger volcanoes. This study demonstrates that broad population-level sampling is required to uncover patterns of diversification within a ubiquitous and long-lived tree species. Hawaiian Metrosideros appears to be a case of incipient radiation in trees and thus should be useful for studies of divergence and the evolution of reproductive isolating barriers at the early stages 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.

Parallel speciation or long-distance dispersal? Lessons from seaweeds (Fucus) in the Baltic Sea

Parallel evolution has been invoked as a forceful mechanism of ecotype and species formation in many animal taxa. However, parallelism may be difficult to separate from recently monophyletically diverged species that are likely to show complex genetic relationships as a result of considerable shared ancestral variation and secondary hybridization in local areas. Thus, species' degrees of reproductive isolation, barriers to dispersal and, in particular, limited capacities for long-distance dispersal will affect demographical structures underlying mechanisms of divergent evolution. Here, we used nine microsatellite DNA markers to study intra- and interspecific genetic diversity of two recently diverged species of brown macroalgae, Fucus radicans (L. Bergström & L. Kautsky) and F. vesiculosus (Linnaeus), in the Baltic Sea. We further performed biophysical modelling to identify likely connectivity patterns influencing the species' genetic structures. For each species, we found intraspecific contrasting patterns of clonality incidence and population structure. In addition, strong genetic differentiation between the two species within each locality supported the existence of two distinct evolutionary lineages (FST  = 0.15-0.41). However, overall genetic clustering analyses across both species' populations revealed that all populations from one region (Estonia) were more genetically similar to each other than to their own taxon from the other two regions (Sweden and Finland). Our data support a hypothesis of parallel speciation. Alternatively, Estonia may be the ancestral source of both species, but is presently isolated by oceanographic barriers to dispersal. Thus, a limited gene flow in combination with genetic drift could have shaped the seemingly parallel structure.

Quantitative traits and mode of speciation in Martinique anoles

We investigate extensive quantitative trait variation (dewlap hue, colour pattern, dorsum hue, body proportions and scalation) in the Martinique anole across eight transects representing nascent parapatric ecological speciation, nascent allopatric speciation and allopatric divergence without sufficient genetic structure to suggest speciation. Quantitative trait divergence can be extremely large between adjacent sets of populations, but with one exception that this is associated with difference in habitat rather than past allopatry. Nascent ecological speciation shows the greatest level of quantitative trait divergence across all character sets including those implicated in natural, as well as sexual selection. The sole example of nascent allopatric speciation is associated with fairly strong quantitative trait divergence among most character sets, but not the set most implicated in natural (rather than sexual) selection. The role of sexual selection in ecological speciation is discussed, both in terms of female choice with assortative mating and male-male competition with condition-dependant sexual signals.

The Role of Parasitism in Adaptive Radiations—When Might Parasites Promote and When Might They Constrain Ecological Speciation?

Research on speciation and adaptive radiation has flourished during the past decades, yet factors underlying initiation of reproductive isolation often remain unknown. Parasites represent important selective agents and have received renewed attention in speciation research. We review the literature on parasite-mediated divergent selection in context of ecological speciation and present empirical evidence for three nonexclusive mechanisms by which parasites might facilitate speciation: reduced viability or fecundity of immigrants and hybrids, assortative mating as a pleiotropic by-product of host adaptation, and ecologically-based sexual selection. We emphasise the lack of research on speciation continuums, which is why no study has yet made a convincing case for parasite driven divergent evolution to initiate the emergence of reproductive isolation. We also point interest towards selection imposed by single versus multiple parasite species, conceptually linking this to strength and multifariousness of selection. Moreover, we discuss how parasites, by manipulating behaviour or impairing sensory abilities of hosts, may change the form of selection that underlies speciation. We conclude that future studies should consider host populations at variable stages of the speciation process, and explore recurrent patterns of parasitism and resistance that could pinpoint the role of parasites in imposing the divergent selection that initiates ecological speciation.