Proteomic evidence of a paedomorphic evolutionary process within a marine snail species: a strategy for adapting to extreme ecological conditions?

The exposed and sheltered ecotypes of the marine snail Littorina saxatilis from European rocky shores are considered a key model system to study adaptation and ecological speciation. Previous studies showed that two ecotypes (RB and SU) of this species in NW Spain have adapted differently to different shore levels and microhabitats. In order to understand how this divergent adaptive process has been accomplished, we followed a quantitative proteomic approach to investigate the proteome variation in a number of different biological factors, that is, ecotype, ontogeny and their interactions. This approach allowed testing the hypothesis that one of the ecotypes has evolved by paedomorphosis, and also whether or not the molecular mechanisms related to ecotype differentiation are set up in early developmental stages. Additionally, the identification of some candidate proteins using mass spectrometry provides some functional insights into these evolutionary processes. Results from this study provided evidence of higher ontogenetic differentiation at proteome level in the RB (metamorphic) than in SU (paedomorphic) ecotype that point to the possibility of juvenile stage retention in this latter ecotype. The level of protein expression (proteome) differences between ecotypes maintained nearly constant from late embryonic stages to adulthood, although some proteins involved in these changes considerably differed in embryonic compared to other ontogenetic stages. Paedomorphosis may be the evolutionary response of the SU ecotype of solving the trade‐off during sexually immaturity that is caused by the evolution of small size arising from adaptation to the wave‐exposed habitat. Some potential candidate genes of adaptation related to energetic metabolism have been identified, providing a promising baseline for future functional analyses.     

Genetic differentiation and estimation of effective population size and migration rates in two sympatric ecotypes of the marine snail Littorina saxatilis

On exposed rocky shores in Galicia (northwest Spain), a striking polymorphism exists between two ecotypes (RB and SU) of Littorina saxatilis that occupy different levels of the intertidal zone and exhibit an incomplete reproductive isolation. The setting has been suggested to represent ongoing sympatric speciation by ecological adaptation of the two ecotypes to their respective habitats. In this article we address whether or not the ecotypes have developed their own population structures in response to the rigors of their corresponding environments and life histories. We analyzed four to five allozymic loci from three surveys of the same sites, spanning a 14-year period. An experimental design including three localities with two transects per locality and three shore levels allowed studying temporal and spatial population structure and estimation of effective population sizes (N(e)), neighborhood sizes (N(n)), and migration rates (m). Genetic differentiation was significantly lower in RB populations (theta(ST) = 0.067) than in SU ones (theta(ST) = 0.124). Mean estimates of N(e), N(n), and m did not differ significantly between ecotypes, but local ecotype differences in migration between the two closest localities (larger migration rates in RB than in SU populations) could explain the pattern in population differentiation.     

Nonallopatric and parallel origin of local reproductive barriers between two snail ecotypes

Theory suggests that speciation is possible without physical isolation of populations (hereafter, nonallopatric speciation), but recent nonallopatric models need the support of irrefutable empirical examples. We collected snails (Littorina saxatilis) from three areas on the NW coast of Spain to investigate the population genetic structure of two ecotypes. Earlier studies suggest that these ecotypes may represent incipient species: a large, thick-shelled 'RB' ecotype living among the barnacles in the upper intertidal zone and a small, thin-shelled 'SU' ecotype living among the mussels in the lower intertidal zone only 10-30 m away. The two ecotypes overlap and hybridize in a midshore zone only 1-3 m wide. Three different types of molecular markers [allozymes, mitochondrial DNA (mtDNA) and microsatellites] consistently indicated partial reproductive isolation between the RB and the SU ecotypes at a particular site. However, each ecotype was related more closely to the other ecotype from the same site than to the same ecotype from another site further along the Galician coast (25-77 km away). These findings supported earlier results based solely on allozyme variation and we could now reject the possibility that selection produced these patterns. The patterns of genetic variation supported a nonallopatric model in which the ecotypes are formed independently at each site by parallel evolution and where the reproductive barriers are a byproduct of divergent selection for body size. We argue that neither our laboratory hybridization experiments nor our molecular data are compatible with a model based on allopatric ecotype formation, secondary overlap and introgression.     

