Functional Diets Modulate lncRNA-Coding RNAs and Gene Interactions in the Intestine of Rainbow Trout Oncorhynchus mykiss

Exploring the Integrated Role of miRNAs and lncRNAs in Regulating the Transcriptional Response to Amino Acids and Insulin-like Growth Factor 1 in Gilthead Sea Bream (Sparus aurata) Myoblasts

In this study, gilthead sea bream (Sparus aurata) fast muscle myoblasts were stimulated with two pro-growth treatments, amino acids (AA) and insulin-like growth factor 1 (Igf-1), to analyze the transcriptional response of mRNAs, microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) and to explore their possible regulatory network using bioinformatic approaches. AA had a higher impact on transcription (1795 mRNAs changed) compared to Igf-1 (385 mRNAs changed). Both treatments stimulated the transcription of mRNAs related to muscle differentiation (GO:0042692) and sarcomere (GO:0030017), while AA strongly stimulated DNA replication and cell division (GO:0007049). Both pro-growth treatments altered the transcription of over 100 miRNAs, including muscle-specific miRNAs (myomiRs), such as miR-133a/b, miR-206, miR-499, miR-1, and miR-27a. Among 111 detected lncRNAs (>1 FPKM), only 30 were significantly changed by AA and 11 by Igf-1. Eight lncRNAs exhibited strong negative correlations with several mRNAs, suggesting a possible regulation, while 30 lncRNAs showed strong correlations and interactions with several miRNAs, suggesting a role as sponges. This work is the first step in the identification of the ncRNAs network controlling muscle development and growth in gilthead sea bream, pointing out potential regulatory mechanisms in response to pro-growth signals.

Identification and characterization of long non-coding RNAs in intestinal immune regulation of largemouth bass, Micropterus salmoides, under acute heat stress

Long noncoding RNAs (lncRNAs), which are RNA molecules that do not code for proteins and have a length exceeding 200 base pairs, have been found to play a crucial role in regulating intestinal immunity. The high mortality of various fish species induced by high temperatures is known to be associated with enteritis. Our investigation demonstrated that acute heat stress was responsible for inducing fish enteritis. However, the specific lncRNAs involved this process remains unknown. In this current study, we utilized intestinal sequencing data from the largemouth bass species Micropterus salmoides under acute heat stress, resulting in a total of 347,351,492 clean reads obtained from six cDNA libraries. A total of 3399 novel lncRNA transcripts originating from 2488 distinct lncRNA genes were successfully identified. Consistent with previous findings in other fish species, these lncRNAs demonstrated comparatively shorter transcript lengths when compared to protein-coding genes. Furthermore, a total of 216 novel lncRNA exhibited differential expression (DE) in the intestine of largemouth bass, meeting the criteria of absolute log2 fold change exceeding 2 and a p-value below 0.05. Additionally, these DE-lncRNAs were found to regulate 210 neighboring genes in a cis-regulatory manner. An examination of GO/KEGG enrichment revealed a notable enrichment of immune regulation (p < 0.05) among these cis-genes, with lncRNA MSTRG.8573.1 playing a significant role in regulating the jak-stat signaling pathway during this process. This study presents a comprehensive inventory of novel DE-lncRNA implicated in the development of enteritis in largemouth bass under acute heat stress. These findings offer valuable insights for future investigations on the regulation of lncRNAs to mitigate heat stress-induced fish enteritis.

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.

Comprehensive identification and expression profiling of immune-related lncRNAs and their target genes in the intestine of turbot (Scophthalmus maximus L.) in response to Vibrio anguillarum infection

