miR210 modulates respiratory burst in Apostichopus japonicus coelomocytes via targeting Toll-like receptor

Regulatory mechanism of miR-722 on C5aR1 and its functions against bacterial inflammation in half-smooth tongue sole (Cynoglossus semilaevis)

MicroRNAs (miRNAs) are small non-coding RNAs involved in various biological processes, including immunity. Previously, we investigated the miRNAs of half-smooth tongue sole (Cynoglossus semilaevis) and found that miR-722 (designated Cse-miR-722) was significantly differentially expressed after infection with Vibrio anguillarum, reflecting its importance in immune response. Our preliminary bioinformatic analysis suggested that Cse-miR-722 could target C5aR1 (designated CsC5aR1), which was known to play crucial roles in complement activation and inflammatory response, as a receptor of C5a. However, the underlying mechanisms of their interactions and specific functions in inflammatory and immune response are still enigmas. In this study, we successfully cloned the precursor sequence of Cse-miR-722 (94 bp) and the full length of CsC5aR1 (1541 bp, protein molecular weight 39 kDa). The target gene of Cse-miR-722 was verified as CsC5aR1 by a dual luciferase reporter assay, and Cse-miR-722 was confirmed to regulate CsC5aR1 at the protein level using quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting. The expression of CsC5aR1 and Cse-miR-722 in liver cells and four immune tissues of half-smooth tongue sole changed significantly after LPS stimulation and infection with V. anguillarum. To explore the functional role of Cse-miR-722 in half-smooth tongue sole, we performed both in vitro and in vivo experiments. Cse-miR-722 was observed to affect phagocytosis and respiratory burst activity of macrophages by regulating CsC5aR1 in half-smooth tongue sole. Furthermore, we found that Cse-miR-722 regulated the expression of CsC5aR1, CsC5a, and the inflammatory factors CsIL1-β, CsIL6, CsIL8, and CsTNF-α both in vitro and in vivo. In addition, Cse-miR-722 reduced mortality and pathological damage. This study clarified the regulatory mechanism of Cse-miR-722 on CsC5aR1 and provided insight into the regulatory roles of Cse-miR-722 in immune responses, laying a theoretical foundation for the feasibility of using miR-722 to prevent and control bacterial diseases in teleost.

Exploring the Role of a Novel Interleukin-17 Homolog from Invertebrate Marine Mussel Mytilus coruscus in Innate Immune Response: Is Negative Regulation by Mc-Novel_miR_145 the Key?

Interleukin-17 (IL-17) represents a class of proinflammatory cytokines involved in chronic inflammatory and degenerative disorders. Prior to this study, it was predicted that an IL-17 homolog could be targeted by Mc-novel_miR_145 to participate in the immune response of Mytilus coruscus. This study employed a variety of molecular and cell biology research methods to explore the association between Mc-novel_miR_145 and IL-17 homolog and their immunomodulatory effects. The bioinformatics prediction confirmed the affiliation of the IL-17 homolog with the mussel IL-17 family, followed by quantitative real-time PCR assays (qPCR) to demonstrate that McIL-17-3 was highly expressed in immune-associated tissues and responded to bacterial challenges. Results from luciferase reporter assays confirmed the potential of McIL-17-3 to activate downstream NF-κb and its targeting by Mc-novel_miR_145 in HEK293 cells. The study also produced McIL-17-3 antiserum and found that Mc-novel_miR_145 negatively regulates McIL-17-3 via western blotting and qPCR assays. Furthermore, flow cytometry analysis indicated that Mc-novel_miR_145 negatively regulated McIL-17-3 to alleviate LPS-induced apoptosis. Collectively, the current results showed that McIL-17-3 played an important role in molluscan immune defense against bacterial attack. Furthermore, McIL-17-3 was negatively regulated by Mc-novel_miR_145 to participate in LPS-induced apoptosis. Our findings provide new insights into noncoding RNA regulation in invertebrate models.

circRNA432 enhances the coelomocyte phagocytosis via regulating the miR-2008-ELMO1 axis in Vibrio splendidus-challenged Apostichopus japonicus

