
==== Front
PeerJ
PeerJ
PeerJ
PeerJ
2167-8359
PeerJ Inc. San Diego, USA

17928
10.7717/peerj.17928
Aquaculture, Fisheries and Fish Science
Genetics
Marine Biology
Population Biology
Genetic investigation of population structure in Atlantic chub mackerel, Scomber colias Gmelin, 1789 along the West African coast
http://orcid.org/0009-0003-0383-879X
Sbiba Salah eddine 12
Quintela María 3
Øyro Johanne 3
Dahle Geir 3
http://orcid.org/0000-0003-4562-045X
Jurado-Ruzafa Alba 4
Iita Kashona 5
Nikolioudakis Nikolaos 3
http://orcid.org/0000-0002-1788-0580
Bazairi Hocein 16
Chlaida Malika 2ma_chlaida@hotmail.com

1 Biodiversity, Ecology and Genome Laboratory, Faculty of Sciences, Mohammed V University in Rabat, Rabat, Morocco
2 Research and Development Unit on Marine Biology, National Institute of Fisheries Research, Casablanca, Morocco
3 Department of Population Genetics, Institute of Marine Research, Bergen, Norway
4 Oceanographic Centre of the Canary Islands, Spanish Institute of Oceanography (IEO-CSIC), Tenerife, Spain
5 National Marine Information and Research Centre (NATMIRC), Ministry of Fisheries and Marine Resources, Swakopmund, Namibia
6 University of Gibraltar, Europa Point Campus, Natural Sciences and Environment Research Hub, Gibraltar, Gibraltar
Jeffery Nicholas
5 9 2024
2024
12 e1792817 10 2023
24 7 2024
© 2024 FAO & IMR
2024
Sbiba et al.
https://creativecommons.org/licenses/by-nc/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits using, remixing, and building upon the work non-commercially, as long as it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.

Sustainable management of transboundary fish stocks hinges on accurate delineation of population structure. Genetic analysis offers a powerful tool to identify potential subpopulations within a seemingly homogenous stock, facilitating the development of effective, coordinated management strategies across international borders. Along the West African coast, the Atlantic chub mackerel (Scomber colias) is a commercially important and ecologically significant species, yet little is known about its genetic population structure and connectivity. Currently, the stock is managed as a single unit in West African waters despite new research suggesting morphological and adaptive differences. Here, eight microsatellite loci were genotyped on 1,169 individuals distributed across 33 sampling sites from Morocco (27.39°N) to Namibia (22.21°S). Bayesian clustering analysis depicts one homogeneous population across the studied area with null overall differentiation (FST = 0.0001ns), which suggests panmixia and aligns with the migratory potential of this species. This finding has significant implications for the effective conservation and management of S. colias within a wide scope of its distribution across West African waters from the South of Morocco to the North-Centre of Namibia and underscores the need for increased regional cooperation in fisheries management and conservation.

Scomber colias
Microsatellites
West Africa
Stock identification
Management
Food security
Food and Agriculture Organisation of the United NationsNorwegian Agency for Development CooperationInstitute of Marine Research in NorwayGlobal Environmental FacilityCanary Current Large Marine Ecosystem ProjectThis study was funded by EAF-Nansen project, a collaboration between the Food and Agriculture Organisation of the United Nations, the Norwegian Agency for Development Cooperation, and the Institute of Marine Research in Norway. This project also received support from Global Environmental Facility, as well as the Canary Current Large Marine Ecosystem project, implemented by the FAO, the United Nations Environment Program, and seven partner countries: Morocco, Mauritania, Senegal, The Gambia, Cape Verde, and Guinea Bissau. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
==== Body
pmcIntroduction

Sustainable stock management and conservation of marine resources require knowledge of the species’ population variability in a bio-geographical context. Failure to take into consideration the underlying components of the population structure of a species (e.g., spatio-temporal mixing of populations) can result in differential exploitation patterns and discrepancies between management measures among countries, which can lead to overexploitation of unique spawning components and subsequent resources declines (e.g., Allendorf et al., 2008; Kerr et al., 2017). As such, a reliable delineation of the biological units is needed to maintain the sustainability of fisheries, especially as climate change and other human pressures on marine ecosystems continue to intensify (Gissi et al., 2021; IPCC, 2022; Ramírez et al., 2022).

Genetic and genomic tools are commonly used to accurately define stock structure of fish (Ward, 2000) and to assess connectivity among marine populations (Bekkevold et al., 2015; Besnier et al., 2014; Dahle et al., 2018b; Hemmer-Hansen et al., 2019; Henriques et al., 2017; Le Moan, Bekkevold & Hemmer-Hansen, 2021). Using genetic and genomic tools to identify stock structure has proved useful for a variety of issues such as updating management plans (Mullins et al., 2018; Quintela et al., 2020; Saha et al., 2017; Westgaard et al., 2017), cost-effective fisheries management enforcement (Glover, 2010; Martinsohn et al., 2019), “real-time” regulation of harvest (Dahle et al., 2018a; Johansen et al., 2018) and harmonizing biological and management units (Aguirre-Sarabia et al., 2021; Leone et al., 2019; Rodríguez-Ezpeleta et al., 2019), which is particularly important for transboundary stocks of highly migratory and commercial important fish. Microsatellites are one of such tool, which consist of short repeating sequences of DNA found in organisms (O’Connell & Wright, 1997; Ellegren, 2004). They are highly variable due to their high mutation rate and are ideal for genetic studies to assess genetic diversity, identify population structure, and conduct parentage analysis (Ramya & Behera, 2023).

Widely distributed species with shared/transboundary distributions (Munro, Van Houtte & Willmann, 2004) such as the Atlantic chub mackerel, Scomber colias Gmelin, 1789, are particularly challenging to manage due to the complexities of their stock structure, but make them a well-suited candidate for genetic analysis. S. colias, a medium-size pelagic schooling species belonging to the Scombridae family, inhabits the continental shelf down to 300 m depth throughout the Atlantic Ocean, the Mediterranean and the Black Seas (Collette, 1986; Collette & Nauen, 1983; Scoles, Graves & Collette, 1988). In the Northeast Atlantic, the species is distributed from southern European waters (Iberian Peninsula) southwards along the African coast to South of Africa, and is also found in the connected Mediterranean and Black Seas (Castro Hernández & Santana Ortega, 2000). The species is economically valuable and targeted by industrial and artisanal fisheries throughout its geographical distribution (Castro Hernández & Santana Ortega, 2000).

Various studies suggested significant differences between Atlantic chub mackerel from the Western and the Eastern Atlantic, supporting the existence of two different populations occurring on either side of the Atlantic Ocean (Costa et al., 2011; Roldán et al., 2000; Scoles, Graves & Collette, 1988). In the NE Atlantic and surrounding waters, the available genetic evidence suggests significant lack of genetic diversity in S. colias indicating the presence of a large panmictic unit (Rodríguez-Ezpeleta et al., 2016). According to Zardoya et al. (2004), the 5′ end of the mitochondrial control region indicates that S. colias behaves as a panmictic population in the Mediterranean Sea and the south of Portugal. Based upon microsatellites, Stroganov et al. (2023) found that in the southern region of Morocco (23°N–21°N), samples show a high genetic homogeneity, suggesting large degree of connectivity. Congeneric species such as S. scombrus reveal population differentiation only at a transatlantic scale (Nesbø et al., 2000), whereas S. japonicus displays weak genetic structure across the NW Pacific (Cheng et al., 2015), with local genetic differentiation seemingly linked to differences in spawning time and migratory behaviour (Zeng, Cheng & Chen, 2012).

