PHYLOGENETIC EVIDENCE FOR MULTIPLE SYMPATRIC ECOLOGICAL DIVERSIFICATION IN A MARINE SNAIL

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

A Simulation Study of the Ecological Speciation Conditions in the Galician Marine Snail Littorina saxatilis

In the last years, the interest in evolutionary divergence at small spatial scales has increased and so did the study of speciation caused by ecologically based divergent natural selection. The evolutionary interplay between gene flow and local adaptation can lead to low-dispersal locally adapted specialists. When this occurs, the evolutionary interplay between gene flow and local adaptation could eventually lead to speciation. The L. saxatilis system consists of two ecotypes displaying a microhabitat-associated intraspecific dimorphism along the wave-exposed rocky shores of Galicia. Despite being a well-known system, the dynamics of the ecotype formation remain unclear and cannot be studied from empirical evidence alone. In this study, individual-based simulations were used to incorporate relevant ecological, spatial, and genetic information, to check different evolutionary scenarios that could evolve non-random mating preferences and finally may facilitate speciation. As main results, we observed the evolution of intermediate values of choice which matches the estimates from empirical data of L. saxatilis in Galician shores and coincides with previous theoretical outcomes. Also, the use of the mating correlation as a proxy for assortative mating led to spuriously inferring greater reproductive isolation in the middle habitat than in the others, which does not happen when directly considering the choice values from the simulations. We also corroborate the well-known fact that the occurrence of speciation is influenced by the strength of selection. Taken together, this means, also according to other L. saxatilis systems, that speciation is not an immediate consequence of local divergent selection and mating preferences, but a fine tuning among several factors including the ecological conditions in the shore levels, the selection strength, the mate choice stringency, and cost to choosiness. The L. saxatilis system could correspond to a case of incomplete reproductive isolation, where the choice intensity is intermediate and local adaptation within the habitat is strong. These results support previous interpretations of the L. saxatilis model system and indicate that further empirical studies would be interesting to test whether the mate choice mechanism functions as a similarity-like mechanism as has been shown in other littorinids.

Genome-Wide Investigation of the Multiple Origins Hypothesis for Deep-Spawning Kokanee Salmon (Oncorhynchus nerka) across its Pan-Pacific Distribution

Salmonids have emerged as important study systems for investigating molecular processes underlying parallel evolution given their tremendous life history variation. Kokanee, the resident form of anadromous sockeye salmon (Oncorhynchus nerka), have evolved multiple times across the species' pan-Pacific distribution, exhibiting multiple reproductive ecotypes including those that spawn in streams, on lake-shores, and at lake depths >50 m. The latter has only been detected in 5 locations in Japan and British Columbia, Canada. Here, we investigated the multiple origins hypothesis for deep-spawning kokanee, using 9721 single nucleotide polymorphisms distributed across the genome analyzed for the vast majority of known populations in Japan (Saiko Lake) and Canada (Anderson, Seton, East Barrière Lakes) relative to stream-spawning populations in both regions. We detected 397 outlier loci, none of which were robustly identified in paired-ecotype comparisons in Japan and Canada independently. Bayesian clustering and principal components analyses based on neutral loci revealed 6 distinct clusters, largely associated with geography or translocation history, rather than ecotype. Moreover, a high level of divergence between Canadian and Japanese populations, and between deep- and stream-spawning populations regionally, suggests the deep-spawning ecotype independently evolved on the 2 continents. On a finer level, Japanese kokanee populations exhibited low estimates of heterozygosity, significant levels of inbreeding, and reduced effective population sizes relative to Canadian populations, likely associated with transplantation history. Along with preliminary evidence for hybridization between deep- and stream-spawning ecotypes in Saiko Lake, these findings should be considered within the context of on-going kokanee fisheries management in Japan.

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.

Parallel introgression, not recurrent emergence, explains apparent elevational ecotypes of polyploid Himalayan snowtrout

The recurrence of similar evolutionary patterns within different habitats often reflects parallel selective pressures acting upon either standing or independently occurring genetic variation to produce a convergence of phenotypes. This interpretation (i.e. parallel divergences within adjacent streams) has been hypothesized for drainage-specific morphological 'ecotypes' observed in polyploid snowtrout (Cyprinidae: Schizothorax). However, parallel patterns of differential introgression during secondary contact are a viable alternative hypothesis. Here, we used ddRADseq (N = 35 319 de novo and N = 10 884 transcriptome-aligned SNPs), as derived from Nepali/Bhutanese samples (N = 48 each), to test these competing hypotheses. We first employed genome-wide allelic depths to derive appropriate ploidy models, then a Bayesian approach to yield genotypes statistically consistent under the inferred expectations. Elevational 'ecotypes' were consistent in geometric morphometric space, but with phylogenetic relationships at the drainage level, sustaining a hypothesis of independent emergence. However, partitioned analyses of phylogeny and admixture identified subsets of loci under selection that retained genealogical concordance with morphology, suggesting instead that apparent patterns of morphological/phylogenetic discordance are driven by widespread genomic homogenization. Here, admixture occurring in secondary contact effectively 'masks' previous isolation. Our results underscore two salient factors: (i) morphological adaptations are retained despite hybridization and (ii) the degree of admixture varies across tributaries, presumably concomitant with underlying environmental or anthropogenic factors.

