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

Integration analysis of miRNA-mRNA uncovers the molecular immune mechanism of macrophage response to Aeromonas veronii infection

Macrophages are innate immunity cells which play pivotal roles in infectious immunity. Aeromonas veronii is a zoonotic agent capable of causing sepsis and poses a serious threat to public health. However, few studies have focused on miRNA-mRNA integration analysis to address the immune mechanisms of macrophage response to A. veronii infection. Herein, we characterized the immunophysiological, biochemical, and transcriptome changes of macrophage under A. veronii infection. We found that macrophages infected with A. veronii released large amounts of cytokines and triggered NLRP3-dependent pyroptosis. Subsequently, 603 differentially expressed miRNAs (DEMIs) and 3693 differentially expressed mRNAs (DEMs) were identified by RNA-seq analysis under A. veronii infection. Moreover, integrated analysis of miRNA-mRNA yielded 66 miRNA-target gene pairs composed of 41 DEMIs and 27 DEMs. We next identified the Toll-like receptor, NOD-like receptor, TNF and NF-κB pathways as necessary for macrophage to respond to A. veronii infection. miR-847 and miR-627 were involved in macrophage response to A. veronii infection by negatively regulating Pannexin-1 and thioredoxin interacting protein (TXNIP). Our findings elucidate the molecular mechanism of macrophage response to A. veronii infection at the miRNA level, providing many candidate miRNAs and mRNAs therapeutic targets for the prevention and treatment of A. veornii infectious diseases.

Transcriptomic analysis of turbot (Scophthalmus maximus) treated with zymosan a reveals that lncRNAs and inflammation-related genes mediate the protection conferred against Aeromonas salmonicida

Aeromonas salmonicida is one of the most harmful pathogens in finfish aquaculture worldwide. Immunostimulants such as β-glucans are used to enhance the immunity of cultured fish. However, their effects on fish physiology are not completely understood. In the present work, we evaluated the effect of a single intraperitoneal (ip) injection of zymosan A on fish survival against A. salmonicida infection. A single administration of this compound protected fish against A. salmonicida challenge and reduce the bacterial load in the head kidney one week after its administration. Transcriptome analyses of head kidney samples revealed several molecular mechanisms involved in the protection conferred by zymosan A and their regulation by long noncoding RNAs. The transcriptome profile of turbot exposed only to zymosan A was practically unaltered one week after ip injection. However, the administration of this immunostimulant induced significant transcriptomic changes once the fish were in contact with the bacteria and increased the survival of the infected turbot. Our results suggest that the restraint of the infection-induced inflammatory response, the management of apoptotic cell death, cell plasticity and cellular processes involving cytoskeleton dynamics support the protective effects of zymosan A. All this information provides insights on the cellular and molecular mechanisms involved in the protective effects of this widely used immunostimulant.

Genome-wide transcriptional profiling and functional analysis of long noncoding RNAs and mRNAs in chicken macrophages associated with the infection of avian pathogenic E. coli

BACKGROUND

Avian pathogenic E. coli (APEC) can cause localized or systemic infections, collectively known as avian colibacillosis, resulting in huge economic losses to poultry industry globally per year. In addition, increasing evidence indicates that long non-coding RNAs (lncRNAs) play a critical role in regulating host inflammation in response to bacterial infection. However, the role of lncRNAs in the host response to APEC infection remains unclear.

RESULTS

Here, we found 816 differentially expressed (DE) lncRNAs and 1,798 DE mRNAs in APEC infected chicken macrophages by RNAseq. The identified DE lncRNA-mRNAs were involved in Toll like receptor signaling pathway, VEGF signaling pathway, fatty acid metabolism, phosphatidylinositol signaling system, and other types of O-glycan biosynthesis. Furthermore, we found the novel lncRNA TCONS_00007391 as an important immune regulator in APEC infection was able to regulate the inflammatory response by directly targeting CD86.

CONCLUSION

These findings provided a better understanding of host response to APEC infection and also offered the potential drug targets for therapy development against APEC infection.

Effects of miR-722 on gene expression and alternative splicing in the liver of half-smooth tongue sole after infection with Vibrio anguillarum

