Novel insights into the response of Atlantic salmon (Salmo salar) to Piscirickettsia salmonis: Interplay of coding genes and lncRNAs during bacterial infection

Transcriptome analyses reveal key roles of alternative splicing regulation in atlantic salmon during the infectious process of Piscirickettsiosis disease

In the Chilean salmon farming industry, infection by Piscirickettsia salmonis is the primary cause of the main bacterial disease known as Piscirickettsiosis, which has an overwhelming economic impact. Although it has been demonstrated that Piscirickettsiosis modifies the expression of numerous salmonids genes, it is yet unknown how alternative splicing (AS) contributes to salmonids bacterial infection. AS, has the potential to create heterogeneity at the protein and RNA levels and has been associated as a relevant molecular mechanism in the immune response of eukaryotes to several diseases. In this study, we used RNA data to survey P. salmonis-induced modifications in the AS of Atlantic salmon and found that P. salmonis infection promoted a substantial number (158,668) of AS events. Differentially spliced genes (DSG) sensitive to Piscirickettsiosis were predominantly enriched in genes involved in RNA processing, splicing and spliceosome processes (e.g., hnRNPm, hnRPc, SRSF7, SRSF45), whereas among the DSG of resistant and susceptible to Piscirickettsiosis, several metabolic and immune processes were found, most notably associated to the regulation of GTPase, lysosome and telomere organization-maintenance. Furthermore, we found that DSG were mostly not differentially expressed (5-7 %) and were implicated in distinct biological pathways. Therefore, our results underpin AS achieving a significant regulatory performance in the response of salmonids to Piscirickettsiosis.

Alternative splicing in Atlantic salmon head kidney and SHK-1 cell line during the Piscirickettsia salmonis infection: A comparative transcriptome survey

Piscirickettsia salmonis, an intracellular bacterium in salmon aquaculture, is a big challenge because it is responsible for 54.2% of Atlantic salmon mortalities. In recent years, the high relevance of Alternative Splicing (AS) as a molecular mechanism associated with infectious conditions and host-pathogen interaction processes, especially in host immune activation, has been observed. Several studies have highlighted the role of AS in the host's immune response during viral, bacterial, and endoparasite infection. In the present study, we evaluated AS transcriptome profiles during P. salmonis infection in the two most used study models, SHK-1 cell line and salmon head kidney tissue. First, the SHK-1 cell line was exposed to P. salmonis infection at 0-, 7-, and 14-days post-infection (dpi). Following, total RNA was extracted for Illumina sequencing. On the other hand, RNA-Seq datasets of Atlantic salmon head kidney infected with the same P. salmonis strayingwase used. For both study models, the highest number of differentially alternative splicing (DAS) events was observed at 7 dpi, 16,830 DAS events derived from 9213 DAS genes in SHK-1 cells, and 13,820 DAS events from 7684 DAS genes in salmon HK. Alternative first exon (AF) was the most abundant AS type in the three infection times analyzed, representing 31% in SHK-1 cells and 228.6 in salmon HK; meanwhile, mutually exclusive exon (MX) was the least abundant. Notably, functional annotation of DAS genes in SHK-1 cells infected with P. salmonis showed a high presence of genes related to nucleotide metabolism. In contrast, the salmon head kidney exhibited many GO terms associated with immune response. Our findings reported the role of AS during P. salmonis infection in Atlantic salmon. These studies would contribute to a better understanding of the molecular bases that support the pathogen-host interaction, evidencing the contribution of AS regulating the transcriptional host response.

Applying genetic technologies to combat infectious diseases in aquaculture

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

The impact of Piscirickettsia salmonis infection on genome-wide DNA methylation profile in Atlantic Salmon

Salmon rickettsial septicaemia (SRS), caused by the bacteria Piscirickettsia salmonis (P. salmonis), is responsible for significant mortality in farmed Atlantic salmon in Chile. Currently there are no effective treatments or preventive measures for this disease, although genetic selection or genome engineering to increase salmon resistance to SRS are promising strategies. The accuracy and efficiency of these strategies are usually influenced by the available biological background knowledge of the disease. The aim of this study was to investigate DNA methylation changes in response to P. salmonis infection in the head kidney and liver tissue of Atlantic salmon, and the interaction between gene expression and DNA methylation in the same tissues. The head kidney and liver methylomes of 66 juvenile salmon were profiled using reduced representation bisulphite sequencing (RRBS), and compared between P. salmonis infected animals (3 and 9 days post infection) and uninfected controls, and between SRS resistant and susceptible fish. Methylation was correlated with matching RNA-Seq data from the same animals, revealing that methylation in the first exon leads to an important repression of gene expression. Head kidney methylation showed a clear response to the infection, associated with immunological processes such as actin cytoskeleton regulation, phagocytosis, endocytosis and pathogen associated pattern receptor signaling. Our results contribute to the growing understanding of the role of methylation in regulation of gene expression and response to infectious diseases and could inform the incorporation of epigenetic markers into genomic selection for disease resistant and the design of diagnostic epigenetic markers to better manage fish health in salmon aquaculture.

Effect of low-dose Piscirickettsia salmonis infection on haematological-biochemical blood parameters in Atlantic salmon (Salmo salar)

Piscirickettsia salmonis is the etiological agent of Piscirickettsiosis, a severe disease that affects Atlantic salmon (Salmo salar) farmed in Chile and many other areas (Norway, Scotland, Ireland, Canada and the USA). This study investigated the effects of low-dose P. salmonis infection (1 × 10²  CFU/ml) on Atlantic salmon. In this study, we challenged fish with an isolated representative of the EM-90 genogroup via intraperitoneal injection for 42 days. Infected fish displayed decreased haematocrit and haemoglobin levels at day 13 post-infection, indicating erythropenia, haemolysis and haemodilution. Conversely, their white blood cell counts increased on days 13 and 21 post-infection. Additionally, their iron levels decreased from day 2 post-infection, indicating iron deficiency and an inability to retrieve stored iron before infection. Their magnesium levels also decreased at day 28 post-infection, possibly due to osmoregulatory problems. Also, we observed an increase in lactate dehydrogenase activity on days 5, 21, and 28 post-infection, suggesting early symptoms of hepatotoxicity. Later analyses determined a decrease in plasma glucose levels from day 2 post-infection. This may be attributed to the hypoxic conditions caused by P. salmonis, leading to an excess utilization of stored carbohydrates. Our results suggest that the blood parameters we studied are useful for monitoring the physiological status of Atlantic salmon infected with P. salmonis.

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.

