10 which environmental factor could lead to a decrease in genetic variation in a population of tuna? Ideas

Contents

Below is information and knowledge on the topic which environmental factor could lead to a decrease in genetic variation in a population of tuna? gather and compiled by the show.vn team. Along with other related topics like: .


NA SNP-panel data suggests stability in bluefin tuna genetic diversity despite centuries of fluctuating catches in the eastern Atlantic and Mediterranean

  • Journal List
  • Scientific Reports
  • PMC8528830

Sci Rep. 2021; 11: 20744.

Adam J. Andrews,corresponding author#1,2 Gregory N. Puncher,corresponding author#1,3 Darío Bernal-Casasola,4 Antonio Di Natale,5 Francesco Massari,1 Vedat Onar,6 Nezir Yaşar Toker,6 Alex Hanke,7 Scott A. Pavey,3 Castrense Savojardo,8 Pier Luigi Martelli,8 Rita Casadio,8 Elisabetta Cilli,2 Arturo Morales-Muñiz,9 Barbara Mantovani,10 Fausto Tinti,1 and Alessia Cariani1

Abstract

Atlantic bluefin tuna (Thunnus thynnus; BFT) abundance was depleted in the late 20th and early 21st century due to overfishing. Historical catch records further indicate that the abundance of BFT in the Mediterranean has been fluctuating since at least the 16th century. Here we build upon previous work on ancient DNA of BFT in the Mediterranean by comparing contemporary (2009–2012) specimens with archival (1911–1926) and archaeological (2nd century BCE–15th century CE) specimens that represent population states prior to these two major periods of exploitation, respectively. We successfully genotyped and analysed 259 contemporary and 123 historical (91 archival and 32 archaeological) specimens at 92 SNP loci that were selected for their ability to differentiate contemporary populations or their association with core biological functions. We found no evidence of genetic bottlenecks, inbreeding or population restructuring between temporal sample groups that might explain what has driven catch fluctuations since the 16th century. We also detected a putative adaptive response, involving the cytoskeletal protein synemin which may be related to muscle stress. However, these results require further investigation with more extensive genome-wide data to rule out demographic changes due to overfishing, and other natural and anthropogenic factors, in addition to elucidating the adaptive drivers related to these.

Subject terms: Ecological genetics, Evolutionary ecology, Marine biology

Introduction

Overfishing has reduced numerous fish populations to remnants of their historical levels1,2, yet we have a poor understanding of what impact this has had on their evolutionary potential and resilience3. This information is crucial to predict future demographic changes and thus promote sustainable fisheries management4,5. Studies of historical marine ecology offer an opportunity to learn and heed these past lessons57. In particular, genetic/genomic studies can infer past history from contemporary samples8, or directly test archaeological and archival samples9 for losses in genetic diversity, population restructuring, or adaptive responses to natural factors e.g., climate, or anthropogenic ones e.g., fisheries-induced evolution (FIE)10. A decade ago, Riccioni et al.11 were the first to investigate temporal demographic changes in the key species Atlantic bluefin tuna (Thunnus thynnus, hereafter BFT) using archival early-20th century samples and microsatellite markers. Here, we build on this work by genotyping archival and archaeological samples to extend investigations into the pre-industrial era, when fishing may have also had the potential to impact BFT.

BFT is a highly migratory pelagic top predator, characterized by its large size (up to 3.3 m in length and 725 kg in weight), slow maturation (between 4 and 8 years)12,13, and inshore migration behaviour, that has made it vulnerable to overfishing. Recent genomic studies14,15 support the delineation of two BFT populations. These are a western Atlantic component that spawns predominantly in the Gulf of Mexico16, and an eastern Atlantic and Mediterranean component that spawns predominantly in the Mediterranean Sea17. Individuals of both populations migrate into the Atlantic Ocean to feed, including as juveniles18, and exhibit high-levels of mixing14,15. The role of additional contemporary and historical spawning areas i.e. the Slope Sea (East of Cape Hatteras, USA)16, the Bay of Biscay19, and the Black Sea20,21, are yet to be clearly defined, especially regarding the Slope Sea where connectivity between populations was observed15.

During the last few years, the eastern Atlantic and Mediterranean population of BFT has recovered to 1970’s levels following heavy overfishing that depleted spawning stock biomass, restructured the population toward younger individuals, and contracted the species range, in the late 20th and early 21st century2225. However, reconstructions of 16th-20th century BFT trap catch records suggest abundance across the Mediterranean has been fluctuating for centuries23,26,27. Pelagic species are particularly susceptible to fluctuations in abundance since dynamic food and environmental conditions drive large variability in recruitment, but fishing magnifies poor recruitment and therefore population declines when large catches occur28. Multiple factors need to be taken into account to be able to interpret trap catch fluctuations as abundance29, though, it appears that catch numbers in the 16th and 18th century may be comparable to those during the industrial fishing of the last 50 years26,27. Hence, fishing appears to have been intense in this period.

The current study investigates genetic variability in eastern Atlantic and Mediterranean BFT prior to both their 21st century population collapse, and record trap catches in the 16th and 18th century, using archived early-20th century specimens, and archaeological remains, respectively. Despite overfished species having an overall lower genetic diversity when contemporary data are compared3, Riccioni et al.11 were unable to detect losses in BFT genetic diversity when comparing contemporary and early-20th century samples. Likewise, no genetic erosion was observed following overfishing in the closely-related albacore (Thunnus alalunga)30, the Pacific herring (Clupea pallasii)31, or the European anchovy (Engraulis encrasicolus)32. Even marine species (e.g., sawfish, Pristis spp.) depleted to between 1 and 5% of their historical biomass appear to have retained genetic diversity33. However, several studies have noted genetic diversity declines or population losses following overfishing in Atlantic cod (Gadus morhua)3436, Atlantic salmon (Salmo salar)37, and Chinook salmon (Oncorhynchus tshawytscha)38. In addition, adaptive responses to size- (in Walleye, Sander vitreus)39, and sex- (in Atlantic salmon)40 selective harvesting, and environmental drivers (in Atlantic cod)41 have also been reported in studies using archival or archaeological fish samples. It remains unclear to what extent the inability of some studies to detect these differences results from the selection of genetic markers with low resolution, or the resilience offered by complex life history traits during times of population decline. A recent whole genome sequencing (WGS) study42, that did not detect genetic erosion or adaptive responses in two Atlantic cod populations following 20th century overfishing, may indicate, however, that the latter is the case for some populations.

Here we test the hypotheses that the genetic diversity of BFT has declined and that their populations restructured following periods of intense fishing in the eastern Atlantic and Mediterranean. Further, we sought to identify adaptive responses that may be related to ecological or environmental conditions. To this end, our objectives were to genotype archaeological and archival specimens on a single nucleotide polymorphic (SNP) panel to (1) characterise their genetic diversity and population structure, (2) to compare those patterns to analogous ones from contemporary groups, and (3) to explore markers under putative selection and identify their associated function, if possible.

Methods

Samples

We collected samples of contemporary, archival and archaeological BFT specimens for analysis as follows: Contemporary reference specimens (GOM: Gulf of Mexico, CMAS: Central Mediterranean Adriatic Sea, CMSI: Central Mediterranean Sicily, EABB: East Atlantic Bay of Biscay, EAGI: Eastern Atlantic Gibraltar, EMLS: Eastern Mediterranean Levantine Sea, WMBA: Western Mediterranean Balearic Islands, WTYR: Western Mediterranean Tyrrhenian Sea, n = 277, Table S1) at each life stage were collected across the species range between 2009 and 2012 (Fig. 1, Table S1) where tissue samples from each specimen were preserved in 96% ethanol or RNAlater (Thermo Fisher Scientific, USA) and stored at − 20 °C until further processing. Archived vertebrae (HBOS: Historical Bosporus, HADR: Historical Adriatic Sea, HION: Historical Ionian Sea, HTYR: Historical Tyrrhenian Sea: n = 147, Table S1) from the Massimo Sella Archive (see11) were collected between 1911 and 1941 (Fig. 1). Archaeological vertebrae (n = 136, Table S1) were retrieved from several excavations (Fig. 1) including 4th–15th century CE Yenikapi (HIST: Historical Istanbul, Turkey)43, 2nd century BCE–5th century CE Baelo Claudia (HBC: Historical Baelo Claudia, Spain)44, 2nd century BCE Tavira (HTAV: Historical Tavira, Portugal), and 4th–2nd century BCE Palacio de Justicia, (HPJ: Historical Palacio de Justicia, Spain)45. See Supplementary Materials 1 for more details on historical samples and their dating.

