RADseq and mate choice assays reveal unidirectional gene flow among three lamprey ecotypes despite weak assortative mating: Insights into the formation and stability of multiple ecotypes in sympatry

Contrasting patterns of genetic and phenotypic divergence of two sympatric congeners, Phragmites australis and P. hirsuta, in heterogeneous habitats

Habitat heterogeneity leads to genome-wide differentiation and morphological and ecological differentiation, which will progress along the speciation continuum, eventually leading to speciation. Phragmites hirsuta and Phragmites australis are sympatric congeners that coexist in saline-alkaline meadow soil (SAS) and sandy soil (SS) habitats of the Songnen Meadow. The results provided genetic evidence for two separate species of reeds. Genetic diversity and spatial genetic structure supported the specialist-generalist variation hypothesis (SGVH) in these two sympatric reed species, suggesting that P. australis is a generalist and P. hirsuta is a habitat specialist. When we compared these different species with respect to phenotypic and genetic variation patterns in different habitats, we found that the phenotypic differentiation of P. australis between the two habitats was higher than that of P. hirsuta. Multiple subtle differences in morphology, genetic background, and habitat use collectively contribute to ecological success for similar congeners. This study provided evidence of the two reed congeners, which should contribute to their success in harsh environments.

Integrating Pool-seq uncertainties into demographic inference

Next-generation sequencing of pooled samples (Pool-seq) is a popular method to assess genome-wide diversity patterns in natural and experimental populations. However, Pool-seq is associated with specific sources of noise, such as unequal individual contributions. Consequently, using Pool-seq for the reconstruction of evolutionary history has remained underexplored. Here we describe a novel Approximate Bayesian Computation (ABC) method to infer demographic history, explicitly modelling Pool-seq sources of error. By jointly modelling Pool-seq data, demographic history and the effects of selection due to barrier loci, we obtain estimates of demographic history parameters accounting for technical errors associated with Pool-seq. Our ABC approach is computationally efficient as it relies on simulating subsets of loci (rather than the whole-genome) and on using relative summary statistics and relative model parameters. Our simulation study results indicate Pool-seq data allows distinction between general scenarios of ecotype formation (single versus parallel origin) and to infer relevant demographic parameters (e.g. effective sizes and split times). We exemplify the application of our method to Pool-seq data from the rocky-shore gastropod Littorina saxatilis, sampled on a narrow geographical scale at two Swedish locations where two ecotypes (Wave and Crab) are found. Our model choice and parameter estimates show that ecotypes formed before colonization of the two locations (i.e. single origin) and are maintained despite gene flow. These results indicate that demographic modelling and inference can be successful based on pool-sequencing using ABC, contributing to the development of suitable null models that allow for a better understanding of the genetic basis of divergent adaptation.

Sperm competition experiments reveal low prezygotic postmating isolation between parasitic and nonparasitic lamprey ecotypes

The role of postmating sexual selection as a potential reproductive barrier in speciation is not well understood. Here, we studied the effects of sperm competition and cryptic female choice as putative postmating barriers in two lamprey ecotypes with a partial reproductive isolation. The European river lamprey Lampetra fluviatilis is anadromous and parasitic of other fish species, whereas the brook lamprey Lampetra planeri is freshwater resident and nonparasitic. We measured sperm traits in both ecotypes and designed sperm competition experiments to test the occurrence of cryptic female choice. We also performed sperm competition experiments either at equal semen volume or equal sperm number to investigate the role of sperm velocity on fertilization success. We observed distinct sperm traits between ecotypes with a higher sperm concentration and a lower sperm velocity for L. planeri compared with L. fluviatilis. The outcomes of sperm competition reflected these differences in sperm traits, and there was no evidence for cryptic female choice irrespective of female ecotype. At equal semen volume, L. planeri males had a higher fertilization success than L. fluviatilis and vice versa at equal sperm number. Our results demonstrate that different sperm traits between ecotypes can influence the male reproductive success and thus gene flow between L. planeri and L. fluviatilis. However, postmating prezygotic barriers are absent and thus cannot explain the partial reproductive isolation between ecotypes.

A SNP marker to discriminate the european brook lamprey (Lampetra planeri), river lamprey (L. fluviatilis) and their hybrids

BACKGROUND

