Population genetics reveals divergent lineages and ongoing hybridization in a declining migratory fish species complex.

Deciphering the effects of historical and recent demographic processes responsible for the spatial patterns of genetic diversity and structure is a key objective in evolutionary and conservation biology. Using population genetic analyses, we investigated the demographic history, the contemporary genetic diversity and structure, and the occurrence of hybridization and introgression of two species of anadromous fish with contrasting life history strategies and which have undergone recent demographic declines, the allis shad (Alosa alosa) and the twaite shad (Alosa fallax). We genotyped 706 individuals from 20 rivers and 5 sites at sea in Southern Europe at thirteen microsatellite markers. Genetic structure between populations was lower for the nearly semelparous species A. alosa, which disperses greater distances compared to the iteroparous species, A. fallax. Individuals caught at sea were assigned at the river level for A. fallax and at the region level for A. alosa. Using an approximate Bayesian computation framework, we inferred that the most likely long term historical divergence scenario between both species and lineages involved historical separation followed by secondary contact accompanied by strong population size decline. Accordingly, we found evidence for contemporary hybridization and bidirectional introgression due to gene flow between both species and lineages. Moreover, our results support the existence of at least one distinct species in the Mediterrannean sea: A. agone in Golfe du Lion area, and another divergent lineage in Corsica. Overall, our results shed light on the interplay between historical and recent demographic processes and life history strategies in shaping population genetic diversity and structure of closely related species. The recent demographic decline of these species' populations and their hybridization should be carefully considered while implementing conservation programs.

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.

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.

Expansions, diversification, and interindividual copy number variations of AID/APOBEC family cytidine deaminase genes in lampreys.

Cytidine deaminases of the AID/APOBEC family catalyze C-to-U nucleotide transitions in mRNA or DNA. Members of the APOBEC3 branch are involved in antiviral defense, whereas AID contributes to diversification of antibody repertoires in jawed vertebrates via somatic hypermutation, gene conversion, and class switch recombination. In the extant jawless vertebrate, the lamprey, two members of the AID/APOBEC family are implicated in the generation of somatic diversity of the variable lymphocyte receptors (VLRs). Expression studies linked CDA1 and CDA2 genes to the assembly of VLRA/C genes in T-like cells and the VLRB genes in B-like cells, respectively. Here, we identify and characterize several CDA1-like genes in the larvae of different lamprey species and demonstrate that these encode active cytidine deaminases. Structural comparisons of the CDA1 variants highlighted substantial differences in surface charge; this observation is supported by our finding that the enzymes require different conditions and substrates for optimal activity in vitro. Strikingly, we also found that the number of CDA-like genes present in individuals of the same species is variable. Nevertheless, irrespective of the number of different CDA1-like genes present, all lamprey larvae have at least one functional CDA1-related gene encoding an enzyme with predicted structural and chemical features generally comparable to jawed vertebrate AID. Our findings suggest that, similar to APOBEC3 branch expansion in jawed vertebrates, the AID/APOBEC family has undergone substantial diversification in lamprey, possibly indicative of multiple distinct biological roles.

Intra- and interspecific densities shape life-history traits in a salmonid population.

Population dynamics can be regulated through intra- and interspecific density dependence. In species with close ecological requirements, interspecific competition for resources may add to intraspecific density, or even exceed its effect; it may impact single or multiple traits. However, the relative impact of intra- and interspecific densities on demographic parameters has been rarely empirically assessed. We analyzed 18 years of capture-mark-recapture data from brown trout (Salmo trutta) coexisting with Atlantic salmon (Salmo salar) during the juvenile freshwater phase in the Oir River (France) to estimate the relative effects of intra- and interspecific density on trout early life. In trout, a species with optional migration, we estimated the migration probability of young-of-the-year trout out of their natal site, survival probability during the first winter, as well as body size, in relation to both intra- and interspecific density. Trout density correlated negatively with body size and with winter survival in resident trout but not with trout migration. Salmon density correlated positively with trout migration, but no impact was detected on trout body size or survival. Our study highlighted contrasting effects of intra- and interspecific density on trout early life, and the need to account for both factors when studying population dynamics in coexisting species. In particular, by affecting trout migration decision, salmon density may drive trout life history.

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.

Higher temperature exacerbates the impact of sediments on embryo performances in a salmonid.

