ABSTRACT
Steeper and sometimes extreme salinity gradients increasingly affect aquatic organisms because of climate change. Hypersalinity habitats demand powerful physiological adaptive strategies. Few teleost species have the capacity to spend their whole life cycle in salinities way over seawater levels. Focusing on the multifunctional gill, we unraveled the tilapia S. melanotheron key strategies to cope with different environmental conditions, ranging from freshwater up to hypersaline habitats. De novo transcriptome assembly based on RNAseq allowed for the analysis of 40,967 annotated transcripts among samples collected in three wild populations at 0, 40 and 80 . A trend analysis of the expression patterns revealed responses across the salinity gradient with different gene pathways involved. Genes linked to ion transport, pH regulation and cell surface receptor signaling were mainly upregulated in the high salinity habitat. We identified tight junction proteins that were critical in high salinity habitats and that were different from the well-known tightening junctional proteins identified and expressed in fresh water. Expression profiles also suggest a change in the vascular tone that could be linked to an osmorespiratory compromise not only in fresh water, but also in high salinity environments. A striking downregulation of genes linked to the immune system and to the heat shock response was observed suggesting an energetic trade-off between immunity and acclimation/adaptation in the hypersaline habitat. The high expression of transcripts coding for immune and heat shock response in the freshwater habitat suggests the establishment of powerful mechanisms to protect gills from environmental threats and to maintain protein integrity. Non-directional expression trends were also detected with an upregulation of genes only in the hypersaline habitat (80 ) or only in the marine habitat (40 ). Unravel physiological strategies in S. melanotheron populations will help to better understand the molecular basis of fish euryhalinity in salinity-contrasted environments.
Subject(s)
Gene Expression Profiling , Gills , Salinity , Tilapia , Transcriptome , Animals , Gills/metabolism , Tilapia/genetics , Tilapia/physiology , SeawaterABSTRACT
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.
ABSTRACT
The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.
Subject(s)
Genome-Wide Association Study , Oncorhynchus mykiss , Animals , Oncorhynchus mykiss/genetics , Genome , Skin Pigmentation/geneticsABSTRACT
Viral nervous necrosis (VNN) is an infectious disease caused by the red-spotted grouper nervous necrosis virus (RGNNV) in European sea bass and is considered a serious concern for the aquaculture industry with fry and juveniles being highly susceptible. To understand the genetic basis for resistance against VNN, a survival phenotype through the challenge test against the RGNNV was recorded in populations from multiple year classes (YC2016 and YC2017). A total of 4,851 individuals from 181 families were tested, and a subset (nâ¼1,535) belonging to 122 families was genotyped using a â¼57K Affymetrix Axiom array. The survival against the RGNNV showed low to moderate heritability with observed scale estimates of 0.18 and 0.25 obtained using pedigree vs. genomic information, respectively. The genome-wide association analysis showed a strong signal of quantitative trait loci (QTL) at LG12 which explained â¼33% of the genetic variance. The QTL region contained multiple genes (ITPK1, PLK4, HSPA4L, REEP1, CHMP2, MRPL35, and SCUBE) with HSPA4L and/or REEP1 genes being highly relevant with a likely effect on host response in managing disease-associated symptoms. The results on the accuracy of predicting breeding values presented 20-43% advantage in accuracy using genomic over pedigree-based information which varied across model types and applied validation schemes.
ABSTRACT
BACKGROUND: In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from 'omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments. RESULTS: We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs. CONCLUSION: Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.
Subject(s)
Bass , Animals , Bass/genetics , Epigenomics , Erythrocytes , Genetic Linkage , ZebrafishABSTRACT
Understanding the genetic underpinnings of fitness trade-offs across spatially variable environments remains a major challenge in evolutionary biology. In Mediterranean gilthead sea bream, first-year juveniles use various marine and brackish lagoon nursery habitats characterized by a trade-off between food availability and environmental disturbance. Phenotypic differences among juveniles foraging in different habitats rapidly appear after larval settlement, but the relative role of local selection and plasticity in phenotypic variation remains unclear. Here, we combine phenotypic and genetic data to address this question. We first report correlations of opposite signs between growth and condition depending on juvenile habitat type. Then, we use single nucleotide polymorphism (SNP) data obtained by Restriction Associated DNA (RAD) sequencing to search for allele frequency changes caused by a single generation of spatially varying selection between habitats. We found evidence for moderate selection operating at multiple loci showing subtle allele frequency shifts between groups of marine and brackish juveniles. We identified subsets of candidate outlier SNPs that, in interaction with habitat type, additively explain up to 3.8% of the variance in juvenile growth and 8.7% in juvenile condition; these SNPs also explained significant fraction of growth rate in an independent larval sample. Our results indicate that selective mortality across environments during early-life stages involves complex trade-offs between alternative growth strategies.
Subject(s)
Gene-Environment Interaction , Genetic Fitness/genetics , Sea Bream/genetics , Selection, Genetic/genetics , Animals , Ecosystem , Environment , Gene Frequency , Multifactorial Inheritance/genetics , Polymorphism, Single Nucleotide/genetics , Quantitative Trait Loci/geneticsABSTRACT
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.
Subject(s)
Polymorphism, Single Nucleotide , Trout/genetics , Animals , Genetics, Population , Microsatellite RepeatsABSTRACT
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.
Subject(s)
Genetics, Population , Hybridization, Genetic , Trout/genetics , Animals , Animals, Wild , Aquaculture , Conservation of Natural Resources , France , Gene Frequency , Genotype , Polymorphism, Single NucleotideABSTRACT
We examined specimens of the macrostigma trout Salmo macrostigma, which refers to big black spots on the flanks, to assess whether it is an example of taxonomic inflation within the brown trout Salmo trutta complex. Using new specimens, publicly available data and a mitogenomic protocol to amplify the control and cytochrome b regions of the mitochondrial genome from degraded museum samples, including one syntype specimen, the present study shows that the macrostigma trout is not a valid species. Our results suggest the occurrence of a distinct evolutionary lineage of S. trutta in North Africa and Sicily. The name of the North African lineage is proposed for this lineage, which was found to be sister to the Atlantic lineage of brown trout, S. trutta.
Subject(s)
Phylogeny , Salmonidae/classification , Animals , Biological Evolution , Salmonidae/genetics , Trout/geneticsABSTRACT
Evaluating biodiversity and understanding the processes involved in diversification are noticeable conservation issues in fishes subject to large, sometimes illegal, ornamental trade purposes. Here, the diversity and evolutionary history of the Neotropical dwarf cichlid genus Apistogramma from several South American countries are investigated. Mitochondrial and nuclear markers are used to infer phylogenetic relationships between 31 genetically identified species. The monophyly of Apistogramma is suggested, and Apistogramma species are distributed into four clades, corresponding to three morphological lineages. Divergence times estimated with the Yule process and an uncorrelated lognormal clock dated the Apistogramma origin to the beginning of the Eocene (≈ 50 Myr) suggesting that diversification might be related to marine incursions. Our molecular dating also suggests that the Quaternary glacial cycles coincide with the phases leading to Apistogramma speciation. These past events did not influence diversification rates in the speciose genus Apistogramma, since diversification appeared low and constant through time. Further characterization of processes involved in recent Apistogramma diversity will be necessary.
Subject(s)
Biodiversity , Cichlids/physiology , Animals , Cytochromes b/genetics , Electron Transport Complex IV/genetics , Haplotypes/genetics , Likelihood Functions , Phylogeny , Polymerase Chain Reaction , South America , Species Specificity , Time FactorsABSTRACT
Describing and explaining the geographic within-species variation in phenotypes ("phenogeography") in the sea over a species distribution range is central to our understanding of a variety of eco-evolutionary topics. However, phenogeographic studies that have a large potential to investigate adaptive variation are overcome by phylogeographic studies, still mainly focusing on neutral markers. How genotypic and phenotypic data could covary over large geographic scales remains poorly understood in marine species. We crossed 75 noninbred sires (five origins) and 26 dams (two origins; each side of a hybrid zone) in a factorial diallel cross in order to investigate geographic variation for early survival and sex ratio in the metapopulation of the European sea bass (Dicentrarchus labrax), a highly prized marine fish species. Full-sib families (N = 1,950) were produced and reared in a common environment. Parentage assignment of 7,200 individuals was performed with seven microsatellite markers. Generalized linear models showed significant additive effects for both traits and pleiotropy between traits. A significant nonadditive genetic effect was detected. Different expression of traits and distinct relative performances were found for reciprocal crosses involving populations located on each side of the main hybrid zone located at the Almeria-Oran front, illustrating asymmetric reproductive isolation. The poor fitness performance observed for the Western Mediterranean population of sea bass is discussed as it represents the main source of seed hatchery production, but also because it potentially illustrates nonadaptive introgression and maladaptation.
ABSTRACT
High-density linkage maps are valuable tools for conservation and eco-evolutionary issues. In salmonids, a complex rediploidization process consecutive to an ancient whole genome duplication event makes linkage maps of prime importance for investigating the evolutionary history of chromosome rearrangements. Here, we developed a high-density consensus linkage map for the brown trout (Salmo trutta), a socioeconomically important species heavily impacted by human activities. A total of 3977 ddRAD markers were mapped and ordered in 40 linkage groups using sex- and lineage-averaged recombination distances obtained from two family crosses. Performing map comparison between S. trutta and its sister species, S. salar, revealed extensive chromosomal rearrangements. Strikingly, all of the fusion and fission events that occurred after the S. salar/S. trutta speciation happened in the Atlantic salmon branch, whereas the brown trout remained closer to the ancestral chromosome structure. Using the strongly conserved synteny within chromosome arms, we aligned the brown trout linkage map to the Atlantic salmon genome sequence to estimate the local recombination rate in S. trutta at 3721 loci. A significant positive correlation between recombination rate and within-population nucleotide diversity (π) was found, indicating that selection constrains variation at linked neutral sites in brown trout. This new high-density linkage map provides a useful genomic resource for future aquaculture, conservation, and eco-evolutionary studies in brown trout.
Subject(s)
Chromosome Mapping , Chromosomes/genetics , Gene Rearrangement/genetics , Genetic Variation , Selection, Genetic , Trout/genetics , Animals , Centromere/metabolism , Genome , Recombination, Genetic/genetics , Synteny/geneticsABSTRACT
The black-chinned tilapia Sarotherodon melanotheron heudelotii Rüppell 1852 (Teleostei, Cichlidae) displays remarkable acclimation capacities. When exposed to drastic changes of salinity, which can be the case in its natural habitat, it develops quick physiological responses and keeps reproducing. The present study focused on the physiological impact of salinity on male reproductive capacities, using gene expression as a proxy of acclimation process. Two series of experimental fish were investigated: the first one was composed of fish maintained in freshwater for several generations and newly acclimated to salinities of 35 and 70, whereas the second one consisted of the descendants of the latter born and were raised under their native salinity. Expression patterns of 43 candidate genes previously identified from the testes of wild males was investigated in the three salinities and two generations. Twenty of them showed significant expression differences between salinities, and their predicted function revealed that most of them are involved in the osmotic tolerance of sperm cells and/or in the maintenance of sperm motility. A high level of expression variation was evidenced, especially for fish maintained in freshwater. In spite of this, gene expression patterns allowed the differentiation between fish raised in freshwater and those maintained in hypersaline water in both generations. Altogether, the results presented here suggest that this high variability of expression is likely to ensure the reproductive success of this species under varying salinities.
ABSTRACT
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.
Subject(s)
Adaptation, Physiological , Bass/genetics , Genetic Speciation , Genome , Animals , Atlantic Ocean , Bass/physiology , Chromosome Mapping , Genetic Variation , Molecular Sequence Data , PhylogenyABSTRACT
Although food deprivation is a major ecological pressure in fishes, there is wide individual variation in tolerance of fasting, whose mechanistic bases are poorly understood. Two thousand individually tagged juvenile European sea bass were submitted to two 'fasting/feeding' cycles each comprising 3 weeks of food deprivation followed by 3 weeks of ad libitum feeding at 25°C. Rates of mass loss during the two fasting periods were averaged for each individual to calculate a population mean. Extreme fasting tolerant (FT) and sensitive (FS) phenotypes were identified that were at least one and a half standard deviations, on opposing sides, from this mean. Respirometry was used to investigate two main hypotheses: (1) tolerance of food deprivation reflects lower mass-corrected routine metabolic rate (RMR) in FT phenotypes when fasting, and (2) tolerance reflects differences in substrate utilisation; FT phenotypes use relatively less proteins as metabolic fuels during fasting, measured as their ammonia quotient (AQ), the simultaneous ratio of ammonia excretion to RMR. There was no difference in mean RMR between FT and FS over 7 days fasting, being 6.70±0.24 mmol h(-1) fish(-1) (mean ± s.e.m., N=18) versus 6.76±0.22 mmol h(-1) fish(-1) (N=17), respectively, when corrected to a body mass of 130 g. For any given RMR, however, the FT lost mass at a significantly lower rate than FS, overall 7-day average being 0.72±0.05 versus 0.90±0.05 g day(-1) fish(-1), respectively (P<0.01, t-test). At 20 h after receiving a ration equivalent to 2% body mass as food pellets, ammonia excretion and simultaneous RMR were elevated and similar in FT and FS, with AQs of 0.105±0.009 and 0.089±0.007, respectively. At the end of the period of fasting, ammonia excretion and RMR had fallen in both phenotypes, but AQ was significantly lower in FT than FS, being 0.038±0.004 versus 0.061±0.005, respectively (P<0.001, t-test). There was a direct linear relationship between individual fasted AQ and rate of mass loss, with FT and FS individuals distributed at opposing lower and upper extremities, respectively. Thus the difference between the phenotypes in their tolerance of food deprivation did not depend upon their routine energy use when fasting. Rather, it depended upon their relative use of tissue proteins as metabolic fuels when fasting, which was significantly lower in FT phenotypes.
Subject(s)
Bass/growth & development , Bass/physiology , Food Deprivation/physiology , Animal Nutritional Physiological Phenomena , Animals , Energy Metabolism/physiology , Gastrointestinal Tract/physiology , Linear ModelsABSTRACT
The bonga shad, Ethmalosa fimbriata, is a West African pelagic species still abundant in most habitats of its distribution range and thought to be only recently affected by anthropogenic pressure (habitat destruction or fishing pressure). Its presence in a wide range of coastal habitats characterised by different hydrodynamic processes, represents a case study useful for evaluating the importance of physical structure of the west African shoreline on the genetic structure of a small pelagic species. To investigate this question, the genetic diversity of E. fimbriata was assessed at both regional and species range scales, using mitochondrial (mt) and nuclear DNA markers. Whereas only three panmictic units were identified with mtDNA at the large spatial scale, nuclear genetic markers (EPIC: exon-primed intron-crossing) indicated a more complex genetic pattern at the regional scale. In the northern-most section of shad's distribution range, up to 4 distinct units were identified. Bayesian inference as well as spatial autocorrelation methods provided evidence that gene flow is impeded by the presence of deep-water areas near the coastline (restricting the width of the coastal shelf), such as the Cap Timiris and the Kayar canyons in Mauritania and Senegal, respectively. The added discriminatory power provided by the use of EPIC markers proved to be essential to detect the influence of more subtle, contemporary processes (e.g. gene flow, barriers, etc.) acting within the glacial refuges identified previously by mtDNA.
Subject(s)
Fishes/genetics , Africa, Western , Animal Migration , Animals , Bayes Theorem , DNA/genetics , DNA/isolation & purification , Ecosystem , Gene Flow , Genetic Variation , Genetics, PopulationABSTRACT
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.
Subject(s)
Gene-Environment Interaction , Microsatellite Repeats , Sea Bream/genetics , Alleles , Amplified Fragment Length Polymorphism Analysis , Animals , Ecosystem , Gene Flow , Gene Frequency , Genetic Variation , Genetics, Population , Genotype , Mediterranean Sea , Promoter Regions, Genetic , Sequence Analysis, DNAABSTRACT
The population genetic structure of sea bass (Dicentrarchus labrax) along a transect from the Atlantic Ocean (AO) to the Eastern Mediterranean (EM) Sea differs from that of most other marine taxa in this area. Three populations (AO, Western Mediterranean [WM], EM) are recognized today, which were originally two allopatric populations. How two ancestral genetic units have evolved into three distinct units has not been addressed yet. Therefore, to investigate mechanisms that lead to the emergence of the central WM population, its current status, and its connectivity with the two parental populations, we applied 20 nuclear loci that were either gene associated or gene independent. Results confirmed the existence of three distinct gene pools, with higher differentiation at two transitional areas, the Almeria-Oran Front (AOF) and of the Siculo-Tunisian Strait (STS), than within any population. Significant linkage disequilibrium and heterozygote excess indicated that the STS is probably another tension zone, as already described for the AOF. Neutrality tests fail to reveal marker loci that could be driven by selection within or among metapopulations, except for locus DLA0068. Collectively, results support that the central WM population arose by trapping two tensions zones at distinct geographic locations of limited connectivity. Population assignment further revealed that WM individuals were more introgressed than individuals from the other two metapopulations. This suggests that this population might result from hybrid swarming, and was or is still seeded by genes received through the filter of each tension zone.
ABSTRACT
Growth-related characters, condition factor, and genetic differentiation were investigated for a single cohort of young-of-the-year (YOY) sole within and among nurseries with differing levels of heavy metals (Cd, Cu and Zn) contamination in the two Charentais Straits, Bay of Biscay, France. Analyses were performed when individuals recruited (May), then after a full summer spent in each nursery (October). Levels of phenotypic and genetic diversity were compared, together with genetic differentiation at a candidate metallothionein (MT) locus and three putatively neutral microsatellite loci. No phenotypic or genetic differentiation was detected among nurseries in May, but significant variation at each phenotypic trait and at the multilocus level in October (P < 0.001). Single locus analysis demonstrated that only the MT locus was significantly differentiated among nurseries, whether corrected for null alleles or not (θ=0.0401 and θ(corr.FreeNA)=0.0326, respectively; P < 0.001). Results indicate that phenotypic differences among YOY sole nurseries present a molecular correlate acting at identical spatio-temporal scales among nurseries, potentially reflecting differential selective pressure among nurseries in response to contamination.