The European River lamprey Lampetra fluviatilis and the brook lamprey L. planeri are two closely related species that are also considered as partially reproductively isolated ecotypes. At the larval stage, they cannot be distinguished morphologically or genetically by mitochondrial DNA. We aimed at developing a molecular tool to identify early life stages of L. fluviatilis and L. planeri.

METHODS AND RESULTS

We first identified Single Nucleotide Polymorphism (SNP) markers with distinct alleles between L. fluviatilis and L. planeri based on RADseq data from 186 individuals collected in France. Then, we developed a quantitative PCR protocol to genotype a species-specific SNP, which was tested on 270 samples including larvae and adults from both species and 17 sites from Western Europe. Results were consistent with morphological identification in all cases except for samples from the Rhône drainage and the Loch Lomond in Scotland.

CONCLUSIONS

This marker allows the identification of larval stages of L. fluviatilis, L. planeri as well as hybrids. The two cases where the marker was not diagnostic might be explained by an isolated and distinct L. planeri lineage in the Rhône drainage, and high levels of admixture among L. fluviatilis, L. planeri and the lake-parasitic form in Loch Lomond.

Drivers of neutral and adaptive differentiation in pike (Esox lucius) populations from contrasting environments

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

An assessment of terminology for intraspecific diversity in fishes, with a focus on "ecotypes" and "life histories"

Understanding and preserving intraspecific diversity (ISD) is important for species conservation. However, ISD units do not have taxonomic standards and are not universally recognized. The terminology used to describe ISD is varied and often used ambiguously. We compared definitions of terms used to describe ISD with use in recent studies of three fish taxa: sticklebacks (Gasterosteidae), Pacific salmon and trout (Oncorhynchus spp., "PST"), and lampreys (Petromyzontiformes). Life history describes the phenotypic responses of organisms to environments and includes biological parameters that affect population growth or decline. Life-history pathway(s) are the result of different organismal routes of development that can result in different life histories. These terms can be used to describe recognizable life-history traits. Life history is generally used in organismal- and ecology-based journals. The terms paired species/species pairs have been used to describe two different phenotypes, whereas in some species and situations a continuum of phenotypes may be expressed. Our review revealed overlapping definitions for race and subspecies, and subspecies and ecotypes. Ecotypes are genotypic adaptations to particular environments, and this term is often used in genetic- and evolution-based journals. "Satellite species" is used for situations in which a parasitic lamprey yields two or more derived, nonparasitic lamprey species. Designatable Units, Evolutionary Significant Units (ESUs), and Distinct Population Segments (DPS) are used by some governments to classify ISD of vertebrate species within distinct and evolutionary significant criteria. In situations where the genetic or life-history components of ISD are not well understood, a conservative approach would be to call them phenotypes.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Fine-Scale Population Genetic Structure and Parapatric Cryptic Species of Kuruma Shrimp (Marsupenaeus japonicus), Along the Northwestern Pacific Coast of China

The kuruma shrimp (Marsupenaeus japonicus) includes two cryptic species, which are distributed mostly allopatrically but co-occur in the northern South China Sea (from Huilai to Beihai). To obtain a better understanding of the fine-scale genetic structure and parapatric diversification of these two varieties in the northwestern Pacific region, we used a genotyping-by-sequencing (GBS) and comparative transcriptomics approach to establish their phylogenetic relationships. Using the GBS technique, we genotyped 28891 SNPs in 160 individuals in the Northwest Pacific. The results supported two highly diverged evolutionary lineages of kuruma shrimp (var. I and II). The ND and XM populations showed complex genetic patterns, which might be affected by the complex environment of the Taiwan Strait. In addition, the migration rates and inbreeding coefficients of XM and BH were much lower than those of the other populations, which might be related to the land-sea changes and complex ocean currents in the Taiwan Strait and Qiongzhou Strait. Based on the synonymous substitution rates (ds) of 2,491 candidate orthologs, we estimated that the divergence time between the two varieties was 0.26~0.69 Mya. Choice and no-choice interbreeding experiments provided support for the biological species concept, by showing the existence of reproductive isolation or incompatibility. In view of these differences between the two Marsupenaeus species, we believe that it is essential and urgent to establish a genetic database for each and reevaluate their ecological suitable conditions in order to improve species-specific culturing techniques. Moreover, this research can serve as a case study for future research on speciation and hybridization.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.