Identification and characterization of miRNAs and lncRNAs of coho salmon (Oncorhynchus kisutch) in normal immune organs

Genome-Wide Identification and Involvement in Response to Biotic and Abiotic Stresses of lncRNAs in Turbot (Scophthalmus maximus)

Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

Regulatory mechanism of LncRNAs in gonadal differentiation of hermaphroditic fish, Monopterus albus

BACKGROUND

Monopterus albus is a hermaphroditic fish with sex reversal from ovaries to testes via the ovotestes in the process of gonadal development, but the molecular mechanism of the sex reversal was unknown.

METHODS

We produced transcriptomes containing mRNAs and lncRNAs in the crucial stages of the gonad, including the ovary, ovotestis and testis. The expression of the crucial lncRNAs and their target genes was detected using qRT‒PCR and in situ hybridization. The methylation level and activity of the lncRNA promoter were analysed by applying bisulfite sequencing PCR and dual-luciferase reporter assays, respectively.

RESULTS

This effort revealed that gonadal development was a dynamic expression change. Regulatory networks of lncRNAs and their target genes were constructed through integrated analysis of lncRNA and mRNA data. The expression and DNA methylation of the lncRNAs MSTRG.38036 and MSTRG.12998 and their target genes Psmβ8 and Ptk2β were detected in developing gonads and sex reversal gonads. The results showed that lncRNAs and their target genes exhibited consistent expression profiles and that the DNA methylation levels were negatively regulated lncRNA expression. Furthermore, we found that Ptk2β probably regulates cyp19a1 expression via the Ptk2β/EGFR/STAT3 pathway to reprogram sex differentiation.

CONCLUSIONS

This study provides novel insight from lncRNA to explore the potential molecular mechanism by which DNA methylation regulates lncRNA expression to facilitate target gene transcription to reprogram sex differentiation in M. albus, which will also enrich the sex differentiation mechanism of teleosts.

Effect of TRAF6-knockout on gene expression and lncRNA expression in Epithelioma papulosum cyprini (EPC) cells

TRAF6 is a key immune gene that plays a significant role in toll-like receptor signal transduction and activates downstream immune genes involved in antiviral immunity in fish. To explore the role of TRAF6 in Epithelioma papulosum cyprini (EPC) cells, we knocked out the TRAF6 gene using the Clustered Regularly Interspaced Short Palindromic Repeats-Cas9 (CRISPR-Cas9) technique and then analyzed the transcriptomes of the knockout cells. In this study, we identified that 232 transcripts were differentially expressed in naive cells. Using the pipeline, we identified 381 novel lncRNAs in EPC cells, 23 of which were differentially expressed. Gene Ontology enrichment analysis demonstrated that differentially expressed genes (DEG) are implicated in various immune processes, such as neutrophil chemotaxis and mitogen-activated protein kinase binding. In addition, the KEGG pathway analysis revealed enrichment in immune-related pathways (Interleukin-17 signaling pathway, cytokine-cytokine receptor interaction, and TNF signaling pathway). Furthermore, the target genes of the differentially expressed lncRNAs were implicated in the negative regulation of interleukin-6 and tumor necrosis factor production. These results indicate that lncRNAs and protein-coding genes participate in the regulation of immune and metabolic processes in fish.

The potential regulatory role of the lncRNA-miRNA-mRNA axis in teleost fish

Research over the past two decades has confirmed that noncoding RNAs (ncRNAs), which are abundant in cells from yeast to vertebrates, are no longer "junk" transcripts but functional regulators that can mediate various cellular and physiological processes. The dysregulation of ncRNAs is closely related to the imbalance of cellular homeostasis and the occurrence and development of various diseases. In mammals, ncRNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), have been shown to serve as biomarkers and intervention targets in growth, development, immunity, and disease progression. The regulatory functions of lncRNAs on gene expression are usually mediated by crosstalk with miRNAs. The most predominant mode of lncRNA-miRNA crosstalk is the lncRNA-miRNA-mRNA axis, in which lncRNAs act as competing endogenous RNAs (ceRNAs). Compared to mammals, little attention has been given to the role and mechanism of the lncRNA-miRNA-mRNA axis in teleost species. In this review, we provide current knowledge about the teleost lncRNA-miRNA-mRNA axis, focusing on its physiological and pathological regulation in growth and development, reproduction, skeletal muscle, immunity to bacterial and viral infections, and other stress-related immune responses. Herein, we also explored the potential application of the lncRNA-miRNA-mRNA axis in the aquaculture industry. These findings contribute to an enhanced understanding of ncRNA and ncRNA-ncRNA crosstalk in fish biology to improve aquaculture productivity, fish health and quality.

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.

