Evaluation of reference genes for quantitative real-time PCR analysis of messenger RNAs and microRNAs in rainbow trout Oncorhynchus mykiss under heat stress

PACAP binds conserved receptors and modulates cytokine gene expression and protein secretion in trout cell lines

Antimicrobial peptides (AMPs) are an alternative to antibiotics for treatment and prevention of infections with a lower risk of bacterial resistance. Pituitary adenylate cyclase activating polypeptide (PACAP) is an outstanding AMP with versatile effects including antimicrobial activity and modulation of immune responses. The objective of this research was to study PACAP immunomodulatory effect on rainbow trout cell lines infected with Aeromonas salmonicida. PACAP from Clarias gariepinus (PACAP1) and a modified PACAP (PACAP5) were tested. RT-qPCR results showed that il1b and il8 expression in RTgutGC was significantly downregulated while tgfb expression was upregulated after PACAP treatment. Importantly, the concentration of IL-1β and IFN-γ increased in the conditioned media of RTS11 cells incubated with PACAP1 and exposed to A. salmonicida. There was a poor correlation between gene expression and protein concentration, suggesting a stimulation of the translation of IL-1β protein from previously accumulated transcripts or the cleavage of accumulated IL-1β precursor. In-silico studies of PACAP-receptor interactions showed a turn of the peptide characteristic of PACAP-PAC1 interaction, correlated with the higher number of interactions observed with this specific receptor, which is also in agreement with the higher PACAP specificity described for PAC1 compared to VPAC1 and VPACA2. Finally, the in silico analysis revealed nine amino acids related to the PACAP receptor-associated functionality.

The effects of Pediococcus acidilactici MA18/5M on growth performance, gut integrity, and immune response using in vitro and in vivo Pacific salmonid models

Microbial management is central to aquaculture's efficiency. Pediococcus acidilactici MA18/5M has shown promising results promoting growth, modulation of the immune response, and disease resistance in many fishes. However, the mechanisms through which this strain confers health benefits in fish are poorly understood, particularly in Pacific salmonid models. Briefly, the aims of this study were to i) assess the protective effects of P. acidilactici MA18/5M by examining gut barrier function and the expression of tight junction (TJ) and immune genes in vitro and in vivo, and ii) to determine the protective effects of this strain against a common saltwater pathogen, Vibrio anguillarum J382. An in vitro model of the salmonid gut was employed utilizing the cell line RTgutGC. Barrier formation and integrity assessed by TEER measurements in RTgutGC, showed a significant decrease in resistance in cells exposed only to V. anguillarum J382 for 24 h, but pre-treatment with P. acidilactici MA18/5M for 48 h mitigated these effects. While P. acidilactici MA18/5M did not significantly upregulate tight junction and immune molecules, pre-treatment with this strain protected against pathogen-induced insults to the gut barrier. In particular, the expression of ocldn was significantly induced by V. anguillarum J382, suggesting that this molecule might play a role in the host response against this pathogen. To corroborate these observations in live fish, the effects of P. acidilactici MA18/5M was evaluated in Chinook salmon reared in real aquaculture conditions. Supplementation with P. acidilactici MA18/5M had no effect on Chinook salmon growth parameters after 10 weeks. Interestingly, histopathological results did not show alterations associated with P. acidilactici MA18/5M supplementation, indicating that this strain is safe to be used in the industry. Finally, the expression pattern of transcripts encoding TJ and immune genes in all the treatments suggest that variation in expression is more likely to be due to developmental processes rather than P. acidilactici MA18/5M supplementation. Overall, our results showed that P. acidilactici MA18/5M is a safe strain for use in fish production, however, to assess the effects on growth and immune response previously observed in other salmonid species, an assessment in adult fish is needed.

Involvement of miR-495 in the skin pigmentation of rainbow trout (Oncorhynchus mykiss) through the regulation of mc1r

MicroRNAs (miRNAs) play crucial roles in skin pigmentation in animals. Rainbow trout (Oncorhynchus mykiss) is a key economic fish species worldwide, and skin color directly affects its economic value. However, the functions of miRNAs in rainbow trout skin pigmentation remain largely unknown. Herein, we overexpressed and silenced miR-495 in vitro and in vivo to investigate its functions. The analysis of spatial and temporal expression patterns suggested that miR-495 is a potential regulator during the process of skin pigmentation. In vitro, mc1r was validated as a direct target for miR-495 by dual-luciferase reporter assay, and overexpression of miR-495 significantly inhibited mc1r expression; in contrast, mc1r and its downstream gene mitf levels were markedly upregulated by decreased miR-495. In vivo, overexpressed miR-495 by injecting agomiR-495 led to a substantial decrease in the expression of mc1r and mitf in dorsal skin and liver, while the opposite results were obtained after miR-495 silencing by antagomiR-495. These findings suggested that miR-495 can target mc1r to regulate rainbow trout skin pigmentation, which provide a potential basis for using miRNAs as target drugs to treat pigmentation disorders and melanoma.

