Divergent natural selection drives the evolution of reproductive isolation in an Australian wildflower

Genomics of hybrid parallel origin in Aquilegia ecalcarata

BACKGROUND

The parallel evolution of similar traits or species provides strong evidence for the role of natural selection in evolution. Traits or species that evolved repeatedly can be driven by separate de novo mutations or interspecific gene flow. Although parallel evolution has been reported in many studies, documented cases of parallel evolution caused by gene flow are scarce by comparison. Aquilegia ecalcarata and A. kansuensis belong to the genus of Aquilegia, and are the closest related sister species. Mutiple origins of A. ecalcarata have been reported in previous studies, but whether they have been driven by separate de novo mutations or gene flow remains unclear.

RESULTS

In this study, We conducted genomic analysis from 158 individuals of two repeatedly evolving pairs of A. ecalcarata and A. kansuensis. All samples were divided into two distinct clades with obvious geographical distribution based on phylogeny and population structure. Demographic modeling revealed that the origin of the A. ecalcarata in the Eastern of China was caused by gene flow, and the Eastern A. ecalcarata occurred following introgression from Western A. ecalcarata population. Analysis of Treemix and D-statistic also revealed that a strong signal of gene flow was detected from Western A. ecalcarata to Eastern A. ecalcarata. Genetic divergence and selective sweep analyses inferred parallel regions of genomic divergence and identified many candidate genes associated with ecologically adaptive divergence between species pair. Comparative analysis of parallel diverged regions and gene introgression confirms that gene flow contributed to the parallel evolution of A. ecalcarata.

CONCLUSIONS

Our results further confirmed the multiple origins of A. ecalcarata and highlighted the roles of gene flow. These findings provide new evidence for parallel origin after hybridization as well as insights into the ecological adaptation mechanisms underlying the parallel origins of species.

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

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

Population genetics and independently replicated evolution of predator-associated burst speed ecophenotypy in mosquitofish

Many species show replicated ecophenotypy due to recurring patterns of natural selection. Based on the presence or absence of pursuit predators, at least 17 species of fish repeatedly differentiated in body shape in a manner that increases burst swimming speed and the likelihood of predator escape. The predator-associated burst speed (PABS) ecophenotype is characterized by a small head and trunk and enlarged caudal region. Mechanisms promoting replicated phenotype-environment association include selection (without evolution), a single instance of adaptive evolution followed by biased habitat occupation, repeated instances of local adaptation, or adaptive phenotypic plasticity. Common garden rearing of mosquitofish, Gambusia affinis, demonstrated a likely heritable basis for PABS phenotypy, but it is unknown whether populations are otherwise genetically distinct or whether replicated ecophenotypy represents a single or replicated instances of adaptation. To genetically characterize the populations and test hypotheses of single or multiple adaptations, we characterized variation in 12 polymorphic DNA microsatellites in the previously studied G. affinis populations. Populations were genetically distinct by multilocus analysis, exhibited high allelic diversity, and were heterozygote deficient, which effects were attributed to G. affinis's shoaling nature and habitat patchiness. Genetic and phenotypic distances among populations were correlated for non-PABS but not PABS morphology. Multilocus analysis demonstrated ecophenotype polyphyly and scattered multivariate genetic structure which support only the replicated-adaptation model. As all of the diverse tests performed demonstrated lack of congruence between patterns of molecular genetic and PABS differentiation, it is likely that divergent natural selection drove multiple instances of adaptive evolution.

Natural hybridization among three Rhododendron species (Ericaceae) revealed by morphological and genomic evidence

BACKGROUND

Natural hybridization can influence the adaptive response to selection and accelerate species diversification. Understanding the composition and structure of hybrid zones may elucidate patterns of hybridization processes that are important to the formation and maintenance of species, especially for taxa that have experienced rapidly adaptive radiation. Here, we used morphological traits, ddRAD-seq and plastid DNA sequence data to investigate the structure of a Rhododendron hybrid zone and uncover the hybridization patterns among three sympatric and closely related species.

RESULTS