Circular RNAs (circRNAs) are a kind of extensive and diverse covalently closed circular endogenous RNA, which exert crucial functions in immune regulation in mammals. However, the functions and mechanisms of circRNAs in invertebrates are largely unclarified. In our previous work, 261 differentially expressed circRNAs including circRNA432 (circ432) were identified from skin ulcer syndrome (SUS) diseased sea cucumber Apostichopus japonicus by RNA-seq. To better address the functional role of sea cucumber circRNAs, circ432 was first found to be significantly induced by Vibrio splendidus challenge and LPS exposure in this study. Knock-down circ432 could depress the V. splendidus-induced coelomocytes phagocytosis. Moreover, circ432 is validated to serve as the sponge of miR-2008, a differential expressed miRNA in SUS-diseased sea cucumbers, by Argonaute 2-RNA immunoprecipitation (AGO2-RIP) assay, luciferase reporter assay and RNA fluorescence in situ hybridization (FISH) in vitro. Engulfment and cell motility protein 1 (AjELMO1) is further demonstrated to be the target of miR-2008, and silencing AjELMO1 inhibits the V. splendidus-induced coelomocytes phagocytosis, and this phenomenon could be further suppressed by supplementing with miR-2008 mimics, suggesting that circ432 might regulate coelomocytes phagocytosis via miR-2008-AjELMO1 axis. We further confirm that the depressed coelomocytes' phagocytosis by circ432 silencing is consistent with the decreased abundance of AjELMO1, and could be recovered by miR-2008 inhibitors transfection. All our results provide the evidence that circ432 is involved in regulating pathogen-induced coelomocyte phagocytosis via sponge miR-2008 and promotes the abundance of AjELMO1. These findings will enrich the regulatory mechanism of phagocytosis in echinoderm and provide theoretical data for SUS disease prevention and control in sea cucumbers.

Phenotypic and functional characterization of two coelomocyte subsets in Apostichopus japonicus

The hemocytes of invertebrates are composed of different cell subsets with different morphologies and structures. Different cell subsets have different immune functions, which play an important role in innate immune response against pathogens. However, the understanding of the classification of Apostichopus japonicus coelomocytes and the molecular basis of immune function of different cell subsets is very limited. In this study, two coelomocyte subpopulations of A. japonicus were isolated by Percoll density gradient centrifugation. They were identified from their morphological and structural characteristics, namely, spherical cells with a size of 10-12 μm spherical in shape and a large number of small granules inside; lymphocyte-like cells with a size of 4-5 μm spherical or oval in shape, and 1-3 filopodia. Functionally, the phagocytic capacity and lysosomal activity in spherical cells were significantly greater than those in lymphocyte-like cells. The results suggest that spherical cells may play a more critical role in the immune responses. Meanwhile, transcriptome sequencing analysis was performed to further clarify the functional differences between the two cell subsets. The data indicated significantly different gene expression patterns in them. Spherical cells tend to participate in immune defense, whereas lymphocyte-like cells tend to participate in energy metabolism. In addition, lymphocyte-like cells may convert oxidative phosphorylation to glycolysis by changing the manner of energy metabolism to quickly adapt to the energy demand of external stimuli. Spherical cells may respond to LPS stimulation through phagocytosis, and their response time is slower than that of lymphocyte-like cells. The expression of genes involved in endocytosis, phagocytosis, and lysosomal and humoral immunity in spherical cells was significantly higher than that in lymphocyte-like cells. These data provide valuable information for understanding the molecular basis of cellular and humoral immunity in A. japonicus.

Differential expression of miRNAs in the body wall of the sea cucumber Apostichopus japonicus under heat stress

MicroRNAs, as one of the post-transcriptional regulation of genes, play an important role in the development process, cell differentiation and immune defense. The sea cucumber Apostichopus japonicus is an important cold-water species, known for its excellent nutritional and economic value, which usually encounters heat stress that affects its growth and leads to significant economic losses. However, there are few studies about the effect of miRNAs on heat stress in sea cucumbers. In this study, high-throughput sequencing was used to analyze miRNA expression in the body wall of sea cucumber between the control group (CS) and the heat stress group (HS). A total of 403 known miRNAs and 75 novel miRNAs were identified, of which 13 miRNAs were identified as significantly differentially expressed miRNAs (DEMs) in response to heat stress. A total of 16,563 target genes of DEMs were predicted, and 101 inversely correlated target genes that were potentially regulated by miRNAs in response to heat stress of sea cucumbers were obtained. Based on these results, miRNA-mRNA regulatory networks were constructed. The expression results of high-throughput sequencing were validated in nine DEMs and four differentially expressed genes (DEGs) by quantitative real-time polymerase chain reaction (RT-qPCR). Moreover, pathway enrichment of target genes suggested that several important regulatory pathways may play an important role in the heat stress process of sea cucumber, including ubiquitin-mediated proteolysis, notch single pathway and endocytosis. These results will provide basic data for future studies in miRNA regulation and molecular adaptive mechanisms of sea cucumbers under heat stress.