In a warming climate, higher temperatures are likely to modulate positively or negatively the effect of other environmental factors on biota, although such interactions are poorly documented. Here, we explore under controlled conditions the combined effects of two common stressors in freshwater ecosystems, higher temperature and sediment load, on the embryonic development of arctic charr (Salvelinus alpinus L.). In the warm treatment, embryos had a lower survival, a longer incubation period and a smaller body size with a bigger yolk sac volume. Our data show a significant interaction between temperature and sediment load with temperature increasing dramatically the negative effects of sediment load on fitness-related traits. In the climate change context, these findings highlight the importance of taking into account different thermal scenarios when examining the effect of environmental or anthropogenic stressors.

Genetic and phenotypic changes in an Atlantic salmon population supplemented with non-local individuals: a longitudinal study over 21 years.

While introductions and supplementations using non-native and potentially domesticated individuals may have dramatic evolutionary effects on wild populations, few studies documented the evolution of genetic diversity and life-history traits in supplemented populations. Here, we investigated year-to-year changes from 1989 to 2009 in genetic admixture at 15 microsatellite loci and in phenotypic traits in an Atlantic salmon (Salmo salar) population stocked during the first decade of this period with two genetically and phenotypically distinct source populations. We detected a pattern of temporally increasing introgressive hybridization between the stocked population and both source populations. The proportion of fish returning to the river after a single winter at sea (versus several ones) was higher in fish assigned to the main source population than in local individuals. Moreover, during the first decade of the study, both single-sea-winter and multi-sea-winter (MSW) fish assigned to the main source population were smaller than local fish. During the second decade of the study, MSW fish defined as hybrids were lighter and smaller than fish from parental populations, suggesting outbreeding depression. Overall, this study suggests that supplementation with non-local individuals may alter not only the genetic diversity of wild populations but also life-history traits of adaptive significance.

Persistent unequal sex ratio in a population of grayling (Salmonidae) and possible role of temperature increase.

In some fishes, water chemistry or temperature affects sex determination or creates sex-specific selection pressures. The resulting population sex ratios are hard to predict from laboratory studies if the environmental triggers interact with other factors, whereas in field studies, singular observations of unusual sex ratios may be particularly prone to selective reporting. Long-term monitoring largely avoids these problems. We studied a population of grayling (Thymallus thymallus) in Lake Thun, Switzerland, that has been monitored since 1948. Samples of spawning fish have been caught about 3 times/week around spawning season, and water temperature at the spawning site has been continuously recorded since 1970. We used scale samples collected in different years to determine the average age of spawners (for life-stage specific analyses) and to identify the cohort born in 2003 (an extraordinarily warm year). Recent tissue samples were genotyped on microsatellite markers to test for genetic bottlenecks in the past and to estimate the genetically effective population size (N(e)). Operational sex ratios changed from approximately 65% males before 1993 to approximately 85% males from 1993 to 2011. Sex ratios correlated with the water temperatures the fish experienced in their first year of life. Sex ratios were best explained by the average temperature juvenile fish experienced during their first summer. Grayling abundance is declining, but we found no evidence of a strong genetic bottleneck that would explain the apparent lack of evolutionary response to the unequal sex ratio. Results of other studies show no evidence of endocrine disruptors in the study area. Our findings suggest temperature affects population sex ratio and thereby contributes to population decline.

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.

Determinants of hierarchical genetic structure in Atlantic salmon populations: environmental factors vs. anthropogenic influences.

Disentangling the effects of natural environmental features and anthropogenic factors on the genetic structure of endangered populations is an important challenge for conservation biology. Here, we investigated the combined influences of major environmental features and stocking with non-native fish on the genetic structure and local adaptation of Atlantic salmon (Salmo salar) populations. We used 17 microsatellite loci to genotype 975 individuals originating from 34 French rivers. Bayesian analyses revealed a hierarchical genetic structure into five geographically distinct clusters. Coastal distance, geological substrate and river length were strong predictors of population structure. Gene flow was higher among rivers with similar geologies, suggesting local adaptation to geological substrate. The effect of river length was mainly owing to one highly differentiated population that has the farthest spawning grounds off the river mouth (up to 900km) and the largest fish, suggesting local adaptation to river length. We detected high levels of admixture in stocked populations but also in neighbouring ones, implying large-scale impacts of stocking through dispersal of non-native individuals. However, we found relatively few admixed individuals suggesting a lower fitness of stocked fish and/or some reproductive isolation between wild and stocked individuals. When excluding stocked populations, genetic structure increased as did its correlation with environmental factors. This study overall indicates that geological substrate and river length are major environmental factors influencing gene flow and potential local adaptation among Atlantic salmon populations but that stocking with non-native individuals may ultimately disrupt these natural patterns of gene flow among locally adapted populations.

