The population genomics of yellowfin tuna (Thunnus albacares) at global geographic scale challenges current stock delineation

Integration of genomic and ecological methods inform management of an undescribed, yet highly exploited, sardine species

Assessing genetic diversity within species is key for conservation strategies in the context of human-induced biotic changes. This is important in marine systems, where many species remain undescribed while being overfished, and conflicts between resource-users and conservation agencies are common. Combining niche modelling with population genomics can contribute to resolving those conflicts by identifying management units and understanding how past climatic cycles resulted in current patterns of genetic diversity. We addressed these issues on an undescribed but already overexploited species of sardine of the genus Harengula. We find that the species distribution is determined by salinity and depth, with a continuous distribution along the Brazilian mainland and two disconnected oceanic archipelagos. Genomic data indicate that such biogeographic barriers are associated with two divergent intraspecific lineages. Changes in habitat availability during the last glacial cycle led to different demographic histories among stocks. One coastal population experienced a 3.6-fold expansion, whereas an island-associated population contracted 3-fold, relative to the size of the ancestral population. Our results indicate that the island population should be managed separately from the coastal population, and that a Marine Protected Area covering part of the island population distribution can support the viability of this lineage.

Development of biomarkers to distinguish different origins of red seabreams (Pagrus major) from Korea and Japan by fatty acid, amino acid, and mineral profiling

Red seabream (Pagrus major) has been one of the most popular fish in East Asia since early times. However, the discharge of nuclear wastewater into the sea following the Fukushima nuclear disaster in Japan has led to violations of the country of origin labeling. Therefore, the aim of the present study was to determine the origin of fish based on fatty acid, amino acid, and mineral analyses, and to develop biomarkers that can discriminate between Japanese and Korean red seabream. To identify the differences between the two groups, 29 fatty acid families, 17 amino acids, and 4 minerals were analyzed in 60 fish samples (standard sample collected in autumn), and fatty acid profiles were analyzed using heatmap with hierarchical clustering analysis and orthogonal projections to latent structures discriminant analysis. The top 10 fatty acids that were different between the two groups were selected from all seasonal fish samples by combining variable importance in projection scores and p-values. According to the receiver operating characteristic curve analysis results, we proposed percentage linoleic acid (C18:2n-6, cis) as a candidate biomarker with excellent sensitivity and specificity. This study introduces a strategy to identify the origins of red seabream using linoleic acid obtained from fatty acid analysis.

What Darwin could not see: island formation and historical sea levels shape genetic divergence and island biogeography in a coastal marine species

Oceanic islands play a central role in the study of evolution and island biogeography. The Galapagos Islands are one of the most studied oceanic archipelagos but research has almost exclusively focused on terrestrial organisms compared to marine species. Here we used the Galapagos bullhead shark (Heterodontus quoyi) and single nucleotide polymorphisms (SNPs) to examine evolutionary processes and their consequences for genetic divergence and island biogeography in a shallow-water marine species without larval dispersal. The sequential separation of individual islands from a central island cluster gradually established different ocean depths between islands that pose barriers to dispersal in H. quoyi. Isolation by resistance analysis suggested that ocean bathymetry and historical sea level fluctuations modified genetic connectivity. These processes resulted in at least three genetic clusters that exhibit low genetic diversity and effective population sizes that scale with island size and the level of geographic isolation. Our results exemplify that island formation and climatic cycles shape genetic divergence and biogeography of coastal marine organisms with limited dispersal comparable to terrestrial taxa. Because similar scenarios exist in oceanic islands around the globe our research provides a new perspective on marine evolution and biogeography with implications for the conservation of island biodiversity.

Phylogeography of sharks and rays: a global review based on life history traits and biogeographic partitions