MiR-7 Regulates Pathogen-Induced Immune Response via PAK1 in the Sea Cucumber Apostichopus japonicus

MicroRNA-7 (miR-7) is a highly conserved short non-coding RNA involved in various bioprocesses via the regulation of multiple target genes. To enrich our knowledge of the functions of miR-7 in innate immune regulation in echinoderms, we first investigated the targeting relationship between miR-7 and PAK1 in the sea cucumber Apostichopus japonicus and then explored the functions of miR-7, the PAK1 gene, and the miR-7/PAK1 axis in the pathogen-induced immune response of A. japonicus. Our results showed that miR-7 can bind to the 3ʹUTR of PAK1 and negatively regulate the expression of PAK1 in A. japonicus. Overexpression and inhibition of miR-7 and inhibition of the expression of PAK1 can alter phagocytosis, cellular agglutination, and lysozyme contents in A. japonicus. Both miR-7 and the PAK1 gene are involved in immune defense against Vibrio splendidus infection; the miR-7/AjPAK1 axis showed immune regulatory function at 48 to 72 h post-infection (hpi) after V. splendidus infection in A. japonicus. In summary, the results of this study established that miR-7 regulates the pathogen-induced immune response by targeting PAK1 in A. japonicus.

Expression characteristics and interaction networks of microRNAs in spleen tissues of grass carp (Ctenopharyngodon idella)

The spleen is an important immune organ in fish. MicroRNAs (miRNAs) have been shown to play an important role in the regulation of immune function. However, miRNA expression profiles and their interaction networks associated with the postnatal late development of spleen tissue are still poorly understood in fish. The grass carp (Ctenopharyngodon idella) is an important economic aquaculture species in China. Here, two small RNA libraries were constructed from the spleen tissue of healthy grass carp at one-year-old and three-year-old. A total of 324 known conserved miRNAs and 9 novel miRNAs were identified by using bioinformatic analysis. Family analysis showed that 23 families such as let-7, mir-1, mir-10, mir-124, mir-8, mir-7, mir-9, and mir-153 were highly conserved between vertebrates and invertebrates. In addition, 14 families such as mir-459, mir-430, mir-462, mir-7147, mir-2187, and mir-722 were present only in fish. Expression analysis showed that the expression patterns of miRNAs in the spleen of one-year-old and three-year-old grass carp were highly consistent, and the percentage of miRNAs with TPM > 100 was above 39%. Twenty significant differentially expressed (SDE) miRNAs were identified. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that these SDE miRNAs were primarily involved in erythrocyte differentiation, lymphoid organ development, immune response, lipid metabolic process, the B cell receptor signaling pathway, the T cell receptor signaling pathway, and the PPAR signaling pathway. In addition, the following miRNA-mRNA interaction networks were constructed: immune and hematopoietic, cell proliferation and differentiation, and lipid metabolism. This study determined the miRNA transcriptome as well as miRNA-mRNA interaction networks in normal spleen tissue during the late development stages of grass carp. The results expand the number of known miRNAs in grass carp and are a valuable resource for better understanding the molecular biology of the spleen development in grass carp.

Identification of a novel RhoA gene in the sea cucumber Apostichopus japonicus and its immune regulatory function via interacting with miR-2012-5p

RhoA (Ras homolog A) protein is a representative member of the Rho GTPase family and is involved in various cellular processes. The function of RhoA in sea cucumbers is unclear. In this study, we hypothesized that RhoA may regulate the innate immune response of Apostichopus japonicus. Our data showed that 1) the complete sequence of RhoA from A. japonicus (named AjRhoA) was 968 bp, with a high level sequence conservation across the echinoderms and other phyla; 2) tissue expression analysis showed that AjRhoA transcripts and protein exhibited higher abundance in coelomocytes, whereas the relative expression of miR-2012-5p was lower in coelomocytes; 3) interactive binding sites and a negative regulatory targeting relationship between AjRhoA and miR-2012-5p were confirmed through a dual-luciferase reporter assay and functional validation in vivo; 4) the relative expression levels of AjRhoA transcripts and protein were upregulated in coelomocytes 4- and 72-hour post infection (hpi) with Vibrio splendidus, whereas miR-2012-5p was expressed in the opposite pattern; 5) both AjRhoA silencing and miR-2012-5p overexpression suppressed the phagocytic capacity of A. japonicus compared with the control at 4 and 72 hpi. Our observations suggest that AjRhoA can regulate the pathogen-induced immune response of A. japonicus through the "AjRhoA-miR-2012-5p" module during the early infection, while miR-2012-5p plays a direct immunomodulatory role as the infection progresses.

