Genetic parameters and genotype-by-environment interaction estimates for growth and feed efficiency related traits in Chinook salmon, Oncorhynchus tshawytscha, reared under low and moderate flow regimes.

BACKGROUND: A genotype-by-environment (G × E) interaction is defined as genotypes responding differently to different environments. In salmonids, G × E interactions can occur in different rearing conditions, including changes in salinity or temperature. However, water flow, an important variable that can influence metabolism, has yet to be considered for potential G × E interactions, although water flows differ across production stages. The salmonid industry is now manipulating flow in tanks to improve welfare and production performance, and expanding sea pen farming offshore, where flow dynamics are substantially greater. Therefore, there is a need to test whether G × E interactions occur under low and higher flow regimes to determine if industry should consider modifying their performance evaluation and selection criteria to account for different flow environments. Here, we used genotype-by-sequencing to create a genomic-relationship matrix of 37 Chinook salmon, Oncorhynchus tshawytscha, families to assess possible G × E interactions for production performance under two flow environments: a low flow regime (0.3 body lengths per second; bl s-1) and a moderate flow regime (0.8 bl s-1). RESULTS: Genetic correlations for the same production performance trait between flow regimes suggest there is minimal evidence of a G × E interaction between the low and moderate flow regimes tested in this study, for Chinook salmon reared from 82.9 ± 16.8 g ( x ¯ ± s.d.) to 583.2 ± 117.1 g ( x ¯ ± s.d.). Estimates of genetic and phenotypic correlations between traits did not reveal any unfavorable trait correlations for size- (weight and condition factor) and growth-related traits, regardless of the flow regime, but did suggest measuring feed intake would be the preferred approach to improve feed efficiency because of the strong correlations between feed intake and feed efficiency, consistent with previous studies. CONCLUSION: This new information suggests that Chinook salmon families do not need to be selected separately for performance across different flow regimes. However, further studies are needed to confirm this across a wider range of fish sizes and flows. This information is key for breeding programs to determine if separate evaluation groups are required for different flow regimes that are used for production (e.g., hatchery, post smolt recirculating aquaculture system, or offshore).

Integrating otolith and genetic tools to reveal intraspecific biodiversity in a highly impacted salmon population.

Intraspecific biodiversity is vital for species persistence in an increasingly volatile world. By embracing methods that integrate information at different spatiotemporal scales, we can directly monitor and reconstruct changes in intraspecific biodiversity. Here we combined genetics and otolith biochronologies to describe the genotypic and phenotypic diversity of Chinook salmon (Oncorhynchus tshawytscha) in the Yuba River, California, comparing cohorts that experienced a range of hydroclimatic conditions. Yuba River salmon have been heavily impacted by habitat loss and degradation, and large influxes of unmarked hatchery fish each year have led to concern about introgression and uncertainty around the viability of its wild populations, particularly the rarer spring-run salmon. Otolith strontium isotopes showed that Yuba River origin fish represented, on average, 42% (range 7%-73%) of spawners across six return years (2009-2011, 2018-2020), with large interannual variability. The remainder of adult Chinook salmon in the river were primarily strays from the nearby Feather River hatchery, and since 2018 from the Mokelumne River hatchery. Among the Yuba-origin spawners, on average, 30% (range 14%-50%) exhibited the spring-run genotype. The Yuba-origin fish also displayed a variety of outmigration phenotypes that differed in the timing and size at which they left the Yuba river. Early-migrating fry dominated the returns (mean 59%, range 33%-89%), and their contribution rates were negatively correlated with freshwater flows. It is unlikely that fry survival rates are elevated during droughts, suggesting that this trend reflects disproportionately low survival of larger later migrating parr, smolts, and yearlings along the migratory corridor in drier years. Otolith daily increments indicated generally faster growth rates in non-natal habitats, emphasizing the importance of continuing upstream restoration efforts to improve in-river growing conditions. Together, these findings show that, despite a long history of habitat degradation and hatchery introgression, the Yuba River maintains intraspecific biodiversity that should be taken into account in future management, restoration, and reintroduction plans. The finding that genotypic spring-run are reproducing, surviving, and returning to the Yuba River every year suggests that re-establishment of an independent population is possible, although hatchery-wild interactions would need to be carefully considered. Integrating methods is critical to monitor changes in key genetic, physiological, and behavioral traits to assess population viability and resilience.