Proteomic analysis of F1 hybrids and intermediate variants in a littorina saxatilis hybrid zone

Proteomic analysis was carried out on the Crab (upper-shore) and Wave (lower-shore) ecotypes of Littorina saxatilis from a hybrid zone at Silleiro Cape, Spain. Proteome profiles of individual snails were obtained. Protein expression in F₁ hybrid snails bred in the laboratory and snails with intermediate shell phenotypes collected from the mid-shore were compared with Crab and Wave ecotypes using analytical approaches used to study dominance. Multivariate analysis over many protein spots showed that the F₁ snails are distinct from both ecotypes but closer to the Wave ecotype. The intermediate snails are highly variable, some closer to the Crab and others to the Wave ecotype. Considered on a protein by protein basis, some proteins are significantly closer in expression to the Crab and others to the Wave ecotype for both F₁ and intermediate snails. Furthermore, a significant majority of proteins were closer in expression to the Wave ecotype for the F₁, consistent with the multivariate analysis. No such significant majority toward either the Crab or Wave ecotype was observed for the intermediate snails. The closer similarity of F₁ and Wave ecotype expression patterns could be the result of similar selective pressures in the similar mid-shore and low-shore environments. For a significantly larger number of proteins, intermediate snails were closer in expression to the ecotype having the lower expression, for both Crab and Wave ecotypes. This is somewhat unexpected as lower expression might be expected to be an indication of impairment of function and lower fitness. Proteomic analysis could be important for the identification of candidate proteins useful for gaining improved understanding of adaptation and barriers to gene flow in hybrid zones.

Population genetic structure of sharpbelly Hemiculter leucisculus (Basilesky, 1855) and morphological diversification along climate gradients in China

Sharpbelly Hemiculter leucisculus (Basilewski, 1855) is a small, widespread, and native cyprinid fish with prominent habitat suitability and high invasive potential and is becoming the dominant species in freshwater ecosystems under intensified environmental disturbances. But how H. leucisculus acclimates to extremely heterogeneous environments remains unclear. In current study, the genetic structure of H. leucisculus was analyzed using Bayesian phylogenetic inference, haplotype network, and STRUCTURE base on cytb gene across 18 populations spanning 20 degrees of latitude and 18 degrees of longitude in China. The morphological diversification of body size and shape for H. leucisculus along the climate gradient was studied. The results showed that the 18 H. leucisculus populations were divided into 3 clusters: one cluster mainly from Huanghe River Basin, another cluster mainly from Yangzi River Basin, and H cluster containing Hainan and Beihai populations. The fish from southern populations were deeper bodied while individuals from northern populations were more slender. Inland individuals were more streamlined while coastal individuals were of deeper body. The partial Mantel test predicts that the potential mechanism underlining the intraspecies morphological diversification along climate gradients is primarily the divergent selection pressures among different environments, while genetic variation had less contribution to morphological differentiation. The formation of the Nanling Mountain Range could drive genetic differentiation between Beihai population and those from Yangzi River Basin. The present results highlight strong selective pressures of climate on widespread species and enrich morphological differentiation basis of acclimation for species with high habitat suitability and invasive potential.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

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.

Population genomics of parallel evolution in gene expression and gene sequence during ecological adaptation

Natural selection often produces parallel phenotypic changes in response to a similar adaptive challenge. However, the extent to which parallel gene expression differences and genomic divergence underlie parallel phenotypic traits and whether they are decoupled or not remains largely unexplored. We performed a population genomic study of parallel ecological adaptation among replicate ecotype pairs of the rough periwinkle (Littorina saxatilis) at a regional geographical scale (NW Spain). We show that genomic changes underlying parallel phenotypic divergence followed a complex pattern of both repeatable differences and of differences unique to specific ecotype pairs, in which parallel changes in expression or sequence are restricted to a limited set of genes. Yet, the majority of divergent genes were divergent either for gene expression or coding sequence, but not for both simultaneously. Overall, our findings suggest that divergent selection significantly contributed to the process of parallel molecular differentiation among ecotype pairs, and that changes in expression and gene sequence underlying phenotypic divergence could, at least to a certain extent, be considered decoupled processes.

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave-adapted and crab-adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double-digested restriction-associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%-15.1% of their divergent loci with a third more-distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid-shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.