Identification of mRNA-miRNA-lncRNA regulatory network associated with the immune response to Aeromonas salmonicides infection in the black rockfish (Sebastes schlegelii)

In aquaculture, Aeromonas salmonicides (A. salmonicida) is a main fish pathogen because of its nearly worldwide distribution, and broad host range. Recently, an increasing number of evidences have uncovered the roles of mRNA-miRNA-lncRNA network in fish diseases. In current study, RNA-seq was conducted in the black rockfish spleen following A. salmonicida infection at 0 h (Sp0 or control) and three different post-infection time-points (2 h: Sp2, 12 h: Sp12 and 24 h: Sp24, respectively) to comprehensively identify differentially expressed (DE) mRNAs, miRNAs and lncRNAs. Enrichment analysis and protein-protein interaction (PPI) analysis of DE mRNAs were performed. Then, expression and correlation analysis for mRNAs and their upstream miRNAs and lncRNAs were conducted. Finally, a total of 1364 mRNAs, 17 miRNAs and 1584 lncRNAs exhibited significantly differential expressions during bacterial infection in the black rockfish spleen. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Amino sugar and nucleotide sugar metabolism, Cell adhesion molecules (CAMs), Neuroactive ligand-receptor interaction, Nicotinate and nicotinamide metabolism, Pentose and glucuronate interconversions, Phagosome, Proteasome, etc. Subsequently, 1091 lncRNA-miRNA-mRNA pathways (323 in Sp2, 609 in Sp12 and 207 in Sp24) were constructed including 400 lncRNAs, 69 miRNAs, and 70 mRNAs. Meanwhile, NLRC3/novel-264/LNC_00116154 pathway demonstrated important immune modulating function in the black rockfish against A. salmonicida infection. Finally, the novel mRNA-miRNA-lncRNA sub-networks were established, among which all mRNAs and ncRNAs possessed significant predictive values for further studies for immune responses in the black rockfish.

PAMPs of Piscirickettsia salmonis Trigger the Transcription of Genes Involved in Nutritional Immunity in a Salmon Macrophage-Like Cell Line

The innate immune system can limit the growth of invading pathogens by depleting micronutrients at a cellular and tissue level. However, it is not known whether nutrient depletion mechanisms discriminate between living pathogens (which require nutrients) and pathogen-associated molecular patterns (PAMPs) (which do not). We stimulated SHK-1 cells with different PAMPs (outer membrane vesicles of Piscirickettsia salmonis "OMVs", protein extract of P. salmonis "TP" and lipopolysaccharides of P. salmonis "LPS") isolated from P. salmonis and evaluated transcriptional changes in nutritional immunity associated genes. Our experimental treatments were: Control (SHK-1 stimulated with bacterial culture medium), OMVs (SHK-1 stimulated with 1μg of outer membrane vesicles), TP (SHK-1 stimulated with 1μg of total protein extract) and LPS (SHK-1 stimulated with 1μg of lipopolysaccharides). Cells were sampled at 15-, 30-, 60- and 120-minutes post-stimulation. We detected increased transcription of zip8, zip14, irp1, irp2 and tfr1 in all three experimental conditions and increased transcription of dmt1 in cells stimulated with OMVs and TP, but not LPS. Additionally, we observed generally increased transcription of ireg-1, il-6, hamp, irp1, ft-h and ft-m in all three experimental conditions, but we also detected decreased transcription of these markers in cells stimulated with TP and LPS at specific time points. Our results demonstrate that SHK-1 cells stimulated with P. salmonis PAMPs increase transcription of markers involved in the transport, uptake, storage and regulation of micronutrients such as iron, manganese and zinc.

Why Does Piscirickettsia salmonis Break the Immunological Paradigm in Farmed Salmon? Biological Context to Understand the Relative Control of Piscirickettsiosis

Piscirickettsiosis (SRS) has been the most important infectious disease in Chilean salmon farming since the 1980s. It was one of the first to be described, and to date, it continues to be the main infectious cause of mortality. How can we better understand the epidemiological situation of SRS? The catch-all answer is that the Chilean salmon farming industry must fight year after year against a multifactorial disease, and apparently only the environment in Chile seems to favor the presence and persistence of Piscirickettsia salmonis. This is a fastidious, facultative intracellular bacterium that replicates in the host’s own immune cells and antigen-presenting cells and evades the adaptive cell-mediated immune response, which is why the existing vaccines are not effective in controlling it. Therefore, the Chilean salmon farming industry uses a lot of antibiotics—to control SRS—because otherwise, fish health and welfare would be significantly impaired, and a significantly higher volume of biomass would be lost per year. How can the ever-present risk of negative consequences of antibiotic use in salmon farming be balanced with the productive and economic viability of an animal production industry, as well as with the care of the aquatic environment and public health and with the sustainability of the industry? The answer that is easy, but no less true, is that we must know the enemy and how it interacts with its host. Much knowledge has been generated using this line of inquiry, however it remains insufficient. Considering the state-of-the-art summarized in this review, it can be stated that, from the point of view of fish immunology and vaccinology, we are quite far from reaching an effective and long-term solution for the control of SRS. For this reason, the aim of this critical review is to comprehensively discuss the current knowledge on the interaction between the bacteria and the host to promote the generation of more and better measures for the prevention and control of SRS.

Global gene expression responses of Atlantic salmon skin to Moritella viscosa

Moritella viscosa is a Gram-negative pathogen that causes large, chronic ulcers, known as winter-ulcer disease, in the skin of several fish species including Atlantic salmon. We used a bath challenge approach to profile the transcriptome responses of M. viscosa-infected Atlantic salmon skin at the lesion (Mv-At) and away from the lesion (Mv-Aw) sites. M. viscosa infection was confirmed through RNA-based qPCR assays. RNA-Seq identified 5212 and 2911 transcripts differentially expressed in the Mv-At compared to no-infection control and Mv-Aw groups, respectively. Also, there were 563 differentially expressed transcripts when comparing the Mv-Aw to control samples. Our results suggest that M. viscosa caused massive and strong, but largely infection site-focused, transcriptome dysregulations in Atlantic salmon skin, and its effects beyond the skin lesion site were comparably subtle. The M. viscosa-induced transcripts of Atlantic salmon were mainly involved in innate and adaptive immune response-related pathways, whereas the suppressed transcripts by this pathogen were largely connected to developmental and cellular processes. As validated by qPCR, M. viscosa dysregulated transcripts encoding receptors, signal transducers, transcription factors and immune effectors playing roles in TLR- and IFN-dependent pathways as well as immunoregulation, antigen presentation and T-cell development. This study broadened the current understanding of molecular pathways underlying M. viscosa-triggered responses of Atlantic salmon, and identified biomarkers that may assist to diagnose and combat this pathogen.