An external file that holds a picture, illustration, etc.
Object name is 41598_2021_99708_Fig1_HTML.jpg

Map of the collection location for samples used in analyses. Historical (archival and archaeological) sample groups (in boldface, denoted with H) use approximate locations and the locations of archaeological sites where fish remains were recovered. Map created using ESRI ArcMap (v.10.6, https://arcgis.com). Only sample groups that were successfully genotyped and analysed are displayed. Numbers (n) represent those included in the final analysis for each sample group.

Contemporary DNA extractions

DNA was isolated from fin (adults) or muscle (juveniles, young-of-the-year) of contemporary samples (Table S1) as part of another study14 using the Wizard®SV96 Genomic DNA Purification Kit (Promega, USA), following the manufacturer’s instructions. Quantification was performed using a NanoDrop 2000 (Thermo Fisher Scientific, USA). Negative controls indicated that no cross-contamination took place between samples.

Ancient DNA extractions

Archival and archaeological samples underwent ancient DNA (aDNA) extraction in sterile, PCR-free conditions at the Ancient DNA Laboratory of the Department of Cultural Heritage (University of Bologna, Ravenna Campus, Italy), as part of another study which investigated their species identification via barcoding46. All bone specimens were sprayed with 1–2% sodium hypochlorite (bleach), left to soak for ten minutes, rinsed with distilled water and air-dried (as per47). Specimens were then mechanically cleaned using sandpaper, and the bleaching process was repeated. After, each specimen was exposed to UV light (254 nm) for 15 min before drilling to obtain ~ 200 mg bone powder. Bones that were too small for drilling were bisected, and their inner matrices were crushed.

Isolation of aDNA was performed using a modified version of Dabney et al.48,49. Briefly, ~ 200 mg of bone powder from each sample was divided in two and placed into separate tubes. After an overnight incubation in EDTA (0.5 M, pH 8.0) and proteinase K, lysates (1000 µl) of each sample were pooled and combined with 3000 µl binding buffer composed of guanidine thiocyanate (5 M), Tween 20 (0.05%), isopropyl alcohol (40% v/v), sodium acetate (90 mM, pH 5.2), and distilled water. This mixture was then centrifuged through a MinElute spin column (Qiagen, Germany), and washed twice with 720 µl PE buffer, before a final elution in 60 µl of distilled water.

The total DNA obtained from each extraction was quantified using a Qubit® dsDNA HS (High Sensitivity) Assay Kit (Thermo Fisher Scientific, USA). Negative controls employed for each batch of samples extracted indicated an undetectable level of contamination (< 500 pg/ml).

DNA genotyping

A total of 273 contemporary samples, and 280 historical (145 archival and 135 archaeological) samples contained sufficient quantities of DNA (100 ng total) for genotyping (Table S1). Samples were genotyped using a 96 SNP-panel we developed from SNP’s identified by two studies14,30 that were polymorphic between contemporary sample groups (see14) or matched with gene functions. To identify protein association we blasted the flanking regions of these loci against sequences for Atlantic cod50, sea bass (Dicentrarchus labrax)51, BFT52,53 and an umbrella set of teleost sequences, on NCBI GenBank (https://blast.ncbi.nlm.nih.gov/Blast.cgi, Table S2) using the blastn option. Queries were considered matches if alignment coverage was > 80% and identity scores were > 80% (Table S2).

SNP genotyping was conducted first using Fluidigm 96.96 Dynamic Array™ Integrated Fluidic Circuits (Probes: SNPtype-FAM:SNPtype-HEX, Passive reference: ROX) on the BioMarkHD034 platform (SGIKER, Spain). Historical samples were re-genotyped at a second facility using the Fluidigm EP1 platform (ABL, Bedford Institute of Oceanography, Canada) to assess genotype error rates. Genotyping employed two negative controls for each run, which confirmed no cross-contamination, and three positive controls (CMAS01, CMAS02, CMAS03), reporting identical genotypes. Similarly, 21 (7.5%) historical samples were extracted and genotyped twice and reported acceptable replicates at 97.8 ± 3.6% accuracy.

Quality control filtering

Prior to analyses, two loci (SNP85, SNP86, Table S2) with low call rates (98–100% missing data) were discarded. Individuals (148 out of 553, 26.7%) and two further loci (SNP45, SNP79) that contained > 10% missing data were then removed. Inconsistencies between the two facilities at the remaining 92 loci were then assessed. The remaining 146 historical individuals were subject to further filtering whereby they were removed if their genotypes were inconsistent between the two facilities at > 5% of loci. This removed a further 21 (14.4%) historical individuals achieved an overall genotyping success of 98.8% at 92 loci. Sample groups that contained a single individual as a result of filtering (HBOS, Table S1) were also removed. Historical duplicate samples resulting from the potential sampling of two or more bone specimens of the same individual were identified and removed by applying the function clonecorrect in the Poppr package54 as implemented in R v.4.0.355. A single clone was evident in the HIST archaeological sample group (Table 1).

Table 1

Genetic diversity, Hardy–Weinberg deviation, and FIS in contemporary and historical (archival and archaeological) sample groups.

Contemporary Historical
GOM CMAS CMSI EABB EAGI EMLS WMBA WMTY HADR HION HTYR HIST HTAV
N 22 39 36 40 29 29 24 40 16 40 35 30 2
aRa 177 178 178 178 178 176 176 177 175 177 177 177 139
Hea 0.346 0.369 0.364 0.359 0.369 0.367 0.369 0.365 0.380 0.358 0.362 0.367 0.353
Hoa 0.332 0.358 0.358 0.344 0.326 0.372 0.377 0.340 0.357 0.344 0.355 0.353 0.305
PHW 0.936 0.932 0.973 0.980 1.000 0.197 0.155 1.000 0.977 0.992 0.902 0.971 0.956
FISa 0.032 0.024 0.033 0.034 0.137 − 0.018 − 0.022 0.066 0.052 0.041 0.023 0.038

Loci evaluation

Deviation from Hardy–Weinberg equilibrium (HWE) was assessed at each locus using the R package Pegas56. Linkage disequilibrium (LD) between loci was tested using the R package Genepop57. Outlier loci were identified using Bayescan58 and OutFLANK v0.259 to obtain a neutral dataset and identify potential adaptive responses. Analysis was run excluding the western Atlantic sample group (GOM) between the following: all sample groups, pooled contemporary and historical sample groups, contemporary sample groups, and historical sample groups. Loci detected as outliers were removed from the dataset prior to demographic analyses and investigated as follows: gene association was inferred from the above blastn searches, and non-synonymous mutations were explored with the Expasy Translation tool as implemented online (https://web.expasy.org/translate). Default settings were used in the analyses. Significance was judged using the False Discovery Rate (FDR)60 approach at the 5% level, as calculated using 999 permutations.

Read More:  10 how does the law of conservation of mass apply to chemical changes Ideas

Genetic diversity

Allelic richness (aR), heterozygosity expected/observed (He, Ho), and the inbreeding coefficient FIS, were calculated for each sample group with the R package Hierfstat61. Significance of heterozygote excess was calculated with Genepop in R using the global excess method and default settings. Differences in aR, He, Ho and FIS between pooled contemporary and historical sample groups were assessed using unpaired t-tests in R. Significance was judged at the 5% level. Effective population size (Ne) estimates were calculated only for samples consistently scored across all 89 neutral loci, as summarised in Table 2. Estimations were calculated using the linkage disequilibrium approach62 as implemented in NeEstimator v2.163 and an allele frequency threshold of 0.01. A random down-sampling to generate and analyse equal-size groups is summarised in Table S3. Because Ne estimates are often unreliable at low sample sizes64, we calculated per locus round-robin estimates of minor allele frequencies in R (as per65) and plotted trajectories between temporal sample groups. We performed a hierarchical analysis of molecular variance (AMOVA) in Poppr, with 10,000 permutations to assess significance. AMOVAs were performed excluding the GOM sample group on the following levels: between periods; between sample groups; between samples (i.e., individuals); and within samples.

Table 2

Effective population size (Ne) and 95% confidence Intervals of contemporary and historical sample groups for samples (n) consistently scored across all 89 neutral loci analysed herein, under two approaches, where separate estimates were made for each sample group and for contemporary and historical pools.

Contemporary Historical
GOM CMAS CMSI EABB EAGI EMLS WMBA WMTY HADR HION HTYR HIST
Separate
n 14 32 27 30 14 16 36 24 24 18
Ne 140 164 391 1454 2515 140 150 825 60 14
CI 56 − ∞ 78–799 111 − ∞ 158 − ∞ 68 − ∞ 46 − ∞ 53 − ∞ 118 − ∞ 40–114 11–18
Pooled
n 32 27 30 4 14 16 36 4 24 24 18
Ne 939 298
CI 497–5465 178–787

Population structure

A discriminant analysis of principal components (DAPC) was performed with the R package Adegenet66 to explore how the historical groups relate to the contemporary reference groups. DAPC is a geometric clustering method free of HWE and LD assumptions, that attempts to maximise the inter-variation between clusters while minimising the intra-variation observed within clusters. DAPC clusters were set a priori to the number of sample groups. We retained 4 discriminant functions and the number of principal components (PC’s) according to the function optim.a.score, based on an initial selection of all PC’s before refinement. Population structuring was also evaluated using STRUCTURE v.2.3.467, which implements a Bayesian clustering method to identify the most likely number of populations (K). We followed the Evano et al.68 method, and thus, we carried out 10 runs per each value of K ranging from 1 to 10. Runs used the locprior and admixture models and assumed correlated allele frequencies. Each run used 500,000 burn-in and Markov Chain Monte Carlo replicates. We estimated the ad hoc statistic ΔK in order to infer the most likely number of populations using STRUCTURE HARVESTER69. CLUMPAK70 was used to merge the 10 runs from the most probable K, and reported similarity scores > 95. We used a hierarchical approach to improve resolution due to the identification of 5 outliers (EAGI 6 & 17, WMTY 52, 57 & 66) in two modern sample groups that constituted two distinct populations at K = 3 in the first run. Hence, these individuals were removed from the dataset and STRUCTURE was run a second time. Pairwise distances between sample groups were calculated with Nei’s estimator of FST71 in the hierfstat R package, using 999 permutations to calculate the respective p-values, which were judged for significance under the FDR approach at the 5% level.

Results

Loci evaluation

Overall, 259 contemporary, and 123 historical (91 archival and 32 archaeological) samples were analysed at 92 loci (Table S1). No loci deviated from HWE or were in LD in more than a single population. BayeScan and OutFLANK both detected three loci (SNP41, SNP43 & SNP93, Table S1, Figure S1) as outliers. Loci SNP41 and SNP43 were outliers between contemporary sample groups and locus SNP93 was an outlier between historical sample groups. Locus SNP41 was identified as a putative adaptive response after being detected as an outlier between pooled contemporary and historical groups. Locus SNP41 was found to be in potential association with the gene SYNM that encodes Synemin, which is an intermediate filament protein. This putative adaptive locus was found to be under selection in all contemporary sample groups except CMSI, comprising a nucleotide mutation (T to A) that was non-synonymous, resulting in the production of glutamine instead of histidine. In contrast, SNP41 was not under selection in a single historical sample group.

Genetic diversity

We found no significant differences in gene diversity aR (p = 0.181, t(11) = 1.426), He (p = 0.923, t(11) = 0.099) and Ho (p = 0.575, t(11) = 0.578) between pooled contemporary and historical groups (Table 1). Heterozygote deficiency was not significant in any sample group (Table 1). Inbreeding (FIS) was rare within all sample groups (Table 1) and was not significantly different between pooled contemporary and historical samples (p = 0.939, t(9) = 0.0791). The dataset lacked power to define reliable estimates of Ne using both methods for each sample group i.e., our CIs contained infinity until they were pooled (Table 2). Randomly excluding samples to create equal size sample groups had minimal influence on estimations (Table S3). Ne estimates were higher for both contemporary sample groups, analysed separately, and the contemporary eastern Atlantic and Mediterranean when pooled (Table 2). Allele trajectories (Figure S1) showed stochastic fluctuations in minor allele frequencies between all sample groups, and no consistent drop-out or over-estimation in all contemporary or historical sample groups, respectively. Within the eastern Atlantic and Mediterranean samples, AMOVAs indicated significant differences in variance within and between samples, and between sample groups, but not between periods (Table 3, Figure S2).

Table 3

Variance of the eastern Atlantic and Mediterranean BFT samples as computed by AMOVAs using a hierarchical approach as indicated by the four levels.

Σ % variance p value
Within samples 31.203 96.010 < 0.001
Between samples 1.203 3.703 0.001
Between sample groups 0.087 0.268 0.003
Between contemporary and historical groups 0.005 0.017 0.306

Population structure

DAPC clustered eastern Atlantic and Mediterranean sample groups together while the western Atlantic (GOM) sample group was substantially separated (Fig. 2). Considerable overlap was observed between contemporary and historical clusters of the eastern Atlantic and Mediterranean. ΔK suggested that the most likely number of populations identified with STRUCTURE was K = 3. All individuals shared mixed membership (q). Separate structuring of the GOM sample group was evident and the historical sample group HIST contained three individuals with this signature (Fig. 3). Overall, no evidence of population structure was evident between contemporary or historical sample groups of the eastern Atlantic and Mediterranean (Figs. 2, 3). Pairwise FST values were significant between the GOM sample group and all others (Table 4). In addition, the sample groups EMLS and WMBA, and EAGI and HIST were significantly different. No other significant differences were observed between contemporary and historical sample groups.

An external file that holds a picture, illustration, etc.
Object name is 41598_2021_99708_Fig2_HTML.jpg

Discriminant analysis of principal components scatterplot showing how historical (archival and archaeological, denoted with H) sample groups relate to contemporary reference populations of the Gulf of Mexico (GOM) and the eastern Atlantic and Mediterranean. DAPC cluster ellipses were set to contain 95% of genotypes. Discriminant analysis (DA) eigenvalues and principal component analysis (PCA) eigenvalues were selected as displayed to avoid overfitting, utilising the optim.a.score approach within the R package adegenet.

An external file that holds a picture, illustration, etc.
Object name is 41598_2021_99708_Fig3_HTML.jpg

STRUCTURE barplot showing membership probabilities (q) for each sample group analysed herein with K = 3 (each represented by a different shade). K = 3 was the most likely number of populations identified by the ΔK method. Historical (archival and archaeological) sample groups are denoted with H).

Table 4

Pairwise FST (below the diagonal) and non-corrected P values (above the diagonal) between contemporary and historical sample groups.

Contemporary Historical
GOM CMAS CMSI EABB EAGI EMLS WMBA WMTY HADR HION HTYR HIST HTAV
GOM 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001
CMAS 0.0154 0.643 0.387 0.065 0.153 0.023 0.127 0.446 0.810 0.351 0.450 0.425
CMSI 0.0128 − 0.0009 0.566 0.053 0.448 0.045 0.364 0.0295 0.602 0.029 0.538 0.221
EABB 0.0118 0.0006 − 0.0004 0.108 0.111 0.174 0.416 0.749 0.635 0.593 0.625 0.702
EAGI 0.0240 0.0041 0.0046 0.0033 0.002 0.124 0.045 0.155 0.249 0.004 0.001 0.615
EMLS 0.0255 0.0027 0.0003 0.0034 0.0096 0.001 0.002 0.194 0.015 0.009 0.398 0.568
WMBA 0.0171 0.0064 0.0051 0.0029 0.0037 0.0109 0.116 0.689 0.273 0.007 0.085 0.103
WMTY 0.0167 0.0026 0.0007 0.0003 0.0047 0.0087 0.0035 0.099 0.022 0.197 0.597 0.413
HADR 0.0228 0.0004 0.0017 − 0.0025 0.0040 0.0037 − 0.0024 0.0048 0.853 0.248 0.876 0.211
HION 0.0129 − 0.0018 − 0.0006 − 0.0008 0.0017 0.0059 0.0016 0.0049 − 0.0040 0.217 0.562 0.238
HTYR 0.0117 0.0016 0.0053 − 0.0004 0.0083 0.0070 0.0087 0.0019 0.0026 0.0017 0.155 0.351
HIST 0.0155 0.0003 − 0.0003 − 0.0007 0.0097 0.0008 0.0042 − 0.0007 − 0.0047 − 0.0004 0.0027 0.594
HTAV 0.0350 0.0037 0.0183 − 0.0127 − 0.0079 − 0.0046 0.0322 0.0035 0.0201 0.0168 0.0085 − 0.0066

p values that were significant after FDR correction are presented in boldface.

Discussion

We found no evidence of genetic diversity loss or population restructuring in contemporary BFT sample groups of the eastern Atlantic and Mediterranean compared with those from the early 20th century CE prior to spawning biomass depletion and species range contraction22,25, and the 4th–15th century CE prior to a significant period of intense trap fishing marked by fluctuating catches26,27. If overfishing had resulted in a genetic bottleneck, we would expect to see significant decreases in minor allele frequencies, allelic richness, and observed heterozygosity72 for contemporary samples compared with historical samples. Therefore, we would also expect to observe an increase in inbreeding and a decrease in effective population size8, which we did not. The impact of overfishing on genetic diversity and allele frequencies has been observed in a variety of studies that directly test archaeological and archival samples3438. At its most extreme, overfishing has been observed to restructure marine fish populations36, yet we found no evidence of genetic restructuring in BFT. Likewise, a recent study found that Atlantic cod had not been impacted by 20th century overfishing at the genomic level42. Our findings are similar to those of Riccioni et al.11 using microsatellite markers, though we did not observe significant sub-structuring within Mediterranean BFT as they did, and this is yet to be resolved to clarify alternative population structure hypotheses18,23,73. No recent genetic study, however, has detected population structure within the eastern Atlantic and Mediterranean BFT14,15,74. Perspectives from threatened populations of other taxa, inform us that a wide range of genomic responses are expected, along a continuous scale from resistance to collapse34,7579, and recovery33,8082. Despite differences between taxa, these data would suggest that there is likely no “one-size fits all” response to the depletion of marine fish populations, according to species life history traits and the extent and rate of overfishing.

The most common explanation for the maintenance of genetic diversity in threatened populations is that gene flow acts as a buffer72. This is plausible for BFT, though its western Atlantic population is smaller than the eastern Atlantic and Mediterranean population (ca. 1/10 the size) and was heavily fished itself since at least the early 19th century83,84. Connectivity with alternative spawning sites (e.g., the Slope Sea16, the Bay of Biscay19, Azores, Canary Islands, Ibero-Moroccan, Gulf of Guinea13,17) remains poorly understood, and the unresolved frequency and duration of spawning at these locations means we cannot assess its effect on gene flow. Likewise, introgression occurs at a low rate between Thunnus species85 but could also act as a buffer. On the other hand, eastern BFT may be resilient towards genetic erosion due to their relatively large population size (enhanced by connectivity between spawning sites within the Mediterranean), and a long life cycle which promotes heavily overlapping generations13,17. In any case, our findings leave us with two possible explanations; either (1) overfishing was not severe enough to cause a genetic bottleneck in BFT, or (2) our observation of significant demographic changes was hindered by the methods we employed.

To address this first point, it is evident that BFT were overfished, at least in the 20th and early 21st century, if not as we suspect between the 16th–19th century. Studies by the management body ICCAT (the International Commission for the Conservation of Atlantic Tunas)25, and independent estimates e.g.,22 suggested that BFT abundance and range declined by 70% and 46–53%, respectively, between 1960 and 2010. However, there is debate on the extent of the population decline, where on one hand, impending population collapse was predicted in 200986, yet on the other, poorly understood population dynamics and incorrect assignments of catches has caused uncertainty in population estimates73,87. Hence, it is difficult to deduce whether we should expect to find evidence of a genetic bottleneck because the recent recovery of the population within just two generations from its lowest point in 200788 could suggest that either the population decline was not that severe, or that overfishing did trigger a severe population decline but BFT is remarkably resilient due to its complex life history traits.

Nonetheless, fishing effort is not the only factor that influences catches and abundance (as shown for the historical trap fishery data29), which one might expect to be reflected in genetic diversity and structure. Climate is likely the largest regulator of recruitment and thus fish abundance89,90 and as a pelagic species, BFT are certainly no exception91,92. Therefore, one might expect to find evidence of fluctuating abundance—and potentially genetic diversity—that is merely exacerbated by fishing28. BFT’s Atlantic distribution varies with Atlantic multidecadal oscillation phases93, and thus gene flow and inbreeding is expected to vary accordingly because connectivity of populations is enhanced in warm years as ranges overlap, as attested by isotope data94. This is notwithstanding time-related effects driven by evolutionary processes i.e., mutation and genetic drift that we might expect to alter allele frequencies over time. Therefore, our observation of homogeneity between contemporary and historical BFT samples is somewhat striking. One might pose the question: at what rate should we expect to observe demographic changes at the genomic level? We analysed moderate sample sizes from 1911 to 1926 (~ 20 generations ago) and the 4th–15th century (~ 100 + generations ago), yet we did not detect time-related effects. Thus, to address this, even at conservative mutation rates lower than those shown for marine fish95, we would expect to observe changes in allele frequencies as a result of genetic drift alone.

Read More:  10 how fast does mars rotate on its axis Ideas

Alternatively, our observations may be explained by our methodological approach. By pre-selecting loci that were polymorphic in contemporary sample groups, our data are subject to an unknown degree of ascertainment bias. Theoretically, ascertainment bias could influence any analysis or inference based on SNP allele frequencies when SNPs are discovered in a limited sample but applied in another context (e.g., our historical sample groups)96. The expectation that this should inflate diversity in the ascertainment sample is a widely accepted hindrance of SNP-panel studies9799. Studies usually correct for this by LD pruning100 or modifying raw genotypes following maximum-likelihood simulations98, however this was not possible herein due to the few loci that were available. Indeed, the effect of ascertainment bias is likely to be exacerbated herein because we analysed few loci. This reduces the likelihood of detecting rare alleles and thereby erodes power97 which is particularly crucial when differentiating marine samples due to high gene flow and low diversity in marine populations101. Theoretically, this might have inflated our estimates of genetic diversity among contemporary samples, and hence genetic diversity was comparatively low in historical samples. This theory is further supported by our AMOVA results and might explain why variance was lower than expected between temporal samples, and why structure was only observed between contemporary sample groups for which SNP discovery was made.

Moreover, our Ne estimate CIs often contained infinity, suggesting that we have little power to make any inferences on Ne. In many cases Ne was strikingly lower in (supposedly unimpacted) historical samples than the empirical rule-of-thumb threshold of Ne (500) proposed to maintain long-term genetic diversity in marine populations72. In any case, Ne is often unreliable when using sample sizes such as ours64 and we caution the interpretation of our results for this reason. Additionally, our sampling strategy may have been limiting. For example, if genetic diversity had decreased following population declines (e.g., between 16th–18th century, and/or during the 20th and early 21st century) but was restored prior to our analogous archival samples of the early 20th century, or 2009–2012 contemporary samples, respectively. Species differ in their rates of genetic recovery according to their life history traits72, and as this rate is unknown in BFT, we cannot rule out this possibility.

Clearly aDNA approaches offer utility to fisheries management because long-term trends are understudied and we lack fisheries-independent indices7,90. However, genome-wide approaches are more likely to provide a better resolution to assess demographic impacts and adaptive responses. Assuming the availability of a reference genome, WGS approaches are increasingly cost-effective9, particularly for shallow sequencing102. This approach may also facilitate the recovery of data from arid Mediterranean specimens which were challenging to genotype herein due to their poor preservation46. Importantly, WGS would reduce ascertainment bias compared with SNP-genotyping. This is crucial where allele frequency distributions are used to infer demographic history, but also to scan for past targets of selection98. We were limited herein to detecting a single adaptive response in BFT: a potential change to the function of the protein synemin, which is a cytoskeletal protein that we speculate might be related to growth changes induced by size selective harvesting (FIE), although this remains to be tested. WGS studies able to detect additional loci under putative selection are ultimately required for the association of this response (and others) with natural or anthropogenic factors, in addition to discounting hitchhiking effects103.

Conclusion

We identify that aDNA preserved within archival and archaeological fish remains has the potential to inform fisheries management by providing novel fisheries-independent baselines with which to observe unstudied long-term demographic and adaptive changes. We found no evidence that genetic diversity decreased or that populations restructured following several centuries of intense fishing, in line with a previous study11. This may hint at BFT’s resilience which has been recently shown by rebounds in abundance25 and a return to previous habitats88. However, we acknowledge limitations in our dataset i.e., few markers and the potential for ascertainment bias, and suggest that future studies might benefit from obtaining WGS data to observe rare alleles and reduce bias. Genome-wide data will be especially necessary to investigate adaptive responses, such as the putative selection on the cytoskeletal protein synemin found herein, and associate these with natural or anthropogenic factors to elucidate the drivers of change.

Supplementary Information

Acknowledgements

We thank two anonymous reviewers whose comments improved the quality of the manuscript. This work was carried out under the provision of the FARB Project (University of Bologna), the ICCAT Atlantic Wide Research Programme for Bluefin Tuna (GBYP) Data Recovery Plan, the International Governance Strategy of Fisheries and Oceans Canada (F5211-180183 and F5211-160166), and the MARES Joint Doctorate programme selected under Erasmus Mundus and coordinated by Ghent University (Grant Number: FPA 2011‐0016). This work is also a contribution to the MSCA SeaChanges ITN and was part-funded by EU Horizon 2020 (Grant Number: 813383).

Author contributions

A.C & F.T. conceived the study. G.N.P, A.C., F.T., B.M., A.M.M., N.Y.T, V.O. and D.B.C. collected tissues for analysis. G.N.P., E.C., A.H., S.A.P. and F.M. conducted the lab work. A.D.N., C.S., P.L.M. and R.C. provided data. A.J.A. performed the statistical analysis. A.J.A., G.N.P., A.C. and F.T. wrote the manuscript. All authors reviewed the manuscript.

Data availability

Flanking region sequences for each locus, and genotypes for all individuals, are attached as supplementary files.

Competing interests

The authors declare no competing interests.

Footnotes

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Adam J. Andrews and Gregory N. Puncher.

Contributor Information

Adam J. Andrews, Email: gro.emoealap@mada.

Gregory N. Puncher, Email: moc.liamg@rehcnupgerg.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-021-99708-9.

References

1. Pauly D, et al. Towards sustainability in world fisheries. Nature. 2002;418:689–695. doi: 10.1038/nature01017. [PubMed] [CrossRef] [Google Scholar]

2. Butchart SHM, et al. Global biodiversity: Indicators of recent declines. Science. 2010;328:1164–1168. doi: 10.1126/science.1187512. [PubMed] [CrossRef] [Google Scholar]

3. Pinsky ML, Palumbi SR. Meta-analysis reveals lower genetic diversity in overfished populations. Mol. Ecol. 2014;23:29–39. doi: 10.1111/mec.12509. [PubMed] [CrossRef] [Google Scholar]

4. Neubauer P, Jensen OP, Hutchings JA, Baum JK. Resilience and recovery of overexploited marine populations. Science. 2013;340:347–349. doi: 10.1126/science.1230441. [PubMed] [CrossRef] [Google Scholar]

5. Lotze, H. K., Hoffmann, R. & Erlandson, J. Lessons from historical ecology and management. In The Sea, Volume 19: Ecosystem-Based Management (Harvard University Press, 2014).

6. Erlandson JM, Rick TC. Archaeology meets marine ecology: The antiquity of maritime cultures and human impacts on marine fisheries and ecosystems. Ann. Rev. Mar. Sci. 2010;2:231–251. doi: 10.1146/annurev.marine.010908.163749. [PubMed] [CrossRef] [Google Scholar]

7. Schwerdtner Máñez, K. et al. The future of the oceans past: Towards a global marine historical research initiative. PLoS ONE9, e101466 (2014). [PMC free article] [PubMed]

8. Palsbøll, P. J., Zachariah Peery, M., Olsen, M. T., Beissinger, S. R. & Bérubé, M. Inferring recent historic abundance from current genetic diversity. Mol. Ecol.22, 22–40 (2013). [PubMed]

9. Oosting T, et al. Unlocking the potential of ancient fish DNA in the genomic era. Evol. Appl. 2019;12:1513–1522. doi: 10.1111/eva.12811. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

10. Heino M, Pauli BD, Dieckmann U. Fisheries-induced evolution. Annu. Rev. Ecol. Evol. Syst. 2015;46:461–480. doi: 10.1146/annurev-ecolsys-112414-054339. [CrossRef] [Google Scholar]

11. Riccioni G, et al. Spatio-temporal population structuring and genetic diversity retention in depleted Atlantic Bluefin tuna of the Mediterranean Sea. Proc. Natl. Acad. Sci. 2010;107:2102–2107. doi: 10.1073/pnas.0908281107. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

12. Cort, J. L. Age and growth of the bluefin tuna (Thunnus thynnus) of the Northeast Atlantic. In Domestication of the bluefin tuna Thunnus thynnus thynnus. Cahiers Options Méditerranéennes (CIHEAM) 45–49 (2003).

13. Mather, F. J., Mason, J. M. & Jones, A. C. Historical document: life history and fisheries of Atlantic bluefin tuna. (1995). NOAA Technical Memorandum NMFS-SEFSC – 370.

14. Puncher GN, et al. Spatial dynamics and mixing of bluefin tuna in the Atlantic Ocean and Mediterranean Sea revealed using next-generation sequencing. Mol. Ecol. Resour. 2018;18:620–638. doi: 10.1111/1755-0998.12764. [PubMed] [CrossRef] [Google Scholar]

15. Rodríguez-Ezpeleta N, et al. Determining natal origin for improved management of Atlantic bluefin tuna. Front. Ecol. Environ. 2019;17:439–444. doi: 10.1002/fee.2090. [CrossRef] [Google Scholar]

16. Richardson DE, et al. Discovery of a spawning ground reveals diverse migration strategies in Atlantic bluefin tuna (Thunnus thynnus) Proc. Natl. Acad. Sci. USA. 2016;113:3299–3304. doi: 10.1073/pnas.1525636113. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

17. Piccinetti, C., Di Natale, A. & Arena, P. Eastern bluefin tuna (Thunnus thynnus, L.) reproduction and reproductive areas and season. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA69, 891–912 (2013).

18. Cort, J. L. & Abaunza, P. The present state of traps and fisheries research in the strait of Gibraltar. In The Bluefin Tuna Fishery in the Bay of Biscay : Its Relationship with the Crisis of Catches of Large Specimens in the East Atlantic Fisheries from the 1960s (eds. Cort, J. L. & Abaunza, P.) 37–78 (Springer International Publishing, 2019).

19. Alemany, F., Tensek, S. & Pagà Garcia, A. ICCAT Atlantic-Wide Research programme for Bluefin Tuna (GBYP) activity report for the Phase 9 and the first part of Phase 10. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA77, 666–700 (2020).

20. MacKenzie, B. R. & Mariani, P. Spawning of bluefin tuna in the Black Sea: historical evidence, environmental constraints and population plasticity. PLoS ONE7, e39998 (2012). [PMC free article] [PubMed]

21. Di Natale, A. The Eastern Atlantic bluefin tuna: Entangled in a big mess, possibly far from a conservation red alert. Some comments after the proposal to include bluefin tuna in CITES Appendix I. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA 65(3), 1004–1043 (2010).

22. Worm B, Tittensor DP. Range contraction in large pelagic predators. Proc. Natl. Acad. Sci. USA. 2011;108:11942–11947. doi: 10.1073/pnas.1102353108. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

23. Fromentin J-M. Lessons from the past: Investigating historical data from bluefin tuna fisheries. Fish Fish. 2009;10:197–216. doi: 10.1111/j.1467-2979.2008.00311.x. [CrossRef] [Google Scholar]

24. Siskey MR, Wilberg MJ, Allman RJ, Barnett BK, Secor DH. Forty years of fishing: Changes in age structure and stock mixing in northwestern Atlantic bluefin tuna (Thunnus thynnus) associated with size-selective and long-term exploitation. ICES J. Mar. Sci. 2016;73:2518–2528. doi: 10.1093/icesjms/fsw115. [CrossRef] [Google Scholar]

25. ICCAT. Report of the 2020 second ICCAT intersessional meeting of the bluefin tuna species group. Online, 20–28 July 2020. SECOND BFT INTERSESSIONAL MEETING – ONLINE 2020 (2020).

26. Ravier C, Fromentin J-M. Long-term fluctuations in the eastern Atlantic and Mediterranean bluefin tuna population. ICES J. Mar. Sci. 2001;58:1299–1317. doi: 10.1006/jmsc.2001.1119. [CrossRef] [Google Scholar]

27. Garcia, A. P. et al. Report on revised trap data recovered by ICCAT GBYP from Phase 1 to Phase 6. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA73, 2074–2098 (2017).

28. Anderson CNK, et al. Why fishing magnifies fluctuations in fish abundance. Nature. 2008;452:835–839. doi: 10.1038/nature06851. [PubMed] [CrossRef] [Google Scholar]

29. Di Natale, A. & Idrissi, M. Factors to be taken into account for a correct reading of tuna trap catch series. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA67, 242–261 (2012).

30. Laconcha, U. et al. New nuclear SNP markers unravel the genetic structure and effective population size of Albacore Tuna (Thunnus alalunga). PLoS ONE10, e0128247 (2015). [PMC free article] [PubMed]

31. Speller, C. F. et al. High potential for using DNA from ancient herring bones to inform modern fisheries management and conservation. PLoS ONE7, e51122 (2012). [PMC free article] [PubMed]

32. Montes, I. et al. No loss of genetic diversity in the exploited and recently collapsed population of Bay of Biscay anchovy (Engraulis encrasicolus, L.). Mar. Biol.163, 98 (2016).

33. Chapman DD, et al. Genetic diversity despite population collapse in a critically endangered marine fish: The smalltooth sawfish (Pristis pectinata) J. Hered. 2011;102:643–652. doi: 10.1093/jhered/esr098. [PubMed] [CrossRef] [Google Scholar]

34. Hutchinson WF, van Oosterhout C, Rogers SI, Carvalho GR. Temporal analysis of archived samples indicates marked genetic changes in declining North Sea cod (Gadus morhua) Proc. Biol. Sci. 2003;270:2125–2132. doi: 10.1098/rspb.2003.2493. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

35. Ólafsdóttir GÁ, Westfall KM, Edvardsson R, Pálsson S. Historical DNA reveals the demographic history of Atlantic cod (Gadus morhua) in medieval and early modern Iceland. Proc. Biol. Sci. 2014;281:20132976. [PMC free article] [PubMed] [Google Scholar]

36. Bonanomi S, et al. Archived DNA reveals fisheries and climate induced collapse of a major fishery. Sci. Rep. 2015;5:15395. doi: 10.1038/srep15395. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

37. Nielsen EE, Hansen MM, Loeschcke V. Analysis of microsatellite DNA from old scale samples of Atlantic salmon Salmo salar : A comparison of genetic composition over 60 years. Mol. Ecol. 1997;6:487–492. doi: 10.1046/j.1365-294X.1997.00204.x. [CrossRef] [Google Scholar]

38. Johnson, B. M., Kemp, B. M. & Thorgaard, G. H. Increased mitochondrial DNA diversity in ancient Columbia River basin Chinook salmon Oncorhynchus tshawytscha. PLoS ONE13, e0190059 (2018). [PMC free article] [PubMed]

39. Bowles E, Marin K, Mogensen S, MacLeod P, Fraser DJ. Size reductions and genomic changes within two generations in wild walleye populations: associated with harvest? Evol. Appl. 2020;13:1128–1144. doi: 10.1111/eva.12987. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

40. Royle, T. C. A. et al. Investigating the sex-selectivity of a middle Ontario Iroquoian Atlantic salmon (Salmo salar) and lake trout (Salvelinus namaycush) fishery through ancient DNA analysis. J. Archaeol. Sci. Rep.31, 102301 (2020).

41. Therkildsen NO, et al. Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod. Mol. Ecol. 2013;22:2424–2440. doi: 10.1111/mec.12260. [PubMed] [CrossRef] [Google Scholar]

42. Pinsky, M. L. et al. Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic. Proc. Natl. Acad. Sci. USA118, (2021). [PMC free article] [PubMed]

43. Onar, V., Pazvant, G. & Armutak, A. Radiocarbon dating results of the animal remains uncovered at Yenikapi Excavations. In Istanbul Archaeological Museums, Proceedings of the 1st Symposium on Marmaray-Metro Salvage Excavations 249–256 (2008).

44. Bernal-Casasola D, Expósito JA, Díaz JJ. The Baelo Claudia paradigm: The exploitation of marine resources in Roman cetariae. J. Marit. Archaeol. 2018;13:329–351. doi: 10.1007/s11457-018-9209-z. [CrossRef] [Google Scholar]

45. Bernal, D. & Monclova, A. Pescar con Arte. Fenicios y romanos en el origen de los aparejos andaluces. Monografías del Proyecto Sagena3, (2011).

46. Puncher GN, et al. Comparison and optimization of genetic tools used for the identification of ancient fish remains recovered from archaeological excavations and museum collections in the Mediterranean region. Int J Osteoarchaeol. 2019;29:365–376. doi: 10.1002/oa.2765. [CrossRef] [Google Scholar]

Read More:  10 all of the following are necessary ingredients for cloud formation except: Ideas

47. Kemp BM, Smith DG. Use of bleach to eliminate contaminating DNA from the surface of bones and teeth. Forensic Sci. Int. 2005;154:53–61. doi: 10.1016/j.forsciint.2004.11.017. [PubMed] [CrossRef] [Google Scholar]

48. Dabney J, et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. USA. 2013;110:15758–15763. doi: 10.1073/pnas.1314445110. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

49. Serventi P, et al. Iron Age Italic population genetics: The Piceni from Novilara (8th–7th century BC) Ann. Hum. Biol. 2018;45:34–43. doi: 10.1080/03014460.2017.1414876. [PubMed] [CrossRef] [Google Scholar]

50. Star B, et al. The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011;477:207–210. doi: 10.1038/nature10342. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

51. Tine M, et al. European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat. Commun. 2014;5:5770. doi: 10.1038/ncomms6770. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

52. Chini V, et al. Genes expressed in bluefin tuna (Thunnus thynnus) liver and gonads. Gene. 2008;410:207–213. doi: 10.1016/j.gene.2007.12.012. [PubMed] [CrossRef] [Google Scholar]

53. Gardner LD, Jayasundara N, Castilho PC, Block B. Microarray gene expression profiles from mature gonad tissues of Atlantic bluefin tuna, Thunnus thynnus in the Gulf of Mexico. BMC Genomics. 2012;13:530. doi: 10.1186/1471-2164-13-530. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

54. Kamvar, Z. N., Tabima, J. F. & Grünwald, N. J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ2, e281 (2014). [PMC free article] [PubMed]

55. Team, R. C. R development core team. RA Lang. Environ. Stat. Comput.55, 275–286 (2013).

56. Paradis E. pegas: An R package for population genetics with an integrated–modular approach. Bioinformatics. 2010;26:419–420. doi: 10.1093/bioinformatics/btp696. [PubMed] [CrossRef] [Google Scholar]

57. Rousset F. genepop’007: A complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 2008;8:103–106. doi: 10.1111/j.1471-8286.2007.01931.x. [PubMed] [CrossRef] [Google Scholar]

58. Foll M, Gaggiotti O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics. 2008;180:977–993. doi: 10.1534/genetics.108.092221. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

59. Whitlock MC, Lotterhos KE. Reliable detection of loci responsible for local adaptation: Inference of a null model through trimming the distribution of FST. Am. Nat. 2015;186:S24–S36. doi: 10.1086/682949. [PubMed] [CrossRef] [Google Scholar]

60. Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. (1995).

61. Goudet J. hierfstat, a package for r to compute and test hierarchical F-statistics. Mol. Ecol. Notes. 2005;5:184–186. doi: 10.1111/j.1471-8286.2004.00828.x. [CrossRef] [Google Scholar]

62. Waples RS. A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci. Conserv. Genet. 2006;7:167–184. doi: 10.1007/s10592-005-9100-y. [CrossRef] [Google Scholar]

63. Do C, et al. NeEstimator v2: Re-implementation of software for the estimation of contemporary effective population size (Ne ) from genetic data. Mol. Ecol. Resour. 2014;14:209–214. doi: 10.1111/1755-0998.12157. [PubMed] [CrossRef] [Google Scholar]

64. Wang J, Santiago E, Caballero A. Prediction and estimation of effective population size. Heredity. 2016;117:193–206. doi: 10.1038/hdy.2016.43. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

65. Jenkins TL, Ellis CD, Triantafyllidis A, Stevens JR. Single nucleotide polymorphisms reveal a genetic cline across the north-east Atlantic and enable powerful population assignment in the European lobster. Evol. Appl. 2019;12:1881–1899. doi: 10.1111/eva.12849. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

66. Jombart T, et al. Package ‘adegenet’ Bioinform. Appl. Note. 2008;24:1403–1405. doi: 10.1093/bioinformatics/btn129. [PubMed] [CrossRef] [Google Scholar]

67. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–959. doi: 10.1093/genetics/155.2.945. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

68. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005;14:2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x. [PubMed] [CrossRef] [Google Scholar]

69. Earl, D. A. & vonHoldt, B. M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361 (2012).

70. Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 2015;15:1179–1191. doi: 10.1111/1755-0998.12387. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

71. Nei, M. Molecular Evolutionary Genetics. (Columbia University Press, 1987). 10.7312/nei-92038.

72. Frankham, R., Scientist Emeritus Jonathan, Briscoe, D. A. & Ballou, J. D. Introduction to Conservation Genetics. (Cambridge University Press, 2002).

73. Di Natale, A. Due to the new scientific knowledge, is it time to reconsider the stock composition of the Atlantic bluefin tuna? Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA75, 1282–1292 (2019).

74. Di Natale, A., Tensek, S. & Pagá García, A. ICCAT Atlantic-wide research programme for bluefin tuna (GBYP) activity report for the last part of phase and the first part of phase (2016–2017). https://www.iccat.int/Documents/CVSP/CV074_2017/n_6/CV074063100.pdf (2017).

75. Leonard JA. Ancient DNA applications for wildlife conservation. Mol. Ecol. 2008;17:4186–4196. doi: 10.1111/j.1365-294X.2008.03891.x. [PubMed] [CrossRef] [Google Scholar]

76. Alter, S. E., Newsome, S. D. & Palumbi, S. R. Pre-whaling genetic diversity and population ecology in eastern Pacific gray whales: Insights from ancient DNA and stable isotopes. PLoS ONE7, e35039 (2012). [PMC free article] [PubMed]

77. Cole TL, et al. Ancient DNA of crested penguins: Testing for temporal genetic shifts in the world’s most diverse penguin clade. Mol. Phylogenet. Evol. 2019;131:72–79. doi: 10.1016/j.ympev.2018.10.025. [PubMed] [CrossRef] [Google Scholar]

78. Dures SG, et al. A century of decline: Loss of genetic diversity in a southern African lion-conservation stronghold. Divers. Distrib. 2019;25:870–879. doi: 10.1111/ddi.12905. [CrossRef] [Google Scholar]

79. Thomas, J. E. et al. Demographic reconstruction from ancient DNA supports rapid extinction of the great auk. Elife8, (2019). [PMC free article] [PubMed]

80. Colson, I. & Hughes, R. N. Rapid recovery of genetic diversity of dogwhelk (Nucella lapillus L.) populations after local extinction and recolonization contradicts predictions from life-history characteristics. Mol. Ecol.13, 2223–2233 (2004). [PubMed]

81. McEachern MB, Van Vuren DH, Floyd CH, May B, Eadie JM. Bottlenecks and rescue effects in a fluctuating population of golden-mantled ground squirrels (Spermophilus lateralis) Conserv. Genet. 2011;12:285–296. doi: 10.1007/s10592-010-0139-z. [CrossRef] [Google Scholar]

82. Jangjoo M, Matter SF, Roland J, Keyghobadi N. Connectivity rescues genetic diversity after a demographic bottleneck in a butterfly population network. Proc. Natl. Acad. Sci. USA. 2016;113:10914–10919. doi: 10.1073/pnas.1600865113. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

83. Porch, C. E., Bonhommeau, S., Diaz, G. A., Haritz, A. & Melvin, G. The journey from overfishing to sustainability for Atlantic bluefin tuna, Thunnus thynnus. In The Future of Bluefin Tunas: Ecology, Fisheries Management, and Conservation 3–44 (2019).

84. Di Natale, A., Macias, D. & Cort, J. L. Atlantic bluefin tuna fisheries: temporal changes in the exploitation pattern, feasibility of sampling, factors that can influence our ability to understand spawning structure and dynamics. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA76, 354–388 (2020).

85. Viñas, J. & Tudela, S. A validated methodology for genetic identification of tuna species (genus Thunnus). PLoS ONE4, e7606 (2009). [PMC free article] [PubMed]

86. MacKenzie BR, Mosegaard H, Rosenberg AA. Impending collapse of bluefin tuna in the northeast Atlantic and Mediterranean. Conserv. Lett. 2009;2:26–35. doi: 10.1111/j.1755-263X.2008.00039.x. [CrossRef] [Google Scholar]

87. Collette BB. Bluefin tuna science remains vague. Science. 2017;358:879–880. [PubMed] [Google Scholar]

88. Nøttestad, L., Boge, E. & Ferter, K. The comeback of Atlantic bluefin tuna (Thunnus thynnus) to Norwegian waters. Fish. Res.231, 105689 (2020).

89. Lehodey P, et al. Climate variability, fish, and fisheries. J. Clim. 2006;19:5009–5030. doi: 10.1175/JCLI3898.1. [CrossRef] [Google Scholar]

90. Kuwae M, et al. Sedimentary DNA tracks decadal-centennial changes in fish abundance. Commun Biol. 2020;3:558. doi: 10.1038/s42003-020-01282-9. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

91. Domingues R, et al. Variability of preferred environmental conditions for Atlantic bluefin tuna (Thunnus thynnus) larvae in the Gulf of Mexico during 1993–2011. Fish. Oceanogr. 2016;25:320–336. doi: 10.1111/fog.12152. [CrossRef] [Google Scholar]

92. Reglero P, et al. Pelagic habitat and offspring survival in the eastern stock of Atlantic bluefin tuna. ICES J. Mar. Sci. 2019;76:549–558. doi: 10.1093/icesjms/fsy135. [CrossRef] [Google Scholar]

93. Faillettaz, R., Beaugrand, G., Goberville, E. & Kirby, R. R. Atlantic Multidecadal Oscillations drive the basin-scale distribution of Atlantic bluefin tuna. Sci. Adv.5, eaar6993 (2019). [PMC free article] [PubMed]

94. Hanke, A. et al. Stock mixing rates of bluefin tuna from Canadian landings: 1975–2015. Collect. Vol. Sci. Pap. ICCAT/Recl. Doc. Sci. CICTA/Colecc. Doc. Cient. CICAA74, 2622–2634 (2017).

95. Fraser DJ, et al. Comparative estimation of effective population sizes and temporal gene flow in two contrasting population systems. Mol. Ecol. 2007;16:3866–3889. doi: 10.1111/j.1365-294X.2007.03453.x. [PubMed] [CrossRef] [Google Scholar]

96. Albrechtsen A, Nielsen FC, Nielsen R. Ascertainment biases in SNP chips affect measures of population divergence. Mol. Biol. Evol. 2010;27:2534–2547. doi: 10.1093/molbev/msq148. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

97. Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R. Ascertainment bias in studies of human genome-wide polymorphism. Genome Res. 2005;15:1496–1502. doi: 10.1101/gr.4107905. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

98. Lachance J, Tishkoff SA. SNP ascertainment bias in population genetic analyses: Why it is important, and how to correct it. BioEssays. 2013;35:780–786. doi: 10.1002/bies.201300014. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

99. Hofreiter M, et al. The future of ancient DNA: Technical advances and conceptual shifts. BioEssays. 2015;37:284–293. doi: 10.1002/bies.201400160. [PubMed] [CrossRef] [Google Scholar]

100. Malomane DK, et al. Efficiency of different strategies to mitigate ascertainment bias when using SNP panels in diversity studies. BMC Genomics. 2018;19:22. doi: 10.1186/s12864-017-4416-9. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

101. Bradbury IR, et al. Evaluating SNP ascertainment bias and its impact on population assignment in Atlantic cod, Gadus morhua. Mol. Ecol. Resour. 2011;11:218–225. doi: 10.1111/j.1755-0998.2010.02949.x. [PubMed] [CrossRef] [Google Scholar]

102. Lou RN, Jacobs A, Wilder A, Therkildsen NO. A beginner’s guide to low-coverage whole genome sequencing for population genomics. Mol. Ecol. 2020 doi: 10.1111/mec.16077. [PubMed] [CrossRef] [Google Scholar]

103. Schlötterer C. Hitchhiking mapping–functional genomics from the population genetics perspective. Trends Genet. 2003;19:32–38. doi: 10.1016/S0168-9525(02)00012-4. [PubMed] [CrossRef] [Google Scholar]


Articles from Scientific Reports are provided here courtesy of Nature Publishing Group


Extra Information About which environmental factor could lead to a decrease in genetic variation in a population of tuna? That You May Find Interested

If the information we provide above is not enough, you may find more below here.

Ancient DNA SNP-panel data suggests stability in bluefin tuna …

Ancient DNA SNP-panel data suggests stability in bluefin tuna ...

  • Author: ncbi.nlm.nih.gov

  • Rating: 4⭐ (753753 rating)

  • Highest Rate: 5⭐

  • Lowest Rate: 1⭐

  • Sumary: Atlantic bluefin tuna (Thunnus thynnus; BFT) abundance was depleted in the late 20th and early 21st century due to overfishing. Historical catch records further indicate that the abundance of BFT in the Mediterranean has been fluctuating since …

  • Matching Result: by AJ Andrews · 2021 · Cited by 2 — The current study investigates genetic variability in eastern … or that overfishing did trigger a severe population decline but BFT is …

  • Intro: Ancient DNA SNP-panel data suggests stability in bluefin tuna genetic diversity despite centuries of fluctuating catches in the eastern Atlantic and Mediterranean Journal List Scientific Reports PMC8528830 Sci Rep. 2021; 11: 20744. Adam J. Andrews,#1,2 Gregory N. Puncher,#1,3 Darío Bernal-Casasola,4 Antonio Di Natale,5 Francesco Massari,1 Vedat Onar,6 Nezir Yaşar Toker,6…
  • Source: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8528830/

Evolution by Natural Selection – University of Hawaii at Manoa

Evolution by Natural Selection - University of Hawaii at Manoa

  • Author: manoa.hawaii.edu

  • Rating: 4⭐ (753753 rating)

  • Highest Rate: 5⭐

  • Lowest Rate: 1⭐

  • Sumary: Without physically counting bacteria or viruses, scientists estimate that there are more than 8.7 million species of organisms living on Earth today. The origin and extinction of so many species has fascinated scientists for thousands of years, since the days…

  • Matching Result: This led him to develop the theory of natural selection, which is the best … Sexual reproduction can also increase genetic variation in a population.

  • Intro: Evolution by Natural Selection | manoa.hawaii.edu/ExploringOurFluidEarthIntroduction to Evolution by Natural Selection Without physically counting bacteria or viruses, scientists estimate that there are more than 8.7 million species of organisms living on Earth today. The origin and extinction of so many species has fascinated scientists for thousands of years, since the…
  • Source: https://manoa.hawaii.edu/exploringourfluidearth/biological/what-alive/evolution-natural-selection

Frequently Asked Questions About which environmental factor could lead to a decrease in genetic variation in a population of tuna?

If you have questions that need to be answered about the topic which environmental factor could lead to a decrease in genetic variation in a population of tuna?, then this section may help you solve it.

What leads to a decline in genetic diversity?

Fragmented and threatened populations are typically exposed to these conditions, which are likely to increase their risk of extinction (Saccheri et al.

Which environmental factor in a population is most likely to increase genetic diversity?

Genetic diversity in a population would most likely increase in response to an increase in food and a decrease in predators.

What kind of environmental element might affect natural selection?

The genetic and chemical diversity of plant populations is influenced by environmental factors such as soil nutrients, temperature, water availability, and light intensity. These environmental factors can exert strong selection pressures and may even influence the evolutionary path of plant populations [3,4].

What elements influence genetic variability?

Mutation, which can produce completely new alleles in a population, random mating, random fertilization, and recombination between homologous chromosomes during meiosis (which rearranges alleles within an organism’s progeny) can all contribute to genetic variation.

Which of the subsequent reduces genetic diversity?

Therefore, stabilizing natural selection is the right response.

Which population’s genetic diversity has decreased?

Because some versions of a gene can be lost due to random chance, and this is more likely to occur when populations are small, small populations tend to lose genetic diversity more quickly than large populations due to stochastic sampling error (i.e., genetic drift).

Which two factors are most responsible for genetic variation in populations?

Gene flow, or the movement of genes between different groups of organisms, mutations, which are changes in the DNA’s gene sequences, and sexual reproduction, which results in the creation of new gene combinations, are the other two sources of genetic variation.

How does genetic variation interact with the environment?

Some harmful environmental factors can change a gene’s nucleotide sequence, such as UV radiation, which can break DNA strands, and either of these effects can change the proteins that are made from a gene, which in turn affects traits.

What is the most significant environmental factor that can influence an animal’s genetic potential?

Environmental elements can affect which genes in an animal are expressed, which in turn affects the animal’s phenotype. These environmental factors include diet, temperature, oxygen levels, humidity, light cycles, and the presence of mutagens.

Which option chooses the best answer to reduce genetic variation?

Answer and explanation: The right response is (4). Gene flow is random and increases genetic variation, whereas artificial selection is selective and decreases it.

What three elements of the environment have the potential to impact gene expression?

By altering the way molecules bind to DNA or the structure of proteins that DNA wraps around, environmental factors such as food, drugs, or exposure to toxins can cause epigenetic changes.

Share this post