Incipient post‐zygotic barrier in a model system of ecological speciation with gene flow

The role of post‐zygotic isolation in nonallopatric ecological speciation is still mostly unknown and information on the nature and strength of these barriers in well‐known speciation models is essential for a deeper understanding of such processes. The Galician ecotypes of the marine snail Littorina saxatilis represent one of the best studied cases of nonallopatric ecological speciation. Here, we test the existence of incipient post‐zygotic isolation by comparing the fertility of male hybrids with that of both pure forms [ridged and banded (RB) and smooth and unbanded (SU) ecotypes]. We analysed the degree of sperm DNA fragmentation (SDF) of individuals morphologically classified as RB, SU and hybrids, sampled from two locations. SDF analyses were chosen to study sperm quality because, in other animal species, SDF rates correlate with important parameters for speciation research, such as fertilization and abortion rates and viability of adult progeny. In the present work, hybrids showed significantly higher SDF rates than RB and SU males in one location and significantly higher variances in both locations. These results suggest the existence of an incipient post‐zygotic barrier, the strength of which may vary across the Galician shore, and highlight the potential of SDF analyses for speciation research.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Case studies and mathematical models of ecological speciation. 3: Ecotype formation in a Swedish snail

Formation of partially reproductively isolated ecotypes in the rough periwinkle, Littorina saxatilis , may be a case of incipient nonallopatric ecological speciation. To better understand the dynamics of ecotype formation, its timescale, driving forces and evolutionary consequences, we developed a spatially explicit, individual-based model incorporating relevant ecological, spatial and mate selection data for Swedish L. saxatilis . We explore the impact of bounded hybrid superiority, ecological scenarios and mate selection systems on ecotype formation, gene flow and the evolution of prezygotic isolation. Our model shows that ecotypes are expected to form rapidly in parapatry under conditions applicable to Swedish L. saxatilis and may proceed to speciation. However, evolution of nonrandom mating had complex behaviour. Ecotype evolution was inhibited by pre-existing mating preferences, but facilitated by the evolution of novel preferences. While in many scenarios positive assortative mating reduced gene flow between ecotypes, in others negative assortative mating arose, preferences were lost after ecotype formation, preferences were confined to one ecotype or the ancestral ecotype became extinct through sexual selection. Bounded hybrid superiority (as observed in nature) enhanced ecotype formation but increased gene flow. Our results highlight that ecotype formation and speciation are distinct processes: factors that contribute to ecotype formation can be detrimental to speciation and vice versa. The complex interactions observed between local adaptation and nonrandom mating imply that generalization from data is unreliable without quantitative theory for speciation.     

Lack of early laboratory postzygotic reproductive isolation between two ecotypes of Littorina saxatilis (Mollusca, Gastropoda) showing strong premating sexual isolation

Two sympatric and divergent adaptive ecotypes of Littorina saxatilis (RB and SU) are known to hybridize showing partial premating isolation in the wild. Previous studies have revealed that morphological intermediate forms (presumably hybrids) present fitness (viability, sexual selection and fecundity) similar to that from pure ecotypes at the mid-shore. However, the absence of postzygotic isolation due to genetic incompatibility cannot be ruled out unless it is measured directly on true F ₁ hybrids. In this study, we overcome this problem and present data on 56 individual crosses including the four possible mating combinations (RB/RB, RB/SU, SU/RB and SU/SU) to compare fertilization and fecundity rates (including young progeny viability) between the four type crosses. Pooled RB female crosses showed apparently larger fertility and fecundity than pooled SU female crosses, probably because of differences in fecundity and laboratory survivorship between ecotypes. However, similar fertilization and fecundity rates were found for both RB and SU females when mated with different male types, supporting the idea that genetic-incompatibility-based postzygotic isolation can be ignored as a major determinant of this polymorphism in nature.     