Behavioral transcriptomic effects of triploidy and probiotic therapy (Bifidobacterium, Lactobacillus, and Lactococcus mixture) on juvenile Chinook salmon (Oncorhynchus tshawytscha).

Aquaculturists use polyploid fish to maximize production albeit with some unintended consequences including compromised behaviors and physiological function. Given benefits of probiotic therapies (e.g., improved immune response, growth, and metabolism), we explored probiotic supplementation (mixture of Bifidobacterium, Lactobacillus, and Lactococcus), to overcome drawbacks. We first examined fish gut bacterial community composition using 16S metabarcoding (via principal coordinate analyses and PERMANOVA) and determined probiotics significantly impacted gut bacteria composition (p = 0.001). Secondly, we examined how a genomic disruptor (triploidy) and diet supplements (probiotics) impact gene transcription and behavioral profiles of hatchery-reared Chinook salmon (Oncorhynchus tshawytscha). Juveniles from four treatment groups (diploid-regular feed, diploid-probiotic feed, triploid-regular feed, and triploid-probiotic feed; n = 360) underwent behavioral assays to test activity, exploration, neophobia, predator evasion, aggression/sociality, behavioral sensitivity, and flexibility. In these fish, transcriptional profiles for genes associated with neural functions (neurogenesis/synaptic plasticity) and biomarkers for stress response and development (growth/appetite) were (i) examined across treatments and (ii) used to describe behavioral phenotypes via principal component analyses and general linear mixed models. Triploids exhibited a more active behavioral profile (p = 0.002), and those on a regular diet had greater Neuropeptide Y transcription (p = 0.02). A growth gene (early growth response protein 1, p = 0.02) and long-term neural development genes (neurogenic differentiation factor, p = 0.003 and synaptysomal-associated protein 25-a, p = 0.005) impacted activity and reactionary profiles, respectively. Overall, our probiotic treatment did not compensate for triploidy. Our research highlights novel applications of behavioral transcriptomics for identifying candidate genes and dynamic, mechanistic associations with complex behavioral repertoires.

Genetic architecture and correlations between the gut microbiome and gut gene transcription in Chinook salmon (Oncorhynchus tshawytscha).

Population divergence through selection can drive local adaptation in natural populations which has implications for the effective restoration of declining and extirpated populations. However, adaptation to local environmental conditions is complicated when both the host and its associated microbiomes must respond via co-evolutionary change. Nevertheless, for adaptation to occur through selection, variation in both host and microbiome traits should include additive genetic effects. Here we focus on host immune function and quantify factors affecting variation in gut immune gene transcription and gut bacterial community composition in early life-stage Chinook salmon (Oncorhynchus tshawytscha). Specifically, we utilized a replicated factorial breeding design to determine the genetic architecture (sire, dam and sire-by-dam interaction) of gut immune gene transcription and microbiome composition. Furthermore, we explored correlations between host gut gene transcription and microbiota composition. Gene transcription was quantified using nanofluidic qPCR arrays (22 target genes) and microbiota composition using 16 S rRNA gene (V5-V6) amplicon sequencing. We discovered limited but significant genetic architecture in gut microbiota composition and transcriptional profiles. We also identified significant correlations between gut gene transcription and microbiota composition, highlighting potential mechanisms for functional interactions between the two. Overall, this study provides support for the co-evolution of host immune function and their gut microbiota in Chinook salmon, a species recognized as locally adapted. Thus, the inclusion of immune gene transcription profile and gut microbiome composition as factors in the development of conservation and commercial rearing practices may provide new and more effective approaches to captive rearing.

Insights from Hi-C data regarding the Pacific salmon louse (Lepeophtheirus salmonis) sex chromosomes.