Downregulation of miR-210 Promoted Apoptosis of Hippocampal Neurons by Negatively Regulating the TLR4/NF-кB1 Signaling Pathway in a Rat Model of Status Epilepticus

PURPOSE

Status epilepticus (SE) is a life-threatening condition causing brain damage, hippocampal necrosis and apoptosis. This study aimed to determine whether microRNA-210 regulates seizure and apoptosis by targeting the TLR4 /NF-κB1 associated signaling pathway.

METHODS

In a pilocarpine-induced epileptic rat model, the expressions of microRNA-210 (miR-210), TLR4, NF-κB1 and caspase-3 were assessed by a quantitative polymerase chain reaction and Western blotting. Tunel detects hippocampal neuron apoptosis.

RESULTS

We found that miR-210, TLR4, NF-κB1 and caspase-3 were upregulated in the hippocampus of the rat model compared with that of control. The knockdown of miR-210 significantly restored the expression levels of TLR4, NF-κB1 and caspase-3 and increased hippocampal apoptosis.

CONCLUSION

These findings showed that the downregulation of miR-210 promoted apoptosis of hippocampal neurons by negatively regulating the TLR4/NF-кB1 signaling pathway.

Integrated analysis of miRNAome and transcriptome reveals miRNA-mRNA network regulation in Vibrio alginolyticus infected thick shell mussel Mytilus coruscus

The thick shell mussel Mytilus coruscus has developed into a model species for studying the interaction between molluscs and environmental stimuli. Herein, integrated analysis of miRNAome and transcriptome was performed to reveal miRNA-mRNA network regulation in Vibrio alginolyticus infected M. coruscus. There have detected some histological abnormalities in digestive gland and gills of V. alginolyticus challenged mussels, ascertaining the effective irritation by the present bacterial strain. A total of 265 novel miRNAs were finally predicted, of which 26 were differentially expressed miRNAs (DEMs). Additionally, 667 differentially expressed genes (DEGs) were detected, which may be potentially associated with innate immune response to V. alginolyticus infection. A regulatory network linked to 22 important pathways and 16 DEMs and 34 OGs was constructed. Some traditional immune-related signaling pathways such as toll-like receptor signaling pathway (TLR) signaling pathway, transforming growth factor-beta (TGF-beta) signaling pathway, peroxisome, phagosome, lysosome, mammalian target of rapamyoin (mTOR) signaling pathway were linked to specific miRNAs and genes in this network. Further, interactional relationship between certain miRNAs and TLR pathway was dissected, which the results predicted that a number of TLRs and TLR-associated signaling genes including TLR1, TLR2, TLR4, TLR6, IRAK1, TRAF6, MAPK, and IL-17 were negatively regulated by novel_miR_11, novel_miR_145, novel_miR_196, novel_miR_5, novel_miR_163 and novel_miR_217 in the TLR pathway. Additionally, interactional relationship between novel_miR_145 and TLR2 was validated by laboratory experiment. The integrated analysis of mRNA and microRNA deep sequencing data exhibited a sophisticated miRNA-mRNA regulation network in M. coruscus in response to V. alginolyticus challenge, which shed a new light on the underlying mechanism of molluscan confronting bacterial infection.

Hypoxia-challenged MSC-derived exosomes deliver miR-210 to attenuate post-infarction cardiac apoptosis

BACKGROUND

Myocardial infarction (MI) is a major cause of death worldwide. Although percutaneous coronary intervention and coronary artery bypass grafting can prolong life, cardiac damage persists. In particular, cardiomyocytes have no regenerative capacity. Mesenchymal stem cells (MSCs) are attractive candidates for the treatment of MI. The manner by which MSCs exert a beneficial effect upon injured cells is a source of continued study.

METHODS