MicroRNAs play crucial roles in various biological processes, including but not limited to differentiation, development, disease, and immunity. However, their immunoregulatory roles in half-smooth tongue sole are lacking. Our previous studies indicated that miR-722 could target C5aR1 to modulate the complement pathway to alleviate inflammatory response and even affect the mortality after the bacterial infection with Vibrio anguillarum. Driven by the purpose of revealing the underlying mechanisms, in this study, we investigated the effects of miR-722 on the gene expression and alternative splicing (AS) in the liver of half-smooth tongue sole after Vibrio anguillarum infection, with the approach of miR-722 overexpression/silencing and subsequent RNA-seq. Among the different comparisons, the I group (miR-722 inhibitor and V. anguillarum) versus blank control (PBS) exhibited the highest number of differentially expressed genes (DEGs), suggesting that the immune response was overactivated after inhibiting the miR-722. In addition, enrichment analyses were performed to reveal the functions of DEGs and differential AS (DAS) genes, reflecting the enrichment of RNA splicing and immune-related pathways including NF-κB and T cell receptor signaling pathway. Comparing the M group (miR-722 mimic and V. anguillarum) with the negative control (random sequence and V. anguillarum), two immune-related genes, cd48 and mapk8, were differentially expressed, of which mapk8 was also differentially spliced, indicating their importance in the immune response. Furthermore, representative gene analysis was performed, suggesting their corresponding functional changes due to AS. To verify the RNA-seq data, quantitative real-time PCR was employed with twenty pairs of primers for DEGs and DAS events. Overall, our results demonstrated that miR-722 could mediate the transcriptome-wide changes of gene expression and AS in half-smooth tongue sole, and provided insights into the regulatory role of miR-722 in immune responses, laying the foundation for further functional analyses and practical applications in aquaculture.

The evolution of NLRC3 subfamily genes in Sebastidae teleost fishes

BACKGROUND

With more than 36,000 valid fish species, teleost fishes constitute the most species-rich vertebrate clade and exhibit extensive genetic and phenotypic variation, including diverse immune defense strategies. NLRC3 subfamily genes, which are specific to fishes, play vital roles in the immune system of teleosts. The evolution of teleosts has been impacted by several whole-genome duplication (WGD) events, which might be a key reason for the expansions of the NLRC3 subfamily, but detailed knowledge of NLRC3 subfamily evolution in the family Sebastidae is still limited.

RESULTS

Phylogenetic inference of NLRC3 subfamily protein sequences were conducted to evaluate the orthology of NLRC3 subfamily genes in black rockfish (Sebastes schlegilii), 13 other fish species from the families Sebastidae, Serranidae, Gasterosteidae and Cyclopteridae, and three species of high vertebrates (bird, reptile and amphibian). WGD analyses were used to estimate expansions and contractions of the NLRC3 subfamily, and patterns of expression of NLRC3 subfamily genes in black rockfish following bacterial infections were used to investigate the functional roles of these genes in the traditional and mucosal immune system of the Sebastidae. Different patterns of gene expansions and contractions were observed in 17 fish and other species examined, and one and two whole-genome duplication events were observed in two members of family Sebastidae (black rockfish and honeycomb rockfish, Sebastes umbrosus), respectively. Subsequently, 179 copy numbers of NLRC3 genes were found in black rockfish and 166 in honeycomb rockfish. Phylogenetic analyses corroborated the conservation and evolution of NLRC3 orthologues between Sebastidae and other fish species. Finally, differential expression analyses provided evidence of the immune roles of NLRC3 genes in black rockfish during bacterial infections and gene ontology analysis also indicated other functional roles.

CONCLUSIONS

We hypothesize that NLRC3 genes have evolved a variety of different functions, in addition to their role in the immune response, as a result of whole genome duplication events during teleost diversification. Importantly, this study had underscored the importance of sampling across taxonomic groups, to better understand the evolutionary patterns of the innate immunity system on which complex immunological novelties arose. Moreover, the results in this study could extend current knowledge of the plasticity of the immune system.

Systematic analysis of circRNA-related ceRNA networks of black rockfish (Sebastes schlegelii) in response to Aeromonas salmonicides infection