Phylogeographic analysis of the East Asian goldenrod (Solidago virgaurea complex, Asteraceae) reveals hidden ecological diversification with recurrent formation of ecotypes

BACKGROUND AND AIMS

The processes and mechanisms underlying lineage diversification are major topics in evolutionary biology. Eurasian goldenrod species of the Solidago virgaurea complex show remarkable morphological and ecological diversity in the Japanese Archipelago, with ecotypic taxa well adapted to specific environments (climate, edaphic conditions and disturbance regimes). The species complex is a suitable model to investigate the evolutionary processes of actively speciating plant groups, due to its ability to evolve in relation to environmental adaptation and its historical population dynamics.

METHODS

Two chloroplast markers, 18 nuclear microsatellite markers and ddRAD-sequencing were used to infer population genetic demography of S. virgaurea complex with its related species/genera.

KEY RESULTS

Our analysis showed that populations in Japan form an evolutionary unit, which was genetically diverged from adjacent continental populations. The phylogenetic structure within the archipelago strongly corresponds to the geography, but interestingly there is no concordance between genetic structure and ecotypic boundaries; neighbouring populations of distinct ecotypes share a genetic background.

CONCLUSIONS

We propose that the traits specific to the ecotypic entities are maintained by natural selection or are very recently generated and have little effect on the genomes, making genome-wide genetic markers unsuitable for detecting ecotypic differentiation. Furthermore, some sporadically distributed taxa (found as rheophytes and alpine plants) were repeatedly generated from a more widespread taxon in geographically distant areas by means of selection. Overall, this study showed that the goldenrod complex has a high ability to evolve, enabling rapid ecological diversification over a recent timeframe.

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.

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow - a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co-occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment-length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host-associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.

Size selection by a gape-limited predator of a marine snail: Insights into magic traits for speciation

The intertidal snail Littorina saxatilis has repeatedly evolved two parallel ecotypes assumed to be wave adapted and predatory shore crab adapted, but the magnitude and targets of predator‐driven selection are unknown. In Spain, a small, wave ecotype with a large aperture from the lower shore and a large, thick‐shelled crab ecotype from the upper shore meet in the mid‐shore and show partial size‐assortative mating. We performed complementary field tethering and laboratory predation experiments; the first set compared the survival of two different size‐classes of the crab ecotype while the second compared the same size‐class of the two ecotypes. In the first set, the large size‐class of the crab ecotype survived significantly better than the small size‐class both on the upper shore and in the laboratory. In the second set, the small size‐class of the crab ecotype survived substantially better than that of the wave ecotype both on the upper shore and in the laboratory. Shell‐breaking predation on tethered snails was almost absent within the lower shore. In the laboratory shore crabs (Pachygrapsus marmoratus) with larger claw heights selected most strongly against the small size‐class of the crab ecotype, whereas those with medium claw heights selected most strongly against the thin‐shelled wave ecotype. Sexual maturity occurred at a much larger size in the crab ecotype than in the wave ecotype. Our results showed that selection on the upper shore for rapid attainment of a size refuge from this gape‐limited predator favors large size, thick shells, and late maturity. Model parameterization showed that size‐selective predation restricted to the upper shore resulted in the evolution of the crab ecotype despite gene flow from the wave ecotype snails living on the lower shore. These results on gape‐limited predation and previous ones showing size‐assortative mating between ecotypes suggest that size may represent a magic trait for the thick‐shelled ecotype.

Historical Biogeography of the Marine Snail Littorina saxatilis Inferred from Haplotype and Shell Morphology Evolution in NW Spain

The marine snail Littorina saxatilis exhibits extreme morphological variation between and within geographical regions and represents an excellent model for assessing local adaptation. Previous studies support the hypothesis of parallel evolution in sympatry of two morphologically different ecotypes (named as RB and SU) that co-inhabit different habitats from Galician rocky shores (NW Spain), and which are interrupted by sheltered areas inhabited by a different morph never studied before (named as SRB). Here, we use morphological and mitochondrial DNA (mtDNA) sequence data to test hypotheses on the origin and diversification of SRB snails and to assess their evolutionary relationships with RB and SU ecotypes. Our results show that the SRB morph displays the largest size and shell elongation and the smallest relative shell aperture, representing an extreme type of the RB vs. SU polymorphism, which has been linked to adaptation to sheltered ecological factors. Phylogenetic analysis shows that the SRB morph shares ancestry with RB and SU ecotypes, rejecting the hypothesis that the SRB morph marks relict populations from which these ecotypes evolved in Galician coasts. Our data support that genetic differentiation among SRB, RB and SU morphs results from a general pattern of restricted gene flow and isolation by distance linked to the colonization of Galician coasts by two independent mtDNA lineages, rather than from a random fragmentation of the initial distributional range. Therefore, the confinement of distinct lineages to specific geographical areas denote evident limits to the distances these snails can disperse. Morphological analysis indicates no association between mtDNA lineage and a specific morphotype, and suggests the independent gain of convergent morphological patterns within each mtDNA lineage in populations occupying contrasting habitats following the colonization of Galician coasts.