Long non-coding RNA (lncRNA) play vital regulatory roles in various biological processes. Intestine is one of the most sensitive organs to environmental and homeostatic disruptions for fish. However, systematic profiles of lncRNAs in the intestine of teleost in responses to pathogen infections is still limited. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection, resulted in dramatic economic loss. Hereinto, the intestinal tissues of turbot were sampled at 0 h, 2 h, 12 h, and 48 h following V. anguillarum infection. The histopathological analysis revealed that the pathological trauma was mainly present in intestinal tunica mucosal epithelium. After high-throughput sequencing and bioinformatic analysis, a total of 9722 lncRNAs and 21,194 mRNAs were obtained, and the average length and exon number of lncRNAs were both less than those of mRNAs. Among which, a set of 158 lncRNAs and 226 mRNAs were differentially expressed (DE-lncRNAs and DEGs) in turbot intestine at three time points, related to many immune-related genes such as complement, interleukin, chemokine, lysosome, and macrophage, indicating their potential critical roles in immune responses. In addition, 2803 and 1803 GO terms were enriched for DEGs and co-expressed target genes of DE-lncRNAs, respectively. Moreover, 127 and 50 KEGG pathways including cell adhesion molecules (CAMs), phagosome, JAK-STAT signaling pathway, cytokine-cytokine receptor interaction, and intestinal immune network for IgA production, were enriched for DEGs and co-expressed target genes of DE-lncRNAs, respectively. Finally, qRT-PCR was conducted to confirm the reliability of sequencing data. The present study will set the foundation for the future exploration of lncRNA functions in teleost in response to bacterial infection.

Full-length transcriptome and analysis of bmp-related genes in Platypharodon extremus

Platypharodon extremus is an endemic species on the Qinghai–Tibet Plateau. As a secondary protected species in China, the basic genomic information of this species has not yet been reported. Here, through third-generation sequencing, the full-length transcriptome of P. extremus was obtained. We identified 323,290 CCS sequences, and a total of 50,083 unigenes were extracted after correction with second-generation sequencing data and the removal of redundant reads. A total of 50,067 transcripts were annotated with the various databases. Based on the sequence information, three members in the bone morphogenetic proteins (bmps) family and their receptors, were identified. We found that the special structures of these proteins (zinc-dependent metalloproteinase domain, CUB domains, EGF-like domains and TGF-β domain) are highly conserved in fish and that they are closely evolutionarily related to the bmps and bmp receptors of Cyprinidae fishes. This is the first study to sequence the full-length transcriptome of P. extremus, which will help us to further understand its biology.

Transcriptome profiling of the early developmental stages in the giant mussel Choromytilus chorus exposed to delousing drugs

The giant mussel Choromytilus chorus is a marine bivalve commonly collected in central - southern Chile from fishery zones shared with the salmon industry. These economically relevant areas are also affected by the use of pesticides for controlling sea lice infestations in salmon aquaculture. Their main target is the sea louse Caligus rogercresseyi. However, other than some physiological impacts, the molecular effects of delousing drugs in non-target species such as C. chorus remain largely understudied. This study aimed to explore the transcriptome modulation of Trochophore and D larvae stages of C. chorus after exposure to azamethiphos and deltamethrin drugs. Herein, RNA-seq analyses and mRNA-lncRNAs molecular interactions were obtained. The most significant changes were found between different larval development stages exposed to delousing drugs. Notably, significant transcriptional variations were correlated with the drug concentrations tested. The biological processes involved in the development, such as cell movement and transcriptional activity, were mainly affected. Long non-coding RNAs (lncRNAs) were also identified in this species, and the transcription activity showed similar patterns with coding mRNAs. Most of the significantly expressed lncRNAs were associated with genes annotated to matrix metalloproteinases, collagenases, and transcription factors. This study suggests that exposure to azamethiphos or deltamethrin drugs can modulate the transcriptome signatures related to the early development of the giant mussel C. chorus.

Exploring Sea Lice Vaccines against Early Stages of Infestation in Atlantic Salmon (Salmo salar)

