Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse.

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna.

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Subtle limits to connectivity revealed by outlier loci within two divergent metapopulations of the deep-sea hydrothermal gastropod Ifremeria nautilei.

Hydrothermal vents form archipelagos of ephemeral deep-sea habitats that raise interesting questions about the evolution and dynamics of the associated endemic fauna, constantly subject to extinction-recolonization processes. These metal-rich environments are coveted for the mineral resources they harbour, thus raising recent conservation concerns. The evolutionary fate and demographic resilience of hydrothermal species strongly depend on the degree of connectivity among and within their fragmented metapopulations. In the deep sea, however, assessing connectivity is difficult and usually requires indirect genetic approaches. Improved detection of fine-scale genetic connectivity is now possible based on genome-wide screening for genetic differentiation. Here, we explored population connectivity in the hydrothermal vent snail Ifremeria nautilei across its species range encompassing five distinct back-arc basins in the Southwest Pacific. The global analysis, based on 10,570 single nucleotide polymorphism (SNP) markers derived from double digest restriction-site associated DNA sequencing (ddRAD-seq), depicted two semi-isolated and homogeneous genetic clusters. Demogenetic modeling suggests that these two groups began to diverge about 70,000 generations ago, but continue to exhibit weak and slightly asymmetrical gene flow. Furthermore, a careful analysis of outlier loci showed subtle limitations to connectivity between neighbouring basins within both groups. This finding indicates that migration is not strong enough to totally counterbalance drift or local selection, hence questioning the potential for demographic resilience at this latter geographical scale. These results illustrate the potential of large genomic data sets to understand fine-scale connectivity patterns in hydrothermal vents and the deep sea.

Systematics of European coastal anchovies (genus Engraulis Cuvier).

Reports of morphological differences between European anchovy (Engraulis cf. encrasicolus) from coastal and marine habitats have long existed in the ichthyologic literature and have given rise to a long-standing debate on their taxonomic status. More recently, molecular studies have confirmed the existence of genetic differentiation between the two anchovy ecotypes. Using ancestry-informative markers, we show that coastal anchovies throughout the Mediterranean share a common ancestry and that substantial genetic differentiation persists in different pairs of coastal/marine populations despite the presence of limited gene flow. On the basis of genetic and ecological arguments, we propose that coastal anchovies deserve a species status of their own (E. maeoticus) and argue that a unified taxonomical framework is critical for future research and management.

Spatial gradients of introgressed ancestry reveal cryptic connectivity patterns in a high gene flow marine fish.

Assessing genetic connectivity among populations in high gene flow species is sometimes insufficient to evaluate demographic connectivity. Genetic differentiation quickly becomes zero as soon as a few dozen migrants are exchanged per generation. This provides little information to determine whether migration can ensure demographic coupling. The resulting difficulties in delineating conservation units for the management of commercially exploited marine fish species are well illustrated in the case of the European sea bass (Dicentrarchus labrax). Previous attempts to assess connectivity patterns in the northeast Atlantic have been hampered by a lack of spatial genetic structure. In contrast, mark-recapture data suggested low migration rates between regional spawning areas. Here, we show how a spatial gradient of introgressed Mediterranean ancestry across the northeast Atlantic reflects cryptic patterns of genetic and demographic connectivity. Using a 1K SNP chip data set in 827 individuals sampled from Portugal to the North Sea, we found null overall genetic differentiation across the northeast Atlantic. We however detected a subtle latitudinal admixture gradient originating at the edge of the contact zone with the Mediterranean sea bass lineage. Two significant breaks in the ancestry gradient at the tip of Galicia and northern Brittany indicated barriers to effective dispersal between demographically distinct units. Moreover, a northward expansion signal in Irish and North Seas was revealed by the surfing of rare Mediterranean alleles at the edge of the species range. Our results show that introgressed ancestry gradients offer a powerful alternative to assess genetic and demographic connectivity when the neutral migration-drift balance is not informative.

Speciation history of European (Anguilla anguilla) and American eel (A. rostrata), analysed using genomic data.

