Expression analysis of microRNAs related to the skin ulceration syndrome of sea cucumber Apostichopus japonicus

MiR-210, not let-7, regulates Akt gene expression against Vibrio splendidus infection in the sea cucumber Apostichopus japonicus

The immune regulatory roles of microRNAs (miRNAs) have recently attracted considerable attention. Bioinformatics prediction revealed that both let-7 and miR-210 provide potential binding sites for the Akt (rac-alpha serine/threonine-protein kinase) gene sequence in the sea cucumber Apostichopus japonicus (termed AjAkt). In this study, we first used a dual-luciferase reporter assay and functional validation techniques to verify the interactions between these two miRNAs (let-7 and miR-210) and AjAkt, and then investigated the functions of the validated miRNA/mRNA pair as part of the innate immune response against Vibrio splendidus infection. We found that AjAkt interacts with miR-210 rather than let-7, and miR-210 negatively regulates the expression of AjAkt. From 8 to 48 h after infection with V. splendidus, opposite trends were observed in the expression levels of miR-210 and AjAkt (mRNA and protein) in coelomocytes, suggesting that the miR-210/AjAkt pair is involved in immune regulation during this period after infection. Both AjAkt silencing and miR-210 overexpression enhanced the phagocytic capacity and reduced the infectivity of A. japonicus after pathogen infection, suggesting that the miR-210/AjAkt pair may regulate the innate immune response of A. japonicus by altering phagocytic capacity. The findings of this study enrich our knowledge of the role of miRNA/mRNA pairs in immune regulation in sea cucumbers and provide insights into the molecular mechanisms of the innate immune response in marine echinoderms.

MicroRNA sequencing analysis reveals injury-induced immune responses of Scylla paramamosain against cheliped autotomy

During pond culture or intensive culture system of crabs (mainly Eriocheir sinensis, Portunus trituberculatus and Scylla paramamosain), high-density farming has typically contributed to a higher limb autotomy level in juvenile animals, especially in S. paramamosain which has a high level of cannibalism. Due to the high limb autotomy level, the survival and growth rates in S. paramamosain farming are restricted, which limit the growth of the mud crab farming industry. MicroRNAs (miRNAs) are small noncoding RNAs that regulate a series of biological processes including innate immune responses by post-transcriptional suppression of their target genes. MiRNAs are believed to be crucial for innate immune process of host wound healing. Many miRNAs have been verified to be required in host immune responses to repair wound and to defense pathogen after tissue damage. However, to our best knowledge, the miRNAs functions of crustacean innate immune reactions against injury induced by limb autotomy have not been studied yet. Here in this study, for the first time, miRNAs involved in the S. paramamosain immune reactions against injury induced by cheliped autotomy were obtained by high-throughput sequencing. A total of 575 miRNAs (518 known miRNAs and 57 novel predicted miRNAs) were obtained, of which 141 differentially expressed microRNAs (93 up-regulated microRNAs and 48 down-regulated microRNAs) were revealed to be modified against cheliped autotomy, and the qPCR results of randomly selected miRNAs confirmed the expression patterns in the miRNAs sequencing data. Numerous immune-related target genes associated with innate immune system were mediated by miRNAs to induce host humoral immune and cellular immune defense to minimize acute physical damage. Furthermore, the genes expression in hemolymph coagulation and melanization pathways, as well as Toll and Imd signaling pathways were mediated by miRNAs to activate host immune responses including melanization and antimicrobial peptides for rapid wound healing and killing invaded pathogens. These results will help to understand injury-induced immune responses in crabs and to develop an effective control strategy of autotomy rate in crabs farming.

Small RNA analysis of Perna viridis after exposure to Prorocentrum lima, a DSP toxins-producing dinoflagellate

