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

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

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Genomic analyses indicate resilience of a commercially and culturally important marine gastropod snail to climate change

Genomic vulnerability analyses are being increasingly used to assess the adaptability of species to climate change and provide an opportunity for proactive management of harvested marine species in changing oceans. Southeastern Australia is a climate change hotspot where many marine species are shifting poleward. The turban snail, Turbo militaris is a commercially and culturally harvested marine gastropod snail from eastern Australia. The species has exhibited a climate-driven poleward range shift over the last two decades presenting an ongoing challenge for sustainable fisheries management. We investigate the impact of future climate change on T. militaris using genotype-by-sequencing to project patterns of gene flow and local adaptation across its range under climate change scenarios. A single admixed, and potentially panmictic, demographic unit was revealed with no evidence of genetic subdivision across the species range. Significant genotype associations with heterogeneous habitat features were observed, including associations with sea surface temperature, ocean currents, and nutrients, indicating possible adaptive genetic differentiation. These findings suggest that standing genetic variation may be available for selection to counter future environmental change, assisted by widespread gene flow, high fecundity and short generation time in this species. We discuss the findings of this study in the content of future fisheries management and conservation.

Genetic monitoring on the world's first MSC eco-labeled common octopus (O. vulgaris) fishery in western Asturias, Spain

Octopus vulgaris (Cuvier, 1797) is a cephalopod species with great economic value. In western Asturias (northwest of Spain), O. vulgaris artisanal fisheries are relatively well monitored and conditionally eco-labeled by the Marine Stewardship Council (MSC). Despite this, the Asturian octopus stocks have not been genetically assessed so far. In order to improve the current fishery plan and contrast the octopus eco-label validity in Asturias, 539 individuals from five regions of the O. vulgaris geographic distribution, including temporal samplings in Asturias, were collected and genotyped at thirteen microsatellite loci. All the samples under analysis were in agreement with Hardy-Weinberg expectations. Spatial levels of genetic differentiation were estimated using F-statistics, multidimensional scaling, and Bayesian analyses. Results suggested that the O. vulgaris consists of at least four genetically different stocks coming from two ancestral lineages. In addition, temporal analyses showed stability in terms of genetic variation and high N_E (> 50) for several generations in different localities within Asturias, pointing out to indeed sustainable fishery exploitation levels. Even though, the current Asturias fishery plan shows no significant genetic damages to the stocks, the regional-specific management plans need systematic genetic monitoring schemes as part of an efficient and preventive regional fishery regulation strategy.

A Window of Vulnerability: Chronic Environmental Stress Does Not Impair Reproduction in the Swordfish Xiphias gladius

Xiphias gladius is an important fishing resource. The Mediterranean stock is affected by overfishing and is declining. In this light, the aim of this study was to evaluate the cross-talk among metabolism, stress response, immune system and reproduction in immature and mature females, coupling histological and transcriptomic approaches. The transcriptome of livers from 3 immature and 3 mature females was analyzed using the Artificial Intelligence RNA-Seq. For the histological analysis, ovary and liver samples were collected from 50 specimens caught during the reproductive season in the Mediterranean Sea. A total of 750 genes were differentially expressed between the livers. The gene ontologtabey analysis showed 91 upregulated and 161 downregulated biological process GO terms. Instead, the KEGG enrichment analysis revealed 15 enriched pathways. Furthermore, the binding occurring between estrogen receptors and aryl hydrocarbon receptor nuclear translocator, upregulated in mature females, could be liable for the inhibition of detoxification pathway. Indeed, at the histological level, mature females showed a higher density and number of melanomacrophage centers, biomarkers of stress. The present findings reveal the cross-talk among response to environmental stressors, metabolism and reproduction, highlighting that mature females invest a lot of energy in reproduction instead of immune response and detoxification.

Three mitochondrial lineages and no Atlantic-Mediterranean barrier for the bogue Boops boops across its widespread distribution

