Next Article in Journal
Two Seas for One Great Diversity: Checklist of the Marine Heterobranchia (Mollusca; Gastropoda) from the Salento Peninsula (South-East Italy)
Next Article in Special Issue
A Mediterranean Monk Seal Pup on the Apulian Coast (Southern Italy): Sign of an Ongoing Recolonisation?
Previous Article in Journal
Final Note of Special Issue “Tardigrades Taxonomy, Biology and Ecology”
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 12
Issue 5
10.3390/d12050170
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Loss of Mitochondrial Genetic Diversity in Overexploited Mediterranean Swordfish (Xiphias gladius, 1759) Population

by
Tommaso Righi
,
Andrea Splendiani
,
Tatiana Fioravanti
,
Elia Casoni
,
Giorgia Gioacchini
,
Oliana Carnevali
and
Vincenzo Caputo Barucchi
*
Dipartimento di Scienze della Vita e dell’Ambiente, Università Politecnica delle Marche, Via Brecce Bianche, 60131 Ancona, Italy
*
Author to whom correspondence should be addressed.
Diversity 2020, 12(5), 170; https://0-doi-org.brum.beds.ac.uk/10.3390/d12050170
Submission received: 4 March 2020 / Revised: 13 April 2020 / Accepted: 24 April 2020 / Published: 26 April 2020
(This article belongs to the Special Issue Historical and Current Diversity Patterns of Mediterranean Marine Species)
Browse Figures
Versions Notes

Abstract

:
Intense and prolonged mortality caused by over-exploitation could drive the decay of genetic diversity which may lead to decrease species’ resilience to environmental changes, thus increasing their extinction risk. Swordfish is a high commercial value species, especially in the Mediterranean Sea, where it is affected by high catch levels. Mediterranean swordfish consist of a population genetically and biologically distinct from Atlantic ones and therefore managed as a separate stock. The last Mediterranean swordfish stock assessment reported that in the last forty years Mediterranean swordfish has been overfished and, to date, it is still subject to overfishing. A comparison between an available mitochondrial sequence dataset and a homologous current sample was carried out to investigate temporal genetic variation in the Mediterranean swordfish population over near twenty years. Our study provides the first direct measure of reduced genetic diversity for Mediterranean swordfish during a short period, as measured both in the direct loss of mitochondrial haplotypes and reduction in haplotype diversity. A reduction of the relative females’ effective population size in the recent sample has been also detected. The possible relationship between fishery activities and the loss of genetic diversity in the Mediterranean swordfish population is discussed.

1. Introduction

The swordfish (Xiphias gladius) is a large pelagic and migratory fish found in the open waters of all oceans, including the Mediterranean Sea [1]. Swordfish migration is complex and multi-directional [1]. Seasonally adults swordfish migrate for reproductive reasons to spawning grounds, while a diel pattern in vertical movement occurring for feeding, reaching deep water (300–600 m) during daylight and staying closer to the sea-surface during the night [2,3]. Swordfish growth is relatively rapid and sexually dimorphic. Females generally reach a larger maximum size than males and large individuals are usually female [2,4]. Mediterranean swordfish consist of a population genetically and biologically distinct from the Atlantic one. The growth parameters are different, and the sexual maturity is reached at younger ages than in the Atlantic [4]. In the Mediterranean, the estimated size at which 50% (L50) of the female population is mature occurs from 131.5 cm [5] and 140 cm [6], corresponding to 3–4-year-old fish [2]. Males reach sexual maturity at smaller sizes and mature specimens have been found at about 90 cm (Lower jaw fork length, LJFL) [7]. Genetic differentiation [8], as well as several differences in biological parameters [2], have been used to identify and manage the Atlantic and Mediterranean swordfish as separated stocks [7]. A popular genetic marker used to examine swordfish populations is the mitochondrial control region (D-loop). As reported for other large pelagic species (Thunnus alalunga [9]; Thunnus thynnus [10] Seriola dumerili [11]) the mitochondrial control region of swordfish can be divided into two divergent clades. Specifically for swordfish, a ubiquitous Clade I widespread in all stocks, and Clade II which decreases in abundance with distance from the Mediterranean Sea [10,12,13] and is absent from the Pacific Ocean [14,15,16]. These phylogeographic associations support the inferred pattern of unidirectional historical gene flow from the Indo–Pacific into the Atlantic and the Atlantic origin of Clade II [10]. A subdivision in two monophyletic subclades to the Atlantic and Mediterranean Clade II lineage confirmed genetic isolation between the two stocks [10,13]. Furthermore, within the Mediterranean Sea Swordfish population sub-structuring was suggested based on a clinal decrease in mitochondrial DNA control region intra-clade genetic variability from Western and the Eastern basins [17]. The Mediterranean population showed the lowest levels of genetic variation compared to any other population worldwide [10,12,13,18,19]. According to Bremer et al. [10], the reduced level of variation in Mediterranean swordfish mtDNA is congruent with multiple faunal extinctions and subsequent re-colonization events occurred during the Pleistocene. However, according to Kotoulas et al. [20], the low genetic variability observed using nuclear markers could be related to the small effective population size of Mediterranean swordfish, probably caused by fishing pressure [19]. Moreover, according to Arocha [4], prolonged exploitation of swordfish in the Mediterranean area may have reduced the numbers of older specimens and consequently explain the lower mean size of the catches and lower size of maturation observed for the Mediterranean population.
Swordfish is a high commercial value species, especially in the Mediterranean Sea, where despite this basin represents less than 10% of the swordfish global range, catch levels are relatively high and similar to those of larger areas such as the North Atlantic [21]. Some characteristics as the presence of a larger spawning area in relation to the area of distribution of the stock, lower abundance of predators and higher recruitment were suggested as the determining factors explaining the higher abundance of swordfish in the Mediterranean Sea. However, in the last forty years, the swordfish stock of the Mediterranean Sea has been overfished and, to date, is still subject to overfishing [7]. The Spawning Stock Biomass (SSB) estimated represents less than 15% of the maximum sustainable yield (BMSY) and the fish mortality caused by harvesting (F) is almost twice the maximum sustainable yield. The 50–70% of the total yearly catches is represented by fish of small size often less than 3 years old, with a high level of immature swordfish reported [7]. The loss of genetic variation caused by population reduction is an important concern in fisheries management. Intense and prolonged mortality caused by overexploitation could drive the decay of genetic diversity across a wide range of marine fishes [22]. Genetic diversity has been defined as the variety of alleles and genotypes present in a population and that is reflected in morphological, physiological and behavioral differences between individuals and populations [23]. The reduction of genetic diversity may lead to decreased species resilience to environmental changes, thus increasing their extinction risk [23,24]. Small and isolated populations are more vulnerable to genetic degradation [23]. Moreover, in marine fishes, effective population size (Ne) can be several orders of magnitude smaller than census size (N) [25,26,27] and species may suffer a loss of genetic diversity under fishing pressure despite large census size.
The most powerful way to detect loss of genetic variation in exploited populations is the examination of genetic samples collected over time [28]. In this study, a dataset of available mitochondrial sequences of Mediterranean swordfish sampled in the 1990s and early 2000s has been compared with homologous mtDNA data obtained from current samples (2016–2018). This comparison was carried out to investigate if a temporal genetic variation occurred in the Mediterranean swordfish, where this species is considered threatened [21] and overexploited.