Considerable research exists on the life history traits, evolutionary history, and environmental factors that shape the population genetic structure of marine organisms, including sharks and rays. Conservation concerns are particularly strong for this group as they are highly susceptible to anthropogenic stressors due to a combination of life history traits including late maturity and low fecundity. Here, we provide a review and synthesis of the global phylogeography of sharks and rays. We examined existing data for 40 species of sharks belonging to 17 genera and 19 species of rays belonging to 11 genera. Median joining haplotype networks were constructed for each species for the mtDNA cytochrome C oxidase subunit I (COI), and an Analysis of Molecular Variance (AMOVA) was conducted to understand patterns of genetic diversity and structure across the three major ocean basins-the Indian, Atlantic and Pacific Oceans. Haplotype networks showed very shallow coalescence in most species, a finding previously reported for marine teleosts. Star topologies were predominant among sharks while complex mutational topologies predominated among rays, a finding we attribute to extremely limited dispersal in the early life history of rays. Population structuring varied amongst species groups, apparently due to differences in life history traits including reproductive philopatry, site fidelity, pelagic habitat, migratory habits, and dispersal ability. In comparison to reef-associated and demersal species, pelagic and semi pelagic species showed lower levels of structure between and within ocean basins. As expected, there is variation between taxa and groups, but there are also some broad patterns that can guide management and conservation strategies.

Genomic assessment reveals signal of adaptive selection in populations of the Spotted rose snapper Lutjanus guttatus from the Tropical Eastern Pacific

Background

The lack of barriers in the marine environment has promoted the idea of panmixia in marine organisms. However, oceanographic conditions and habitat characteristics have recently been linked to genetic structure in marine species. The Tropical Eastern Pacific (TEP) is characterized by dynamic current systems and heterogeneous oceanographic conditions. The Gulf of Panama (part of the equatorial segment for the TEP) is influenced by a complex current system and heterogeneous environment, which has been shown to limit the gene flow for shoreline species. Next Generation Sequencing (NGS) has contributed to detect genetic differences in previously reported panmictic species by the assessment of loci associated with selection and to understand how selection acts affects marine populations. Lutjanus guttatus is a species distributed in the TEP for which previous studies using mitochondrial data recovered a panmictic pattern along its distributional range. In this study, we used SNP data of L. guttatus individuals sampled along its range to evaluate population genetic structure and investigate whether oceanographic factors influence the species’ genetic architecture. Finally, we assessed the role of adaptive selection by evaluating the contribution of outlier and neutral loci to genetic divergence.

Methods

The RADcap method was used to obtain 24 million paired reads for 123 individuals of L. guttatus covering nearly all its distributional area. Genetic variation was assessed using both spatial and non-spatial methods by comparing three different data sets: (i) a Combined Loci (CL dataset = 2003 SNPs); a search for putative loci under selection allowed the evaluation of (ii) Neutral Loci (NL dataset = 1858 SNPs) and (iii) Outlier Loci (OL dataset = 145 SNPs). We used the estimating effective migration surface (EEMS) approach to detect possible barriers to gene flow.

Results

Genetic differences were found in the OL dataset, showing two clusters (Northern and Southern), whereas NL showed no differences. This result may be related to the Selection-Migration balance model. The limit between the Northern and Southern groups was in the Gulf of Panama, which has been previously identified as a barrier to gene flow for other species, mainly due to its heterogeneous oceanographic conditions. The results suggest that selection plays an important role in generating genetic differences in Lutjanus guttatus. A migration corridor was detected that coincides with the Costa Rica Coastal Current that flows from Central America to the Gulf of California, allowing the homogenization of the northern population. In the Southern cluster, a migration corridor was observed with the OL from Panama to Colombia, which could be associated with the currents found in the Gulf of Panama. Genetic variation found in the OL of Lutjanus guttatus highlights the usefulness of NGS data in evaluating the role of selection in population differentiation.

Stepping up to genome scan allows stock differentiation in the worldwide distributed blue shark Prionace glauca

The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome-wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo-west Pacific samples (Southern population). Significant pairwise F_ST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies.

Lack of genetic differentiation in yellowfin tuna has conservation implications in the Eastern Pacific Ocean