Comparative Analysis of Salmon Cell Lines and Zebrafish Primary Cell Cultures Infection with the Fish Pathogen Piscirickettsia salmonis

Piscirickettsia salmonis is the etiologic agent of piscirickettsiosis, a disease that causes significant losses in the salmon farming industry. In order to unveil the pathogenic mechanisms of P. salmonis, appropriate molecular and cellular studies in multiple cell lines with different origins need to be conducted. Toward that end, we established a cell viability assay that is suitable for high-throughput analysis using the alamarBlue reagent to follow the distinct stages of the bacterial infection cycle. Changes in host cell viability can be easily detected using either an absorbance- or fluorescence-based plate reader. Our method accurately tracked the infection cycle across two different Atlantic salmon-derived cell lines, with macrophage and epithelial cell properties, and zebrafish primary cell cultures. Analyses were also carried out to quantify intracellular bacterial replication in combination with fluorescence microscopy to visualize P. salmonis and cellular structures in fixed cells. In addition, dual gene expression analysis showed that the pro-inflammatory cytokines IL-6, IL-12, and TNFα were upregulated, while the cytokines IL1b and IFNγ were downregulated in the three cell culture types. The expression of the P. salmonis metal uptake and heme acquisition genes, together with the toxin and effector genes ospD3, ymt, pipB2 and pepO, were upregulated at the early and late stages of infection regardless of the cell culture type. On the other hand, Dot/Icm secretion system genes as well as stationary state and nutrient scarcity-related genes were upregulated only at the late stage of P. salmonis intracellular infection. We propose that these genes encoding putative P. salmonis virulence factors and immune-related proteins could be suitable biomarkers of P. salmonis infection. The infection protocol and cell viability assay described here provide a reliable method to compare the molecular and cellular changes induced by P. salmonis in other cell lines and has the potential to be used for high-throughput screenings of novel antimicrobials targeting this important fish intracellular pathogen.

Integrated analysis of mRNA and long non-coding RNA expression profiles reveals the potential roles of lncRNA-mRNA network in carp macrophage immune regulation

Long non-coding RNAs (lncRNAs) have emerged as a hot topic in research as mounting evidence has indicated their transcriptional or post-transcriptional regulatory potential in multiple biological processes. Previous studies have revealed the involvement of lncRNAs in the immunoregulation of mammalian macrophages by changing mRNA expression; however, studies on the lncRNAs in fish macrophages and their potential roles in the immune system remain unknown. Primary macrophages were isolated from the head kidney (HK) of red common carp (Cyprinus carpio) and high-throughput lncRNA-mRNA sequencing was performed using the Illumina HiSeq platform. The results revealed that the most highly expressed mRNAs in primary HK macrophages were mainly involved in immune-related signal pathways. Furthermore, the most enriched immune-related GO term and KEGG pathway of the mRNAs were "immune system development" and "chemokine signaling pathway," respectively. A total of 20,333 lncRNAs, composed of 10,512 known and 9821 novel lncRNAs, were identified, and functional enrichment analysis of the lncRNA-mRNA network indicated that the expressed lncRNAs in primary HK macrophages could be associated with the regulation of multiple immune-related signaling pathways. In addition, the expressions of several selected lncRNAs and their related mRNAs were determined in carp macrophages following a 6-h exposure to lipopolysaccharide (LPS) and Poly(I: C), the results of which confirmed the co-expression regulation of lncRNAs and target mRNAs in the immune response of carp macrophages. These results suggest the correlative of the lncRNA-mRNA network in fish macrophage immune response, which may further affect the cross-talk of various signaling pathways by interaction with other network genes. Here, we provided fundamental data about the transcriptome profiles of primary HK macrophages from red common carp by analysis of the lncRNA-mRNA network, and ultimately suggest the potential roles of lncRNA-mRNA networks in immune regulation in teleost fish.

Integrated analysis of lncRNA and mRNA in liver of Megalobrama amblycephala post Aeromonas hydrophila infection

Background

As non-coding RNA molecules of more than 200 bp in length, long non-coding RNAs (lncRNAs) play a variety of roles in biological processes, including regulating the immune responses to bacterial infections. In recent years, there have been many in-depth studies on mammalian lncRNAs, but the relevant studies in fish are very limited. Meanwhile, since lncRNAs are not conserved among species, it is difficult to apply the existing results directly to unstudied species.

Results

To obtain the information of lncRNAs in Megalobrama amblycephala, one of the most economically important freshwater fish in China, also to better understand the biological significance of lncRNAs in the immunity system, the fish liver at 0, 4, 12, 24, and 72 h post Aeromonas hydrophila infection (hpi) were obtained for lncRNA-sequencing (lncRNA-seq). A total of 14,849 lncRNAs were identified, and 2196 lncRNAs showed significant differences at different time points post A. hydrophila infection. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that the target genes of the differentially expressed lncRNAs were enriched in several pathways related to immune such as apoptosis, inflammation, and immune response. Time-specific modules were then identified, using weighted correlation network analysis (WGCNA), and 28 modules significantly correlated with different time point after infection were found. Furthermore, four immune-related genes and six lncRNAs in the time-specific modules were subsequently verified by RT-qPCR.

Conclusions

The above findings reveal the discovery of widespread differentially expressed lncRNAs in the M. amblycephala liver post A. hydrophila infection, suggesting that lncRNAs might participate in the regulation of host response to bacterial infection, enriching the information of lncRNAs in teleost and providing a resources basis for further studies on the immune function of lncRNAs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07969-5.

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.

Transcriptomic analysis reveals a Piscirickettsia salmonis-induced early inflammatory response in rainbow trout skeletal muscle

Skeletal muscle is the most abundant tissue in teleosts and is essential for movement and metabolism. Recently, it has been described that skeletal muscle can express and secrete immune-related molecules during pathogen infection. However, the role of this tissue during infection is poorly understood. To determine the immunocompetence of fish skeletal muscle, juvenile rainbow trout (Oncorhynchus mykiss) were challenged with Piscirickettsia salmonis strain LF-89. P. salmonis is the etiological agent of piscirickettsiosis, a severe disease that has caused major economic losses in the aquaculture industry. This gram-negative bacterium produces a chronic systemic infection that involves several organs and tissues in salmonids. Using high-throughput RNA-seq, we found that 60 transcripts were upregulated in skeletal muscle, mostly associated with inflammatory response and positive regulation of interleukin-8 production. Conversely, 141 transcripts were downregulated in association with muscle filament sliding and actin filament-based movement. To validate these results, we performed in vitro experiments using rainbow trout myotubes. In myotubes coincubated with P. salmonis strain LF-89 at an MOI of 50, we found increased expression of the proinflammatory cytokine il1b and the pattern recognition receptor tlr5s 8 and 12 h after infection. These results demonstrated that fish skeletal muscle is an immunologically active organ that can implement an early immunological response against P. salmonis.

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.

