Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. ...The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species’ range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species’ natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (FCT = 0.0011; P‐value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean FST = 0.00185; CI: 0.0007–0.0021, P‐value < 0.0002), providing strong evidence for weak, albeit fine‐scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest FST. Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid ‘high‐grading bias’, 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine‐scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.
Investigating how environmental features shape the genetic structure of populations is crucial for understanding how they are potentially adapted to their habitats, as well as for sound management. ...In this study, we assessed the relative importance of spatial distribution, ocean currents and sea surface temperature (SST) on patterns of putatively neutral and adaptive genetic variation among American lobster from 19 locations using population differentiation (PD) approaches combined with environmental association (EA) analyses. First, PD approaches (using bayescan, arlequin and outflank) found 28 outlier SNPs putatively under divergent selection and 9770 neutral SNPs in common. Redundancy analysis revealed that spatial distribution, ocean current‐mediated larval connectivity and SST explained 31.7% of the neutral genetic differentiation, with ocean currents driving the majority of this relationship (21.0%). After removing the influence of spatial distribution, no SST were significant for putatively neutral genetic variation whereas minimum annual SST still had a significant impact and explained 8.1% of the putatively adaptive genetic variation. Second, EA analyses (using Pearson correlation tests, bayescenv and lfmm) jointly identified seven SNPs as candidates for thermal adaptation. Covariation at these SNPs was assessed with a spatial multivariate analysis that highlighted a significant temperature association, after accounting for the influence of spatial distribution. Among the 505 candidate SNPs detected by at least one of the three approaches, we discovered three polymorphisms located in genes previously shown to play a role in thermal adaptation. Our results have implications for the management of the American lobster and provide a foundation on which to predict how this species will cope with climate change.
Adaptive genomics may help predicting how a species will respond to future environmental changes. Genomic signatures of local adaptation in marine organisms are often driven by environmental ...selective agents impacting the physiology of organisms. With one of the highest salinity level, the Mediterranean Sea provides an excellent model to investigate adaptive genomic divergence underlying salinity adaptation. In the present study, we combined six genome scan methods to detect potential genomic signal of selection in the striped red mullet (Mullus surmuletus) populations distributed across a wide salinity gradient. We then blasted these outlier sequences on published fish genomic resources in order to identify relevant potential candidate genes for salinity adaptation in this species.
Altogether, the six genome scan methods found 173 outliers out of 1153 SNPs. Using a blast approach, we discovered four candidate SNPs belonging to three genes potentially implicated in adaptation of M. surmuletus to salinity. The allele frequency at one of these SNPs significantly increases with salinity independently from the effect of longitude. The gene associated to this SNP, SOCS2, encodes for an inhibitor of cytokine and has previously been shown to be expressed under osmotic pressure in other marine organisms. Additionally, our results showed that genome scan methods not correcting for spatial structure can still be an efficient strategy to detect potential footprints of selection, when the spatial and environmental variation are confounded, and then, correcting for spatial structure in a second step represents a conservative method.
The present outcomes bring evidences of potential genomic footprint of selection, which suggest an adaptive response of M. surmuletus to salinity conditions in the Mediterranean Sea. Additional genomic data such as sequencing of a full-genome and transcriptome analyses of gene expression would provide new insights regarding the possibility that some striped red mullet populations are locally adapted to their saline environment.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Marine populations are typically characterized by weak genetic differentiation due to the potential for long‐distance dispersal favouring high levels of gene flow. However, strong directional ...advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (FST = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R2 = .49) than geographic distance (IBD) (R2 = .18), and directional processes played an important role in shaping fine‐scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine 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.
Interactions between environmental factors and complex life‐history characteristics of marine organisms produce the genetic diversity and structure observed within species. Our main goal was to test ...for genetic differentiation among eastern oyster populations from the coastal region of Canadian Maritimes against expected genetic homogeneity caused by historical events, taking into account spatial and environmental (temperature, salinity, turbidity) variation. This was achieved by genotyping 486 individuals originating from 13 locations using RADSeq. A total of 11,321 filtered SNPs were used in a combination of population genomics and environmental association analyses. We revealed significant neutral genetic differentiation (mean FST = 0.009) between sampling locations, and the occurrence of six major genetic clusters within the studied system. Redundancy analyses (RDAs) revealed that spatial and environmental variables explained 3.1% and 4.9% of the neutral genetic variation and 38.6% and 12.2% of the putatively adaptive genetic variation, respectively. These results indicate that these environmental factors play a role in the distribution of both neutral and putatively adaptive genetic diversity in the system. Moreover, polygenic selection was suggested by genotype–environment association analysis and significant correlations between additive polygenic scores and temperature and salinity. We discuss our results in the context of their conservation and management implications for the eastern oyster.
Marine reserves are being established worldwide to conserve and manage marine resources, and are also often assumed to conserve the evolutionary potential of marine populations. However, comparisons ...of genomic patterns inside and outside marine reserves in multiple species are scarce. Here, we aim to fill this gap by 1) comparing genomic variation in protected versus unprotected areas in a network of eight reserves in the Mediterranean Sea, and 2) disentangling the effect of seascape, space and marine reserves on this genomic variation in 1297 individuals from four species genotyped at more than 10 000 SNPs each. We report homogeneous patterns of local genetic diversity within and outside reserves in three of the four species (the white seabream Diplodus sargus, the striped red mullet Mullus surmuletus and the European spiny lobster Palinurus elephas). Surprisingly, the comber Serranus cabrilla shows lower genetic diversity in protected areas, highlighting an apparent conservation paradox that may be due to increased predation by protected fishes in the reserves. Multivariate analyses reveal that seascape factors (salinity, temperature and chlorophyll) and space exert a significant influence on genomic variation in some cases, while protection status has no detectable effect. Nevertheless, four marine reserves (Cabo de Palos, Cerbère‐Banyuls, Cap de Creus and Illes Columbretes) harbor singular genetic variation in one or two species. This pattern, observed in the four species, suggests that reserves contribute to preserve genetic variation locally. Our study illustrates that we can document key aspects of the genetic variation of marine species and their interaction with seascape factors, space and protection status in multiple species. It stresses the need for a multi‐species approach to inform marine conservation planning, opening up new perspectives at the community level.
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., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates.
Minor allele frequency correlation comparing GBS and Rapture.
The FishMed database provides traits, phylogeny, current and projected species distribution of Mediterranean fishes, and associated sea surface temperature (SST) from the regional oceanic model ...NEMOMED8. Data for the current geographical distributions of 635 Mediterranean fish species were compiled from a published expert knowledge atlas of fishes of the northern Atlantic and the Mediterranean (FNAM) edited between 1984 and 1986 and from an updated exotic fish species list. Two future sets of projected species distributions were obtained for the middle and end of the 21st century by using an ensemble forecasting approach for 288 coastal Mediterranean fish species based on SST according to the IPPC/SRES A2 scenario implemented with the Mediterranean climatic model NEMOMED8. The functional part of the database encompasses 12 biological and ecological traits (maximal and common lengths, vertical distribution, habitat, migration type, mode of reproduction, sex shift, semelparity, diet type (larvae and adults), social behavior, species origin, and depth) for the 635 fish species. To build the phylogeny we inferred the timing and geographic origins of Mediterranean teleost species diversity using nucleotide sequences collected from GenBank including 62% of Mediterranean teleost species plus nine outgroups. Maximum likelihood Bayesian phylogenetic and dating analyses were calibrated using 20 fossil species. An additional 124 fish species were grafted onto the chronogram according to their taxonomic affinity to obtain a phylogenetic tree including 498 species. Finally we also present the associated SST data for the observed period (1961-1980) and for the middle (2040-2059) and the end of the 21st century (2080-2099) obtained from NEMOMED8 according to the IPCC A2 scenario. The FishMed database might be of interest in the context of global anthropogenic changes as coastal Mediterranean ecosystems are currently recognized as one of the most impacted ecosystems on earth.