Marine species exhibiting wide distributional ranges are frequently subdivided into discrete genetic units over limited spatial scales. This is often due to specific life-history traits or oceanographic barriers that prevent gene flow. Fine-scale sampling studies revealed distinct phylogeographic patterns in the northeastern Atlantic and the Mediterranean, ranging from panmixia to noticeable population genetic structure. Here, we used mitochondrial sequence data to analyse connectivity in the bogue Boops boops throughout most of its widespread distribution. Our results identified the existence of three clades, one comprising specimens from the Azores and eastern Atlantic/Mediterranean, another with individuals from the Canary Islands, Madeira and Cape Verde archipelagos, and the third with samples from Mauritania only. One of the branches of the northern subtropical gyre (Azores Current) that drifts towards the Gulf of Cádiz promotes a closer connection between the Azores, southern Portugal and the Mediterranean B. boops populations. The Almería-Oran Front, widely recognised as an oceanographic barrier for many organisms to cross the Atlantic-Mediterranean divide, does not seem to affect the dispersal of this benthopelagic species. The southward movement of the Cape Verde Frontal Zone during the winter, combined with the relatively short duration of the pelagic larval stage of B. boops, may be potential factors for preventing the connectivity between the Atlantic oceanic archipelagos and Mauritania shaping the genetic signature of this species.

A new GTSeq resource to facilitate multijurisdictional research and management of walleye Sander vitreus

Conservation and management professionals often work across jurisdictional boundaries to identify broad ecological patterns. These collaborations help to protect populations whose distributions span political borders. One common limitation to multijurisdictional collaboration is consistency in data recording and reporting. This limitation can impact genetic research, which relies on data about specific markers in an organism's genome. Incomplete overlap of markers between separate studies can prevent direct comparisons of results. Standardized marker panels can reduce the impact of this issue and provide a common starting place for new research. Genotyping-in-thousands (GTSeq) is one approach used to create standardized marker panels for nonmodel organisms. Here, we describe the development, optimization, and early assessments of a new GTSeq panel for use with walleye (Sander vitreus) from the Great Lakes region of North America. High genome-coverage sequencing conducted using RAD capture provided genotypes for thousands of single nucleotide polymorphisms (SNPs). From these markers, SNP and microhaplotype markers were chosen, which were informative for genetic stock identification (GSI) and kinship analysis. The final GTSeq panel contained 500 markers, including 197 microhaplotypes and 303 SNPs. Leave-one-out GSI simulations indicated that GSI accuracy should be greater than 80% in most jurisdictions. The false-positive rates of parent-offspring and full-sibling kinship identification were found to be low. Finally, genotypes could be consistently scored among separate sequencing runs >94% of the time. Results indicate that the GTSeq panel that we developed should perform well for multijurisdictional walleye research throughout the Great Lakes region.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Advancing the protection of marine life through genomics

The rapid growth in genomic techniques provides the potential to transform how we protect, manage, and conserve marine life. Further, solutions to boost the resilience of marine species to climate change and other disturbances that characterize the Anthropocene require transformative approaches, made more effective if guided by genomic data. Although genetic techniques have been employed in marine conservation for decades and the availability of genomic data is rapidly expanding, widespread application still lags behind other data types. This Essay reviews how genetics and genomics have been utilized in management initiatives for ocean conservation and restoration, highlights success stories, and presents a pathway forward to enhance the uptake of genomic data for protecting our oceans.

Population structure of lake whitefish (Coregonus clupeaformis) from the Mississippian lineage in North America

The lake whitefish ( Coregonus clupeaformis) is a commercially valuable freshwater species with a broad distribution in North America. Some phylogeographic work has been done on this species, but little is known about genetic population subdivision among populations of the widely dispersed Mississippian lineage. We used 3,173 single nucleotide polymorphisms in 508 lake whitefish from 22 different lakes to examine population structure across central Canada and the United States. Bayesian clustering, ordination, and fixation indices identified population subdivision that largely reflected geographic distance and hydrological connectivity, with greater differentiation between lakes that are farther apart. Population subdivision was hierarchical, with greater differentiation between Canadian provinces and less differentiation based on river basins within provincial boundaries. Interestingly, isolation by distance alone was not sufficient to account for all of the observed genetic differentiation among populations. We conclude that important components of lake whitefish genetic diversity are present at different spatial scales, and that populations within the Mississippian lineage have differentiated widely across their range.

