Revealing the evolutionary history and contemporary population structure of Pacific salmon in the Fraser River through genome resequencing.

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 (O. kisutch), and sockeye salmon (O. nerka) from the Fraser River at moderate coverage (∼16x). 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.1x 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, three Fraser River genetic groups were identified for each species using principal component and admixture analyses, which is 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.

An amplicon panel for high-throughput and low-cost genotyping of Pacific oyster.

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 = 190). 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.

The prokaryotic and eukaryotic microbiome of Pacific oyster spat is shaped by ocean warming but not acidification.

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 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 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 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.

Population genomics of harbour seal Phoca vitulina from northern British Columbia through California and comparison to the Atlantic subspecies.

The harbour seal Phoca vitulina is a ubiquitous pinniped species found throughout coastal waters of the Northern Hemisphere. Harbour seal impacts on ecosystem dynamics may be significant due to their high abundance and food web position. Two subspecies exist in North America, P. v. richardii in the Pacific Ocean and P. v. vitulina in the Atlantic. Strong natal philopatry of harbour seals can result in fine-scale genetic structure and isolation by distance. Management of harbour seals is expected to benefit from improved resolution of seal population structure and dynamics. Here, we use genotyping-by-sequencing to genotype 146 harbour seals from the eastern Pacific Ocean (i.e. British Columbia (BC), Oregon and California) and the western Atlantic Ocean (i.e. Québec, Newfoundland and Labrador). Using 12,742 identified variants, we confirm the recently identified elevated genetic diversity in the eastern Pacific relative to the western Atlantic and greatest differentiation between the subspecies. Further, we demonstrate that this is independent of reference genome bias or other potential technical artefacts. Coast-specific analyses with 8933 and 3828 variants in Pacific and Atlantic subspecies, respectively, identify divergence between BC and Oregon-California, and between Québec and Newfoundland-Labrador. Unexpected PCA outlier clusters were observed in two populations due to cryptic relatedness of individuals; subsequently, closely related samples were removed. Admixture analysis indicates an isolation-by-distance signature where Oregon seals contained some of the BC signature, whereas California did not. Additional sampling is needed in the central and north coast of BC to determine whether a discrete separation of populations exists within the region.

Biochemical and Cellular Characterization of the Function of Fluorophosphonate-Binding Hydrolase H (FphH) in Staphylococcus aureus Support a Role in Bacterial Stress Response.

The development of new treatment options for bacterial infections requires access to new targets for antibiotics and antivirulence strategies. Chemoproteomic approaches are powerful tools for profiling and identifying novel druggable target candidates, but their functions often remain uncharacterized. Previously, we used activity-based protein profiling in the opportunistic pathogen Staphylococcus aureus to identify active serine hydrolases termed fluorophosphonate-binding hydrolases (Fph). Here, we provide the first characterization of S. aureus FphH, a conserved, putative carboxylesterase (referred to as yvaK in Bacillus subtilis) at the molecular and cellular level. First, phenotypic characterization of fphH-deficient transposon mutants revealed phenotypes during growth under nutrient deprivation, biofilm formation, and intracellular survival. Biochemical and structural investigations revealed that FphH acts as an esterase and lipase based on a fold well suited to act on a small to long hydrophobic unbranched lipid group within its substrate and can be inhibited by active site-targeting oxadiazoles. Prompted by a previous observation that fphH expression was upregulated in response to fusidic acid, we found that FphH can deacetylate this ribosome-targeting antibiotic, but the lack of FphH function did not infer major changes in antibiotic susceptibility. In conclusion, our results indicate a functional role of this hydrolase in S. aureus stress responses, and hypothetical functions connecting FphH with components of the ribosome rescue system that are conserved in the same gene cluster across Bacillales are discussed. Our atomic characterization of FphH will facilitate the development of specific FphH inhibitors and probes to elucidate its physiological role and validity as a drug target.

Genomic diversity of wild and cultured Yesso scallop Mizuhopecten yessoensis from Japan and Canada.

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.

Chinook salmon diversity contributes to fishery stability and trade-offs with mixed-stock harvest.