The sea louse Caligus rogercresseyi genome has opened the opportunity to apply the reverse vaccinology strategy for identifying antigens with potential effects on lice development and its application in sea lice control. This study aimed to explore the efficacy of three sea lice vaccines against the early stage of infestation, assessing the transcriptome modulation of immunized Atlantic salmon. Therein, three experimental groups of Salmo salar (Atlantic salmon) were vaccinated with the recombinant proteins: Peritrophin (prototype A), Cathepsin (prototype B), and the mix of them (prototype C), respectively. Sea lice infestation was evaluated during chalimus I-II, the early-infective stages attached at 7-days post infestation. In parallel, head kidney and skin tissue samples were taken for mRNA Illumina sequencing. Relative expression analyses of genes were conducted to identify immune responses, iron transport, and stress responses associated with the tested vaccines during the early stages of sea lice infection. The vaccine prototypes A, B, and C reduced the parasite burden by 24, 44, and 52% compared with the control group. In addition, the RNA-Seq analysis exhibited a prototype-dependent transcriptome modulation. The high expression differences were observed in genes associated with metal ion binding, molecular processes, and energy production. The findings suggest a balance between the host’s inflammatory response and metabolic process in vaccinated fish, increasing their transcriptional activity, which can alter the early host–parasite interactions. This study uncovers molecular responses produced by three vaccine prototypes at the early stages of infestation, providing new knowledge for sea lice control in the salmon aquaculture.

Comparative transcriptome analysis reveals immunoregulation mechanism of lncRNA-mRNA in gill and skin of large yellow croaker (Larimichthys crocea) in response to Cryptocaryon irritans infection

Background

Cryptocaryonosis caused by Cryptocaryon irritans is one of the major diseases of large yellow croaker (Larimichthys crocea), which lead to massive economic losses annually to the aquaculture industry of L. crocea. Although there have been some studies on the pathogenesis for cryptocaryonosis, little is known about the innate defense mechanism of different immune organs of large yellow croaker.

Results

In order to analyze the roles of long non-coding RNAs and genes specifically expressed between immune organs during the infection of C. irritans, in this study, by comparing transcriptome data from different tissues of L. crocea, we identified tissue-specific transcripts in the gills and skin, including 507 DE lncRNAs and 1592 DEGs identified in the gills, and 110 DE lncRNAs and 1160 DEGs identified in the skin. Furthermore, we constructed transcriptome co-expression profiles of L. crocea gill and skin, including 7,503 long noncoding RNAs (lncRNAs) and 23,172 protein-coding genes. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that the DEGs and the target genes of the DE lncRNAs in the gill were specifically enriched in several pathways related to immune such as HIF-1 signaling pathway. The target genes of DE lncRNAs and DEGs in the skin are specifically enriched in the complement and coagulation cascade pathways. Protein–protein interaction (PPI) network analysis identified 3 hub genes including NFKBIA, TNFAIP3 and CEBPB, and 5 important DE lncRNAs including MSTRG.24134.4, MSTRG.3038.5, MSTRG.27019.3, MSTRG.26559.1, and MSTRG.10983.1. The expression patterns of 6 randomly selected differentially expressed immune-related genes were validated using the quantitative real-time PCR method.

Conclusions

In short, our study is helpful to explore the potential interplay between lncRNAs and protein coding genes in different tissues of L. crocea post C. irritans and the molecular mechanism of pathogenesis for cryptocaryonosis.

Highlights

Skin and gills are important sources of pro-inflammatory molecules, and their gene expression patterns are tissue-specific after C. irritans infection.

15 DEGs and 5 DE lncRNAs were identified as hub regulatory elements after C. irritans infection

The HIF-1 signaling pathway and the complement and coagulation cascade pathway may be key tissue-specific regulatory pathways in gills and skin, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08431-w.

The Emerging Role of Long Non-Coding RNAs in Development and Function of Gilthead Sea Bream (Sparus aurata) Fast Skeletal Muscle