Yellowfin tuna, Thunnus albacares, is an important global fishery and of particular importance in the Eastern Pacific Ocean (EPO). According to the 2019 Inter-American Tropical Tuna Commission (IATTC) assessment, yellowfin tuna within the EPO is a single stock, and is being managed as one stock. However, previous studies indicate site fidelity, or limited home ranges, of yellowfin tuna which suggests the potential for multiple yellowfin tuna stocks within the EPO, which was supported by a population genetic study using microsatellites. If numerous stocks are present, management at the wrong spatial scales could cause the loss of minor yellowfin tuna populations in the EPO. In this study we used double digestion RADseq to assess the genetic structure of yellowfin tuna in the EPO. A total of 164 yellowfin tuna from Cabo San Lucas, México, and the Galápagos Islands and Santa Elena, Ecuador, were analysed using 18,011 single nucleotide polymorphisms. Limited genetic differentiation (F_ST = 0.00058–0.00328) observed among the sampling locations (México, Ecuador, Peru, and within Ecuador) is consistent with presence of a single yellowfin tuna population within the EPO. Our findings are consistent with the IATTC assessment and provide further evidence of the need for transboundary cooperation for the successful management of this important fishery throughout the EPO.

Fisheries genomics of snapper ( Chrysophrys auratus ) along the west Australian coast

The efficacy of fisheries management strategies depends on stock assessment and management actions being carried out at appropriate spatial scales. This requires understanding of spatial and temporal population structure and connectivity, which is challenging in weakly structured and highly connected marine populations. We carried out a population genomics study of the heavily exploited snapper (Chrysophrys auratus) along ~2600 km of the Australian coastline, with a focus on Western Australia (WA). We used 10,903 filtered SNPs in 341 individuals from eight sampling locations to characterize population structure and connectivity in snapper across WA and to assess if current spatial scales of stock assessment and management agree with evidence from population genomics. Our dataset also enabled us to investigate temporal stability in population structure as well as connectivity between WA and its nearest, eastern jurisdictional neighbour. As expected for a species influenced by the extensive ocean boundary current in the region, low genetic differentiation and high connectivity were uncovered across WA. However, we did detect strong isolation by distance and genetic discontinuities in the mid‐west and south‐east. The discontinuities correlate with boundaries between biogeographic regions, influenced by on‐shelf oceanography, and the sites of important spawning aggregations. We also detected temporal instability in genetic structure at one of our sites, possibly due to interannual variability in recruitment in adjacent regions. Our results partly contrast with the current spatial management of snapper in WA, indicating the likely benefits of a review. This study supports the value of population genomic surveys in informing the management of weakly structured and wide‐ranging marine fishery resources.

Population structure of blackfin tuna (Thunnus atlanticus) in the western Atlantic Ocean inferred from microsatellite loci

The blackfin tuna, Thunnus atlanticus, is a small tropical tuna exploited by recreational and commercial fisheries in various parts of its range. Information on stock structure is needed to develop management plans for this species but is currently lacking. In this work, 470 blackfin tuna from nine geographic populations were assayed at 13 homologous microsatellite markers to provide a first assessment of stock structure across the species range. The overall divergence among locality samples was very low (overall FST = 0.0004) indicating high connectivity of blackfin tuna across their range. No clear grouping of localities in differentiated units was inferred but structuring followed a weak isolation by distance pattern (r = 0.16, P = 0.032). Pairwise exact tests and spatial analysis of molecular variance suggested divergence of the sample collected offshore Baía Formosa (Brazil) possibly reflecting reproductive isolation of Brazilian populations from those in the Caribbean region and further north. Further study of the status of Brazilian populations and the transition between this region and the Caribbean is warranted. Cryptic subdivision within the Northern Hemisphere part of the range is possible and should be evaluated using increased marker density and a more comprehensive geographic coverage.

Genomic Stock Structure of the Marine Teleost Tarakihi (Nemadactylus macropterus) Provides Evidence of Potential Fine-Scale Adaptation and a Temperature-Associated Cline Amid Panmixia