Chimeric Protein IPath® with Chelating Activity Improves Atlantic Salmon's Immunity against Infectious Diseases

Infection processes displayed by pathogens require the acquisition of essential inorganic nutrients and trace elements from the host to survive and proliferate. Without a doubt, iron is a crucial trace metal for all living organisms and also a pivotal component in the host–parasite interactions. In particular, the host reduces the iron available to face the infectious disease, increasing iron transport proteins’ expression and activating the heme synthesis and degradation pathways. Moreover, recent findings have suggested that iron metabolism modulation in fish promotes the immune response by reducing cellular iron toxicity. We hypothesized that recombinant proteins related to iron metabolism could modulate the fish’s immune system through iron metabolism and iron-responsive genes. Here a chimeric iron transport protein (IPath^®) was bioinformatically designed and then expressed in a recombinant bacterial system. The IPath^® protein showed a significant chelating activity under in vitro conditions and biological activity. Taking this evidence, a vaccine candidate based on IPath^® was evaluated in Atlantic salmon challenged with three different fish pathogens. Experimental trials were conducted using two fish groups: one immunized with IPath^® and another injected with adjutant as the control group. After 400 accumulated thermal units (ATUs), two different infection trials were performed. In the first one, fish were infected with the bacterium Aeromonas salmonicida, and in a second trial, fish were exposed to the ectoparasite Caligus rogercresseyi and subsequently infected with the intracellular bacterium Piscirickettsia salmonis. Fish immunized with IPath^® showed a significant delay in the mortality curve in response to A. salmonicida and P. salmonis infections. However, no significant differences between infected and control fish groups were observed at the end of the experiment. Notably, sea lice burden reduction was observed in vaccinated Atlantic salmon. Transcriptional analysis evidenced a high modulation of iron-homeostasis-related genes in fish vaccinated with IPath^® compared to the control group during the infection. Moreover, increasing expression of Atlantic salmon IgT was associated with IPath^® immunization. This study provides evidence that the IPath^® protein could be used as an antigen or booster in commercial fish vaccines, improving the immune response against relevant pathogens for salmon aquaculture.

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.

Identification and Characterization of Long Non-coding RNAs in the Intestine of Olive Flounder (Paralichthys olivaceus) During Edwardsiella tarda Infection

Long non-coding RNAs (lncRNAs) play widespread roles in fundamental biological processes, including immune responses. The olive flounder (Paralichthys olivaceus), an important economical flatfish widely cultured in Japan, Korea, and China, is threatened by infectious pathogens, including bacteria, viruses, and parasites. However, the role of lncRNAs in the immune responses of this species against pathogen infections is not well-understood. Therefore, in this study, we aimed to identify lncRNAs in the intestine of olive flounder and evaluate their differential expression profiles during Edwardsiella tarda infection, which is an important zoonotic and intestinal pathogen. A total of 4,445 putative lncRNAs were identified, including 3,975 novel lncRNAs and 470 annotated lncRNAs. These lncRNAs had shorter lengths and fewer exons compared with mRNAs. In total, 115 differentially expressed lncRNAs (DE-lncRNAs) were identified during E. tarda infection. To validate the expression pattern of lncRNAs, six DE-lncRNAs were randomly selected for quantitative real-time PCR. The co-located and co-expressed mRNAs of DE-lncRNAs were predicted, which were used to conduct the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The target genes of DE-lncRNAs enriched numerous immune-related processes and exhibited a strong correlation with immune-related signaling pathways. To better understand the extensive regulatory functions of lncRNAs, the lncRNA–miRNA–mRNA regulatory networks were constructed, and two potential competing endogenous RNA (ceRNA) networks, LNC_001979-novel_171-Potusc2 and LNC_001979-novel_171-Podad1, were preliminarily identified from the intestine of olive flounders for the first time. In conclusion, this study provides an invaluable annotation and expression profile of lncRNAs in the intestine of olive flounder infected with E. tarda; this forms a basis for further studies on the regulatory function of lncRNAs in the intestinal mucosal immune responses of olive flounder.

The Dynamic Immune Response of Yellow Catfish (Pelteobagrus fulvidraco) Infected With Edwardsiella ictaluri Presenting the Inflammation Process

Edwardsiella ictaluri is a highly destructive pathogen in cultured yellow catfish, thus it was very necessary to study the immune response of yellow catfish against bacterial infection. In this study, RNA-Seq technology was used to study the immune response in two distinct tissues of yellow catfish at eight different time points (h) after E. ictaluri infection. The number of differentially expressed genes (DEGs) in the spleen and liver was low at 3 h and 6 h post-infection, respectively. Afterwards, the most number of DEGs in the spleen was detected at 72 h, while the number of DEGs in the liver maintained a high level from 24 h to 120 h. The GO and KEGG enrichment analyses of DEGs at different time points uncovered that cytokines were continuously transcribed at 6 h to 120 h; whereas the liver is the main organ that secretes the components of the complement system, and metabolic regulation was activated from 12 h to 120 h. Moreover, an overview of the inflammation response of yellow catfish was exhibited including pattern-recognition receptors, inflammatory cytokines, chemokines, complements, and inflammation-related signal pathways. The similar expression tendency of nine genes by qRT-PCR validated the accuracy of transcriptome analyses. The different transcriptomic profiles obtained from the spleen and liver will help to better understand the dynamic immune response of fish against bacterial infection, and will provide basic information for establishing effective measures to prevent and control diseases in fish.

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.

Global Proteomic Profiling of Piscirickettsia salmonis and Salmon Macrophage-Like Cells during Intracellular Infection