Male discrimination of female mucous trails permits assortative mating in a marine snail species

Recent research has shown the potential for nonallopatric speciation, but we lack an adequate understanding of the mechanisms of prezygotic barriers and how these evolve in the presence of gene flow. The marine snail Littorina saxatilis has distinct ecotypes in different shore microhabitats. Ecotypes hybridize in contact zones, but gene flow is impeded by assortative mating. Earlier studies have shown that males and females of the same ecotype copulate for longer than mates of different ecotype. Here we report a new mechanism that further contributes to reproductive isolation between ecotypes in the presence of gene flow. This mechanism is linked to the ability of males to track potential partners by following their mucous trail. We show that cliff ecotype males follow the trails of females of the same ecotype for longer than females of the alternate (boulder) ecotype. In addition, cliff males are more likely to follow the mucous trail in the correct direction if the trail is laid by a cliff-female. The capacity to discriminate the ecotype of female mucous trails combined with differential copulation times creates a strong prezygotic reproductive barrier between ecotypes of L. saxatilis that reduces gene flow from cliff to boulder ecotypes by >/=80%.     

Habitat-dependent Ecotype Micro-distribution at the Mid-shore in Natural Populations of <Emphasis Type="Italic">Littorina saxatilis</Emphasis>

On the NW coast of Spain, there is a hybrid zone in natural populations of L. saxatilis, where the largest ecotype is adapted to the desiccation and heat stress of the upper shore, whereas the smaller one is adapted to the wave impacts of the lower shore. The two ecotypes meet and hybridise in a mid-shore area producing intermediate individuals or hybrids. It has been postulated that this hybrid zone is maintained by a selection-gradient model that assumes habitat-dependent selection for the different environments of upper and lower shore areas. In this study, we focus on the hybrid (mid-shore) area, where a transition occurs between upper-shore and lower-shore environments, in order to ascertain whether the habitat-dependent selection is maintained at the micro-habitat scale. We present data on snail density at several levels of the mid-shore (varying accordingly in the relative frequency of mussels and barnacles) obtained for three seasons in three consecutive years and three localities. In the mid-shore, the abundance of the RB ecotype increases with increasing tidal height, whereas the abundance of the SU ecotype increases with decreasing tidal height, suggesting that habitat-dependent selection is maintaining the ecotype micro-distribution.     

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.     

Proteomic analysis of Arabidopsis thaliana ecotypes with contrasted root architecture in response to phosphate deficiency

Owing to a weak availability in soil, plants have developed numerous morphological, physiological and biochemical adaptations to acquire phosphate (Pi). Identification and characterisation of key genes involved in the initial steps of Pi-signalling might provide clues about the regulation of the complex Pi deficiency adaptation mechanism. A two-dimensional gel electrophoresis approach was performed to investigate proteome responses to Pi starvation in Arabidopsis. Two ecotypes were selected according to contrasting responses of their root system architecture to low availability of Pi. Thirty protein spots were shown to be affected by Pi deficiency. Fourteen proteins appeared to be up-regulated and ten down-regulated with ecotype Be-0, wheras only thirteen proteins were observed as down-regulated for ecotype Ll-0. Furthermore, systematic and opposite responses to Pi deficiency were observed between the two ecotypes. The sequences of these 30 differentially expressed protein spots were identified using mass spectrometry, and most of the proteins were involved in oxidative stress, carbohydrate and proteins metabolism. The results suggested that the modulation of alcohol dehydrogenase, malic enzyme and aconitate hydratase may contribute to the contrasted adaptation strategy to Pi deficiency of Be-0 and Ll-0 ecotypes. A focus on aconitate hydratase highlighted a complex reverse response of the pattern of corresponding spots between the two ecotypes. This protein, also potentially involved in iron homeostasis, was speculated to contribute, at least indirectly, to the root architecture response of these ecotypes.     

Immigrant and extrinsic hybrid seed inviability contribute to reproductive isolation between forest and dune ecotypes of Epipactis helleborine (Orchidaceae)