Genotyping-in-Thousands by sequencing of archival fish scales reveals maintenance of genetic variation following a severe demographic contraction in kokanee salmon

Historical DNA analysis of archival samples has added new dimensions to population genetic studies, enabling spatiotemporal approaches for reconstructing population history and informing conservation management. Here we tested the efficacy of Genotyping-in-Thousands by sequencing (GT-seq) for collecting targeted single nucleotide polymorphism genotypic data from archival scale samples, and applied this approach to a study of kokanee salmon (Oncorhynchus nerka) in Kluane National Park and Reserve (KNPR; Yukon, Canada) that underwent a severe 12-year population decline followed by a rapid rebound. We genotyped archival scales sampled pre-crash and contemporary fin clips collected post-crash, revealing high coverage (> 90% average genotyping across all individuals) and low genotyping error (< 0.01% within-libraries, 0.60% among-libraries) despite the relatively poor quality of recovered DNA. We observed slight decreases in expected heterozygosity, allelic diversity, and effective population size post-crash, but none were significant, suggesting genetic diversity was retained despite the severe demographic contraction. Genotypic data also revealed the genetic distinctiveness of a now extirpated population just outside of KNPR, revealing biodiversity loss at the northern edge of the species distribution. More broadly, we demonstrated GT-seq as a valuable tool for collecting genome-wide data from archival samples to address basic questions in ecology and evolution, and inform applied research in wildlife conservation and fisheries management.

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

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

Long-lived marine species may be resilient to environmental variability through a temporal portfolio effect

Maintenance of genetic variation may provide resilience of populations to natural environmental variability. We used Pacific ocean perch (POP; Sebastes alutus) to test for the maintenance of adaptive variation across overlapping generations. POP are a long-lived species characterized by widespread larval dispersal in their first year and a longevity of over 100 years. In order to understand how early marine dispersal affects POP survival and population structure, we used restriction site-associated DNA sequencing (RADseq) to obtain 11,146 single-nucleotide polymorphisms (SNPs) from 401 young-of-the-year (YOY) POP collected during surveys conducted in 2014 (19 stations) and 2015 (4 stations) in the eastern Gulf of Alaska. Population clustering analysis showed that the POP samples represented four distinct ancestral populations mixed throughout the sampling area. Based on prior work on larval dispersal of POP, these larvae are most likely from distinct parturition locations that are mixing during their pelagic dispersal life stage. Latent factor mixed models revealed that POP larvae face significant selection during their first year at sea, which is specific to the year of their birth. Thus each adult cohort's genetic composition is heavily influenced by the environmental conditions experienced during their first year at sea. Long-lived species relying on broadcast spawning strategies may therefore be uniquely resilient to environmental variability by maintaining a portfolio of cohort-specific adaptive genotypes, and age truncation due to overfishing of older cohorts may have detrimental effect on the population viability.

Sex matters: Otolith shape and genomic variation in deacon rockfish (Sebastes diaconus)

Little is known about intraspecific variation within the deacon rockfish (Sebastes diaconus), a recently described species found in the northeast Pacific Ocean. We investigated population structure among fish sampled from two nearshore reefs (Siletz Reef and Seal Rock) and one offshore site (Stonewall Bank) within a <50-km2 area off the Oregon coast. Fish from the three sample sites exhibited small but statistically significant differences based on genetic variation at >15,000 neutral loci, whether analyzed independently or classified into nearshore and offshore groups. Male and females were readily distinguished using genetic data and 92 outlier loci were associated with sex, potentially indicating differential selection between males and females. Morphometric results indicated that there was significant secondary sexual dimorphism in otolith shape, but further sampling is required to disentangle potential confounding influence of age. This study is the first step toward understanding intraspecific variation within the deacon rockfish and the potential management implications. Since differentiation among the three sample sites was small, we consider the results to be suggestive of a single stock. However, future studies should evaluate how the stock is affected by differences in sex, age, and gene flow between the nearshore and offshore environments.