Diarrheic shellfish poisoning toxins (DSP toxins) are a set of the most important phycotoxins produced by some dinoflagellates. Studies have shown that DSP toxins have various toxicities such as genotoxicity, cytotoxicity, and immunotoxicity to bivalve mollusks. However, these toxicities appear decreasing with exposure time and concentration of DSP toxins. The underlying mechanism involved remains unclear. In this study, small RNA sequencing was performed in the digestive gland of the mussel Perna viridis after exposure to DSP toxins-producing dinoflagellate Prorocentrum lima for different time periods. The potential roles of miRNAs in response and detoxification to DSP toxins in the mussel were analyzed. Small RNA sequencing of 12 samples from 72 individuals was conducted by BGISEQ-500. A total of 123 mature miRNAs were identified, including 90 conserved miRNAs and 33 potential novel miRNAs. After exposure to P. lima, multiple important miRNAs displayed some alterations. Further miRNA target prediction revealed some important genes involved in cytoskeleton, apoptosis, complement system and immune stress. qPCR demonstrated that miR-71_5, miR-750_1 and novel_mir4 were significantly up-regulated at 6 h after exposure to P. lima, while miR-100_2 was significantly down-regulated after 96 h of exposure. Accordingly, putative target genes of these differentially expressed miRNAs experienced some changes. After 6 h of DSP toxins exposure, NHLRC2 and C1q-like were significantly down-regulated. After 96 h of DSP toxins exposure, NHLRC2 was significantly up-regulated. It is reasonable to speculate that the mussel P. viridis might respond to DSP toxins through miR-750_1, novel_mir4 and miR-71_5 regulating the expression of relevant target genes involved in apoptosis, cytoskeleton, and immune response, etc. This study might provide new clues to uncover the toxic response of bivalve to DSP toxins and lay a foundation for revealing the roles of miRNAs in the environmental adaptation in shellfish.

Comparative Mucous miRomics in Cynoglossus semilaevis Related to Vibrio harveyi Caused Infection

Epidermal mucus is an important barrier and regulating mediator in fish. MicroRNAs (miRNAs) are proved to be involved in various biological processes, also as promising biomarkers for disease diagnosis. Vibrio harveyi has long been a noticeable bacterial pathogen in Cynoglossus semilaevis aquaculture. To find the evidence whether there are indicating miRNAs in mucus and whether the miRNAs are related to infections caused by V. harveyi, miRNA profiles of mucus from V. harveyi infected fish and healthy controls were screened by small RNA sequencing and verified by quantitative real-time PCR. This is the first report about miRNA profiling of flatfish mucus, aiming at illustrating the pathogenesis of V. harveyi caused infection and developing disease-related biomarkers. The results revealed significant differences in expression levels of some miRNAs between infected fish and healthy ones. Three hundred differentially expressed miRNAs were obtained after filtering through FC > 2 or FC < 0.5 and most of the differential miRNAs were downregulated. After verification through qRT-PCR, four unique miRNAs, dre-miR-451, dre-miR-184, dre-miR-205-5p > ssa-miR-205b-5p, and dre-miR-181a-5p > ssa-miR-181a-5p, were identified as V. harveyi infection-related signatures, consistent with sequencing trend. The expression levels of these four miRNAs in the infected fish were all significantly lower than controls. These miRNAs in mucus could be used to differentiate diseased and healthy fish in a non-invasive way with practical value for large-scale disease screening. They also provided new insights into the mechanism underlying the bacterial infections in fish.

Sea cucumbers in a high temperature and low dissolved oxygen world: Roles of miRNAs in the regulation of environmental stresses

The exacerbation of global warming has driven changes in environmental factors, including water temperature and oxygen concentration. The sea cucumber Apostichopus japonicus, an economically important aquatic animal, is constantly and directly challenged by heat and hypoxia. In this study, 12 small RNA libraries were constructed for this species, and a total of 21, 26 and 22 differentially expressed (DE) miRNAs were clarified in A. japonicus under thermal (26 °C), hypoxic (2 mg/L) and the combined stresses. Comparative miRNA sequencing analysis and real-time PCR were used to identify and validate the representative miRNAs, including Aja-miR-novel-299, Aja-let-7b-3p, Aja-miR-71b-5p, Aja-miR-novel-13218 and Aja-miR-2004 in response to high temperature, and Aja-miR-92b-3p, Aja-miR-210-5p and Aja-miR-novel-26331 in response to oxygen limitation. GO and KEGG pathway analysis revealed that the potential target genes of DE-miRNAs involved in biosynthesis, metabolism, immunity, cell growth and death, translation and signaling transduction. Key DE-miRNAs with potentially targeted genes associated with heat shock and hypoxia response were also determined. These results may help explaining the role of miRNA regulation in stress resistance, as well as the potential molecular regulation mechanism of the echinoderm A. japonicus in the context of global warming.