Speciation in the ocean could differ from terrestrial environments due to fewer barriers to gene flow. Hence, sympatric speciation might be common, with American and European eel being candidates for exemplifying this. They show disjunct continental distributions on both sides of the Atlantic, but spawn in overlapping regions of the Sargasso Sea from where juveniles are advected to North American, European and North African coasts. Hybridization and introgression are known to occur, with hybrids almost exclusively observed in Iceland. Different speciation scenarios have been suggested, involving either vicariance or sympatric ecological speciation. Using RAD sequencing and whole-genome sequencing data from parental species and F1 hybrids, we analysed speciation history based on the joint allele frequency spectrum (JAFS) and pairwise sequentially Markovian coalescent (PSMC) plots. JAFS supported a model involving a split without gene flow 150,000-160,000 generations ago, followed by secondary contact 87,000-92,000 generations ago, with 64% of the genome experiencing restricted gene flow. This supports vicariance rather than sympatric speciation, likely associated with Pleistocene glaciation cycles and ocean current changes. Whole-genome PSMC analysis of F1 hybrids from Iceland suggested divergence 200,000 generations ago and indicated subsequent gene flow rather than strict isolation. Finally, simulations showed that results from both approaches (JAFS and PSMC) were congruent. Hence, there is strong evidence against sympatric speciation in North Atlantic eels. These results reiterate the need for careful consideration of cases of possible sympatric speciation, as even in seemingly barrier-free oceanic environments palaeoceanographic factors may have promoted vicariance and allopatric speciation.

Development of a large SNPs resource and a low-density SNP array for brown trout (Salmo trutta) population genetics.

BACKGROUND: The brown trout (Salmo trutta) is an economically and ecologically important species for which population genetic monitoring is frequently performed. The most commonly used genetic markers for this species are microsatellites and mitochondrial markers that lack replicability among laboratories, and a large genome coverage. An alternative that may be particularly efficient and universal is the development of small to large panels of Single Nucleotide Polymorphism markers (SNPs). Here, we used Restriction site Associated DNA sequences (RADs) markers to identify a set of 12,204 informative SNPs positioned on the brown trout linkage map and suitable for population genetics studies. Then, we used this novel resource to develop a cost-effective array of 192 SNPs (96 × 2) evenly spread on this map. This array was tested for genotyping success in five independent rivers occupied by two main brown trout evolutionary lineages (Atlantic -AT- and Mediterranean -ME-) on a total of 1862 individuals. Moreover, inference of admixture rate with domestic strains and population differentiation were assessed for a small river system (the Taurion River, 190 individuals) and results were compared to a panel of 13 microsatellites. RESULTS: A high genotyping success was observed for all rivers (< 1% of non-genotyped loci per individual), although some initially used SNP failed to be amplified, probably because of mutations in primers, and were replaced. These SNPs permitted to identify patterns of isolation-by-distance for some rivers. Finally, we found that microsatellite and SNP markers yielded very similar patterns for population differentiation and admixture assessments, with SNPs having better ability to detect introgression and differentiation. CONCLUSIONS: The novel resources provided here opens new perspectives for universality and genome-wide studies in brown trout populations.

Polygenic selection drives the evolution of convergent transcriptomic landscapes across continents within a Nearctic sister species complex.

In contrast to the plethora of studies focusing on the genomic basis of adaptive phenotypic divergence, the role of gene expression during speciation has been much less investigated and consequently less understood. Yet, the convergence of differential gene expression patterns between closely related species-pairs might reflect the role of natural selection during the process of ecological speciation. Here, we test for intercontinental convergence in differential transcriptional signatures between limnetic and benthic sympatric species-pairs of Lake Whitefish (Coregonus clupeaformis) and its sister lineage, the European Whitefish (Coregonus lavaretus), using six replicated sympatric species-pairs (two in North America, two in Norway and two in Switzerland). We characterized both sequence variation in transcribed regions and differential gene expression between sympatric limnetic and benthic species across regions and continents. Our first finding was that differentially expressed genes (DEG) between limnetic and benthic whitefish tend to be enriched in shared polymorphism among sister lineages. We then used both genotypes and covariation in expression in order to infer polygenic selection at the gene level. We identified parallel outliers and DEG involving genes primarily overexpressed in limnetic species relative to the benthic species. Our analysis finally revealed the existence of shared genomic bases underlying parallel differential expression across replicated species-pairs from both continents, such as a cis-eQTL affecting the pyruvate kinase expression level involved in glycolysis. Our results are consistent with a long-standing role of natural selection in maintaining trans-continental diversity at phenotypic traits involved in ecological speciation between limnetic and benthic whitefishes.

The role of recombination on genome-wide patterns of local ancestry exemplified by supplemented brook charr populations.

