Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin

Parallel evolution despite low genetic diversity in three-spined sticklebacks

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south.

Cast away in the Adriatic: Low degree of parallel genetic differentiation in three-spined sticklebacks

The three-spined stickleback (Gasterosteus aculeatus) has repeatedly and independently adapted to freshwater habitats from standing genetic variation (SGV) following colonization from the sea. However, in the Mediterranean Sea G. aculeatus is believed to have gone extinct, and thus the spread of locally adapted alleles between different freshwater populations via the sea since then has been highly unlikely. This is expected to limit parallel evolution, that is the extent to which phylogenetically related alleles can be shared among independently colonized freshwater populations. Using whole genome and 2b-RAD sequencing data, we compared levels of genetic differentiation and genetic parallelism of 15 Adriatic stickleback populations to 19 Pacific, Atlantic and Caspian populations, where gene flow between freshwater populations across extant marine populations is still possible. Our findings support previous studies suggesting that Adriatic populations are highly differentiated (average F_ST  ≈ 0.45), of low genetic diversity and connectivity, and likely to stem from multiple independent colonizations during the Pleistocene. Linkage disequilibrium network analyses in combination with linear mixed models nevertheless revealed several parallel marine-freshwater differentiated genomic regions, although still not to the extent observed elsewhere in the world. We hypothesize that current levels of genetic parallelism in the Adriatic lineages are a relic of freshwater adaptation from SGV prior to the extinction of marine sticklebacks in the Mediterranean that has persisted despite substantial genetic drift experienced by the Adriatic stickleback isolates.

Threespine Stickleback in Lake Constance: The Ecology and Genomic Substrate of a Recent Invasion

Invasive species can be powerful models for studying contemporary evolution in natural environments. As invading organisms often encounter new habitats during colonization, they will experience novel selection pressures. Threespine stickleback (Gasterosteus aculeatus complex) have recently colonized large parts of Switzerland and are invasive in Lake Constance. Introduced to several watersheds roughly 150 years ago, they spread across the Swiss Plateau (400–800 m a.s.l.), bringing three divergent hitherto allopatric lineages into secondary contact. As stickleback have colonized a variety of different habitat types during this recent range expansion, the Swiss system is a useful model for studying contemporary evolution with and without secondary contact. For example, in the Lake Constance region there has been rapid phenotypic and genetic divergence between a lake population and some stream populations. There is considerable phenotypic variation within the lake population, with individuals foraging in and occupying littoral, offshore pelagic, and profundal waters, the latter of which is a very unusual habitat for stickleback. Furthermore, adults from the lake population can reach up to three times the size of adults from the surrounding stream populations, and are large by comparison to populations globally. Here, we review the historical origins of the threespine stickleback in Switzerland, and the ecomorphological variation and genomic basis of its invasion in Lake Constance. We also outline the potential ecological impacts of this invasion, and highlight the interest for contemporary evolution studies.

Multiple waves of freshwater colonization of the three-spined stickleback in the Japanese Archipelago

Background

The three-spined stickleback (Gasterosteus aculeatus) is a remarkable system to study the genetic mechanisms underlying parallel evolution during the transition from marine to freshwater habitats. Although the majority of previous studies on the parallel evolution of sticklebacks have mainly focused on postglacial freshwater populations in the Pacific Northwest of North America and northern Europe, we recently use Japanese stickleback populations for investigating shared and unique features of adaptation and speciation between geographically distant populations. However, we currently lack a comprehensive phylogeny of the Japanese three-spined sticklebacks, despite the fact that a good phylogeny is essential for any evolutionary and ecological studies. Here, we conducted a phylogenomic analysis of the three-spined stickleback in the Japanese Archipelago.

Results

We found that freshwater colonization occurred in multiple waves, each of which may reflect different interglacial isolations. Some of the oldest freshwater populations from the central regions of the mainland of Japan (hariyo populations) were estimated to colonize freshwater approximately 170,000 years ago. The next wave of colonization likely occurred approximately 100,000 years ago. The inferred origins of several human-introduced populations showed that introduction occurred mainly from nearby habitats. We also found a new habitat of the three-spined stickleback sympatric with the Japan Sea stickleback (Gasterosteus nipponicus).

Conclusions