Tarakihi (Nemadactylus macropterus) is an important fishery species with widespread distribution around New Zealand and off the southern coasts of Australia. However, little is known about whether the populations are locally adapted or genetically structured. To address this, we conducted whole-genome resequencing of 175 tarakihi from around New Zealand and Tasmania (Australia) to obtain a dataset of 7.5 million genome-wide and high-quality single nucleotide polymorphisms (SNPs). Variant filtering, FST-outlier analysis, and redundancy analysis (RDA) were used to evaluate population structure, adaptive structure, and locus-environment associations. A weak but significant level of neutral genetic differentiation was found between tarakihi from New Zealand and Tasmania (FST = 0.0054–0.0073, P ≤ 0.05), supporting the existence of at least two separate reproductive stocks. No clustering was detected among the New Zealand populations (ΦST < 0.001, P = 0.77). Outlier-based, presumably adaptive variation suggests fine-scale adaptive structure between locations around central New Zealand off the east (Wairarapa, Cape Campbell, and Hawke’s Bay) and the west coast (Tasman Bay/Golden Bay and Upper West Coast of South Island). Allele frequencies from 55 loci were associated with at least one of six environmental variables, of which 47 correlated strongly with yearly mean water temperature. Although genes associated with these loci are linked to various functions, the most common functions were integral components of membrane and cilium assembly. Projection of the RDA indicates the existence of a latitudinal temperature cline. Our work provides the first genomic insights supporting panmixia of tarakihi in New Zealand and evidence of a genomic cline that appears to be driven by the temperature gradients, together providing crucial information to inform the stock assessment of this species, and to widen the insights of the ecological drivers of adaptive variation in a marine species.

Are genomic updates of well-studied species worth the investment for conservation? A case study of the Critically Endangered Magdalena River turtle

As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used RADseq to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment.

Transporter-interfering chemicals inhibit P-glycoprotein of yellowfin tuna (Thunnus albacares)

Marine pollutants bioaccumulate at high trophic levels of marine food webs and are transferred to humans through consumption of apex species. Yellowfin tuna (Thunnus albacares) are marine predators, and one of largest commercial fisheries in the world. Previous studies have shown that yellowfin tuna can accumulate high levels of persistent organic pollutants, including Transporter Interfering Chemicals (TICs), which are chemicals shown to bind to mammalian xenobiotic transporters and interfere with their function. Here, we examined the extent to which these same compounds might interfere with the activity of the yellowfin tuna (Thunnus albacares) ortholog of this transporter. To accomplish this goal we identified, expressed, and functionally assayed tuna ABCB1. The results demonstrated a common mode of vertebrate ABCB1 interaction with TICs that predicts effects across these species, based on high conservation of specific interacting residues. Importantly several TICs showed potent inhibition of Ta-ABCB1, such as the organochlorine pesticides Endrin (EC₅₀ = 1.2 ± 0.2 μM) and Mirex (EC₅₀ = 2.3 ± 0.9 μM). However, unlike the effects observed on mouse ABCB1, low concentrations of the organochlorine pesticide TICs p,p'-DDT and its metabolite p,p'-DDD co-stimulated verapamil-induced Ta-ABCB1 ATPase activity possibly suggesting a low transport activity for these ligands in tuna. These results provide a mechanistic basis for understanding the potential vulnerability of tuna to these ubquitous pollutants.

Genomic assessment of global population structure in a highly migratory and habitat versatile apex predator, the tiger shark (Galeocerdo cuvier)

Understanding the population dynamics of highly mobile, widely distributed, oceanic sharks, many of which are overexploited, is necessary to aid their conservation management. We investigated the global population genomics of tiger sharks (Galeocerdo cuvier), a circumglobally distributed, apex predator displaying remarkable behavioral versatility in its diet, habitat use (near coastal, coral reef, pelagic), and individual movement patterns (spatially resident to long-distance migrations). We genotyped 242 tiger sharks from 10 globally distributed locations at more than 2000 single nucleotide polymorphisms. Although this species often conducts massive distance migrations, the data show strong genetic differentiation at both neutral (FST=0.125-0.144) and candidate outlier loci (FST=0.570-0.761) between western Atlantic and Indo-Pacific sharks, suggesting the potential for adaptation to the environments specific to these oceanic regions. Within these regions, there was mixed support for population differentiation between northern and southern hemispheres in the western Atlantic, and none for structure within the Indian Ocean. Notably, the results demonstrate a low level of population differentiation of tiger sharks from the remote Hawaiian archipelago compared to sharks from the Indian Ocean (FST=0.003-0.005, P<0.01). Given concerns about biodiversity loss and marine ecosystem impacts caused by overfishing of oceanic sharks in the midst of rapid environmental change, our results suggest it imperative that international fishery management prioritize conservation of the evolutionary potential of the highly genetically differentiated Atlantic and Indo-Pacific populations of this unique apex predator. Furthermore, we suggest targeted management attention to tiger sharks in the Hawaiian archipelago based on a precautionary biodiversity conservation perspective.

