A transcriptome analysis focusing on splenic immune-related mciroRNAs of rainbow trout upon Aeromonas salmonicida subsp. salmonicida infection

Integrative analysis of miRNA and mRNA expression associated with the immune response in the intestine of rainbow trout (Oncorhynchus mykiss) infected with infectious hematopoietic necrosis virus

Rainbow trout (Oncorhynchus mykiss), an economically important cold-water fish cultured worldwide, suffers from infectious hematopoietic necrosis virus (IHNV) infection, resulting in huge financial losses. In order to understand the immune response of rainbow trout during virus infection, we explored trout intestine transcriptome profiles following IHNV challenge, and identified 3355 differentially expressed genes (DEGs) and 80 differentially expressed miRNAs (DEMs). Transcriptome analysis revealed numerous DEGs involved in immune responses, such as toll-like receptor 3 (TLR3), toll-like receptor 7/8 (TLR7/8), tripartite motif-containing 25 (TRIM25), DExH-Box helicase 58 (DHX58), interferon-induced with helicase C domain 1 (IFIH1), interferon regulatory factor 3 (IRF3/7), signal transducer and activator of transcription 1 (STAT1) and heat shock protein 90-alpha 1 (HSP90A1). Integrated analysis identified five key miRNAs (miR-19-y, miR-181-z, miR-203-y, miR-143-z and miR-206-y) targeting at least two important immune genes (TRIM25, DHX58, STAT1, TLR7/8 and HSP90A1). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that DEGs and target genes were significantly enriched in various immune-related terms including immune system process, binding, cell part and pathways of Toll-like receptor signalling, RIG-I-like receptor signalling, NOD-like receptor signalling, JAK-STAT signalling, PI3K-Akt signalling, NF-kappa B signalling, IL-17 signalling and AGE-RAGE signalling. In addition, protein-protein interaction networks (PPI) was used to display highly interactive DEG networks involving eight immune-related pathways. The expression trends of 12 DEGs and 10 DEMs were further verified by quantitative real-time PCR, which confirmed the reliability of the transcriptome sequencing results. This study expands our understanding of the immune response of rainbow trout infected with IHNV, and provides valuable resources for future studies on the immune molecular mechanism and disease resistance breeding.

Time-course transcriptome analyses of spleen in rainbow trout (Oncorhynchus mykiss) post-Flavobacterium psychrophilum infection

Flavobacterium psychrophilum, the etiological agent of bacterial coldwater disease and rainbow trout fry syndrome, causes considerable losses in salmonid aquaculture globally. Systemic F. psychrophilum infections in rainbow trout (Oncorhynchus mykiss) lead to a range of clinical signs, including ulcerative lesions in the skin and muscle and splenitis. Previous studies offered an integrative analysis of the skeletal muscle response to F. psychrophilum infection in rainbow trout. However, little is known about the molecular mechanism of immune response in the spleen, which is an important immune organ of rainbow trout. Here, we investigated the time-course splenic transcriptome profiles in uninfected rainbow trout (CK) and F. psychrophilum–infected rainbow trout at day 3 and day 7 (D3, D7) by RNA-seq analyses. Among the 7,170 differentially expressed genes (DEGs) in the three comparisons (D3 vs. CK, D7 vs. CK, D3 vs. D7), 1,286 DEGs showed consistent upregulation or downregulation at D3 and D7 and were associated with pattern recognition, acute-phase response, complement cascade, chemokine and cytokine signaling, and apoptosis. The Real time quantitative PCR (RT-qPCR) analysis of eight DEGs confirmed the accuracy of the RNA-Sequencing (RNA-seq) data. Our results reflected a general process from pathogen recognition to inflammatory cytokine generation and delineated a putative Toll-like receptor signaling pathway in rainbow trout spleen, following F. psychrophilum infection. Taken together, these results provide new insights into the molecular mechanism of the immune response to F. psychrophilum infection and are a valuable resource for future research on the prevention and control of bacterial coldwater disease during salmon culture.

Expression characteristics and interaction networks of microRNAs in spleen tissues of grass carp (Ctenopharyngodon idella)

