MtDNA markers reveal the existence of allopatric evolutionary lineages in the threatened lampreys Lampetra fluviatilis (L.) and Lampetra planeri (Bloch) in the Iberian glacial refugium

Seascape genomics reveals limited dispersal and suggests spatially varying selection among European populations of sea lamprey (Petromyzon marinus)

Sea lamprey Petromyzon marinus is an anadromous and semelparous fish without homing behaviors. Despite being a freshwater, free-living organism for a large part of their life cycle, its adulthood is spent as a parasite of marine vertebrates. In their native European range, while it is well-established that sea lampreys comprise a single nearly-panmictic population, few studies have further explored the evolutionary history of natural populations. Here, we performed the first genome-wide characterization of sea lamprey's genetic diversity in their European natural range. The objectives were to investigate the connectivity among river basins and explore evolutionary processes mediating dispersal during the marine phase, with the sequencing of 186 individuals from 8 locations spanning the North Eastern Atlantic coast and the North Sea with double-digest RAD-sequencing, obtaining a total of 30,910 bi-allelic SNPs. Population genetic analyses reinforced the existence of a single metapopulation encompassing freshwater spawning sites within the North Eastern Atlantic and the North Sea, though the prevalence of private alleles at northern latitudes suggested some limits to the species' dispersal. Seascape genomics suggested a scenario where oxygen concentration and river runoffs impose spatially varying selection across their distribution range. Exploring associations with the abundance of potential hosts further suggested that hake and cod could also impose selective pressures, although the nature of such putative biotic interactions was unresolved. Overall, the identification of adaptive seascapes in a panmictic anadromous species could contribute to conservation practices by providing information for restoration activities to mitigate local extinctions on freshwater sites.

Riverscape genetics in brook lamprey: genetic diversity is less influenced by river fragmentation than by gene flow with the anadromous ecotype

Understanding the effect of human-induced landscape fragmentation on gene flow and evolutionary potential of wild populations has become a major concern. Here, we investigated the effect of riverscape fragmentation on patterns of genetic diversity in the freshwater resident European brook lamprey (Lampetra planeri) that has a low ability to pass obstacles to migration. We tested the hypotheses of (i) asymmetric gene flow following water current and (ii) an effect of gene flow with the closely related anadromous river lamprey (L. fluviatilis) ecotype on L. planeri genetic diversity. We genotyped 2472 individuals, including 225 L. fluviatilis, sampled from 81 sites upstream and downstream barriers to migration, in 29 western European rivers. Linear modelling revealed a strong positive relationship between genetic diversity and the distance from the river source, consistent with expected patterns of decreased gene flow into upstream populations. However, the presence of anthropogenic barriers had a moderate effect on spatial genetic structure. Accordingly, we found evidence for downstream-directed gene flow, supporting the hypothesis that barriers do not limit dispersal mediated by water flow. Downstream L. planeri populations in sympatry with L. fluviatilis displayed consistently higher genetic diversity. We conclude that genetic drift and slight downstream gene flow drive the genetic make-up of upstream L. planeri populations whereas gene flow between ecotypes maintains higher levels of genetic diversity in L. planeri populations sympatric with L. fluviatilis. We discuss the implications of these results for the design of conservation strategies of lamprey, and other freshwater organisms with several ecotypes, in fragmented dendritic river networks.

Phylogeography of the European brook lamprey (Lampetra planeri) and the European river lamprey (Lampetra fluviatilis) species pair based on mitochondrial data

The European river lamprey Lampetra fluviatilis and the European brook lamprey Lampetra planeri (Block 1784) are classified as a paired species, characterized by notably different life histories but morphological similarities. Previous work has further shown limited genetic differentiation between these two species at the mitochondrial DNA level. Here, we expand on this previous work, which focused on lamprey species from the Iberian Peninsula in the south and mainland Europe in the north, by sequencing three mitochondrial marker regions of Lampetra individuals from five river systems in Ireland and five in southern Italy. Our results corroborate the previously identified pattern of genetic diversity for the species pair. We also show significant genetic differentiation between Irish and mainland European lamprey populations, suggesting another ichthyogeographic district distinct from those previously defined. Finally, our results stress the importance of southern Italian L. planeri populations, which maintain several private alleles and notable genetic diversity.

Genetic variation and population structure among larval Lethenteron spp. within the Yukon River drainage, Alaska

The absence of information on genetic variation and population structure of lampreys Lethenteron spp. in the eastern part of their distribution limits our understanding of the migration ecology and spatial population genetic structure of the species. We examined genetic variation within and among three aggregations of Lethenteron spp. larvae in the Yukon River drainage, Alaska, using microsatellite genotypes. A total of 120 larval lampreys were genotyped at eight microsatellite loci. Global F_ST was 0.053 (95% CI 0.021-0.086), while pairwise F_ST values ranged from 0.048-0.057. Model-based Bayesian clustering analyses with sample locality priors (LOCPRIOR) identified three distinct, but admixed, genetic clusters that corresponded with the three aggregations. Estimates of contemporary gene flow indicate substantial reciprocal migration among sites consistent with no or low-fidelity natal homing. These results are largely in agreement with previous reports of historic and contemporary gene flow among Lethenteron spp. in other parts of their geographic distribution.

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.