Demography of swordfish (Xiphias gladius Linneus) populations from the coasts of Turkey, based on mitochondrial DNA and microsatellites

The genetic diversity of the Mediterranean swordfish (Xiphias gladius Linneus) has not been explored extensively at its easternmost range so far. In this study, modern X. gladius samples from the eastern part of the Mediterranean basin, north of the Aegean Sea (Aegean-2013, n = 26) and the Mediterranean coast of Turkey (N.Levantine-2013, n = 42) were studied genetically, along with ancient samples from Yenikapı excavation (n = 6). Partial mitochondrial DNA control region sequences (entire sequences, clade I and clade II) were evaluated spatially and temporally together with previously published sequences (Alvarado Bremer et al., Molecular Phylogenetics and Evolution, 2005, 36, 169-187; Viñas et al., ICES Journal of Marine Science, 2010, 67, 1222-1229; Righi et al., Diversity, 2020, 12, 170) from the rest of the Mediterranean Sea. Pair-wise F_ST and pair-wise AMOVA tests showed that, in general, groups of eastern populations and western Mediterranean populations have not genetically differed from each other significantly nearly in the past 20 years. Therefore, the results direct reconsideration of previous descriptions of population sub-structure within the Mediterranean and support high gene flow throughout the region. On the contrary, the results of this study confirmed the existence of genetic diversity differences between western and eastern Mediterranean, with eastern being low. One-tailed permutation tests revealed that θ, which is directly proportional to long-term female effective population size (Ne), decreased significantly (P < 0.05) in both regions over the past two decades. On the Turkish coasts, θ is not significantly different from that of the nearly contemporary eastern Mediterranean population. Nonetheless, θ of the ancient sample was consistently and significantly (P < 0.001) higher than those of the eastern and western Mediterranean populations in clade I and clade II. Furthermore, it contains two mitochondrial haplotypes that are not observed in modern samples, suggesting that the Ne of X. gladius in the eastern was high in Byzantium times. Eight microsatellite loci were also genotyped in modern samples. The microsatellite-based present Ne estimate of the pooled Aegean-2013 and N.Levantine-2013 populations was lower than 1000 according to the upper limit of 95% c.i. and possibly even lower than 100 according to the mean of posterior distribution of the present Ne estimate calculated by the software package MSVAR. These alarming genetic signals for the sustainability of X. gladius on the coasts of Turkey are in agreement with the nearly collapsing X. gladius fisheries as depicted also in the fisheries statistics. Overall, congruent with the previous studies, the data presented here show that sustainability of the X. gladius population in Mediterranean is under major threat. Therefore, X. gladius around the Turkish coasts need an immediate stringent action and management plan.

Adaptive markers distinguish North and South Pacific Albacore amid low population differentiation

Albacore (Thunnus alalunga) support an economically valuable global fishery, but surprisingly little is known about the population structure of this highly migratory species. Physical tagging data suggest that Albacore from the North and South Pacific Ocean are separate stocks, but results from previous genetic studies did not support this two stock hypothesis. In addition, observed biological differences among juveniles suggest that there may be population substructure in the North Pacific. We used double-digest restriction site-associated DNA sequencing to assess population structure among 308 Albacore caught in 12 sample areas across the Pacific Ocean (10 North, 2 South). Since Albacore are highly migratory and spawning areas are unknown, sample groups were not assumed to be equivalent to populations and the genetic data were analyzed iteratively. We tested for putatively adaptive differences among groups and for genetic variation associated with sex. Results indicated that Albacore in the North and South Pacific can be distinguished using 84 putatively adaptive loci, but not using the remaining 12,788 presumed neutral sites. However, two individuals likely represent F1 hybrids between the North and South Pacific populations, and 43 Albacore potentially exhibit lower degrees of mixed ancestry. In addition, four or five cross-hemisphere migrants were potentially identified. No genetic evidence was found for population substructure within the North Pacific, and no loci appeared to distinguish males from females. Potential functions for the putatively adaptive loci were identified, but an annotated Albacore genome is required for further exploration. Future research should try to locate spawning areas so that life history, demography, and genetic population structure can be linked and spatiotemporal patterns can be investigated.