Ecological Specialization of Two Photobiont-Specific Maritime Cyanolichen Species of the Genus Lichina

All fungi in the class Lichinomycetes are lichen-forming and exclusively associate with cyanobacteria. Two closely related maritime species of the genus Lichina (L. confinis and L. pygmaea) show similar distribution ranges in the Northeast Atlantic, commonly co-occurring at the same rocky shores but occupying different littoral zones. By means of 16S rRNA and phycocyanin operon markers we studied a) the phylogenetic relationships of cyanobionts associated with these species, b) the match of divergence times between both symbionts, and c) whether Lichina species differ in photobiont association and in how geography and ecology affect selectivity. The cyanobionts studied are closely related to both marine and freshwater strains of the genus Rivularia. We found evidence of a high specificity to particular cyanobiont lineages in both species: Lichina pygmaea and L. confinis incorporate specific lineages of Rivularia that do not overlap at the haplotype nor the OTU levels. Dating divergences of the fungal and cyanobacterial partners revealed an asynchronous origin of both lineages. Within each fungal species, selectivity varied across the studied area, influenced by environmental conditions (both atmospheric and marine), although patterns were highly correlated between both lichen taxa. Ecological speciation due to the differential association of photobionts to each littoral zone is suspected to have occurred in marine Lichina.

What can be learnt from a snail?

The marine snail Littorina saxatilis is a common inhabitant of intertidal shores of the north Atlantic. It is amazingly polymorphic and forms reproductively isolated ecotypes in microhabitats where crabs are either present and wave action is less furious, or where waves are strong and crabs are absent. Decades of research have unveiled much of the ecological and demographic context of the formation of crab‐ and wave‐ecotype snails showing important phenotypic differences being inherited, differential selection being strong over adjacent microhabitats, local dispersal being restricted, and long‐distance transports of individuals being rare. In addition, strong assortative mating of ecotypes has been shown to include a component of male mate preference based on female size. Several studies support ecotypes being diverged locally and under gene flow in a parallel and highly replicated fashion. The high level of replication at various levels of independence (from local to pan‐European scale) provides excellent opportunities to investigate the detailed mechanisms of microevolution, including the formation of barriers to gene flow. Current investigations benefit from a draft reference genome and an integration of genomic approaches, modelling and experiments to unveil molecular and ecological components of speciation and their interactions.

Do the same genes underlie parallel phenotypic divergence in different Littorina saxatilis populations?

Parallel patterns of adaptive divergence and speciation are cited as powerful evidence for the role of selection driving these processes. However, it is often not clear whether parallel phenotypic divergence is underlain by parallel genetic changes. Here, we asked about the genetic basis of parallel divergence in the marine snail Littorina saxatilis, which has repeatedly evolved coexisting ecotypes adapted to either crab predation or wave action. We sequenced the transcriptome of snails of both ecotypes from three distant geographical locations (Spain, Sweden and United Kingdom) and mapped the reads to the L. saxatilis reference genome. We identified genomic regions potentially under divergent selection between ecotypes within each country, using an outlier approach based on F(ST) values calculated per locus. In line with previous studies indicating that gene reuse is generally common, we expected to find extensive sharing of outlier loci due to recent shared ancestry and gene flow between at least two of the locations in our study system. Contrary to our expectations, we found that most outliers were country specific, suggesting that much of the genetic basis of divergence is not shared among locations. However, we did find that more outliers were shared than expected by chance and that differentiation of shared outliers is often generated by the same SNPs. We discuss two mechanisms potentially explaining the limited amount of sharing we observed. First, a polygenic basis of divergent traits might allow for multiple distinct molecular mechanisms generating the same phenotypic patterns. Second, additional, location-specific axes of selection that we did not focus on in this study may produce distinct patterns of genetic divergence within each site.

Species and gene divergence in Littorina snails detected by array comparative genomic hybridization

Background