Reproductive isolation caused by divergent natural selection arising from differences between ecological environments or ecological interactions represents a major mechanism contributing to speciation, but its relative importance is poorly known. In this study, controlled reciprocal crossings and seed germination experiments were combined with genetic and morphometric analyses to test the hypothesis that previously described differences in mycorrhizal communities between forest and dune ecotypes of Epipactis helleborine were sufficiently strong to create complete reproductive isolation between the two ecotypes. Molecular analyses using 770 SNP markers showed that the two ecotypes were genetically distinct and that populations of the dune ecotype were genetically impoverished compared to populations of the forest ecotype. Morphologically, the two ecotypes were also significantly different, with plants of the dune ecotype generally being smaller than plants of the forest ecotype. The results further showed that immigrant seeds had a significantly lower probability of protocorm formation than native seeds, indicating strong immigrant inviability. Although both ecotypes were able to cross easily and to produce a large number of viable seeds, hybrid seeds showed significantly lower protocorm formation than pure seeds, further contributing to reproductive isolation. Overall, these results indicate that interfertile populations of a widespread orchid adapting to contrasting environments diverge as a consequence of concurrent selection acting against immigrants and hybrids and suggest that mycorrhizal fungi can play a role in the early stages of plant speciation.     

Proteomic characterization of Phragmites communis in ecotypes of swamp and desert dune

Phragmites communis Trin. (common reed) is a recognized model plant for studying its adaptation to contrasting and harsh environments. To understand the inherent molecular basis for its remarkable resistance to combined stresses, we performed a comprehensive proteomic analysis of the leaf proteins from two ecotypes, i.e. swamp and desert dune, naturally growing in the desert region of northwestern China. First, a proteome reference map of Phragmites was established based on the swamp ecotype. Proteins were resolved by 2‐D/SDS‐PAGE and identified by MALDI‐TOF/TOF MS. In total, 177 spots were identified corresponding to 51 proteins. The major proteins identified are proteins involved in photosynthesis, glutathione and ascorbic acid metabolism as well as protein synthesis and quality control. Second, the 2‐DE profiles of the two ecotypes were compared quantitatively via DIGE analysis. Compared with swamp ecotype, 51 proteins spots are higher‐expressed and 58 protein spots are lower‐expressed by twofold or more in desert dune ecotype. Major differences were found for the proteins involved in light reaction of photosynthesis, protein biosynthesis and quality control and antioxidative reactions. The physiological significance of such differences is discussed in the context of a flow of complex events in relation to plant adaptation to combined environmental stresses.     

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.     

Differences in timing of migration and response to sexual signalling drive asymmetric hybridization across a migratory divide

Ecological traits and sexual signals may both contribute to the process of ecological speciation. Here we investigate the roles of an ecological trait, seasonal migratory behaviour and a sexual trait, song, in restricting or directing gene flow across a migratory divide in the Swainson's thrush (Catharus ustulatus). We show that short-distance migratory ecotypes wintering in Central America arrive earlier at the breeding grounds than long-distance migratory ecotypes wintering primarily in South America, providing the potential for some premating isolation. Playback experiments suggest that early- and late-arriving forms recognize each other as competitors, but that the early-arriving form responds more aggressively to a broader spectrum of stimuli. Genetic analysis suggests that hybridization occurs more often between males of the early-arriving ecotype and females of the late-arriving ecotype. Together our results suggest that differences in arrival times may reduce the temporal coincidence of mate choice, but asymmetry in response to heterotypic song may hinder complete divergence. These data provide further insight into the roles of ecological traits and sexual signals during the incipient stages of speciation.     

Recent stable insertion of mitochondrial DNA into an Arabidopsis polyubiquitin gene by nonhomologous recombination.