Comparative transcriptome analysis of long non coding RNA (lncRNA) in RTG-2 cells infected by infectious hematopoietic necrosis virus

Infectious hematopoietic necrosis virus (IHNV) is the vital pathogen that has caused the great economic loss in salmonid fisheries. To date, there is limited information concerning the changes of lncRNAs in RTG-2 cells infected by IHNV. In this study, a comparative transcriptome analysis of lncRNAs was performed in RTG-2 cells with and without IHNV infection to determine their changes and the effects on IHNV infection. The results showed that IHNV infection significantly changed the expression levels of lncRNAs and mRNAs, including 3693 differentially expressed lncRNAs (DE-lncRNAs) and 3503 differentially expressed mRNAs (DE-mRNAs) respectively. These DE-lncRNAs and DE-mRNAs induced by IHNV were mostly associated with immune response, RNA processing, and viral diseases related pathways. Further analysis found that some DE-lncRNAs might participate in the regulation of extracellular matrix metabolism, apoptosis, lipid synthesis, autophagy, and immune responses referring to the functions of their target genes. Afterwards, 349 co-expression relationships were constructed by 223 DE-lncRNAs and 271 DE-mRNAs, of which LTCONS_00146935 was the pivotal node in the interaction networks, and was together with its target genes modulated the immune responses under the IHNV infection. RT-qPCR results showed that the changes of the selected immune-related DEGs were in consistent with the RNA-seq data, suggesting that the sequencing data was relatively reliable. In summary, this is the first study to determine the changes and interactions of lncRNA-mRNA in RTG-2 cells under the IHNV infection. The results provided the valuable information concerning the lncRNAs in salmonid fish, which will benefit for future study on uncovering the roles of lncRNAs-mRNAs during the viral infection.

Genome-wide integrated analysis reveals functions of lncRNA-miRNA-mRNA interactions in Atlantic salmon challenged by Aeromonas salmonicida

Aeromonas salmonicida (A. salmonicida) is a pathogenic bacterium that causes serious problems in the global Atlantic salmon aquaculture industry. In this study, we comprehensively analyzed the profiles of lncRNAs, miRNAs and mRNAs in gills of Atlantic salmon at high-dose A. salmonicida infection (3.06 × 10⁸ CFU/mL), low-dose A. salmonicida infection (3.06 × 10⁵ CFU/mL), and a PBS (100 μL) control. We identified 65 differentially expressed lncRNAs, 41 miRNAs, and 512 mRNAs between the control group and infection groups. Functional analysis showed that these genes were significantly enriched in the p53 signaling pathway, Wnt signaling pathway, mTOR signaling pathway, JAK-STAT signaling pathway, and Toll-like receptor signaling pathway. In addition, we predicted key genes in immune-related pathways and constructed a lncRNA-miRNA-mRNA network based on whole transcriptomic analysis. We further predicted three lncRNA-miRNA-mRNA axes as potential novel biomarkers in regulating the immune response of Atlantic salmon against A. salmonicida infection.

Temporal Gene Expression Signature of Plasma Extracellular Vesicles-MicroRNAs from Post-Smolt Coho Salmon Challenged with Piscirickettsia salmonis

Piscirickettsiosis is the most important bacterial disease in the Chilean salmon industry, which has borne major economic losses due to failure to control it. Cells use extracellular vesicles (EVs) as an inter-cellular communicators to deliver several factors (e.g., microRNAs) that may regulate the responses of other cells. However, there is limited knowledge about the identification and characterization of EV-miRNAs in salmonids or the effect of infections on these. In this study, Illumina sequencing technology was used to identify Coho salmon plasma EV-miRNAs upon Piscirickettsia salmonis infection at four different time points. A total of 118 novels and 188 known EV-miRNAs, including key immune teleost miRNAs families (e.g., miR-146, miR-122), were identified. A total of 245 EV-miRNAs were detected as differentially expressed (FDR < 5%) in terms of control, with a clear down-regulation pattern throughout the disease. KEGG enrichment results of EV-miRNAs target genes showed that they were grouped mainly in cellular, stress, inflammation and immune responses. Therefore, it is hypothesized that P. salmonis could potentially benefit from unbalanced modulation response of Coho salmon EV-miRNAs in order to promote a hyper-inflammatory and compromised immune response through the suppression of different key immune host miRNAs during the course of the infection, as indicated by the results of this study.