Salmon lice, Lepeophtheirus salmonis (family Caligidae), are ectoparasites that have negatively impacted the salmon aquaculture industry and vulnerable wild salmon populations. Researchers have studied salmon lice to better understand their biology to develop effective control strategies. In this study, we updated the chromosome-level reference genome assembly of the Pacific subspecies of L. salmonis using Hi-C data. The previous version placed contigs/scaffolds using an Atlantic salmon louse genetic map. By utilizing Hi-C data from Pacific salmon lice, we were able to properly assign locations to contigs/scaffolds previously unplaced or misplaced. This resulted in a more accurate genome assembly and a more comprehensive characterization of the sex chromosome unique to females (W). We found evidence that the same ZW-ZZ mechanism is common in both Atlantic and Pacific subspecies of salmon lice using PCR assays. The W chromosome was approximately 800 kb in size, which is ∼30 times smaller than the Z chromosome (24 Mb). The W chromosome contained 61 annotated genes, including 32 protein-coding genes, 27 long noncoding RNA (lncRNA) genes, and 2 pseudogenes. Among these 61 genes, 39 genes shared homology to genes found on other chromosomes, while 20 were unique to the W chromosome. Two genes of interest on the W chromosome, prohibitin-2 and kinase suppressor of ras-2, were previously identified as potential sex-linked markers in the salmon louse. However, we prioritized the 20 unique genes on the W chromosome as sex-determining candidates. This information furthers our understanding of the biology of this ectoparasite and may help in the development of more effective management strategies.

Invasion status of hatchery-origin pink salmon in an unstocked river at the Shiretoko World Natural Heritage Site in northern Japan.

Hatchery fish and their offspring (including hatchery-wild hybrids) have lower reproductive success than wild fish. Thus, the straying of hatchery fish may negatively impact wild populations, depending on the number of wild salmon returning and hatchery strays. We investigated the straying status of hatchery-origin pink salmon (Oncorhynchus gorbuscha), which have a higher straying rate than other salmonids, in an unstocked river at the Shiretoko World Natural Heritage Site, Japan. The hatchery strays accounted for 40.0% and 19.0% of the total samples in 2021 and 2022, respectively. These results indicate that hatchery pink salmon have invaded unstocked rivers and potentially genetically affect wild populations.

Salmon sperm DNA increases sample recovery from cotton swabs.

Forensic samples with low DNA template amounts are difficult to analyze and interpret. There is a large body of research demonstrating that adding carrier nucleic acid to storage tubes, solid phase extractions, or filtering devices can improve yields of target DNA. However, the addition of carrier nucleic acid to sampling substrates, like cotton swabs, has not yet been attempted. In this proof-of-concept study, carrier nucleic acids in the form of either Poly (A) RNA or salmon sperm DNA were spotted onto cotton swabs, followed by human genomic DNA, to determine if introducing the carrier prior to sample collection would increase recovery from the swabs post-extraction. Extracts were also evaluated to determine whether adding the carrier nucleic acids to human DNA would interfere with downstream forensic DNA analysis processes such as real-time PCR quantitation, PCR amplification of STR loci, or capillary electrophoresis. The RNA carrier did not improve human sample recovery from cotton swabs. The extraction efficiency of human DNA from cotton swabs was increased when the DNA carrier was applied to the swabs prior to sample deposition, and the scale of the increase depended on the amount of carrier DNA used. When applying the salmon sperm DNA carrier to cotton swabs, with each increase from no carrier to 0.001-1-10 µg, human DNA recovery went from ∼29 % to ∼50 % to ∼75 % to ∼100 %. Additionally, no inhibitory effects from the carrier DNA were observed post-extraction with quantitation or in the DNA profile after amplification. Therefore, salmon sperm DNA carrier will increase human DNA yield from cotton swabs without negative effects on downstream forensic DNA profiling methods, with the optimal carrier amount being 10 µg.

Rescue of Infectious Salmon Anemia Virus (ISAV) from Cloned cDNA.