Variation among populations in life history and intrinsic population characteristics (i.e., population diversity) helps maintain resilience to environmental change and dampen interannual variability in ecosystem services. As a result, ecological variation, and the processes that generate it, is considered central to strategies for managing risks to ecosystems in an increasingly variable and uncertain world. However, characterizing population diversity is difficult, particularly in large and remote regions, which often prevents its formal consideration in management advice. We combined genetic stock identification of archived scale and tissue samples with state-space run-reconstruction models to estimate migration timing and annual return abundance for eight geographically and genetically distinct Chinook salmon populations within the Canadian portion of the Yukon River. We found that among-population variation in migration timing and return abundances resulted in aggregate return migrations that were 2.1 times longer and 1.4 times more stable than if they had composed a single homogeneous population. We then fit state-space spawner-recruitment models to the annual return abundances to characterize among-population diversity in intrinsic productivity and population size and their consequences for the fisheries they support. Productivity and carrying capacity varied among populations by approximately 2.4-fold (2.9 to 6.9 recruits per spawner) and three-fold (8800 to 27,000 spawners), respectively. This diversity implies an equilibrium trade-off between harvesting of the population aggregate and the conservation of individual populations whereby the harvest rate predicted to maximize aggregate harvests comes at the cost of overfishing ~40% of the populations but with a relatively low risk of extirpating the weakest ones. Our findings illustrate how population diversity in one of the largest salmon-producing river basins in the world contributes to fishery stability and food security in a region where salmon have high cultural and subsistence value. More generally, our work demonstrates the utility of molecular analyses of archived biological material for characterizing diversity in biological systems and its benefits and consequences for trade-offs in decision-making.

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.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors.

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

The development and psychometric testing of three instruments that measure person-centred caring as three concepts - Personalization, participation and responsiveness.

