miRNAs Predicted to Regulate Host Anti-viral Gene Pathways in IPNV-Challenged Atlantic Salmon Fry Are Affected by Viral Load, and Associated With the Major IPN Resistance QTL Genotypes in Late Infection

Whole-Transcriptome Analysis Reveals Potential CeRNA Regulatory Mechanism in Takifugu rubripes against Cryptocaryon irritans Infection

Cryptocaryon irritans (C. irritans) is a proto-ciliate parasite that infects marine fishes, including the cultured species Takifugu rubripes (T. rubripes), causing disease and potential mortality. In host organisms, infection by parasites triggers an immune response that is modulated by regulatory elements including proteins and non-coding RNAs. In this study, the whole transcriptome RNA sequencing of T. rubripes gill tissue before and after infection with C. irritans was performed to reveal the competitive endogenous RNA (ceRNA) regulatory network. Histomorphology revealed gill segment swelling and parasitic invasion in the infected group. The analysis identified 18 differentially expressed miRNAs (DEMs), 214 lncRNAs (DELs), 2501 genes (DEGs), and 7 circRNAs (DECs) in the infected group. Gene Ontology (GO) enrichment analysis revealed that these genes were notably enriched in the Wnt signaling pathway and mTOR signaling pathway. The co-expression networks (lncRNA/circRNA-miRNA-mRNA) were constructed based on correlation analysis of the differentially expressed RNAs. Further analysis suggested that the LOC105418663-circ_0000361-fru-miR-204a-fzd3a ceRNA axis was potentially involved in the regulation of immune responses against C. irritans infection. Finally, the expression levels of DEG, DEL, and DEM were validated. This study reveals the regulatory mechanism of a candidate ceRNA network, providing insights into the potential mechanism of T. rubripes' infection with C. irritans.

Characterization of the impact of dietary immunostimulant CpG on the expression of mRNA biomarkers involved in the immune responses in Atlantic salmon (Salmo salar)

Infectious diseases have significantly impacted Atlantic salmon aquaculture worldwide. Modulating fish immunity with immunostimulant-containing functional feeds could be an effective strategy in mitigating disease problems. Previously, we characterized the impact of polyriboinosinic polyribocytidylic acid (pIC) and formalin-killed typical Aeromonas salmonicida bacterin on miRNA expression in Atlantic salmon fed a commercial diet with and without immunostimulant CpG. A set of miRNA biomarkers of Atlantic salmon head kidney responding to pIC and/or bacterin immune stimulations was identified (Xue et al., 2019) [1]. Herein, we report a complementary qPCR study that investigated the impact of the pIC, bacterin and dietary CpG on the expression of immune-relevant mRNAs (n = 31) using the same samples as in the previous study (Xue et al., 2019) [1]. Twenty-six of these genes were predicted target transcripts of the pIC- and/or bacterin-responsive miRNAs identified in the earlier study. The current data showed that pIC and/or bacterin stimulations significantly modulated the majority of the qPCR-analyzed genes involved in various immune pathways. Some genes responded to both stimulations (e.g. tnfa, il10rb, ifng, irf9, cxcr3, campb) while others appeared to be stimulation specific [e.g. irf3, irf7a, il1r1, mxa, mapk3 (pIC only); clra (bacterin only)]. A. salmonicida bacterin stimulation produced a strong inflammatory response (e.g. higher expression of il1b, il8a and tnfa), while salmon stimulated with pIC showed robust interferon responses (both type I and II). Furthermore, the current data indicated significant down-regulation of immune-relevant transcripts (e.g. tlr9, irf5, il1r1, hsp90ab1, itgb2) by dietary immunostimulant CpG, especially among pre-injection and PBS-injected fish. Together with our prior miRNA study, the present research provided complementary information on Atlantic salmon anti-viral and anti-bacterial immune responses and on how dietary CpG may modulate these responses.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Expression Analysis of Moritella viscosa-Challenged Atlantic Salmon Identifies Disease-Responding Genes, MicroRNAs and Their Predicted Target Genes and Pathways