The piscine orthomyxovirus called infectious salmon anemia virus (ISAV) is one of the most important emerging pathogens affecting the salmon industry worldwide. The first reverse genetics system for ISAV, which allows the generation of recombinant ISA virus (rISAV), is an important tool for the characterization and study of this virus. The plasmid-based reverse genetics system for ISAV includes the use of a novel fish promoter, the Atlantic salmon internal transcribed spacer region 1 (ITS-1). The salmon, viral, and mammalian genetic elements included in the pSS-URG vectors allow the expression of the eight viral RNA segments. In addition to four cytomegalovirus (CMV)-based vectors that express the four proteins of the ISAV ribonucleoprotein complex, the eight pSS-URG vectors allowed the generation of infectious rISAV in salmon cells.

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

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

ISOTHERMAL RECOMBINANT POLYMERASE AMPLIFICATION AND CRIPSR(CAS12A) ASSAY DETECTION OF RENIBACTERIUM SALMONINARUM AS AN EXAMPLE FOR WILDLIFE PATHOGEN DETECTION IN ENVIRONMENTAL DNA SAMPLES.

Improving rapid detection methods for pathogens is important for research as we collectively aim to improve the health of ecosystems globally. In the northern hemisphere, the success of salmon (Oncorhynchus spp.) populations is vitally important to the larger marine, aquatic, and terrestrial ecosystems they inhabit. This has led to managers cultivating salmon in hatcheries and aquaculture to bolster their populations, but young salmon face many challenges, including diseases such as bacterial kidney disease (BKD). Early detection of the BKD causative agent, Renibacterium salmoninarum, is useful for managers to avoid outbreaks in hatcheries and aquaculture stocks to enable rapid treatment with targeted antibiotics. Isothermal amplification and CRIPSR-Cas12a systems may enable sensitive, relatively rapid, detection of target DNA molecules from environmental samples compared to quantitative PCR (qPCR) and culture methods. We used these technologies to develop a sensitive and specific rapid assay to detect R. salmoninarum from water samples using isothermal recombinase polymerase amplification (RPA) and an AsCas12a RNA-guided nuclease detection. The assay was specific to R. salmoninarum (0/10 co-occurring or closely related bacteria detected) and sensitive to 0.0128 pg/µL of DNA (approximately 20-40 copies/µL) within 10 min of Cas activity. This assay successfully detected R. salmoninarum environmental DNA in 14/20 water samples from hatcheries with known quantification for the pathogen via previous qPCR (70% of qPCR-positive samples). The RPA-CRISPR/AsCas12a assay had a limit of detection (LOD) of >10 copies/µL in the hatchery water samples and stochastic detection below 10 copies/µL, similar to but slightly higher than the qPCR assay. This LOD enables 37 C isothermal detection, potentially in the field, of biologically relevant levels of R. salmoninarum in water. Further research is needed to develop easy-to-use, cost-effective, sensitive RPA/CRISPR-AsCas12a assays for rapidly detecting low concentrations of wildlife pathogens in environmental samples.

Tidal gradients, fine-scale homing and a potential cryptic ecotype of wild spawning pink salmon (Oncorhynchus gorbuscha).

The homing behaviour of salmon is a remarkable natural phenomenon, critical for shaping the ecology and evolution of populations yet the spatial scale at which it occurs is poorly understood. This study investigated the spatial scale and mechanisms driving homing as depicted by spawning site-choice behaviour in pink salmon (Oncorhynchus gorbuscha) in Prince William Sound, Alaska. Molecular pedigree analyses of over 30,000 adult spawners in four streams revealed that pink salmon exhibit fine-scale site fidelity within a stream, returning to within <100 m of their parents. Homing behaviours were driven in part by a salinity gradient between intertidal and freshwater environments, with individuals incubated in freshwater environments more than twice as likely to spawn upstream of tidal influence than those incubated in the intertidal. Our findings challenge the traditional view that pink salmon populations are genetically and phenotypically homogenous due to their short freshwater residency as juveniles and high rates of dispersal as returning adults (i.e. straying). This study has important implications for rates of inbreeding, local adaptation and gene flow within populations, and is particularly relevant to the management of salmon hatcheries, given the high incidence of hatchery-origin pink salmon, reared in freshwater hatchery environments, that stray into wild populations of Prince William Sound.

Diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) have altered microRNA responses in immune tissues after infection with Vibrio anguillarum.

Production of sterile fishes through artificial retention of a third set of chromosomes (triploidy) is a sustainable alternative for aquaculture since it reduces escapee pressure on wild populations. However, these fishes have reduced survival in stressful conditions and in response to infection. In this study, the impact of Vibrio anguillarum infection on diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) was investigated to identify if there was any significant immune regulation by microRNAs (miRNA). Small RNAs from hindgut, head kidney, and spleen were sequenced to determine if miRNA transcript abundance was altered due to ploidy and infection in nine-month old full-sibling diploids and triploids. All three tissues had differentially expressed miRNA prior to infection, indicating subtle changes in epigenetic regulation due to increased ploidy. Additionally, miRNA were altered by infection, but there was only a difference in spleen miRNA expression between diploids and triploids at three days of infection. Furthermore, one miRNA (ssa-miR-2188-3p) was confirmed as having an altered response to infection in triploids compared to diploids, implicating potential immune dysregulation due to increased ploidy. The miRNAs identified in this study are predicted to target immune pathways, providing evidence for their importance in regulating responses to pathogens. This study is the first to investigate how increased ploidy alters miRNA expression in response to infection. Additionally, it provides evidence for epigenetic dysregulation in triploid fishes, which may contribute to their poor performance in response to stress.

CRISPR-Cas- induced IRF3 and MAVS knockouts in a salmonid cell line disrupt PRR signaling and affect viral replication.

Background: Interferon (IFN) responses are critical in the resolution of viral infections and are actively targeted by many viruses. They also play a role in inducing protective responses after vaccination and have been successfully tested as vaccine adjuvants. IFN responses are well conserved and function very similar in teleosts and mammals. Like in mammals, IFN responses in piscine cells are initiated by intracellular detection of the viral infection by different pattern recognition receptors. Upon the recognition of viral components, IFN responses are rapidly induced to combat the infection. However, many viruses may still replicate and be able to inhibit or circumvent the IFN response by different means. Methods: By employing CRISPR Cas9 technology, we have disrupted proteins that are central for IFN signaling in the salmonid cell line CHSE-214. We successfully generated KO clones for the mitochondrial antiviral signaling protein MAVS, the transcription factors IRF3 and IRF7-1, as well as a double KO for IRF7-1/3 using an optimized protocol for delivery of CRISPR-Cas ribonucleoproteins through nucleofection. Results: We found that MAVS and IRF3 KOs inhibited IFN and IFN-stimulated gene induction after intracellular poly I:C stimulation as determined through gene expression and promoter activation assays. In contrast, the IRF7-1 KO had no clear effect. This shows that MAVS and IRF3 are essential for initiation of intracellular RNA-induced IFN responses in CHSE-214 cells. To elucidate viral interference with IFN induction pathways, the KOs were infected with Salmon alphavirus 3 (SAV3) and infectious pancreatic necrosis virus (IPNV). SAV3 infection in control and IRF7-1 KO cells yielded similar titers and no cytopathic effect, while IRF3 and MAVS KOs presented with severe cytopathic effect and increased titers 6 days after SAV 3 infection. In contrast, IPNV yields were reduced in IRF3 and MAVS KOs, suggesting a dependency on interactions between viral proteins and pattern recognition receptor signaling components during viral replication. Conclusion: Aside from more insight in this signaling in salmonids, our results indicate a possible method to increase viral titers in salmonid cells.

Regulation of host gene expression by gastrointestinal tract microbiota in Chinook Salmon (Oncorhynchus tshawytscha).