Piscirickettsiasalmonis is an intracellular bacterial fish pathogen that causes piscirickettsiosis, a disease with numerous negative impacts in the Chilean salmon farming industry. Although transcriptomic studies of P. salmonis and its host have been performed, dual host-pathogen proteomic approaches during infection are still missing. Considering that gene expression does not always correspond with observed phenotype, and bacteriological culture studies inadequately reflect infection conditions, to improve the existing knowledge for the pathogenicity of P. salmonis, we present here a global proteomic profiling of Salmon salar macrophage-like cell cultures infected with P. salmonis LF-89. The proteomic analyses identified several P. salmonis proteins from two temporally different stages of macrophages infection, some of them related to key functions for bacterial survival in other intracellular pathogens. Metabolic differences were observed in early-stage infection bacteria, compared to late-stage infections. Virulence factors related to membrane, lipopolysaccharide (LPS) and surface component modifications, cell motility, toxins, and secretion systems also varied between the infection stages. Pilus proteins, beta-hemolysin, and the type VI secretion system (T6SS) were characteristic of the early-infection stage, while fimbria, upregulation of 10 toxins or effector proteins, and the Dot/Icm type IV secretion system (T4SS) were representative of the late-infection stage bacteria. Previously described virulence-related genes in P. salmonis plasmids were identified by proteomic assays during infection in SHK-1 cells, accompanied by an increase of mobile-related elements. By comparing the infected and un-infected proteome of SHK-1 cells, we observed changes in cellular and redox homeostasis; innate immune response; microtubules and actin cytoskeleton organization and dynamics; alteration in phagosome components, iron transport, and metabolism; and amino acids, nucleoside, and nucleotide metabolism, together with an overall energy and ATP production alteration. Our global proteomic profiling and the current knowledge of the P. salmonis infection process allowed us to propose a model of the macrophage-P. salmonis interaction.

Pharmacological iron-chelation as an assisted nutritional immunity strategy against Piscirickettsia salmonis infection

Salmonid Rickettsial Septicaemia (SRS), caused by Piscirickettsia salmonis, is a severe bacterial disease in the Chilean salmon farming industry. Vaccines and antibiotics are the current strategies to fight SRS; however, the high frequency of new epizootic events confirms the need to develop new strategies to combat this disease. An innovative opportunity is perturbing the host pathways used by the microorganisms to replicate inside host cells through host-directed antimicrobial drugs (HDAD). Iron is a critical nutrient for P. salmonis infection; hence, the use of iron-chelators becomes an excellent alternative to be used as HDAD. The aim of this work was to use the iron chelator Deferiprone (DFP) as HDAD to treat SRS. Here, we describe the protective effect of the iron chelator DFP over P. salmonis infections at non-antibiotic concentrations, in bacterial challenges both in vitro and in vivo. At the cellular level, our results indicate that DFP reduced the intracellular iron content by 33.1% and P. salmonis relative load during bacterial infections by 78%. These findings were recapitulated in fish, where DFP reduced the mortality of rainbow trout challenged with P. salmonis in 34.9% compared to the non-treated group. This is the first report of the protective capacity of an iron chelator against infection in fish, becoming a potential effective host-directed therapy for SRS and other animals against ferrophilic pathogens.

Potential Involvement of lncRNAs in the Modulation of the Transcriptome Response to Nodavirus Challenge in European Sea Bass (Dicentrarchus labrax L.)

Long noncoding RNAs (lncRNAs) are being increasingly recognised as key modulators of various biological mechanisms, including the immune response. Although investigations in teleosts are still lagging behind those conducted in mammals, current research indicates that lncRNAs play a pivotal role in the response of fish to a variety of pathogens. During the last several years, interest in lncRNAs has increased considerably, and a small but notable number of publications have reported the modulation of the lncRNA profile in some fish species after pathogen challenge. This study was the first to identify lncRNAs in the commercial species European sea bass. A total of 12,158 potential lncRNAs were detected in the head kidney and brain. We found that some lncRNAs were not common for both tissues, and these lncRNAs were located near coding genes that are primarily involved in tissue-specific processes, reflecting a degree of cellular specialisation in the synthesis of lncRNAs. Moreover, lncRNA modulation was analysed in both tissues at 24 and 72 h after infection with nodavirus. Enrichment analysis of the neighbouring coding genes of the modulated lncRNAs revealed many terms related to the immune response and viral infectivity but also related to the stress response. An integrated analysis of the lncRNAs and coding genes showed a strong correlation between the expression of the lncRNAs and their flanking coding genes. Our study represents the first systematic identification of lncRNAs in European sea bass and provides evidence regarding the involvement of these lncRNAs in the response to nodavirus.

The regulation mechanism of lncRNAs and mRNAs in sea cucumbers under global climate changes: Defense against thermal and hypoxic stresses

The aquatic environment can be greatly impacted by thermal and hypoxic stresses, particularly caused by intensified global warming. Hence, there is an urgency to understand the response mechanisms of marine organisms to adverse environment. Although long non-coding RNAs (lncRNAs) are involved in many biological processes, their roles in stress responses still remain unclear. Here, differentially expressed (DE) lncRNAs and mRNAs were identified as responses to environmental stresses in the economically important sea cucumber, Apostichopus japonicus, and their potential roles were explored. Based on a total of 159, 355 and 495 significantly upregulated genes and 230, 518 and 647 significantly downregulated genes identified in the thermal, hypoxic and combination thermal + hypoxic stress treatments, respectively, we constructed DE-lncRNA-mRNA coexpression networks. Among the networks, eight shared pairs were identified from the three treatments, and based on the connectivity degree, MSTRG.27265, MSTRG.19729 and MSTRG.95524 were shown to be crucial lncRNAs. Among all the significantly changed lncRNAs identified by RT-qPCR and sequencing data, binding sites were found in four other lncRNAs (MSTRG.34610, MSTRG.10941, MSTRG.81281 and MSTRG.93731) with Aja-miR-2013-3p, a key miRNA that responds to hypoxia in sea cucumbers. The hypoxia-inducible factor (HIF-1α) was also shown as the possible targeted mRNA of Aja-miR-2013-3p. As indicated by a dual-luciferase reporter assay system, "HIF-1α gene/Aja-miR-2013-3p/MSTRG.34610" network and the "HIF-1α gene/Aja-miR-2013-3p/MSTRG.10941" network may play important roles in sea cucumbers under environmental stresses. Moreover, environmental stress altered the expression of multiple lncRNAs and mRNAs, thus affecting various biological processes in A. japonicus, including immunity, energy metabolism and the cell cycle. At the molecular level, more comprehensive responses were elicited by the combined thermal/hypoxic stress treatment than by individual stresses alone in sea cucumbers. This study lays the groundwork for future research on molecular mechanisms of echinoderm responses to thermal and hypoxic stress in the context of global climate changes.

Identification and Evaluation of Long Noncoding RNAs in Response to Handling Stress in Red Cusk-Eel (Genypterus chilensis) via RNA-seq