Moritella viscosa is a bacterial pathogen causing winter-ulcer disease in Atlantic salmon. The lesions on affected fish lead to increased mortality, decreased fish welfare, and inferior meat quality in farmed salmon. MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional regulation by guiding the miRNA-induced silencing complex to specific mRNA transcripts (target genes). The goal of this study was to identify miRNAs responding to Moritella viscosa in salmon by investigating miRNA expression in the head-kidney and the muscle/skin from lesion sites caused by the pathogen. Protein coding gene expression was investigated by microarray analysis in the same materials. Seventeen differentially expressed guide-miRNAs (gDE-miRNAs) were identified in the head-kidney, and thirty-nine in lesion sites, while the microarray analysis reproduced the differential expression signature of several thousand genes known as infection-responsive. In silico target prediction and enrichment analysis suggested that the gDE-miRNAs were predicted to target genes involved in immune responses, hemostasis, angiogenesis, stress responses, metabolism, cell growth, and apoptosis. The majority of the conserved gDE-miRNAs (e.g., miR-125, miR-132, miR-146, miR-152, miR-155, miR-223 and miR-2188) are known as infection-responsive in other vertebrates. Collectively, the findings indicate that gDE-miRNAs are important post-transcriptional gene regulators of the host response to bacterial infection.

Integrated analysis of immune parameters, miRNA-mRNA interaction, and immune genes expression in the liver of rainbow trout following infectious hematopoietic necrosis virus infection

Rainbow trout (Oncorhynchus mykiss) is an important economical cold-water fish worldwide. However, infection with infectious hematopoietic necrosis virus (IHNV) has severely restricted the development of aquaculture and caused huge economic losses. Currently, little is known about the immune defense mechanisms of rainbow trout against IHNV. In this study, we detected the changes of immune parameters over different post-infection periods (6-, 12-, 24-, 48-, 72-, 96-, 120-, and 144 hours post-infection (hpi)), mRNA and miRNA expression profiles under 48 hpi (T48L) compared to control (C48L), and key immune-related genes expression patterns in rainbow trout liver following IHNV challenge through biochemical methods, RNA-seq, and qRT-PCR, and the function of miR-330-y was verified by overexpression and silencing in vitro and in vivo. The results revealed that alkaline phosphatase (AKP), alanine aminotransferase (ALT), catalase (CAT), and total superoxide dismutase (T-SOD) activities, and lysozyme (LZM) content showed significant peaks at 48 hpi, whereas malondialdehyde (MDA) content and aspartate aminotransferase (AST) activity decreased continuously during infection, and acid phosphatase (ACP) activity varied slightly. From RNA-seq, a total of 6844 genes and 86 miRNAs were differentially expressed, and numerous immune-related differentially expressed genes (DEGs) involved in RIG-I-like receptor signaling pathway, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, cytokine-cytokine receptor interaction, and antigen processing and presentation were significantly upregulated in T48Lm group, including IFIH1, DHX58, MAVS, TRAF3, IRF3, IRF7, MX1, TLR3, TLR8, MYD88, NOD1, NOD2, IL-8, CXCR1, CD209, CD83, and TAP1. Integrated analysis identified seven miRNAs (miR-425-x, miR-185-x, miR-338-x, miR-330-y, miR-361-x, miR-505-y, and miR-191-x) that target at least three key immune-related DEGs. Expression analysis showed that IFIH1, DHX58, IRF3, IRF7, MX1, TLR3, TLR8, and MYD88 showed a marked increase after 24 hpi during infection. Further research confirmed TAP1 as one of the targets of miR-330-y, overexpression of miR-330-y with mimics or agomir significantly reduced the expression levels of TAP1, IRF3, and IFN, and the opposite effects were obtained by inhibitor. These results facilitate in-depth understanding of the immune mechanisms in rainbow trout against IHNV.

Differential Expression of miRNAs and Their Predicted Target Genes Indicates That Gene Expression in Atlantic Salmon Gill Is Post-Transcriptionally Regulated by miRNAs in the Parr-Smolt Transformation and Adaptation to Sea Water

Smoltification (parr-smolt transformation) is a complex developmental process consisting of developmental changes that lead to remodeling of the Atlantic salmon gill. Here, the expression changes of miRNAs and mRNAs were studied by small-RNA sequencing and microarray analysis, respectively, to identify miRNAs and their predicted targets associated with smoltification and subsequent sea water adaptation (SWA). In total, 18 guide miRNAs were identified as differentially expressed (gDE miRNAs). Hierarchical clustering analysis of expression changes divided these into one cluster of 13 gDE miRNAs with decreasing expression during smoltification and SWA that included the miRNA-146, miRNA-30 and miRNA-7132 families. Another smaller cluster that showed increasing expression consisted of miR-101a-3p, miR-193b-5p, miR-499a-5p, miR-727a-3p and miR-8159-5p. The gDE miRNAs were predicted to target 747 of the genes (DE mRNAs), showing expression changes in the microarray analysis. The predicted targets included genes encoding NKA-subunits, aquaporin-subunits, cystic fibrosis transmembrane conductance regulator and the solute carrier family. Furthermore, the predicted target genes were enriched in biological processes associated with smoltification and SWA (e.g., immune system, reactive oxygen species, stress response and extracellular matrix organization). Collectively, the results indicate that remodeling of the gill involves the post-transcriptional regulation of gene expression by the characterized gDE miRNAs.