Differences in gut microbiome composition are linked with health, disease and ultimately host fitness; however, the molecular mechanisms underlying that relationship are not well characterized. Here, we modified the fish gut microbiota using antibiotic and probiotic feed treatments to address the effect of host microbiome on gene expression patterns. Chinook salmon (Oncorhynchus tshawytscha) gut gene expression was evaluated using whole transcriptome sequencing (RNA-Seq) on hindgut mucosa samples from individuals treated with antibiotic, probiotic and control diets to determine differentially expressed (DE) host genes. Fifty DE host genes were selected for further characterization using nanofluidic qPCR chips. We used 16S rRNA gene metabarcoding to characterize the rearing water and host gut microbiome (bacterial) communities. Daily administration of antibiotics and probiotics resulted in significant changes in fish gut and aquatic microbiota as well as more than 100 DE genes in the antibiotic and probiotic treatment fish, relative to healthy controls. Normal microbiota depletion by antibiotics mostly led to downregulation of different aspects of immunity and upregulation of apoptotic process. In the probiotic treatment, genes related to post-translation modification and inflammatory responses were up-regulated relative to controls. Our qPCR results revealed significant effects of treatment (antibiotic and probiotic) on rabep2, aifm3, manf, prmt3 gene transcription. Moreover, we found significant associations between members of Lactobacillaceae and Bifidobacteriaceae with host gene expression patterns. Overall, our analysis showed that the microbiota had significant impacts on many host signalling pathways, specifically targeting immune, developmental and metabolic processes. Our characterization of some of the molecular mechanisms involved in microbiome-host interactions will help develop new strategies for preventing/ treating microbiome disruption-related diseases.

Divergent RNA viruses infecting sea lice, major ectoparasites of fish.

Sea lice, the major ectoparasites of fish, have significant economic impacts on wild and farmed finfish, and have been implicated in the decline of wild salmon populations. As blood-feeding arthropods, sea lice may also be reservoirs for viruses infecting fish. However, except for two groups of negative-strand RNA viruses within the order Mononegavirales, nothing is known about viruses of sea lice. Here, we used transcriptomic data from three key species of sea lice (Lepeophtheirus salmonis, Caligus clemensi, and Caligus rogercresseyi) to identify 32 previously unknown RNA viruses. The viruses encompassed all the existing phyla of RNA viruses, with many placed in deeply branching lineages that likely represent new families and genera. Importantly, the presence of canonical virus-derived small interfering RNAs (viRNAs) indicates that most of these viruses infect sea lice, even though in some cases their closest classified relatives are only known to infect plants or fungi. We also identified both viRNAs and PIWI-interacting RNAs (piRNAs) from sequences of a bunya-like and two qin-like viruses in C. rogercresseyi. Our analyses showed that most of the viruses found in C. rogercresseyi occurred in multiple life stages, spanning from planktonic to parasitic stages. Phylogenetic analysis revealed that many of the viruses infecting sea lice were closely related to those that infect a wide array of eukaryotes with which arthropods associate, including fungi and parasitic tapeworms, implying that over evolutionary time there has been cross-phylum and cross-kingdom switching of viruses between arthropods and other eukaryotes. Overall, this study greatly expands our view of virus diversity in crustaceans, identifies viruses that infect and replicate in sea lice, and provides evidence that over evolutionary time, viruses have switched between arthropods and eukaryotic hosts in other phyla and kingdoms.

The effects of host quantitative genetic architecture on the gut microbiota composition of Chinook salmon (Oncorhynchus tshawytscha).

The microbiota consists of microbes living in or on an organism and has been implicated in host health and function. Environmental and host-related factors were shown to shape host microbiota composition and diversity in many fish species, but the role of host quantitative architecture across populations and among families within a population is not fully characterized. Here, Chinook salmon were used to determine if inter-population differences and additive genetic variation within populations influenced the gut microbiota diversity and composition. Specifically, hybrid stocks of Chinook salmon were created by crossing males from eight populations with eggs from an inbred line created from self-fertilized hermaphrodite salmon. Based on high-throughput sequencing of the 16S rRNA gene, significant gut microbial community diversity and composition differences were found among the hybrid stocks. Furthermore, additive genetic variance components varied among hybrid stocks, indicative of population-specific heritability patterns, suggesting the potential to select for specific gut microbiota composition for aquaculture purposes. Determining the role of host genetics in shaping their gut microbiota has important implications for predicting population responses to environmental changes and will thus impact conservation efforts for declining populations of Chinook salmon.