Effects of different mating scenarios on embryo viability in brown trout.

Mating with attractive or dominant males is often predicted to offer indirect genetic benefits to females, but it is still largely unclear how important such non-random mating can be with regard to embryo viability. We sampled a natural population of adult migratory brown trout (Salmo trutta), bred them in vitro in a half-sib breeding design to separate genetic from maternal environmental effects, raised 2098 embryos singly until hatching, and exposed them experimentally to different levels of pathogen stress at a late embryonic stage. We found that the embryos' tolerance to the induced pathogen stress was linked to the major histocompatibility complex (MHC) of their parents, i.e. certain MHC genotypes appeared to provide better protection against infection than others. We also found significant additive genetic variance for stress tolerance. Melanin-based dark skin patterns revealed males with 'good genes', i.e. embryos fathered by dark coloured males had a high tolerance to infection. Mating with large and dominant males would, however, not improve embryo viability when compared to random mating. We used simulations to provide estimates of how mate choice based on MHC or melanin-based skin patterns would influence embryos' tolerance to the experimentally induced pathogen stress.

Cytidine deaminase 2 is required for VLRB antibody gene assembly in lampreys.

The antibodies of jawless vertebrates consist of leucine-rich repeat arrays encoded by somatically assembled VLRB genes. It is unknown how the incomplete germline VLRB loci are converted into functional antibody genes during B lymphocyte development in lampreys. In Lampetra planeri larvae lacking the cytidine deaminase CDA2 gene, VLRB assembly fails, whereas the T lineage-associated VLRA and VLRC antigen receptor gene assemblies occur normally. Thus, CDA2 acts in a B cell lineage-specific fashion to support the somatic diversification of VLRB antibody genes. CDA2 is closely related to activation-induced cytidine deaminase (AID), which is essential for the elaboration of immunoglobulin gene repertoires in jawed vertebrates. Our results thus identify a convergent mechanism of antigen receptor gene assembly and diversification that independently evolved in the two sister branches of vertebrates.

Parallel changes in genetic diversity and species diversity following a natural disturbance.

We examined the spatial and temporal variation of species diversity and genetic diversity in a metacommunity comprising 16 species of freshwater gastropods. We monitored species abundance at five localities of the Ain river floodplain in southeastern France, over a period of four years. Using 190 AFLP loci, we monitored the genetic diversity of Radix balthica, one of the most abundant gastropod species of the metacommunity, twice during that period. An exceptionally intense drought occurred during the last two years and differentially affected the study sites. This allowed us to test the effect of natural disturbances on changes in both genetic and species diversity. Overall, local (alpha) diversity declined as reflected by lower values of gene diversity H(S) and evenness. In parallel, the among-sites (beta) diversity increased at both the genetic (F(ST)) and species (F(STC)) levels. These results suggest that disturbances can lead to similar changes in genetic and community structure through the combined effects of selective and neutral processes.

Viability of brown trout embryos positively linked to melanin-based but negatively to carotenoid-based colours of their fathers.

'Good-genes' models of sexual selection predict significant additive genetic variation for fitness-correlated traits within populations to be revealed by phenotypic traits. To test this prediction, we sampled brown trout (Salmo trutta) from their natural spawning place, analysed their carotenoid-based red and melanin-based dark skin colours and tested whether these colours can be used to predict offspring viability. We produced half-sib families by in vitro fertilization, reared the resulting embryos under standardized conditions, released the hatchlings into a streamlet and identified the surviving juveniles 20 months later with microsatellite markers. Embryo viability was revealed by the sires' dark pigmentation: darker males sired more viable offspring. However, the sires' red coloration correlated negatively with embryo survival. Our study demonstrates that genetic variation for fitness-correlated traits is revealed by male colour traits in our study population, but contrary to predictions from other studies, intense red colours do not signal good genes.

Male dominance linked to size and age, but not to 'good genes' in brown trout (Salmo trutta).

