Phylogeographical analysis reveals multiple conservation units in brook lampreys Lampetra planeri of Portuguese streams

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.

Transcriptional Basis of Copper-Induced Olfactory Impairment in the Sea Lamprey, a Primitive Invasive Fish

Olfaction mediates behaviors necessary for survival and reproduction in fishes. Anthropogenic inputs of contaminants into aquatic environments, specifically copper, are known to disrupt a broad range of olfactory-mediated behaviors and can cause long-lasting damage even at low concentrations that have profound impacts on the biology of aquatic organisms. The sea lamprey (Petromyzon marinus) is a primitive fish species invasive to the North American Great Lakes that relies on olfaction to navigate during natal homing and in mate choice during reproduction. To investigate effects of copper on sea lamprey olfaction and the potential for maintenance of olfactory function during copper exposure, we exposed juvenile sea lamprey to environmentally ecologically relevant copper concentrations (0, 5, 10 and 30 µg/L) for 24 hr and characterized gene transcription response in olfactory tissue (i.e., peripheral olfactory organ and olfactory bulb) and forebrain using whole transcriptome sequencing. Copper exposure induced a pattern of positive dose-dependent transcriptional response. Expression changes primarily reflected up-regulation of genes involved in apoptosis and wound healing. Unlike higher vertebrates, genes specifically related to the olfactory senses of the sea lamprey, e.g., olfactory receptors, exhibited little transcriptional response to copper exposure, suggesting the mechanism of copper-induced olfactory impairment is through necrosis of the olfactory bulb and not copper-selective inhibition of olfactory receptors. Fully two-thirds of the differentially expressed genes at higher doses of copper have no known function and thus represent important candidates for further study of the responses to copper-induced olfactory injury. Our results shed light on the evolution of vertebrate olfactory repair mechanisms and have important implications for the conservation and management of both invasive and native populations of lamprey.

The mitogenome of the southern Brook Lamprey, Ichthyomyzon gagei (Cyclostomata: Petromyzontidae)

The southern brook lamprey (Ichthyomyzon gagei) is a non-parasitic lamprey endemic to the southeastern US. Here, we report the complete mitogenome of this basal vertebrate and found its genomic organization to be similar to that of other reported lamprey mitogenomes. Nucleotide sequence identities for individual proteins range from 90% to 94% when compared with the congeneric species I. fossor and I. unicuspis. Finally, phylogenetic analysis revealed I. gagei to be sister to these other species of Ichthyomyzon. These genomic data provide a baseline for future investigations regarding the molecular evolution of basal vertebrates.

Inferring the demographic history underlying parallel genomic divergence among pairs of parasitic and nonparasitic lamprey ecotypes

Understanding the evolutionary mechanisms generating parallel genomic divergence patterns among replicate ecotype pairs remains an important challenge in speciation research. We investigated the genomic divergence between the anadromous parasitic river lamprey (Lampetra fluviatilis) and the freshwater-resident nonparasitic brook lamprey (Lampetra planeri) in nine population pairs displaying variable levels of geographic connectivity. We genotyped 338 individuals with RAD sequencing and inferred the demographic divergence history of each population pair using a diffusion approximation method. Divergence patterns in geographically connected population pairs were better explained by introgression after secondary contact, whereas disconnected population pairs have retained a signal of ancient migration. In all ecotype pairs, models accounting for differential introgression among loci outperformed homogeneous migration models. Generating neutral predictions from the inferred divergence scenarios to detect highly differentiated markers identified greater proportions of outliers in disconnected population pairs than in connected pairs. However, increased similarity in the most divergent genomic regions was found among connected ecotype pairs, indicating that gene flow was instrumental in generating parallelism at the molecular level. These results suggest that heterogeneous genomic differentiation and parallelism among replicate ecotype pairs have partly emerged through restricted introgression in genomic islands.

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.

Molecular diversity of Germany's freshwater fishes and lampreys assessed by DNA barcoding

This study represents the first comprehensive molecular assessment of freshwater fishes and lampreys from Germany. We analysed COI sequences for almost 80% of the species mentioned in the current German Red List. In total, 1056 DNA barcodes belonging to 92 species from all major drainages were used to (i) build a reliable DNA barcode reference library, (ii) test for phylogeographic patterns, (iii) check for the presence of barcode gaps between species and (iv) evaluate the performance of the barcode index number (BIN) system, available on the Barcode of Life Data Systems. For over 78% of all analysed species, DNA barcodes are a reliable means for identification, indicated by the presence of barcode gaps. An overlap between intra- and interspecific genetic distances was present in 19 species, six of which belong to the genus Coregonus. The Neighbour-Joining phenogram showed 60 nonoverlapping species clusters and three singleton species, which were related to 63 separate BIN numbers. Furthermore, Barbatula barbatula, Leucaspius delineatus, Phoxinus phoxinus and Squalius cephalus exhibited remarkable levels of cryptic diversity. In contrast, 11 clusters showed haplotype sharing, or low levels of divergence between species, hindering reliable identification. The analysis of our barcode library together with public data resulted in 89 BINs, of which 56% showed taxonomic conflicts. Most of these conflicts were caused by the use of synonymies, inadequate taxonomy or misidentifications. Moreover, our study increased the number of potential alien species in Germany from 14 to 21 and is therefore a valuable groundwork for further faunistic investigations.

Complete mitochondrial genome of the river lamprey, Lampetra japonica (Petromyzontiformes, Petromyzonidae)

The complete mitochondrial genome was sequenced from the river lamprey, Lampetra japonica. The genome sequence was 16,267 bp in size, and the gene order and contents were identical with the genera Lampetra reissneri. However, L. japonica has two control regions with different mitochondrial genomic organization compared to other previously reported fish mitochondrial genomes. Of 13 protein-coding genes (PCGs), 1 gene (CO1) had incomplete stop codon. The base composition of L. japonica showed anti-G bias (4.65%) on the third position of PCGs.