Long non-coding RNAs (lncRNAs) are an emerging group of ncRNAs that can modulate gene expression at the transcriptional or translational levels. In the present work, previously published transcriptomic data were used to identify lncRNAs expressed in gilthead sea bream skeletal muscle, and their transcription levels were studied under different physiological conditions. Two hundred and ninety lncRNAs were identified and, based on transcriptomic differences between juveniles and adults, a total of seven lncRNAs showed potential to be important for muscle development. Our data suggest that the downregulation of most of the studied lncRNAs might be linked to increased myoblast proliferation, while their upregulation might be necessary for differentiation. However, with these data, as it is not possible to propose a formal mechanism to explain their effect, bioinformatic analysis suggests two possible mechanisms. First, the lncRNAs may act as sponges of myoblast proliferation inducers microRNAs (miRNAs) such as miR-206, miR-208, and miR-133 (binding energy MEF < -25.0 kcal). Secondly, lncRNA20194 had a strong predicted interaction towards the myod1 mRNA (ndG = -0.17) that, based on the positive correlation between the two genes, might promote its function. Our study represents the first characterization of lncRNAs in gilthead sea bream fast skeletal muscle and provides evidence regarding their involvement in muscle development.

Differentially expressed lncRNAs involved in immune responses of Haliotis diversicolor and H. discus hannai challenged with Vibrio parahaemolyticus

Although many studies have shown that lncRNA, a non-coding RNA with a length of more than 200 bases, is involved in various biological functions, including the immune process, stress process, and cell development process. However, the function of lncRNA in abalone, especially in immunity, has been rarely studied. H. discus hannai and H. diversicolor are two main aquaculture abalone, and their growth is easily affected by the main pathogen Vibrio parahaemolyticus. Through rigorous screening procedures for transcripts in this study, we found that lncRNAs were 34,240, 23,022 in Haliotis diversicolor and H. discus hannai injected with V. parahaemolyticus, respectively. We also identified the unique and common lncRNAs and mRNAs of two abalone species for the first time; the shared lncRNAs and mRNAs in Haliotis diversicolor and H. discus hannai were 2352 and 13,165, respectively. Then gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the differentially expressed target genes of common and unique lncRNAs has shown that common lncRNAs could be widely involved in the biological processes of stress and cell development in both abalone species. In contrast, unique lncRNAs are linked to the Toll-like receptor, NF-kappaB signaling pathway of H. diversicolor, and pattern recognition receptors and lectins immune-related pathways of H. discus hannai. The co-expression network shows that some immune-related genes, such as INFK1, INFK2, CASP2, CASP8, IRAK1, lectin C, were closely related to lncRNAs. Further, we identified the targeted relationship between some immune-related genes and lncRNAs by qRT-PCR, through which we showed that the expression trend between targeted genes, such as INFK1 and Lnc7057, lectin C and Lnc6943, Lnc5637, and PLCG1 and Lnc1692, were consistent. In general, our results showed that lncRNA expression was induced in the two species of abalone after being infected with V. parahaemolyticus, and lncRNA was involved in the immune response of abalone by targeting coding genes.

Identification and characterization of immune-related lncRNAs and lncRNA-miRNA-mRNA networks of Paralichthys olivaceus involved in Vibrio anguillarum infection

BACKGROUND

Long non-coding RNAs (lncRNAs) structurally resemble mRNAs and exert crucial effects on host immune defense against pathogen infection. Japanese flounder (Paralichthys olivaceus) is an economically important marine fish susceptible to Vibrio anguillarum infection. To date, study on lncRNAs in flounder is scarce.

RESULTS

Here, we reported the first systematic identification and characterization of flounder lncRNAs induced by V. anguillarum infection at different time points. A total of 2,368 lncRNAs were identified, 414 of which were differentially expressed lncRNAs (DElncRNAs) that responded significantly to V. anguillarum infection. For these DElncRNAs, 3,990 target genes (named DETGs) and 42 target miRNAs (named DETmiRs) were identified based on integrated analyses of lncRNA-mRNA and lncRNA-miRNA expressions, respectively. The DETGs were enriched in a cohort of functional pathways associated with immunity. In addition to modulating mRNAs, 36 DElncRNAs were also found to act as competitive endogenous RNAs (ceRNAs) that regulate 37 DETGs through 16 DETmiRs. The DETmiRs, DElncRNAs, and DETGs formed ceRNA regulatory networks consisting of 114 interacting DElncRNAs-DETmiRs-DETGs trinities spanning 10 immune pathways.