Screening and Stability Analysis of Reference Genes for Gene Expression Normalization in Hybrid Yellow Catfish (Pelteobagrus fulvidraco ♀ × Pelteobagrus vachelli ♂) Fed Diets Containing Different Soybean Meal Levels

In this study, we screened the expression stability of six reference genes (18S rRNA, β-actin, GAPDH, EF1a, B2M, and HPRT1) in hybrid yellow catfish (n = 6), considering the SBM levels, sampling time points, and different tissues. Four different statistical programs, BestKeeper, NormFinder, Genorm, and Delta Ct, combined with a method that comprehensively considered all results, were used to evaluate the expression stability of these reference genes systematically. The results showed that SBM levels significantly impacted the expression stability of most of the reference genes studied and that this impact was time-, dose-, and tissue-dependent. The expression stability of these six reference genes varied depending on tissue, sampling time point, and SBM dosage. Additionally, more variations were found among different tissues than among different SBM levels or sampling time points. Due to its high expression, 18S rRNA was excluded from the list of candidate reference genes. β-actin and GAPDH in the liver and β-actin, HPRT1 and EF1a in the intestine were the most stable reference genes when SBM levels were considered. HPRT1, and EF1a in tissues sampled at 2 W and EF1a and β-actin in tissues sampled at 4 and 6 W were proposed as two stable reference genes when different tissues were considered. When the sampling time points were considered, β-actin, EF1a, and HPRT1 were the top three stable reference genes in the intestine. In contrast, β-actin and B2M are the most stable reference genes in the liver. In summary, β-actin, EF1a, and HPRT1 were the more stable reference genes in this study. The stability of reference genes depends on the tissues, sampling time points, and SBM diet levels in hybrid yellow catfish. Therefore, attention should be paid to these factors before selecting suitable reference genes for normalizing the target genes.

miR-330 targeting BCO2 is involved in carotenoid metabolism to regulate skin pigmentation in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

MicroRNAs (miRNAs) play a critical role in regulating skin pigmentation. As a key economic trait, skin color directly affects the market value of rainbow trout (Oncorhynchus mykiss), however, the regulatory mechanism of most miRNAs in fish skin color is still unclear.

RESULTS

In this study, the full-length cDNA sequence of β-carotene oxygenase 2 (BCO2, a key regulator of carotenoid metabolism) from the rainbow trout was obtained using rapid-amplification of cDNA ends (RACE) technology, and qRT-PCR was used to investigate the differential expression of miR-330 and BCO2 in 14 developmental stages and 13 tissues between wild-type rainbow trout (WTrt) and yellow mutant rainbow trout (YMrt). Additionally, the function of miR-330 was verified by overexpression and silencing in vitro and in vivo. The results showed that the complete cDNA sequence of BCO2 was 2057 bp with a 1707 bp ORF, encoding a 568 amino acid protein having a molecular weight of 64.07 kD. Sequence alignment revealed that higher conservation of BCO2 protein amongst fishes than amongst other vertebrates, which was further confirmed by phylogenetic analysis. The analysis of spatial and temporal expression patterns suggested that BCO2 and miR-330 were abundantly expressed from fertilized-stage to multi-cell as well as in the dorsal and ventral skin of WTrt and YMrt, and their expression patterns were opposite in most of the same periods and tissues. In vitro, luciferase reporter assay confirmed that BCO2 was a direct target of miR-330, and transfection of miR-330 mimics into rainbow trout liver cells resulted in a decrease in the expression of BCO2; conversely, miR-330 inhibitor had the opposite effect to the miR-330 mimics. In vivo, miR-330 agomir significantly decreased BCO2 expression in dorsal skin, tail fin, and liver. Furthermore, overexpression of miR-330 could suppress cell proliferation and induce apoptosis.

CONCLUSION

Our results showed that miR-330 is involved in the regulation of skin pigmentation in rainbow trout by targeting BCO2 and shows its promise as a potential molecular target to assist the selection of rainbow trout with better skin color patterns.

Comprehensive analysis of immune parameters, mRNA and miRNA profiles, and immune genes expression in the gill of rainbow trout infected with infectious hematopoietic necrosis virus (IHNV)