2. Materials and Methods

2.1. Sampling and DNA Extraction

A total of 287 swordfish were collected from six areas within the eastern, central and western Mediterranean regions, trying to obtain a representative coverage of the basin (For details see Table 1 and Figure 1).
Samples were achieved at the fishing landing of commercial catch by longline or trap bycatch (only in case of Sardinian samples) from May to October in three years 2016–2018. For each individual, a piece of the caudal fin was collected and was stored in ethanol absolute and kept at −20 °C until DNA extraction. Total genomic DNA was extracted using specific cartridge 401 in the MagCore® automated Nucleic Acid extractor (RBC Bioscience, New Taipei City, Taiwan) following the specific protocol. Historical mitochondrial data were obtained from a re-analysis of the dataset available from Viñas [17] (Table S2). That work included specimens sampled in several Mediterranean areas (Figure 1) from 1992 to 1995 and the early 2000s (see [11] for sample details), allowing us to evaluate the mitochondrial genetic variability between two complementary sample groups of Mediterranean swordfish over a 20-year period.

2.2. Mitochondrial DNA Amplification and Analysis

PCR-RFLP (Restrictin Fragment Length Polymorphism) analysis was performed to identify the two mitochondrial lineages. The restriction enzymes were selected by virtual restriction analysis (see Supplementary material S3 for virtual analysis details). A 360 bps long portion of the hypervariable L-domain of the mitochondrial control region was amplified using primers L15998 (5’-TAC CCC AAA CTC CCA AAG CTA-3’) [12] and SWO 5’ CCC TGT GAA ATA TGC TGG TTG 3’ (designed in this study). The amplification was carried out in 25 μL reaction volume containing: 1 × MyTaq reaction buffer (BioLine GmbH, Luckenwalde, Germany), 10 pmol of each primer, 0.5 U of BIOTAQ DNA-Polymerase (BioLine GmbH, Luckenwalde, Germany) and approximately 80 ng of the isolated DNA. For thermal cycler Viñas et al. [9] has been followed. Amplicons were double-digested with selected PacI and VspI restriction enzymes (Thermo Fisher Scientific, Waltham, Massachusetts, MA, USA). The digestion mix reaction was carried out in a total volume of 20μL with 4 μL of amplified DNA, 2 μL of Buffer G 10X (Thermo Fisher Scientific, Waltham, Massachusetts, MA, USA), and 0.3 μL [10 U/μL] of each enzyme. The digestion took place at 37 °C overnight. Restriction profiles were visualized in 2% agarose gel. Subsequently, to screen mitochondrial genetic variability a PCR-SSCP (Single Strand Conformation Polimorphism) analysis was performed. 5 μL amplicon was added to 4 μL loading buffer (98% formamide, 10 mM, EDTA (0.5 M, pH 8), 0.05% bromophenol blue, 0.05% xylene cyanol), heated to 95 °C for 5 min and immediately chilled on ice. Vertical electrophoresis run was performed in a nondenaturing polyacrylamide gel (8% acrylamide/polyacrylamide [49:1], 10% glycerol) at 5 W for 12 h at room temperature, with 1 × TBE as the running buffer. To identify different haplotypes Sanger sequencing of D-loop was performed on a subsample for each SSCP pattern on an ABI 3730XL DNA (Applied Biosystems, Foster City, CA, USA). To validate the sensitivity of SSCP technique in the identification of different haplotypes, all individuals belonging to the Clade II lineage, which showed a fewer number of SSCP morphs, were sequenced. Sequences were aligned using ClustalW [29] and checked by eye in BioEdit [30]. Analysis of sequence variability was analyzed pooling all samples together and, following Viñas [17] grouping samples in Western (Balearic, Sardinian, Tyrrhenian) and Eastern (Sicilian, Adriatic and Aegean), according to the basin of origin. However, no genetic differentiation was detected by AMOVA (Analysis of Molecular VAriance) between the western and eastern group therefore we reported the only comparison among whole historical and contemporary datasets.
The number of haplotypes (N), and polymorphic site, haplotypic (h) and nucleotide diversities (π) were computed in Arlequin v. 3.5 [31]. The analysis was performed for both dataset pooling Clades I and II together, and for each clade separately. To account for the sampling variance of haplotypic richness, rarefaction curves [32] were generated for both samples. Rarefaction curves were generated using iNEXT [33,34]. Median-joining (MJ) networks for the entire dataset were built using the software Network v 4.510 [35] with the default settings.
The extent of genetic differentiation between historical and contemporary datasets was investigated with an analysis of molecular variance, AMOVA [36], and the exact tests of population differentiation implemented in Arlequin [31]. The AMOVA significance levels were determined by 10,000 permutations while the exact test was conducted by using 10,000 Markov chain steps and 1000 dememorization steps.Finally, temporal changes in the relative effective number of females (Nef) was estimated comparing θ for each dataset. The θ is expected to be equal to 2 Nef μ for neutral mutations in mtDNA, where μ is the mutation rate per generation and Nef is the effective population size of females. Considering μ to be constant between two temporal samples, we can treat differences in θ among samples as differences in relative female effective population sizes. Population size parameter (θ) was inferred for each dataset using the maximum likelihood (ML) inference in LAMARC 2.1 [37] using Watterson’s estimator as starting points, default model of evolution was used. A Markov Chain Monte Carlo was run for ten short chains and three long chains with 1000 and 10,000 recorded genealogies, respectively, after discarding the first 10,000 genealogies. One of every 20 reconstructed genealogies was sampled for both short and long chains.