Detecting Cynoglossus semilaevis infected with Vibrio harveyi using micro RNAs from mucous exosomes

Exosomes are important mediators of vesicle transportation and contain microRNAs (miRNAs) that mediate transcriptional gene knockout and silencing in biological processes. Moreover, exosomic miRNAs are promising biomarkers for disease diagnosis and physiological status indication in many species, including fish. The impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis aquaculture is becoming more and more serious as the industry expands. To overcome this challenge, miRNAs in mucous exosomes were screened by small RNA sequencing and verified by quantitative real-time PCR to develop biomarkers. This is the first capture of exosomes from flatfish mucus coupled with miRNA profiling. The results revealed significant differences in expression levels of some miRNAs between infected and healthy fish. Three unique miRNAs were identified for V. harveyi infection diagnosis; expression levels of dre-miR-205-5p and dre-miR-205-5p in infected fish were significantly lower than controls, while dre-miR-100-5p expression was higher. These miRNAs in mucous exosomes could be used to differentiate diseased and healthy fish in an early screening method with practical value for breeding disease-resistant C. semilaevis.

Differential expression of microRNAs in mud crab Scylla paramamosain in response to white spot syndrome virus (WSSV) infection

Till date numerous microRNAs (miRNAs) have been discovered from various organisms, including fish, shellfish and crustaceans. The miRNAs are known to regulate immune functions in crustaceans, but little is known about the role of miRNAs against viral infection in crab. We performed small RNA sequencing to characterize the differentially expressed miRNAs in WSSV infected Scylla paramamosain, in comparison to that in uninfected crab, at 2 h and 12 h post infection. In total, 24 host miRNAs were up-regulated and 25 host miRNAs were down-regulated in response to WSSV infection at 2 h post infection. And 27 host miRNAs were up-regulated and 30 host miRNAs were down-regulated in response to WSSV infection at 12 h post infection. Further, the gene ontology analysis revealed that many signaling pathways were mediated by these miRNAs. The integral component of membrane is the most important biological process and endocytosis pathway is the most important pathway, which indicates that endocytosis is very important for WSSV infection. This study is one important attempt at characterizing crab miRNAs that response to WSSV infection, and will help unravel the miRNA pathways involved in antiviral immunity of crab.

Insights into the DNA methylation of sea cucumber Apostichopus japonicus in response to skin ulceration syndrome infection

DNA methylation is an important epigenetic modification that regulates gene expression in many biological processes, including immune response. In this study, whole-genome bisulfite sequencing (WGBS) was carried out on healthy body wall (HB) and skin ulceration syndrome (SUS) infected body wall (SFB) to gain insights into the epigenetic regulatory mechanism in sea cucumber Apostichopus japonicus. After comparison, a total of 116,522 differentially methylated regions (DMRs) were obtained including 67,269 hyper-methylated and 49,253 hypo-methylated DMRs (p < 0.05, FDR < 0.001). GO enrichment analysis indicated that regulation of DNA-templated transcription (GO: 0006355), where DNA methylation occurred, was the most significant term in the biology process. The integration of methylome and transcriptome analysis revealed that 10,499 DMRs were negatively correlated with 496 differentially expressed genes (DEGs). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that these DEGs were enriched in the phosphoinositide 3-kinase-protein kinase B (PI3K/Akt)/mammalian target of rapamycin (mTOR) signaling pathway. Interestingly, two serine/threonine-protein kinases, nemo-like kinase (NLK) and mTOR, were highlighted after functional analysis. The variations of methylation in these two genes were associated with SUS infection and immune regulation. They regulated gene expression at different levels and showed interaction during response process. The validation of methylation sites showed high consistency between pyrosequencing and WGBS. WGBS analysis not only revealed the changes of DNA methylation, but also presented important information about the regulation of key genes after SUS infection in A. japonicus.