Assessing the immediate and long-term evolutionary consequences of human-mediated hybridization is of major concern for conservation biology. Several studies have documented how selection in interaction with recombination modulates introgression at a genome-wide scale, but few have considered the dynamics of this process within and among chromosomes. Here, we used an exploited freshwater fish, the brook charr (Salvelinus fontinalis), for which decades of stocking practices have resulted in admixture between wild populations and an introduced domestic strain, to assess both the temporal dynamics and local chromosomal variation in domestic ancestry. We provide a detailed picture of the domestic ancestry patterns across the genome using about 33,000 mapped single nucleotide polymorphisms genotyped in 611 individuals from 24 supplemented populations. For each lake, we distinguished early- and late-generation hybrids using information regarding admixture tracts. To assess the selective outcomes following admixture we then evaluated the relationship between recombination and admixture proportions at three different scales: the whole genome, chromosomes and within 2-Mb windows. This allowed us to detect a wide range of evolutionary mechanisms varying along the genome, as reflected by the finding of favoured or disfavoured introgression of domestic haplotypes. Among these, the main factor modulating local ancestry was probably the presence of deleterious recessive mutations in the wild populations, which can be efficiently hidden to selection in the presence of long admixture tracts. Overall, our results emphasize the relevance of taking into consideration local ancestry information to assess both the temporal and the chromosomal variation in local admixture ancestry toward better understanding post-hybridization evolutionary outcomes.

Model-based demographic inference of introgression history in European whitefish species pairs'.

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi-permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome-wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post-glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.

Population-specific variations of the genetic architecture of sex determination in wild European sea bass Dicentrarchus labrax L.

Polygenic sex determination (PSD) may show variations in terms of genetic and environmental components between populations of fish species exposed/adapted to different environments. The European sea bass (Dicentrarchus labrax) is an interesting model, combining both a PSD system and a genetic subdivision into an Atlantic and a Mediterranean lineage, with genetic substructures within the Mediterranean Sea. Here, we produced experimental progeny crosses (N = 927) from broodstock sampled in four wild populations (North Atlantic, NAT; Western Mediterranean, WEM; North-Eastern Mediterranean, NEM; South-Eastern Mediterranean, SEM). We found less females than males in the progeny, both in the global dataset (32.5%) and within each paternal group (from 25.1% for NEM to 39.0% for WEM), with significant variation among populations, dams, and sires. Sex, body weight (BW), and body length (BL) showed moderate heritability (0.52 ± 0.17, 0.46 ± 0.17, 0.34 ± 0.15, respectively) and sex was genetically correlated with BW and BL, with rAsex/BW = 0.69 ± 0.12 and rA sex/BL = 0.66 ± 0.13. A weighted GWAS performed both on the global dataset and within each paternal group revealed a different genetic architecture of sex determination between Atlantic and Mediterranean populations (9 QTLs found in NAT, 7 in WEM, 5 in NEM, and 4 in SEM, with a cumulated variance explained of 27.04%, 21.87%, 15.89%, and 12.10%, respectively) and a more similar genetic architecture among geographically close populations compared to geographically distant populations, consistent with the hypothesis of a population-specific evolution of polygenic sex determination systems in different environments.

Genomic consequences of a recent three-way admixture in supplemented wild brown trout populations revealed by local ancestry tracts.

Understanding the evolutionary consequences of human-mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make-up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double-digested restriction site-associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple-way admixture with population genomic data.

Genomic and geographic footprints of differential introgression between two divergent fish species (Solea spp.).

Investigating gene flow between closely related species and its variation across the genome is important to understand how reproductive barriers shape genome divergence before speciation is complete. An efficient way to characterize differential gene flow is to study how the genetic interactions that take place in hybrid zones selectively filter gene exchange between species, leading to heterogeneous genome divergence. In the present study, genome-wide divergence and introgression patterns were investigated between two sole species, Solea senegalensis and Solea aegyptiaca, using restriction-associated DNA sequencing (RAD-Seq) to analyze samples taken from a transect spanning the hybrid zone. An integrative approach combining geographic and genomic clines methods with an analysis of individual locus introgression accounting for the demographic history of divergence was conducted. Our results showed that the two sole species have come into secondary contact postglacially, after experiencing a prolonged period (ca. 1.1 to 1.8 Myrs) of allopatric separation. Secondary contact resulted in the formation of a tension zone characterized by strong reproductive isolation, which only allowed introgression in a limited fraction of the genome. We found multiple evidence for a preferential direction of introgression in the S. aegyptiaca genetic background, indicating a possible recent or ongoing movement of the hybrid zone. Deviant introgression signals found in the opposite direction suggested that S. senegalensis could have possibly undergone adaptive introgression that has not yet spread throughout the entire species range. Our study thus illustrates the varied outcomes of genetic interactions between divergent gene pools that recently met after a long history of divergence.

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.