3. Results

The SSCP analysis of 360 bps portion of mtDNA D-loop from 287 contemporary Mediterranean swordfish revealed 36 different morphs. For each morph, a subsample has been sequenced for a total of 172 swordfish, corresponding to 60% of the total sample. The sequencing of all individuals belonging to the Clade II corroborated the SSCP profile interpretation. The sequence length ranged from 290 to 297 bps due to the presence of 1–3 tandemly repeated copies of the motif TACA near the 5’ end. Sequence alignment revealed 49 polymorphic sites, defining 36 distinct haplotypes (Table 2).
According to Bremer et al. [12] the mtDNA haplotypes of swordfish clustered into two highly divergent clades, Clade I and Clade II. The MJ network (Figure 2) identified the two mtDNA clades.
Clade I included four centroids (Genbank n. AY650768, AY650763, AY650778, AY650781) represented respectively by 61, 55, 14 and 25 fish, while Clade II was characterized by a single star-like formation featuring one major centroid (AY650762) with a total of 73 individuals. These five centroids haplotypes were observed with high frequencies in all Mediterranean localities investigated in this study, as already observed in previous works [10,17]. The other haplotypes, instead, were represented in very few numbers of the individuals often representing singletons. In this study seventeen, previously unidentified haplotypes were sampled. Fifteen of which belong to Clade I and two belong to Clade II and all were observed at low frequencies (Table S1). Sequences are available from the NCBI database with accession numbers MN652595-MN652611.
Comparison among historical dataset from Viñas et al. [17] and the present dataset, respectively composed by 251 and 287 swordfish, revealed the same percentage of individuals belonging to Clade I (about 66%) and Clade II (about 33%). However, dataset comparison indicates a reduction in the number of haplotypes from 93 in the historical to 36 in the contemporary samples (Table 3).
Haplotype reduction was confirmed by the rarefaction curve that reported the highest haplotype richness values in the older sample (Figure 3).
Considering the two clades separately, the observed haplotype number decreased from 60 to 26 for Clade I and from 33 to 10 for Clade II. Nineteen haplotypes were shared between datasets (11 for Clade I and 8 for Clade II). Moreover, in the modern sample, a shift towards the fixation of the main five haplotypes (i.e., centroids) was also observed (Figure 2, Table S2).
A reduction of haplotype diversity was observed between historical (h = 0.947 + 0.007 s.d.) and modern datasets (h = 0.844 + 0.01 s.d.). Haplotype diversity decrease was more evident considering the two clades separately. For Clade I the diversity in historical samples (h = 0.91) was greater than the contemporary (h = 0.793) as well as for Clade II where diversity was h = 0.861 in the historical sample and h = 0.418 in the modern sample. Levels of nucleotide diversity were similar throughout the two databases. All estimates of haplotype and nucleotide diversity are given in Table 3. The AMOVA analysis identifies significant genetic differentiation between older and recent datasets. Low levels of differentiation (FST = 0.018, p < 0.001) were detected between samples pooling clades and considering Clade I singularly, while historical and modern samples were moderately differentiated (FST = 0.121, p < 0.001) considering only Clade II (Table 4).
The older sample significantly differed from the modern one when compared through an exact test of sample differentiation. The effective population size parameters of swordfish populations were estimated for each mtDNA CR-I clade separately. The comparison of the θ values showed a reduction of the relative females’ effective population size in the recent sample (θ = 0.033, 95% C.I. 0.028–0.044) compared to the historical one (θ = 0.142, 95% C.I. 0.121–0.188) (Table 3). Ratios between the historical and current θ suggested a reduction of approximately three-quarters of the female population size. Estimates were congruent considering both the clades separately and pooled them together.

4. Discussion

