Anja Marie Westram

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Shared and nonshared genomic divergence in parallel ecotypes of Littorina saxatilis at a local scale

Parallel speciation occurs when selection drives repeated, independent adaptive divergence that reduces gene flow between ecotypes. Classical examples show parallel speciation originating from shared genomic variation, but this does not seem to be the case in the rough periwinkle (Littorina saxatilis) that has evolved considerable phenotypic diversity across Europe, including several distinct ecotypes. Small ‘wave’ ecotype snails inhabit exposed rocks and experience strong wave action, while thick‐shelled, ‘crab’ ecotype snails are larger and experience crab predation on less exposed shores. Crab and wave ecotypes appear to have arisen in parallel, and recent evidence suggests only marginal sharing of molecular variation linked to evolution of similar ecotypes in different parts of Europe. However, the extent of genomic sharing is expected to increase with gene flow and more recent common ancestry. To test this, we used de novo RAD‐sequencing to quantify the extent of shared genomic divergence associated with phenotypic similarities amongst ecotype pairs on three close islands (<10 km distance) connected by weak gene flow (Nm ~ 0.03) and with recent common ancestry (<10 000 years). After accounting for technical issues, including a large proportion of null alleles due to a large effective population size, we found ~8–28% of positive outliers were shared between two islands and ~2–9% were shared amongst all three islands. This low level of sharing suggests that parallel phenotypic divergence in this system is not matched by shared genomic divergence despite a high probability of gene flow and standing genetic variation.

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.

Spatial Distribution of Cryptic Species Diversity in European Freshwater Amphipods (Gammarus fossarum) as Revealed by Pyrosequencing

In order to understand and protect ecosystems, local gene pools need to be evaluated with respect to their uniqueness. Cryptic species present a challenge in this context because their presence, if unrecognized, may lead to serious misjudgement of the distribution of evolutionarily distinct genetic entities. In this study, we describe the current geographical distribution of cryptic species of the ecologically important stream amphipod Gammarus fossarum (types A, B and C). We use a novel pyrosequencing assay for molecular species identification and survey 62 populations in Switzerland, plus several populations in Germany and eastern France. In addition, we compile data from previous publications (mainly Germany). A clear transition is observed from type A in the east (Danube and Po drainages) to types B and, more rarely, C in the west (Meuse, Rhone, and four smaller French river systems). Within the Rhine drainage, the cryptic species meet in a contact zone which spans the entire G. fossarum distribution range from north to south. This large-scale geographical sorting indicates that types A and B persisted in separate refugia during Pleistocene glaciations. Within the contact zone, the species rarely co-occur at the same site, suggesting that ecological processes may preclude long-term coexistence. The clear phylogeographical signal observed in this study implies that, in many parts of Europe, only one of the cryptic species is present.

Diversity and Distribution of Freshwater Amphipod Species in Switzerland (Crustacea: Amphipoda)

Amphipods are key organisms in many freshwater systems and contribute substantially to the diversity and functioning of macroinvertebrate communities. Furthermore, they are commonly used as bioindicators and for ecotoxicological tests. For many areas, however, diversity and distribution of amphipods is inadequately known, which limits their use in ecological and ecotoxicological studies and handicaps conservation initiatives. We studied the diversity and distribution of amphipods in Switzerland (Central Europe), covering four major drainage basins, an altitudinal gradient of>2,500 m, and various habitats (rivers, streams, lakes and groundwater). We provide the first provisional checklist and detailed information on the distribution and diversity of all amphipod species from Switzerland. In total, we found 29 amphipod species. This includes 16 native and 13 non-native species, one of the latter (Orchestia cavimana) reported here for the first time for Switzerland. The diversity is compared to neighboring countries. We specifically discuss species of the genus Niphargus, which are often receiving less attention. We also found evidence of an even higher level of hidden diversity, and the potential occurrence of further cryptic species. This diversity reflects the biogeographic past of Switzerland, and suggests that amphipods are ideally suited to address questions on endemism and adaptive radiations, post-glaciation re-colonization and invasion dynamics as well as biodiversity-ecosystem functioning relationships in aquatic systems.

Are cryptic host species also cryptic to parasites? Host specificity and geographical distribution of acanthocephalan parasites infecting freshwater Gammarus