Insights Into miRNA-mRNA Regulatory Mechanisms of Cold Adaptation in Gymnocypris eckloni: Ubiquitin-Mediated Proteolysis Is Pivotal for Adaptive Energy Metabolism

This study aimed to understand cold stress adaptations mechanism in fish. Thus, the transcriptional response to cold conditions in Gymnocypris eckloni was evaluated using RNA-seq and microRNA (miRNA)-seq analyses. Low-temperature (LT) group G. eckloni was cultivated outdoors in waters cooled to 2–4°C for 3 weeks, while individuals in the control temperature (CT) group were exposed to 14–16°C. Significantly different responses were observed in both mRNA and miRNA expression profiles, with more mRNAs (1,833 and 1,869 mRNAs were up- and downregulated, respectively) and fewer miRNAs (15 and 6 were up- and downregulated, respectively) observed in the LT group individuals relative to the CT group individuals. A miRNA-mRNA network involved in the regulation of G. eckloni responses to cold stress was constructed; this network included ubiquitin-mediated proteolysis, protein processing, and oxidative phosphorylation. These results provided new insights into mechanisms of cold tolerance by fish, including decreased metabolic activity in addition to proteolysis.

Dynamic mRNA and miRNA expression of the head during early development in bighead carp (Hypophthalmichthys nobilis)

Background

Head of fish species, an exquisitely complex anatomical system, is important not only for studying fish evolution and development, but also for economic values. Currently, although some studies have been made on fish growth and body shapes, very limited information is available on the molecular mechanism of head development.

Results

In this study, RNA sequencing (RNA–Seq) and small RNA sequencing (sRNA–Seq) technologies were used to conduct integrated analysis for the head of bighead carp at different development stages, including 1, 3, 5, 15 and 30 Dph (days post hatch). By RNA-Seq data, 26 pathways related to growth and bone formation were identified as the main physiological processes during early development. Coupling this to sRNA–Seq data, we picked out six key pathways that may be responsible for head development, namely ECM receptor interaction, TNF signaling pathway, osteoclast differentiation, PI3K–Akt signaling pathway, Neuroactive ligand–receptor interaction and Jak–STAT signaling pathway. Totally, 114 important candidate genes from the six pathways were obtained. Then we found the top 20 key genes according to the degree value by cytohubba, which regulated cell growth, skeletal formation and blood homeostasis, such as pik3ca, pik3r1, egfr, vegfa, igf1 and itga2b. Finally, we also acquired 19 key miRNAs playing multiple roles in the perfection of various tissues in the head (such as brain, eye and mouth) and mineralization of head bone system, such as let–7e, miR–142a–5p, miR–144–3p, miR–23a–3p and miR–223.

Conclusions

Results of this study will be informative for genetic mechanisms of head development and also provide potential candidate targets for the interaction regulation during early growth in bighead carp.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08387-x.

Characterization of miRNAs in Embryonic, Larval, and Adult Lumpfish Provides a Reference miRNAome for Cyclopterus lumpus

Simple Summary

Lumpfish (Cyclopterus lumpus) is an emergent aquaculture species, and its miRNA repertoire is still unknown. miRNAs are critical post-transcriptional modulators of teleost gene expression. Therefore, a lumpfish reference miRNAome was characterized by small RNA sequencing and miRDeep analysis of samples from different organs and developmental stages. The resulting miRNAome, an essential reference for future expression analyses, consists of 443 unique mature miRNAs from 391 conserved and eight novel miRNA genes. Enrichment of specific miRNAs in particular organs and developmental stages indicates that some conserved lumpfish miRNAs regulate organ and developmental stage-specific functions reported in other teleosts.

Abstract