Genetic population structure of the Blackspot seabream (Pagellus bogaraveo): contribution of mtDNA control region to fisheries management

Marine fisheries management models have traditionally considered biological parameters and geopolitical boundaries. The result is the existence of fisheries management units that do not match genetic populations. However, this panorama is changing with the contribution of genetic and genomic data. Pagellus bogaraveo is a commercially important sparid in the northeast Atlantic, with three stock components being considered by ICES: the Celtic Sea and Bay of Biscay, Atlantic Iberian waters and the Azores. The northern stock collapsed (1975-1985) and is essential to characterize the genetic makeup of the species, particularly in the Iberian Peninsula, where it is managed as a single stock. The mitochondrial control region was used to screen the intraspecific diversity and population structure of individuals from six locations across the species range. The genetic diversity found is similar among sites, and there is differentiation between the Azores and the remaining locations.

A review of the fisheries, life history and stock structure of tropical tuna (skipjack Katsuwonus pelamis, yellowfin Thunnus albacares and bigeye Thunnus obesus) in the Indian Ocean

Skipjack (Katsuwonus pelamis), yellowfin (Thunnus albacares) and bigeye (Thunnus obesus) tuna are the target species of tropical tuna fisheries in the Indian Ocean, with high commercial value in the international market. High fishing pressure over the past three decades has raised concerns about their sustainability. Understanding life history strategies and stock structure is essential to determine species resilience and how they might respond to exploitation. Here we provide a comprehensive review of available knowledge on the biology, ecology, and stock structure of tropical tuna species in the Indian Ocean. We describe the characteristics of Indian Ocean tropical tuna fisheries and synthesize skipjack, yellowfin, and bigeye tuna key life history attributes such as biogeography, trophic ecology, growth, and reproductive biology. In addition, we evaluate the available literature about their stock structure using different approaches such as analysis of fisheries data, genetic markers, otolith microchemistry and tagging, among others. Based on this review, we conclude that there is a clear lack of ocean basin-scale studies on skipjack, yellowfin and bigeye tuna life history, and that regional stock structure studies indicate that the panmictic population assumption of these stocks should be investigated further. Finally, we identify specific knowledge gaps that should be addressed with priority to ensure a sustainable and effective management of these species.

Connectivity and population structure of albacore tuna across southeast Atlantic and southwest Indian Oceans inferred from multidisciplinary methodology

Albacore tuna (Thunnus alalunga) is an important target of tuna fisheries in the Atlantic and Indian Oceans. The commercial catch of albacore is the highest globally among all temperate tuna species, contributing around 6% in weight to global tuna catches over the last decade. The accurate assessment and management of this heavily exploited resource requires a robust understanding of the species’ biology and of the pattern of connectivity among oceanic regions, yet Indian Ocean albacore population dynamics remain poorly understood and its level of connectivity with the Atlantic Ocean population is uncertain. We analysed morphometrics and genetics of albacore (n = 1,874) in the southwest Indian (SWIO) and southeast Atlantic (SEAO) Oceans to investigate the connectivity and population structure. Furthermore, we examined the species’ dispersal potential by modelling particle drift through major oceanographic features. Males appear larger than females, except in South African waters, yet the length–weight relationship only showed significant male–female difference in one region (east of Madagascar and Reunion waters). The present study produced a genetic differentiation between the southeast Atlantic and southwest Indian Oceans, supporting their demographic independence. The particle drift models suggested dispersal potential of early life stages from SWIO to SEAO and adult or sub-adult migration from SEAO to SWIO.

Population genomics reveals a mismatch between management and biological units in green abalone (Haliotis fulgens)