In this study, for the first time, a temporal approach has been used to study recent changes in mitochondrial genetic diversity in the threatened Mediterranean swordfish. Comparison between temporal samples allows us to identify short-term changes that happened in the last few decades. Our results revealed that Mediterranean swordfish has undergone genetic depletion during the last twenty years.
The number of mitochondrial haplotypes observed in the recent sample (36 in total; 26 Clade I and 10 Clade II) was largely reduced than in previous studies [10,17]. The PCR-SSCP is a valuable tool to detect nucleotidic variation in DNA because even a single nucleotide change may alter the conformation of ssDNA. Despite this high resolving power, there is evidence that not all nucleotide modifications result in conformation changes [38]. To exclude loss of haplotypes due to the SSCP analysis we sequenced all the individuals belonging to the Clade II and all sequences confirmed the SSCP interpretation. Therefore, we consider the SSCP technique used in this paper sensitive enough to detect each single nucleotide polymorphism.
The temporal comparison based on the mtDNA genetic variability highlighted as approximately one-third of the haplotypes was lost in the last twenty years. This reduction was noticeable as a lower haplotype diversity observed on modern samples (h = 0.84, this study) compared to previous ones (h = 0.94) [10,17]. When the two mitochondrial clades were analysed separately, the reduction of genetic diversity was not uniformly distributed between clades and the reduction in genetic diversity was more evident for the Clade II. Diversity comparison indicates approximately half less haplotype diversity in the contemporary Clade II samples than in the historical. Genetic differentiation between contemporary and historical swordfish was corroborated by both AMOVA and exact test results (Table 4).
The reduction in haplotype diversity with an increment in the proportion of identical haplotype and the loss of rare haplotypes observed in current Mediterranean swordfish suggest that population has undergone a demographic reduction [28]. As far as it’s known, nowadays, the Mediterranean swordfish stock is considered threatened by overfishing [7]. In fact, the spawning stock biomass (SSB) estimated is less than 15% of the maximum sustainable yield (BMSY). From 2007–2010 the mesopelagic longline was introduced in Mediterranean swordfish fishery, substituting the traditional surface longline [7]. The two types of gears differ mainly for the depth of displacement, with the mesopelagic reaching deeper depths (50–800 m) than the traditional longline (15–60 m). Increment in fish size catch using the mesopelagic longline has revealed the presence of a bulk of large spawners, that found their refuge in the deep [39]. However, after only six years the mean swordfish size exploited shifted towards smaller sizes, in a situation very similar to that recorded in the past for surface driftnets and surface longline [40]. Thus, it is possible to assume that the mesopelagic gear catching bigger swordfish may have increased the impact on adults, affecting the spawning stock [39,40,41]. Specifically, mtDNA is maternally inherited therefore, the analysis of mtDNA variability reflects only the history of the Mediterranean swordfish female portion. Swordfish females reach larger sizes than males, outnumbering males in larger size classes (LJFL > 140 cm) [5,39,42]. Therefore, an increase in the average size may be related to an increase in the number of mature females exploited. Furthermore, females mature at larger age and size than males, about two years later. The high percentage (50–70%) of small fish, often still immature, reported by the annual catch estimates, may affect the number of recruited new females in the spawning stock. Comparison of θs values supports the scenario of a recent population size decrease in Mediterranean swordfish stock, suggesting a reduction of approximately three-quarters of effective female population size (Table 3). Thus, the fishery may have altered the sex ratio, reducing the number of females in the breeding population, therefore leading to the loss in mitochondrial diversity.
Overexploitation and environmental degradation were identified as the main causes of reductions in stocks and extinction of marine species [43,44]. Human activities may reduce the genetic variability of the population in an extremely short period [45]. Overexploitation has been related to a decline in genetic diversity across a wide range of marine fish species [22]. The genetic diversity estimated in harvested and unharvested fish stocks over 140 species, resulted to be lower in the exploited than in the non-harvest stocks [22]. However, a limited number of studies have assessed the temporal reduction of effective population size (Ne) and the loss of genetic variation related to harvesting due to lack of historical data as well as samples. Some of these studies, using DNA extracted from archived otoliths or scales, found a loss of neutral genetic diversity in some fish species. A significant reduction in genetic diversity was described in New Zealand Hoplostethus atlanticus over only 6 years, during which time the exploitation reduced the biomass by 70% [46]. Marked genetic changes were also detected using microsatellite in the Gadus morhua Flamborough Head population across a period during which the population exhibited a decline in SSB related to high levels of exploitation [27]. Hauser [25] detected a significant decline in genetic diversity in a New Zealand snapper population over the 50 years since the onset of exploitation. Moreover, in Adriatic European anchovy (Engraulis encrasicolus) a significant reduction of mean genetic parameters have been detected and related to the stock collapse in 1987 [47]. The number of fish in a population (census population size, N) is often much larger than the fish that are reproducing and helping to maintain genetic diversity (the genetically effective population size, Ne). Hauser [25] using the Ne/N estimate in snapper suggested that fish stocks of several million individuals may be in danger of losing genetic variability in the long term. Considering that, globally 31% of fish populations are fished unsustainably, with an additional 60% fully fished [48], fishing may have already caused a considerable loss of overall biodiversity.
This study reports the first direct measure of reduction in genetic diversity for Mediterranean swordfish during a short period, as measured both in the direct loss of mitochondrial haplotypes and reductions in haplotype diversity. Genetic changes observed in this study suggest that reduction in the SSB of Mediterranean swordfish caused significant changes in genetic composition. The significant loss of mitochondrial diversity in Mediterranean swordfish over a short period is alarming because the rapid loss of genetic diversity has been shown to have harmful effects. Reduction in genetic diversity, coupled with low population size [20] and reproductive isolation [2], can disfavour Mediterranean swordfish recovery. This may result in genetic drift and inbreeding depression and reduction in fitness limiting long-term adaptability [23,24], particularly if abundance remains low and diversity continues to decay [22]. This result underlines the urgent necessity to re-evaluate management strategies of the swordfish in the Mediterranean Sea.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-2818/12/5/170/s1, Table S1. Haplotype frequencies observed have been reported for each sample locality, Table S2. GenBank accession numbers of the mtDNA CR haplotypes, clade assignment and frequencies observed in each temporal sample, Table S3. Predicted fragment sizes (bp) of a 290 bps long portion of the mitochondrial control region virtually double-digested with VspI and PacI, Table S4. Variable sites in the control-region segment of all 110 haplotypes of swordfish mtDNA detected in Viñas et al., 2010 and this work (in bold). Numbers above sites refer to nucleotide position on the L-strand of the control region of the swordfish. Identity with sequence 1 is indicated by dots and deletions of 1 residue with dashes.