MicroRNAs (miRNAs) are endogenous small RNA molecules involved in the post-transcriptional regulation of protein expression by binding to the mRNA of target genes. They are key regulators in teleost development, maintenance of tissue-specific functions, and immune responses. Lumpfish (Cyclopterus lumpus) is becoming an emergent aquaculture species as it has been utilized as a cleaner fish to biocontrol sea lice (e.g., Lepeophtheirus salmonis) infestation in the Atlantic Salmon (Salmo salar) aquaculture. The lumpfish miRNAs repertoire is unknown. This study identified and characterized miRNA encoding genes in lumpfish from three developmental stages (adult, embryos, and larvae). A total of 16 samples from six different adult lumpfish organs (spleen, liver, head kidney, brain, muscle, and gill), embryos, and larvae were individually small RNA sequenced. Altogether, 391 conserved miRNA precursor sequences (discovered in the majority of teleost fish species reported in miRbase), eight novel miRNA precursor sequences (so far only discovered in lumpfish), and 443 unique mature miRNAs were identified. Transcriptomics analysis suggested organ-specific and age-specific expression of miRNAs (e.g., miR-122-1-5p specific of the liver). Most of the miRNAs found in lumpfish are conserved in teleost and higher vertebrates, suggesting an essential and common role across teleost and higher vertebrates. This study is the first miRNA characterization of lumpfish that provides the reference miRNAome for future functional studies.

Differential Transcriptomic Response of Rainbow Trout to Infection with Two Strains of IPNV

The IPN virus (IPNV) causes a highly contagious disease that affects farmed salmonids. IPNV isolates have been phylogenetically classified into seven genogroups, of which two are present in Chile, genogroups 1 and 5. This study aimed to compare the transcriptomic response of rainbow trout fry challenged with two Chilean isolates of IPNV, RTTX (genogroup 1), and ALKA (genogroup 5). Tissue samples from challenged individuals and controls were taken at 1, 7, and 20 days post-challenge and analyzed by RNA-Seq. The results revealed that infection with RTTX elicited a greater modulation of the trout transcriptome compared to ALKA infection, generating a greater number of highly differentially expressed genes in relation to the control fish. Gene Ontology enrichment indicated that functions related to the inflammatory and immune responses were modulated in fish challenged with both isolates throughout the trial, but with different regulation patterns. On day 1 post challenge, these functions were activated in those challenged with ALKA, but suppressed in RTTX-challenged fish. These results suggest that rainbow trout exhibit a differential transcriptomic response to infection with the two genetically distinct IPNV isolates, especially at early times post-infection.

Comparison of Gene Expression Between Resistant and Susceptible Families Against VPAHPND and Identification of Biomarkers Used for Resistance Evaluation in Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) has caused a heavy loss to shrimp aquaculture since its outbreak. Vibrio parahaemolyticus (VP_AHPND) is regarded as one of the main pathogens that caused AHPND in the Pacific white shrimp Litopenaeus vannamei. In order to learn more about the mechanism of resistance to AHPND, the resistant and susceptible shrimp families were obtained through genetic breeding, and comparative transcriptome approach was used to analyze the gene expression patterns between resistant and susceptible families. A total of 95 families were subjected to VP_AHPND challenge test, and significant variations in the resistance of these families were observed. Three pairs of resistant and susceptible families were selected for transcriptome sequencing. A total of 489 differentially expressed genes (DEGs) that presented in at least two pairwise comparisons were screened, including 196 DEGs highly expressed in the susceptible families and 293 DEGs in the resistant families. Among these DEGs, 16 genes demonstrated significant difference in all three pairwise comparisons. Gene set enrichment analysis (GSEA) of all 27,331 expressed genes indicated that some energy metabolism processes were enriched in the resistant families, while signal transduction and immune system were enriched in the susceptible families. A total of 32 DEGs were further confirmed in the offspring of the detected families, among which 19 genes were successfully verified. The identified genes in this study will be useful for clarifying the genetic mechanism of shrimp resistance against Vibrio and will further provide molecular markers for evaluating the disease resistance of shrimp in the breeding program.

PUF60 of Japanese flounder is regulated by pol-miR-novel_395 and involved in pathogen infection, autophagy, and apoptosis

