Juan Galindo

Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can ‘go genomic’, it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.

Adaptation genomics: the next generation

Understanding the genetics of how organisms adapt to changing environments is a fundamental topic in modern evolutionary ecology. The field is currently progressing rapidly because of advances in genomics technologies, especially DNA sequencing. The aim of this review is to first briefly summarise how next generation sequencing (NGS) has transformed our ability to identify the genes underpinning adaptation. We then demonstrate how the application of these genomic tools to ecological model species means that we can start addressing some of the questions that have puzzled ecological geneticists for decades such as: How many genes are involved in adaptation? What types of genetic variation are responsible for adaptation? Does adaptation utilise pre-existing genetic variation or does it require new mutations to arise following an environmental change?

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.

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.

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 FST 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 FST partitions: differences between ecotypes (overall and local), differences between localities and micro‐geographical differences within ecotypes. The magnitude of FST differed between candidate and noncandidate loci for all partitions except in the case of microgeographical 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.

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.

Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.

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.

Advances in Ecological Speciation: an integrative approach

The role of natural selection in promoting reproductive isolation has received substantial renewed interest within the last two decades. As a consequence, the study of ecological speciation has become an extremely productive research area in modern evolutionary biology. Recent innovations in sequencing technologies offer an unprecedented opportunity to study the mechanisms involved in ecological speciation. Genome scans provide significant insights but have some important limitations; efforts are needed to integrate them with other approaches to make full use of the sequencing data deluge. An international conference ‘Advances in Ecological Speciation’ organized by the University of Porto (Portugal) aimed to review current progress in ecological speciation. Using some of the examples presented at the conference, we highlight the benefits of integrating ecological and genomic data and discuss different mechanisms of parallel evolution. Finally, future avenues of research are suggested to advance our knowledge concerning the role of natural selection in the establishment of reproductive isolation during ecological speciation.