Migration patterns could be potentially impacting population structure. Thus, in both African and European waters, the species would migrate from the wintering areas (mainly located off Mauritania, South Portugal and the inner part of the Bay of Biscay) towards northern areas in summer and, in the case of the Bay of Biscay, also towards the western Iberian Peninsula (ICES, 2021). Recent studies investigating the spatial variability of life-history parameters of S. colias suggested the existence of latitudinal trends with a mixing zone between Atlantic African, Mediterranean and Atlantic Iberian population components (e.g., Domínguez-Petit et al., 2022). In the most recent study using otolith shape analysis, Sbiba et al. (2024) pointed out the existence of two populations, a northern one in Morocco from Larache (36°N) to Tarfaya (28°N), and a southern one from Tarfaya to Senegal (14°N). Thus, two stocks have been outlined according to biological traits: a northern one ranging from North Morocco to Cape Bojador in South Morocco (26°07′37″N, 14°29′57″W), and a southern one between Cape Bojador and the South of Senegal (https://firms.fao.org/firms/resource/10100/en). However, the lack of both robust evidence of stock identity (ICES, 2021) as well as information on population connectivity led to a joint assessment of both stocks in this area since 2003, thus highlighting the need for stock identity information in the NE Atlantic part of the species range (FAO, 2020). Morphometric differences on S. colias were also detected across five sampling sites in Morocco (Bouzzammit & El Ouizgani, 2019) as well as across three different locations in Tunisia (Allaya et al., 2016). Moreover, Erguden et al. (2009) used morphometric and meristic analyses to distinguish stocks of S. colias in seven locations throughout the Black, Marmara, Aegean, and Northeastern Mediterranean Seas, revealing the existence of two distinct groups. Overall, these studies also highlight the high potential of phenotypic plasticity of this species.

The migratory behaviour of S. colias, along with the high dispersal potential of its early life stages (eggs and larvae) within the open Atlantic Ocean is expected to favour connectivity and act as a homogenization factor. In addition, the West African coast is affected by several oceanic currents (see Fig. 1), i.e., the Canary, Guinea, and Benguela Currents that can influence both dispersal and migration patterns, thus amplifying the homogenization effects (White et al., 2010). Consequently, a pattern of genetic panmixia is foreseen on S. colias along the west coast of Africa. To verify this hypothesis, eight microsatellites were used to investigate the population structure of this species sampled from south of Morocco (~28°N) to Namibia (~22°S) at 33 locations.

10.7717/peerj.17928/fig-1 Figure 1 Map of sampling locations of Scomber colias in the West coast of Africa.

Detailed information per sample can be found in Table 1. The arrows, based upon Fig. 1 in Reid et al. (2016), indicate some of the major currents affecting the area and their relative temperature: Canary Current (CC), North Equatorial Current (NEC), Equatorial Counter Current (ECC), Guinea Current (GC), Angolan Current (AnC), and South Equatorial Current (SEC).

10.7717/peerj.17928/table-1 Table 1 Summary statistics of Scomber colias genotyped at eight microsatellites.

Country	Sample	Latitude	Longitude	N	No alleles	Ar (N > 20)	No private alleles	H o	u He	F IS	Dev HWE (FDR)	Dev LD (FDR)	
Morocco	MOR_1	27.39	−13.38	51	128	12.1	0	0.791 ± 0.023	0.879 ± 0.013	0.090 ± 0.028	3 (2)	2 (0)	
	MOR_2	25.98	−14.87	58	135	12.4	2	0.803 ± 0.033	0.881 ± 0.015	0.081 ± 0.032	1 (1)	3 (1)	
	MOR_3	24.71	−15.8	45	130	12.6	0	0.780 ± 0.040	0.885 ± 0.011	0.108 ± 0.047	5 (4)	6 (2)	
	MOR_4	24.4	−15.93	41	126	12.5	0	0.832 ± 0.028	0.885 ± 0.011	0.048 ± 0.029	5 (3)	7 (1)	
	MOR_5	24.02	−16.28	44	117	11.8	0	0.763 ± 0.041	0.878 ± 0.014	0.122 ± 0.045	1 (1)	2 (0)	
	MOR_6	23.32	−16.6	40	126	12.4	1	0.794 ± 0.037	0.873 ± 0.016	0.079 ± 0.036	3 (3)	0 (0)	
	MOR_7	22.35	−17	51	126	12.4	1	0.823 ± 0.039	0.889 ± 0.009	0.064 ± 0.043	2 (1)	10 (3)	
	MOR_8	21.06	−17.31	38	121	12.9	0	0.841 ± 0.031	0.891 ± 0.012	0.043 ± 0.032	4 (4)	3 (0)	
	MOR_9	21.09	−17.42	39	125	12.7	1	0.811 ± 0.051	0.884 ± 0.013	0.072 ± 0.054	2 (2)	1 (0)	
Mauritania	MAU_1	20.22	−17.58	41	124	12.5	0	0.803 ± 0.037	0.882 ± 0.013	0.078 ± 0.039	3 (1)	3 (2)	
	MAU_2	19.46	−16.98	41	122	12.4	0	0.843 ± 0.041	0.884 ± 0.012	0.034 ± 0.043	4 (3)	0 (0)	
	MAU_3	19.05	−16.66	34	118	12.8	0	0.777 ± 0.055	0.881 ± 0.014	0.106 ± 0.056	2 (2)	1 (0)	
	MAU_4	18.08	−16.43	40	123	12.5	1	0.839 ± 0.037	0.879 ± 0.012	0.035 ± 0.034	1 (1)	0 (0)	
	MAU_5	16.73	−16.66	44	122	11.9	0	0.812 ± 0.036	0.872 ± 0.016	0.059 ± 0.036	1 (0)	0 (0)	
	MAU_6	17.41	−16.27	23	101	12.1	0	0.790 ± 0.041	0.866 ± 0.025	0.062 ± 0.052	2 (0)	1 (0)	
Senegal	SEN_1	14.56	−17.35	44	129	12.2	4	0.778 ± 0.035	0.869 ± 0.015	0.095 ± 0.033	2 (1)	1 (0)	
	SEN_2	13.16	−17.41	40	124	12.5	0	0.831 ± 0.023	0.882 ± 0.013	0.045 ± 0.024	3 (2)	2 (0)	
Gambia	GAM	13.25	−17.47	22	102	12.4	1	0.806 ± 0.052	0.877 ± 0.022	0.061 ± 0.052	2 (2)	1 (1)	
Guinea-Bissau	GUIBIS	10.92	−16.86	5	46	n.a.	0	0.725 ± 0.053	0.850 ± 0.023	0.034 ± 0.087	0 (0)	0 (0)	
Guinea	GUI_1	9.38	−15.61	29	106	11.8	0	0.773 ± 0.038	0.869 ± 0.021	0.095 ± 0.036	3 (3)	1 (0)	
	GUI_2	9.51	−15.22	24	102	12.1	0	0.802 ± 0.037	0.871 ± 0.014	0.058 ± 0.043	2 (1)	0 (0)	
Sierra Leone	SIER	8.06	−14.09	11	71	n.a.	1	0.839 ± 0.054	0.866 ± 0.019	−0.011 ± 0.054	0 (0)	0 (0)	
Liberia	LIB	4.66	−8.74	33	116	12.7	0	0.814 ± 0.048	0.894 ± 0.010	0.075 ± 0.053	1 (1)	1 (0)	
Ivory Coast	IVCO	4.81	−5.13	23	96	11.5	1	0.828 ± 0.027	0.877 ± 0.014	0.034 ± 0.029	3 (2)	3 (0)	
Ghana	GHA_1	5.19	−0.6	51	142	13.0	1	0.820 ± 0.025	0.886 ± 0.013	0.065 ± 0.025	2 (1)	2 (0)	
	GHA_2	4.59	−2.19	43	124	12.3	0	0.806 ± 0.034	0.883 ± 0.012	0.077 ± 0.036	2 (1)	0 (0)	
	GHA_3	4.92	−0.89	30	116	12.9	2	0.806 ± 0.032	0.881 ± 0.015	0.070 ± 0.032	2 (1)	0 (0)	
	GHA_4	4.86	−0.92	35	114	11.9	0	0.798 ± 0.041	0.877 ± 0.015	0.076 ± 0.047	0 (0)	0 (0)	
Gambia	GAB_1	−2.41	9.06	38	137	13.6	1	0.819 ± 0.023	0.884 ± 0.013	0.061 ± 0.019	0 (0)	1 (0)	
	GAB_2	−3.45	9.87	24	104	12.2	0	0.813 ± 0.036	0.879 ± 0.014	0.055 ± 0.042	2 (1)	1 (0)	
Angola	ANG_1	−10.64	13.52	27	111	12.5	0	0.792 ± 0.041	0.875 ± 0.019	0.077 ± 0.047	2 (1)	0 (0)	
	ANG_2	−15.06	12.02	38	114	12.2	0	0.768 ± 0.057	0.873 ± 0.016	0.109 ± 0.061	1 (1)	0 (0)	
Namibia	NAM	−22.21	13.17	22	96	11.6	0	0.733 ± 0.049	0.861 ± 0.019	0.129 ± 0.051	1 (0)	0 (0)	
Note:

Sampling sites per country, geographic coordinates (decimal degrees), number of individuals (N), total number of alleles, allelic richness based on a minimum sample of 20 diploid individuals, number of private alleles, observed heterozygosity Ho (mean ± SE), unbiased expected heterozygosity uHE (mean ± SE), inbreeding coefficient FIS (mean ± SE), number of deviations from Hardy-Weinberg equilibrium (HWE) at α = 0.05, number of deviations from Linkage Disequilibrium (LD) at α = 0.05 with False Discovery Rate (FDR) correction.

Materials and Methods

Sampling

Fish were collected during the scientific surveys on board the R/V Dr. Fridtjof Nansen carried out in 2017 and 2019 in the area stretching from Morocco (~28°N) to Namibia (~22°S) (Fig. 1). Sampling permits were obtained through the “Application of the Ecosystem Approach to Fisheries management considering climate and pollution impacts” (GCP/GLO/690/NOR). Pelagic trawl hauls were performed after echosounders target identification of fish schools and, when possible, 30 individuals were collected per trawling station (Table 1). The sampling procedure reflected the species’ distribution (Baird, 1977; Scoles, Graves & Collette, 1988), whose southern limit of highest abundance seems to be placed at low latitudes in the southern hemisphere. Thus, some 75.5% of the individuals were collected north to Ghana whereas the southernmost countries such as Angola and Namibia only accounted for 7.5% of the individuals in line with the species biomass levels observed in those zones. Fin clips from 1,848 fish were collected and preserved in ethanol prior to DNA extraction.

DNA isolation and genotyping

DNA was extracted using the Qiagen DNeasy 96 Blood & Tissue Kit in (Qiagen, Hilden, Germany) 96-well plates; each of which contained two or more negative controls. DNA concentration was quantified using a NanoDrop 8000 (Thermo Fisher Scientific, Waltham, MA, USA). The molecular markers used in the current study were initially developed for both S. colias (one) and S. japonicus (seven) (Catanese et al., 2010; Chen et al., 2017; Zeng & Cheng, 2012) as microsatellite cross-species amplification has proven a successful tool for fisheries management and conservation (e.g., Maduna et al., 2014) formerly applied to Scomber species (Tang et al., 2009; Zeng, Cheng & Chen, 2012). Eight microsatellite loci isolated from S. colias (Sco2_1) and S. japonicus (SJ78, SJNT19, SJT5, SJT53, SJT122, SJT182, SJT199), respectively (Table S1 in Supporting Information), were genotyped in three multiplexed reactions. PCR amplifications were performed in a final volume of 10 μL containing approximately 20 ng DNA template, 1 × buffer, 2–3 mM MgCl2, 0.2 mM dNTPs, 0.2–0.5 μM of each primer and 1U GoTaq polymerase in an Applied Biosystems Gene Amp PCR Systems 2700 thermal cycler. PCR profiles consisted of 15 min denaturation at 95 °C followed by 27–32 cycles of 30 s denaturation at 94 °C, 90 s at annealing temperature (Table S1 in Supporting Information) and 90 s extension at 72 °C with a final step of 7 min at 72 °C. PCR products were analysed on an ABI 3730 Genetic Analyser and the 500LIZTM size standard (Applied Biosystems, Waltham, MA, USA) was used to accurately determine the size of the fragments and allelic variation. Scoring was performed using the software Fragment Profiler 1.2 (Amersham Biosciences, Buckinghamshire, England). Automatically binned alleles were manually checked by two researchers prior to exporting data for statistical analyses.

As it was not always possible to obtain a minimum of thirty good-quality genotyped individuals per site, nearby sampling sites within countries were merged to achieve significant sampling sizes (Hale, Burg & Steeves, 2012). Sampling sizes ranged from 22 individuals (GAM in Gambia) to 58 (MOR_2 in Morocco). However, in Sierra Leone (N = 11), and Guinea-Bissau (N = 5) (Table 1), sampling sizes were too low and therefore these samples were excluded from some of the analyses.

Statistical analysis

Loci were screened for null alleles, large allele dropouts and potential scoring errors due to stuttering bands using the software Micro-Checker v.2.2.3 (Van Oosterhout et al., 2004). The frequency of null allele(s) was estimated with the maximum likelihood method using the EM algorithm of Dempster, Laird & Rubin (1977) implemented in the software Genepop 7 (Rousset, 2008). To assess the effect of including loci with possible null allele(s) on population differentiation estimates, the software FreeNA (Chapuis & Estoup, 2006) was used to calculated both uncorrected and corrected FST values. Confidence intervals (95%) of null frequencies were based on 1,000 bootstraps.

The statistical power of this microsatellites set to detect genetic differentiation was estimated using POWSIM 4.1 (Ryman & Palm, 2006), which uses χ2 and Fisher tests to assess whether the observed data set carries enough statistical power to detect a Nei’s FST value significantly larger than zero. The analyses were conducted maintaining effective population size constant (Ne = 15,000) following Henriques et al. (2017), and changing the number of generations of drift since divergence (t: 15, 35, 100, 160, 200, 350 and 500 generations) to model different levels of FST (0.0005, 0.0012, 0.0033, 0.0053, 0.0066, 0.0117, and 0.0165, respectively). The probability of rejecting the null hypotheses of no population differentiation (type I error) was also estimated (FST = 0, t = 0). Statistical power was determined as the proportion of tests indicating statistical significance (p < 0.05), with 1,000 replicates. In addition, to assess if this suite of twelve microsatellites would accurately discriminate between individuals in a population, the genotype accumulation curve was built using the function genotype curve in the R (R Core Team, 2020) package poppr (Kamvar, Tabima & Grünwald, 2014) by randomly sampling x loci without replacement and counting the number of observed multilocus genotypes (MLGs).

The total number of alleles, number of private alleles, and allelic richness per sample was calculated using MSA 4.05 (Dieringer & Schlötterer, 2003). The observed (Ho) and unbiased expected heterozygosity (uHE), inbreeding coefficient (FIS) as well as the number of deviations from Hardy-Weinberg expectations (HWE) were computed per sample using GenAlEx v6.1 (Peakall & Smouse, 2006). Linkage disequilibrium (LD) between pairwise loci per sample was computed using the program Genepop (Rousset, 2008). The false discovery rate (FDR) correction of Benjamini & Hochberg (1995) was applied to p-values to control for Type I errors.

In marine fish species, loci carrying signatures of locally divergent selection might function as powerful markers to evaluate spatially explicit genetic structure, as well as to outline fisheries stocks for sustainable management (Russello et al., 2012). Thus, two approaches were combined to identify loci eventually departing from neutrality: BayeScan 2.1 (Foll & Gaggiotti, 2008) and Arlequin v.3.5.1.2 (Excoffier, Laval & Schneider, 2005). In BayeScan, sample size was set to 10,000 and thinning interval to 50; and loci with a posterior probability over 0.99, corresponding to a Bayes Factor >2 (i.e., “decisive selection” (Foll & Gaggiotti, 2006), were retained as outliers. In Arlequin, analysis was simulated based on 1,000 demes with 50,000 simulations under a hierarchal island model.

Pairwise FST (Weir & Cockerham, 1984) was computed using GENEPOP 4.0.6 (Rousset, 2008), whereas population differentiation was tested via G exact tests calculated using the following Monte Carlo Markov Chain (MCMC) parameters: 10,000 steps of dememorisation, and 5,000 iterations for 100 batches, also with GENEPOP 4.0.6.

The Bayesian clustering approach implemented in STRUCTURE v.2.3.4 (Pritchard, Stephens & Donnelly, 2000), and conducted using the software ParallelStructure (Besnier & Glover, 2013), was used to identify genetic groups under a model assuming admixture and correlated allele frequencies, both using and without using geographic LOCPRIORS to assist the clustering. Ten runs with a burn-in period consisting of 100,000 replications and a run length of 1,000,000 MCMC iterations were performed for K = 1 to K = 10 clusters. To determine the number of genetic groups, STRUCTURE output was analysed using two approaches: (a) the ad hoc summary statistic ΔK of Evanno, Regnaut & Goudet (2005) and (b) the Puechmaille (2016) four statistics (MedMedK, MedMeanK, MaxMedK and MaxMeanK), specially recommended for uneven sampling sizes, both conducted using StructureSelector (Li & Liu, 2018). Finally, the ten runs for the selected Ks were averaged with CLUMPP v.1.1.1 (Jakobsson & Rosenberg, 2007) using the FullSearch algorithm and the G′ pairwise matrix similarity statistic, and graphically displayed using barplots. Furthermore, the relationship among samples was also examined using the Discriminant Analysis of Principal Components (DAPC) (Jombart, Devillard & Balloux, 2010) implemented in the R (R Core Team, 2020) package adegenet (Jombart, 2008) in which groups were defined a priori using geographically explicit samples. To avoid overfitting, both the optimal number of principal components and discriminant functions to be retained were determined through cross validation using the xvalDapc function from adegenet (Jombart & Collins, 2015; Miller, Cullingham & Peery, 2020).

The relationship between genetic (FST) and geographic (km) distance was examined to investigate if it followed the expectations of an “Isolation by Distance” pattern (IBD), i.e., increasing genetic differentiation with geographic distance as a result of restricted gene flow and drift (Rousset, 1997; Slatkin, 1993; Wright, 1943). A two-tailed Mantel (1967) test was conducted using PASSaGE v2 (Rosenberg & Anderson, 2011) and significance was assessed via 10,000 permutations. The matrix of pairwise shortest distance by water was created by calculating least-cost distances via seas (avoiding landmasses) between sampling sites using the lc.dist function from the R (R Core Team, 2020) package marmap v1.0 (Pante & Simon-Bouhet, 2013).

Results

Quality control

Some 679 individuals out of the 1,848 genotyped ones were discarded due to exceeding the threshold of acceptance of missing genotypic data, which was set at 25% or due to unreliable scoring. From the 1,169 retained individuals, no missing markers were reported for 84% of the fish, whereas only 4% of the individuals displayed the maximum missing data allowed (see distribution of missing data per sample and locus in Fig. S1 in Supporting Information).

Data validation conducted with Micro-Checker did not find evidence of scoring error due to stuttering or large allelic dropout but suggested the possibility of null alleles in two of the loci, which showed a significant deficit of heterozygotes (p < 0.0001) after global Hardy-Weinberg tests. The frequency of null alleles per sample in these loci took a maximum value of 0.216 (locus STJ5 in ANG_2), whereas it was <0.05 in 71% of the cases (see Table S2 in Supporting information). As global and pairwise divergence calculated with and without using the ENA correction implemented in FreeNA displayed similar values (i.e., a difference of 0.00045 in the overall FST values); the two loci suggestive of displaying null alleles were retained for further analyses.

The simulation analyses revealed that the dataset carried enough statistical power to detect genetic differentiation as low as FST = 0.0066 in 80% of the cases, whereas only in 11% of the cases when levels of differentiation were five-fold lower (FST = 0.0012). The probability of rejecting the null hypothesis of panmixia if true (type I error) was estimated at 5% (Fisher test) and thus considered appropriate (Fig. S2 in Supporting Information). In addition, the number of microsatellite markers with the capacity of differentiating unique individuals can be inferred from the plateau of the genotype accumulation curve and this was achieved with some 37% of the markers used (Fig. S3 in Supporting Information).

Summary statistics

The eight loci showed high rates of amplification success (93–100%) and suitable levels of polymorphism for population genetic analysis (100% polymorphic loci per sample). Allelic richness based on a minimum sample size of 20 diploid individuals took similar values throughout all the geographic range and moved between 11.5 (IVCO) and 13.6 (GAB_1), whereas the number of alleles for the same set of samples ranged between 96 and 137, respectively (Table 1). A total of 17 private alleles were found; they were distributed across all loci but SJT122 and ranged from 1 (Sco2-1, SJ78) to 4 (SJNT19, SJT199 and SJT53) and across 12 samples ranging from 1 (MOR_6, MOR_7, MOR_9, MAU_4, GAM, SIER, IVCO,GHA_1, GAB_1) to 4 (SEN_1). Ho and uHE took similar values across samples and in the range 0.73–0.9. uHE was consistently higher than HO in all samples; in consequence, multilocus test detected heterozygote deficit in every sample but SIER (p-val = 0.257). The three loci displaying deviations from HWE in a larger number of samples after correcting for FDR were: SJT5 (in 20 samples), and Sco2-1 and SJ78 (in six). Both SJT5 and Sco2-1 exhibited homozygote excess and likely presence of null alleles.

Genetic differentiation

No locus was found to be deviate from neutral expectations according to Arlequin, whereas BayeScan flagged all of them as putative candidates to balancing selection (Table S3 in Supporting Information). Due to the lack of consensus of neutrality between methods, the full set of microsatellites was retained for subsequent analyses.

Overall, no significant population structuring was detected across the study area as indicated by a global FST of 0.001 (p = 0.587). Pairwise FST ranged from 0.000 to 0.011 (Table 2) with no significant differentiation whatsoever across the investigated geographic range.

10.7717/peerj.17928/table-2 Table 2 Tests of genetic differentiation of S. colias samples collected from Northeast Atlantic waters of West Africa between Morocco and Namibia.

Pairwise FST values (bottom diagonal) and significance of G exact tests (upper diagonal).

	MOR
_1	MOR
_2	MOR
_3	MOR
_4	MOR
_5	MOR
_6	MOR
_7	MOR
_8	MOR
_9	MAU
_1	MAU
_2	MAU
_3	MAU
_4	MAU
_5	MAU
_6	SEN
_1	SEN
_2	GAM	GUI
_1	GUI
_2	LIB	IVCO	GHA
_1	GHA
_2	GHA
_3	GHA
_4	GAB
_1	GAB
_2	ANG
_1	ANG
_2	NAM	
MOR
_1		1.000	1.000	0.468	0.454	1.000	0.538	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.574	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MOR
_2	0.000		1.000	0.504	0.504	1.000	0.490	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.599	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MOR
_3	0.000	0.000		0.512	0.490	1.000	0.498	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.550	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MOR
_4	0.000	0.001	0.000		0.244	0.483	0.264	0.483	0.493	0.458	0.501	0.516	0.524	0.521	0.514	0.526	0.493	0.524	0.486	0.504	0.499	0.281	0.472	0.495	0.498	0.515	0.500	0.515	0.518	0.488	0.529	
MOR
_5	0.000	0.000	0.000	0.000		0.489	0.246	0.493	0.507	0.503	0.509	0.488	0.506	0.499	0.532	0.480	0.496	0.559	0.502	0.537	0.514	0.301	0.493	0.504	0.529	0.490	0.451	0.551	0.527	0.471	0.521	
MOR
_6	0.000	0.000	0.000	0.000	0.000		0.556	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.546	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MOR
_7	0.000	0.000	0.002	0.000	0.000	0.002		0.482	0.477	0.496	0.513	0.523	0.500	0.487	0.559	0.526	0.479	0.589	0.515	0.548	0.543	0.326	0.498	0.503	0.544	0.482	0.503	0.573	0.496	0.507	0.593	
MOR
_8	0.000	0.001	0.000	0.000	0.000	0.001	0.002		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.526	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MOR
_9	0.000	0.000	0.001	0.000	0.000	0.000	0.000	0.003		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.510	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MAU
_1	0.000	0.000	0.000	0.002	0.000	0.001	0.004	0.002	0.001		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.514	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MAU
_2	0.000	0.003	0.002	0.002	0.001	0.000	0.000	0.000	0.000	0.004		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.540	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MAU
_3	0.000	0.000	0.000	0.000	0.003	0.001	0.002	0.000	0.000	0.000	0.002		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.493	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MAU
_4	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.540	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MAU
_5	0.000	0.000	0.000	0.001	0.003	0.000	0.004	0.000	0.001	0.001	0.000	0.000	0.000		1.000	1.000	1.000	1.000	1.000	1.000	1.000	0.551	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
MAU
_6	0.002	0.001	0.008	0.005	0.006	0.004	0.005	0.010	0.007	0.005	0.006	0.004	0.008	0.002		1.000	1.000	1.000	1.000	1.000	1.000	0.484	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
SEN
_1	0.000	0.000	0.001	0.001	0.004	0.002	0.003	0.002	0.002	0.000	0.004	0.000	0.002	0.000	0.001		1.000	1.000	1.000	1.000	1.000	0.570	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
SEN
_2	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.001	0.002	0.000	0.000	0.001	0.005	0.001		1.000	1.000	1.000	1.000	0.548	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
GAM	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.002	0.000	0.001	0.000	0.000	0.000	0.000	0.000	0.000		1.000	1.000	1.000	0.499	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
GUI
_1	0.000	0.000	0.000	0.000	0.000	0.000	0.001	0.001	0.002	0.001	0.003	0.001	0.000	0.000	0.000	0.003	0.000	0.000		1.000	1.000	0.496	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
GUI
_2	0.000	0.003	0.004	0.006	0.004	0.005	0.008	0.003	0.007	0.001	0.005	0.007	0.003	0.002	0.003	0.003	0.007	0.000	0.010		1.000	0.508	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
LIB	0.000	0.003	0.001	0.001	0.002	0.003	0.004	0.001	0.003	0.002	0.004	0.000	0.005	0.004	0.010	0.007	0.000	0.001	0.006	0.005		0.521	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
IVCO	0.000	0.001	0.004	0.000	0.000	0.001	0.002	0.000	0.000	0.001	0.000	0.000	0.000	0.000	0.004	0.002	0.000	0.000	0.000	0.004	0.005		0.591	0.531	0.478	0.510	0.538	0.501	0.483	0.496	0.499	
GHA
_1	0.000	0.000	0.001	0.000	0.002	0.006	0.000	0.000	0.002	0.000	0.003	0.003	0.001	0.003	0.005	0.000	0.002	0.000	0.003	0.004	0.005	0.003		1.000	1.000	1.000	1.000	1.000	1.000	1.000	1.000	
GHA
_2	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.004	0.000	0.000	0.000	0.000	0.002	0.002	0.000	0.000		1.000	1.000	1.000	1.000	1.000	1.000	1.000	
GHA
_3	0.000	0.000	0.004	0.001	0.000	0.001	0.000	0.004	0.003	0.002	0.004	0.003	0.001	0.002	0.000	0.000	0.000	0.000	0.000	0.002	0.005	0.002	0.002	0.001		1.000	1.000	1.000	1.000	1.000	1.000	
GHA
_4	0.000	0.001	0.000	0.001	0.000	0.000	0.004	0.002	0.002	0.002	0.004	0.002	0.000	0.002	0.009	0.004	0.003	0.000	0.003	0.006	0.002	0.004	0.001	0.000	0.005		1.000	1.000	1.000	1.000	1.000	
GAB
_1	0.000	0.001	0.001	0.001	0.001	0.001	0.000	0.001	0.002	0.005	0.001	0.003	0.000	0.004	0.008	0.005	0.000	0.001	0.001	0.008	0.004	0.000	0.000	0.001	0.002	0.002		1.000	1.000	1.000	1.000	
GAB
_2	0.000	0.000	0.003	0.000	0.000	0.000	0.000	0.000	0.001	0.003	0.000	0.006	0.000	0.003	0.008	0.002	0.000	0.000	0.001	0.004	0.005	0.000	0.000	0.000	0.000	0.005	0.000		1.000	1.000	1.000	
ANG
_1	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.001	0.000	0.000	0.000	0.000	0.000	0.001	0.006	0.003	0.000	0.000	0.000	0.008	0.006	0.000	0.000	0.000	0.002	0.000	0.000	0.000		1.000	1.000	
ANG
_2	0.000	0.001	0.004	0.006	0.003	0.002	0.003	0.001	0.003	0.004	0.000	0.004	0.000	0.001	0.008	0.003	0.002	0.000	0.001	0.003	0.010	0.000	0.001	0.000	0.001	0.002	0.000	0.000	0.000		1.000	
NAM	0.000	0.000	0.000	0.001	0.002	0.001	0.000	0.000	0.000	0.001	0.001	0.000	0.000	0.000	0.006	0.001	0.000	0.000	0.000	0.011	0.003	0.001	0.003	0.000	0.001	0.004	0.000	0.000	0.000	0.000		

STRUCTURE conducted both with and without LOCPRIORS reported a decreasing trend of LnP(D) across subsequent values of K with the highest average likelihood at K = 1 (Figs. S4A and S4B in Supporting Information). Evanno test, which by definition cannot yield K = 1, reported K = 2 (no priors) and K = 5 (LOCPRIORS), respectively but with extremely little support as both ΔK values were low (ΔK < 7). Puechmaille (2016)’s four statistics (MedMedK, MedMeanK, MaxMedK and MaxMeanK) implemented in StructureSelector confirmed one single genetic group in the model using priors, and K = 2 when not using them (Figs. S4A and S4B in Supporting Information). The barplot for K = 2 for the model without using LOCPRIORS revealed inferred ancestry to cluster ranging around 0.5 for each individual (Fig. S5 in Supporting Information), which is consistent with one single genetic unit.

Cross validation determined that 200 was the optimal number of PCs to be retained for the DAPC analysis using the 33 geographically explicit locations. In agreement with STRUCTURE, the corresponding DAPC plot revealed a large overlap of the individuals sampled across all the studied geographic range (Fig. 2) with none of the axes reaching 7% of the variation. Although some some individuals of SEN_1 partially deviated from the cloud, SEN_1 was not significantly different from any of the remaining samples according to pairwise FST.

10.7717/peerj.17928/fig-2 Figure 2 Genetic differentiation among Scomber colias sampled in the West African coast and genotyped at eight microsatellite markers depicted using Discriminant Analysis of Principal Components (DAPC) after retaining 200 principal components.

Individuals from different samples are represented by coloured dots, and name labels are placed on the centroids of each geographically-explicit site.

Mantel test revealed lack of correlation between genetic differentiation (pairwise FST) and geographic distance along the latitudinal gradient ranging from Morocco to Namibia (rxy = 0.012, p = 0.401).

Discussion

Investigating transboundary fish stocks with genetics is crucial for sustainable fisheries management. By identifying population structure in important commercial fish species, we can prevent overexploitation and allow for coordinated management strategies across national borders. This is the first attempt to describe the population structure of the Atlantic chub mackerel over most of its African distribution. The main result of the present study is the lack of geographically explicit structure in the coastline stretching from Morocco (27.39°N) to Namibia (22.21°S) based upon eight microsatellite markers.

Genetic diversity

Many marine fishes exhibit high levels of genetic diversity due to historically large effective population sizes (Hauser et al., 2002; Hauser & Carvalho, 2008). All samples analysed in the current study showed high levels of genetic diversity in terms of allelic richness (11.5 to 13.6) and heterozygosity (≥0.7) in agreement with results obtained for other medium and small pelagic species, e.g., S. australasicus in the western North Pacific (Tzeng et al., 2009), European sprat Sprattus sprattus (Glover et al., 2011), European anchovy Engraulis encrasicolus (Zarraonaindia et al., 2009) or European pilchard Sardina pilchardus (Baibai et al., 2012) in the North-East Atlantic. Overall, all samples (with the exception of the one from Sierra Leone represented by just eleven individuals) deviated from HWE expectations showing a deficit of heterozygotes. Heterozygote deficiency is not new to Scomber literature as it has been shown in microsatellite studies of S. colias in the East Atlantic Ocean and the Mediterranean Sea (Medina-Alcaraz, 2016) as well as in S. japonicus along the Chinese coast (Cheng, Zhu & Chen, 2014; Zeng, Cheng & Chen, 2012). Heterozygote deficiency in non-selfing, diploid populations is relatively common (Waldman & McKinnon, 1993) and in largely outcrossing species, it has commonly been interpreted as indirect evidence of the mixing of differentiated gene pools (Wahlund effect). However, selection and the presence of null alleles are common cause of heterozygote deficiency. In agreement with the latter, in the current study, two out of the eight genotyped loci were likely to display null alleles according to Micro-Checker assessment. Finally, and given the S. colias stocks are considered to be overexploited along West Africa (FAO, 2020), the third hypothesis to explain heterozygote deficiency could be overexploitation as overfishing reduces population size drastically, leading to genetic bottlenecks and loss of alleles (Smith, 1994).

Genetic population structure

In this study, the lack of genetic differentiation in S. colias revealed by eight microsatellites in the geographic area from Morocco (27.39°N) to Namibia (22.21°S) supporting a theory of panmixia would be a foreseen outcome due to pertinent ecological features of pelagic fish, such as high fecundity, large effective population sizes that impose limitations to genetic drift as well as a notable dispersal capacity as planktonic larvae and free-swimming adults that promotes gene flow (e.g., Carlsson et al., 2004; Martínez et al., 2006; Waples, 1987, 1998). The combination of both features acts as a major homogenizing force hindering genetic differentiation and eventually leading to panmixia (Canales-Aguirre et al., 2018; De Bie et al., 2012; Roy et al., 2014; Weersing & Toonen, 2009).

S. colias is a highly migratory species with individuals capable of covering extensive distances during their annual migrations favoured by a body shape well fitted for swift motion and free swimming mode of life (Collette, 1999; Uriarte et al., 2001). In addition, its planktonic larval stages also display considerable dispersal capabilities as pelagic larvae can be transported over vast distances by ocean currents and eddies (Hernández-León, Gómez & Arístegui, 2007) effectively leading to genetic homogeneity among populations. In agreement, other studies have revealed a notable absence of genetic differentiation on S. colias in the Northeast Atlantic and adjacent waters suggesting a large panmictic unit (Rodríguez-Ezpeleta et al., 2016; Stroganov et al., 2023; Zardoya et al., 2004). Likewise, congeneric species such as S. scombrus reveal population differentiation only at a transatlantic scale (Nesbø et al., 2000), whereas S. japonicus displays weak genetic structure across the NW Pacific (Cheng et al., 2015).

The hydrological conditions along the West African coast play a crucial role in various ecological processes and are influenced by factors such as ocean currents, eddies, temperature gradients, and upwelling phenomena. The region experiences the influence of several major ocean currents (see Fig. 1), including the Canary Current, the Guinea Current, and the Benguela Current, which affect the distribution of nutrients and plankton, which in turn impacts the abundance and distribution of marine life, including fish populations and seem to shape genetic differentiation of coastal species in this region (Alpers et al., 2013; Barton, Field & Roy, 2013; Mittelstaedt, 1991; Nielsen et al., 2018; Pelegrí & Peña-Izquierdo, 2015; Shannon, 2001). In consequence, distant sites may be connected by a strong current between them, so, the absence of the genetic structure that can result from the long-distance dispersal of the early development stages (Vásquez et al., 2013; White et al., 2010).

However, regional differences in various life-history traits with latitudinal trends have been detected in both European and NW African waters, with the Strait of Gibraltar as changing point (ICES, 2021). Several studies have addressed the delineation of S. colias population structure in the East Atlantic using phenotypic characters, and some latitudinal trends seem to occur in length composition and life-history traits from the Strait of Gibraltar northwards into European waters and southwards, through African waters (Domínguez-Petit et al., 2022; ICES, 2021; Landa et al., 2022; Perrotta, Carvalho & Isidro, 2005). Differences include size at first maturity, which seems to be larger in the South of Africa than in the North (ICES, 2021). Likewise, the spawning period, a crucial evolutionary driver in pelagic species (e.g., Lamichhaney et al., 2017; Petrou et al., 2021), corresponds to January-March in the region comprising Morocco to Senegal, in contrast to June-September in the area from Ghana to South Africa (ICES, 2021). Regarding the Atlantic coast of North America, the population structure of the Atlantic chub mackerel has not been extensively studied; however, differences in growth (Daley & Leaf, 2019), habitat suitability (Chen et al., 2009), spawning seasons (Weber & McClatchie, 2011), size at maturity (Cerna & Plaza, 2014) and morphology (Erguden et al., 2009) suggest that different sub-stocks might exist. The latitudinal trends observed in the East Atlantic, with the Strait of Gibraltar as inflection point (ICES, 2021), indicate that significant differences might be expected if individuals from the whole East Atlantic distribution of S. colias (including Mediterranean samples) were jointly analysed. This study is lacking samples from North of Morocco as well as from the Iberian Peninsula, which should be included in future genetic assessments.

Management implications

The lack of geographically explicit structure found here suggests the existence of one single putative genetic stock for S. colias inhabiting the West African waters from Morocco (27.39°N) to Namibia (22.21°S). Panmixia in a transboundary stock can present both opportunities and challenges. First, if a transboundary stock truly acts as a single panmictic population, then a single stock assessment might be sufficient, and the complexity of modelling population dynamics and setting catch quotas becomes simplified. In this case, all countries sharing the stock would benefit from implementing consistent management practices and should work together to set catch quotas and implement regulations that benefit the entire population.

However, a determination of panmixia based on a limited dataset might mask the presence of subtle subpopulations with distinct spawning grounds or migration patterns. In this case, management might need to consider specific population components, which if unnoticed could lead to inadequate management strategies and potential overexploitation of vulnerable subpopulations.

Although the sampling strategy of the present study comprises most of the West African geographical distribution of the species, this did not include the entire NE Atlantic distribution and the seasonal migrations could not be considered in the study. Therefore, although genetic divergences are not noticeable, some kind of local adaptations seem to exist, more when seasonal migration processes of the species have been described in the study area (García, 1982). From this perspective, and considering the current global warming scenario in which the northwards expansion of the species across the East Atlantic Ocean from regions of higher abundance off northwest Africa to the waters of the Atlantic Iberian and the Mediterranean Sea is a fact (Jurado-Ruzafa et al., 2024), fisheries management of Atlantic chub mackerel should account not only with punctual assessment statuses of political-based stocks but consider all the variables mentioned along the present section.

Although our data suggests that a nearly panmictic model is the most plausible scenario for the genetic population pattern of S. colias, implying high level of genetic connectivity across the Western African countries ranging from 28°N to 22°S. However, phenotypic variability and life history of S. colias, which follows a latitudinal gradient, support the evidence of regional demographic units as suggested by Sbiba et al. (2024) based upon otolith shape. Therefore, relying solely on genetic data may not provide a complete snapshot of the demographic connectivity, especially in wild marine fish populations characterized by high migration rates and large population sizes. It is crucial to consider demographic connectivity in addition to genetic connectivity when managing an important commercial species to ensure effective and sustainable measures. Finally, it cannot be neglected that genetic differentiation could have eluded the scrutiny of the set of markers used here as well as the possible impact of ascertainment bias driven by markers isolated from sampled from the Pacific ocean (Zeng & Cheng, 2012). Future studies would benefit from the genomic tools currently available such as the genome assembly of this species recently published (Machado et al., 2022).

Conclusions

The findings of this study shed light on the population genetics of Scomber colias along the West African coast. The study indicates that the species has very low genetic population differentiation and large migratory (dispersal) potential. These results have significant implications for the management of this crucial pelagic species in this region. The study emphasizes the need for cooperation efforts and further research to fully comprehend the underlying mechanisms of population connectivity and genetic homogenization. Advanced genomic tools, such as whole-genome sequencing in combination with phenotypic and environmental information could provide a more comprehensive understanding of any complex population structure or adaptive potential of S. colias. Such an understanding is important for developing more rational effective management strategies and ensuring the sustainability of this shared stock. It could also help refine our understanding of this species’ evolutionary trajectories and guide sustainable management approaches at fine scale.

Supplemental Information

10.7717/peerj.17928/supp-1 Supplemental Information 1 Supplemental figures and tables.

10.7717/peerj.17928/supp-2 Supplemental Information 2 Allele frequency of the eight loci genotyped on Scomber colias.

Alleles identified by their size (bp) in the X-axis and number of individuals presenting each allele in the Y-axis.

We thank the captain and crew of R.V. “Dr. Fridtjof Nansen” for their invaluable assistance and contribution to this study, which made a significant impact on the research’s success. Likewise, we thank Alejandro Mateos-Rivera for his advice in the laboratory. The authors would like to extend their gratitude to the INRH for nominating Salah eddine Sbiba to embark on the R/V Dr. Fridtjof Nansen. This opportunity has significantly contributed to the success of our research. Disclaimer: The views expressed in this publication are those of the author(s) and do not necessarily reflect the views or policies of the Food and Agriculture Organization of the United Nations.

Additional Information and Declarations

Competing Interests

Author Contributions

Field Study Permissions

Data Availability

The authors declare that they have no competing interests.

Salah eddine Sbiba conceived and designed the experiments, performed the experiments, analyzed the data, authored or reviewed drafts of the article, and approved the final draft.

María Quintela analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Johanne Øyro performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Geir Dahle performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Alba Jurado-Ruzafa conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Kashona Iita conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Nikolaos Nikolioudakis conceived and designed the experiments, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Hocein Bazairi conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Malika Chlaida conceived and designed the experiments, authored or reviewed drafts of the article, and approved the final draft.

The following information was supplied relating to field study approvals (i.e., approving body and any reference numbers):

The EAF-Nansen Programme approved the study (GCP/GLO/690/NOR).

The following information was supplied regarding data availability:

Genotype data is available at Zenodo and Havforskningsinstituttet:

Sbiba, S. E. (2024). Microsatellite data on Scomber colias in NW Africa [Data set]. Zenodo. https://doi.org/10.5281/zenodo.11518098

https://imr.brage.unit.no/imr-xmlui/handle/11250/3072605
==== Refs
References

Aguirre-Sarabia et al. (2021) Aguirre-Sarabia I Díaz-Arce N Pereda-Agirre I Mendibil I Urtizberea A Gerritsen HD Burns F Holmes I Landa J Coscia I Quincoces I Santurtún M Zanzi A Martinsohn JT Rodríguez-Ezpeleta N Evidence of stock connectivity, hybridization, and misidentification in white anglerfish supports the need of a genetics-informed fisheries management framework Evolutionary Applications 2021 14 9 2221 2230 10.1111/eva.13278 34603494
Allaya et al. (2016) Allaya H ben faleh A Rebaya M Zrelli S Hajjej G Hattour A Jean-pierre Q Trabelsi M Identification of Atlantic chub mackerel Scomber colias population through the analysis of body shape in Tunisian waters Cahiers de Biologie Marine 2016 57 195 207 10.21411/CBM.A.6C4E0868
Allendorf et al. (2008) Allendorf FW England PR Luikart G Ritchie PA Ryman N Genetic effects of harvest on wild animal populations Trends in Ecology & Evolution 2008 23 6 327 337 10.1016/j.tree.2008.02.008 18439706
Alpers et al. (2013) Alpers W Brandt P Lazar A Dagorne D Sow B Faye S Hansen MW Rubino A Poulain PM Brehmer P A small-scale oceanic eddy off the coast of West Africa studied by multi-sensor satellite and surface drifter data Remote Sensing of Environment 2013 129 C8 132 143 10.1016/j.rse.2012.10.032
Baibai et al. (2012) Baibai T Oukhattar L Quinteiro JV Mesfioui A Rey-Méndez M soukri A First global approach: morphological and biological variability in a genetically homogeneous population of the European pilchard, Sardina pilchardus (Walbaum, 1792) in the North Atlantic coast Reviews in Fish Biology and Fisheries 2012 22 1 63 80 10.1007/s11160-011-9223-9
Baird (1977) Baird D Age, growth and aspects of reproduction of the mackerel, Scomber japonicus in South African waters (Pisces: Scombridae) African Zoology 1977 12 347 362 10.1080/00445096.1977.11447580
Barton, Field & Roy (2013) Barton ED Field DB Roy C Canary current upwelling: more or less? Progress in Oceanography 2013 116 1 167 178 10.1016/j.pocean.2013.07.007
Bekkevold et al. (2015) Bekkevold D Helyar SJ Limborg MT Nielsen EE Hemmer-Hansen J Clausen LAW Carvalho GR FishPopTrace Consortium Gene-associated markers can assign origin in a weakly structured fish, Atlantic herring ICES Journal of Marine Science 2015 72 6 1790 1801 10.1093/icesjms/fsu247
Benjamini & Hochberg (1995) Benjamini Y Hochberg Y Controlling the false discovery rate: a practical and powerful approach to multiple testing Journal of the Royal Statistical Society: Series B (Methodological) 1995 57 1 289 300 10.1111/j.2517-6161.1995.tb02031.x
Besnier & Glover (2013) Besnier F Glover KA ParallelStructure: a R package to distribute parallel runs of the population genetics program STRUCTURE on multi-core computers PLOS ONE 2013 8 e70651 10.1371/journal.pone.0070651 23923012
Besnier et al. (2014) Besnier F Kent M Skern-Mauritzen R Lien S Malde K Edvardsen RB Taylor S Ljungfeldt LER Nilsen F Glover KA Human-induced evolution caught in action: SNP-array reveals rapid amphi-atlantic spread of pesticide resistance in the salmon ecotoparasite Lepeophtheirus salmonis BMC Genomics 2014 15 1 937 10.1186/1471-2164-15-937 25344698
Bouzzammit & El Ouizgani (2019) Bouzzammit N El Ouizgani H Morphometric and meristic variation in the Atlantic chub mackerel Scomber colias Gmelin, 1789 from the Moroccan coast Indian Journal of Fisheries 2019 66 2 8 15 10.21077/ijf.2019.66.2.78488-02
Canales-Aguirre et al. (2018) Canales-Aguirre CB Ferrada-Fuentes S Galleguillos R Oyarzun FX Buratti CC Hernández CE High genetic diversity and low-population differentiation in the Patagonian sprat (Sprattus fuegensis) based on mitochondrial DNA Mitochondrial DNA Part A 2018 29 8 1148 1155 10.1080/24701394.2018.1424841
Carlsson et al. (2004) Carlsson J McDowell JR Díz-Jaimes P Carlsson JEL Boles SB Gold JR Graves JE Microsatellite and mitochondrial DNA analyses of Atlantic bluefin tuna (Thunnus thynnus thynnus) population structure in the Mediterranean Sea Molecular Ecology 2004 13 11 3345 3356 10.1111/j.1365-294X.2004.02336.x 15487994
Castro Hernández & Santana Ortega (2000) Castro Hernández JJ Santana Ortega AT Bonfil R Minow J Synopsis of biological data on the chub mackerel (Scomber japonicus Houttuyn, 1782) FAO Fisheries Synopsis 2000 Rome FAO 77
Catanese et al. (2010) Catanese G Funes V Perez L Infante C Microsatelites for Scomber colias 2010 https://www.researchgate.net/publication/281005067
Cerna & Plaza (2014) Cerna F Plaza G Life history parameters of chub mackerel (Scomber japonicus) from two areas off Chile Bulletin of Marine Science 2014 90 3 833 848 10.5343/bms.2013.1077
Chapuis & Estoup (2006) Chapuis MP Estoup A Microsatellite null alleles and estimation of population differentiation Molecular Biology and Evolution 2006 24 3 621 631 10.1093/molbev/msl191 17150975
Chen et al. (2009) Chen X Li G Feng B Tian S Habitat suitability index of chub mackerel (Scomber japonicus) from July to September in the East China Sea Journal of Oceanography 2009 65 1 93 102 10.1007/s10872-009-0009-9
Chen et al. (2017) Chen FF Ma CY Yan LP Zhang H Wang W Zhang Y Ma LB Isolation and characterization of polymorphic microsatellite markers for the chub mackerel (Scomber japonicus) and cross-species amplification in the blue mackerel (S. australasicus) Genetics and Molecular Research 2017 16 3 gmr16039712 10.4238/gmr16039712
Cheng et al. (2015) Cheng J Yanagimoto T Song N Gao TX Population genetic structure of chub mackerel Scomber japonicus in the Northwestern Pacific inferred from microsatellite analysis Molecular Biology Reports 2015 42 2 373 382 10.1007/s11033-014-3777-2 25366174
Cheng, Zhu & Chen (2014) Cheng Q Zhu Y Chen X High polymorphism and moderate differentiation of chub mackerel, Scomber japonicus (Perciformes: Scombridae), along the coast of China revealed by fifteen novel microsatellite markers Conservation Genetics 2014 15 5 1021 1035 10.1007/s10592-014-0596-x
Collette (1986) Collette BB Whitehead PJP Bauchot ML Hureau JC Nielsen J Tortonese E Scombridae Fishes of the North-Eastern Atlantic and the Mediterranean 1986 Paris, France United Nations Educational Scientific and Cultural Organization 981 997
Collette (1999) Collette B Séret B Sire JY Mackerels, molecules, and morphology Proceedings of the 5th Indo-Pacific Fish Conference 1999 Nouméa Soc. Fr. Ichtyol 149 164
Collette & Nauen (1983) Collette BB Nauen CE Scombrids of the world: an annotated and illustrated catalogue of tunas, mackerels, bonitos, and related species known to date 1983 Roma Food and Agriculture Organization of the United Nations
Costa et al. (2011) Costa G Cavallero S D’Amelio S Paggi L Santamaría M Perera C Santos M Khadem M Helminth parasites of the Atlantic chub mackerel, Scomber colias Gmelin, 1789 from Canary Islands, Central North Atlantic, with comments on their relations with other Atlantic regions Acta Parasitologica 2011 56 1 98 104 10.1017/S0022149X17000104
Dahle et al. (2018a) Dahle G Johansen T Westgaard J-I Aglen A Glover KA Genetic management of mixed-stock fisheries real-time: the case of the largest remaining cod fishery operating in the Atlantic in 2007–2017 Fisheries Research 2018a 205 6 77 85 10.1016/j.fishres.2018.04.006
Dahle et al. (2018b) Dahle G Quintela M Johansen T Westgaard J-I Besnier F Aglen A Jørstad KE Glover KA Analysis of coastal cod (Gadus morhua L.) sampled on spawning sites reveals a genetic gradient throughout Norway’s coastline BMC Genetics 2018b 19 1 42 10.1186/s12863-018-0625-8 29986643
Daley & Leaf (2019) Daley TT Leaf RT Age and growth of Atlantic chub mackerel (Scomber colias) in the Northwest Atlantic Journal of Northwest Atlantic Fishery Science 2019 50 1 12 10.2960/J.v50.m717
De Bie et al. (2012) De Bie T De Meester L Brendonck L Martens K Goddeeris B Ercken D Hampel H Denys L Vanhecke L Van der Gucht K Van Wichelen J Vyverman W Declerck SAJ Body size and dispersal mode as key traits determining metacommunity structure of aquatic organisms Ecology Letters 2012 15 7 740 747 10.1111/j.1461-0248.2012.01794.x 22583795
Dempster, Laird & Rubin (1977) Dempster AP Laird NM Rubin DB Maximum likelihood from incomplete data via the EM algorithm Journal of the Royal Statistical Society Series B (Methodological) 1977 39 1 38 10.1111/J.2517-6161.1977.TB01600.X
Dieringer & Schlötterer (2003) Dieringer D Schlötterer C Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data sets Molecular Ecology Notes 2003 3 1 167 169 10.1046/j.1471-8286.2003.00351.x
Domínguez-Petit et al. (2022) Domínguez-Petit R Navarro MR Cousido-Rocha M Tornero J Ramos F Jurado-Ruzafa A Nunes C Hernández C Silva AV Landa J Spatial variability of life-history parameters of the Atlantic chub mackerel (Scomber colias), an expanding species in the northeast Atlantic Scientia Marina 2022 86 4 e048 10.3989/scimar.05296.048
Ellegren (2004) Ellegren H Microsatellites: simple sequences with complex evolution Nature Reviews Genetics 2004 5 435 445 10.1038/nrg1348
Erguden et al. (2009) Erguden D Oztürk B Erdoğan Z Turan C Morphologic structuring between populations of chub mackerel Scomber japonicus in the Black, Marmara, Aegean, and northeastern Mediterranean Seas Fisheries Science 2009 75 1 129 135 10.1007/s12562-008-0032-6
Evanno, Regnaut & Goudet (2005) Evanno G Regnaut S Goudet J Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study Molecular Ecology 2005 14 2611 2620 10.1111/j.1365-294X.2005.02553.x 15969739
Excoffier, Laval & Schneider (2005) Excoffier L Laval G Schneider S Arlequin ver. 3.0: an integrated software package for population genetics data analysis Evolutionary Bioinformatics Online 2005 1 47 50 10.1177/117693430500100003
FAO (2020) FAO Report of the working group on the assessment of small pelagic fish of Northwest Africa Casablanca, Morocco, 8–13 July 2019. Rapport de groupe de travail sur l’évaluation des petits pêlagiques au large de l’Afrique Nord-Occidentale Casablanca, Maroc, 8-13 juillet 2019. Fishery Committee for the Eastern Central Atlantic (CECAF)/Comité des pêches pour l’Atlantique Centre-Est (COPACE). RFAO Fisheries and Aquaculture Report No. 1309/FAO, Rapport sur les pêches et l’aquaculture no 1309. Rome 2020
Foll & Gaggiotti (2006) Foll M Gaggiotti O Identifying the environmental factors that determine the genetic structure of populations Genetics 2006 174 2 875 891 10.1534/genetics.106.059451 16951078
Foll & Gaggiotti (2008) 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 2 977 993 10.1534/genetics.108.092221 18780740
García (1982) García S Distribution, migration and spawning of the main 1982 fish resources in the northern CECAF area FAO, CECAF/ECAF Series 82/(E) 1982 9
Gissi et al. (2021) Gissi E Manea E Mazaris AD Fraschetti S Almpanidou V Bevilacqua S Coll M Guarnieri G Lloret-Lloret E Pascual M Petza D Rilov G Schonwald M Stelzenmüller V Katsanevakis S A review of the combined effects of climate change and other local human stressors on the marine environment Science of the Total Environment 2021 755 3 142564 10.1016/j.scitotenv.2020.142564 33035971
Glover (2010) Glover KA Forensic identification of fish farm escapees: the Norwegian experience Aquaculture Environment Interactions 2010 1 1 10 10.3354/aei00002
Glover et al. (2011) Glover KA Skaala Ø Limborg M Kvamme C Torstensen E Microsatellite DNA reveals population genetic differentiation among sprat (Sprattus sprattus) sampled throughout the Northeast Atlantic, including Norwegian fjords ICES Journal of Marine Science 2011 68 10 2145 2151 10.1093/icesjms/fsr153
Hale, Burg & Steeves (2012) Hale ML Burg TM Steeves TE Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies PLOS ONE 2012 7 9 e45170 10.1371/journal.pone.0045170 22984627
Hauser et al. (2002) Hauser L Adcock GJ Smith PJ Bernal Ramírez JH Carvalho GR Loss of microsatellite diversity and low effective population size in an overexploited population of New Zealand snapper (Pagrus auratus) Proceedings of the National Academy of Sciences 2002 99 18 11742 11747 10.1073/pnas.172242899
Hauser & Carvalho (2008) Hauser L Carvalho GR Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts Fish and Fisheries 2008 9 4 333 362 10.1111/j.1467-2979.2008.00299.x
Hemmer-Hansen et al. (2019) Hemmer-Hansen J Hüssy K Baktoft H Huwer B Bekkevold D Haslob H Herrmann J-P Hinrichsen H-H Krumme U Mosegaard H Nielsen EE Reusch TBH Storr-Paulsen M Velasco A von Dewitz B Dierking J Eero M Genetic analyses reveal complex dynamics within a marine fish management area Evolutionary Applications 2019 12 830 844 10.1111/eva.1276010.1111/eva.12760 30976313
Henriques et al. (2017) Henriques R Nielsen ES Durholtz D Japp D von der Heyden S Genetic population sub-structuring of kingklip (Genypterus capensis—Ophidiidiae), a commercially exploited demersal fish off South Africa Fisheries Research 2017 187 86 95 10.1016/j.fishres.2016.11.007
Hernández-León, Gómez & Arístegui (2007) Hernández-León S Gómez M Arístegui J Mesozooplankton in the canary current system: the coastal-ocean transition zone Progress in Oceanography 2007 74 2–3 397 421 10.1016/j.pocean.2007.04.010
ICES (2021) ICES Second workshop on atlantic chub mackerel (Scomber colias) ICES Scientific Reports 2021 3 236 10.17895/ices.pub.8142
IPCC (2022) IPCC Pörtner HO Roberts DC Tignor MMB Poloczanska ES Mintenbeck K Alegría A Craig M Langsdorf S Löschke S Möller V Okem A Rama B Climate change 2022: impacts, adaptation and vulnerability Working Group II Contribution to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change 2022 UK and New York, USA Cambridge 3056
Jakobsson & Rosenberg (2007) Jakobsson M Rosenberg NA CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure Bioinformatics 2007 23 14 1801 1806 10.1093/bioinformatics/btm233 17485429
Johansen et al. (2018) Johansen T Westgaard J-I Seliussen BB Nedreaas K Dahle G Glover KA Kvalsund R Aglen A Real-time genetic monitoring of a commercial fishery on the doorstep of an MPA reveals unique insights into the interaction between coastal and migratory forms of the Atlantic cod ICES Journal of Marine Science 2018 75 3 1093 1104 10.1093/icesjms/fsx224
Jombart (2008) Jombart T adegenet: a R package for the multivariate analysis of genetic markers Bioinformatics 2008 24 11 1403 1405 10.1093/bioinformatics/btn129 18397895
Jombart & Collins (2015) Jombart T Collins C A tutorial for discriminant analysis of principal components (DAPC) using adegenet 2.0.0 2015 https://adegenet.r-forge.r-project.org/files/tutorial-dapc.pdf
Jombart, Devillard & Balloux (2010) Jombart T Devillard S Balloux F Discriminant analysis of principal components: a new method for the analysis of genetically structured populations BMC Genetics 2010 11 1 94 10.1186/1471-2156-11-94 20950446
Jurado-Ruzafa et al. (2024) Jurado-Ruzafa A Vasconcelos J Otero-Ferrer JL Navarro MR Massaro A Hernández C Tuset VM Phenotypic response of a geographically expanding species, Scomber colias: clues in the fish otolith shape Estuarine, Coastal and Shelf Science 2024 305 7 108880 10.1016/j.ecss.2024.108880
Kamvar, Tabima & Grünwald (2014) Kamvar ZN Tabima JF Grünwald NJ Poppr: an R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction PeerJ 2014 2 e281 10.7717/peerj.281 24688859
Kerr et al. (2017) Kerr LA Hintzen NT Cadrin SX Clausen LW Dickey-Collas M Goethel DR Hatfield EMC Kritzer JP Nash RDM Lessons learned from practical approaches to reconcile mismatches between biological population structure and stock units of marine fish ICES Journal of Marine Science 2017 74 6 1708 1722 10.1093/icesjms/fsw188
Lamichhaney et al. (2017) Lamichhaney S Fuentes-Pardo AP Rafati N Ryman N McCracken GR Bourne C Singh R Ruzzante DE Andersson L Parallel adaptive evolution of geographically distant herring populations on both sides of the North Atlantic Ocean Proceedings of the National Academy of Sciences 2017 114 E3452 10.1073/pnas.1617728114
Landa et al. (2022) Landa J Domínguez-Petit R Navarro R Silva AV Tornero J Parra C Ramos F Nunes C Jurado-Ruzafa A Growth variability and demographic structure of Northeast Atlantic chub mackerel (Scomber colias) in southern European Atlantic waters 2022 ISMS Las Palmas de Gran Canaria, 6–8 July
Le Moan, Bekkevold & Hemmer-Hansen (2021) Le Moan A Bekkevold D Hemmer-Hansen J Evolution at two time frames: ancient structural variants involved in post-glacial divergence of the European plaice (Pleuronectes platessa) Heredity 2021 126 6 668 683 10.1038/s41437-021-00432-x 33531657
Leone et al. (2019) Leone A Álvarez P García D Saborido-Rey F Rodríguez-Ezpeleta N Genome-wide SNP based population structure in European hake reveals the need for harmonizing biological and management units ICES Journal of Marine Science 2019 76 7 2260 2266 10.1093/icesjms/fsz161
Li & Liu (2018) Li YL Liu JX StructureSelector: a web-based software to select and visualize the optimal number of clusters using multiple methods Molecular Ecology Resources 2018 18 1 176 177 10.1111/1755-0998.12719 28921901
Machado et al. (2022) Machado AM Gomes-Dos-Santos A Fonseca MM da Fonseca RR Veríssimo A Felício M Capela R Alves N Santos M Salvador-Caramelo F Domingues M Ruivo R Froufe E Castro LFC A genome assembly of the Atlantic chub mackerel (Scomber colias): a valuable teleost fishing resource GigaByte 2022 2022 3 gigabyte40 10.46471/gigabyte.40 36824513
Maduna et al. (2014) Maduna SN Rossouw C Roodt-Wilding R Bester-van der Merwe AE Microsatellite cross-species amplification and utility in southern African elasmobranchs: a valuable resource for fisheries management and conservation BMC Research Notes 2014 7 1 352 10.1186/1756-0500-7-352 24915745
Mantel (1967) Mantel N The detection of disease of clustering and a generalized regression approach Cancer Research 1967 27 209 220 6018555
Martinsohn et al. (2019) Martinsohn JT Raymond P Knott T Glover KA Nielsen EE Eriksen LB Ogden R Casey J Guillen J DNA-analysis to monitor fisheries and aquaculture: too costly? Fish and Fisheries 2019 20 2 391 401 10.1111/faf.12343
Martínez et al. (2006) Martínez P González EG Castilho R Zardoya R Genetic diversity and historical demography of Atlantic bigeye tuna (Thunnus obesus) Molecular Phylogenetics and Evolution 2006 39 2 404 416 10.1016/j.ympev.2005.07.022 16188460
Medina-Alcaraz (2016) Medina-Alcaraz C Diversidad genética y estructura poblacional de caballa (Scomber colias, Gmelin, 1789) en aguas del Atlántico y del Mediterráneo PhD thesis. Universidad de Las Palmas de Gran Canaria 2016
Miller, Cullingham & Peery (2020) Miller JM Cullingham CI Peery RM The influence of a priori grouping on inference of genetic clusters: simulation study and literature review of the DAPC method Heredity 2020 125 5 269 280 10.1038/s41437-020-0348-2 32753664
Mittelstaedt (1991) Mittelstaedt E The ocean boundary along the northwest African coast: circulation and oceanographic properties at the sea surface Progress in Oceanography 1991 26 4 307 355 10.1016/0079-6611(91)90011-A
Mullins et al. (2018) Mullins RB McKeown NJ Sauer WHH Shaw PW Genomic analysis reveals multiple mismatches between biological and management units in yellowfin tuna (Thunnus albacares) ICES Journal of Marine Science 2018 75 6 2145 2152 10.1093/icesjms/fsy102
Munro, Van Houtte & Willmann (2004) Munro G Van Houtte A Willmann R The conservation and management of shared fish stocks: legal and economic aspects 2004 In: FAO, editor. Rome: FAO Fisheries Technical Paper, 69
Nesbø et al. (2000) Nesbø CL Rueness EK Iversen SA Skagen DW Jakobsen KS Phylogeography and population history of Atlantic mackerel (Scomber scombrus L.): a genealogical approach reveals genetic structuring among the eastern Atlantic stocks Proceedings Biological Sciences 2000 267 1440 281 292 10.1098/rspb.2000.0998 10714883
Nielsen et al. (2018) Nielsen ES Henriques R Toonen RJ Knapp ISS Guo B von der Heyden S Complex signatures of genomic variation of two non-model marine species in a homogeneous environment BMC Genomics 2018 19 1 347 10.1186/s12864-018-4721-y 29743012
O’Connell & Wright (1997) O’Connell M Wright JM Microsatellite DNA in fishes Reviews in Fish Biology and Fisheries 1997 7 331 363 10.1023/A:1018443912945
Pante & Simon-Bouhet (2013) Pante E Simon-Bouhet B marmap: a package for importing, plotting and analyzing bathymetric and topographic data in R PLOS ONE 2013 8 9 e73051 10.1371/journal.pone.0073051 24019892
Peakall & Smouse (2006) Peakall R Smouse PE GenAlEx 6: genetic analysis in excel. Population genetic software for teaching and research Molecular Ecology Notes 2006 6 1 288 295 10.1111/j.1471-8286.2005.01155.x
Pelegrí & Peña-Izquierdo (2015) Pelegrí JL Peña-Izquierdo J Eastern boundary currents off North-West Africa. Oceanographic and biological features in the Canary Current Large Marine Ecosystem Technical Series: 115 2015 Paris, France IOC-UNESCO
Perrotta, Carvalho & Isidro (2005) Perrotta RG Carvalho N Isidro E Comparative study on growth of chub mackerel (Scomber japonicus Houttuyn, 1782) from three different regions: NW Mediterranean, NE and SW Atlantic Revista de Investigación y Desarrollo Pesquero 2005 17 67 79
Petrou et al. (2021) Petrou EL Fuentes-Pardo AP Rogers LA Orobko M Tarpey C Jiménez-Hidalgo I Moss ML Yang D Pitcher TJ Sandell T Lowry D Ruzzante DE Hauser L Functional genetic diversity in an exploited marine species and its relevance to fisheries management Proceedings of the Royal Society B: Biological Sciences 2021 288 1945 20202398 10.1098/rspb.2020.2398
Pritchard, Stephens & Donnelly (2000) Pritchard JK Stephens M Donnelly P Inference of population structure using multilocus genotype data Genetics 2000 155 2 945 959 10.1093/genetics/155.2.945 10835412
Puechmaille (2016) Puechmaille SJ The program structure does not reliably recover the correct population structure when sampling is uneven: subsampling and new estimators alleviate the problem Molecular Ecology Resources 2016 16 3 608 627 10.1111/1755-0998.12512 26856252
Quintela et al. (2020) Quintela M Kvamme C Bekkevold D Nash RDM Jansson E Sørvik AG Taggart JB Skaala Ø Dahle G Glover KA Genetic analysis redraws the management boundaries for the European sprat Evolutionary Applications 2020 13 8 1906 1922 10.1111/eva.12942 32908594
Ramírez et al. (2022) Ramírez F Shannon LJ van der Lingen CD Julià L Steenbeek J Coll M Climate and fishing simultaneously impact small pelagic fish in the oceans around the southernmost tip of Africa Frontiers in Marine Science 2022 9 45 10.3389/fmars.2022.1031784
Ramya & Behera (2023) Ramya VL Behera BK Behera BK Molecular markers and their application in fisheries and aquaculture Biotechnological Tools in Fisheries and Aquatic Health Management 2023 Singapore Springer Nature Singapore 115 150
R Core Team (2020) R Core Team R: a language and environment for statistical computing 2020 Vienna R Foundation for Statistical Computing http://www.R-project.org/
Reid et al. (2016) Reid K Hoareau TB Graves JE Potts WM dos Santos SMR Klopper AW Bloomer P Secondary contact and asymmetrical gene flow in a cosmopolitan marine fish across the Benguela upwelling zone Heredity 2016 117 5 307 315 10.1038/hdy.2016.51 27436525
Rodríguez-Ezpeleta et al. (2016) Rodríguez-Ezpeleta N Bradbury IR Mendibil I Álvarez P Cotano U Irigoien X Population structure of Atlantic mackerel inferred from RAD-seq-derived SNP markers: effects of sequence clustering parameters and hierarchical SNP selection Molecular Ecology Resources 2016 16 4 991 1001 10.1111/1755-0998.12518 26936210
Rodríguez-Ezpeleta et al. (2019) Rodríguez-Ezpeleta N Díaz-Arce N Walter Iii JF Richardson DE Rooker JR Nøttestad L Hanke AR Franks JS Deguara S Lauretta MV Addis P Varela JL Fraile I Goñi N Abid N Alemany F Oray IK Quattro JM Sow FN Itoh T Karakulak FS Pascual-Alayón PJ Santos MN Tsukahara Y Lutcavage M Fromentin J-M Arrizabalaga H Determining natal origin for improved management of Atlantic bluefin tuna Frontiers in Ecology and the Environment 2019 17 8 439 444 10.1002/fee.2090
Roldán et al. (2000) Roldán MI Perrotta RG Cortey M Pla C Molecular and morphologic approaches to discrimination of variability patterns in chub mackerel, Scomber japonicus Journal of Experimental Marine Biology and Ecology 2000 253 1 63 74 10.1016/S0022-0981(00)00244-6 11018237
Rosenberg & Anderson (2011) Rosenberg MS Anderson CD PASSaGE: pattern analysis, spatial statistics and geographic exegesis. Version 2 Methods in Ecology and Evolution 2011 2 3 229 232 10.1111/j.2041-210X.2010.00081.x
Rousset (1997) Rousset F Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance Genetics 1997 145 1219 1228 10.1093/genetics/145.4.1219 9093870
Rousset (2008) Rousset F GENEPOP′007: a complete re-implementation of the genepop software for windows and Linux Molecular Ecology Resources 2008 8 1 103 106 10.1111/j.1471-8286.2007.01931.x 21585727
Roy et al. (2014) Roy D Hardie DC Treble MA Reist JD Ruzzante DE Evidence supporting panmixia in Greenland halibut (Reinhardtius hippoglossoides) in the Northwest Atlantic Canadian Journal of Fisheries and Aquatic Sciences 2014 71 5 763 774 10.1139/cjfas-2014-0004
Russello et al. (2012) Russello MA Kirk SL Frazer KK Askey PJ Detection of outlier loci and their utility for fisheries management Evolutionary Applications 2012 5 1 39 52 10.1111/j.1752-4571.2011.00206.x 25568028
Ryman & Palm (2006) Ryman N Palm S POWSIM: a computer program for assessing statistical power when testing for genetic differentiation Molecular Ecology Notes 2006 6 3 600 602 10.1111/j.1471-8286.2006.01378.x
Saha et al. (2017) Saha A Johansen T Hedeholm R Nielsen EE Westgaard J-I Hauser L Planque B Cadrin SX Boje J Geographic extent of introgression in Sebastes mentella and its effect on genetic population structure Evolutionary Applications 2017 10 1 77 90 10.1111/eva.12429 28035237
Sbiba et al. (2024) Sbiba SE Correia AT Nikolioudakis N Bazairi H Chlaida M Insights into the stock structure of Scomber colias Gmelin, 1789 along the Northwest African coast inferred from otolith shape analysis Fisheries Research 2024 272 106955 10.1016/j.fishres.2024.106955
Scoles, Graves & Collette (1988) Scoles D Graves J Collette B Global phylogeography of mackerels of the genus Scomber Fishery Bulletin- National Oceanic and Atmospheric Administration 1988 96 823 842
Shannon (2001) Shannon LV Steele JH Benguela current Encyclopedia of Ocean Sciences 2001 Oxford Academic Press 255 267
Slatkin (1993) Slatkin M Isolation by distance in equilibrium and non-equilibrium populations Evolution Education and Outreach 1993 47 1 264 279 10.1111/j.1558-5646.1993.tb01215.x
Smith (1994) Smith PJ Genetic diversity of marine fisheries resources: possible impacts of fishing FAO Fisheries Technical Paper 1994 Rome FAO 53
Stroganov et al. (2023) Stroganov AN Nikitenko AI Rakitskaya TA Belyaev VA Study of population genetic structure of Atlantic mackerel Scomber colias Gmelin, 1789 (Scombridae) in the Central East Atlantic Russian Journal of Genetics 2023 59 1 66 72 10.1134/S1022795423010118
Tang et al. (2009) Tang CY Tzeng CH Chen CS Chiu TS Microsatellite DNA markers for population-genetic studies of blue mackerel (Scomber australasicus) and cross-specific amplification in S. japonicus Molecular Ecology Resources 2009 9 3 824 827 10.1111/j.1755-0998.2008.02278.x 21564757
Tzeng et al. (2009) Tzeng CH Chen CS Tang PC Chiu TS Microsatellite and mitochondrial haplotype differentiation in blue mackerel (Scomber australasicus) from the western North Pacific ICES Journal of Marine Science 2009 66 5 816 825 10.1093/icesjms/fsp120
Uriarte et al. (2001) Uriarte A Alvarez P Iversen SA Molloy J Villamor B Martíns MM Myklevoll S Spatial pattern of migration and recruitment of North East Atlantic mackerel 2001 ICES Annual Science Conference Session O: Application of Mark Recapture Experiments to Stock Assessment
Van Oosterhout et al. (2004) Van Oosterhout C Hutchinson WF Wills DPM Shipley P MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data Molecular Ecology Notes 2004 4 3 535 538 10.1111/j.1471-8286.2004.00684.x
Vásquez et al. (2013) Vásquez S Correa-Ramírez M Parada C Sepúlveda A The influence of oceanographic processes on jack mackerel (Trachurus murphyi) larval distribution and population structure in the southeastern Pacific Ocean ICES Journal of Marine Science 2013 70 6 1097 1107 10.1093/icesjms/fst065
Waldman & McKinnon (1993) Waldman B McKinnon J Thornhill N Inbreeding and outbreeding in fishes, amphibians and reptiles The Natural History of Inbreeding and Outbreeding 1993 Chicago University of Chicago Press 250 282
Waples (1987) Waples RS A multispecies approach to the analysis of gene flow in marine shore fishes Evolution Education and Outreach 1987 41 2 385 400 10.1111/j.1558-5646.1987.tb05805.x
Waples (1998) Waples RS Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species Journal of Heredity 1998 89 5 438 450 10.1093/jhered/89.5.438
Ward (2000) Ward RD Genetics in fisheries management Hydrobiologia 2000 420 191 201 10.1007/978-94-017-2184-4_18
Weber & McClatchie (2011) Weber E McClatchie S Effect of environmental conditions on the distribution of Pacific mackerel (Scomber japonicus) larvae in the California Current system Fishery Bulletin 2011 110 85 97
Weersing & Toonen (2009) Weersing K Toonen RJ Population genetics, larval dispersal, and demographic connectivity in marine systems Marine Ecology Progress Series 2009 393 1 12 10.3354/meps08287
Weir & Cockerham (1984) Weir BS Cockerham C Estimating F-statistics for the analysis of population structure Evolution 1984 38 6 1358 1370 10.2307/2408641 28563791
Westgaard et al. (2017) Westgaard J-I Saha A Kent M Hansen HH Knutsen H Hauser L Cadrin SX Albert OT Johansen T Genetic population structure in Greenland halibut (Reinhardtius hippoglossoides) and its relevance to fishery management Canadian Journal of Fisheries and Aquatic Sciences 2017 74 4 475 485 10.1139/cjfas-2015-0430
White et al. (2010) White C Selkoe KA Watson J Siegel DA Zacherl DC Toonen RJ Ocean currents help explain population genetic structure Proceedings of the Royal Society B: Biological Sciences 2010 277 1688 1685 1694 10.1098/rspb.2009.2214
Wright (1943) Wright S Isolation by distance Genetics 1943 28 114 138 10.1093/genetics/28.2.114 17247074
Zardoya et al. (2004) Zardoya R Castilho R Grande C Favre-Krey L Caetano S Marcato S Krey G Patarnello T Differential population structuring of two closely related fish species, the mackerel (Scomber scombrus) and the chub mackerel (Scomber japonicus), in the Mediterranean Sea Molecular Ecology 2004 13 7 1785 1798 10.1111/j.1365-294X.2004.02198.x 15189203
Zarraonaindia et al. (2009) Zarraonaindia I Pardo MA Iriondo M Manzano C Estonba A Microsatellite variability in European anchovy (Engraulis encrasicolus) calls for further investigation of its genetic structure and biogeography ICES Journal of Marine Science 2009 66 10 2176 2182 10.1093/icesjms/fsp187
Zeng & Cheng (2012) Zeng L Cheng Q Thirty novel microsatellite markers for the coastal pelagic fish, Scomber japonicus (Scombridae) Journal of Genetics 2012 91 S1 64 68 10.1007/s12041-012-0153-7
Zeng, Cheng & Chen (2012) Zeng L Cheng Q Chen X Microsatellite analysis reveals the population structure and migration patterns of Scomber japonicus (Scombridae) with continuous distribution in the East and South China Seas Biochemical Systematics and Ecology 2012 42 83 93 10.1016/j.bse.2012.02.014