MicroRNAs (miRNAs) are evolutionary conserved, non-coding small RNAs that have been shown to regulate diverse biological processes including immunity. In a previous study, a novel miRNA of Japanese flounder (Paralichthys olivaceus), pol-miR-novel_395, was found to be responsive in expression to the infection of the bacterial pathogen Edwardsiella tarda. In the present study, we examined the regulation and immune effect of pol-miR-novel_395 and its target gene. We found that pol-miR-novel_395 expression was regulated by E. tarda and megalocytivirus, and pol-miR-novel_395 targeted the gene of PUF60 (poly (U)-binding-splicing factor 60 kDa) of flounder (named PoPUF60). Constitutive expression of PoPUF60 occurred in relatively high levels in the heart and liver of flounder. Bacterial infection upregulated PoPUF60 expression, whereas viral infection downregulated PoPUF60 expression. Interference with PoPUF60 expression or overexpression of pol-miR-novel_395 in flounder cells strongly potentiated E. tarda infection. Consistently, in vivo knockdown of PoPUF60 enhanced bacterial dissemination in the tissues of flounder but blocked viral replication, whereas in vivo overexpression of PoPUF60 inhibited bacterial dissemination but facilitated viral replication. Additionally, pol-miR-novel_395 and PoPUF60 were involved in the process of autophagy and apoptosis. Collectively, these results indicated that PoPUF60 and pol-miR-novel_395 play an important role in pathogen infection, autophagy, and apoptosis.

MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3'UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome

Complete 3′UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3′UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3′UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3′UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3′UTR motifs.

Transcriptome Profiling of Atlantic Salmon Adherent Head Kidney Leukocytes Reveals That Macrophages Are Selectively Enriched During Culture

The Atlantic salmon (Salmo salar) is an economically important fish, both in aquaculture and in the wild. In vertebrates, macrophages are some of the first cell types to respond to pathogen infection and disease. While macrophage biology has been characterized in mammals, less is known in fish. Our previous work identified changes in the morphology, phagocytic ability, and miRNA profile of Atlantic salmon adherent head kidney leukocytes (HKLs) from predominantly “monocyte-like” at Day 1 of in vitro culture to predominantly “macrophage-like” at Day 5 of culture. Therefore, to further characterize these two cell populations, we examined the mRNA transcriptome profile in Day 1 and Day 5 HKLs using a 44K oligonucleotide microarray. Large changes in the transcriptome were revealed, including changes in the expression of macrophage and immune-related transcripts (e.g. csf1r, arg1, tnfa, mx2), lipid-related transcripts (e.g. fasn, dhcr7, fabp6), and transcription factors involved in macrophage differentiation and function (e.g. klf2, klf9, irf7, irf8, stat1). The in silico target prediction analysis of differentially expressed genes (DEGs) using miRNAs known to change expression in Day 5 HKLs, followed by gene pathway enrichment analysis, supported that these miRNAs may be involved in macrophage maturation by targeting specific DEGs. Elucidating how immune cells, such as macrophages, develop and function is a key step in understanding the Atlantic salmon immune system. Overall, the results indicate that, without the addition of exogenous factors, the adherent HKL cell population differentiates in vitro to become macrophage-like.

Expression of ssa-miR-155 during ISAV infection in vitro: Putative role as a modulator of the immune response in Salmo salar

Multiple cellular components are involved in pathogen-host interaction during viral infection; in this context, the role of miRNAs have become highly relevant. We assessed the expression of selected miRNAs during an in vitro infection of a Salmo salar cell line with Infectious Salmon Anemia Virus (ISAV), the causative agent of a severe disease by the same name. Salmon orthologs for miRNAs that regulate antiviral responses were measured using RT-qPCR in an in vitro time-course assay. We observed a modulation of specific miRNAs expression, where ssa-miR-155-5p was differentially over-expressed. Using in silico analysis, we identified the putative mRNA targets for ssa-miR-155-5p, finding a high prevalence of hosts immune response-related genes; moreover, several mRNAs involved in the viral infective process were also identified as targets for this miRNA. Our results suggest a relevant role for miR-155-5p in Salmo salar during an ISAV infection as a regulator of the immune response to the virus.

A de novo Full-Length mRNA Transcriptome Generated From Hybrid-Corrected PacBio Long-Reads Improves the Transcript Annotation and Identifies Thousands of Novel Splice Variants in Atlantic Salmon

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

Modulation of hepatic miRNA expression in Atlantic salmon (Salmo salar) by family background and dietary fatty acid composition