Transcriptional profiles of early stage red sea urchins (Mesocentrotus franciscanus) reveal differential regulation of gene expression across development

The red sea urchin, Mesocentrotus franciscanus, is an ecologically important kelp forest species that also serves as a valuable fisheries resource. In this study, we have assembled and annotated a developmental transcriptome for M. franciscanus that represents eggs and six stages of early development (8- to 16-cell, morula, hatched blastula, early gastrula, prism and early pluteus). Characterization of the transcriptome revealed distinct patterns of gene expression that corresponded to major developmental and morphological processes. In addition, the period during which maternally-controlled transcription was terminated and the zygotic genome was activated, the maternal-to-zygotic transition (MZT), was found to begin during early cleavage and persist through the hatched blastula stage, an observation that is similar to the timing of the MZT in other sea urchin species. The presented developmental transcriptome will serve as a useful resource for investigating, in both an ecological and fisheries context, how the early developmental stages of this species respond to environmental stressors.

Genomic and transcriptomic approaches to study immunology in cyprinids: What is next?

Accelerated by the introduction of Next-Generation Sequencing (NGS), a number of genomes of cyprinid fish species have been drafted, leading to a highly valuable collective resource of comparative genome information on cyprinids (Cyprinidae). In addition, NGS-based transcriptome analyses of different developmental stages, organs, or cell types, increasingly contribute to the understanding of complex physiological processes, including immune responses. Cyprinids are a highly interesting family because they comprise one of the most-diversified families of teleosts and because of their variation in ploidy level, with diploid, triploid, tetraploid, hexaploid and sometimes even octoploid species. The wealth of data obtained from NGS technologies provides both challenges and opportunities for immunological research, which will be discussed here. Correct interpretation of ploidy effects on immune responses requires knowledge of the degree of functional divergence between duplicated genes, which can differ even between closely-related cyprinid fish species. We summarize NGS-based progress in analysing immune responses and discuss the importance of respecting the presence of (multiple) duplicated gene sequences when performing transcriptome analyses for detailed understanding of complex physiological processes. Progressively, advances in NGS technology are providing workable methods to further elucidate the implications of gene duplication events and functional divergence of duplicates genes and proteins involved in immune responses in cyprinids. We conclude with discussing how future applications of NGS technologies and analysis methods could enhance immunological research and understanding.

Characterization of natural variation in North American Atlantic Salmon populations (Salmonidae: Salmo salar) at a locus with a major effect on sea age

Age at maturity is a key life‐history trait of most organisms. In anadromous salmonid fishes such as Atlantic Salmon (Salmo salar), age at sexual maturity is associated with sea age, the number of years spent at sea before the spawning migration. For the first time, we investigated the presence of two nonsynonymous vgll3 polymorphisms in North American Atlantic Salmon populations that relate to sea age in European salmon and quantified the natural variation at these and two additional candidate SNPs from two other genes. A targeted resequencing assay was developed and 1,505 returning adult individuals of size‐inferred sea age and sex from four populations were genotyped. Across three of four populations sampled in Québec, Canada, the late‐maturing component (MSW) of the population of a given sex exhibited higher proportions of SNP genotypes 54Thr_vgll3 and 323Lys_vgll3 compared to early‐maturing fish (1SW), for example, 85% versus 53% of females from Trinité River carried 323Lys_vgll3 (n_MSW = 205 vs. n_1SW = 30; p < .001). However, the association between vgll3 polymorphism and sea age was more pronounced in females than in males in the rivers we studied. Logistic regression analysis of vgll3 SNP genotypes revealed increased probabilities of exhibiting higher sea age for 54Thr_vgll3 and 323Lys_vgll3 genotypes compared to alternative genotypes, depending on population and sex. Moreover, individuals carrying the heterozygous vgll3 SNP genotypes were more likely (>66%) to be female. In summary, two nonsynonymous vgll3 polymorphisms were confirmed in North American populations of Atlantic Salmon and our results suggest that variation at those loci correlates with sea age and sex. Our results also suggest that this correlation varies among populations. Future work would benefit from a more balanced sampling and from adding data on juvenile riverine life stages to contrast our data.