The spleen is an important immune organ in fish. MicroRNAs (miRNAs) have been shown to play an important role in the regulation of immune function. However, miRNA expression profiles and their interaction networks associated with the postnatal late development of spleen tissue are still poorly understood in fish. The grass carp (Ctenopharyngodon idella) is an important economic aquaculture species in China. Here, two small RNA libraries were constructed from the spleen tissue of healthy grass carp at one-year-old and three-year-old. A total of 324 known conserved miRNAs and 9 novel miRNAs were identified by using bioinformatic analysis. Family analysis showed that 23 families such as let-7, mir-1, mir-10, mir-124, mir-8, mir-7, mir-9, and mir-153 were highly conserved between vertebrates and invertebrates. In addition, 14 families such as mir-459, mir-430, mir-462, mir-7147, mir-2187, and mir-722 were present only in fish. Expression analysis showed that the expression patterns of miRNAs in the spleen of one-year-old and three-year-old grass carp were highly consistent, and the percentage of miRNAs with TPM > 100 was above 39%. Twenty significant differentially expressed (SDE) miRNAs were identified. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that these SDE miRNAs were primarily involved in erythrocyte differentiation, lymphoid organ development, immune response, lipid metabolic process, the B cell receptor signaling pathway, the T cell receptor signaling pathway, and the PPAR signaling pathway. In addition, the following miRNA-mRNA interaction networks were constructed: immune and hematopoietic, cell proliferation and differentiation, and lipid metabolism. This study determined the miRNA transcriptome as well as miRNA-mRNA interaction networks in normal spleen tissue during the late development stages of grass carp. The results expand the number of known miRNAs in grass carp and are a valuable resource for better understanding the molecular biology of the spleen development in grass carp.

Transcriptional Signatures of Immune, Neural, and Endocrine Functions in the Brain and Kidney of Rainbow Trout (Oncorhynchus mykiss) in Response to Aeromonas salmonicida Infection

Rainbow trout (Oncorhynchus mykiss) serves as one of the most important commercial fish with an annual production of around 800,000 tonnes. However, infectious diseases, such as furunculosis caused by Aeromonas salmonicida infection, results in great economic loss in trout culture. The brain and kidney are two important organs associated with “sickness behaviors” and immunomodulation in response to disease. Therefore, we worked with 60 trout and investigated transcriptional responses and enrichment pathways between healthy and infected trout. We observed that furunculosis resulted in the activation of toll-like receptors with neuroinflammation and neural dysfunction in the brain, which might cause the “sickness behaviors” of infected trout including anorexia and lethargy. We also showed the salmonid-specific whole genome duplication contributed to duplicated colony stimulating factor 1 (csf-1) paralogs, which play an important role in modulating brain immunomodulation. Enrichment analyses of kidneys showed up-regulated immunomodulation and down-regulated neural functions, suggesting an immune-neural interaction between the brain and kidney. Moreover, the kidney endocrine network was activated in response to A. salmonicida infection, further convincing the communications between endocrine and immune systems in regulating internal homeostasis. Our study provided a foundation for pathophysiological responses of the brain and kidney in response to furunculosis and potentially offered a reference for generating disease-resistant trout strains.

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.

Identification and characterization of caspases genes in rainbow trout (Oncorhynchus mykiss) and their expression profiles after Aeromonas salmonicida and Vibrio anguillarum infection

Caspases are highly conserved cysteine-dependent aspartyl-specific proteases that play an important role in regulating cell death and inflammation. However, the caspase genes have not been systematically studied in rainbow trout (Oncorhynchus mykiss). Rainbow trout experienced 4 rounds (4R) of genome duplication in the evolutionary history. Thereby an increased numbers of paralogs are observed in trout, probably with more complicated gene functions. We identified 18 caspase genes in rainbow trout, including two inflammatory caspases (casp1a, casp1b), six apoptosis executioner caspases (casp3, casp3a1, casp3a2, casp3b, casp6, and casp7), nine apoptosis initiator caspases (casp2a, casp2b, casp8, casp9a, casp9b, casp10a, casp10b, casp20a, and casp20b) and one uncategorized caspase gene (casp17). To investigate the potentially physiological functions of caspase genes, we challenged the rainbow trout with Aeromonas salmonicida (A. salmonicida) and Vibrio anguillarum (V. anguillarum). Results showed that the CASP3-regulated intrinsic apoptosis was activated after A. salmonicida infection, while the CASP8 and CASP6-regulated extrinsic apoptosis exerted the greatest effect on trout challenged with V. anguillarum. In response to V. anguillarum infection, the data of RNA-Seq further showed the casp8 was tightly integrated with the significantly enriched Gene Ontology terms and functional pathways, including apoptosis regulation, pathogen detection and immunomodulation. Our study provides a foundation for the physiological functions and regulatory network of the caspase genes in teleosts.

Molecular insight into regulation of miRNAs in the spleen of zebrafish (Danio rerio) upon pathogenic Streptococcus parauberis infection