Our results show that the hybrid zone is complex, where bi-directional hybridization takes place among the three sympatric parental species: R. spinuliferum, R. scabrifolium, and R. spiciferum. Hybrids between R. spinuliferum and R. spiciferum (R. ×duclouxii) comprise multiple hybrid classes and a high proportion of F₁ generation hybrids, while a novel hybrid taxon between R. spinuliferum and R. scabrifolium dominated the F₂ generation, but no backcross individuals were detected. The hybrid zone showed basically coincident patterns of population structure between genomic and morphological data.

CONCLUSIONS

Natural hybridization exists among the three Rhododendron species in the hybrid zone, although patterns of hybrid formation vary between hybrid taxa, which may result in different evolutionary outcomes. This study represents a unique opportunity to dissect the ecological and evolutionary mechanisms associated with adaptive radiation of Rhododendron species in a biodiversity hotspot.

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.

Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio

The independent and repeated adaptation of populations to similar environments often results in the evolution of similar forms. This phenomenon creates a strong correlation between phenotype and environment and is referred to as parallel evolution. However, we are still largely unaware of the dynamics of parallel evolution, as well as the interplay between phenotype and genotype within natural systems. Here, we examined phenotypic and genotypic parallel evolution in multiple parapatric Dune‐Headland coastal ecotypes of an Australian wildflower, Senecio lautus. We observed a clear trait‐environment association in the system, with all replicate populations having evolved along the same phenotypic evolutionary trajectory. Similar phenotypes have arisen via mutational changes occurring in different genes, although many share the same biological functions. Our results shed light on how replicated adaptation manifests at the phenotypic and genotypic levels within populations, and highlight S. lautus as one of the most striking cases of phenotypic parallel evolution in nature.

Divergence in flowering time is a major component contributing to reproductive isolation between two wild rice species (Oryza rufipogon and O. nivara)

It is of critical importance for our understanding of speciation process to determine the forms of reproductive isolation and their relative importance in species divergence. Oryza nivara and O. rufipogon are direct ancestors of Asian cultivated rice and a progenitor-daughter species pair. Investigating the reproductive isolation between them provides insights into plant speciation and helps understanding of the rice domestication. Here, we quantitatively measured the major components of reproductive isolation between the two species based on common garden and crossing experiments for three pairs of sympatric populations in Nepal, Cambodia and Laos. We revealed significant differences in the flowering times between species pairs, with O. nivara flowering much earlier than O. rufipogon. A very weak reduction in seed set but no reduction in F1 viability and fertility were detected for the crosses between species relative to those within species. Moreover, we detected asymmetrical compatibility between species and found that emasculation significantly decreased pollination success in O. nivara but not in O. rufipogon. Our study demonstrates that the divergence between O. nivara and O. rufipogon is maintained almost entirely by the difference in flowering times and suggests that differential flowering times contribute to both habitat preferences and reproductive isolation between species.

Loss of ecologically important genetic variation in late generation hybrids reveals links between adaptation and speciation

Adaptation to contrasting environments occurs when advantageous alleles accumulate in each population, but it remains largely unknown whether these same advantageous alleles create genetic incompatibilities that can cause intrinsic reproductive isolation leading to speciation. Identifying alleles that underlie both adaptation and reproductive isolation is further complicated by factors such as dominance and genetic interactions among loci, which can affect both processes differently and obscure potential links between adaptation and speciation. Here, we use a combination of field and glasshouse experiments to explore the connection between adaptation and speciation while accounting for dominance and genetic interactions. We created a hybrid population with equal contributions from four contrasting ecotypes of Senecio lautus (Asteraceae), which produced hybrid genomes both before (F1 hybrid generation) and after (F4 hybrid generation) recombination among the parental ecotypes. In the glasshouse, plants in the second generation (F2 hybrid generation) showed reduced fitness as a loss of fertility. However, fertility was recovered in subsequent generations, suggesting that genetic variation underlying the fitness reduction was lost in subsequent generations. To quantify the effects of losing genetic variation at the F2 generation on the fitness of later generation hybrids, we used a reciprocal transplant to test for fitness differences between parental ecotypes, and F1 and F4 hybrids in all four parental habitats. Compared to the parental ecotypes and F1 hybrids, variance in F4 hybrid fitness was lower, and lowest in habitats that showed stronger native-ecotype advantage, suggesting that stronger natural selection for the native ecotype reduced fitness variation in the F4 hybrids. Fitness trade-offs that were present in the parental ecotypes and F1 hybrids were absent in the F4 hybrid. Together, these results suggest that the genetic variation lost after the F2 generation was likely associated with both adaptation and intrinsic reproductive isolation among ecotypes from contrasting habitats.

