Pierre-Alexandre Gagnaire

European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation

The European sea bass (Dicentrarchus labrax) is a temperate-zone euryhaline teleost of prime importance for aquaculture and fisheries. This species is subdivided into two naturally hybridizing lineages, one inhabiting the north-eastern Atlantic Ocean and the other the Mediterranean and Black seas. Here we provide a high-quality chromosome-scale assembly of its genome that shows a high degree of synteny with the more highly derived teleosts. We find expansions of gene families specifically associated with ion and water regulation, highlighting adaptation to variation in salinity. We further generate a genome-wide variation map through RAD-sequencing of Atlantic and Mediterranean populations. We show that variation in local recombination rates strongly influences the genomic landscape of diversity within and differentiation between lineages. Comparing predictions of alternative demographic models to the joint allele-frequency spectrum indicates that genomic islands of differentiation between sea bass lineages were generated by varying rates of introgression across the genome following a period of geographical isolation.

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.

The genetic consequences of spatially varying selection in the panmictic American eel (Anguilla rostrata).

Our understanding of the genetic basis of local adaptation has recently benefited from the increased power to identify functional variants associated with environmental variables at the genome scale. However, it often remains challenging to determine whether locally adaptive alleles are actively maintained at intermediate frequencies by spatially varying selection. Here, we evaluate the extent to which this particular type of balancing selection explains the retention of adaptive genetic variation in the extreme situation of perfect panmixia, using the American eel (Anguilla rostrata) as a model. We first conducted a genome scan between two samples from opposite ends of a latitudinal environmental gradient using 454 sequencing of individually tagged cDNA libraries. Candidate SNPs were then genotyped in 992 individuals from 16 sampling sites at different life stages of the same cohort (including larvae from the Sargasso Sea, glass eels, and 1-year-old individuals) as well as in glass eels of the following cohort. Evidence for spatially varying selection was found at 13 loci showing correlations between allele frequencies and environmental variables across the entire species range. Simulations under a multiple-niche Levene’s model using estimated relative fitness values among genotypes rarely predicted a stable polymorphic equilibrium at these loci. Our results suggest that some genetic-by-environment interactions detected in our study arise during the progress toward fixation of a globally advantageous allele with spatially variable effects on fitness.

Population genetics of the American eel (Anguilla rostrata): FST = 0 and North Atlantic Oscillation effects on demographic fluctuations of a panmictic species

We performed population genetic analyses on the American eel (Anguilla rostrata) with three main objectives. First, we conducted the most comprehensive analysis of neutral genetic population structure to date to revisit the null hypothesis of panmixia in this species. Second, we used this data to provide the first estimates of contemporary effective population size (Ne) and to document temporal variation in effective number of breeders (Nb) in American eel. Third, we tested for statistical associations between temporal variation in the North Atlantic Oscillation (NAO), the effective number of breeders and two indices of recruit abundance. A total of 2142 eels from 32 sampling locations were genotyped with 18 microsatellite loci. All measures of differentiation were essentially zero, and no evidence for significant spatial or temporal genetic differentiation was found. The panmixia hypothesis should thus be accepted for this species. Nb estimates varied by a factor of 23 among 12 cohorts, from 473 to 10 999. The effective population size Ne was estimated at 10 532 (95% CI, 9312–11 752). This study also showed that genetically based demographic indices, namely Nb and allelic richness (Ar), can be used as surrogates for the abundance of breeders and recruits, which were both shown to be positively influenced by variation during high (positive) NAO phases. Thus, long‐term genetic monitoring of American glass eels at several sites along the North American Atlantic coast would represent a powerful and efficient complement to census monitoring to track demographic fluctuations and better understand their causes.

Using neutral, selected, and hitchhiker loci to assess connectivity of marine populations in the genomic era.

Estimating the rate of exchange of individuals among populations is a central concern to evolutionary ecology and its applications to conservation and management. For instance, the efficiency of protected areas in sustaining locally endangered populations and ecosystems depends on reserve network connectivity. The population genetics theory offers a powerful framework for estimating dispersal distances and migration rates from molecular data. In the marine realm, however, decades of molecular studies have met limited success in inferring genetic connectivity, due to the frequent lack of spatial genetic structure in species exhibiting high fecundity and dispersal capabilities. This is especially true within biogeographic regions bounded by well‐known hotspots of genetic differentiation. Here, we provide an overview of the current methods for estimating genetic connectivity using molecular markers and propose several directions for improving existing approaches using large population genomic datasets. We highlight several issues that limit the effectiveness of methods based on neutral markers when there is virtually no genetic differentiation among samples. We then focus on alternative methods based on markers influenced by selection. Although some of these methodologies are still underexplored, our aim was to stimulate new research to test how broadly they are applicable to nonmodel marine species. We argue that the increased ability to apply the concepts of cline analyses will improve dispersal inferences across physical and ecological barriers that reduce connectivity locally. We finally present how neutral markers hitchhiking with selected loci can also provide information about connectivity patterns within apparently well‐mixed biogeographic regions. We contend that one of the most promising applications of population genomics is the use of outlier loci to delineate relevant conservation units and related eco‐geographic features across which connectivity can be measured.