CONCLUSIONS

This study provides a comprehensive picture of lncRNAs involved in V. anguillarum infection. The identified lncRNAs and ceRNA networks add new insights into the anti-bacterial immunity of flounder.

Transcriptome Profiling of Long Non-coding RNAs During the Atlantic Salmon Smoltification Process

For salmon aquaculture, one of the most critical phase is the parr-smolt transformation. Studies around this process have mainly focused on physiological changes and the Na⁺/K⁺-ATPase activity during the osmoregulatory activity. However, understanding how the salmon genome regulates the parr-smolt transformation, specifically the molecular mechanisms involved, remains uncovered. This study aimed to explore the transcriptional modulation of long non-coding RNAs (lncRNAs), as key molecular regulators, during the freshwater (FW) to seawater (SW) transfer in Atlantic salmon. Transcriptome sequencing was performed from gill samples of Atlantic salmon adapted from FW to SW through gradual salinity changes from 0 to 30 PSU. The results showed that most transcripts differently modulated were downregulated in all salinity conditions. Relevant biological processes were associated with growth, collagen formation, immune response, metabolism, and heme transport. Notably, 2864 putative lncRNAs were identified in Atlantic salmon gills differently expressed during fish smoltification. The highest number of lncRNAs differently modulated was observed at 30 PSU. Correlation expression analysis suggests putative regulatory roles of lncRNAs with smoltification-related genes. Herein, co-localization of Na⁺/K⁺-ATPase, growth hormone receptor, and thyroid hormone receptor genes with lncRNAs differentially expressed suggest putative regulatory mechanisms in the Atlantic salmon genome. The lncRNAs can be used as novel biomarkers for the fish smoltification process. Here, the lncRNA_145326 and lncRNA_18762 are putatively related to the parr-smolt transfer in Atlantic salmon. This study is the first description of lncRNAs with putative regulatory roles in Atlantic salmon during the SW adaptation.

Characterization of a novel lncRNA (SETD3-OT) in turbot (Scophthalmus maximus L.)

LncRNAs have been demonstrated to play pivotal roles in virous biological processes, especially the gene expression regulation, including transcriptional regulation, posttranscriptional control and epigenetic processes. However, most of the current studies of lncRNAs are still limited in mammalian species, the investigations of functional roles of lncRNAs in teleost species are still lacking. In current study, we identified a novel lncRNA (SETD3-OT) in turbot, with 2,504 bp full-length obtained by 5' and 3' RACE, located in turbot chromosome 17, ranged from 20,933,835 to 20,936,302 bp. In addition, 8 neighboring genes of SETD3-OT were identified within 100 kbp in genome location. From the annotation of the neighboring adjacent genes, SETD3-OT might involve in regulation of cell apoptosis and cycle, the immune cell development, and the immune response against infection, and its expression pattern is similar to majority of the neighboring genes following Aeromonas salmonicida challenge. Intriguingly, SETD3-OT showed significant high expression levels in mucosal surfaces (intestine, gill and skin), and was dramatically down-regulated in these mucosal tissues following Vibrio anguillarum challenge, especially in gill and skin. In addition, SETD3-OT was distributed in nucleus, it might regulate the neighboring genes in cis or in trans. Taken together, our results provide insights for lncRNA in fish innate immunity, further studies should be conduct to explore the detailed molecular mechanism of the gene regulation between SETD3-OT and its neighboring genes.

Long Noncoding RNA Signatures Induced by Toll-Like Receptor 7 and Type I Interferon Signaling in Activated Human Plasmacytoid Dendritic Cells

