J. Grahame

Differential gene exchange between parapatric morphs of Littorina saxatilis detected using AFLP markers

Speciation requires the acquisition of reproductive isolation, and the circumstances under which this could evolve are of great interest. Are new species formed after the acquisition of generalized incompatibility arising between physically separated populations, or may they arise as a result of the action of disruptive selection beginning with the divergence of a rather restricted set of gene loci? Here we apply the technique of amplified fragment length polymorphism (AFLP) analysis to an intertidal snail whose populations display a cline in shell shape across vertical gradients on rocky shores. We compare the FST values for 306 AFLP loci with the distribution of FST estimated from a simulation model using values of mutation and migration derived from the data. We find that about 5% of these loci show greater differentiation than expected, providing evidence of the effects of selection across the cline, either direct or indirect through linkage. This is consistent with expectations from nonallopatric speciation models that propose an initial divergence of a small part of the genome driven by strong disruptive selection while divergence at other loci is prevented by gene flow. However, the pattern could also be the result of differential introgression after secondary contact.

Sympatric, parapatric or allopatric: the most important way to classify speciation?

The most common classification of modes of speciation begins with the spatial context in which divergence occurs: sympatric, parapatric or allopatric. This classification is unsatisfactory because it divides a continuum into discrete categories, concentrating attention on the extremes, and it subordinates other dimensions on which speciation processes vary, such as the forces driving differentiation and the genetic basis of reproductive isolation. It also ignores the fact that speciation is a prolonged process that commonly has phases in different spatial contexts. We use the example of local adaptation and partial reproductive isolation in the intertidal gastropod Littorina saxatilis to illustrate the inadequacy of the spatial classification of speciation modes. Parallel divergence in shell form in response to similar environmental gradients in England, Spain and Sweden makes this an excellent model system. However, attempts to demonstrate ‘incipient’ and ‘sympatric’ speciation involve speculation about the future and the past. We suggest that it is more productive to study the current balance between local adaptation and gene flow, the interaction between components of reproductive isolation and the genetic basis of differentiation.

Adaptation to a steep environmental gradient and an associated barrier to gene exchange in Littorina saxatilis

Abstract Steep environmental gradients offer important opportunities to study the interaction between natural selection and gene flow. Allele frequency clines are expected to form at loci under selection, but unlinked neutral alleles may pass easily across these clines unless a generalized barrier evolves. Here we consider the distribution of forms of the intertidal gastropod Littorina saxatilis, analyzing shell shape and amplified fragment length polymorphism (AFLP) loci on two rocky shores in Britain. On the basis of previous work, the AFLP loci were divided into differentiated and undifferentiated groups. On both shores, we have shown a sharp cline in allele frequencies between the two morphs for differentiated AFLP loci. This is coincident with a habitat transition on the shore where the two habitats (cliff and boulder field) are immediately contiguous. The allele frequency clines coincide with a cline in shell morphology. In the middle of the cline, linkage disequilibrium for the differentiated loci rises in accordance with expectation. The clines are extremely narrow relative to dispersal, probably as a result of both strong selection and habitat choice. An increase in FST for undifferentiated AFLPs between morphs, relative to within‐morph comparisons, is consistent with there being a general barrier to gene flow across the contact zone. These features are consistent either with an episode of allopatric divergence followed by secondary contact or with primary, nonallopatric divergence. Further data will be needed to distinguish between these alternatives

Sequence differentiation in regions identified by a genome scan for local adaptation.

Genome scans using large numbers of randomly selected markers have revealed a small proportion of loci that deviate from neutral expectations and so may mark genomic regions that contribute to local adaptation. Measurements of sequence differentiation and identification of genes in these regions is important but difficult, especially in organisms with limited genetic information available. We have followed up a genome scan in the marine gastropod, Littorina saxatilis, by searching a bacterial artificial chromosome library with differentiated and undifferentiated markers, sequencing four bacterial artificial chromosomes and then analysing sequence variation in population samples for fragments at, and close to the original marker polymorphisms. We show that sequence differentiation follows the patterns expected from the original marker frequencies, that differentiated markers identify independent and highly localized sites and that these sites fall outside coding regions. Two differentiated loci are characterized by insertions of putative transposable elements that appear to have increased in frequency recently and which might influence expression of downstream genes. These results provide strong candidate loci for the study of local adaptation in Littorina. They demonstrate an approach that can be applied to follow up genome scans in other taxa and they show that the genome scan approach can lead rapidly to candidate genes in nonmodel organisms.

Ecological genetics in the North Atlantic: environmental gradients and adaptation at specific loci.