The red cusk-eel (Genypterus chilensis) is a native species with strong potential to support Chilean aquaculture diversification. Under commercial conditions, fish are exposed to several stressors. To date, little is known about the mechanism involved in the stress response of red cusk-eel, and there is no information related to the regulation mediated by long noncoding RNAs (lncRNAs). The objective of this work was to identify for the first time the lncRNAs in the transcriptome of G. chilensis and to evaluate the differential expression levels of lncRNAs in the liver, head kidney, and skeletal muscle in response to handling stress. We used previously published transcriptome data to identify the lncRNAs by applying a series of filters based on annotation information in several databases to discard coding sequences. We identified a total of 14,614 putative lncRNAs in the transcriptome of red cusk-eel, providing a useful lncRNA reference resource to be used in future studies. We evaluated their differential expression in response to handling stress in the liver, head kidney, and skeletal muscle, identifying 112, 323, and 108 differentially expressed lncRNAs, respectively. The results suggest that handling stress in red cusk-eel generate an altered metabolic status in liver, altered immune response in head kidney, and skeletal muscle atrophy through an important coding and noncoding gene network. This is the first study that identifies lncRNAs in Genypterus genus and that evaluates the relation between handling stress and lncRNAs in teleost fish, thereby providing valuable information regarding noncoding responses to stress in Genypterus species.

Comparative modulation of lncRNAs in wild-type and rag1-heterozygous mutant zebrafish exposed to immune challenge with spring viraemia of carp virus (SVCV)

Although the modulation of immune-related genes after viral infection has been widely described in vertebrates, the potential implications of non-coding RNAs (ncRNAs), especially long non-coding RNAs (lncRNAs), in immunity are still a nascent research field. The model species zebrafish could serve as a useful organism for studying the functionality of lncRNAs due to the numerous advantages of this teleost, including the existence of numerous mutant lines. In this work, we conducted a whole-transcriptome analysis of wild-type (WT) and heterozygous rag1 mutant (rag1^+/-) zebrafish after infection with the pathogen spring viraemia of carp virus (SVCV). WT and rag1^+/- zebrafish were infected with SVCV for 24 h. Kidney samples were sampled from infected and uninfected fish for transcriptome sequencing. From a total of 198,540 contigs, 12,165 putative lncRNAs were identified in zebrafish. Most of the putative lncRNAs were shared by the two zebrafish lines. However, by comparing the lncRNA profiles induced after SVCV infection in WT and rag1^+/- fish, most of the lncRNAs that were significantly induced after viral challenge were exclusive to each line, reflecting a highly differential response to the virus. Analysis of the neighboring genes of lncRNAs that were exclusively modulated in WT revealed high representation of metabolism-related terms, whereas those from rag1^+/- fish showed enrichment in terms related to the adaptive immune response, among others. On the other hand, genes involved in numerous antiviral processes surrounded commonly modulated lncRNAs, as expected. These results clearly indicate that after SVCV infection in zebrafish, the expression of an array of lncRNAs with functions in different aspects of immunity is induced.

Catching the complexity of salmon-louse interactions

The study of host-parasite relationships is an integral part of the immunology of aquatic species, where the complexity of both organisms has to be overlayed with the lifecycle stages of the parasite and immunological status of the host. A deep understanding of how the parasite survives in its host and how they display molecular mechanisms to face the immune system can be applied for novel parasite control strategies. This review highlights current knowledge about salmon and sea louse, two key aquatic animals for aquaculture research worldwide. With the aim to catch the complexity of the salmon-louse interactions, molecular information gleaned through genomic studies are presented. The host recognition system and the chemosensory receptors found in sea lice reveal complex molecular components, that in turn, can be disrupted through specific molecules such as non-coding RNAs.

Pathogenicity of the bacterium New Zealand rickettsia-like organism (NZ-RLO2) in Chinook salmon Oncorhynchus tshawytscha smolt

Farmed New Zealand Chinook salmon Oncorhynchus tshawytscha Walbaum have been found to be infected by rickettsia-like organisms (NZ-RLO). While these Gram-negative intra-cellular bacteria are closely related to Piscirickettsia salmonis, a significant pathogen for farmed salmon globally, the pathogenicity of NZ-RLO is unknown. The aim of the present study was to determine if one strain, NZ-RLO2, causes disease in Chinook salmon. Post-smolt salmon were inoculated with NZ-RLO2 by intraperitoneal injection at high, medium and low doses and observed for 30 d. All fish in the high and medium dosed groups died by the end of the study and 63% of the low dose group died within 30 d of inoculation. Necropsy revealed the fish inoculated with NZ-RLO2 had internal multifocal haemorrhages. The most consistent histological finding in fish inoculated with NZ-RLO2 was neutrophilic and necrotizing pancreatitis and steatitis with intra-cytoplasmic organisms often visible within areas of inflammation. Other histological lesions included multifocal hepatic necrosis, haematopoietic cell necrosis and splenic and renal lymphoid depletion. The presence of NZ-RLO2 within the inoculated fish was confirmed by replication in cell culture and qPCR. The results suggest NZ-RLO2 can cause disease in Chinook salmon and therefore could be a significant pathogen in farmed Chinook salmon.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

Transcriptome analysis is a powerful tool that enables a deep understanding of complicated physiological pathways, including immune responses. RNA sequencing (RNA-Seq)-based transcriptome analysis and various bioinformatics tools have also been used to study non-model animals, including aquaculture species for which reference genomes are not available. Rapid developments in these techniques have not only accelerated investigations into the process of pathogenic infection and defense strategies in fish, but also used to identify immunity-related genes in fish. These findings will contribute to fish immunotherapy for the prevention and treatment of bacterial infections through the design of more specific and effective immune stimulants, adjuvants, and vaccines. Until now, there has been little information regarding the universality and diversity of immune reactions against pathogenic infection in fish. Therefore, one of the aims of this paper is to introduce the RNA-Seq technique for examination of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of immune responses against bacteria, based on our previous findings in largemouth bass (Micropterus salmoides) against Nocardia seriolae, gray mullet (Mugil cephalus) against Lactococcus garvieae, orange-spotted grouper (Epinephelus coioides) against Vibrio harveyi, and koi carp (Cyprinus carpio) against Aeromonas sobria, using RNA-seq techniques. We demonstrated that only 39 differentially expressed genes (DEGs) were present in all species. However, the number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs in largemouth bass against N. seriolae, 819 DEGs in mullets against L. garvieae, 909 in groupers against V. harveyi, and 1471 in carps against A. sobria. The DEGs in different fish species were also representative of specific immune-related pathways. The results of this study will enhance our understanding of the immune responses of fish, and will aid in the development of effective vaccines, therapies, and disease-resistant strains.

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.