Molecular evidence for stress, inflammation and structural changes in non-specific ulcers in skin of farmed Chinook salmon (Oncorhynchus tshawytscha).

Fish skin is critical to physical defence against pathogens and there is a need to understand the physiological processes impacting ulcers and their healing. Ulcers have been reported in farmed Chinook salmon in New Zealand. This study investigated stress, immune and structural gene expression in farmed Chinook salmon skin with and without ulcers from two sites in New Zealand sampled from February (higher temperature, late summer) to May (lower temperature, late autumn). Skin samples taken adjacent to non-specific ulcers in May and control fish in February demonstrated upregulation of heat shock protein 70 relative to control fish in May. Anterior gradient 2 expression was upregulated in fish with ulcers relative to control fish (both February and May), suggesting increased mucous cell activity. Based on the results of this study, fish with non-specific ulcers showed evidence of stress, inflammation, re-epithelisation, and delayed healing near the ulcer site, elucidating the importance of these processes in the pathogenesis of non-specific ulcers in farmed chinook salmon.

Functional validation of transposable element-derived cis-regulatory elements in Atlantic salmon.

Transposable elements (TEs) are hypothesized to play important roles in shaping genome evolution following whole-genome duplications (WGDs), including rewiring of gene regulation. In a recent analysis, duplicate gene copies that had evolved higher expression in liver following the salmonid WGD ∼100 million years ago were associated with higher numbers of predicted TE-derived cis-regulatory elements (TE-CREs). Yet, the ability of these TE-CREs to recruit transcription factors (TFs) in vivo and impact gene expression remains unknown. Here, we evaluated the gene-regulatory functions of 11 TEs using luciferase promoter reporter assays in Atlantic salmon (Salmo salar) primary liver cells. Canonical Tc1-Mariner elements from intronic regions showed no or small repressive effects on transcription. However, other TE-CREs upstream of transcriptional start sites increased expression significantly. Our results question the hypothesis that TEs in the Tc1-Mariner superfamily, which were extremely active following WGD in salmonids, had a major impact on regulatory rewiring of gene duplicates, but highlights the potential of other TEs in post-WGD rewiring of gene regulation in the Atlantic salmon genome.

Genomic variation across Chinook salmon populations reveals effects of a duplication on migration alleles and supports fine scale structure.

The distribution of ecotypic variation in natural populations is influenced by neutral and adaptive evolutionary forces that are challenging to disentangle. This study provides a high-resolution portrait of genomic variation in Chinook salmon (Oncorhynchus tshawytscha) with emphasis on a region of major effect for ecotypic variation in migration timing. With a filtered data set of ~13 million single nucleotide polymorphisms (SNPs) from low-coverage whole genome resequencing of 53 populations (3566 barcoded individuals), we contrasted patterns of genomic structure within and among major lineages and examined the extent of a selective sweep at a major effect region underlying migration timing (GREB1L/ROCK1). Neutral variation provided support for fine-scale structure of populations, while allele frequency variation in GREB1L/ROCK1 was highly correlated with mean return timing for early and late migrating populations within each of the lineages (r2  = .58-.95; p < .001). However, the extent of selection within the genomic region controlling migration timing was much narrower in one lineage (interior stream-type) compared to the other two major lineages, which corresponded to the breadth of phenotypic variation in migration timing observed among lineages. Evidence of a duplicated block within GREB1L/ROCK1 may be responsible for reduced recombination in this portion of the genome and contributes to phenotypic variation within and across lineages. Lastly, SNP positions across GREB1L/ROCK1 were assessed for their utility in discriminating migration timing among lineages, and we recommend multiple markers nearest the duplication to provide highest accuracy in conservation applications such as those that aim to protect early migrating Chinook salmon. These results highlight the need to investigate variation throughout the genome and the effects of structural variants on ecologically relevant phenotypic variation in natural species.