Array comparative genomic hybridization (aCGH) is commonly used to screen different types of genetic variation in humans and model species. Here, we performed aCGH using an oligonucleotide gene-expression array for a non-model species, the intertidal snail Littorina saxatilis. First, we tested what types of genetic variation can be detected by this method using direct re-sequencing and comparison to the Littorina genome draft. Secondly, we performed a genome-wide comparison of four closely related Littorina species: L. fabalis, L. compressa, L. arcana and L. saxatilis and of populations of L. saxatilis found in Spain, Britain and Sweden. Finally, we tested whether we could identify genetic variation underlying “Crab” and “Wave” ecotypes of L. saxatilis.

Results

We could reliably detect copy number variations, deletions and high sequence divergence (i.e. above 3%), but not single nucleotide polymorphisms. The overall hybridization pattern and number of significantly diverged genes were in close agreement with earlier phylogenetic reconstructions based on single genes. The trichotomy of L. arcana, L. compressa and L. saxatilis could not be resolved and we argue that these divergence events have occurred recently and very close in time. We found evidence for high levels of segmental duplication in the Littorina genome (10% of the transcripts represented on the array and up to 23% of the analyzed genomic fragments); duplicated genes and regions were mostly the same in all analyzed species. Finally, this method discriminated geographically distant populations of L. saxatilis, but we did not detect any significant genome divergence associated with ecotypes of L. saxatilis.

Conclusions

The present study provides new information on the sensitivity and the potential use of oligonucleotide arrays for genotyping of non-model organisms. Applying this method to Littorina species yields insights into genome evolution following the recent species radiation and supports earlier single-gene based phylogenies. Genetic differentiation of L. saxatilis ecotypes was not detected in this study, despite pronounced innate phenotypic differences. The reason may be that these differences are due to single-nucleotide polymorphisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-687) contains supplementary material, which is available to authorized users.

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.

Species tree of a recent radiation: the subfamily Delphininae (Cetacea, Mammalia)

Lineages undergoing rapid radiations provide exceptional opportunities for studying speciation and adaptation, but also represent a challenge for molecular systematics because retention of ancestral polymorphisms and the occurrence of hybridization can obscure relationships among lineages. Dolphins in the subfamily Delphininae are one such case. Non-monophyly, rapid speciation events, and discordance between morphological and molecular characters have made the inference of phylogenetic relationships within this subfamily very difficult. Here we approach this problem by applying multiple methods intended to estimate species trees using a multi-gene dataset for the Delphininae (Sousa, Sotalia, Stenella, Tursiops, Delphinus and Lagenodelphis). Incongruent gene trees obtained indicate that incomplete lineage sorting and possibly hybridization are confounding the inference of species history in this group. Nonetheless, using coalescent-based methods, we have been able to extract an underlying species-tree signal from divergent histories of independent genes. This is the first time a molecular study provides support for such relationships. This study further illustrates how methods of species-tree inference can be very sensitive both to the characteristics of the dataset and the evolutionary processes affecting the evolution of the group under study.

Explosive radiation of the genus Schizopera on a small subterranean island in Western Australia (Copepoda:Harpacticoida): unravelling the cases of cryptic speciation, size differentiation and multiple invasions

A previously unsurveyed calcrete aquifer in the Yilgarn region of Western Australia revealed an unprecedented diversity of copepods, representing 67% of that previously recorded in this whole region. Especially diverse was the genus Schizopera, with up to four morphospecies per bore and a significant size difference between them. Aims of this study were to: (1) survey the extent of this diversity using morphological and molecular tools; (2) derive a molecular phylogeny based on COI; and (3) investigate whether high diversity is a result of an explosive radiation, repeated colonisations, or both, size differentiation is a result of parallel evolution or different phylogeny, and whether Schizopera is a recent invasion in inland waters. More than 300 samples were analysed and the COI fragment successfully amplified by PCR from 43 specimens. Seven species and one subspecies are described as new, and three possible cryptic species were detected. Reconstructed phylogenies reveal that both explosive radiation and multiple colonisations are responsible for this richness, and that Schizopera is probably a recent invasion in these habitats. No evidence for parallel evolution was found, interspecific size differentiation being a result of different phylogeny. Sister species have parapatric distributions and show niche partitioning in the area of overlap.

Parallel Ecological Speciation in Plants?

Populations that have independently evolved reproductive isolation from their ancestors while remaining reproductively cohesive have undergone parallel speciation. A specific type of parallel speciation, known as parallel ecological speciation, is one of several forms of evidence for ecology's role in speciation. In this paper we search the literature for candidate examples of parallel ecological speciation in plants. We use four explicit criteria (independence, isolation, compatibility, and selection) to judge the strength of evidence for each potential case. We find that evidence for parallel ecological speciation in plants is unexpectedly scarce, especially relative to the many well-characterized systems in animals. This does not imply that ecological speciation is uncommon in plants. It only implies that evidence from parallel ecological speciation is rare. Potential explanations for the lack of convincing examples include a lack of rigorous testing and the possibility that plants are less prone to parallel ecological speciation than animals.

