Search

Maintaining genetic diversity in cultured shellfish can be challenging due to high variance in individual reproductive success, founder effects, and rapid genetic drift, but is important to retain adaptive potential and avoid inbreeding depression. To support broodstock management and selective breeding in cultured Pacific oysters (Crassostrea (Magallana) gigas), we developed an amplicon panel targeting 592 genomic regions and SNP variants with an average of 50 amplicons per chromosome. Target SNPs were selected based on elevated observed heterozygosity or differentiation in Pacific oyster populations in British Columbia, Canada. The use of the panel for parentage applications was evaluated using multiple generations of oysters from a breeding program on Vancouver Island, Canada (n = 181) and families selected for Ostreid herpesvirus-1 resistance from the Molluscan Broodstock Program in Oregon, USA (n = 136). Population characterization was evaluated using wild, naturalized, farmed, or hatchery oysters sampled throughout the Northern Hemisphere (n = 189). Technical replicates showed high genotype concordance (97.5%; n = 68 replicates). Parentage analysis found suspected pedigree and sample handling errors, demonstrating the panel's value for quality control in breeding programs. Suspected null alleles were identified and found to be largely population dependent, suggesting population-specific variation impacting target amplification. Null alleles were identified using existing data without the need for pedigree information, and once they were removed, assignment rates increased to 93.0 and 86.0% of possible assignments in the two breeding program datasets. A pipeline for analyzing the amplicon sequence data from sequencer output, amplitools , is also provided. [ABSTRACT FROM AUTHOR]