Aeromonas salmonicides is a type of Gram-negative bacteria and has become the main fish pathogen in aquaculture because of its characteristics of worldwide distribution, broad host range and potentially devastating impacts. In the past years, studies have been focused to explore the regulatory roles of circRNA-miRNA-mRNA network in fish diseases. However, there are only few systematic studies linked to the anti-bacterial roles of circRNA-related ceRNA networks in the spleen immune system of black rockfish (Sebastes schlegelii). In this study, the whole-transcriptome sequencing (RNA-seq) was conducted in the black rockfish spleen with A. salmonicida challenging. The differentially expressed (DE) circRNAs were identified comprehensively for the following enrichment analysis. Interactions of miRNA-circRNA pairs and miRNA-mRNA pairs were predicted for the construction of circRNA-related ceRNA regulatory networks. Then, protein-protein interaction (PPI) analysis of mRNAs from these ceRNA networks were conducted. Finally, a total number of 39 circRNAs exhibited significantly differential expressions during A. salmonicida infection in the black rockfish spleen in 4338 identified circRNAs from 12 samples in 4 libraries. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Endocytosis, FoxO signaling pathway, Jak-STST signaling pathway, Herpes simplex infection, etc. Subsequently, 290 circRNA-miRNA-mRNA pathways (91 at 2 h, 142 at 12 h and 65 at 24 h) were constructed including 31 circRNAs, 50 miRNAs, and 156 mRNAs. In conclusion, the circRNA-related ceRNA networks were established, which will provide some novel insights in molecular mechanistic investigations of anti-bacterial immune response in teleost. Also, these findings will propose significant predictive values for the development of methods of treatment and prevention in black rockfish after bacterial infection in the future.

Molecular characterization, expression analysis and function identification of TNFα in black rockfish (Sebastes schlegelii)

TNFα, as a pro-inflammatory cytokine, plays an important role in inflammation and immune homeostasis maintaining. However, the knowledge about the immune functions of teleost TNFα against bacterial infections is still limited. In this study, the TNFα was characterized from black rockfish (Sebastes schlegelii). The bioinformatics analyses showed the evolutionary conservations in sequence and structure. The expression levels of Ss_TNFα mRNA were significantly up-regulated in the spleen and intestine after Aeromonas salmonicides and Edwardsiella tarda infections, and dramatically down-regulated in PBLs after LPS and poly I:C stimulations. Meanwhile, the extremely up-regulated expressions of other inflammatory cytokines (especially for IL-1β and IL17C) were observed in the intestine and spleen after bacterial infection and down-regulations were obtained in PBLs. The significant regulation with expression patterns of Ss_TNFα and other inflammatory cytokine mRNAs illustrated the variations of immunity in different tissues and cells of black rockfish. The regulated functions of Ss_TNFα in the up/downstream signaling pathways were preliminarily verified on the transcription and translation levels. Subsequently, in vitro knockdown of Ss_TNFα in the intestine cells of black rockfish confirmed the important immune roles of Ss_TNFα. Finally, the apoptotic analyses were conducted in PBLs and intestine cells of black rockfish. The rapid increases of the apoptotic rates were obtained in both PBLs and intestine cells after treatment with rSs_TNFα, but distinct apoptotic rates at the early and late stages of apoptosis were observed between these two types of cells. The results of apoptotic analyses suggested that Ss_TNFα could trigger apoptosis of different cells in different strategies in black rockfish. Overall, the findings in this study indicated the important roles of Ss_TNFα in the immune system of black rockfish during pathogenic infection, as well as the potential function on biomarker for monitoring the health status.

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

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

LncRNA-ANAPC2 and lncRNA-NEFM positively regulates the inflammatory response via the miR-451/npr2/ hdac8 axis in grass carp

In grass carp (Ctenopharyngodon idella), septicemia is a systemic inflammatory response to bacterial infection and could be leaded to lethality. MiR-451 involved in septicemia progression has been reported. However, the underlying mechanism of miR-451 in septicemia induced inflammatory response remains to be revealed. In the present study, miR-451 was highly expressed in Aeromonas hydrophila induced CIK cells, opposite to lncRNA-ANAPC2 and lncRNA-NEFM expression. Furthermore, we found that miR-451 interacted with lncRNA-ANAPC2 and lncRNA-NEFM, also targeted the 3' UTR of npr2 and hdac8. In CIK cells, the inhibition of npr2 and hdac8 were down-regulated by lncRNA-ANAPC2 and lncRNA-NEFM knockdown, while downstream proinflammatory factors were inhibited. In a word, this study indicates that lncRNA-ANAPC2 and lncRNA-NEFM regulation the LPS-induced progression of inflammatory response by modulating miR-451/npr2/hdac8 axis.