Sequence analysis of a newly identified polyubiquitin gene (UBQ13) from the Columbia ecotype of Arabidopsis thaliana revealed that the gene contained a 3.9-kb insertion in the coding region. All subclones of the 3.9-kb insert hybridized to isolated mitochondrial DNA. The insert was found to consist of at least two, possibly three, distinct DNA segments from the mitochondrial genome. A 590-bp region of the insert is nearly identical to the Arabidopsis mitochondrial nad1 gene. UBQ13 restriction fragments in total cellular DNA from ecotypes Ler, No-0, Be-0, WS, and RLD were identified and, with the exception of Be-0, their sizes were equivalent to that predicted from the corresponding ecotype Columbia UBQ13 restriction fragment without the mitochondrial insert. Isolation by polymerase chain reaction and sequence determination of UBQ13 sequences from the other ecotypes showed that all lacked the mitochondrial insert. All ecotypes examined, except Columbia, contain intact open reading frames in the region of the insert, including four ubiquitin codons which Columbia lacks. This indicates that the mitochondrial DNA in UBQ13 in ecotype Columbia is the result of an integration event that occurred after speciation of Arabidopsis rather than a deletion event that occurred in all ecotypes except Columbia. This stable movement of mitochondrial DNA to the nucleus is so recent that there are few nucleotide changes subsequent to the transfer event. This allows for precise analysis of the sequences involved and elucidation of the possible mechanism. The presence of intron sequences in the transferred nucleic acid indicates that DNA was the transfer intermediate. The lack of sequence identity between the integrating sequence and the target site, represented by the other Arabidopsis ecotypes, suggests that integration occurred via nonhomologus recombination. This nuclear/organellar gene transfer event is strikingly similar to the experimentally accessible process of nuclear integration of introduced heterologous DNA.     

Do floral and ecogeographic isolation allow the co‐occurrence of two ecotypes of Anacamptis papilionacea (Orchidaceae)?

Ecotypes are relatively frequent in flowering plants and considered central in ecological speciation as local adaptation can promote the insurgence of reproductive isolation. Without geographic isolation, gene flow usually homogenizes the allopatrically generated phenotypic and ecological divergences, unless other forms of reproductive isolation keep them separated. Here, we investigated two orchid ecotypes with marked phenotypic floral divergence that coexist in contact zones. We found that the two ecotypes show different ecological habitat preferences with one being more climatically restricted than the other. The ecotypes remain clearly morphologically differentiated both in allopatry and in sympatry and differed in diverse floral traits. Despite only slightly different flowering times, the two ecotypes achieved floral isolation thanks to different pollination strategies. We found that both ecotypes attract a wide range of insects, but the ratio of male/female attracted by the two ecotypes was significantly different, with one ecotype mainly attracts male pollinators, while the other mainly attracts female pollinators. As a potential consequence, the two ecotypes show different pollen transfer efficiency. Experimental plots with pollen staining showed a higher proportion of intra‐ than interecotype movements confirming floral isolation between ecotypes in sympatry while crossing experiments excluded evident postmating barriers. Even if not completely halting the interecotypes pollen flow in sympatry, such incipient switch in pollination strategy between ecotypes may represent a first step on the path toward evolution of sexual mimicry in Orchidinae.     

Genetic variation for shell traits in a direct-developing marine snail involved in a putative sympatric ecological speciation process

Populations of the marine gastropod Littorina saxatilis from exposed rocky shores of NW Spain provide one of the few putative cases of sympatric ecological speciation. Two ecotypes with large differences in shell morphology and strong assortative mating are living at different vertical levels of the shore separated by a few meters. It has been hypothesized that shell size is the main determinant for the reproductive isolation observed between the ecotypes, and that several shell shape traits are subject to divergent natural selection and are responsible for the adaptation of each ecotype to its respective habitat. Using embryos extracted from wild females we obtain estimates of genetic variation for shell size and shape and compare them with those from neutral molecular markers. Estimates of heritability are significantly larger for the ecotype found in the upper shore than for that in the lower shore, in concordance with a similar result observed for heterozygosity of neutral markers. The large genetic differentiation between ecotypes for the shell traits, contrasting the smaller close to neutral differentiation between populations of the same ecotype, supports the implication of the traits in adaptation.