These Japanese stickleback systems differ from those in the Pacific Northwest of North America and northern Europe in terms of divergence time and history. Stickleback populations in the Japanese Archipelago offer valuable opportunities to study diverse evolutionary processes in historical and contemporary timescales.

Admixture between old lineages facilitated contemporary ecological speciation in Lake Constance stickleback

Ecological speciation can sometimes rapidly generate reproductively isolated populations coexisting in sympatry, but the origin of genetic variation permitting this is rarely known. We previously explored the genomics of very recent ecological speciation into lake and stream ecotypes in stickleback from Lake Constance. Here, we reconstruct the origin of alleles underlying ecological speciation by combining demographic modelling on genome-wide single nucleotide polymorphisms, phenotypic data and mitochondrial sequence data in the wider European biogeographical context. We find that parallel differentiation between lake and stream ecotypes across replicate lake-stream ecotones resulted from recent secondary contact and admixture between old East and West European lineages. Unexpectedly, West European alleles that introgressed across the hybrid zone at the western end of the lake, were recruited to genomic islands of differentiation between ecotypes at the eastern end of the lake. Our results highlight an overlooked outcome of secondary contact: ecological speciation facilitated by admixture variation.

Ecological speciation can proceed rapidly, but the origin of genetic variation facilitating it has remained elusive. Here, the authors show that secondary contact and introgression between deeply diverged lineages of stickleback fish facilitated rapid ecological speciation into lake and stream ecotypes in Lake Constance.

Phylogeography and Conservation Genetics of the Ibero-Balearic Three-Spined Stickleback (Gasterosteus aculeatus)

Genetic isolation and drift may imperil peripheral populations of wide-ranging species more than central ones. Therefore, information about species genetic variability and population structure is invaluable for conservation managers. The Iberian populations of three-spined stickleback lie at the southwestern periphery of the European distribution of Gasterosteus aculeatus. This teleost is a protected species in Portugal and Spain and local extinctions have been reported in both countries during the last decades. Our objectives were (i) to determine whether the Iberian populations of G. aculeatus are unique or composed of any of the major evolutionary lineages previously identified and (ii) to assess the evolutionary potential of these peripheral populations. We genotyped 478 individuals from 17 sites at 10 polymorphic microsatellite loci to evaluate the genetic variability and differentiation of the Ibero-Balearic populations. We also sequenced 1,165 bp of the mitochondrial genome in 331 of those individuals in order to complement the estimates of genetic diversity in the Ibero-Balearic region. We predicted the evolutionary potential of the different sites analysed based on the contribution of each of them to total allelic/mitochondrial diversity. An intraspecific phylogeny at European level was reconstructed using our data from the mitochondrial cytochrome b gene (755 bp) and published sequences. The so-called Transatlantic, European and Mediterranean mitochondrial lineages were found to be present in the Ibero-Balearic region. Their phylogeography suggests a history of multiple colonisations. The nuclear results show, however, a strong correlation between population structure and drainage system. The following basins should be prioritised by conservation policies in order to preserve those populations with the highest evolutionary potential: the Portuguese Vouga and Tagus as well as the Spanish Majorca and Limia. Maintenance of their connectivity, control of exotic species and monitoring of habitat properties are strongly recommended in those areas. Genetic variation alone cannot, however, ensure the persistence of these peripheral southern populations of G. aculeatus. On the one hand, the analysis of a historical sample from Eastern Spain (Penyscola) revealed no genetic erosion, which suggests a fairly sudden extinction of that population. On the other hand, the reintroduction program implemented in the Valencian Community has mostly failed despite our finding of similar level of genetic diversity between the wild source and the captive-bred released individuals.

Displaced phylogeographic signals from Gyrodactylus arcuatus, a parasite of the three-spined stickleback Gasterosteus aculeatus, suggest freshwater glacial refugia in Europe