The North Atlantic intertidal community provides a rich set of organismal and environmental material for the study of ecological genetics. Clearly defined environmental gradients exist at multiple spatial scales: there are broad latitudinal trends in temperature, meso-scale changes in salinity along estuaries, and smaller scale gradients in desiccation and temperature spanning the intertidal range. The geology and geography of the American and European coasts provide natural replication of these gradients, allowing for population genetic analyses of parallel adaptation to environmental stress and heterogeneity. Statistical methods have been developed that provide genomic neutrality tests of population differentiation and aid in the process of candidate gene identification. In this paper, we review studies of marine organisms that illustrate associations between an environmental gradient and specific genetic markers. Such highly differentiated markers become candidate genes for adaptation to the environmental factors in question, but the functional significance of genetic variants must be comprehensively evaluated. We present a set of predictions about locus-specific selection across latitudinal, estuarine, and intertidal gradients that are likely to exist in the North Atlantic. We further present new data and analyses that support and contradict these simple selection models. Some taxa show pronounced clinal variation at certain loci against a background of mild clinal variation at many loci. These cases illustrate the procedures necessary for distinguishing selection driven by internal genomic vs. external environmental factors. We suggest that the North Atlantic intertidal community provides a model system for identifying genes that matter in ecology due to the clarity of the environmental stresses and an extensive experimental literature on ecological function. While these organisms are typically poor genetic and genomic models, advances in comparative genomics have provided access to molecular tools that can now be applied to taxa with well-defined ecologies. As many of the organisms we discuss have tight physiological limits driven by climatic factors, this synthesis of molecular population genetics with marine ecology could provide a sensitive means of assessing evolutionary responses to climate change.

An EST-based genome scan using 454 sequencing in the marine snail Littorina saxatilis

Genome scans have been used in the studies of ecological speciation to find genomic regions (‘outlier loci’) showing reduced gene flow between divergent populations/species. High‐throughput sequencing (‘454’) offers new opportunities in this field via transcriptome sequencing. Divergent ecotypes of the marine gastropod Littorina saxatilis represent a good example of incipient ecological speciation. We performed a 454‐based genome scan between H and M ecotypes of L. saxatilis from the British Isles using cDNA of pooled individuals. Allele frequencies were calculated for 2454 single nucleotide polymorphisms (SNPs), within 572 contigs, and 7% of loci were detected as outliers. Functional annotation of the contigs containing outlier SNPs showed that they included shell matrix and muscle proteins (lithostathine, mucin, titin), proteins involved in energetic metabolism (arginine kinase, NADH dehydrogenase) and reverse transcriptases. Follow‐up investigations into these proteins and unannotated outliers will be a promising route in the study of ecological speciation in L. saxatilis.

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.

Partial sequence of the mitochondrial genome of Littorina saxatilis: relevance to gastropod phylogenetics.

Abstract. A 8022 base pair fragment from the mitochondrial DNA of the prosobranch gastropod Littorina saxatilis has been sequenced and shown to contain the complete genes for 12 transfer RNAs and five protein genes (CoII, ATPase 6, ATPase 8, ND1, ND6), two partial protein genes (CoI and cyt b), and two ribosomal RNAs (small and large subunits). The order of these constituent genes differs from those of other molluscan mitochondrial gene arrangements. Only a single rearrangement involving a block of protein coding genes and three tRNA translocations are necessary to produce identical gene orders between L. saxatilis and K. tunicata. However, only one gene boundary is shared between the L. saxatilis gene order and that of the pulmonate gastropod Cepaea nemoralis. This extends the observation that there is little conservation of mitochrondrial gene order amongst the Mollusca and suggests that radical mitochondrial DNA gene rearrangement has occurred on the branch leading to the pulmonates.

Reproductive effort and r-and K-selection in two species of Lacuna (Gastropoda: Prosobranchia)

The growth and spawning of Lacuna pallidula and L. vincta were measured in the laboratory over 7 months. In both cases, there were significant differences between the mean number of eggs per batch or weight of egg batches from females of the same species and also between rates of spawning by females of the same species. The reproductive effort of L. vincta, estimated by the ratio total spawn weight: body weight and by the time taken to exceed body weight in cumulative spawn, is approximately twice that of L. pallidula. L. vincta has a long-lived planktotrophic larva, while L. pallidula has direct, lecithotrophic development, and therefore in this instance planktotrophic development seems to be more expensive to the parent than lecithotrophic development. Published work on two Pacific seastars leads to the opposite conclusion, and it is suggested that the paradox can be resolved in terms of r-K-selection theory. In both cases, the r-selected species has a higher reproductive effort, notwithstanding that the snail is planktotrophic and the seastar lecithotrophic in development.

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 FST 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.