Search

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.

Quantifying gene flow between lineages at different stages of the speciation continuum is central to understanding speciation.Heliconiusbutterflies have undergone an adaptive radiation in wing colour patterns driven partly by natural selection for local mimicry. They are also known for assortative mating based on wing pattern. Therefore, wing patterns may facilitate speciation. Here, we take a multifaceted approach to explore speciation and species boundaries between closely-related speciesH. hecaleandH. ismenius. We focus our research in geographic regions where the two are mimetic and contrast this with geographic regions where they do not mimic each other. To examine population history and patterns of gene flow, we developed a four population model accounting for linked selection. This model suggests that the two species have remained isolated for the majority of their history, yet with a small amount of gene exchange. Accordingly, signatures of genomic introgression were small except at a major wing pattern allele. Additionally, behavioural assays suggest that other cues determine strong sexual isolation. Tests for chemical differentiation between species identified major differences in putative pheromones which likely mediate mate choice and the maintenance of species differences. Our results show that strong divergence and adaptive introgression at wing pattern loci may not reveal the wholeHeliconiusspeciation history. In a clade where assortative mating is readily triggered by wing pattern differences, we show that this alone may not lead to speciation and that the accumulation of other barriers to gene flow are important in the completion of speciation.

The Anthropocene threatens worldwide biodiversity, with a potential to induce species decline and range contraction. In this context, it is important to quantify species’ adaptive potential and how demographic changes may impact selection efficacy and the burden of deleterious mutations in different populations. Here we show that key evolutionary processes, including variation in effective population size through postglacial change in demography and recombination rates have affected the efficacy of selection and impacted the load in Coho salmon (Oncorhynchus kisutch), a widely distributed salmonid species on the west coast of North America. Using whole genome resequencing (30x coverage) data from populations at different latitudinal distances from their southern glacial refugium, we found support for postglacial gene surfing, with reducedNeat the range recolonization front, thus inducing both a reduction of the adaptive potential and a surf of deleterious alleles in populations evolving at lowNe. This inference was robust to various proxy of the load, namely ratios of non-synonymous to synonymous diversity, the number of missense and “loss of function” mutations. In addition, comparing residual tetrasomic and re-diploidizing regions of the salmon genome, we found support for a prime role of recombination rates in shaping the within-genome variation of the load. Overall, our empirical results are remarkably consistent with expectations under the nearly neutral theory of molecular evolution. We discuss the fundamental and applied implications of these findings for evolutionary and conservation genomics.

In marine species experiencing intense fishing pressures, knowledge of genetic structure and local adaptation represent a critical information to assist sustainable management. In this study, we performed a landscape genomics analysis in the American lobster to investigate the issues pertaining to the consequences of making use of putative adaptive loci to reliably infer population structure and thus more rigorously delineating biological management units in marine exploited species. Toward this end, we genotyped 14,893 single nucleotide polymorphism (SNPs) in 4,190 lobsters sampled across 96 sampling sites distributed along 1000 km in the northwest Atlantic in both Canada and the USA. As typical for most marine species, we observed a weak, albeit highly significant genetic structure. We also found that adaptive genetic variation allows detecting fine-scale population structure not resolved by neutral genetic variation alone. Using the recent genome assembly of the American lobster, we were able to map and annotate several SNPs located in functional genes potentially implicated in adaptive processes such as thermal stress response, salinity tolerance and growth metabolism pathways. Taken together, our study indicates that weak population structure in high gene flow systems can be resolved at various spatial scales, and that putatively adaptive genetic variation can substantially enhance the delineation of biological management units of marine exploited species.

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter‐ and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated two chromosome-level genome assemblies and sequenced 24 RNA-seq libraries to better annotate the coho salmon genome assemblies. We also resequenced the genomes of 83 coho salmon across their North American range to identify nucleotide variants, characterize the broad effects of isolation-by-distance using a genome-wide association analysis approach, and understand the demographic histories of these salmon by modeling population size from genome-wide data. We observed that more than 13% of all SNPs were associated with latitude (before multiple test correction), likely an affect of isolation-by-distance. From demographic history modeling, we estimated that the SNP latitudinal gradient likely developed as recently as 8,000 years ago. In addition, we identified four genes each harboring multiple SNPs associated with latitude; all of these SNPs were also predicted to modify the function of the gene. Three of these genes have roles in cell junction maintenance and may be involved in osmoregulation. This signifies that ocean salinity may have been a factor influencing coho salmon recolonization after the last glaciation period – generating the current pattern of variation in these three genes.

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 brook lamprey (Lampetra planeri) that has a low ability to pass obstacles to migration. We also tested the hypotheses of i) asymmetric gene flow following water current and ii) admixture with the closely related anadromousL. fluviatilisecotype having a positive effect onL. planerigenetic diversity. We genotyped 2472 individuals, including 225L. fluviatilis, sampled in 81 sites upstream and downstream from barriers to migration, in 29 West-European rivers. Linear modelling revealed a strong positive relationship between the distance to the source and genetic diversity, 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 following water flow. DownstreamL. planeripopulations in sympatry withL. fluviatilisdisplayed consistently higher genetic diversity. We conclude that genetic drift and slight downstream gene flow mainly drive the genetic make up of upstreamL. planeripopulations whereas admixture between ecotypes maintains higher levels of genetic diversity inL. planeripopulations sympatric withL. 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.