We examined the global mitochondrial phylogeography of Gyrodactylus arcuatus, a flatworm ectoparasite of three-spined stickleback Gasterosteus aculeatus. In accordance with the suggested high divergence rate of 13%/million years, the genetic variation of the parasite was high: haplotype diversity h=0.985 and nucleotide diversity π=0.0161. The differentiation among the parasite populations was substantial (Φst=0.759), with two main allopatric clades (here termed Euro and North) accounting for 54% of the total genetic variation. The diversity center of the Euro clade was in the Baltic Sea, while the North clade was spread across the Barents and White Seas. A single haplotype within the North clade was found in the western and eastern Pacific Ocean. Divergence of main clades was estimated to be circa 200 thousand years ago. Each main clade was further divided into six distinct subclades, estimated to have diverged in isolation since 135 thousand years ago. This second division corresponds approximately to the Eemian interglacial predating the last glacial maximum. A demographic expansion of the subclades is associated with colonisation of northern Europe since the last glacial maximum, circa 15-40 thousand years ago. The parasite phylogeny is most likely explained by sequential isolated bottlenecks and expansions in numerous allopatric refugia. The postglacial intermingling and high variation in the marine parasite populations, separately in the Baltic and Barents Seas, suggest low competition of divergent parasite matrilines, coupled with a large population size and high rate of dispersal of hosts. The genetic contribution of the assumed refugial fish populations maintaining the parasite during the last glacial maximum was not detected among the marine sticklebacks, which perhaps were infected after range expansion.

Comparative analysis of Japanese three-spined stickleback clades reveals the Pacific Ocean lineage has adapted to freshwater environments while the Japan Sea has not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.

Unique mitochondrial DNA lineages in Irish stickleback populations: cryptic refugium or rapid recolonization?

Repeated recolonization of freshwater environments following Pleistocene glaciations has played a major role in the evolution and adaptation of anadromous taxa. Located at the western fringe of Europe, Ireland and Britain were likely recolonized rapidly by anadromous fishes from the North Atlantic following the last glacial maximum (LGM). While the presence of unique mitochondrial haplotypes in Ireland suggests that a cryptic northern refugium may have played a role in recolonization, no explicit test of this hypothesis has been conducted. The three-spined stickleback is native and ubiquitous to aquatic ecosystems throughout Ireland, making it an excellent model species with which to examine the biogeographical history of anadromous fishes in the region. We used mitochondrial and microsatellite markers to examine the presence of divergent evolutionary lineages and to assess broad-scale patterns of geographical clustering among postglacially isolated populations. Our results confirm that Ireland is a region of secondary contact for divergent mitochondrial lineages and that endemic haplotypes occur in populations in Central and Southern Ireland. To test whether a putative Irish lineage arose from a cryptic Irish refugium, we used approximate Bayesian computation (ABC). However, we found no support for this hypothesis. Instead, the Irish lineage likely diverged from the European lineage as a result of postglacial isolation of freshwater populations by rising sea levels. These findings emphasize the need to rigorously test biogeographical hypothesis and contribute further evidence that postglacial processes may have shaped genetic diversity in temperate fauna.

Microscale vicariance and diversification of Western Balkan caddisflies linked to karstification

The karst areas in the Dinaric region of the Western Balkan Peninsula are a hotspot of freshwater biodiversity. Many investigators have examined diversification of the subterranean freshwater fauna in these karst systems. However, diversification of surface-water fauna remains largely unexplored. We assessed local and regional diversification of surface-water species in karst systems and asked whether patterns of population differentiation could be explained by dispersal-diversification processes or allopatric diversification following karst-related microscale vicariance. We analyzed mitochondrial cytochrome c oxidase subunit I (mtCOI) sequence data of 4 caddisfly species (genus Drusus) in a phylogeographic framework to assess local and regional population genetic structure and Pliocene/Pleistocene history. We used BEAST software to assess the timing of intraspecific diversification of the target species. We compared climate envelopes of the study species and projected climatically suitable areas during the last glacial maximum (LGM) to assess differences in the species climatic niches and infer potential LGM refugia. The haplotype distribution of the 4 species (324 individuals from 32 populations) was characterized by strong genetic differentiation with few haplotypes shared among populations (16%) and deep divergence among populations of the 3 endemic species, even at local scales. Divergence among local populations of endemics often exceeded divergence among regional and continental clades of the widespread D. discolor. Major divergences among regional populations dated to 2.0 to 0.5 Mya. Species distribution model projections and genetic structure suggest that the endemic species persisted in situ and diversified locally throughout multiple Pleistocene climate cycles. The pattern for D. discolor was different and consistent with multiple invasions into the region. Patterns of population genetic structure and diversification were similar for the 3 regional endemic Drusus species and consistent with microscale vicariance after the onset of intensified karstification in the Dinaric region. Karstification may induce microscale vicariance of running surface-water habitats and probably promotes allopatric fragmentation of stream insects at small spatial scales.