Rainbow trout (Oncorhynchus mykiss) is a species of cold-water fish with important economic values, widely cultivated worldwide. However, the outbreak of infectious hematopoietic necrosis virus (IHNV) caused the large-scale death of rainbow trout and seriously restricted the development of the trout farming industry. In this study, the changes of immune parameters in different periods (6-, 12-, 24-, 48-, 72-, 96-, 120-, and 144 h post-infection (hpi)), transcriptome profiles of 48 hpi (T48G) compared to control (C48G), and key immune-related genes expression patterns were measured in rainbow trout gill following IHNV challenge through biochemical methods, RNA sequencing (RNA-seq), and quantitative real-time polymerase chain reaction (qRT-PCR). The results showed that alkaline phosphatase (AKP), acid phosphatase (ACP), total superoxide dismutase (T-SOD), alanine aminotransferase (ALT), and aspartate aminotransferase (AST) activities, as well as lysozyme (LZM) and malonaldehyde (MDA) content decreased and then increased during infection, and remained at a high level after 48 hpi (P < 0.05), whereas catalase (CAT) activity showed a significant peak at 48 hpi (P < 0.05). The mRNA and miRNA analysis identified 4343 differentially expressed genes (DEGs) and 11 differentially expressed miRNAs (DEMs), and numerous immune-related DEGs involved in the Toll-like receptor signaling pathway, apoptosis, DNA replication, p53 signaling, RIG-I-like receptor signaling pathway, and NOD-like receptor signaling pathway and expression were significantly up-regulated in T48Gm group, including tlr3, tlr7, tlr8, traf3, ifih1, trim25, dhx58, ddh58, hsp90a.1, nlrc3, nlrc5, socs3, irf3, irf7, casp7, mx1, and vig2. The integrated analysis identified several important miRNAs (ola-miR-27d-3p_R+5, gmo-miR-124-3-5p, ssa-miR-301a-5p_L+2, and ssa-miR-146d-3p) that targeted key immune-related DEGs. Expression analysis showed that tlr3, tlr7, traf3, ifih1, dhx58, hap90a.1, irf3, irf7, and mx1 genes increased and then decreased during infection, and peaked at 72 hpi (P < 0.05). However, trim25 expression peaked at 96 hpi (P < 0.05). This study contributes to understanding immune response of rainbow trout against IHNV infection, and provides new insights into the immune regulation mechanisms and disease resistance breeding studies.

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

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

Integrated analysis of immune parameters, miRNA-mRNA interaction, and immune genes expression in the liver of rainbow trout following infectious hematopoietic necrosis virus infection

Rainbow trout (Oncorhynchus mykiss) is an important economical cold-water fish worldwide. However, infection with infectious hematopoietic necrosis virus (IHNV) has severely restricted the development of aquaculture and caused huge economic losses. Currently, little is known about the immune defense mechanisms of rainbow trout against IHNV. In this study, we detected the changes of immune parameters over different post-infection periods (6-, 12-, 24-, 48-, 72-, 96-, 120-, and 144 hours post-infection (hpi)), mRNA and miRNA expression profiles under 48 hpi (T48L) compared to control (C48L), and key immune-related genes expression patterns in rainbow trout liver following IHNV challenge through biochemical methods, RNA-seq, and qRT-PCR, and the function of miR-330-y was verified by overexpression and silencing in vitro and in vivo. The results revealed that alkaline phosphatase (AKP), alanine aminotransferase (ALT), catalase (CAT), and total superoxide dismutase (T-SOD) activities, and lysozyme (LZM) content showed significant peaks at 48 hpi, whereas malondialdehyde (MDA) content and aspartate aminotransferase (AST) activity decreased continuously during infection, and acid phosphatase (ACP) activity varied slightly. From RNA-seq, a total of 6844 genes and 86 miRNAs were differentially expressed, and numerous immune-related differentially expressed genes (DEGs) involved in RIG-I-like receptor signaling pathway, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, cytokine-cytokine receptor interaction, and antigen processing and presentation were significantly upregulated in T48Lm group, including IFIH1, DHX58, MAVS, TRAF3, IRF3, IRF7, MX1, TLR3, TLR8, MYD88, NOD1, NOD2, IL-8, CXCR1, CD209, CD83, and TAP1. Integrated analysis identified seven miRNAs (miR-425-x, miR-185-x, miR-338-x, miR-330-y, miR-361-x, miR-505-y, and miR-191-x) that target at least three key immune-related DEGs. Expression analysis showed that IFIH1, DHX58, IRF3, IRF7, MX1, TLR3, TLR8, and MYD88 showed a marked increase after 24 hpi during infection. Further research confirmed TAP1 as one of the targets of miR-330-y, overexpression of miR-330-y with mimics or agomir significantly reduced the expression levels of TAP1, IRF3, and IFN, and the opposite effects were obtained by inhibitor. These results facilitate in-depth understanding of the immune mechanisms in rainbow trout against IHNV.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss)

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