Copy number variants (CNVs) are a major component of genotypic and phenotypic variation in genomes. To date, our knowledge of genotypic variation and evolution has largely been acquired by means of single nucleotide polymorphism (SNPs) analyses. Until recently, the adaptive role of structural variants (SVs) and particularly that of CNVs has been overlooked in wild populations, partly due to their challenging identification. Here, we document the usefulness of Rapture, a derived reduced-representation shotgun sequencing approach, to detect and investigate copy number variants (CNVs) alongside SNPs in American lobster (Homarus americanus) populations. We conducted a comparative study to examine the potential role of SNPs and CNVs in local adaptation by sequencing 1,141 lobsters from 21 sampling sites within the southern Gulf of St. Lawrence, which experiences the highest yearly thermal variance of the Canadian marine coastal waters. Our results demonstrated that CNVs accounts for higher genetic differentiation than SNP markers. Contrary to SNPs, for which no significant genetic-environment association was found, 48 CNV candidates were significantly associated with the annual variance of sea surface temperature, leading to the genetic clustering of sampling locations despite their geographic separation. Altogether, we provide a strong empirical case that CNVs putatively contribute to local adaptation in marine species and unveil stronger spatial signal of population structure than SNPs. Our study provides the means to study CNVs in non-model species and highlights the importance of considering structural variants alongside SNPs to enhance our understanding of ecological and evolutionary processes shaping adaptive population structure.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species homing behavior may promote the establishment of local adaptation. We genotyped 7,829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. Results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidates SNP associated with long distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

As climate change threatens species' persistence, predicting the potential for species to adapt to rapidly changing environments is imperative for the development of effective conservation strategies. Eco-evolutionary individual-based models (IBMs) can be useful tools for achieving this objective. We performed a literature review to identify studies that apply these tools in marine systems. Our survey suggested that this is an emerging area of research fuelled in part by developments in modelling frameworks that allow simulation of increasingly complex ecological, genetic and demographic processes. The studies we identified illustrate the promise of this approach and advance our understanding of the capacity for adaptation to outpace climate change. These studies also identify limitations of current models and opportunities for further development. We discuss three main topics that emerged across studies: (i) effects of genetic architecture and non-genetic responses on adaptive potential; (ii) capacity for gene flow to facilitate rapid adaptation; and (iii) impacts of multiple stressors on persistence. Finally, we demonstrate the approach using simple simulations and provide a framework for users to explore eco-evolutionary IBMs as tools for understanding adaptation in changing seas.