Author Contributions

Conceptualization, T.R.; Methodology, T.R.; Formal analysis, T.R.; Investigation, T.R.; Resources, G.G. and V.C.B.; Data curation, E.C., T.R.; Writing—Original draft preparation, T.R.; Writing—Review and editing, A.S., T.F. and T.R.; Visualization, T.R.; Supervision, V.C.B.; Project administration, G.G.; Funding acquisition, O.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Ministry of Agriculture, Food and Forestry Policies (MIPAAF), note 6775, Art.36 Paragraph 1 Reg (UE9 n 508/2014) to O.C.

Acknowledgments

We truly thank the fishermen, Pesca Pronta Import-Export s.r.l. for their extremely precious collaboration during sampling activities. The authors wish to thank all the staff at OCEANIS s.r.l. (Ercolano, Naples) for the support in sampling activities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Palko, B.J.; Beardsley, G.L.; Richards, W.J. Synopsis of the biology of the swordfish, Xiphias gladius (L.). U.S. Dep. Commer. NOAA Tech. Rep. NMFS Circ. 1981, 441, 21. [Google Scholar]
  2. Neilson, J.; Arocha, F.; Cass-Calay, S.; Mejuto, J.; Ortiz, M.; Scott, G.; Smith, C.; Travassos, P.; Tserpes, G.; Andrushchenko, I. The recovery of Atlantic swordfish: The comparative roles of the regional fisheries management organization and species biology. Rev. Fish. Sci. 2013, 21, 59–97. [Google Scholar] [CrossRef]
  3. Abascal, F.J.; Mejuto, J.; Quintans, M.; García-Cortés, B.; Ramos-Cartelle, A. Tracking of the broadbill swordfish, Xiphias gladius, in the central and eastern North Atlantic. Fish. Res. 2015, 162, 20–28. [Google Scholar] [CrossRef]
  4. Arocha, F. Swordfish Reproduction in the Atlantic Ocean: An Overview. GCR 2007, 19. [Google Scholar] [CrossRef] [Green Version]
  5. Marisaldi, L.; Basili, D.; Candelma, M.; Sesani, V.; Pignalosa, P.; Gioacchini, G.; Carnevali, O. Maturity assignment based on histology-validated macroscopic criteria: Tackling the stock decline of the Mediterranean swordfish (Xiphias gladius). Aquat. Conserv. Mar. Freshw. Ecosyst. 2020, 30, 303–314. [Google Scholar] [CrossRef]
  6. Macías, D.; Hattour, A.; De la Serna, J.M.; Gómez-Vives, M.J.; Godoy, D. Reproductive characteristics of swordfish (Xiphias gladius) caught in the southwestern Mediterranean during 2003. Collect. Vol. Sci. Pap. ICCAT 2005, 58, 454–469. [Google Scholar]
  7. Anon. Report for Biennial Period, 2018–19 PART I (2018)—Volume 2; International Commission for the Conservation of Atlantic Tunas: Madrid, Spain, 2019; pp. 182–192. [Google Scholar]
  8. Bremer, J.A.; Mejuto, J.; Gómez-Márquez, J.; Pla-Zanuy, C.; Viñas, J.; Marques, C.; Hazin, F.; Griffiths, M.; Ely, B.; García-Cortés, B. Genetic population structure of Atlantic swordfish: Current status and future directions. Collect. Vol. Sci. Pap. ICCAT 2007, 61, 107–118. [Google Scholar]
  9. Viñas, J.; Bremer, J.A.; Pla, C. Inter-oceanic genetic differentiation among albacore (Thunnus alalunga) populations. Mar. Biol. 2004, 145, 225–232. [Google Scholar] [CrossRef]
  10. Bremer, J.R.A.; Viñas, J.; Mejuto, J.; Ely, B.; Pla, C. Comparative phylogeography of Atlantic bluefin tuna and swordfish: The combined effects of vicariance, secondary contact, introgression, and population expansion on the regional phylogenies of two highly migratory pelagic fishes. Mol. Phylogenetics Evol. 2005, 36, 169–187. [Google Scholar] [CrossRef]
  11. Šegvić-Bubić, T.; Marrone, F.; Grubišić, L.; Izquierdo-Gomez, D.; Katavić, I.; Arculeo, M.; Brutto, S.L. Two seas, two lineages: How genetic diversity is structured in Atlantic and Mediterranean greater amberjack Seriola dumerili Risso, 1810 (Perciformes, Carangidae). Fish. Res. 2016, 179, 271–279. [Google Scholar] [CrossRef]
  12. Bremer, J.R.A.; Mejuto, J.; Greig, T.W.; Ely, B. Global population structure of the swordfish (Xiphias gladius L.) as revealed by analysis of the mitochondrial DNA control region. J. Exp. Mar. Biol. Ecol. 1996, 197, 295–310. [Google Scholar] [CrossRef]
  13. Rosel, P.E.; Block, B.A. Mitochondrial control region variability and global population structure in the swordfish, Xiphias gladius. Mar. Biol. 1996, 125, 11–22. [Google Scholar] [CrossRef]
  14. Chow, S.; Okamoto, H.; Uozumi, Y.; Takeuchi, Y.; Takeyama, H. Genetic stock structure of the swordfish (Xiphias gladius) inferred by PCR-RFLP analysis of the mitochondrial DNA control region. Mar. Biol. 1997, 127, 359–367. [Google Scholar] [CrossRef]
  15. Muths, D.; Grewe, P.; Jean, C.; Bourjea, J. Genetic population structure of the Swordfish (Xiphias gladius) in the southwest Indian Ocean: Sex-biased differentiation, congruency between markers and its incidence in a way of stock assessment. Fish. Res. 2009, 97, 263–269. [Google Scholar] [CrossRef] [Green Version]
  16. Reeb, C.A.; Arcangeli, L.; Block, B.A. Structure and migration corridors in Pacific populations of the Swordfish Xiphius gladius, as inferred through analyses of mitochondrial DNA. Mar. Biol. 2000, 136, 1123–1131. [Google Scholar] [CrossRef]
  17. Viñas, J.; Pérez-Serra, A.; Vidal, O.; Alvarado Bremer, J.R.; Pla, C. Genetic differentiation between eastern and western Mediterranean swordfish revealed by phylogeographic analysis of the mitochondrial DNA control region. ICES J. Mar. Sci. 2010, 67, 1222–1229. [Google Scholar] [CrossRef] [Green Version]
  18. Bremer, J.R.A.; Baker, A.J.; Mejuto, J. Mitochondrial DNA control region sequences indicate extensive mixing of swordfish (Xiphias gladius) populations in the Atlantic Ocean. Can. J. Fish. Aquat. Sci. 1995, 52, 1720–1732. [Google Scholar] [CrossRef]
  19. Pujolar, J.M.; Roldan, M.I.; Pla, C. A genetic assessment of the population structure of swordfish (Xiphias gladius) in the Mediterranean Sea. J. Exp. Mar. Biol. Ecol. 2002, 276, 19–29. [Google Scholar] [CrossRef]
  20. Kotoulas, G.; Mejuto, J.; Antoniou, A.; Kasapidis, P.; Tserpes, G.; Piccinetti, C.; Peristeraki, P.; Garcia-Cortes, B.; Oikonomaki, K.; De la Serna, J.M. Global genetic structure of swordfish (Xiphias gladius) as revealed by microsatellite DNA markers. Collect. Vol. Sci. Pap. ICCAT 2007, 61, 79–88. [Google Scholar]
  21. Xiphias gladius (errata version published in 2016). The IUCN Red List of Threatened Species 2011: e.T23148A88828055. Available online: https://0-dx-doi-org.brum.beds.ac.uk/10.2305/IUCN.UK.2011-2.RLTS.T23148A9422329.en. (accessed on 17 December 2017).
  22. Pinsky, M.L.; Palumbi, S.R. Meta-analysis reveals lower genetic diversity in overfished populations. Mol. Ecol. 2014, 23, 29–39. [Google Scholar] [CrossRef]
  23. Frankham, R. Genetics and extinction. Biol. Conserv. 2005, 126, 131–140. [Google Scholar] [CrossRef]
  24. Spielman, D.; Brook, B.W.; Frankham, R. Most species are not driven to extinction before genetic factors impact them. Proc. Natl. Acad. Sci. USA 2004, 101, 15261–15264. [Google Scholar] [CrossRef] [Green Version]
  25. Hauser, L.; Adcock, G.J.; Smith, P.J.; Ramírez, J.H.B.; Carvalho, G.R. Loss of microsatellite diversity and low effective population size in an overexploited population of New Zealand snapper (Pagrus auratus). Proc. Natl. Acad. Sci. USA 2002, 99, 11742–11747. [Google Scholar] [CrossRef] [Green Version]
  26. Turner, T.F.; Wares, J.P.; Gold, J.R. Genetic effective size is three orders of magnitude smaller than adult census size in an abundant, estuarine-dependent marine fish (Sciaenops ocellatus). Genetics 2002, 162, 1329–1339. [Google Scholar]
  27. Hutchinson, W.F.; Oosterhout, C.; Rogers, S.I.; Carvalho, G.R. Temporal analysis of archived samples indicates marked genetic changes in declining North Sea cod (Gadus morhua). Proc. R. Soc. Lond. Ser. B Biol. Sci. 2003, 270, 2125–2132. [Google Scholar] [CrossRef] [Green Version]
  28. Allendorf, F.W.; England, P.R.; Luikart, G.; Ritchie, P.A.; Ryman, N. Genetic effects of harvest on wild animal populations. Trends Ecol. Evol. 2008, 23, 327–337. [Google Scholar] [CrossRef] [PubMed]
  29. Thompson, J.D.; Gibson, T.J.; Higgins, D.G. Multiple sequence alignment using ClustalW and ClustalX. Curr. Protoc. Bioinform. 2003, 1, 2–3. [Google Scholar] [CrossRef] [PubMed]
  30. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  31. Excoffier, L.; Lischer, H.E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 565–567. [Google Scholar] [CrossRef] [PubMed]
  32. Heck, K.L., Jr.; van Belle, G.; Simberloff, D. Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size. Ecology 1975, 56, 1459–1461. [Google Scholar] [CrossRef]
  33. Hsieh, T.C.; Ma, K.H.; Chao, A. iNEXT: An R package for rarefaction and extrapolation of species diversity (H ill numbers). Methods Ecol. Evol. 2016, 7, 1451–1456. [Google Scholar] [CrossRef]
  34. Chao, A.; Gotelli, N.J.; Hsieh, T.C.; Sander, E.L.; Ma, K.H.; Colwell, R.K.; Ellison, A.M. Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies. Ecol. Monogr. 2014, 84, 45–67. [Google Scholar] [CrossRef] [Green Version]
  35. Bandelt, H.-J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef] [PubMed]
  36. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 1992, 131, 479–491. [Google Scholar] [PubMed]
  37. Kuhner, M.K. LAMARC 2.0: Maximum likelihood and Bayesian estimation of population parameters. Bioinformatics 2006, 22, 768–770. [Google Scholar] [CrossRef] [Green Version]
  38. Hayashi, K.; Yandell, D.W. How sensitive is PCR-SSCP? Hum. Mutat. 1993, 2, 338–346. [Google Scholar] [CrossRef]
  39. Garibaldi, F. Effects of the introduction of the mesopelagic longline on catches and size structure of swordfish in the Ligurian sea (western Mediterranean). Collect. Vol. Sci. Pap. ICCAT 2015, 71, 2006–2014. [Google Scholar]
  40. Garibaldi, F.; Lanteri, L. Notes about a tagged/recaptured swordfish in the Ligurian Sea (western Mediterranean). Collect. Vol. Sci. Pap. ICCAT 2017, 74, 1354–1361. [Google Scholar]
  41. Bertolino, F.; Camolese, C.; Dell’Aquila, M.; Mariani, A.; Valastro, M. Swordfish (Xiphias gladius L.) fisheries using mesopelagic longline in the Mediterranean Sea by Italian fishing fleet. Collect. Vol. Sci. Pap. ICCAT 2015, 71, 2073–2078. [Google Scholar]
  42. Abid, N.; Laglaoui, A.; Arakrak, A.; Bakkali, M. The reproductive biology of swordfish (Xiphias gladius) in the Strait of Gibraltar. J. Mar. Biol. Assoc. U. K. 2019, 99, 649–659. [Google Scholar] [CrossRef]
  43. Dulvy, N.K.; Sadovy, Y.; Reynolds, J.D. Extinction vulnerability in marine populations. Fish Fish. 2003, 4, 25–64. [Google Scholar] [CrossRef]
  44. Dulvy, N.K.; Jennings, S.; Rogers, S.I.; Maxwell, D.L. Threat and decline in fishes: An indicator of marine biodiversity. Can. J. Fish. Aquat. Sci. 2006, 63, 1267–1275. [Google Scholar] [CrossRef]
  45. Kenchington, E.L. 14 The Effects of Fishing on Species and Genetic Diversity. Responsible Fish. Mar. Ecosyst. 2003, 235. [Google Scholar]
  46. Smith, P.J.; Francis, R.; McVeagh, M. Loss of genetic diversity due to fishing pressure. Fish. Res. 1991, 10, 309–316. [Google Scholar] [CrossRef]
  47. Ruggeri, P.; Splendiani, A.; Di Muri, C.; Fioravanti, T.; Santojanni, A.; Leonori, I.; De Felice, A.; Biagiotti, I.; Carpi, P.; Arneri, E. Coupling demographic and genetic variability from archived collections of European anchovy (Engraulis encrasicolus). PLoS ONE 2016, 11, e0151507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. FAO. The state of world fisheries and aquaculture 2016. Contributing to Food Security and Nutrition for all; FAO: Rome, Italy, 2016. [Google Scholar]