Dual RNA-Seq Uncovers Metabolic Amino Acids Dependency of the Intracellular Bacterium Piscirickettsia salmonis Infecting Atlantic Salmon

High-throughput sequencing technologies have offered the possibility to understand the complexity of the transcriptomic responses of an organism during a wide variety of biological scenarios, such as the case of pathogenic infections. Recently, the simultaneous sequencing of both pathogen and host transcriptomes (dual RNA-seq) during the infection has become a promising approach to uncover the complexity of the host-pathogen interactions. In this study, through a double rRNA depletion and RNA sequencing protocols, we simultaneously analyzed the transcriptome of the intracellular bacterium Piscirickettsia salmonis and its host the Atlantic salmon (Salmo salar) during the course of the infection. Beyond canonical host immune-related response and pathogen virulent factors, both bacteria and host displayed a large number of genes associated with metabolism and particularly related with the amino acid metabolism. Notably, genome-wide comparison among P. salmonis genomes and different fish pathogens genomes revealed a lack of the biosynthetic pathway for several amino acids such as valine, leucine, and isoleucine. To support this finding, in vitro experiments evidenced that when these amino acids are restricted the bacterial growth dynamics is significantly affected. However, this condition is phenotypically reversed when the amino acids are supplemented in the bacterial growth medium. Based on our results, a metabolic dependency of P. salmonis on S. salar amino acids is suggested, which could imply novel mechanisms of pathogenesis based on the capacity to uptake nutrients from the host. Overall, dual transcriptome sequencing leads to the understanding of host-pathogen interactions from a different perspective, beyond biological processes related to immunity.

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5115-x) contains supplementary material, which is available to authorized users.

Comparative analysis of long non-coding RNAs in Atlantic and Coho salmon reveals divergent transcriptome responses associated with immunity and tissue repair during sea lice infestation

The increasing capacity of transcriptomic analysis by high throughput sequencing has highlighted the presence of a large proportion of transcripts that do not encode proteins. In particular, long non-coding RNAs (lncRNAs) are sequences with low coding potential and conservation among species. Moreover, cumulative evidence has revealed important roles in post-transcriptional gene modulation in several taxa. In fish, the role of lncRNAs has been scarcely studied and even less so during the immune response against sea lice. In the present study we mined for lncRNAs in Atlantic salmon (Salmo salar) and Coho salmon (Oncorhynkus kisutch), which are affected by the sea louse Caligus rogercresseyi, evaluating the degree of sequence conservation between these two fish species and their putative roles during the infection process. Herein, Atlantic and Coho salmon were infected with 35 lice/fish and evaluated after 7 and 14 days post-infestation (dpi). For RNA sequencing, samples from skin and head kidney were collected. A total of 5658/4140 and 3678/2123 lncRNAs were identified in uninfected/infected Atlantic and Coho salmon transcriptomes, respectively. Species-specific transcription patterns were observed in exclusive lncRNAs according to the tissue analyzed. Furthermore, neighbor gene GO enrichment analysis of the top 100 highly regulated lncRNAs in Atlantic salmon showed that lncRNAs were localized near genes related to the immune response. On the other hand, in Coho salmon the highly regulated lncRNAs were localized near genes involved in tissue repair processes. This study revealed high regulation of lncRNAs closely localized to immune and tissue repair-related genes in Atlantic and Coho salmon, respectively, suggesting putative roles for lncRNAs in salmon against sea lice infestation.

Differential lncRNA expression profiles reveal the potential roles of lncRNAs in antiviral immune response of Crassostrea gigas

Long noncoding RNAs (lncRNAs) may play widespread roles in various biological processes. However, systematic profiles of lncRNAs in the biological responses of Pacific Oyster (Crassostrea gigas) to pathogen infection have not yet been demonstrated. Here, we have conducted an exhaustive comparative transcriptome analysis using a bioinformatics approach to exam the functions of lncRNAs response to Ostreid herpesvirus 1μVar (OsHV-1μVar) challenge. In total, 101 differentially expressed lncRNAs (DE-lncRNA) during OsHV-1μVar infections were identified. Compared with differentially expressed mRNAs (DE-mRNA), DE-lncRNAs are shorter in terms of overall length but longer in terms of exon length. These lncRNAs shared similar characteristics with previously reported invertebrate lncRNAs, such as relatively low GC content, low exon number and low sequence conservation, but low expression level were not observed. 20 DE-lncRNAs are typically co-expressed with their neighboring genes annotated as GO terms (GO: 0044237), indicating that these lncRNAs are involved in binding and cellular process functions in cis mode. The weighted gene co-expression network (WGCNA) analysis resulted in 15 modules. The highlighted blue module was specifically demonstrated a co-expression relationship between 14 DE-lncRNAs and 17 immune-related DE-mRNAs (IR-DE-mRNA). Three hub lncRNAs within this module were co-expressed with one hub IR-DE-mRNA involved in fibrinogen-related protein. It was speculated that lncRNAs is extensively involved in oyster antiviral innate immune system. The present study will facilitate subsequently experimental studies to unravel the function of lncRNAs in marine invertebrate response to pathogen infection.

Under control: The innate immunity of fish from the inhibitors' perspective

The innate immune response involves a concerted network of induced gene products, preformed immune effectors, biochemical signalling cascades and specialised cells. However, the multifaceted activation of these defensive measures can derail or overshoot and, if left unchecked, overwhelm the host. A plenty of regulatory devices therefore mediate the fragile equilibrium between pathogen defence and pathophysiological manifestations. Over the past decade in particular, an almost complete set of teleostean sequences orthologous to mammalian immunoregulatory factors has been identified in various fish species, which prove the remarkable conservation of innate immune-control concepts among vertebrates. This review will present the current knowledge on more than 50 teleostean regulatory factors (plus additional fish-specific paralogs) that are of paramount importance for controlling the clotting cascade, the complement system, pattern-recognition pathways and cytokine-signalling networks. A special focus lies on those immunoregulatory features that have emerged as potential biomarker genes in transcriptome-wide research studies. Moreover, we report on the latest progress in elucidating control elements that act directly with immune-gene-encoding nucleic acids, such as transcription factors, hormone receptors and micro- and long noncoding RNAs. Investigations into the function of teleostean inhibitory factors are still mainly based on gene-expression profiling or overexpression studies. However, in support of structural and in-vitro analyses, evidence from in-vivo trials is also available and revealed many biochemical details on piscine immune regulation. The presence of multiple gene copies in fish adds a degree of complexity, as it is so far hardly understood if they might play distinct roles during inflammation. The present review addresses this and other open questions that should be tackled by fish immunologists in future.