Oaks are dominant forest tree species widely distributed across the northern hemisphere, where they constitute natural resources of economic, ecological, social and historical value. Hybridization and adaptive introgression have long been thought to be major evolutionary drivers of the ecological success of oaks. However, the maintenance of species barriers remains a key question, given the extent of interspecific gene flow. We used a demographically-explicit genome scan approach combining approximate Bayesian computation and backward simulations to scan the genomes of four European white oak species for reproductive barriers. We first confirmed earlier findings of interspecific contact since the last glacial maximum. We then made use of the tremendous genetic variation among these four species (31 million SNPs) to identify genomic regions enriched in SNPs deviating from demographic expectations under a hypothesis of recent secondary contact. The outlier-enriched regions detected were of restricted size and widely distributed, over all the chromosomes. We identified 215 regions, containing 227 oak genes. A Literature-based functional annotation of these genes highlighted important functions driving the reproductive isolation of the four species. These functions included tolerance to biotic and abiotic constraints, and intrinsic mating barriers. The functions of the annotated genes were highly consistent with the ecological preferences of the four species. These genes may, therefore, be considered candidate genes in studies of species responses to past and ongoing environmental changes.

Unraveling genetic population structure is challenging in species potentially characterized by large population size and high dispersal rates, often resulting in weak genetic differentiation. Genotyping a large number of samples can improve the detection of subtle genetic structure, but this may substantially increase sequencing cost and downstream bioinformatics computational time. To overcome this challenge, alternative, cost‐effective sequencing approaches, namely Pool‐seq and Rapture, have been developed. We empirically measured the power of resolution and congruence of these two methods in documenting weak population structure in nonmodel species with high gene flow comparatively to a conventional genotyping‐by‐sequencing (GBS) approach. For this, we used the American lobster (Homarus americanus) as a case study. First, we found that GBS, Rapture, and Pool‐seq approaches gave similar allele frequency estimates (i.e., correlation coefficient over 0.90) and all three revealed the same weak pattern of population structure. Yet, Pool‐seq data showed FST estimates three to five times higher than GBS and Rapture, while the latter two methods returned similar FST estimates, indicating that individual‐based approaches provided more congruent results than Pool‐seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost‐effective in project involving a relatively small numbers of genotyped individuals (e.g.

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

High-throughput sequencing technologies have provided an unprecedented opportunity to study the different evolutionary forces that have shaped present-day patterns of genetic diversity, with important implications for many directions in plant biology research. To manage such massive quantities of sequencing data, biologists, however, need new additional skills in informatics and statistics. In this chapter, our objective is to introduce population genomics methods to beginners following a learning-by-doing strategy in order to help the reader to analyze the sequencing data by themselves. Conducted analyses cover several main areas of evolutionary biology, such as an initial description of the evolutionary history of a given species or the identification of genes targeted by natural or artificial selection. In addition to the practical advices, we performed re-analyses of two cases studies with different kind of data: a domesticated cereal (African rice) and a non-domesticated tree species (sessile oak). All the code needed to replicate this work is publicly available on github ( https://github.com/ThibaultLeroyFr/Intro2PopGenomics/ ).

High-throughput sequencing technologies have provided an unprecedented opportunity to study the different evolutionary forces that have shaped present-day patterns of genetic diversity, with important implications for many directions in plant biology research. To manage such massive quantities of sequencing data, biologists, however, need new additional skills in informatics and statistics. In this chapter, our objective is to introduce population genomics methods to beginners following a learning-by-doing strategy in order to help the reader to analyze the sequencing data by themselves. Conducted analyses cover several main areas of evolutionary biology, such as an initial description of the evolutionary history of a given species or the identification of genes targeted by natural or artificial selection. In addition to the practical advices, we performed re-analyses of two cases studies with different kind of data: a domesticated cereal (African rice) and a non-domesticated tree species (sessile oak). All the code needed to replicate this work is publicly available on github ( https://github.com/ThibaultLeroyFr/Intro2PopGenomics/ ).

International audience; Dispersal is a central process in ecology and evolution. At the individual level, the three stages of the dispersal process (i.e., emigration, transience and immigration) are affected by complex interactions between phenotypes and environmental factors. Condition‐ and context‐dependent dispersal have far‐reaching consequences, both for the demography and the genetic structuring of natural populations and for adaptive processes. From an applied point of view, dispersal also deeply affects the spatial dynamics of populations and their ability to respond to land‐use changes, habitat degradation and climate change. For these reasons, dispersal has received considerable attention from ecologists and evolutionary biologists. Demographic and genetic methods allow quantifying non‐effective (i.e., followed or not by a successful reproduction) and effective (i.e., with a successful reproduction) dispersal and to investigate how individual and environmental factors affect the different stages of the dispersal process. Over the past decade, demographic and genetic methods designed to quantify dispersal have rapidly evolved but interactions between researchers from the two fields are limited. We here review recent developments in both demographic and genetic methods to study dispersal in wild animal populations. We present their strengths and limits, as well as their applicability depending on study objectives and population characteristics. We propose a unified framework allowing researchers to combine methods and select the more suitable tools to address a broad range of important topics about the ecology and evolution of dispersal and its consequences on animal population dynamics and genetics.

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. Here, we investigated the demographic history of Atlantic salmon across the entire species range using 2035 anadromous individuals from North America and Eurasia. By combining results from admixture graphs, geo-genetic maps, and an Approximate Bayesian Computation (ABC) framework, we validated previous hypotheses pertaining to secondary contact between European and Northern American populations, but also identified secondary contacts in European populations from different glacial refugia. We further identified the major sources of admixture from the southern range of North America into more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflect the spatial redistribution of ancestral variation across the entire North American range. Results also support a role for linked selection and differential introgression that likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere. We conclude that studies between partially isolated populations should systematically include heterogeneity in selective and introgressive effects among loci to perform more rigorous demographic inferences of the divergence process.

International audience; Many angiosperms prevent inbreeding through a self-incompatibility (SI) system, but the loss of SI has been frequent in their evolutionary history. The loss of SI may often lead to an increase in the selfing rate, with the purging of inbreeding depression and the ultimate evolution of a selfing syndrome, where plants have smaller flowers with reduced pollen and nectar production. In this study, we used approximate Bayesian computation (ABC) to estimate the timing of divergence between populations of the plant Linaria cavanillesii that differ in SI status and in which SI is associated with low inbreeding depression but not with a transition to full selfing or a selfing syndrome. Our analysis suggests that the mixed-mating self-compatible (SC) population may have begun to diverge from the SI populations around 2810 generation ago, a period perhaps too short for the evolution of a selfing syndrome. We conjecture that the SC population of L. cavanillesii is at an intermediate stage of transition between outcrossing and selfing.

Population genetic theory states that adaptation most frequently occurs from standing genetic variation, which results from the interplay between different evolutionary processes including mutation, chromosomal rearrangements, drift, gene flow and selection. To date, empirical work focusing on the contribution of standing genetic variation to local adaptation in the presence of high gene flow has been limited to a restricted number of study systems. Marine organisms are excellent biological models to address this issue since many species have to cope with variable environmental conditions acting as selective agents despite high dispersal abilities. In this study, we examined how, demographic history, standing genetic variation linked to chromosomal rearrangements and shared polymorphism among glacial lineages contribute to local adaptation to environmental conditions in the marine fish, the capelin (Mallotus villosus). We used a comprehensive dataset of genome-wide single nucleotide polymorphisms (25,904 filtered SNPs) genotyped in 1,359 individuals collected from 31 spawning sites in the northwest Atlantic (North America and Greenland waters). First, we reconstructed the history of divergence among three glacial lineages and showed that they diverged from 3.8 to 1.8 MyA. Depending on the pair of lineages considered, historical demographic modelling provided evidence for divergence with gene flow and secondary contacts, shaped by barriers to gene flow and linked selection. We next identified candidate loci associated with reproductive isolation of these lineages. Given the absence of physical or geographic barriers, we thus propose that these lineages may represent three cryptic species of capelin. Within each of these, our analyses provided evidence for largeNeand high gene flow at both historical and contemporary time scales among spawning sites. Furthermore, we detected a polymorphic chromosomal rearrangement leading to the coexistence of three haplogroups within the Northwest Atlantic lineage, but absent in the other two clades. Genotype-environment associations revealed molecular signatures of local adaptation to environmental conditions prevailing at spawning sites. Altogether, our study shows that standing genetic variation associated with both chromosomal rearrangements and ancestral polymorphism contribute to local adaptation in the presence of high gene flow.

Gene flow has tremendous importance for local adaptation, by influencing the fate of de novo mutations, maintaining standing genetic variation and driving adaptive introgression. Furthermore, structural variation as chromosomal rearrangements may facilitate adaptation despite high gene flow. However, our understanding of the evolutionary mechanisms impending or favouring local adaptation in the presence of gene flow is still limited to a restricted number of study systems. In this study, we examined how demographic history, shared ancestral polymorphism, and gene flow among glacial lineages contribute to local adaptation to sea conditions in a marine fish, the capelin (Mallotus villosus). We first assembled a 490-Mbp draft genome of M. villosus to map our RAD sequence reads. Then, we used a large data set of genome-wide single nucleotide polymorphisms (25,904 filtered SNPs) genotyped in 1,310 individuals collected from 31 spawning sites in the northwest Atlantic. We reconstructed the history of divergence among three glacial lineages and showed that they probably diverged from 3.8 to 1.8 million years ago and experienced secondary contacts. Within each lineage, our analyses provided evidence for large Ne and high gene flow among spawning sites. Within the Northwest Atlantic lineage, we detected a polymorphic chromosomal rearrangement leading to the occurrence of three haplogroups. Genotype-environment associations revealed molecular signatures of local adaptation to environmental conditions prevailing at spawning sites. Our study also suggests that both shared polymorphisms among lineages, resulting from standing genetic variation or introgression, and chromosomal rearrangements may contribute to local adaptation in the presence of high gene flow.

A thorough reconstruction of historical processes is essential for a comprehensive understanding the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded in postglacial time, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

International audience; The particular combinations of alleles that define haplotypes along individual chromosomes can be determined with increasing ease and accuracy by using current sequencing technologies. Beyond allele frequencies, haplotype data collected in population samples contain information about the history of allelic associations in gene genealogies, and this is of tremendous potential for conservation genomics. We provide an overview of how haplotype information can be used to assess historical demography, gene flow, selection, and the evolutionary outcomes of hybridization across different timescales relevant to conservation issues. We address technical aspects of applying such approaches to nonmodel species. We conclude that there is much to be gained by integrating haplotype-based analyses in future conservation genomics studies.

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. Here, we investigated the demographic history of Atlantic salmon across the entire species range using 2035 anadromous individuals from North America and Eurasia. By combining results from admixture graphs, geo-genetic maps, and an Approximate Bayesian Computation (ABC) framework, we validated previous hypotheses pertaining to secondary contact between European and Northern American populations, but also identified secondary contacts in European populations from different glacial refugia. We further identified the major sources of admixture from the southern range of North America into more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflect the spatial redistribution of ancestral variation across the entire North American range. Results also support a role for linked selection and differential introgression that likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere. We conclude that studies between partially isolated populations should systematically include heterogeneity in selective and introgressive effects among loci to perform more rigorous demographic inferences of the divergence process.

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. In the Northern hemisphere in particular, species genetic makeup has been largely influenced by severe climatic oscillations of the Quaternary Period. Here, we investigated the demographic history of Atlantic Salmon across the entire species range using 2035 anadromous individuals from 77 sampling sites from North America and Eurasia genotyped at 4,656 SNPs. By combining results from admixture graphs, geogenetic maps and an approximate Bayesian computation framework, we validate previous hypotheses pertaining to secondary contact between European and Northern American populations, but also demonstrate that European populations from different glacial refugia have been exchanging alleles in contemporary times. We further identify the major sources of admixture from the southern range of North America to more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflects the spatial redistribution of ancestral variation across the entire American range. Results also point to a role for linked selection in the form of background selection and or positive hitchhiking. Altogether, differential introgression and linked selective effects likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere Therefore we conclude that such heterogeneity among loci should be systematically integrated into demographic inferences of the divergence process, even between incompletely reproductively isolated populations.

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.

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.

International audience; Left-right asymmetries in the epithalamic region of the brain are widespread across vertebrates, but their magnitude and laterality varies among species. Whether these differences reflect independent origins of forebrain asymmetries or taxa-specific diversifications of an ancient vertebrate feature remains unknown. Here we show that the catshark Scyliorhinus canicula and the lampreys Petromyzon marinus and Lampetra planeri exhibit conserved molecular asymmetries between the left and right developing habenulae. Long-term pharmacological treatments in these species show that nodal signalling is essential to their generation, rather than their directionality as in teleosts. Moreover, in contrast to zebrafish, habenular left-right differences are observed in the absence of overt asymmetry of the adjacent pineal field. These data support an ancient origin of epithalamic asymmetry, and suggest that a nodal-dependent asymmetry programme operated in the forebrain of ancestral vertebrates before evolving into a variable trait in bony fish.

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.

cited By 4; International audience; Understanding dispersal capacities for migratory species is crucial for their management. By coupling otolith microchemistry and microsatellite genetic analyses, we provided information on snapshot and long-term dispersal capacity of Allis shad (Alosa alosa), an anadromous clupeid in decline throughout its distribution range. The allocation of natal origin was obtained from water chemistry and signatures in otoliths of juveniles and spawners within a Bayesian model. The majority of adults were assigned to a source river with high degrees of confidence; only 4% were undetermined. Otolith natal origins were used to define a population baseline by grouping individuals from the same natal river and not from the same sampling location as usually done. While A. alosa exhibited a high level of natal site fidelity, this species showed weak genetic structure, which supported the evidence of a substantial flow of strayers between river basins in the vicinity or at longer distances. However, long distance straying was probable but not frequent. In a context of global change, straying would be a key mechanism to drive dispersal and allow resilience of Allis shad populations.

ISI Document Delivery No.: CB8UWTimes Cited: 0Cited Reference Count: 5Cited References: Aprahamian MW, 2002, R & D Technical Report W1-014/TR,Rougemont, Quentin Besnard, Anne-Laure Bagliniere, Jean-Luc Launey, SophieEuropean Regional Development Fund; French Ministry of Ecology and Sustainable DevelopmentWe thank our colleagues from the CRESCO laboratory (MNHN Dinard) and I. Lebel and Y. Abdallah from Migrateurs Rhone Mediterrannee for their help in collecting the samples. This study received funding from the European Regional Development Fund (Transnational program Interreg IV, Atlantic Aquatic Resource Conservation Project) and the French Ministry of Ecology and Sustainable Development (project 'Amphihalins, coherence du reseau N2000 en mer').SpringerDordrechtBiodiversity & Conservation; Genetics & Heredity; Allis shad Alosa alosa and twaite shad Alosa fallax are two related anadromous European clupeid species which have undergone a sharp decline in the past decades. We describe the development of 13 microsatellite markers for both species, based on 454 pyrosequencing technology. The number of alleles per locus ranged from 2 to 11 in A. alosa and from 4 to 13 in A. fallax. Gene diversity ranged from 0.364 to 0.800 for A. alosa and 0.226 to 0.849 for A. fallax. These markers will help in conservation genetics studies such as assessing the extent of population decline, defining appropriate conservation units and monitoring reintroduction programs.