Long noncoding RNAs (lncRNAs) exhibit highly lineage-specific expression and act through diverse mechanisms to exert control over a wide range of cellular processes. lncRNAs can function as potent modulators of innate immune responses through control of transcriptional and posttranscriptional regulation of mRNA expression and processing. Recent studies have demonstrated that lncRNAs participate in the regulation of antiviral responses and autoimmune disease. Despite their emerging role as immune mediators, the mechanisms that govern lncRNA expression and function have only begun to be characterized. In this study, we explore the role of lncRNAs in human plasmacytoid dendritic cells (pDCs), which are critical sentinel sensors of viral infection and contribute to the development of autoimmune disease. Using genome-wide sequencing approaches, we dissect the contributions of Toll-like receptor 7 (TLR7) and type I interferon (IFN-I) in the regulation of coding and noncoding RNA expression in CAL-1 pDCs treated with R848 or IFNβ. Functional enrichment analysis reveals both the unique and synergistic roles of TLR7 and IFN-I signaling in the orchestration of pDC function. These observations were consistent with primary cell immune responses elicited by detection of viral infection. We identified and characterized the conditional TLR7- and IFN-I-dependent regulation of 588 lncRNAs. Dysregulation of these lncRNAs could significantly alter pDC maturation, IFN-I and inflammatory cytokine production, antigen presentation, costimulation or tolerance cues, turnover, or localization, all consequential events during viral infection or IFN-I-driven autoimmune diseases such as systemic lupus erythematosus. These findings demonstrate the differential responsiveness of lncRNAs to unique immune stimuli, uncover regulatory mechanisms of lncRNA expression, and reveal a novel and tractable platform for the study of lncRNA expression and function.

Differential Transcriptomic and Metabolomic Responses in the Liver of Nile Tilapia (Oreochromis niloticus) Exposed to Acute Ammonia

Ammonia is toxic to aquatic animal. Currently, only limited works were reported on the responses of aquatic animals after ammonia exposure using "omics" technologies. Tilapia suffers from the stress of ammonia-nitrogen during intensive recirculating aquaculture. Optimizing ammonia stress tolerance has become an important issue in tilapia breeding. The molecular and biochemical mechanisms of ammonia-nitrogen toxicity have not been understood comprehensively in tilapia yet. In this study, using RNA-seq and gas chromatograph system coupled with a Pegasus HT time-of-flight mass spectrometer (GC-TOF-MS) techniques, we investigated differential expressed genes (DEGs) and metabolomes in the liver at 6 h post-challenges (6 hpc) and 24 h post-challenges (24 hpc) under high concentration of ammonia-nitrogen treatment. We detected 2258 DEGs at 6 hpc and 315 DEGs at 24 hpc. Functional enrichment analysis indicated that DEGs were significantly associated with cholesterol biosynthesis, steroid and lipid metabolism, energy conservation, and mitochondrial tissue organization. Metabolomic analysis detected 31 and 36 metabolites showing significant responses to ammonia-nitrogen stress at 6 and 24 hpc, respectively. D-(Glycerol 1-phosphate), fumaric acid, and L-malic acid were found significantly down-regulated at both 6 and 24 hpc. The integrative analysis of transcriptomics and metabolomics suggested considerable alterations and precise control of gene expression at both physiological and molecular levels in response to the stress of ammonia-nitrogen in tilapia.

Transcriptomic response of rainbow trout (Oncorhynchus mykiss) skeletal muscle to Flavobacterium psychrophilum

Flavobacterium psychrophilum is the etiologic agent of rainbow trout fry syndrome (RTFS). This pathogen infects a wide variety of salmonid species during freshwater stages, causing significant losses in the aquaculture industry. Rainbow trout (Oncorhynchus mykiss) infected with F. psychrophilum, presents as the main external clinical sign ulcerative lesions and necrotic myositis in skeletal muscle. We previously reported the in vitro cytotoxic activity of F. psychrophilum on rainbow trout myoblast, however little is known about the molecular mechanisms underlying the in vivo pathogenesis in skeletal muscle. In this study, we examined the transcriptomic profiles of skeletal muscle tissue of rainbow trout intraperitoneally challenged with low infection dose of F. psychrophilum. Using high-throughput RNA-seq, we found that 233 transcripts were up-regulated, mostly associated to ubiquitin mediated proteolysis and apoptosis. Conversely, 189 transcripts were down-regulated, associated to skeletal muscle contraction. This molecular signature was consistent with creatine kinase activity in plasma and oxidative damage in skeletal muscle. Moreover, the increased caspase activity suggests as a whole skeletal muscle atrophy induced by F. psychrophilum. This study offers an integrative analysis of the skeletal muscle response to F. psychrophilum infection and reveals unknown aspects of its pathogenesis in rainbow trout.