Do North Atlantic eels show parallel patterns of spatially varying selection?

BACKGROUND: The two North Atlantic eel species, the European and the American eel, represent an ideal system in which to study parallel selection patterns due to their sister species status and the presence of ongoing gene flow. A panel of 80 coding-gene SNPs previously analyzed in American eel was used to genotype European eel individuals (glass eels) from 8 sampling locations across the species distribution. We tested for single-generation signatures of spatially varying selection in European eel by searching for elevated genetic differentiation using FST-based outlier tests and by testing for significant associations between allele frequencies and environmental variables. RESULTS: We found signatures of possible selection at a total of 11 coding-gene SNPs. Candidate genes for local selection constituted mainly genes with a major role in metabolism as well as defense genes. Contrary to what has been found for American eel, only 2 SNPs in our study correlated with differences in temperature, which suggests that other explanatory variables may play a role. None of the genes found to be associated with explanatory variables in European eel showed any correlations with environmental factors in the previous study in American eel. CONCLUSIONS: The different signatures of selection between species could be due to distinct selective pressures associated with the much longer larval migration for European eel relative to American eel. The lack of parallel selection in North Atlantic eels could also be due to most phenotypic traits being polygenic, thus reducing the likelihood of selection acting on the same genes in both species.

Comparative genomics reveals adaptive protein evolution and a possible cytonuclear incompatibility between European and American Eels.

During the early stages of speciation, interspecific gene flow may be impeded by deleterious epistatic interactions in hybrids, which maintain parental allelic combinations at the speciation genes. The resulting semipermeable nature of the barrier to interspecific gene flow provides a valuable framework to identify the genes involved in hybrid mortality or sterility, as well as the evolutionary mechanisms that initially caused their divergence. The two Atlantic eels Anguilla anguilla and A. rostrata are partially isolated sister species that naturally hybridize, but whose genetic basis of postzygotic isolation remains unknown. We collected high-throughput sequencing data from the transcriptomes of 58 individuals and discovered 94 genes showing differentially fixed mutations between species. Evidence for positive selection at nuclear diagnostic genes was obtained using multilocus extensions of the McDonald-Kreitman test with polymorphism data from each species. In contrast, mitochondrial protein-coding genes experienced strong purifying selection and mostly diverged at synonymous sites, except for the mt-atp6 gene, which showed an atypically high nonsynonymous to synonymous rate ratio. Nuclear-encoded protein interactors of the mt-atp6 gene in the ATP synthase complex were significantly overrepresented in the list of nuclear diagnostic genes. Further analysis of resequencing data showed that positive selection has operated at both the mt-atp6 gene and its nuclear interactor atp5c1. These findings suggest that a cytonuclear incompatibility caused by a disruption of normal ATP synthase function in hybrids contributes to partial reproductive isolation between European and American eels.

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.

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.

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.

Microsatellite length variation in candidate genes correlates with habitat in the gilthead sea bream Sparus aurata.

The genetic basis and evolutionary implications of local adaptation in high gene flow marine organisms are still poorly understood. In several Mediterranean fish species, alternative migration patterns exist between individuals entering coastal lagoons that offer favourable conditions for growth and those staying in the sea where environmental conditions are less subject to rapid and stressful change. Whether these coexisting strategies are phenotypically plastic or include a role for local adaptation through differential survival needs to be determined. Here, we explore the genetic basis of alternate habitat use in western Mediterranean populations of the gilthead sea bream (Sparus aurata). Samples from lagoonal and open-sea habitats were typed for three candidate gene microsatellite loci, seven anonymous microsatellites and 44 amplified fragment length polymorphism markers to test for genotype-environment associations. While anonymous markers globally indicated high levels of gene flow across geographic locations and habitats, non-neutral differentiation patterns correlated with habitat type were found at two candidate microsatellite loci located in the promoter region of the growth hormone and prolactin genes. Further analysis of these two genes revealed that a mechanism based on habitat choice alone could not explain the distribution of genotype frequencies at a regional scale, thus implying a role for differential survival between habitats. We also found an association between allele size and habitat type, which, in the light of previous studies, suggests that polymorphisms in the proximal promoter region could influence gene expression by modulating transcription factor binding, thus providing a potential explanatory link between genotype and growth phenotype in nature.