A METHODOLOGY FOR THE STOCHASTIC MODELING AND SIMULATION OF SYMPATRIC SPECIATION BY SEXUAL SELECTION

In the evolution literature, sympatric speciation is the origin of two, or more, species from a single local population. Many models have been developed to study the role of ecological competition and sexual selection in sympatric speciation.

In this paper we propose a methodology for systematically deriving efficient computational models to study speciation in populations evolving with overlapping generations. As a particular case, we consider sympatric speciation by sexual selection and we follow an individual based approach: a population is represented as a set of individuals that can mate and survive according to given probabilities.

We use our methodology to construct four different models for sympatric speciation, based on male traits and female preferences. These models differ in the genotypical representation of the individuals. Results of simulations in the different models are shown and discussed.

The study of the models show that sympatric speciation by sexual selection is unlikely, also with a favorable distribution of genotypes in the initial population.

Mechanical reproductive isolation facilitates parallel speciation in western North American scincid lizards

Mechanical reproductive barriers have been dismissed as a major driver of animal speciation, yet the extent to which such barriers cause reproductive isolation in most animal groups is largely unknown and rarely tested. In this study, we used hierarchical Bayesian modeling of mate compatibility experiments to show that body size divergence in lizards of the Plestiodon skiltonianus complex contributes to reproductive isolation in at least three ways: males preferably court females that are more similar in size, females reject males that are highly divergent in size, and the size difference of a male and female in copula constrains the ability to align the genitalia for intromission. We used a predictive model to estimate the contributions of behavioral and mechanical barriers to reproductive isolation between populations with differing degrees of size divergence. This model shows that the mechanical barrier is more important than behavioral barriers at small and intermediate degrees of size divergence, suggesting that it acts earlier during speciation when body morphology is more similar between diverging lineages. Given that correlated divergence in size and ecology is common in animals, similar constraints imposed by the geometry of the mating posture may apply to a variety of major animal lineages and merit further attention in speciation research.

The Littorina sequence database (LSD)--an online resource for genomic data

We present an interactive, searchable expressed sequence tag database for the periwinkle snail Littorina saxatilis, an upcoming model species in evolutionary biology. The database is the result of a hybrid assembly between Sanger and 454 sequences, 1290 and 147,491 sequences respectively. Normalized and non-normalized cDNA was obtained from different ecotypes of L. saxatilis collected in the UK and Sweden. The Littorina sequence database (LSD) contains 26,537 different contigs, of which 2453 showed similarity with annotated proteins in UniProt. Querying the LSD permits the selection of the taxonomic origin of blast hits for each contig, and the search can be restricted to particular taxonomic groups. The database allows access to UniProt annotations, blast output, protein family domains (PFAM) and Gene Ontology. The database will allow users to search for genetic markers and identifying candidate genes or genes for expression analyses. It is open for additional deposition of sequence information for L. saxatilis and other species of the genus Littorina. The LSD is available at http://mbio-serv2.mbioekol.lu.se/Littorina/.

Adaptive traits are maintained on steep selective gradients despite gene flow and hybridization in the intertidal zone

Gene flow among hybridizing species with incomplete reproductive barriers blurs species boundaries, while selection under heterogeneous local ecological conditions or along strong gradients may counteract this tendency. Congeneric, externally-fertilizing fucoid brown algae occur as distinct morphotypes along intertidal exposure gradients despite gene flow. Combining analyses of genetic and phenotypic traits, we investigate the potential for physiological resilience to emersion stressors to act as an isolating mechanism in the face of gene flow. Along vertical exposure gradients in the intertidal zone of Northern Portugal and Northwest France, the mid-low shore species Fucus vesiculosus, the upper shore species Fucus spiralis, and an intermediate distinctive morphotype of F. spiralis var. platycarpus were morphologically characterized. Two diagnostic microsatellite loci recovered 3 genetic clusters consistent with prior morphological assignment. Phylogenetic analysis based on single nucleotide polymorphisms in 14 protein coding regions unambiguously resolved 3 clades; sympatric F. vesiculosus, F. spiralis, and the allopatric (in southern Iberia) population of F. spiralis var. platycarpus. In contrast, the sympatric F. spiralis var. platycarpus (from Northern Portugal) was distributed across the 3 clades, strongly suggesting hybridization/introgression with both other entities. Common garden experiments showed that physiological resilience following exposure to desiccation/heat stress differed significantly between the 3 sympatric genetic taxa; consistent with their respective vertical distribution on steep environmental clines in exposure time. Phylogenetic analyses indicate that F. spiralis var. platycarpus is a distinct entity in allopatry, but that extensive gene flow occurs with both higher and lower shore species in sympatry. Experimental results suggest that strong selection on physiological traits across steep intertidal exposure gradients acts to maintain the 3 distinct genetic and morphological taxa within their preferred vertical distribution ranges. On the strength of distributional, genetic, physiological and morphological differences, we propose elevation of F. spiralis var. platycarpus from variety to species level, as F. guiryi.

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F(ST) outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F(ST) partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F(ST) differed between candidate and noncandidate loci for all partitions except in the case of micro-geographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.