European Lampreys: New Insights on Postglacial Colonization, Gene Flow and Speciation

Ice ages are known to be the most dominant palaeoclimatic feature occurring on Earth, producing severe climatic oscillations and consequently shaping the distribution and the population structure of several species. Lampreys constitute excellent models to study the colonization of freshwater systems, as they commonly appear in pairs of closely related species of anadromous versus freshwater resident adults, thus having the ability to colonize new habitats, through the anadromous species, and establish freshwater resident derivates. We used 10 microsatellite loci to investigate the spatial structure, patterns of gene flow and migration routes of Lampetra populations in Europe. We sampled 11 populations including the migratory L. fluviatilis and four resident species, L. planeri, L. alavariensis, L. auremensis and L. lusitanica, the last three endemic to the Iberian Peninsula. In this southern glacial refugium almost all sampled populations represent a distinct genetic cluster, showing high levels of allopatric differentiation, reflecting long periods of isolation. As result of their more recent common ancestor, populations from northern Europe are less divergent among them, they are represented by fewer genetic clusters, and there is evidence of strong recent gene flow among populations. These previously glaciated areas from northern Europe may have been colonized from lampreys expanding out of the Iberian refugia. The pair L. fluviatilis/L. planeri is apparently at different stages of speciation in different locations, showing evidences of high reproductive isolation in the southern refugium, and low differentiation in the north.

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater-resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size-assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.

Contrasting population genetic structure among freshwater-resident and anadromous lampreys: the role of demographic history, differential dispersal and anthropogenic barriers to movement

The tendency of many species to abandon migration remains a poorly understood aspect of evolutionary biology that may play an important role in promoting species radiation by both allopatric and sympatric mechanisms. Anadromy inherently offers an opportunity for the colonization of freshwater environments, and the shift from an anadromous to a wholly freshwater life history has occurred in many families of fishes. Freshwater-resident forms have arisen repeatedly among lampreys (within the Petromyzontidae and Mordaciidae), and there has been much debate as to whether anadromous lampreys, and their derived freshwater-resident analogues, constitute distinct species or are divergent ecotypes of polymorphic species. Samples of 543 European river lamprey Lampetra fluviatilis (mostly from anadromous populations) and freshwater European brook lamprey Lampetra planeri from across 18 sites, primarily in the British Isles, were investigated for 13 polymorphic microsatellite DNA loci, and 108 samples from six of these sites were sequenced for 829 bp of mitochondrial DNA (mtDNA). We found contrasting patterns of population structure for mtDNA and microsatellite DNA markers, such that low diversity and little structure were seen for all populations for mtDNA (consistent with a recent founder expansion event), while fine-scale structuring was evident for nuclear markers. Strong differentiation for microsatellite DNA loci was seen among freshwater-resident L. planeri populations and between L. fluviatilis and L. planeri in most cases, but little structure was evident among anadromous L. fluviatilis populations. We conclude that postglacial colonization founded multiple freshwater-resident populations with strong habitat fidelity and limited dispersal tendencies that became highly differentiated, a pattern that was likely intensified by anthropogenic barriers.

Different stocks of brook lamprey in Spain and their origin from Lampetra fluviatilis at two distinct times and places

Using the mitochondrial non-coding region I, it was shown that the two Spanish Lampetra planeri populations (Cares-Deva and Olabidea-Ugarana) correspond to different genetic units. The Cares-Deva population is probably a recent offshoot of Lampetra fluviatilis, being the Olabidea-Ugarana population less diverse and of older origin.

Genetic diversity, endemism and phylogeny of lampreys within the genus Lampetra sensu stricto (Petromyzontiformes: Petromyzontidae) in western North America

Phylogenetic structure of four Lampetra species from the Pacific drainage of North America (western brook lamprey Lampetra richardsoni, Pacific brook lamprey Lampetra pacifica, river lamprey Lampetra ayresii and Kern brook lamprey Lampetra hubbsi) and unidentified Lampetra specimens (referred to as Lampetra sp.) from 36 locations was estimated using the mitochondrial cytochrome b gene. Maximum parsimony and Bayesian inferences did not correspond with any taxonomic scheme proposed to date. Rather, although L. richardsoni (from Alaska to California) and L. ayresii (from British Columbia to California) together constituted a well-supported clade distinct from several genetically divergent Lampetra populations in Oregon and California, these two species were not reciprocally monophyletic. The genetically divergent populations included L. pacifica (from the Columbia River basin) and L. hubbsi (from the Kern River basin) and four Lampetra sp. populations in Oregon (Siuslaw River and Fourmile Creek) and California (Kelsey and Mark West Creeks). These four Lampetra sp. populations showed genetic divergence between 2.3 and 5.7% from any known species (and up to 8.0% from each other), and may represent morphologically cryptic and thus previously undescribed species. A fifth population (from Paynes Creek, California) may represent a range extension of L. hubbsi into the Upper Sacramento River.