Transcriptomic Analysis of Potential "lncRNA-mRNA" Interactions in Liver of the Marine Teleost Cynoglossus semilaevis Fed Diets With Different DHA/EPA Ratios

Long non-coding RNAs (lncRNA) have emerged as important regulators of lipid metabolism and have been shown to play multifaceted roles in controlling transcriptional gene regulation, but very little relevant information has been available in fish, especially in non-model fish species. With a feeding trial on a typical marine teleost tongue sole C. semilaevis followed by transcriptomic analysis, the present study investigated the possible involvement of lncRNA in hepatic mRNA expression in response to different levels of dietary DHA and EPA, which are two most important fatty acids for marine fish. An 80-day feeding trial was conducted in a flow-through seawater system, and in this trial three experimental diets differing basically in DHA/EPA ratio, i.e., 0.61 (D/E-0.61), 1.46 (D/E-1.46), and 2.75 (D/E-2.75), were randomly assigned to 9 tanks of experimental fish. A total of 124.04 G high quality genome-wide clean data about coding and non-coding transcripts was obtained in the analysis of hepatic transcriptome. Compared to diet D/E-0.61, D/E-1.46 up-regulated expression of 178 lncRNAs and 2629 mRNAs, and down-regulated that of 47 lncRNAs and 3059 mRNAs, while D/E-2.75 resulted in much less change in gene expression. The co-expression and co-localization analysis of differentially expressed (DE) lncRNA and mRNA among dietary groups were then conducted. The co-expressed DE lncRNA and mRNA were primarily enriched in GO terms such as Metabolic process, Intracellular organelle, Catalytic activity, and Oxidoreductase activity, as well as in KEGG pathways such as Ribosome and Oxidative phosphorylation. Overlap of co-expression and co-localization analysis, i.e., lncRNA–mRNA matches “XR_523541.1–solute carrier family 16, member 5 (slc16a5)” and “LNC_000285–bromodomain adjacent to zinc finger domain 2A (baz2a),” were observed in all inter-group comparisons, indicating that they might crucially mediate the effects of dietary DHA and EPA on hepatic gene expression in tongue sole. In conclusion, this was the first time in marine teleost to investigate the possible lncRNA–mRNA interactions in response to dietary fatty acids. The results provided novel knowledge of lncRNAs in non-model marine teleost, and will serve as important resources for future studies that further investigate the roles of lncRNAs in lipid metabolism of marine teleost.

Identifying a Long QTL Cluster Across chrLG18 Associated with Salt Tolerance in Tilapia Using GWAS and QTL-seq

Understanding the genetic mechanism of osmoregulation is important for the improvement of salt tolerance in tilapia. In our previous study, we have identified a major quantitative trait locus (QTL) region located at 23.0 Mb of chrLG18 in a Nile tilapia line by QTL-seq. However, the conservation of these QTLs in other tilapia populations or species is not clear. In this study, we successfully investigated the QTLs associated with salt tolerance in a mass cross population from the GIFT line of Nile tilapia (Oreochromis niloticus) using a ddRAD-seq-based genome-wide association study (GWAS) and in a full-sib family from the Malaysia red tilapia strain (Oreochromis spp) using QTL-seq. Our study confirmed the major QTL interval that is located at nearly 23.0 Mb of chrLG18 in Nile tilapia and revealed a long QTL cluster across chrLG18 controlling for the salt-tolerant trait in both red tilapia and Nile tilapia. This is the first GWAS analysis on salt tolerance in tilapia. Our finding provides important insights into the genetic architecture of salinity tolerance in tilapia and supplies a basis for fine mapping QTLs, marker-assisted selection, and further detailed functional analysis of the underlying genes for salt tolerance in tilapia.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

Non-coding RNAs (ncRNAs) are functional RNA molecules that are transcribed from DNA but not translated into proteins. ncRNAs function as key regulators of gene expression and chromatin modification. Recently, the functional role of ncRNAs in teleost fish has been extensively studied. Teleost fish are a highly diverse group among the vertebrate lineage. Fish are also important in terms of aquatic ecosystem, food production and human life, being the source of animal proteins worldwide and models of biomedical research. However, teleost fish always suffer from the invasion of infectious pathogens including viruses and bacteria, which has resulted in a tremendous economic loss to the fishing industry worldwide. Emerging evidence suggests that ncRNAs, especially miRNAs and lncRNAs, may serve as important regulators in cytokine and chemokine signaling, antigen presentation, complement and coagulation cascades, and T cell response in teleost fish. In this review, we summarize current knowledge and understanding of the roles of both miRNAs and lncRNAs in immune regulation in teleost fish. Molecular mechanism insights into the function of ncRNAs in fish immune response may contribute to the development of potential biomarkers and therapeutic targets for the prevention and treatment of fish diseases.

Comparative Transcriptome Analysis Reveals Growth-Related Genes in Juvenile Chinese Sea Cucumber, Russian Sea Cucumber, and Their Hybrids

Heterosis is important for sea cucumber breeding, but its molecular mechanism remains largely unexplored. In this study, parental lines of Apostichopus japonicus from Russia (R) and China (C) were used to construct hybrids (CR and RC) by reciprocal crossing. We examined the transcriptional profiles of the hybrids (CR and RC) and the purebreds (CC and RR) at different developmental times. A total of 60.27 Gb of clean data was obtained, and 176,649 unigenes were identified, of which 50,312 unigenes were annotated. A total of 414,536 SNPs were identified. A total of 7011 differentially expressed genes (DEGs) were obtained between the purebreds and hybrids at 45 days after fertilization (DAF), and a total of 8218 DEGs were obtained between the purebreds and hybrids at 75 DAF. In addition, a total of 7652 DEGs were obtained between 45 DAF and 75 DAF. The significant DEGs were mainly involved in the MAPK and FOXO signaling pathways, especially in the Ras-Raf-MEK1/2-ERK module, which may be a key regulator of development and growth in juvenile A. japonicus. In addition, we also identified key growth-related genes, such as fgfs, igfs, megfs and hgfs, which were upregulated in the hybrids (RC and CR); these genes may play important roles in heterosis in A. japonicus. Our study provides fundamental information on the molecular mechanisms underlying heterosis in sea cucumber and might suggest strategies for the selection of rapidly growing strains of sea cucumber in aquaculture.

Expression of Antisense Long Noncoding RNAs as Potential Regulators in Rainbow Trout with Different Tolerance to Plant-Based Diets

Reformulation of aquafeeds in salmonid diets to include more plant proteins is critical for sustainable aquaculture. However, increasing plant proteins can lead to stunted growth and enteritis. Toward an understanding of the regulatory mechanisms behind plant protein utilization, directional RNA sequencing of liver tissues from a rainbow trout strain selected for growth on an all plant-protein diet and a control strain, both fed a plant diet for 12 weeks, were utilized to construct long noncoding RNAs. Antisense long noncoding RNAs were selected for differential expression and functional analyses since they have been shown to have regulatory actions within a genome. A total of 142 unique antisense long noncoding RNAs were differentially expressed between strains, 60 of which could be mapped to a gene. Genes underlying these noncoding RNAs are indicated in lipid metabolism and immunity. Six noncoding transcripts were also found to overlap with differentially expressed protein-coding genes, all of which were co-expressed. Associating variation in regulatory elements between rainbow trout strains with differing tolerance to plant-protein diets will assist in future studies toward increased gains throughout carnivorous aquaculture.