Effective fishery management strategies should be based on stock delimitation and knowledge of the spatial scale at which species are distributed. However, a mismatch often occurs between biological and management units of fishery resources. The green abalone (Haliotis fulgens) supports an important artisanal fishery in the west coast of the Baja California Peninsula (BCP), Mexico, which has shown a declining tendency despite the several management measures. Thus, the aim of this study was to characterize the spatial patterns of neutral genomic variation of green abalone along the BCP to test whether the genomic structure patterns support the current green abalone management areas. To test this hypothesis, a set of 2,170 putative neutral single nucleotide polymorphisms discovered by a double digest restriction-site associated DNA approach was used on 10 locations along the BCP. The results revealed a population structure with three putative groups: Guadalupe Island and northern and southern BCP locations. The contemporary gene flow might be explained by local oceanographic features, where it is bidirectional within the southern region but with a predominant southward flow from the northern region. These findings indicated that the administrative areas did not match the biological units of H. fulgens fishery; hence, the stock assessment and management areas should be revised.

Mediterranean swordfish (Xiphias gladius Linnaeus, 1758) population structure revealed by microsatellite DNA: genetic diversity masked by population mixing in shared areas

Background

The Mediterranean swordfish stock is overfished and considered not correctly managed. Elucidating the patterns of the Mediterranean swordfish population structure constitutes an essential prerequisite for effective management of this fishery resource. To date, few studies have investigated intra-Mediterranean swordfish population structure, and their conclusions are controversial.

Methods

A panel of 20 microsatellites DNA was used to investigate fine-scale population structuring of swordfish from six main fishing areas of the Mediterranean Sea.

Results

This study provides evidence to reject the hypothesis of a single swordfish population within the Mediterranean Sea. DAPC analysis revealed the presence of three genetic clusters and a high level of admixture within the Mediterranean Sea. Genetic structure was supported by significant F_ST values while mixing was endorsed by the heterozygosity deficit observed in sampling localities indicative of a possible Wahlund effect, by sampling admixture individuals. Overall, our tests reject the hypothesis of a single swordfish population within the Mediterranean Sea. Homing towards the Mediterranean breeding areas may have generated a weak degree of genetic differentiation between populations even at the intra-basin scale.

Population Genomics of an Anadromous Hilsa Shad Tenualosa ilisha Species across Its Diverse Migratory Habitats: Discrimination by Fine-Scale Local Adaptation

The migration of anadromous fish in heterogenic environments unceasingly imposes a selective pressure that results in genetic variation for local adaptation. However, discrimination of anadromous fish populations by fine-scale local adaptation is challenging because of their high rate of gene flow, highly connected divergent population, and large population size. Recent advances in next-generation sequencing (NGS) have expanded the prospects of defining the weakly structured population of anadromous fish. Therefore, we used NGS-based restriction site-associated DNA (NextRAD) techniques on 300 individuals of an anadromous Hilsa shad (Tenualosa ilisha) species, collected from nine strategic habitats, across their diverse migratory habitats, which include sea, estuary, and different freshwater rivers. The NextRAD technique successfully identified 15,453 single nucleotide polymorphism (SNP) loci. Outlier tests using the F_ST OutFLANK and pcadapt approaches identified 74 and 449 SNPs (49 SNPs being common), respectively, as putative adaptive loci under a divergent selection process. Our results, based on the different cluster analyses of these putatively adaptive loci, suggested that local adaptation has divided the Hilsa shad population into two genetically structured clusters, in which marine and estuarine collection sites were dominated by individuals of one genetic cluster and different riverine collection sites were dominated by individuals of another genetic cluster. The phylogenetic analysis revealed that all the riverine populations of Hilsa shad were further subdivided into the north-western riverine (turbid freshwater) and the north-eastern riverine (clear freshwater) ecotypes. Among all of the putatively adaptive loci, only 36 loci were observed to be in the coding region, and the encoded genes might be associated with important biological functions related to the local adaptation of Hilsa shad. In summary, our study provides both neutral and adaptive contexts for the observed genetic divergence of Hilsa shad and, consequently, resolves the previous inconclusive findings on their population genetic structure across their diverse migratory habitats. Moreover, the study has clearly demonstrated that NextRAD sequencing is an innovative approach to explore how dispersal and local adaptation can shape genetic divergence of non-model anadromous fish that intersect diverse migratory habitats during their life-history stages.