Mapping phenotypic, expression and transmission ratio distortion QTL using RAD markers in the Lake Whitefish (Coregonus clupeaformis)

The evolution of reproductive isolation in an ecological context may involve multiple facets of species divergence on which divergent selection may operate. These include variation in quantitative phenotypic traits, regulation of gene expression, and differential transmission of particular allelic combinations. Thus, an integrative approach to the speciation process involves identifying the genetic basis of these traits, in order to understand how they are affected by divergent selection in nature and how they ultimately contribute to reproductive isolation. In the Lake Whitefish (Coregonus clupeaformis), dwarf and normal species pairs sympatrically occur in several North American postglacial lakes. The limnetic dwarf whitefish distinguishes from its normal benthic relative by numerous life history, behavioural, morphological and gene expression traits, in relation with the exploitation of distinct ecological niches. Here, we have applied the RAD‐Sequencing method to a hybrid backcross family to reconstruct a high‐density genetic linkage map and perform QTL mapping in the Lake Whitefish. The 3061 cM map encompassed 3438 segregating RAD markers distributed over 40 linkage groups, for an average resolution of 0.89 cM. We mapped phenotypic and expression QTL underlying ecologically important traits as well as transmission ratio distortion QTL, and identified genomic regions harbouring clusters of such QTL. A narrow genomic region strongly associated with sex determination was also evidenced. Positional and functional information revealed in this study will be useful in ongoing population genomic studies to illuminate our understanding of the genomic architecture of reproductive isolation between whitefish species pairs.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs.

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.

WITHIN-POPULATION STRUCTURE HIGHLIGHTED BY DIFFERENTIAL INTROGRESSION ACROSS SEMIPERMEABLE BARRIERS TO GENE FLOW IN ANGUILLA MARMORATA

In the marine environment, differential gene exchange between partially reproductively isolated taxa can result in introgression that extends over long distances due to high larval dispersal potential. However, the degree to which this process contributes to interlocus variance of genetic differentiation within introgressed populations remains unclear. Using a genome‐scan approach in the Indo‐Pacific eel Anguilla marmorata, we investigated the degree of interpopulation genetic differentiation, the rate of introgression, and within‐population genetic patterns at 858 AFLP markers genotyped in 1117 individuals. Three divergent populations were identified based on clustering analysis. Genetic assignments of individuals revealed the existence of different types of hybrids that tended to co‐occur with parental genotypes in three population contact zones. Highly variable levels of genetic differentiation were found between populations across the AFLP markers, and reduced rates of introgression were shown at some highly differentiated loci. Gene flow across semipermeable genetic barriers was shown to generate spatial introgression patterns at some loci which define within‐population structure over long distances. These results suggest that differential introgression in subdivided populations may be relevant when interpreting spatial variation patterns displayed by outlying loci in other marine fish populations.

Nonparallelism in MHCIIβ diversity accompanies nonparallelism in pathogen infection of lake whitefish (Coregonus clupeaformis) species pairs as revealed by next-generation sequencing.

Major histocompatibility (MHC) immune system genes may evolve in response to pathogens in the environment. Because they also may affect mate choice, they are candidates for having great importance in ecological speciation. Here, we use next‐generation sequencing to test the general hypothesis of parallelism in patterns of MHCIIβ diversity and bacterial infections among five dwarf and normal whitefish sympatric pairs. A second objective was to assess the functional relationships between specific MHCIIβ alleles and pathogens in natural conditions. Each individual had between one and four alleles, indicating two paralogous loci. In Cliff Lake, the dwarf ecotype was monomorphic for the most common allele. In Webster Lake, the skew in the allelic distribution was towards the same allele but in the normal ecotype, underscoring the nonparallel divergence among lakes. Our signal of balancing selection matched putative peptide binding region residues in some cases, but not in others, supporting other recent findings of substantial functional differences in fish MHCIIβ compared with mammals. Individuals with fewer alleles were less likely to be infected; thus, we found no evidence for the heterozygote advantage hypothesis. MHCIIβ alleles and pathogenic bacteria formed distinct clusters in multivariate analyses, and clusters of certain alleles were associated with clusters of pathogens, or sometimes the absence of pathogens, indicating functional relationships at the individual level. Given that patterns of MHCIIβ and bacteria were nonparallel among dwarf and normal whitefish pairs, we conclude that pathogens driving MHCIIβ evolution did not play a direct role in their parallel phenotypic evolution.

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.