This study finds significant differences in hepatic fatty acid composition between four groups of Atlantic salmon (Salmo salar) consisting of offspring from families selected for high and low capacities to express the delta 6 desaturase isomer b and fed diets with 10% or 75% fish oil. The results demonstrated that hepatic lipid metabolism was affected by experimental conditions (diet/family). The fatty acid composition in the four groups mirrored the differences in dietary composition, but it was also associated with the family groups. Small RNA sequencing followed by RT-qPCR identified 12 differentially expressed microRNAs (DE miRNAs), with expression associated with family groups (miR-146 family members, miR-200b, miR-214, miR-221, miR-125, miR-135, miR-137, miR_nov_1), diets (miR-203, miR-462) or both conditions. All the conserved DE miRNAs have been reported as associated with lipid metabolism in other vertebrates. In silico predictions revealed 37 lipid metabolism pathway genes, including desaturases, transcription factors and key enzymes in the synthesis pathways as putative targets (e.g., srebp-1 and 2, Δ6fad_b and c, hmdh, elovl4 and 5b, cdc42). RT-qPCR analysis of selected target genes showed expression changes that were associated with diet and with family groups (d5fad, d6fad_a, srebp-1). There was a reciprocal difference in the abundance of ssa-miR-203a-3p and srebp-1 in one group comparison, whereas other predicted targets did not reveal any evidence of being negatively regulated by degradation. More experimental studies are needed to validate and fully understand the predicted interactions and how the DE miRNAs may participate in the regulation of hepatic lipid metabolism.

Characterization of miRNAs in Extracellular Vesicles Released From Atlantic Salmon Monocyte-Like and Macrophage-Like Cells

Cell-derived extracellular vesicles (EVs) participate in cell-cell communication via transfer of molecular cargo including genetic material like miRNAs. In mammals, it has previously been established that EV-mediated transfer of miRNAs can alter the development or function of immune cells, such as macrophages. Our previous research revealed that Atlantic salmon head kidney leukocytes (HKLs) change their morphology, phagocytic ability and miRNA profile from primarily “monocyte-like” at Day 1 to primarily “macrophage-like” at Day 5 of culture. Therefore, we aimed to characterize the miRNA cargo packaged in EVs released from these two cell populations. We successfully isolated EVs from Atlantic salmon HKL culture supernatants using the established Vn96 peptide-based pull-down. Isolation was validated using transmission electron microscopy, nanoparticle tracking analysis, and Western blotting. RNA-sequencing identified 19 differentially enriched (DE) miRNAs packaged in Day 1 versus Day 5 EVs. Several of the highly abundant miRNAs, including those that were DE (e.g. ssa-miR-146a, ssa-miR-155 and ssa-miR-731), were previously identified as DE in HKLs and are associated with macrophage differentiation and immune response in other species. Interestingly, the abundance relative of the miRNAs in EVs, including the most abundant miRNA (ssa-miR-125b), was different than the miRNA abundance in HKLs, indicating selective packaging of miRNAs in EVs. Further study of the miRNA cargo in EVs derived from fish immune cells will be an important next step in identifying EV biomarkers useful for evaluating immune cell function, fish health, or response to disease.

Characterization of Differentially Expressed miRNAs and Their Predicted Target Transcripts during Smoltification and Adaptation to Seawater in Head Kidney of Atlantic Salmon

Smoltification and early seawater phase are critical developmental periods with physiological and biochemical changes in Atlantic salmon that facilitates survival in saltwater. MicroRNAs (miRNAs) are known to have important roles in development, but whether any miRNAs are involved in regulation of gene expression during smoltification and the adaption to seawater is largely unknown. Here, small RNA sequencing of materials from head kidney before, during smoltification and post seawater transfer were used to study expression dynamics of miRNAs, while microarray analysis was applied to study mRNA expression dynamics. Comparing all timepoints, 71 miRNAs and 2709 mRNAs were identified as differentially expressed (DE). Hierarchical clustering analysis of the DE miRNAs showed three major clusters with characteristic expression changes. Eighty-one DE mRNAs revealed negatively correlated expression patterns to DE miRNAs in Cluster I and III. Furthermore, 42 of these mRNAs were predicted as DE miRNA targets. Gene enrichment analysis of negatively correlated target genes showed they were enriched in gene ontology groups hormone biosynthesis, stress management, immune response, and ion transport. The results strongly indicate that post-transcriptional regulation of gene expression by miRNAs is important in smoltification and sea water adaption, and this study identifies several putative miRNA-target pairs for further functional studies.