BACKGROUND: Males that are successful in intra-sexual competition are often assumed to be of superior quality. In the mating system of most salmonid species, intensive dominance fights are common and the winners monopolise most mates and sire most offspring. We drew a random sample of mature male brown trout (Salmo trutta) from two wild populations and determined their dominance hierarchy or traits linked to dominance. The fish were then stripped and their sperm was used for in vitro fertilisations in two full-factorial breeding designs. We recorded embryo viability until hatching in both experiments, and juvenile survival during 20 months after release into a natural streamlet in the second experiment. Since offspring of brown trout get only genes from their fathers, we used offspring survival as a quality measure to test (i) whether males differ in their genetic quality, and if so, (ii) whether dominance or traits linked to dominance reveal 'good genes'. RESULTS: We found significant additive genetic variance on embryo survival, i.e. males differed in their genetic quality. Older, heavier and larger males were more successful in intra-sexual selection. However, neither dominance nor dominance indicators like body length, weight or age were significantly linked to genetic quality measured as embryo or juvenile survival. CONCLUSION: We found no evidence that females can improve their offspring's genetic viability by mating with large and dominant males. If there still were advantages of mating with dominant males, they may be linked to non-genetic benefits or to genetic advantages that are context dependent and therefore possibly not revealed under our experimental conditions - even if we found significant additive genetic variation for embryo viability under such conditions.

[Morphological variation of the spermatheca in the garden snail Cantareus aspersus]. / Variation morphologique de la spermathèque chez l'escargot terrestre Cantareus aspersus.

Spermathecal morphology is known to play an important role in postcopulatory sexual selection of many invertebrates. In helicid land snails, the spermatheca is subdivided into tubules, whose number is sometimes subject to a strong inter-individual variation. Significance of this variation for postcopulatory sexual selection is unknown, but it might be related to cryptic female choice. In the present work, we have investigated the fine multi-tubular structure of the sperm storage organ in Cantareus aspersus. We found between 3 and 13 tubules per individual in a single population, which represents a degree of variation rarely observed in helicid land snails.

Reconstructing the demographic history of divergence between European river and brook lampreys using approximate Bayesian computations.

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey (L. planeri) show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.

Understanding admixture patterns in supplemented populations: a case study combining molecular analyses and temporally explicit simulations in Atlantic salmon.

Genetic admixture between wild and introduced populations is a rising concern for the management of endangered species. Here, we use a dual approach based on molecular analyses of samples collected before and after hatchery fish introduction in combination with a simulation study to obtain insight into the mechanisms of admixture in wild populations. Using 17 microsatellites, we genotyped pre- and post-stocking samples from four Atlantic salmon populations supplemented with non-native fish to estimate genetic admixture. We also used individual-based temporally explicit simulations based on realistic demographic and stocking data to predict the extent of admixture. We found a low admixture by hatchery stocks within prestocking samples but moderate to high values in post-stocking samples (from 12% to 60%). The simulation scenarios best fitting the real data suggested a 10-25 times lower survival of stocked fish relative to wild individuals. Simulations also suggested relatively high dispersal rates of stocked and wild fish, which may explain some high levels of admixture in weakly stocked populations and the persistence of indigenous genotypes in heavily stocked populations. This study overall demonstrates that combining genetic analyses with simulations can significantly improve the understanding of admixture mechanisms in wild populations.

Changes in the genetic structure of Atlantic salmon populations over four decades reveal substantial impacts of stocking and potential resiliency.

While the stocking of captive-bred fish has been occurring for decades and has had substantial immediate genetic and evolutionary impacts on wild populations, its long-term consequences have only been weakly investigated. Here, we conducted a spatiotemporal analysis of 1428 Atlantic salmon sampled from 1965 to 2006 in 25 populations throughout France to investigate the influence of stocking on the neutral genetic structure in wild Atlantic salmon (Salmo salar) populations. On the basis of the analysis of 11 microsatellite loci, we found that the overall genetic structure among populations dramatically decreased over the period studied. Admixture rates among populations were highly variable, ranging from a nearly undetectable contribution from donor stocks to total replacement of the native gene pool, suggesting extremely variable impacts of stocking. Depending on population, admixture rates either increased, remained stable, or decreased in samples collected between 1998 and 2006 compared to samples from 1965 to 1987, suggesting either rising, long-lasting or short-term impacts of stocking. We discuss the potential mechanisms contributing to this variability, including the reduced fitness of stocked fish and persistence of wild locally adapted individuals.