The emerging regulatory roles of noncoding RNAs in immune function of fish: MicroRNAs versus long noncoding RNAs

The genome could be considered as raw data expressed in proteins and various types of noncoding RNAs (ncRNAs). However, a large portion of the genome is dedicated to ncRNAs, which in turn represent a considerable amount of the transcriptome. ncRNAs are modulated on levels of type and amount whenever any physiological process occurs or as a response to external modulators. ncRNAs, typically forming complexes with other partners, are key molecules that influence diverse cellular processes. Based on the knowledge of mammalian biology, ncRNAs are known to regulate and control diverse trafficking pathways and cellular activities. Long noncoding RNAs (lncRNAs) notably have diverse and more regulatory roles than microRNAs. Expanding these studies on fish has derived the same conclusion with relevance to other species, including invertebrates, explored the potentials to harness such types of RNA to further understand the biology of such organisms, and opened gates for applying recent technologies, such as RNA interference and delivering micromolecules as microRNAs to living cells and possibly to target organs. These technologies should improve aquaculture productivity and fish health, as well as help understand fish biology.

Temporal genome-wide DNA methylation signature of post-smolt Pacific salmon challenged with Piscirickettsia salmonis

Piscirickettsiosis is the most important bacterial disease in the Chilean salmon industry, which has sorted several efforts to its control, generating enormous economic losses. Epigenetic alterations, such as DNA methylation, can play a relevant role in the modulation of the metazoans response to pathogens. Bacterial disease may activate global and local immune responses generating intricate responses with significant biological impact in the host. However, it is scarcely understood how bacterial infections influence fish epigenetic alterations. In the present study, we utilized Pacific salmon and Piscirickettsiosis as model, to gain understanding into the dynamics of DNA methylation among fish-bacterial infection interactions. A genome-wide analysis of DNA methylation patterns in female spleen tissue of Pacific salmon was achieved by reduced representation bisulphite sequencing from a time course design. We determined 2,251, 1,918, and 2,516 differentially methylated regions DMRs among infected and control Pacific salmon in 1 dpi, 5 dpi, and 15 dpi, respectively. The mean methylation difference per DMR among control and infected groups was of ~35%, with an oscillatory pattern of hypo, hyper, and hypomethylation across the disease. DMCs, among the control and infected group, showed that they were statistically enriched in intergenic regions and depleted in exons. Functional annotation of the DMR genes demonstrated three KEGG principal categories, associated directly with the host response to pathogens infections. Our results provide the first evidence of epigenetic variation in fish provoked by bacterial infection and demonstrate that this variation can be modulated across the disease.

Identification and characterization of miRNA expression profiles across five tissues in giant panda

microRNA (miRNA) is a small endogenous noncoding RNA molecule that plays multiple roles in regulating most biological processes. However, for China's national treasure giant panda, a world-famous rare and protected species, reports of its miRNA have been found only in blood and breast milk. To explore the miRNA expression differences between different giant panda tissues, here, we generated the miRNA profiles of five tissues (heart, liver, spleen, lung, and kidney) from four giant pandas with Illumina Hiseq 2500 platform, and filtered the differentially expressed miRNAs (DEmiRs) in each tissue, predicted the target genes of miRNA from each tissue based on the DEmiRs. Then, the GO and KEGG enrichment analysis were conducted using the target genes predicted from DEmiRs in each tissue. The RNA-seq generated an average of 0.718 GB base per sample. A total of 1,256 known miRNAs and 12 novel miRNAs were identified, and there were 215, 131, 185, 83, and 126 tissue-specific DEmiRs filtered in the heart, liver, spleen, lung, and kidney, respectively, including miR-1b-5p, miR-122-5p, miR-143, miR-126-5p, and miR-10b-5p, respectively. The predicted target genes, including MYL2, LRP5, MIF, CFD, and PEBP1 in the heart, liver, spleen, lung, and kidney, respectively, were closely associated with tissue-specific biological functions. The enrichment analysis results of target genes showed tissue-specific characteristics, such as the significantly enriched GO terms extracellular matrix in the heart and insulin-like growth factor binding in the liver. The miRNA profiles of the heart, liver, spleen, lung, and kidney of giant panda have been reported in this study, it reveals the miRNA expression differences between different tissues of the giant panda, and provides valuable genetic resources for the further related molecular genetic research of the rare and protected species giant panda and other mammals.