Selection and validation of reliable housekeeping genes to evaluate Piscirickettsia salmonis gene expression

Piscirickettsia salmonis is a highly aggressive facultative intracellular bacterium that challenges the sustainability of Chilean salmon production. Due to the limited knowledge of its biology, there is a need to identify key molecular markers that could help define the pathogenic potential of this bacterium. We think a model system should be implemented that efficiently evaluates the expression of putative bacterial markers by using validated, stable, and highly specific housekeeping genes to properly select target genes, which could lead to identifying those responsible for infection and disease induction in naturally infected fish. Here, we selected a set of validated reference or housekeeping genes for RT-qPCR expression analyses of P. salmonis under different growth and stress conditions, including an in vitro infection kinetic. After a thorough screening, we selected sdhA as the most reliable housekeeping gene able to represent stable and highly specific host reference genes for RT-qPCR-driven P. salmonis analysis.

Epigenetic considerations in aquaculture

Epigenetics has attracted considerable attention with respect to its potential value in many areas of agricultural production, particularly under conditions where the environment can be manipulated or natural variation exists. Here we introduce key concepts and definitions of epigenetic mechanisms, including DNA methylation, histone modifications and non-coding RNA, review the current understanding of epigenetics in both fish and shellfish, and propose key areas of aquaculture where epigenetics could be applied. The first key area is environmental manipulation, where the intention is to induce an ‘epigenetic memory’ either within or between generations to produce a desired phenotype. The second key area is epigenetic selection, which, alone or combined with genetic selection, may increase the reliability of producing animals with desired phenotypes. Based on aspects of life history and husbandry practices in aquaculture species, the application of epigenetic knowledge could significantly affect the productivity and sustainability of aquaculture practices. Conversely, clarifying the role of epigenetic mechanisms in aquaculture species may upend traditional assumptions about selection practices. Ultimately, there are still many unanswered questions regarding how epigenetic mechanisms might be leveraged in aquaculture.

MicroRNA-based transcriptomic responses of Atlantic salmon during infection by the intracellular bacterium Piscirickettsia salmonis

MicroRNAs (miRNAs) are small non-coding RNAs that have emerged as key regulators in diverse biological processes across taxa. However, despite the importance of these transcripts, little is known about their role during the immune response in salmonids. Because of this, we use deep sequencing technologies to explore the microRNA-based transcriptomic response of the Atlantic salmon (Salmo salar) to the intracellular bacteria Piscirickettsia salmonis, one of the main threats to salmon aquaculture in Chile. Hence, 594 different miRNAs were identified from head kidney and spleen transcriptomic data. Among them, miRNA families mir-181, mir-143 and mir-21 were the most abundant in control groups, while after infection with P. salmonis, mir-21, mir-181 and mir-30 were the most predominant families. Furthermore, transcriptional analysis revealed 84 and 25 differentially expressed miRNAs in head kidney and spleen respectively, with an overlapping response of 10 miRNAs between the analyzed tissues. Target prediction, coupled with GO enrichment analysis, revealed that the possible targets of the most regulated miRNAs were genes involved in the immune response, such as cortisol metabolism, chemokine-mediated signaling pathway and neutrophil chemotaxis genes. Among these, predicted putative target genes such as C-C motif chemokine 19-like, stromal cell-derived factor 1-like, myxovirus resistance protein 2 and hepcidin-1 were identified. Overall, our results suggest that miRNA expression in co-modulation with transcription activity of target genes is related to putative roles of non-coding RNAs in the immune response of Atlantic salmon against intracellular bacterial pathogens.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

In the shadow: The emerging role of long non-coding RNAs in the immune response of Atlantic salmon

The genomic era has increased the research effort to uncover how the genome of an organism, and specifically the transcriptome, is modulated after interplaying with pathogenic microorganisms and ectoparasites. However, the ever-increasing accessibility of sequencing technology has also evidenced regulatory roles of long non-coding RNAs (lncRNAs) related to several biological processes including immune response. This study reports a high-confidence annotation and a comparative transcriptome analysis of lncRNAs from several tissues of Salmo salar infected with the most prevalent pathogens in the Chilean salmon aquaculture such as the infectious salmon anemia (ISA) virus, the intracellular bacterium Piscirickettsia salmonis and the ectoparasite copepod Caligus rogercresseyi. Our analyses showed that lncRNAs are widely modulated during infection. However, this modulation is pathogen-specific and highly correlated with immuno-related genes associated with innate immune response. These findings represent the first discovery for the widespread differential expression of lncRNAs in response to infections with different types of pathogens in Atlantic salmon, suggesting that lncRNAs are pivotal player during the fish immune response.

Combining 'omics and microscopy to visualize interactions between the Asian citrus psyllid vector and the Huanglongbing pathogen Candidatus Liberibacter asiaticus in the insect gut

Huanglongbing, or citrus greening disease, is an economically devastating bacterial disease of citrus. It is associated with infection by the gram-negative bacterium Candidatus Liberibacter asiaticus (CLas). CLas is transmitted by Diaphorina citri, the Asian citrus psyllid (ACP). For insect transmission to occur, CLas must be ingested during feeding on infected phloem sap and cross the gut barrier to gain entry into the insect vector. To investigate the effects of CLas exposure at the gut-pathogen interface, we performed RNAseq and mass spectrometry-based proteomics to analyze the transcriptome and proteome, respectively, of ACP gut tissue. CLas exposure resulted in changes in pathways involving the TCA cycle, iron metabolism, insecticide resistance and the insect's immune system. We identified 83 long non-coding RNAs that are responsive to CLas, two of which appear to be specific to the ACP. Proteomics analysis also enabled us to determine that Wolbachia, a symbiont of the ACP, undergoes proteome regulation when CLas is present. Fluorescent in situ hybridization (FISH) confirmed that Wolbachia and CLas inhabit the same ACP gut cells, but do not co-localize within those cells. Wolbachia cells are prevalent throughout the gut epithelial cell cytoplasm, and Wolbachia titer is more variable in the guts of CLas exposed insects. CLas is detected on the luminal membrane, in puncta within the gut epithelial cell cytoplasm, along actin filaments in the gut visceral muscles, and rarely, in association with gut cell nuclei. Our study provides a snapshot of how the psyllid gut copes with CLas exposure and provides information on pathways and proteins for targeted disruption of CLas-vector interactions at the gut interface.