Senecio as a model system for integrating studies of genotype, phenotype and fitness

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

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.

Natural selection drives leaf divergence in experimental populations of Senecio lautus under natural conditions

Leaf morphology is highly variable both within and between plant species. This study employs a combination of common garden and reciprocal transplant experiments to determine whether differences in leaf shape between Senecio lautus ecotypes has evolved as an adaptive response to divergent ecological conditions.We created a synthetic population of hybrid genotypes to segregate morphological variation between three ecotypes and performed reciprocal transplants where this hybrid population was transplanted into the three adjacent native environments. We measured nine leaf morphology traits across the experimental and natural populations at these sites.We found significant divergence in multivariate leaf morphology toward the native character in each environment, suggesting environmental conditions at each site exert selective pressure that results in a phenotypic shift toward the local phenotype of the wild populations.These associations suggest that differences in leaf morphology between S. lautus ecotypes have arisen as a result of divergent selection on leaf shape or associated traits that confer an adaptive advantage in each environment, which has led to the formation of morphologically distinct ecotypes.

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.

Environmentally induced development costs underlie fitness tradeoffs

Local adaptation can lead to genotype-by-environment interactions, which can create fitness tradeoffs in alternative environments, and govern the distribution of biodiversity across geographic landscapes. Exploring the ecological circumstances that promote the evolution of fitness tradeoffs requires identifying how natural selection operates and during which ontogenetic stages natural selection is strongest. When organisms disperse to areas outside their natural range, tradeoffs might emerge when organisms struggle to reach key life history stages, or alternatively, die shortly after reaching life history stages if there are greater risks of mortality associated with costs to developing in novel environments. We used multiple populations from four ecotypes of an Australian native wildflower (Senecio pinnatifolius) in reciprocal transplants to explore how fitness tradeoffs arise across ontogeny. We then assessed whether the survival probability for plants from native and foreign populations was contingent on reaching key developmental stages. We found that fitness tradeoffs emerged as ontogeny progressed when native plants were more successful than foreign plants at reaching seedling establishment and maturity. Native and foreign plants that failed to reach seedling establishment died at the same rate, but plants from foreign populations died quicker than native plants after reaching seedling establishment, and died quicker regardless of whether they reached sexual maturity or not. Development rates were similar for native and foreign populations, but changed depending on the environment. Together, our results suggest that natural selection for environment-specific traits early in life history created tradeoffs between contrasting environments. Plants from foreign populations were either unable to develop to seedling establishment, or they suffered increased mortality as a consequence of reaching seedling establishment. The observation of tradeoffs together with environmentally dependent changes in development rate suggest that foreign environments induce organisms to develop at a rate different from their native habitat, incurring consequences for lifetime fitness and population divergence.

Evolution of Genetic Variance during Adaptive Radiation

Genetic correlations between traits can concentrate genetic variance into fewer phenotypic dimensions that can bias evolutionary trajectories along the axis of greatest genetic variance and away from optimal phenotypes, constraining the rate of evolution. If genetic correlations limit adaptation, rapid adaptive divergence between multiple contrasting environments may be difficult. However, if natural selection increases the frequency of rare alleles after colonization of new environments, an increase in genetic variance in the direction of selection can accelerate adaptive divergence. Here, we explored adaptive divergence of an Australian native wildflower by examining the alignment between divergence in phenotype mean and divergence in genetic variance among four contrasting ecotypes. We found divergence in mean multivariate phenotype along two major axes represented by different combinations of plant architecture and leaf traits. Ecotypes also showed divergence in the level of genetic variance in individual traits and the multivariate distribution of genetic variance among traits. Divergence in multivariate phenotypic mean aligned with divergence in genetic variance, with much of the divergence in phenotype among ecotypes associated with changes in trait combinations containing substantial levels of genetic variance. Overall, our results suggest that natural selection can alter the distribution of genetic variance underlying phenotypic traits, increasing the amount of genetic variance in the direction of natural selection and potentially facilitating rapid adaptive divergence during an adaptive radiation.