The northern pike Esox lucius is a freshwater fish with low genetic diversity but ecological success throughout the Northern Hemisphere. Here, we generate an annotated chromosome-level genome assembly of 941 Mbp in length with 25 chromosome-length scaffolds. We then genotype 47 northern pike from Alaska through New Jersey at a genome-wide scale and characterize a striking decrease in genetic diversity along the sampling range. Individuals west of the North American Continental Divide have substantially higher diversity than those to the east (e.g. Interior Alaska and St. Lawrence River have on average 181 and 64K heterozygous SNPs per individual, or a heterozygous SNP every 5.2 and 14.6 kbp, respectively). Individuals clustered within each population with strong support, with numerous private alleles observed within each population. Evidence for recent population expansion was observed for a Manitoba hatchery and the St. Lawrence population (Tajima's D = −1.07 and −1.30, respectively). Several chromosomes have large regions with elevated diversity, including LG24, which holds amhby , the ancestral sex determining gene. As expected amhby was largely male-specific in Alaska and the Yukon and absent southeast to these populations, but we document some amhby (−) males in Alaska and amhby (+) males in the Columbia River, providing evidence for a patchwork of presence of this system in the western region. These results support the theory that northern pike recolonized North America from refugia in Alaska and expanded following deglaciation from west to east, with probable founder effects resulting in loss of both neutral and functional diversity (e.g. amhby). [ABSTRACT FROM AUTHOR]

Pacific oyster Crassostrea gigas, endemic to coastal Asia, has been translocated globally throughout the past century, resulting in self‐sustaining introduced populations (naturalized). Oyster aquaculture industries in many parts of the world depend on commercially available seed (hatchery‐farmed) or naturalized/wild oysters to move onto a farm (naturalized‐farmed). It is therefore important to understand genetic variation among populations and farm types. Here, we genotype naturalized/wild populations from France, Japan, China, and most extensively in coastal British Columbia, Canada. We also genotype cultured populations from throughout the Northern Hemisphere to compare with naturalized populations. In total, 16,942 markers were identified using double‐digest RAD‐sequencing in 182 naturalized, 112 hatchery‐farmed, and 72 naturalized‐farmed oysters (n = 366). Consistent with previous studies, very low genetic differentiation was observed around Vancouver Island (mean FST = 0.0019) and low differentiation between countries in the Japan–Canada–France historical translocation lineage (France–Canada FST = 0.0024; Japan–Canada FST = 0.0060). Chinese populations were more differentiated (China–Japan FST = 0.0241). Hatchery‐propagated populations had higher interindividual relatedness suggesting family structure. Within‐population inbreeding was not detected on farms, but nucleotide diversity and polymorphism rate were lower in one farm population. Moving oysters from nature onto farms did not result in strong within‐generation selection. Private alleles at substantial frequency were identified in several hatchery populations grown in BC, suggesting nonlocal origins. Tests of selection identified outlier loci consistent with selective differences associated with domestication, in some cases consistently identified in multiple farms. Top outlier candidates were nearby genes involved in calcium signaling and calmodulin activity. Implications of potential introgression from hatchery‐farmed oysters depend on whether naturalized populations are valued as a locally adapted resource or as an introduced, invasive species. Given the value of the industry in BC and the challenges the industry faces (e.g., climate change, crop losses, biotic stressors), this remains an important question. [ABSTRACT FROM AUTHOR]

Salmon lice Lepeophtheirus salmonis are an ecologically and economically important parasite of wild and farmed salmon. In Scotland, Norway, and Eastern Canada, L. salmonis have developed resistance to emamectin benzoate (EMB), one of the few parasiticides available for salmon lice. Drug resistance mechanisms can be complex, potentially differing among populations and involving multiple genes with additive effects (i.e., polygenic resistance). Indicators of resistance development may enable early detection and countermeasures to avoid the spread of resistance. Here, we collect sensitive Pacific L. salmonis and sensitive and resistant Atlantic L. salmonis from salmon farms, propagate in laboratory (F1), expose to EMB in bioassays, and evaluate either baseline (Atlantic only) or induced transcriptomic differences between populations. In all populations, induced responses were minor and a cellular stress response was not identified. Pacific lice did not upregulate any genes in response to EMB, but downregulated degradative enzymes and transport proteins at 50 ppb EMB. Baseline differences between sensitive and now resistant Atlantic lice were much greater than responses to exposures. All resistant lice overexpressed degradative enzymes, and resistant males, the most resistant group, overexpressed collagenases to the greatest extent. These results indicate an accumulation of baseline expression differences related to resistance. [ABSTRACT FROM AUTHOR]

Divergent life history strategies occur in steelhead or rainbow trout Oncorhynchus mykiss, and many populations produce both migrant (anadromous fish that move to the ocean after rearing) and resident (do not migrate and remain in fresh water) individuals. Mechanisms leading to each type are only partially understood; while the general tendency of a population is heritable, individual tendency may be plastic, influenced by local environment. Steelhead hatchery programmes aim to mitigate losses in wild stocks by producing trout that will migrate to the ocean and not compete with wild trout for limited freshwater resources. To increase our understanding of gill function in these migratory or resident phenotypes, here we compare gill transcriptome profiles of hatchery-released fish either at the release site (residents) or five river kilometres downstream while still in full fresh water (migrants). To test whether any of these genes can be used as predictive markers for smoltification, we compared these genes between migrant-like and undifferentiated trout while still in the hatchery in a common environment (prerelease). Results confirmed the gradual process of smoltification, and the importance of energetics, gill remodelling and ion transport capacity for migrants. Additionally, residents overexpressed transcripts involved in antiviral defences, potentially for immune surveillance via dendritic cells in the gills. The best smoltification marker candidate was protein s100a4, expression of which was highly correlated with Na+, K+ ATPase ( NKA) activity and smolt-like morphology in pre- and postrelease trout gills. [ABSTRACT FROM AUTHOR]

Background: Salmon species vary in susceptibility to infections with the salmon louse (Lepeophtheirus salmonis). Comparing mechanisms underlying responses in susceptible and resistant species is important for estimating impacts of infections on wild salmon, selective breeding of farmed salmon, and expanding our knowledge of fish immune responses to ectoparasites. Herein we report three L. salmonis experimental infection trials of co-habited Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha, profiling hematocrit, blood cortisol concentrations, and transcriptomic responses of the anterior kidney and skin to the infection. Results: In all trials, infection densities (lice per host weight (g)) were consistently highest on chum salmon, followed by Atlantic salmon, and lowest in pink salmon. At 43 days post-exposure, all lice had developed to motile stages, and infection density was uniformly low among species. Hematocrit was reduced in infected Atlantic and chum salmon, and cortisol was elevated in infected chum salmon. Systemic transcriptomic responses were profiled in all species and large differences in response functions were identified between Atlantic and Pacific (chum and pink) salmon. Pink and chum salmon up-regulated acute phase response genes, including complement and coagulation components, and down-regulated antiviral immune genes. The pink salmon response involved the largest and most diverse iron sequestration and homeostasis mechanisms. Pattern recognition receptors were up-regulated in all species but the active components were often species-specific. C-type lectin domain family 4 member M and acidic mammalian chitinase were specifically up-regulated in the resistant pink salmon. Conclusions: Experimental exposures consistently indicated increased susceptibility in chum and Atlantic salmon, and resistance in pink salmon, with differences in infection density occurring within the first three days of infection. Transcriptomic analysis suggested candidate resistance functions including local inflammation with cytokines, specific innate pattern recognition receptors, and iron homeostasis. Suppressed antiviral immunity in both susceptible and resistant species indicates the importance of future work investigating co-infections of viral pathogens and lice. [ABSTRACT FROM AUTHOR]

The Fraser River once supported massive salmon returns. However, over the last century, the largest returns have consistently been less than half of the recorded historical maximum. There is substantial interest from surrounding communities and governments to increase salmon returns for both human use and functional ecosystems. To generate resources for this endeavor, we resequenced genomes of Chinook (Oncorhynchus tshawytscha), coho (Oncorhynchus kisutch), and sockeye salmon (Oncorhynchus nerka) from the Fraser River at moderate coverage (∼16×). A total of 954 resequenced genomes were analyzed, with 681 collected specifically for this study from tissues sampled between 1997 and 2021. An additional 273 were collected from previous studies. At the species level, Chinook salmon appeared to have 1.6–2.1× more SNPs than coho or sockeye salmon, respectively. This difference may be attributable to large historical declines of coho and sockeye salmon. At the population level, 3 Fraser River genetic groups were identified for each species using principal component and admixture analyses. These were consistent with previous research and supports the continued use of these groups in conservation and management efforts. Environmental factors and a migration barrier were identified as major factors influencing the boundaries of these genetic groups. Additionally, 20 potentially adaptive loci were identified among the genetic groups. This information may be valuable in new management and conservation efforts. Furthermore, the resequenced genomes are an important resource for contemporary genomics research on Fraser River salmon and have been made publicly available. [ABSTRACT FROM AUTHOR]

Pacific oysters (Magallana gigas, a.k.a. Crassostrea gigas), the most widely farmed oysters, are under threat from climate change and emerging pathogens. In part, their resilience may be affected by their microbiome, which, in turn, may be influence d by ocean warming and acidification. To understand these impacts, we exposed earlydevelopment Pacific oyster spat to different temperatures (18°C and 24°C) and pCO2 levels (800, 1,600, and 2,800 µatm) in a fully crossed design for 3 weeks. Under all conditions, the microbiome changed over time, with a large decrease in the relative abundance of potentially pathogenic ciliates (Uronema marinum) in all treatments with time. The microbiome composition differed significantly with temperature, but not acidification, indicating that Pacific oyster spat microbiomes can be altered by ocean warming but is resilient to ocean acidification in our experiments. Microbial taxa differed in relative abundance with temperature, implying different adaptive strategies and ecological specializations among microorganisms. Additionally, a small proportion (~0.2% of the total taxa) of the relatively abundant microbial taxa were core constituents (>50% occurrence among samples) across different temperatures, pCO2 levels, or time. Some taxa, including A4b bacteria and members of the family Saprospiraceae in the phyla Chloroflexi (syn. Chloroflexota) and Bacteroidetes (syn. Bacteroidota), respectively, as well as protists in the genera Labyrinthula and Aplanochytrium in the class Labyrinthulomycetes, and Pseudoperkinsus tapetis in the class Ichthyosporea were core constituents across temperatures, pCO2 levels, and time, suggesting that they play an important, albeit unknown, role in maintaining the structural and functional stability of the P acific oyster spat microbiome in response to ocean warming and acidification. These findings highlight the flexibility of the spat microbiome to environmental changes. IMPORTANCE Pacific oysters are the most economically important and widely farmed species of oyster, and their production depends on healthy oyster spat. In turn, spat health and productivity are affected by the associated microbiota; yet, studies have not scrutinized the effects of temperature and pCO2 on the prokaryotic and eukaryotic microbiomes of spat. Here, we show that both the prokaryotic and, for the first time, eukaryotic microbiome of Pacific oyster spat are surprisingly resilient to changes in acidification, but sensitive to ocean warming. The findings have potential implications for oyster survival amid climate change and underscore the need to understand temperature and pCO2 effects on the microbiome and the cascading effects on oyster health and productivity. [ABSTRACT FROM AUTHOR]

Pacific oysters ( Magallana gigas, a.k.a. Crassostrea gigas ), the most widely farmed oysters, are under threat from climate change and emerging pathogens. In part, their resilience may be affected by their microbiome, which, in turn, may be influenced by ocean warming and acidification. To understand these impacts, we exposed early-development Pacific oyster spat to different temperatures (18°C and 24°C) and p CO 2 levels (800, 1,600, and 2,800 µatm) in a fully crossed design for 3 weeks. Under all conditions, the microbiome changed over time, with a large decrease in the relative abundance of potentially pathogenic ciliates ( Uronema marinum ) in all treatments with time. The microbiome composition differed significantly with temperature, but not acidification, indicating that Pacific oyster spat microbiomes can be altered by ocean warming but is resilient to ocean acidification in our experiments. Microbial taxa differed in relative abundance with temperature, implying different adaptive strategies and ecological specializations among microorganisms. Additionally, a small proportion (~0.2% of the total taxa) of the relatively abundant microbial taxa were core constituents (>50% occurrence among samples) across different temperatures, p CO 2 levels, or time. Some taxa, including A4b bacteria and members of the family Saprospiraceae in the phyla Chloroflexi (syn. Chloroflexota ) and Bacteroidetes (syn. Bacteroidota ), respectively, as well as protists in the genera Labyrinthula and Aplanochytrium in the class Labyrinthulomycetes , and Pseudoperkinsus tapetis in the class Ichthyosporea were core constituents across temperatures, p CO 2 levels, and time, suggesting that they play an important, albeit unknown, role in maintaining the structural and functional stability of the Pacific oyster spat microbiome in response to ocean warming and acidification. These findings highlight the flexibility of the spat microbiome to environmental changes.IMPORTANCEPacific oysters are the most economically important and widely farmed species of oyster, and their production depends on healthy oyster spat. In turn, spat health and productivity are affected by the associated microbiota; yet, studies have not scrutinized the effects of temperature and p CO 2 on the prokaryotic and eukaryotic microbiomes of spat. Here, we show that both the prokaryotic and, for the first time, eukaryotic microbiome of Pacific oyster spat are surprisingly resilient to changes in acidification, but sensitive to ocean warming. The findings have potential implications for oyster survival amid climate change and underscore the need to understand temperature and p CO 2 effects on the microbiome and the cascading effects on oyster health and productivity., Competing Interests: The authors declare no conflict of interest.

The Yesso scallop Mizuhopecten yessoensis is an important aquaculture species that was introduced to Western Canada from Japan to establish an economically viable scallop farming industry. This highly fecund species has been propagated in Canadian aquaculture hatcheries for the past 40 years, raising questions about genetic diversity and genetic differences among hatchery stocks. In this study, we compare cultured Canadian and wild Japanese populations of Yesso scallop using double-digest restriction site-associated DNA (ddRAD) sequencing to genotype 21,048 variants in 71 wild-caught scallops from Japan, 65 scallops from the Vancouver Island University breeding population, and 37 scallops obtained from a commercial farm off Vancouver Island, British Columbia. The wild scallops are largely comprised of equally unrelated individuals, whereas cultured scallops are comprised of multiple families of related individuals. The polymorphism rate estimated in wild scallops was 1.7%, whereas in the cultured strains, it ranged between 1.35 and 1.07%. Interestingly, heterozygosity rates were highest in the cultured populations, which is likely due to shellfish hatchery practices of crossing divergent strains to gain benefits of heterosis and to avoid inbreeding. Evidence of founder effects and drift was observed in the cultured strains, including high genetic differentiation between cultured populations and between cultured populations and the wild population. Cultured populations had effective population sizes ranging from 9 to 26 individuals whereas the wild population was estimated at 25,048-56,291 individuals. Further, a depletion of low-frequency variants was observed in the cultured populations. These results indicate significant genetic diversity losses in cultured scallops in Canadian breeding programs., Competing Interests: Conflicts of interest BJGS is affiliated with Sutherland Bioinformatics. The author has no competing financial interests to declare. The other authors declare that no competing interests exist., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Identifying early gene expression responses to hypoxia ( i.e. , low dissolved oxygen) as a tool to assess the degree of exposure to this stressor is crucial for salmonids, because they are increasingly exposed to hypoxic stress due to anthropogenic habitat change, e.g. , global warming, excessive nutrient loading, and persistent algal blooms. Our goal was to discover and validate gill gene expression biomarkers specific to the hypoxia response in salmonids across multi-stressor conditions. Gill tissue was collected from 24 freshwater juvenile Chinook salmon ( Oncorhynchus tshawytscha ), held in normoxia [dissolved oxygen (DO) > 8 mg L -1 ] and hypoxia (DO = 4‒5 mg L -1 ) in 10 and 18° temperatures for up to six days. RNA-sequencing (RNA-seq) was then used to discover 240 differentially expressed genes between hypoxic and normoxic conditions, but not affected by temperature. The most significantly differentially expressed genes had functional roles in the cell cycle and suppression of cell proliferation associated with hypoxic conditions. The most significant genes (n = 30) were selected for real-time qPCR assay development. These assays demonstrated a strong correlation (r = 0.88; P < 0.001) between the expression values from RNA-seq and the fold changes from qPCR. Further, qPCR of the 30 candidate hypoxia biomarkers was applied to an additional 322 Chinook salmon exposed to hypoxic and normoxic conditions to reveal the top biomarkers to define hypoxic stress. Multivariate analyses revealed that smolt stage, water salinity, and morbidity status were relevant factors to consider with the expression of these genes in relation to hypoxic stress. These hypoxia candidate genes will be put into application screening Chinook salmon to determine the identity of stressors impacting the fish., (Copyright © 2020 Akbarzadeh et al.)

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco ( Coregonus artedi ), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts., (Copyright © 2020 Blumstein et al.)

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes ( e.g. , translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr., (Copyright © 2019 Sutherland et al.)

Migration is a ubiquitous life history trait with profound evolutionary and ecological consequences. Recent developments in telemetry and genomics, when combined, can bring significant insights on the migratory ecology of nonmodel organisms in the wild. Here, we used this integrative approach to document dispersal, gene flow and potential for local adaptation in anadromous Arctic Char from six rivers in the Canadian Arctic. Acoustic telemetry data from 124 tracked individuals indicated asymmetric dispersal, with a large proportion of fish (72%) tagged in three different rivers migrating up the same short river in the fall. Population genomics data from 6,136 SNP markers revealed weak, albeit significant, population differentiation (average pairwise F ST  = 0.011) and asymmetric dispersal was also revealed by population assignments. Approximate Bayesian computation simulations suggested the presence of asymmetric gene flow, although in the opposite direction to that observed from the telemetry data, suggesting that dispersal does not necessarily lead to gene flow. These observations suggested that Arctic Char home to their natal river to spawn, but may overwinter in rivers with the shortest migratory route to minimize the costs of migration in nonbreeding years. Genome scans and genetic-environment associations identified 90 outlier markers putatively under selection, 23 of which were in or near a gene. Of these, at least four were involved in muscle and cardiac function, consistent with the hypothesis that migratory harshness could drive local adaptation. Our study illustrates the power of integrating genomics and telemetry to study migrations in nonmodel organisms in logistically challenging environments such as the Arctic., (© 2017 John Wiley & Sons Ltd.)

Whole-genome duplication (WGD) can have large impacts on genome evolution, and much remains unknown about these impacts. This includes the mechanisms of coping with a duplicated sex determination system and whether this has an impact on increasing the diversity of sex determination mechanisms. Other impacts include sexual conflict, where alleles having different optimums in each sex can result in sequestration of genes into nonrecombining sex chromosomes. Sex chromosome development itself may involve sex-specific recombination rate ( i.e. , heterochiasmy), which is also poorly understood. The family Salmonidae is a model system for these phenomena, having undergone autotetraploidization and subsequent rediploidization in most of the genome at the base of the lineage. The salmonid master sex determining gene is known, and many species have nonhomologous sex chromosomes, putatively due to transposition of this gene. In this study, we identify the sex chromosome of Brook Charr Salvelinus fontinalis and compare sex chromosome identities across the lineage (eight species and four genera). Although nonhomology is frequent, homologous sex chromosomes and other consistencies are present in distantly related species, indicating probable convergence on specific sex and neo-sex chromosomes. We also characterize strong heterochiasmy with 2.7-fold more crossovers in maternal than paternal haplotypes with paternal crossovers biased to chromosome ends. When considering only rediploidized chromosomes, the overall heterochiasmy trend remains, although with only 1.9-fold more recombination in the female than the male. Y chromosome crossovers are restricted to a single end of the chromosome, and this chromosome contains a large interspecific inversion, although its status between males and females remains unknown. Finally, we identify quantitative trait loci (QTL) for 21 unique growth, reproductive, and stress-related phenotypes to improve knowledge of the genetic architecture of these traits important to aquaculture and evolution., (Copyright © 2017 Sutherland et al.)