Diversity 12 00170 g001 550
Figure 1. Map of the Mediterranean Sea with the Swordfish sampling sites. Blue triangles represent historical sampling localities from Viñas et al. [17] and orange dots represent sampling locations of current Swordfish samples.
Figure 1. Map of the Mediterranean Sea with the Swordfish sampling sites. Blue triangles represent historical sampling localities from Viñas et al. [17] and orange dots represent sampling locations of current Swordfish samples.
Diversity 12 00170 g001
Diversity 12 00170 g002 550
Figure 2. MJ network for the complete dataset of Mediterranean swordfish. White slices indicate haplotype observed by Viñas [17] and black slices recent sample (this work). Small red circles correspond to missing (or hypothetical) haplotypes.
Figure 2. MJ network for the complete dataset of Mediterranean swordfish. White slices indicate haplotype observed by Viñas [17] and black slices recent sample (this work). Small red circles correspond to missing (or hypothetical) haplotypes.
Diversity 12 00170 g002
Diversity 12 00170 g003 550
Figure 3. Size-based rarefaction (solid curves) and extrapolation (dashed curves) with 95% confidence intervals (shaded areas) comparing ancient Mediterranean swordfish (blue) and contemporary (orange) samples.
Figure 3. Size-based rarefaction (solid curves) and extrapolation (dashed curves) with 95% confidence intervals (shaded areas) comparing ancient Mediterranean swordfish (blue) and contemporary (orange) samples.
Diversity 12 00170 g003
Table
Table 1. Sampling details for swordfish analyzed in this study.
Table 1. Sampling details for swordfish analyzed in this study.
Sampling Area (FAO Fishing Area/Geographical Subarea)Sampling DatenMean LJFL (cm) (s.d.)
Balearic Sea (GSA 5, 6)July 2016–August 2016–September 201684122.1 (22,2)
September 2019
Southern Sicily (GSA 15, 16)July 201659133.7 (32.0)
June 2017–July 2017
June 2018
Aegean Sea (GSA 22)August 201617124.2 (21.0)
South Adriatic Sea (GSA 18)September 201662113.5 (16.6)
Tyrrhenian Sea (GSA 10)May 201716120.8 (23.8)
Sardinian Sea (GSA 11.2)June 2017–October 2017 49125.0 (32.6)
May 2018–Jun 2018
Total 287
Table
Table 2. Variable sites in the control-region segment of the 36 haplotypes of swordfish mtDNA detected in the modern sample. Numbers above sites refer to nucleotide position on the L-strand of the control region of the swordfish. Identity with sequence 1 is indicated by dots and deletions of 1 residue with dashes.
Table 2. Variable sites in the control-region segment of the 36 haplotypes of swordfish mtDNA detected in the modern sample. Numbers above sites refer to nucleotide position on the L-strand of the control region of the swordfish. Identity with sequence 1 is indicated by dots and deletions of 1 residue with dashes.
Genbank Accession Number2912171819202122232446606777839192101110111118120122123125126127132134135138158159160161165166172173178180186187190194217257262
AY650768GTATACA----CAACATTAGTTAACTGTTATCATTTAAATCC-TGGTAG
MN652597.......TACA..........C.........T..........-......
AY650778.......----..........C.........T..........-......
AY650860.......----..........C.........TG.........-......
MN652598A......----...............................-......
MN652599.......----.........................G.....-......
EU827759.......----....G..........................-......
MN652600.......TACA.....................G.........-......
MN652601.......----..............C......G.........-......
AY650836.......----.....................G.........-......
AY650821.......----.....................C.........-......
MN652602.......----......................C........-......
AY650858.......----..G............................-......
MN652603.......----..........................G....-......
MN652604.......----.........C...T.................-......
MN652605.......----.............T.................-......
MN652606.......----..................G............-......
MN652607.......----...T...........................-......
MN652608.-.....----.................C..T........GT-.A....
AY650781.......----.................C..T........GT-.A....
MN652609.......----...............A.C..T........GT-.A....
AY650763..-....----...T.........T.A.C..T.......C..-.A...A
MN652595..-....----.G.T.........T.A.C..T.......C..-.A...A
AY650855..-....TACA...T.........T.A.C..T.......C..-.A.C.A
MN652596..-....----...T.........T.A.C..T.......C..-.....A
AY650861..-....----...T.........T.A.C..T......GC..-......
AY650809...--------T..T...TC..GG..ACCG.TG.C......TC....GA
MN652610...--------T..T...TC..GG..ACCG.TG.C......TCC...GA
AY650762...--------T..T...TC..GG..ACCG.TG.C......TC..A.GA
AY650814...--------T..T...TC..GG..ACCG.TG.CC.....TC..A.GA
AY650805...--------T..T...TC..GG..ACCG.TG.C..G...TC..A.GA
AY650761...--------T..T..CTC..GG..ACCG.TG.C......TC..A.GA
AY650829...--------T..T...TT..GG..ACCG.TG.C......TC..A.GA
EU827771...--------T..T...TC..GG..A.CGCTG.C......TC..A.GA
MN652611...--------T..T.C.TC..GG..ACCG.TG.C......TC..A..A
EU827787...--------T..T...TT..GG..ACCG.TG.C......TC..A.G.
Table
Table 3. Summary of molecular diversity indices for Mediterranean swordfish temporal samples for pooled clades and each clade separately. The number of individuals (N) and haplotypes (H), the haplotype diversity (h), the nucleotide diversity (π) and effective population size parameter (θ) are shown.
Table 3. Summary of molecular diversity indices for Mediterranean swordfish temporal samples for pooled clades and each clade separately. The number of individuals (N) and haplotypes (H), the haplotype diversity (h), the nucleotide diversity (π) and effective population size parameter (θ) are shown.
TitleNHh (s.d.)π (s.d.)θ (95% C.I.)
Viñas et al., 2010
Clade I + II251930.946 (0.007)0.160 (0.080)0.142 (0.121–0.188)
Clade I156600.910 (0.014)0.034 (0.020)0.076 (0.060–0.107)
Clade II95330.861 (0.031)0.006 (0.004)0.044 (0.0.29–0.062)
This study
Clade I + II286360.844 (0.010)0.130 (0.070)0.033 (0.028–0.044)
Clade I190260.793 (0.018)0.029 (0.015)0.019 (0.015–0.030)
Clade II96100.418 (0.063)0.003 (0.002)0.010 (0.006–0.017)
Table
Table 4. AMOVA results of Mediterranean swordfish temporal sample comparison.
Table 4. AMOVA results of Mediterranean swordfish temporal sample comparison.
Source of VariationVariance ComponentPercentage of VariationFixation Index (ФST)p-Value
Whole dataset
Among temporal samples0.0081.770.01770.0000
Within samples0.44698.23
Clade I
Among temporal samples0.0081.850.01840.0005
Within samples0.42398.15
Clade II
Among temporal samples0.04412.100.1210.0000
Within samples0.31987.90