Molecular characterization and functional analysis of MyD88 from the tropical sea cucumber, Holothuria leucospilota

In this study, a myeloid differentiation factor 88 (MyD88) named as HLMyD88 was identified from the sea cucumber Holothuria leucospilota. The full-length cDNA of HLMyD88 is 4797 bp in size, containing a 227 bp 5'-untranslated region (UTR), a 3721 bp 3'-UTR and an 849 bp open reading frame (ORF) encoding a protein of 282 amino acids with a deduced molecular weight of 32.25 kDa HLMyD88 contains an N-terminal death domain and a C-terminal Toll/interluekin-1 receptor (TIR) domain with three highly conserved sequence motifs named as Box 1, Box 2 and Box 3. The results of luciferase reporter assay showed that over-expressed HLMyD88 in HEK293T cells could activate the transcription factors nuclear factor-κB (NF-κB) and activator protein 1 (AP-1). Additionally, the secretion of proinflammatory cytokines IL-1β and TNF-α in the HEK293T cells was increased by over-expressed HLMyD88, indicating the potential role of HLMyD88 in the innate immunity of sea cucumber. Moreover, we further confirmed that over-expressed HLMyD88 could also induce apoptosis.

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.

Modulation of Atlantic salmon miRNome response to sea louse infestation

MicroRNAs are non-coding RNA that plays a crucial role in post-transcriptional regulation and immune system regulation. On other hand, sea lice are prevalent parasites that affect salmon farming, generating different degrees of immune suppression depending on the salmon and sea louse species. Caligus rogercresseyi for example, which affects the salmon industry in Chile, decreases Th1 response, macrophage activation, TLR-mediated response and iron regulation in infected fish. In this study, we explore Atlantic salmon miRNome during infestation by C. rogercresseyi. Using small RNA sequencing, we annotated 1718 miRNAs for skin and head kidney from infected Atlantic salmon. The most abundant families identified were mir-10, mir-21, mir-30, mir-181 and let7. Significant differences were found between tissue, with 1404 annotated miRNA in head kidney and 529 in skin. Differential analysis of transcript expression indicated that at an early stage of infestation miRNA expression was higher in head kidney than in skin tissue, revealing tissue-specific expression patterns. In parallel, miRNA target prediction using 3'UTRs from highly regulated immune-related genes and iron metabolism showed that mir-140-4 and mir-181a-2-5 modulate the expression of TLR22 and Aminolevulinic acid synthase, respectively. This study contributes knowledge about the immune response of Atlantic salmon during infestation with sea lice.

Profiling of cellular microRNA responses during the early stages of KSHV infection

Kaposi's sarcoma-associated herpesvirus (KSHV) causes a variety of cancers, including Kaposi's sarcoma (KS), primary effusion lymphoma (PEL), and multicentric Castleman disease (MCD). Host cellular microRNAs (miRNAs) play important post-transcriptional regulatory roles in gene expression and can greatly influence virus-host cell interactions. This study investigated cellular miRNA expression profiles operating in response to early stages of KSHV infection of human Burkitt lymphoma B cells (BJAB). We employed deep sequencing to analyze miRNA expression in KSHV-infected BJAB cells 15 min post infection (PI) and compared this to uninfected BJAB cells. A total of 32 known miRNAs and 28 novel miRNA candidates were differentially expressed in KSHV-infected compared to uninfected BJAB cells. Interestingly, miRNA expression profiles during the early stages of viral infection yielded comparable results when UV-inactivated KSHV was used. The deep sequencing results were further confirmed by performing real-time reverse transcription PCR. The target genes predicted to be regulated by both the known and novel miRNAs are mainly involved in assisting virus entry, inducing critical cell signaling, initiating transcription of immediate early genes, promoting latent infection, and modulating the host immune response. For the first time, we provide insight into the host cellular miRNA expression profiles in response to early stages of KSHV infection of human B cells. Furthermore, this study offers a valuable basis for further investigation on the roles of cellular miRNAs in the KSHV entry process.