Many parasites infect multiple host species. In coevolving host-parasite interactions, theory predicts that parasites should be adapted to locally common hosts, which could lead to regional shifts in host preferences. We studied the interaction between freshwater Gammarus (Crustacea, Amphipoda) and their acanthocephalan parasites using a large-scale field survey and experiments, combined with molecular identification of cryptic host and parasite species. Gammarus pulex is a common host for multiple species of Acanthocephala in Europe but, in Switzerland, is less common than two cryptic members of the Gammarus fossarum species complex (type A and type B). We found that natural populations of these cryptic species were frequently infected by Pomphorhynchus tereticollis and Polymorphus minutus. Four additional parasite species occurred only locally. Parasites were more common in G. fossarum type B than in type A. Infection experiments using several host and parasite sources confirmed consistently lower infection rates in G. pulex than in G. fossarum type A, suggesting a general difference in susceptibility between the two species. In conclusion, we could show that cryptic host species differ in their interactions with parasites, but that these differences were much less dramatic than differences between G. fossarum (type A) and G. pulex. Our data suggest that the acanthocephalans in Switzerland have adapted to the two most common Gammarus species in this region where host species frequencies differ from near-by regions in Europe.

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations – ‘outlier loci’ (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false‐positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated FST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated FST. Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of FST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased FST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated FST in geographically distant locations. We discuss possible explanations for this unexpected result.

Hidden biodiversity in an ecologically important freshwater amphipod: differences in genetic structure between two cryptic species.

Cryptic species, i.e. species that are morphologically hard to distinguish, have been detected repeatedly in various taxa and ecosystems. In order to evaluate the importance of this finding, we have to know in how far cryptic species differ in various aspects of their biology. The amphipod Gammarus fossarum is a key invertebrate in freshwater streams and contains several cryptic species. We examined the population genetic structure, genetic diversity and demographic history of two of them (type A and type B) using microsatellite markers and asked whether they show significant differences. We present results of population genetic analyses based on a total of 37 populations from the headwaters of two major European drainages, Rhine and Rhone. We found that, in both species, genetic diversity was geographically structured among and within drainages. For type A in the Rhine and type B in the Rhone, we detected significant patterns of isolation by distance. The increase of genetic differentiation with geographical distance, however, was much higher in type A than in type B. This result indicates substantial interspecific differences in population history and/or the extent of current gene flow between populations. In the Rhine, type B does not show evidence of isolation by distance, and population differentiation is relatively low across hundreds of kilometres. The majority of these populations also show signatures of recent bottlenecks. These patterns are consistent with a recent expansion of type B into the Rhine drainage. In summary, our results suggest considerable and previously unrecognized interspecific differences in the genetic structure of these cryptic keystone species.

Intrinsic incompatibilities evolving as a by-product of divergent ecological selection: Considering them in empirical studies on divergence with gene flow

The possibility of intrinsic barriers to gene flow is often neglected in empirical research on local adaptation and speciation with gene flow, for example when interpreting patterns observed in genome scans. However, we draw attention to the fact that, even with gene flow, divergent ecological selection may generate intrinsic barriers involving both ecologically selected and other interacting loci. Mechanistically, the link between the two types of barriers may be generated by genes that have multiple functions (i.e., pleiotropy), and/or by gene interaction networks. Because most genes function in complex networks, and their evolution is not independent of other genes, changes evolving in response to ecological selection can generate intrinsic barriers as a by‐product. A crucial question is to what extent such by‐product barriers contribute to divergence and speciation—that is whether they stably reduce gene flow. We discuss under which conditions by‐product barriers may increase isolation. However, we also highlight that, depending on the conditions (e.g., the amount of gene flow and the strength of selection acting on the intrinsic vs. the ecological barrier component), the intrinsic incompatibility may actually destabilize barriers to gene flow. In practice, intrinsic barriers generated as a by‐product of divergent ecological selection may generate peaks in genome scans that cannot easily be interpreted. We argue that empirical studies on divergence with gene flow should consider the possibility of both ecological and intrinsic barriers. Future progress will likely come from work combining population genomic studies, experiments quantifying fitness and molecular studies on protein function and interactions.