After the isolation and identification of exosomes from MSCs, the expression of miR-210 was determined by microarray chip. Subsequently, gain- and loss-function approaches were conducted to detect the role of exosomes and exosomal-miR-210 in cell proliferation and apoptosis of cardiomyocytes, as well as the MI in vivo. Dual-Luciferase Report Gene System was used to demonstrate the target gene of miR-210.

RESULTS

We tested the hypothesis that MSC-derived exosomes transfer specific miRNA to protect cardiomyocytes from apoptotic cell death. Interestingly, direct cardiac injection of MSC exosomes reduced infarct size and improved heart function after coronary ligation. In vitro, the MSC exosomes enhanced cardiomyocyte survival to hypoxia. Confirmation of exosome uptake in myocytes was confirmed. Dual-luciferase reporter assay implicated miR-210 as a mediator of the therapeutic effect and AIFM3 as a downstream target. Treatment with miR-210 overexpressing MSC exosomes improved myocyte protection to both in vitro and in vivo stress. Furthermore, the endogenous and exogenous miR-210 had the same therapeutic effects.

CONCLUSION

These results demonstrated that the beneficial effects offered by MSC-exosomes transplantation after MI are at least partially because of excreted exosome containing mainly miR-210.

MicroRNAs involved in innate immunity regulation in the sea cucumber: A review

The sea cucumber is one of the most economically significant echinoderms. The immunity against exogenous stimulation of sea cucumber is of great academic and economic importance. MicroRNAs (miRNAs) are a class of short endogenous non-coding RNAs (ncRNAs) that are considered as vital regulators of both innate and adaptive immune responses in most eukaryotes. In sea cucumbers, some miRNAs (such as miR-133, miR-137, and miR-2008, among others) that participate in the regulation of innate immunity have been recently identified and characterized. This review focuses on those known miRNAs and their corresponding target genes that participate in the regulation of the complement system, Toll-like receptor (TLR) pathway, reactive oxygen species (ROS) production and apoptosis pathways in sea cucumbers. Moreover, we cover immune-related miRNA investigations in sea cucumbers that provide insights into developing more miRNA-based biomarkers and therapeutic strategies for sea cucumber diseases.

HBV infection suppresses the expression of inflammatory macrophage miR‑210

It has been previously reported that hepatitis B e‑antigen (HBeAg) induces microRNA (miR)‑155 expression and promotes liver injury by increasing inflammatory cytokine production in macrophages. Moreover, it was previously demonstrated that miR‑210 alleviates lipopolysaccharide‑stimulated proinflammatory cytokine production in macrophages. In addition, accumulating evidence suggests that miR‑210 is able to suppress hepatitis B virus (HBV) replication in HepG2.2.15 cells. However, it remains unclear whether miR‑210, similar to miR‑155, affects the progress of hepatitis B by regulating macrophage function. Reverse transcription‑quantitative polymerase chain reaction analysis was used to detect miR‑210 levels in serum and cells. HBV‑associated antigens stimulated different types of macrophages and facilitated the observation of the effects of these antigens on miR‑210 expression in macrophages. Co‑culture of peripheral blood monocytes from healthy controls and the serum of patients with chronic hepatitis B (CHB) was conducted to evaluate the effect of HBV‑associated elements in the serum on the expression of the macrophage miR‑210 in vivo. It was observed that miR‑210 expression levels were decreased in the peripheral blood monocytes (PBMs) and serum of patients with CHB and negatively associated with serum alanine aminotransferase and aspartate aminotransferase, but not other clinical parameters including hepatitis B surface antigen (HBsAg), HBeAg, anti‑HBe antibody (HBeAb) and hepatitis B core antibody (HBcAb) and HBV‑DNA. Notably, it was demonstrated that miR‑210 expression was not affected by treatment with HBV‑associated antigens in different types of macrophages. Notably, the serum of patients with CHB was able to markedly downregulate the miR‑210 expression of PBMs in healthy controls. These findings suggested that, unlike the induction of miR‑155 by HBeAg, there may be certain other elements, apart from HBV‑associated antigens, regulating miR‑210 levels in the serum and PBMs of patients with CHB that affect macrophage activation.

MicroRNAs Profiling Identifies miR-125a and Its Target Gene Wnt2 in Skins of Different Haired Rabbits