Share and Cite

MDPI and ACS Style

Righi, T.; Splendiani, A.; Fioravanti, T.; Casoni, E.; Gioacchini, G.; Carnevali, O.; Caputo Barucchi, V. Loss of Mitochondrial Genetic Diversity in Overexploited Mediterranean Swordfish (Xiphias gladius, 1759) Population. Diversity 2020, 12, 170. https://0-doi-org.brum.beds.ac.uk/10.3390/d12050170

AMA Style

Righi T, Splendiani A, Fioravanti T, Casoni E, Gioacchini G, Carnevali O, Caputo Barucchi V. Loss of Mitochondrial Genetic Diversity in Overexploited Mediterranean Swordfish (Xiphias gladius, 1759) Population. Diversity. 2020; 12(5):170. https://0-doi-org.brum.beds.ac.uk/10.3390/d12050170

Chicago/Turabian Style

Righi, Tommaso, Andrea Splendiani, Tatiana Fioravanti, Elia Casoni, Giorgia Gioacchini, Oliana Carnevali, and Vincenzo Caputo Barucchi. 2020. "Loss of Mitochondrial Genetic Diversity in Overexploited Mediterranean Swordfish (Xiphias gladius, 1759) Population" Diversity 12, no. 5: 170. https://0-doi-org.brum.beds.ac.uk/10.3390/d12050170

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop