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

Transcriptome Analysis of Stinging Catfish (Heteropneustes fossilis) Kidney Reveals its Role During Adaptation to Hypertonic Environment

The facultative air-breather stinging catfish Heteropneustes fossilis can easily adapt to anisotonic environment up to a level of 6 ppt. The teleost kidney is an essential osmoregulatory organ and plays a vital role in immune system maintenance. We took kidney samples from control and 100 mM NaCl-treated stinging catfish to study the molecular responses to salinity stress by looking at the DEGs. We used high-throughput sequencing followed by de novo transcriptome assembly to analyze the transcriptome of the fish. A total of 89,813 unigene transcripts with 47.28% GC content and 901 bp average length were produced during the final contig assembly. A BLAST search revealed that 45,855 unigenes had a considerable amount of similarity to sequences found in the Swiss-Prot and Ref-Seq databases. Comparative transcriptome analysis revealed that 1660 genes were differentially expressed (693 were upregulated and 967 downregulated) in the kidneys of fish treated with NaCl compared to control kidneys. Several GO pathway terms, including "inflammatory response" and "neutrophil chemotaxis" in the biological process category, "cytoplasm" and "integral component of membrane" in the cellular component category, and "metal ion binding" and "DNA binding," in the molecular function category, were found to be significantly over-represented by enrichment and functional analysis of the DEGs. In addition, there was an overrepresentation of KEGG pathways like "chemokine signaling pathway," "cytokine-cytokine receptor interaction," and "metabolic pathways." The current work demonstrates that tolerance to hypertonic environmental stress necessitates the activation of many molecular pathways involved in sustaining the physiological systems in air-breathing stinging catfish. Adaptation is achieved through the enrichment of genes related to metabolic pathways, immune system response, stress response such as apoptosis and protein folding, and transporter of molecules involved in ionic and osmoregulation.

Rainbow Trout (Oncorhynchus mykiss) Pre-Smolts Treated with 11-Deoxycorticosterone Regulate Liver Carbohydrate Metabolism and Gill Osmoregulation

Smoltification is stressful for salmonids, and cortisol is one of the central endocrine regulators for seawater adaptation. It has been established that cortisol plays both mineralocorticoid and glucocorticoid functions by MR and GR, respectively, since the aldosterone hormone is absent. Recently, investigations have proposed that the 11-deoxycorticosterone (DOC) mineralocorticoid precursor might support cortisol effects, but this mechanism remains unclear. Hence, we assessed the early effects of DOC on rainbow trout pre-smolts, the key smoltification stage, via metabolic and transcriptomic approaches. Thirty-six juveniles (~120 g) were treated for 3 h with DOC (1 mg/kg) and/or mineralocorticoid (eplerenone) or glucocorticoid (mifepristone) receptor antagonists (n = 6 for each group). DOC decreased plasma glucose and pyruvate and increased phosphate and liver glycogen. DOC also downregulated carbohydrate metabolism-related genes in the liver. Finally, gill RNA-seq analysis presented 1660 differentially expressed transcripts in DOC versus vehicle, 1022 for eplerenone + DOC versus DOC and 3324 for mifepristone + DOC versus DOC. The enrichment analysis mainly revealed the upregulation of ion transmembrane transport and carbohydrate metabolism and the downregulation of stress and innate immune responses. This suggests a significant role of DOC in liver carbohydrate metabolism and gill osmoregulation of pre-smolts through both receptors. Hence, this could contribute to improving animal welfare monitoring during smoltification by featuring novel and potential biomarkers.

Transcriptome Profiling of Gene Expression in Atlantic Salmon (Salmo salar) at Early Stage of Development

For Atlantic salmon development, the most critical phase is the early development stage from egg to fry through alevin. However, the studies investigating the early development of Atlantic salmon based on RNA-seq are scarce and focus only on one stage of development. Therefore, using the RNA-seq technology, the assessment of different gene expressions of various early development stages (egg, alevin, and fry) was performed on a global scale. Over 22 GB of clean data was generated from 9 libraries with three replicates for each stage with over 90% mapping efficiency. A total of 5534 genes were differentially expressed, among which 19, 606, and 826 genes were specifically expressed in each stage, respectively. The transcriptome analysis showed that the number of differentially expressed genes (DEGs) increased as the Atlantic salmon progressed in development from egg to fry stage. In addition, gene ontology enrichment demonstrated that egg and alevin stages are characterized by upregulation of genes involved in spinal cord development, neuron projection morphogenesis, axonogenesis, and cytoplasmic translation. At the fry stage, upregulated genes were enriched in the muscle development process (muscle cell development, striated muscle cell differentiation, and muscle tissue development), immune system (defense response and canonical NF-kappaB signal transduction), as well as epidermis development. These results suggest that the early development of Atlantic salmon is characterized by a dynamic shift in gene expression and DEGs between different stages, which provided a solid foundation for the investigation of Atlantic salmon development.

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.

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.

Whole Transcriptome Analysis Reveals the Global Molecular Responses of mRNAs, lncRNAs, miRNAs, circRNAs, and Their ceRNA Networks to Salinity Stress in Hong Kong Oysters, Crassostrea hongkongensis

The Hong Kong oyster, Crassostrea hongkongensis, is an estuarine bivalve with remarkable commercial value in South China, and the increase of salinity in estuaries during the dry season has posed a major threat to the oyster farming. To explore the global transcriptional response to salinity stress, a whole-transcriptome analysis was performed with the gills of oysters in 6‰, 18‰, and 30‰ filtered seawater. Overall, 2243, 194, 371, and 167 differentially expressed mRNAs (DEmRNAs), differentially expressed long non-coding RNAs (DElncRNAs), differentially expressed circular RNAs (DEcircRNAs), and differentially expressed microRNAs (DEmiRNAs) were identified, respectively. Based on GO enrichment and KEGG pathway analysis, these important DEmRNAs, DElncRNAs, DEcircRNAs, and DEmiRNAs were predicted to be mainly involved in amino acids metabolism, microtubule movement, and immune defense. This demonstrated the complexity of dynamic transcriptomic profiles of C. hongkongensis in response to salinity fluctuation. The regulatory relationships of DEmiRNAs-DEmRNAs, DElncRNAs-DEmiRNAs, and DEcircRNAs-DEmiRNAs were also predicted, and finally, a circRNA-associated competing endogenous RNA (ceRNA) network was constructed, consisting of six DEcircRNAs, eight DEmiRNAs, and five DEmRNAs. The key roles of taurine and hypotaurine metabolism and phenylalanine metabolism were highlighted in this ceRNA network, which was consistent with the major contribution of free amino acids to intracellular osmolality and cell volume regulation. Collectively, this study provides comprehensive data, contributing to the exploration of coding and non-coding RNAs in C. hongkongensis salinity response. The results would benefit the understanding of the response mechanism of bivalves against salinity fluctuation, and provide clues for genetic improvement of C. hongkongensis with hyper-salinity tolerance.

Transcriptome Signatures of Atlantic Salmon-Resistant Phenotypes against Sea Lice Infestation Are Associated with Tissue Repair

Salmon aquaculture is constantly threatened by pathogens that impact fish health, welfare, and productivity, including the sea louse Caligus rogercresseyi. This marine ectoparasite is mainly controlled through delousing drug treatments that have lost efficacy. Therein, strategies such as salmon breeding selection represent a sustainable alternative to produce fish with resistance to sea lice. This study explored the whole-transcriptome changes in Atlantic salmon families with contrasting resistance phenotypes against lice infestation. In total, 121 Atlantic salmon families were challenged with 35 copepodites per fish and ranked after 14 infestation days. Skin and head kidney tissue from the top two lowest (R) and highest (S) infested families were sequenced by the Illumina platform. Genome-scale transcriptome analysis showed different expression profiles between the phenotypes. Significant differences in chromosome modulation between the R and S families were observed in skin tissue. Notably, the upregulation of genes associated with tissue repairs, such as collagen and myosin, was found in R families. Furthermore, skin tissue of resistant families showed the highest number of genes associated with molecular functions such as ion binding, transferase, and cytokine activity, compared with the susceptible. Interestingly, lncRNAs differentially modulated in the R/S families are located near genes associated with immune response, which are upregulated in the R family. Finally, SNPs variations were identified in both salmon families, where the resistant ones showed the highest number of SNPs variations. Remarkably, among the genes with SPNs, genes associated with the tissue repair process were identified. This study reported Atlantic salmon chromosome regions exclusively expressed in R or S Atlantic salmon families' phenotypes. Furthermore, due to the presence of SNPs and high expression of tissue repair genes in the resistant families, it is possible to suggest mucosal immune activation associated with the Atlantic salmon resistance to sea louse infestation.

Comprehensive Relationship Analysis of the Long Noncoding RNAs (lncRNAs) and the Target mRNAs in Response to the Infection of Edwardsiella anguillarum in European eel (Anguilla anguilla) Inoculated with Freund's Adjuvant

Freund's complete adjuvant (FCA) and incomplete adjuvant (FIA), generally applied in subunit fishery vaccine, have not been explored on the molecular mechanism of the non-specific immune enhancement. As long noncoding RNAs (lncRNAs) play vital regulating roles in various biological activities, in this study, we examined the genome-wide expression of transcripts in the liver of European eel (Anguilla anguilla, Aa) inoculated with FCA and FIA (FCIA) to elucidate the regulators of lncRNAs in the process of Edwardsiella anguillarum (Ea) infection and Aa anti-Ea infection using strand-specific RNA-seq. After eels were challenged by Ea at 28 days post the first inoculation (dpi), compared to the control uninfected eels (Li group), the control infected eels (Con_Li group) showed severe bleeding, hepatocyte atrophy, and thrombi formed in the hepatic vessels of the liver, although eels inoculated with FCIA (FCIA_Li group) also formed slight thrombi in the hepatic vessels. Compared to the FCIA_Li group, there was about 10 times colony-forming unit (cfu) in the Con_Li group per 100 μg liver tissue, and the relative percent survival (RPS) of eels was 50% in FCIA_Li vs Con_Li. Using high-throughput transcriptomics, differential expressed genes (DEGs) and transcripts were identified and the results were verified using fluorescence real-time polymerase chain reaction (qRT-PCR). Interactions between the differential expressed lncRNAs (DE-lncRNAs) and the target DEGs were explored using Cytoscape according to their co-expression and co-location relationship. We found 13,499 lncRNAs (10,176 annotated and 3423 novel lncRNAs) between 3 comparisons of Con_Li vs Li, FCIA_Li vs Li, and FCIA_Li vs Con_Li, of which 111, 110, and 129 DE-lncRNAs were ascertained. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs targeted by DE-lncRNAs revealed these DEGs mainly involved in single-organism cellular process in BP, membrane in CC and binding in MF, and KEGG pathways showed that the target DEGs in co-expression and co-location enriched in cell adhesion molecules. Finally, 118 DE-lncRNAs target 1161 DEGs were involved in an interaction network of 8474 co-expression and 333 co-location-related links, of which 16 DE-lncRNAs play vital roles in anti-Ea infection. Taken together, the interaction networks revealed that DE-lncRNAs underlies the process of Ea infection and Aa anti-Ea infection.

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.

Whole-Genome Transcript Expression Profiling Reveals Novel Insights into Transposon Genes and Non-Coding RNAs during Atlantic Salmon Seawater Adaptation

The growing amount of genome information and transcriptomes data available allows for a better understanding of biological processes. However, analysis of complex transcriptomic experimental designs involving different conditions, tissues, or times is relevant. This study proposes a novel approach to analyze complex data sets combining transcriptomes and miRNAs at the chromosome-level genome. Atlantic salmon smolts were transferred to seawater under two strategies: (i) fish group exposed to gradual salinity changes (GSC) and (ii) fish group exposed to a salinity shock (SS). Gills, intestine, and head kidney samples were used for total RNA extraction, followed by mRNA and small RNA illumina sequencing. Different expression patterns among the tissues and treatments were observed through a whole-genome transcriptomic approach. Chromosome regions highly expressed between experimental conditions included a great abundance of transposable elements. In addition, differential expression analysis showed a greater number of transcripts modulated in response to SS in gills and head kidney. miRNA expression analysis suggested a small number of miRNAs involved in the smoltification process. However, target analysis of these miRNAs showed a regulatory role in growth, stress response, and immunity. This study is the first to evidence the interplaying among mRNAs and miRNAs and the structural relationship at the genome level during Atlantic salmon smoltification.

Gill Transcriptome Analysis Revealed the Difference in Gene Expression Between Freshwater and Seawater Acclimated Guppy (Poecilia reticulata)

Guppy (Poecilia reticulata) can adapt to a wide range of salinity changes. To investigate the gene expression changes in the guppy exposed to seawater, we characterized its gill transcriptome using RNA sequencing. Experimental fish were exposed to salinity increase from 0 to 30‰ within 4 days, while control fish were cultured in freshwater (0‰ salinity). Seven days after salinity exposure, the gills were sampled and the mortality within 2 weeks was recorded. No significant difference in the cumulative mortality at the second week was found between the two groups. Transcriptomic analysis identified 3477 differentially expressed genes (DEGs), including 1067 upregulated and 2410 downregulated genes. These DEGs were enriched in several biological processes, including ion transport, ion homeostasis, ATP biosynthetic process, metabolic process, and immune system process. Oxidative phosphorylation was the most activated pathway. DEGs involved in the pathway "endoplasmic reticulum (ER)-mediated phagocytosis," "starch and sucrose metabolism," and "steroid biosynthesis" were mainly downregulated; chemokines and interleukins involved in "cytokine-cytokine receptor interaction" were differentially expressed. The present results suggested that oxidative phosphorylation had essential roles in osmoregulation in the gills of seawater acclimated guppy, during which the decline in the expression of genes encoding V-ATPases and calreticulin had a negative effect on the phagocytosis and immune response. Besides, several metabolic processes including "starch and sucrose metabolism" and "steroid biosynthesis" were affected. This study elucidates transcriptomic changes in osmotic regulation, metabolism, and immunity in seawater acclimated guppy.