MicroRNAs (miRNAs) play critical roles in the control of skin and hair follicle development, epidermal homeostasis and pigmentation. However, the roles of miRNAs in the skins of rabbits with different hair types are unclear. In this study, we profiled miRNAs in the skins of long and short haired rabbits by Illumina deep sequencing. The dataset was compared with known mammalian miRNAs in miRBase 21.0. In total, 118 miRNAs were found to be differentially expressed between the two different rabbit types, of which 94 were upregulated, and 24 were downregulated in the skin of short haired vs. long haired rabbits. The expression levels of five randomly selected miRNAs detected by quantitative real-time PCR indicated that the expression patterns were consistent with Illumina sequencing results. What's more, bioinformatics analysis showed that miR-125a might target Wnt2, an important modulator for hair follicle development. To test whether Wnt2 is a target of miR-125a, luciferase reporter vector (pMir-report-Wnt2-3'-UTR-WT) and its substitution mutant (pMir-report-Wnt2-3'-UTR-MUT) were constructed. Co-transfection and reporter enzyme assays showed that compared with control (miR-125a NC transfection), miR-125a mimics transfection significantly inhibited the reporter luciferase activities expressed by pMir-report-Wnt2-3'-UTR-WT, while transfection of miR-125a inhibitors increased reporter enzyme activities. RT-PCR and Simple Western analysis found that Wnt2 mRNA and protein levels were induced or repressed by miR-125a overexpression or inhibition, respectively. Moreover, the mRNA expression levels of genes in Wnt signaling pathway, such as CTNNB1, LEF-1, PPARD and TGFB1, were also significantly changed (P < 0.05), consistent with Wnt2. It indicated that the regulation of Wnt2 expression by miRNAs may depend on the transcriptional degradation. The results will help to further understand the role of miRNAs in hair follicle development and the genetic mechanism underlying hair length phenotype.

A novel p38 MAPK gene in the sea cucumber Apostichopus japonicus (Ajp38) is associated with the immune response to pathogenic challenge

The p38 mitogen-activated protein kinase (MAPK), an important component of the MAPK signal cascade, is activated by extracellular stimuli, such as environmental stress and pathogenic infection. To clarify the function of p38 MAPKs in echinoderms, we used transcriptome database mining and rapid amplification of cDNA ends (RACE) to identify a novel p38 MAPK gene in the sea cucumber Apostichopus japonicus (here designated Ajp38). The full-length cDNA of Ajp38 was 2231 bp, including an open reading frame encoding 356 amino acid residues. Our sequence analysis indicated that the predicted Ajp38 protein contained the dual phosphorylation site Thr-Gly-Tyr (TGY) and was similar to the p38 homolog in sea urchins. Quantitative real-time PCR analysis showed that Ajp38 was ubiquitously expressed in all examined tissues of healthy adult A. japonicus, with the highest level of expression identified in the coelomocytes. Ajp38 mRNA expression was significantly upregulated in the coelomocytes 4, 12, and 72 h post in vivo infection with Vibrio splendidus. Our results provide more information about the characteristics and immune functions of the p38 homolog in sea cucumbers.

In-depth profiling of miRNA regulation in the body wall of sea cucumber Apostichopus japonicus during skin ulceration syndrome progression

MicroRNAs (miRNAs) are small non-coding RNAs that mediate mRNA degradation or translation repression. Previous study showed that the expression of miRNAs was significantly changed in the body wall of sea cucumber Apostichopus japonicus after skin ulceration syndrome (SUS) infection, which is a dynamic process. However, the critical miRNAs from body wall that involved in different infection stages of SUS remain unknown. In this study, four cDNA libraries were constructed with the body wall from healthy and three SUS-infected stages of A. japonicus. A total of 248 conserved miRNAs and five novel miRNAs were identified through Illumina HiSeq 2000 platform. Compared to the control, 238 miRNAs showed significant differential expression at three stages of SUS progression. Totally, 3149 miRNA-mRNA pairs were identified by target prediction and 314 miRNA-mRNA pairs showed negative correlation. It is noteworthy that 15 miRNAs and four mRNAs were located at the crucial positions of the network built with the anti-correlated miRNA-mRNA pairs. GO and KEGG enrichment analysis indicated that the predicted targets were involved in many immune-related processes. Deep analysis of miR-31c-5p, miR-29b-3p, NF-kB, mucin 2 and titin showed that they may play important roles in the pathogens attachment and recognition, signaling transduction and lesions repair of A. japonicus after SUS infection. These results would be useful for further investigating the potential roles of critical miRNAs and mRNAs in A. japonicus immune regulation.