Historical and contemporary population genetic connectivity of the European short-snouted seahorse Hippocampus hippocampus and implications for management

This first genetic study of Hippocampus hippocampus covers the species' entire geographic range and employs two mtDNA markers (control region and cytochrome b) to establish patterns of population structuring. A total of 255 specimens from 21 locations were used to obtain 89 concatenated haplotypes. The common haplotype was present in all but one population, however, most haplotypes were unique. The haplotype network had a star-like construction, suggesting expansion from a bottleneck event. F(ST) and AMOVA revealed population subdivision into three geographic regions (English Channel + Bay of Biscay, Mediterranean Sea + Atlantic Ocean Iberian coast + Macaronesian Islands, and West Africa) with barriers to gene flow indentified at Cape Finisterre and the Cape Verde frontal zone. Neutrality tests and nested clade analysis suggest a complex demographic history, with both historic events and contemporary processes shaping patterns of genetic differentiation. The genetic population subdivision detected in this study indicates that H. hippocampus should be managed as three separate units. This is especially pertinent as H. hippocampus populations within the West African region are the only ones known to be specifically targeted for exploitation.

Phylogeographic analysis reveals a deep lineage split within North Atlantic Littorina saxatilis

Phylogeographic studies provide critical insight into the evolutionary histories of model organisms; yet, to date, range-wide data are lacking for the rough periwinkle Littorina saxatilis, a classic example of marine sympatric speciation. Here, we use mitochondrial DNA (mtDNA) sequence data to demonstrate that L. saxatilis is not monophyletic for this marker, but is composed of two distinct mtDNA lineages (I and II) that are shared with sister species Littorina arcana and Littorina compressa. Bayesian coalescent dating and phylogeographic patterns indicate that both L. saxatilis lineages originated in the eastern North Atlantic, around the British Isles, at approximately 0.64 Ma. Both lineages are now distributed broadly across the eastern, central and western North Atlantic, and show strong phylogeographic structure among regions. The Iberian Peninsula is genetically distinct, suggesting prolonged isolation from northeastern North Atlantic populations. Western North Atlantic populations of L. saxatilis lineages I and II predate the last glacial maximum and have been isolated from eastern North Atlantic populations since that time. This identification of two distinct, broadly distributed mtDNA lineages further complicates observed patterns of repeated incipient ecological speciation in L. saxatilis, because the sympatric origins of distinct ecotype pairs on eastern North Atlantic shores may be confounded by admixture of divergent lineages.

"Same same but different": replicated ecological speciation at White Sands

Understanding the factors that promote or inhibit species formation remains a central focus in evolutionary biology. It has been difficult to make generalities about the process of ecological speciation in particular given that each example is somewhat idiosyncratic. Here we use a case study of replicated ecological speciation in the same selective environment to assess factors that account for similarities and differences across taxa in progress towards ecological speciation. We study three different species of lizards on the gypsum sand dunes of White Sands, New Mexico, and present evidence that all three fulfill the essential factors for ecological speciation. We use multilocus nuclear data to show that progress toward ecological speciation is unequal across the three species. We also use morphometric data to show that traits other than color are likely under selection and that selection at White Sands is both strong and multifarious. Finally, we implicate geographic context to explain difference in progress toward speciation in the three species. We suggest that evaluating cases from the natural world that are "same same but different" can reveal the mechanisms of ecological speciation.

Insights into the role of differential gene expression on the ecological adaptation of the snail Littorina saxatilis

Background

In the past 40 years, there has been increasing acceptance that variation in levels of gene expression represents a major source of evolutionary novelty. Gene expression divergence is therefore likely to be involved in the emergence of incipient species, namely, in a context of adaptive radiation. In this study, a genome-wide expression profiling approach (cDNA-AFLP), validated by quantitative real-time polymerase chain reaction (qPCR) were used to get insights into the role of differential gene expression on the ecological adaptation of the marine snail Littorina saxatilis. This gastropod displays two sympatric ecotypes (RB and SU) which are becoming one of the best studied systems for ecological speciation.

Results