MicroRNAs (miRNAs) constitute a group of small non-coding RNAs (~22 nucleotides) and one of their main functions is to regulate the immune responses. Gram-positive bacterium, Streptococcus parauberis is the main causative agent of "Streptococcosis" in wide range of fish species. In this study, we performed high throughput sequencing analysis to identify the miRNA profile against S. parauberis infection in the spleen of zebrafish (Danio rerio). Overall, 349 known and 151 novel miRNAs were discovered. Among them, 12 known miRNAs (dre-miR-34b, dre-miR-135a, dre-miR-200b-5p, dre-miR-146b, dre-miR-31, dre-miR-17a-3p, dre-miR-222a-3p, dre-miR-731, dre-miR-301b-3p and dre-miR-30a-3p) and 9 novel miRNAs were differentially expressed (DE) in the spleen of S. parauberis challenged zebrafish. The identified 12 DE miRNAs were predicted to regulate 721 target genes. We confirmed the miRNA expression results by validating selected known and novel DE miRNAs using qRT-PCR. Gene Ontology (GO), Kyoto Encyclopedia of Genes (KEGG) pathway analysis and miRNA-mRNA interactions implies that specific target genes of DE miRNAs are associated with immune responses. The enriched pathways included Toll-like receptor (TLR), C-type lectin, NOD-like receptor, and RIG-I-like receptor signaling pathways, etc. Especially, dre-miR-200b-5p, dre-miR-146b, dre-miR-731, dre-miR-222a-3p, and dre-miR-34b were able to target potential immune-related genes such as il10, irak1, traf6, hspa8 and ikbke upon S. parauberis challenge. Thus, overall results could lay a foundation to understand the underlying immune regulatory role of miRNAs in response to pathogenic S. parauberis infection in teleosts.

Integrated analysis of mRNA-miRNA expression in Tilapia infected with Tilapia lake virus (TiLV) and identifies primarily immuneresponse genes

We investigated differential gene expression in Tilapia infected with the Tilapia Lake virus (TiLV).We used high-throughput sequencing to identify mRNAs and miRNAs involved in TiLV infection progression We identified 25,359 differentially expressed genes that included 863 new genes. We identified 1770, 4142 and 4947 differently expressed genes comparing non-infected controls with 24 and 120 h infections and between the infected groups, respectively. These genes were enriched to 291 GO terms and 62 KEGG pathways and included immune system progress and virion genes. High-throughput miRNA sequencing identified 316 conserved miRNAs, 525 known miRNAs and 592 novel miRNAs. Furthermore, 138, 198 and 153 differently expressed miRNAs were found between the 3 groups listed above, respectively. Target prediction revealed numerous genes including erythropoietin isoform X2, double-stranded RNA-specific adenosine deaminase isoform X1, bone morphogenetic protein 4 and tapasin-related protein that are involved in immune responsiveness. Moreover, these target genes overlapped with differentially expressed mRNAs obtained from RNA-seq. These target genes were significantly enriched to GO terms and KEGG pathways including immune system progress, virion and Wnt signaling pathways. Expression patterns of differentially expressed mRNA and miRNAs were validated in 20 mRNA and 19 miRNAs by qRT-PCR. We also were able to construct a miRNA-mRNA target network that can further understand the molecular mechanisms on the pathogenesis of TiLV and guide future research in developing effective agents and strategies to combat TiLV infections in Tilapia.

Dietary Immunostimulant CpG Modulates MicroRNA Biomarkers Associated with Immune Responses in Atlantic Salmon (Salmo salar)

MicroRNAs (miRNAs) are key regulators in fish immune responses. However, no study has previously characterized the impact of polyriboinosinic polyribocytidylic acid (pIC) and formalin-killed typical Aeromonas salmonicida (ASAL) on miRNA expression in Atlantic salmon fed a commercial diet with and without immunostimulant CpG. To this end, first, we performed small RNA deep sequencing and qPCR analyses to identify and confirm pIC- and/or ASAL-responsive miRNAs in the head kidney of salmon fed a control diet. DESeq2 analyses identified 12 and 18 miRNAs differentially expressed in pIC and ASAL groups, respectively, compared to the controls. Fifteen of these miRNAs were studied by qPCR; nine remained significant by qPCR. Five miRNAs (miR-27d-1-2-5p, miR-29b-2-5p, miR-146a-5p, miR-146a-1-2-3p, miR-221-5p) were shown by qPCR to be significantly induced by both pIC and ASAL. Second, the effect of CpG-containing functional feed on miRNA expression was investigated by qPCR. In pre-injection samples, 6 of 15 miRNAs (e.g., miR-181a-5-3p, miR-462a-3p, miR-722-3p) had significantly lower expression in fish fed CpG diet than control diet. In contrast, several miRNAs (e.g., miR-146a-1-2-3p, miR-192a-5p, miR-194a-5p) in the PBS- and ASAL-injected groups had significantly higher expression in CpG-fed fish. Multivariate statistical analyses confirmed that the CpG diet had a greater impact on miRNA expression in ASAL-injected compared with pIC-injected fish. This study identified immune-relevant miRNA biomarkers that will be valuable in the development of diets to combat infectious diseases of salmon.