Evaluation of qPCR reference genes for taimen (Hucho taimen) under heat stress

As a powerful and attractive method for detecting gene expression, qRT-PCR has been broadly used in aquaculture research. Understanding the biology of taimen (Hucho taimen) has drawn increasing interest because of its ecological and economic value. Stable reference genes are required for the reliable quantification of gene expression, but such genes have not yet been optimized for taimen. In this study, the stability levels of 10 commonly used candidate reference genes were evaluated using geNorm, NormFinder, BestKeeper, and RefFinder. The expression levels of the 10 genes were detected using 240 samples from 48 experimental groups consisting of 40 individuals treated under four heat-stress conditions (18, 20, 22, and 24 °C) for 24 h and 26 °C for 4, 24, 48, and 72 h. Six tissues (blood, heart, brain, gill, skin, and liver) were collected from each individual. Ribosomal protein S29 (RPS29) and ribosomal protein L19 (RPL19) were the most stable genes among all of the samples, whereas 28S ribosomal RNA (28S rRNA), attachment region binding protein (ARBP), and 18S ribosomal RNA (18S rRNA) were the least stable. These results were verified by an expression analysis of taimen heat-stress genes (heat shock protein 60, hsp60, and heat shock protein 70, hsp70). In conclusion, RPS29 and RPL19 are the optimal reference genes for qRT-PCR analyses of taimen, irrespective of the tissue and experimental conditions. These results allow the reliable study of gene expression in taimen.

Genome-wide identification of Hsp70/110 genes in rainbow trout and their regulated expression in response to heat stress

Heat shock proteins (Hsps) play an important role in many biological processes. However, as a typical cold water fish, the systematic identification of Hsp70/110 gene family of rainbow trout (Oncorhynchus mykiss) has not been reported, and the role of Hsp70/110 gene in the evolution of rainbow trout has not been described systematically. In this study, bioinformatics methods were used to analyze the Hsp70/110 gene family of rainbow trout. A total of 16 hsp70/110 genes were identified and classified into ten subgroups. The 16 Hsp70/110 genes were all distributed on chromosomes 2, 4, 8 and 13. The molecular weight is ranged from 78.93 to 91.39 kD. Gene structure and motif composition are relatively conserved in each subgroup. According to RNA-seq analysis of rainbow trout liver and head kidney, a total of four out of 16 genes were significantly upregulated in liver under heat stress, and a total of seven out of 16 genes were significantly upregulated in head kidney. RT-qPCR was carried out on these gene, and the result were consistent with those of RNA-seq. The significantly regulated expressions of Hsp70/110 genes under heat stress indicats that Hsp70/110 genes are involved in heat stress response in rainbow trout. This systematic analysis provided valuable information about the diverse roles of Hsp70/110 in the evolution of teleost, which will contribute to the functional characterization of Hsp70/110 genes in further research.

Genome-Wide Identification of hsp90 Gene in Rainbow Trout (Oncorhynchus mykiss) and Their Regulated Expression in Response to Heat Stress

In this study, we analyzed the gene structure, chemical characterizations, chromosome locations, evolutionary relationship, and expression profile of hsp90 genes with online database. In addition, the expression levels of hsp90s were also investigated under heat stress by quantitative real-time (qRT)-PCR. A total of eight hsp90 genes were identified from the rainbow trout genome. They were all distributed on chromosomes 2, 4, 8, and 13. The molecular weight ranged from 78.93 to 91.39 kDa, and the isoelectric point ranged from 4.84 to 4.96. The eight hsp90 genes were clustered into six subfamilies (A, B, C, D, E, and F). Genetic structure and conserved domain analysis revealed that all eight hsp90 genes had only one exon, and motif 1-motif 10 was shared by most genes. According to RNA-seq analysis of rainbow trout liver and head kidney, a total of seven out of eight genes were significantly upregulated under heat stress, and qRT-PCR was carried out on these seven genes; the expression levels of these genes were significantly upregulated under heat stress. The significantly regulated expressions of hsp90 genes under heat stress indicated that hsp90 genes are involved in heat stress response in rainbow trout. This study provides a theoretical basis for further study on the role of hsp90 in the heat stress tolerance of rainbow trout.