Genetic investigation of population structure in Atlantic chub mackerel, Scomber colias Gmelin, 1789 along the West African coast.

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 (F ST = 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.

Leveraging the genetic diversity of trout in the rivers of the British Isles and northern France to understand the movements of sea trout (Salmo trutta L.) around the English Channel.

Populations of anadromous brown trout, also known as sea trout, have suffered recent marked declines in abundance due to multiple factors, including climate change and human activities. While much is known about their freshwater phase, less is known about the species' marine feeding migrations. This situation is hindering the effective management and conservation of anadromous trout in the marine environment. Using a panel of 95 single nucleotide polymorphism markers we developed a genetic baseline, which demonstrated strong regional structuring of genetic diversity in trout populations around the English Channel and adjacent waters. Extensive baseline testing showed this structuring allowed high-confidence assignment of known-origin individuals to region of origin. This study presents new data on the movements of anadromous trout in the English Channel and southern North Sea. Assignment of anadromous trout sampled from 12 marine and estuarine localities highlighted contrasting results for these areas. The majority of these fisheries are composed predominately of stocks local to the sampling location. However, there were multiple cases of long-distance movements of anadromous trout, with several individuals originating from rivers in northeast England being caught in the English Channel and southern North Sea, in some cases more than 1000 km from their natal region. These results have implications for the management of sea trout in inshore waters around the English Channel and southern North Sea.

Temporal allele frequency changes in large-effect loci reveal potential fishing impacts on salmon life-history diversity.

Fishing has the potential to influence the life-history traits of exploited populations. However, our understanding of how fisheries can induce evolutionary genetic changes remains incomplete. The discovery of large-effect loci linked with ecologically important life-history traits, such as age at maturity in Atlantic salmon (Salmo salar), provides an opportunity to study the impacts of temporally varying fishing pressures on these traits. A 93-year archive of fish scales from wild Atlantic salmon catches from the northern Baltic Sea region allowed us to monitor variation in adaptive genetic diversity linked with age at maturity of wild Atlantic salmon populations. The dataset consisted of samples from both commercial and recreational fisheries that target salmon on their spawning migration. Using a genotyping-by-sequencing approach (GT-seq), we discovered strong within-season allele frequency changes at the vgll3 locus linked with Atlantic salmon age at maturity: fishing in the early season preferentially targeted the vgll3 variant linked with older maturation. We also found within-season temporal variation in catch proportions of different wild Atlantic salmon subpopulations. Therefore, selective pressures of harvesting may vary depending on the seasonal timing of fishing, which has the potential to cause evolutionary changes in key life-history traits and their diversity. This knowledge can be used to guide fisheries management to reduce the effects of fishing practices on salmon life-history diversity. Thus, this study provides a tangible example of using genomic approaches to infer, monitor and help mitigate human impacts on adaptively important genetic variation in nature.

Harnessing the power of regional baselines for broad-scale genetic stock identification: A multistage, integrated, and cost-effective approach.

In mixed-stock fishery analyses, genetic stock identification (GSI) estimates the contribution of each population to a mixture and is typically conducted at a regional scale using genetic baselines specific to the stocks expected in that region. Often these regional baselines cannot be combined to produce broader geographical baselines due to non-overlapping populations and genetic markers. In cases where the mixture contains stocks spanning across a wide area, a broad-scale baseline is created, but often at the cost of resolution. Here, we introduce a new GSI method to harness the resolution capabilities of baselines developed for regional applications in the analysis of mixtures containing individuals from a broad geographic range. This method employs a multistage framework that allows disparate baselines to be used in a single integrated process that produces estimates along with the propagated errors from each stage. All individuals in the mixture sample are required to be genotyped for all genetic markers in the baselines used by this model, but the baselines do not require overlap in genetic markers or populations representing the broad-scale or regional baselines. We demonstrate the utility of our integrated multistage model using a synthesized data set made up of Chinook salmon, Oncorhynchus tshawytscha, from the North Bering Sea of Alaska. The results show an improved accuracy for estimates using an integrated multistage framework, compared to the conventional framework of using separate hierarchical steps. The integrated multistage framework allows GSI of a wide geographic area without first developing a large scale, high-resolution genetic baseline or dividing a mixture sample into smaller regions beforehand. This approach is more cost-effective than updating range-wide baselines with all regionally important markers.

Modeling ocean distributions and abundances of natural- and hatchery-origin Chinook salmon stocks with integrated genetic and tagging data.

Background: Considerable resources are spent to track fish movement in marine environments, often with the intent of estimating behavior, distribution, and abundance. Resulting data from these monitoring efforts, including tagging studies and genetic sampling, often can be siloed. For Pacific salmon in the Northeast Pacific Ocean, predominant data sources for fish monitoring are coded wire tags (CWTs) and genetic stock identification (GSI). Despite their complementary strengths and weaknesses in coverage and information content, the two data streams rarely have been integrated to inform Pacific salmon biology and management. Joint, or integrated, models can combine and contextualize multiple data sources in a single statistical framework to produce more robust estimates of fish populations. Methods: We introduce and fit a comprehensive joint model that integrates data from CWT recoveries and GSI sampling to inform the marine life history of Chinook salmon stocks at spatial and temporal scales relevant to ongoing fisheries management efforts. In a departure from similar models based primarily on CWT recoveries, modeled stocks in the new framework encompass both hatchery- and natural-origin fish. We specifically model the spatial distribution and marine abundance of four distinct stocks with spawning locations in California and southern Oregon, one of which is listed under the U.S. Endangered Species Act. Results: Using the joint model, we generated the most comprehensive estimates of marine distribution to date for all modeled Chinook salmon stocks, including historically data poor and low abundance stocks. Estimated marine distributions from the joint model were broadly similar to estimates from a simpler, CWT-only model but did suggest some differences in distribution in select seasons. Model output also included novel stock-, year-, and season-specific estimates of marine abundance. We observed and partially addressed several challenges in model convergence with the use of supplemental data sources and model constraints; similar difficulties are not unexpected with integrated modeling. We identify several options for improved data collection that could address issues in convergence and increase confidence in model estimates of abundance. We expect these model advances and results provide management-relevant biological insights, with the potential to inform future mixed-stock fisheries management efforts, as well as a foundation for more expansive and comprehensive analyses to follow.

Exploring morphometric frontiers: A comprehensive study of otolith growth patterns in brown comber Serranus hepatus (Linnaeus, 1758).

Otoliths are widely employed in marine sciences to gain insights into fish growth, age, migrations, and population structure. This study investigates the relationships between morphometric measurements, otolith characteristics, and length size patterns in the brown comber (Serranus hepatus) from the Gulf of Cádiz, a species discarded in artisanal trawl fisheries. Our findings reveal significant changes in otolith shape indices as fish grow, with symmetry observed between left and right otolith measurements. Otolith size is found to be related to fish size, supporting its use in estimating body length at different life stages. Otolith shape analysis has potential applications in stock identification, detecting catch misreporting, and studying marine predator diets. Combining otolith shape analysis with other data types can clarify relationships among taxa and inform spatial management strategies, contributing to the long-term sustainability of fish populations and the assessment of the impact of management strategies on fish size and growth. This study enhances our understanding of the broader implications of morphometric and otolith analyses in fisheries research and supports the development of more sustainable fisheries management practices.

Parasites of Sardinella maderensis (Lowe, 1838) (Actinopterygii: Clupeidae) and Their Potential as Biological Tags for Stock Identification along the Coast of West Africa.

This study is the first to provide information on the parasite fauna of Sardinella maderensis along the coasts of Benin and Ghana, and the first to investigate the potential use of parasites as biological tags in fish population studies in the area. It may thus serve as a starting point for upcoming studies. From February to June 2021, a total of 200 S. maderensis were sampled from the fishing port of Cotonou (Benin) and the Elmina landing site (Ghana). The prevalence and abundance of each parasite were recorded. The following are the outcomes of this study: Parasite species, such as Parahemiurus merus, Mazocraeoides sp. and Hysterothylacium fortalezae, were recorded along the coasts of Benin and Ghana, while Anisakis sp(p). and Tentacularia coryphaenae were only recorded along the coast of Benin. Parahemiurus merus was the most prevalent and abundant among all the parasites recorded. Anisakis sp(p). and T. coryphaenae were selected as having potential in the stock identification of S. maderensis. Both parasites were only recorded along the coast of Benin at a low prevalence. As a result, examinations of more S. maderensis from each location for these parasites may justify their use in stock identification studies.

Spatiotemporal Variability of Trace Elements Fingerprints in Otoliths of Japanese Eel (Anguilla japonica) and Its Use in Tracing Geographic Origin.

To secure traceability along supply chains of foodstuffs, the spatiotemporal variability of trace elements' fingerprints (TEF) in fish otoliths provides a powerful tool to determine and discriminate the origin. Spatiotemporal variability of TEF was examined in a commercially important seafood, Japanese eel (Anguilla japonica), by means of laser ablation-inductively coupled plasma mass spectrometry (LA-ICPMS). Six elemental ratios (Na:Ca, Mg:Ca, P:Ca, K:Ca, Sr:Ca, and Ba:Ca) were determined in the otoliths of specimens originating from four aquaculture farms to examine the spatial variability and from one wild habitat over three years to examine the temporal variation. Significant temporal variation was found in Mg:Ca and Sr:Ca ratios; however, discriminant function analysis showed a lower temporal variation (50%) for the three years. Spatial variations were significant in Sr:Ca and Ba:Ca ratios, and discriminant function analysis showed high (80%) spatial variation among the four farms. Otolith TEF in the Japanese eel showed specific spatial variation among aquaculture farms but intangible temporal variation, suggesting the otolith TEF reflect each aquaculture environment. The present study shows that otolith TEF can be a reliable tool to discriminate the geographic origin of the Japanese eel.

A baseline for the genetic stock identification of Atlantic herring, Clupea harengus, in ICES Divisions 6.a, 7.b-c.

Atlantic herring in International Council for Exploration of the Sea (ICES) Divisions 6.a, 7.b-c comprises at least three populations, distinguished by temporal and spatial differences in spawning, which have until recently been managed as two stocks defined by geographical delineators. Outside of spawning the populations form mixed aggregations, which are the subject of acoustic surveys. The inability to distinguish the populations has prevented the development of separate survey indices and separate stock assessments. A panel of 45 single-nucleotide polymorphisms, derived from whole-genome sequencing, were used to genotype 3480 baseline spawning samples (2014-2021). A temporally stable baseline comprising 2316 herring from populations known to inhabit Division 6.a was used to develop a genetic assignment method, with a self-assignment accuracy greater than 90%. The long-term temporal stability of the assignment model was validated by assigning archive (2003-2004) baseline samples (270 individuals) with a high level of accuracy. Assignment of non-baseline samples (1514 individuals) from Divisions 6.a, 7.b-c indicated previously unrecognized levels of mixing of populations outside of the spawning season. The genetic markers and assignment models presented constitute a 'toolbox' that can be used for the assignment of herring caught in mixed survey and commercial catches in Division 6.a into their population of origin with a high level of accuracy.

Exploratory stock identification through morphometric trait analysis of the estuarine round herring, Gilchristella aestuaria (Pisces: Clupeidae).

As a first step in stock identification of estuarine round herring Gilchristella aestuaria, a morphometric analysis of multiple standardised traits and body condition was conducted on samples from sites across the distribution range of the species, including freshwater and estuarine populations. Multivariate analyses, including a principal component analysis (PCA), revealed that sites on the same river system and sites in close geographic proximity tend to share morphometric traits. Most of the variation in PCA was due to caudal fin length (CFL). Single traits were analysed with generalised additive models with river system location as the smooth term. CFL was strongly related to river systems, and to a lesser extent salinity and turbidity. Unlike previous local-scale studies, this broad-scale study did not support the notion of ecotypes based on eye diameter. Condition indices like body depth and relative weight were related to climatic conditions, salinity and turbidity, more than to river system. The findings agree with previous studies on phylogenetic history and limited gene flow in G. aestuaria, and suggest variation in environmental productivity that affects body condition in different populations. Subsequent research should examine temporal changes in traits like CFL and body condition based on long-term seasonal sampling. Safeguarding this potential resource may be enabled by adopting appropriate management methods based on delineation of stocks with different levels of productivity and connectivity, before widespread fishing of this species is promoted.

Genetic assignment of fisheries bycatch reveals disproportionate mortality among Alaska Northern Fulmar breeding colonies.

Global fisheries kill millions of seabirds annually through bycatch, but little is known about population-level impacts, particularly in species that form metapopulations. U.S. North Pacific groundfish fisheries catch thousands of Northern Fulmars (Fulmarus glacialis rodgersii) each year, making fulmars the most frequently caught seabird in federally managed U.S. fisheries. Here, we used genetic stock identification to assign 1,536 fulmars sampled as bycatch to one of four Alaska breeding colonies and quantified the similarity of bycatch locations at sea among colonies. We found disproportionately high bycatch from the Pribilof Islands (6% of metapopulation, 23% of bycatch), and disproportionately low bycatch from Chagulak Island (34% of metapopulation, 14% of bycatch). Overlap between fisheries and colony-specific foraging areas diverge more during the summer breeding season, leading to greater differences in bycatch susceptibility. Contemporary and historical gene flow likely contributes to low genetic differentiation among colonies (FST = 0.003-0.01), yet these values may not represent present connectivity. Our findings illustrate how genetic stock identification can link at-sea threats to colonies and inform management to reduce bycatch from impacted colonies.

In-field genetic stock identification of overwintering coho salmon in the Gulf of Alaska: Evaluation of Nanopore sequencing for remote real-time deployment.

Genetic stock identification (GSI) from genotyping-by-sequencing of single nucleotide polymorphism (SNP) loci has become the gold standard for stock of origin identification in Pacific salmon. The sequencing platforms currently applied require large batch sizes and multiday processing in specialized facilities to perform genotyping by the thousands. However, recent advances in third-generation single-molecule sequencing platforms, such as the Oxford Nanopore minION, provide base calling on portable, pocket-sized sequencers and promise real-time, in-field stock identification of variable batch sizes. Here we evaluate utility and comparability to established GSI platforms of at-sea stock identification of coho salmon (Oncorhynchus kisutch) using targeted SNP amplicon sequencing on the minION platform during a high-sea winter expedition to the Gulf of Alaska. As long read sequencers are not optimized for short amplicons, we concatenate amplicons to increase coverage and throughput. Nanopore sequencing at-sea yielded data sufficient for stock assignment for 50 out of 80 individuals. Nanopore-based SNP calls agreed with Ion Torrent-based genotypes in 83.25%, but assignment of individuals to stock of origin only agreed in 61.5% of individuals, highlighting inherent challenges of Nanopore sequencing, such as resolution of homopolymer tracts and indels. However, poor representation of assayed salmon in the queried baseline data set contributed to poor assignment confidence on both platforms. Future improvements will focus on lowering turnaround time and cost, increasing accuracy and throughput, as well as augmentation of the existing baselines. If successfully implemented, Nanopore sequencing will provide an alternative method to the large-scale laboratory approach by providing mobile small batch genotyping to diverse stakeholders.

Genetic stock identification in Perna viridis (Linnaeus1758) from the Indian Peninsula by using microsatellite markers.

BACKGROUND: Perna viridis (Linnaeus, 1758), the Asian green mussel, is native to the Asia-Pacific region. The species is extensively distributed in the Indian subcontinent and is a candidate species for aquaculture in the Southeast Asian region. Availability of genetic information on wild populations is essential for the effective conservation and management of Perna species. The present study assessed the genetic variation and population structure across the distribution range of this species from the Indian peninsula by using microsatellite markers to determine the genetic structuring among the species. METHODS: A total of 15 microsatellite loci with M13 labeling were used for the genetic characterization of P. viridis along Indian waters. Genotyped data were analyzed using analytical software to determine the genetic stocks and understand the genetic variability across the populations. RESULTS: We identified 15 polymorphic markers to understand the genetic stocks and variability across Perna populations. The mean value of the observed heterozygosity (Hobs: 0.741) for all populations was closer to the expected heterozygosity (Hexp: 0.75). The pairwise Fst values between the west and east coasts of India varied significantly, indicating the existence of significant genetic structure between the populations. CONCLUSIONS: Genetic stock identification using software analysis exhibited two distinct stocks, one along the west coast (Arabian Sea) and another along the east coast (Bay of Bengal). Bottleneck analysis indicated the genetic stability of species in the wild. P. viridis is a commercially vital species in Indian peninsular regions. The present study suggests the adoption of stock-specific relaying programs of the species from Indian waters in future studies.

Russian honey bee genotype identification through enhanced marker panel set.

Russian honey bees (RHB) are a breeding population developed by USDA-ARS as an effort to provide Varroa-resistant honey bees to beekeepers. The selection strategy for this breeding population was the first in honey bees to incorporate genetic stock identification (GSI). The original GSI approach has been in use for over a decade, and though effective, novel technologies and analytical approaches recently developed provide an opportunity for improvement. Here we outline a novel genotyping assay that capitalizes on the markers used in the GSI as well as new loci recently identified in a whole genome pooled study of commercial honey bee stocks. Our approach utilizes a microfluidic platform and machine learning analyses to arrive at an accurate, high throughput assay. This novel approach provides an improved tool that can be readily incorporated into breeding decisions towards healthier more productive bees.

Parentage-based tagging combined with genetic stock identification is a cost-effective and viable replacement for coded-wire tagging in large-scale assessments of marine Chinook salmon fisheries in British Columbia, Canada.

Wild Pacific salmon, including Chinook salmon Oncorhynchus tshawytscha, have been supplemented with hatchery propagation for over 50 years in support of increased ocean harvest, mitigation for hydroelectric development, and conservation of threatened populations. In Canada, the Wild Salmon Policy for Pacific salmon was established with the goal of maintaining and restoring healthy and diverse Pacific salmon populations, making conservation of wild salmon and their habitats the highest priority for resource management decision-making. For policy implementation, a new approach to the assessment and management of Chinook salmon and the associated hatchery production and fisheries management are needed. Implementation of genetic stock identification (GSI) and parentage-based tagging (PBT) for marine fisheries assessment may overcome problems associated with coded-wire tag-based (CWT) assessment and management of Chinook salmon fisheries, providing at a minimum information equivalent to that derived from the CWT program. GSI and PBT were used to identify Chinook salmon sampled in 2018 and 2019 marine fisheries (18,819 individuals genotyped) in British Columbia to specific conservation units (CU), populations, and broodyears. Individuals were genotyped at 391 single nucleotide polymorphisms via direct sequencing of amplicons. Very high accuracy of assignment to population and age (>99.5%) via PBT was observed for 1994 Chinook salmon of ages 2-4 years, with a 105,722-individual, 380-population baseline available for assignment. Application of a GSI-PBT system of identification to individuals in 2019 fisheries provided high-resolution estimates of stock composition, catch, and exploitation rate by CU or population, with fishery exploitation rates directly comparable to those provided by CWTs for 13 populations. GSI and PBT provide an alternate, cheaper, and more effective method in the assessment and management of Canadian-origin Chinook salmon relative to CWTs, and an opportunity for a genetics-based system to replace the current CWT system for salmon assessment.

Mixed-stock analysis using Rapture genotyping to evaluate stock-specific exploitation of a walleye population despite weak genetic structure.

Mixed-stock analyses using genetic markers have informed fisheries management in cases where strong genetic differentiation occurs among local spawning populations, yet many fisheries are supported by multiple, weakly differentiated stocks. Freshwater fisheries exemplify this problem, with many populations supported by multiple stocks of young evolutionary age and isolated across small spatial scales. Consequently, attempts to conduct genetic mixed-stock analyses of inland fisheries have often been unsuccessful. Advances in genomic sequencing offer the ability to discriminate among populations with weak population structure, providing the necessary resolution to conduct mixed-stock assignment among previously indistinguishable stocks. We used genomic data to conduct a mixed-stock analysis of eastern Lake Erie's commercial and recreational walleye (Sander vitreus) fisheries and estimate the relative harvest of weakly differentiated stocks (pairwise F ST < 0.01). Using RAD-capture (Rapture), we sequenced and genotyped individuals from western and eastern basin local spawning stocks at 12,081 loci with 95% reassignment accuracy, which was not possible in the past using microsatellite markers. A baseline assessment of 395 walleye from 11 spawning stocks identified three reporting groups and refined previous assessments of gene flow among walleye stocks. Genetic assignment of 1,075 walleye harvested in eastern Lake Erie's recreational and commercial fisheries indicated that western basin stocks constituted the majority of harvest during the peak walleye fishing season (July-September), whereas eastern basin individuals comprised much of the early season harvest (May-June). Clear spatial structure in harvest composition existed; catches in more easterly sites contained more individuals of eastern basin origin than did more westerly sites. Our study provides important stock contribution estimates for Lake Erie fishery management and demonstrates the utility of genomic data to facilitate mixed-stock analysis in exploited fish populations having weak population structure or limited existing genetic resources.

An integrated model of seasonal changes in stock composition and abundance with an application to Chinook salmon.

Population-specific spatial and temporal distribution data are necessary to identify mechanisms regulating abundance and to manage anthropogenic impacts. However the distributions of highly migratory species are often difficult to resolve, particularly when multiple populations' movements overlap. Here we present an integrated model to estimate spatially-stratified, seasonal trends in abundance and population composition, using data from extensive genetic sampling of commercial and recreational Chinook salmon (Oncorhynchus tshawytscha) fisheries in southern British Columbia. We use the model to estimate seasonal changes in population-specific standardized catch per unit effort (a proxy for abundance) across six marine regions, while accounting for annual variability in sampling effort and uncertain genetic stock assignment. We also share this model as an R package stockseasonr for application to other regions and species. Even at the relatively small spatial scales considered here, we found that patterns in seasonal abundance differed among regions and stocks. While certain locations were clearly migratory corridors, regions within the Salish Sea exhibited diverse, and often weak, seasonal patterns in abundance, emphasizing that they are important, year-round foraging habitats. Furthermore, we found evidence that stocks with similar freshwater life histories and adult run timing, as well as relatively proximate spawning locations, exhibited divergent distributions. Our findings highlight subtle, but important differences in how adult Chinook salmon use marine habitats. Down-scaled model outputs could be used to inform ecosystem-based management efforts by resolving the degree to which salmon overlap with other species of concern, as well as specific fisheries. More broadly, variation in stock-specific abundance among regions indicates efforts to identify mechanisms driving changes in size-at-maturity and natural mortality should account for distinct marine distributions.

Stock differentiation of the greater lizardfish Saurida tumbil (Teleostei: Synodontidae) collected along the western coast of the Arabian Gulf and Sea of Oman using meristic characters.

Meristic variation among stocks of greater lizardfish Saurida tumbil through the western coasts of the Arabian Gulf and Sea of Oman was examined using meristic characters. Statistical analysis of meristic traits proposed that there is constrained migration of populations of greater lizardfish along the western coast of the Arabian Gulf and Sea of Oman. Overlapping of the two samples from the northern part of the Arabian Gulf (Iraq-Kuwait waters), three samples from the middle region of the Arabian Gulf (Bahrain-Qatar-Saudi Arabia) and two samples from the southern part of the Arabian Gulf/Sea of Oman (United Arab Emirates-Sultanate of Oman) suggested that there are three self-recruiting populations in the studied area. Inspection of the role of each meristic trait variable to Canonical discriminant analysis showed that changes among samples appeared to be linked with the pattern of distribution of water temperature and configuration of current in both the Arabian Gulf and Sea of Oman areas.

Discrimination of red porgy Pagrus pagrus (Sparidae) potential stocks in the south-western Atlantic by otolith shape analysis.

Otolith shape analysis is a powerful method for fish stock identification. We compared the otolith shape of Pagrus pagrus (Linnaeus 1758) along with its distribution in four south-western Atlantic regions where it is commercially fished: Rio de Janeiro, Rio Grande do Sul in southern Brazil, the Argentine-Uruguayan Common Fishing Zone (UA) and the Argentinian Exclusive Fishing Zone (AR). Otolith shapes were compared by Elliptical Fourier and Wavelet coefficients among specimens in a size range with similar otoliths, morphometric parameters and ages. Four potential stocks were identified: one in the AR, a second along the UA which included specimens from southern Brazil with well-marked opaque bands in its otoliths (MRS), the third in southern Brazil with faint or absent opaque bands in its otoliths (FRS) and the fourth along Rio de Janeiro. The difference in the otolith shape among regions followed differences reported using other stock identification techniques. The similarity between otoliths from UA and MRS (ANOVA-like, P > 0.01) can be explained by seasonal short-range migrations. Otoliths shape differences between MRS and FRS (ANOVA-like, P < 0.01) suggest that P. pagrus does not form a homogeneous group in southern Brazil.

Estimating the contribution of Greenland Halibut (Reinhardtius hippoglossoides) stocks to nurseries by means of genotyping-by-sequencing: Sex and time matter.

Identification of stocks and quantification of their relative contribution to recruitment are major objectives toward improving the management and conservation of marine exploited species. Next-generation sequencing allows for thousands of genomic markers to be analyzed, which provides the resolution needed to address these questions in marine species with weakly differentiated populations. Greenland Halibut (Reinhardtius hippoglossoides) is one of the most important exploited demersal species throughout the North Atlantic, in particular in the Gulf of St. Lawrence, Canada. There, two nurseries are known, the St. Lawrence Estuary and the northern Anticosti Island, but their contribution to the renewal of stocks remains unknown. The goals of this study were (a) to document the genetic structure and (b) to estimate the contribution of the different identified breeding stocks to nurseries. We sampled 100 juveniles per nursery and 50 adults from seven sites ranging from Saguenay Fjord to offshore Newfoundland, with some sites sampled over two consecutive years in order to evaluate the temporal stability of the contribution. Our results show that after removing sex-linked markers, the Estuary/Gulf of St. Lawrence represents a single stock which is genetically distinct from the Atlantic around Newfoundland (F ST = 0.00146, p-value = .001). Population assignment showed that recruitment in both nurseries is predominantly associated with the St. Lawrence stock. However, we found that the relative contribution of both stocks to the nurseries is temporally variable with 1% contribution of the Newfoundland stock one year but up to 33% for the second year, which may be caused by year-to-year variation in larval transport into the Gulf of St. Lawrence. This study serves as a model for the identification of stocks for fisheries resources in a context where few barriers to dispersal occur, in addition to demonstrating the importance of considering sex-linked markers and temporal replicates in studies of population genomics.