Among the 99 transcripts shared between ecotypes, 12.12% showed significant differential expression. At least 4% of these transcripts still displayed significant differences after correction for multiple tests, highlighting that gene expression can differ considerably between subpopulations adapted to alternative habitats in the face of gene flow. One of the transcripts identified was Cytochrome c Oxidase subunit I (COI). In addition, 6 possible reference genes were validated to normalize and confirm this result using qPCR. α-Tubulin and histone H3.3 showed the more stable expression levels, being therefore chosen as the best option for normalization. The qPCR analysis confirmed a higher COI expression in SU individuals.

Conclusions

At least 4% of the transcriptome studied is being differentially expressed between ecotypes living in alternative habitats, even when gene flow is still substantial between ecotypes. We could identify a candidate transcript of such ecotype differentiation: Cytochrome c Oxidase Subunit I (COI), a mitochondrial gene involved in energy metabolism. Quantitative PCR was used to confirm the differences found in COI and its over-expression in the SU ecotype. Interestingly, COI is involved in the oxidative phosphorylation, suggesting an enhanced mitochondrial gene expression (or increased number of mitochondria) to improve energy supply in the ecotype subjected to the strongest wave action.

Intron sequences of arginine kinase in an intertidal snail suggest an ecotype-specific selective sweep and a gene duplication

Many species with restricted gene flow repeatedly respond similarly to local selection pressures. To fully understand the genetic mechanisms behind this process, the phylogeographic history of the species (inferred from neutral markers) as well as the loci under selection need to be known. Here we sequenced an intron in the arginine kinase gene (Ark), which shows strong clinal variation between two locally adapted ecotypes of the flat periwinkle, Littorina fabalis. The 'small-sheltered' ecotype was almost fixed for one haplotype, H1, in populations on both sides of the North Sea, unlike the 'large-moderately exposed ecotype', which segregated for ten different haplotypes. This contrasts with neutral markers, where the two ecotypes are equally variable. H1 could have been driven to high frequency in an ancestral population and then repeatedly spread to sheltered habitats due to local selection pressures with the colonization of both sides of the North Sea, after the last glacial maximum (~18 000 years ago). An alternative explanation is that a positively selected mutation, in or linked to Ark, arose after the range expansion and secondarily spread through sheltered populations throughout the distribution range, causing this ecotype to evolve in a concerted fashion. Also, we were able to sequence up to four haplotypes consistently from some individuals, suggesting a gene duplication in Ark.

Repeated evolution of reproductive isolation in a marine snail: unveiling mechanisms of speciation

Distinct ecotypes of the snail Littorina saxatilis, each linked to a specific shore microhabitat, form a mosaic-like pattern with narrow hybrid zones in between, over which gene flow is 10-30% of within-ecotype gene flow. Multi-locus comparisons cluster populations by geographic affinity independent of ecotype, while loci under selection group populations by ecotype. The repeated occurrence of partially reproductively isolated ecotypes and the conflicting patterns in neutral and selected genes can either be explained by separation in allopatry followed by secondary overlap and extensive introgression that homogenizes neutral differences evolved under allopatry, or by repeated evolution in parapatry, or in sympatry, with the same ecotypes appearing in each local site. Data from Spain, the UK and Sweden give stronger support for a non-allopatric model of ecotype formation than for an allopatric model. Several different non-allopatric mechanisms can, however, explain the repeated evolution of the ecotypes: (i) parallel evolution by new mutations in different populations; (ii) evolution from standing genetic variation; and (iii) evolution in concert with rapid spread of new positive mutations among populations inhabiting similar environments. These models make different predictions that can be tested using comprehensive phylogenetic information combined with candidate loci sequencing.

It's about time: divergence, demography, and the evolution of developmental modes in marine invertebrates

Differences in larval developmental mode are predicted to affect ecological and evolutionary processes ranging from gene flow and population bottlenecks to rates of population recovery from anthropogenic disturbance and capacity for local adaptation. The most powerful tests of these predictions use comparisons among species to ask how phylogeographic patterns are correlated with the evolution and loss of prolonged planktonic larval development. An important and largely untested assumption of these studies is that interspecific differences in population genetic structure are mainly caused by differences in dispersal and gene flow (rather than by differences in divergence times among populations or changes in effective population sizes), and that species with similar patterns of spatial genetic variation have similar underlying temporal demographic histories. Teasing apart these temporal and spatial patterns is important for understanding the causes and consequences of evolutionary changes in larval developmental mode. New analytical methods that use the coalescent history of allelic diversity can reveal these temporal patterns, test the strength of traditional population-genetic explanations for variation in spatial structure based on differences in dispersal, and identify strongly supported alternative explanations for spatial structure based on demographic history rather than on gene flow alone. We briefly review some of these recent analytical developments, and show their potential for refining ideas about the correspondence between the evolution of larval developmental mode, population demographic history, and spatial genetic variation.