AIM: To develop and test the psychometric properties of three instruments that measure Person-centred Caring: as Personalization, Participation and Responsiveness. DESIGN: A three-phase mixed methods design used two frameworks: content validity determination and quantification; consensus-based standards for selection of health measurement instruments. METHODS: A narrative literature review identified the domain definition. A systematic review of instruments provided the basis for item pools, which were refined by focus groups (N = 4) of multidisciplinary staff and service users (N = 25) and cognitive interviews (N = 11) with service users. Scale content validity indexes were calculated. Three cross-sectional surveys were conducted between April 2015 and June 2016. The instruments' psychometric properties tested included factor structure, internal consistency and construct validity. Convergent validity was tested, hypothesizing that: Personalization related to relational empathy; Participation related to empowerment; and Responsiveness related to trust. RESULTS: Scale content validity indexes were ≥0.96 in all instruments. Response rates were 24% (N = 191), 15% (N = 108) and 19% (N = 124). Two factors were revealed for the Personalization and Responsiveness instruments and one factor for the Participation instrument. All had acceptable: reliability (Cronbach's Alpha >0.7); construct validity (>50%); and convergent validity (Spearman's correlation coefficient >0.25, p < 0.05). CONCLUSION: This study composed definitions and instruments that reflect the multidisciplinary teams' caring behaviours, which have acceptable reliability and validity in the community population. Further psychometric testing of Participation and Responsiveness instruments should be undertaken with a larger sample. IMPACT: The instruments can be used to monitor the variability of multidisciplinary teams' caring behaviours; research effective interventions to improve caring behaviours; and increase understanding of the impact of caring on health outcomes.

Relative genomic impacts of translocation history, hatchery practices, and farm selection in Pacific oyster Crassostrea gigas throughout the Northern Hemisphere.

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 F ST = 0.0019) and low differentiation between countries in the Japan-Canada-France historical translocation lineage (France-Canada F ST = 0.0024; Japan-Canada F ST = 0.0060). Chinese populations were more differentiated (China-Japan F ST = 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.

Identification of Hypoxia-Specific Biomarkers in Salmonids Using RNA-Sequencing and Validation Using High-Throughput qPCR.

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.

Comparative Genomic Analyses and a Novel Linkage Map for Cisco (Coregonus artedi) Provide Insights into Chromosomal Evolution and Rediploidization Across Salmonids.

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.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis.

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.

Genomics and telemetry suggest a role for migration harshness in determining overwintering habitat choice, but not gene flow, in anadromous Arctic Char.

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 FST  = 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.

Effects of the vertically transmitted microsporidian Facilispora margolisi and the parasiticide emamectin benzoate on salmon lice (Lepeophtheirus salmonis).

BACKGROUND: Microsporidia are highly specialized, parasitic fungi that infect a wide range of eukaryotic hosts from all major taxa. Infections cause a variety of damaging effects on host physiology from increased stress to death. The microsporidian Facilispora margolisi infects the Pacific salmon louse (Lepeophtheirus salmonis oncorhynchi), an economically and ecologically important ectoparasitic copepod that can impact wild and cultured salmonids. RESULTS: Vertical transmission of F. margolisi was demonstrated by using PCR and in situ hybridization to identify and localize microsporidia in female L. salmonis and their offspring. Spores and developmental structures of F. margolisi were identified in 77% of F1 generation copepods derived from infected females while offspring from uninfected females all tested negative for the microsporidia. The transcriptomic response of the salmon louse to F. margolisi was profiled at both the copepodid larval stage and the pre-adult stage using microarray technology. Infected copepodids differentially expressed 577 transcripts related to stress, ATP generation and structural components of muscle and cuticle. The infection also impacted the response of the copepodid to the parasiticide emamectin benzoate (EMB) at a low dose of 1.0 ppb for 24 h. A set of 48 transcripts putatively involved in feeding and host immunomodulation were up to 8-fold underexpressed in the F. margolisi infected copepodids treated with EMB compared with controls or either stressor alone. Additionally, these infected lice treated with EMB also overexpressed 101 transcripts involved in stress resistance and signalling compared to the other groups. In contrast, infected pre-adult lice did not display a stress response, suggesting a decrease in microsporidian virulence associated with lice maturity. Furthermore, copepodid infectivity and moulting was not affected by the microsporidian infection. CONCLUSIONS: This study demonstrated that F. margolisi is transmitted vertically between salmon louse generations and that biological impacts of infection differ depending on the stage of the copepod host. The infection caused significant perturbations of larval transcriptomes and therefore must be considered in future studies in which impacts to host development and environmental factors are assessed. Fitness impacts are probably minor, although the interaction between pesticide exposure and microsporidian infection merits further study.

Sex Chromosome Evolution, Heterochiasmy, and Physiological QTL in the Salmonid Brook Charr Salvelinus fontinalis.

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.

Sex matters in massive parallel sequencing: Evidence for biases in genetic parameter estimation and investigation of sex determination systems.

Using massively parallel sequencing data from two species with different life history traits, American lobster (Homarus americanus) and Arctic Char (Salvelinus alpinus), we highlight how an unbalanced sex ratio in the samples and a few sex-linked markers may lead to false interpretations of population structure and thus to potentially erroneous management recommendations. Here, multivariate analyses revealed two genetic clusters separating samples by sex instead of by expected spatial variation: inshore and offshore locations in lobster, or east and west locations in Arctic Char. To further investigate this, we created several subsamples artificially varying the sex ratio in the inshore/offshore and east/west groups and then demonstrated that significant genetic differentiation could be observed despite panmixia in lobster, and that FST values were overestimated in Arctic Char. This pattern was due to 12 and 94 sex-linked markers driving differentiation for lobster and Arctic Char, respectively. Removing sex-linked markers led to nonsignificant genetic structure in lobster and a more accurate estimation of FST in Arctic Char. The locations of these markers and putative identities of genes containing or nearby the markers were determined using available transcriptomic and genomic data, and this provided new information related to sex determination in both species. Given that only 9.6% of all marine/diadromous population genomic studies to date have reported sex information, we urge researchers to collect and consider individual sex information. Sex information is therefore relevant for avoiding unexpected biases due to sex-linked markers as well as for improving our knowledge of sex determination systems in nonmodel species.

Genomic patterns of diversity and divergence of two introduced salmonid species in Patagonia, South America.

Invasive species have become widespread in aquatic environments throughout the world, yet there are few studies that have examined genomic variation of multiple introduced species in newly colonized environments. In this study, we contrast genomic variation in two salmonid species (anadromous Chinook Salmon, Oncorhynchus tshawytscha, 11,579 SNPs and resident Brook Charr Salvelinus fontinalis, 13,522 SNPs) with differing invasion success after introduction to new environments in South America relative to populations from their native range in North America. Estimates of genetic diversity were not significantly different between introduced and source populations for either species, indicative of propagule pressure that has been shown to maintain diversity in founding populations relative to their native range. Introduced populations also demonstrated higher connectivity and gene flow than those in their native range. Evidence for candidate loci under divergent selection was observed, but was limited to specific introduced populations and was not widely evident. Patterns of genomic variation were consistent with general dispersal potential of each species and therefore also the notion that life history variation may contribute to both invasion success and subsequent genetic structure of these two salmonids in Patagonia.

Enhanced transcriptomic responses in the Pacific salmon louse Lepeophtheirus salmonis oncorhynchi to the non-native Atlantic Salmon Salmo salar suggests increased parasite fitness.

BACKGROUND: Outcomes of infections with the salmon louse Lepeophtheirus salmonis vary considerably among its natural hosts (Salmo, Oncorhynchus spp.). Host-parasite interactions range from weak to strong host responses accompanied by high to low parasite abundances, respectively. Parasite behavioral studies indicate that the louse prefers the host Atlantic Salmon (Salmo salar), which is characterized by a weak immune response, and that this results in enhanced parasite reproduction and growth rates. Furthermore, parasite-derived immunosuppressive molecules (e.g., proteases) have been detected at higher amounts in response to the mucus of Atlantic Salmon relative to Coho Salmon (Oncorhynchus kisutch). However, the host-specific responses of the salmon louse have not been well characterized in either of the genetically distinct sub-species that occur in the Atlantic and Pacific Oceans. RESULTS: We assessed and compared the transcriptomic feeding response of the Pacific salmon louse (L. salmonis oncorhynchi,) while parasitizing the highly susceptible Atlantic Salmon and Sockeye Salmon (Oncorhynchus nerka) or the more resistant Coho Salmon (Oncorhynchus kisutch) using a 38 K oligonucleotide microarray. The response of the louse was enhanced both in the number of overexpressed genes and in the magnitude of expression while feeding on the non-native Atlantic Salmon, compared to either Coho or Sockeye Salmon. For example, putative virulence factors (e.g., cathepsin L, trypsin, carboxypeptidase B), metabolic enzymes (e.g., cytochrome B, cytochrome C), protein synthesis enzymes (e.g., ribosomal protein P2, 60S ribosomal protein L7), and reproduction-related genes (e.g., estrogen sulfotransferase) were overexpressed in Atlantic-fed lice, indicating heightened parasite fitness with this host species. In contrast, responses in Coho- or Sockeye-fed lice were more similar to those of parasites deprived of a host. To test for host acclimation by the parasite, we performed a reciprocal host transfer experiment and determined that the exaggerated response to Atlantic Salmon was independent of the initial host species, confirming our conclusion that the Pacific salmon louse exhibits an enhanced response to Atlantic Salmon. CONCLUSIONS: This study characterized global transcriptomic responses of Pacific salmon lice during infection of susceptible and resistant hosts. Similar parasite responses during infection of Coho or Sockeye Salmon, despite differences in natural immunity to infection between these host species, indicate that host susceptibility status alone does not drive the parasite response. We identified an enhanced louse response after feeding on Atlantic Salmon, characterized by up-regulation of virulence factors, energy metabolism and reproductive-associated transcripts. In contrast, the responses of lice infecting Coho or Sockeye Salmon were weaker, with reduced expression of virulence factors. These observations indicate that the response of the louse is independent of host susceptibility and suggest that co-evolutionary host-parasite relationships may influence contemporary host-parasite interactions. This research improves our understanding of the susceptibility of Atlantic Salmon and may assist in the development of novel control measures against the salmon louse.