Identification and Characterization of MicroRNAs in the Liver of Blunt Snout Bream (Megalobrama amblycephala) Infected by Aeromonas hydrophila

MicroRNAs (miRNAs) are small RNA molecules that play key roles in regulation of various biological processes. In order to better understand the biological significance of miRNAs in the context of Aeromonas hydrophila infection in Megalobrama amblycephala, small RNA libraries obtained from fish liver at 0 (non-infection), 4, and 24 h post infection (poi) were sequenced using Illumina deep sequencing technology. A total of 11,244,207, 9,212,958, and 7,939,157 clean reads were obtained from these three RNA libraries, respectively. Bioinformatics analysis identified 171 conserved miRNAs and 62 putative novel miRNAs. The existence of ten randomly selected novel miRNAs was validated by RT-PCR. Pairwise comparison suggested that 61 and 44 miRNAs were differentially expressed at 4 and 24 h poi, respectively. Furthermore, the expression profiles of nine randomly selected miRNAs were validated by qRT-PCR. MicroRNA target prediction, gene ontology (GO) annotation, and Kyoto Encylopedia of Genes and Genomes (KEGG) analysis indicated that a variety of biological pathways could be affected by A. hydrophila infection. Additionally, transferrin (TF) and transferrin receptor (TFR) genes were confirmed to be direct targets of miR-375. These results will expand our knowledge of the role of miRNAs in the immune response of M. amblycephala to A. hydrophila infection, and facilitate the development of effective strategies against A. hydrophila infection in M. amblycephala.

RNA sequencing analysis to capture the transcriptome landscape during skin ulceration syndrome progression in sea cucumber Apostichopus japonicus

Background

Sea cucumber Apostichopus japonicus is an important economic species in China, which is affected by various diseases; skin ulceration syndrome (SUS) is the most serious. In this study, we characterized the transcriptomes in A. japonicus challenged with Vibrio splendidus to elucidate the changes in gene expression throughout the three stages of SUS progression.

Results

RNA sequencing of 21 cDNA libraries from various tissues and developmental stages of SUS-affected A. japonicus yielded 553 million raw reads, of which 542 million high-quality reads were generated by deep-sequencing using the Illumina HiSeq™ 2000 platform. The reference transcriptome comprised a combination of the Illumina reads, 454 sequencing data and Sanger sequences obtained from the public database to generate 93,163 unigenes (average length, 1,052 bp; N50 = 1,575 bp); 33,860 were annotated. Transcriptome comparisons between healthy and SUS-affected A. japonicus revealed greater differences in gene expression profiles in the body walls (BW) than in the intestines (Int), respiratory trees (RT) and coelomocytes (C). Clustering of expression models revealed stable up-regulation as the main pattern occurring in the BW throughout the three stages of SUS progression. Significantly affected pathways were associated with signal transduction, immune system, cellular processes, development and metabolism. Ninety-two differentially expressed genes (DEGs) were divided into four functional categories: attachment/pathogen recognition (17), inflammatory reactions (38), oxidative stress response (7) and apoptosis (30). Using quantitative real-time PCR, twenty representative DEGs were selected to validate the sequencing results. The Pearson’s correlation coefficient (R) of the 20 DEGs ranged from 0.811 to 0.999, which confirmed the consistency and accuracy between these two approaches.

Conclusions

Dynamic changes in global gene expression occur during SUS progression in A. japonicus. Elucidation of these changes is important in clarifying the molecular mechanisms associated with the development of SUS in sea cucumber.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2810-3) contains supplementary material, which is available to authorized users.