Transcriptome Analysis and Discovery of Genes Involved in Immune Pathways from Coelomocytes of Sea Cucumber (Apostichopus japonicus) after Vibrio splendidus Challenge

The nuclear receptor coactivator 4 regulates ferritinophagy induced by vibrio splendidus in coelomocytes of Apostichopus japonicus

Iron homeostasis is vital for the host's defense against pathogenic invasion and the ferritinophagy is a crucial mechanism in maintaining intracellular iron homeostasis by facilitating the degradation and recycling of stored iron. The nuclear receptor coactivator 4 (NCOA4) serves as a ferritinophagy receptor, facilitating the binding and delivery of ferritin to the autophagosome and lysosome. However, NCOA4 of the sea cucumber Apostichopus japonicus (AjNCOA4) has not been reported until now. In this study, we identified and characterized AjNCOA4 in A. japonicus. This gene encodes a polypeptide containing 597 amino acids with an open reading frame of 1794 bp. The inferred amino acid sequence of AjNCOA4 comprises an ARA70 domain. Furthermore, a multiple sequence alignment demonstrated varying degrees of sequence homology between AjNCOA4 from A. japonicus and other NCOA4 orthologs. The phylogenetic tree of NCOA4 correlates with the established timeline of metazoan evolution. Expression analysis revealed that AjNCOA4 is expressed in all tested tissues, including the body wall, muscle, intestine, respiratory tree, and coelomocytes. Following challenge with Vibrio splendidus, the coelomocytes exhibited a significant increase in AjNCOA4 mRNA levels, peaking at 24 h. We successfully obtained recombinant AjNCOA4 protein through prokaryotic expression and prepared a specific polyclonal antibody. Immunofluorescence and co-immunoprecipitation experiments demonstrated an interaction between AjNCOA4 and AjFerritin in coelomocytes. RNA interference-mediated knockdown of AjNCOA4 expression resulted in elevated iron ion levels in coelomocytes. Bacterial stimulation enhanced ferritinophagy in coelomocytes, while knockdown of AjNCOA4 reduced the occurrence of ferritinophagy. These findings suggest that AjNCOA4 modulates ferritinophagy induced by V. splendidus in coelomocytes of A. japonicus.

Cloning and characterization of an aerolysin gene from a marine pathogen Vibrio splendidus

Vibrio splendidus is one of the main pathogens caused diseases with a diversity of marine cultured animals, especially the skin ulcer syndrome in Apostichopus japonicus. However, limited virulence factors have been identified in V. splendidus. In this study, one aerAVs gene coding an aerolysin of V. splendidus was cloned and conditionally expressed in Escherichia coli. The haemolytic activity of the recombinant AerAVs was analyzed. Western blotting was used to study of the secretion pathway of proaerolysin, and it showed that the proaerolysin was secreted via both outer membrane vehicles and classical secretion pathways. Since no active protein of aerolysin was obtained, one aerolysin surface displayed bacterium DH5α/pAT-aerA was constructed, and its haemolytic activity and virulence were determined. The results showed that the AerAVs displayed on the surface showed obvious haemolytic activity and cytotoxic to the coelomocyte of A. japonicus. Artificial immerse infection separately using the DH5α/pAT or DH5α/pAT-aerA was conducted. The result showed that the mortality percent of sea cucumber A. japonicus challenged with DH5α/pAT-aerA was 38.89 % higher than that challenged with the control strain DH5α/pAT, and earlier death occurred. Combined all the results indicates that aerolysin with the haemolytic activity and cytotoxic activity is a virulence factor of V. splendidus.

Vibrio splendidus AJ01 Promotes Pathogenicity via L-Glutamic Acid

Vibrio splendidus is a pathogen that infects a wide range of hosts, especially the sea cucumber species Apostichopus japonicus. Previous studies showed that the level of L-glutamic acid (L-Glu) significantly increased under heat stress, and it was found to be one of the best carbon sources used by V. splendidus AJ01. In this study, the effects of exogenous L-Glu on the coelomocyte viability, tissue status, and individual mortality of sea cucumbers were analyzed. The results showed that 10 mM of L-Glu decreased coelomocyte viability and increased individual mortality, with tissue rupture and pyknosis, while 0.1 mM of L-Glu slightly affected the survival of sea cucumbers without obvious damage at the cellular and tissue levels. Transcriptomic analysis showed that exogenous L-Glu upregulated 343 and downregulated 206 genes. Gene Ontology (GO) analysis showed that differentially expressed genes (DEGs) were mainly enriched in signaling and membrane formation, while a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that DEGs were significantly enriched in the upregulated endocytosis and downregulated lysosomal pathways. The coelomocyte viability further decreased by 20% in the simultaneous presence of exogenous L-Glu and V. splendidus AJ01 compared with that in the presence of V. splendidus AJ01 infection alone. Consequently, a higher sea cucumber mortality was also observed in the presence of exogenous L-Glu challenged by V. splendidus AJ01. Real-time reverse transcriptase PCR showed that L-Glu specifically upregulated the expression of the fliC gene coding the subunit protein of the flagellar filament, promoting the swimming motility activity of V. splendidus. Our results indicate that L-Glu should be kept in a state of equilibrium, and excess L-Glu at the host-pathogen interface prompts the virulence of V. splendidus via the increase of bacterial motility.

Ajpacifastin-like is involved in the immune response of Apostichopus japonicus challenged by Vibrio splendidus

Pacifastin proteins are previously found to regulate the phenoloxidase system in invertebrates and arthropods. In this study, the immune response that was regulated by Ajpacifastin-like in the sea cucumber Apostichopus japonicus was determined. RNA interference was used to knock down the expression of the Ajpacifastin-like gene in A. japonicus, followed by challenge with Vibrio splendidus, and the colony count showed that the survival of V. splendidus in the si-Ajpacifastin group increased 4.64-fold compared to that of the control group. The purified recombinant Ajpacifastin-like showed an inhibitory effect on the extracellular protease activity of the supernatant collected from the V. splendidus culture. Consequently, a comparative transcriptome analysis of the coelomocytes from the control group and the si-Ajpacifastin group was performed to explore the global regulatory effect of the Ajpacifastin-like. A total of 1486 differentially expressed genes (DEGs) were identified, including 745 upregulated genes and 741 downregulated genes. GO enrichment showed that the DEGs were mainly enriched in translation, cytosolic ribosomal subunit and structural constituent of ribosome. KEGG analysis showed that the DEGs were significantly enriched in the retinoic acid-inducible gene I (RIG-I)-like receptor signaling pathway, antigen processing and presentation, toll-like receptor signaling pathway, mitogen-activated protein kinase signaling pathway, nuclear factor-kappa B signaling pathway and other immune-related pathways. Furthermore, real-time reverse transcriptase PCR was used to determine the RNA levels of six DEGs, i.e., cathepsinB, CYLD, caspase8, TRAF6, hsp90 and FADD, to verify the RNA-seq results. Overall, our results specified the immune response and pathways of A. japonicus in which Ajpacifastin-like was involved in.

The lysosome-phagosome pathway mediates immune regulatory mechanisms in Mesocentrotus nudus against Vibrio coralliilyticus infection

Sea urchins are a popular model species for studying invertebrate diseases. The immune regulatory mechanisms of the sea urchin Mesocentrotus nudus during pathogenic infection are currently unknown. This study aimed to reveal the potential molecular mechanisms of M. nudus during resistance to Vibrio coralliilyticus infection by integrative transcriptomic and proteomic analyses. Here, we identified a total of 135,868 unigenes and 4,351 proteins in the four infection periods of 0 h, 20 h, 60 h and 100 h in M. nudus. In the I20, I60 and I100 infection comparison groups, 10,861, 15,201 and 8,809 differentially expressed genes (DEGs) and 2,188, 2,386 and 2,516 differentially expressed proteins (DEPs) were identified, respectively. We performed an integrated comparative analysis of the transcriptome and proteome throughout the infection phase and found very a low correlation between transcriptome and proteome changes. KEGG pathway analysis revealed that most upregulated DEGs and DEPs were involved in immune strategies. Notably, "lysosome" and "phagosome" activated throughout the infection process, could be considered the two most important enrichment pathways at the mRNA and protein levels. The significant increase in phagocytosis of infected M. nudus coelomocytes further demonstrated that the lysosome-phagosome pathway played an important immunological role in M. nudus resistance to pathogenic infection. Key gene expression profiles and protein‒protein interaction analysis revealed that cathepsin family and V-ATPase family genes might be key bridges in the lysosome-phagosome pathway. In addition, the expression patterns of key immune genes were verified using qRT‒PCR, and the different expression trends of candidate genes reflected, to some extent, the regulatory mechanism of immune homeostasis mediated by the lysosome-phagosome pathway in M. nudus against pathogenic infection. This work will provide new insights into the immune regulatory mechanisms of sea urchins under pathogenic stress and help identify key potential genes/proteins for sea urchin immune responses.

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.

Comparative analysis for immune response of coelomic fluid from coelom and polian vesicle in Apostichopus japonicus to Vibrio splendidus infection

The polian vesicle and coelom of sea cucumber Apostichopus japonicus were full of coelomic fluid in which many types of coelomocytes with different functions were suspended. Our previous work has indicated the differences of coelomocytes between two sites mainly in subtype proportion, non-specific immune enzymes activities and several immune-related genes expression levels in healthy A. japonicus. However, the functional similarities and differences of coelomic fluid in two sites including the coelom and polian vesicle after pathogenic infection still remain unclear. Here, we investigated the changes of the total coelomocyte density (TCD) and differential coelomocyte density (DCD) after pathogen infection by Vibrio splendidus in coelom and polian vesicle. After infected by V. splendidus, the TCD in the coelom and polian vesicle rapidly declined at 12 h, and then the TCD in the coelom showed a stably ascending trend, while the TCD in the polian vesicle reached a peak at 24 h post infection (hpi), and then showed a continuously decline trend from 24 hpi to 72 hpi followed by a slow elevation until recovering the normal level from 72 hpi to 96 hpi. Then the activities of acidic phosphatase (ACP), alkaline phosphatase (AKP), catalase (CAT) and superoxide dismutase (SOD) were determined to evaluate the response of cell-free coelomic fluid to V. splendidus infection. The activities of ACP, AKP and CAT showed similar trends in the coelom and polian vesicle. The SOD activity significantly increased in the polian vesicle, whereas it exhibited a decreasing trend in the coelom. Finally, the expression profiles of nine immune-related genes including Aj-MyD88, Aj-IRAK4, Aj-i-Lys, Aj-Rel, Aj-p50, Aj-DMBT1, Aj-CDC, Aj-Rrp15 and Aj-Fibrinogen C were detected after V. splendidus challenge. The results suggested all the detected genes were significantly up-regulated both in the coelom and polian vesicle, and the expression levels of these genes in two sites shared similar trends except Aj-MyD88 and Aj-DMBT1. This research provides a new insight into the differentially immune roles of coelomic fluid and coelomocytes in polian vesicle and coelom response to bacterial infections and supplements comprehensive resources for better understanding the innate immune response of A. japonicus.

Transcriptomic responses to air exposure stress in coelomocytes of the sea cucumber, Apostichopus japonicus

During rearing in hatcheries and transportation to restocking sites, sea cucumbers are often exposed to air for several hours, which may depress their non-specific immunity and lead to mass mortality. We performed transcriptome analysis of Apostichopus japonicus coelomocytes after air exposure to identify stress-related genes and pathways. After exposure to air for 1 h, individuals were re-submerged in aerated seawater and coelomocytes were collected at 0, 1, 4, and 16 h (B, H1, H4, and H16, respectively). We identified 6148 differentially expressed genes, of which 3216 were upregulated and 2932 were downregulated. Many genes involved in the immune response, antioxidant defense, and apoptosis were highly induced in response to air exposure. Enrichment analysis of Gene Ontology terms showed that the most abundant terms in the biological process category were oxidation-reduction process, protein folding and phosphorylation, and receptor-mediated endocytosis for the comparison of H1 vs. B, H4 vs. H1, and H16 vs. H4, respectively. Kyoto Eecyclopedia of Genes and Genomes enrichment analysis showed that six pathways related to the metabolism of proteins, fats, and carbohydrates were shared among the three comparisons. These results indicated that sea cucumbers regulate the expression of genes related to the antioxidant system and energy metabolism to resist the negative effects of air exposure stress. These findings may be applied to optimize juvenile sea cucumber production, and facilitate molecular marker-assisted selective breeding of an anoxia-resistant strain.

Transcriptome analysis of growth variation in early juvenile stage sandfish Holothuria scabra

The sandfish Holothuria scabra is a high-value tropical sea cucumber species representing a major mariculture prospect across the Indo-Pacific. Advancements in culture technology, rearing, and processing present options for augmenting capture production, stock restoration, and sustainable livelihood activities from hatchery-produced sandfish. Further improvements in mariculture production may be gained from the application of genomic technologies to improve performance traits such as growth. In this study, we performed de novo transcriptome assembly and characterization of fast- and slow-growing juvenile H. scabra from three Philippine populations. Analyses revealed 66 unigenes that were consistently differentially regulated in fast-growing sandfish and found to be associated with immune response and metabolism. Further, we identified microsatellite and single nucleotide polymorphism markers potentially associated with fast growth. These findings provide insight on potential genomic determinants underlying growth regulation in early juvenile sandfish which will be useful for further functional studies.

Identification and localization of growth factor genes in the sea cucumber, Holothuria scabra

The sea cucumber Holothuria scabra is both an economically important species in Asian countries and an emerging experimental model for research studies in regeneration and medicinal bioactives. Growth factors and their receptors are known to be key components that guide tissue repair and renewal, yet validation of their presence in H. scabra has not been established. We performed a targeted in silico search of H. scabra transcriptome data to elucidate conserved growth factor family and receptor genes. In total, 42 transcripts were identified, of which 9 were validated by gene cloning and sequencing. The H. scabra growth factor genes, such as bone morphogenetic protein 2A (BMP 2A), bone morphogenetic protein 5-like (BMP5-like), neurotrophin (NT) and fibroblast growth factor 18 (FGF18), were selected for further analyses, including phylogenetic comparison and spatial gene expression using RT-PCR and in situ hybridization. Expression of all genes investigated were widespread in multiple tissues. However, BMP 2A, BMP5-like and NT were found extensively in the radial nerve cord cells, while FGF18 was highly expressed in connective tissue layer of the body wall. Our identification and expression analysis of the H. scabra growth factor genes provided the molecular information of growth factors in this species which may ultimately complement the research in regenerative medicine.

Transcriptome analysis reveals roles of polian vesicle in sea cucumber Apostichopus japonicus response to Vibrio splendidus infection

Polian vesicle is originally regarded as a hematopoietic and inflammatory response organ in sea cucumber by the operations of cell depletion and heterogeneous cells injection, respectively. In the present study, to reveal the role and immune mechanisms of polian vesicle in response to pathogen, Vibrio splendidus, we first performed a comparative transcriptome analysis for the cells from polian vesicle wall in V. splendidus-challenged Apostichopus japonicus through RNA high-throughput sequencing technology. Briefly, 465,356,848 clean reads were obtained after cleaning up low-quality reads in total. Approximately 73% of the sequenced reads could be aligned to the reference genome of A. japonicus. The DEGs of CG (control group) vs TG 24 h (24 h post-infection group), CG vs TG 72 h (72 h post-infection group) and TG 24 h vs TG 72 h were 3762, 1391 and 3258, respectively. Gene Ontology (GO) annotation assay revealed that those genes associated with the processes such as cell process, cell, binding and catalytic activity were significantly induced in all three groups post V. splendidus infection. KEGG enrichment analysis suggested the DEGs in TG 24 h were enriched in Toll-like receptor (TLR) signaling pathway, complement and coagulation cascades, antigen processing and presentation and IL-17 signaling pathway compared with that in CG, while the pathways including ribosome biogenesis in eukaryotes, DNA replication, and cell cycle related with cell proliferation were mainly enriched in TG 72 h than that of CG. Furthermore, six important DEGs were chosen and showed the consistent expression patterns with the results of RNA-seq by qPCR. Overall, our analysis towards the current data demonstrates that polian vesicle may play an essential role in the regulation of immune response in A. japonicus and provide new insights into hematopoietic function of polian vesicle in response to pathogen infection.

mTORC2/Rictor is essential for coelomocyte endocytosis in Apostichopus japonicus

Endocytosis plays an important role in the immune defence systems of invertebrates through the interaction between the mechanical target of rapamycin complex 2 (mTORC2) and the AGC kinase family. Rictor is the most important unique subunit protein of mTORC2 and is thought to regulate almost all functions of mTORC2, including endocytosis. In the present study, a novel invertebrate Rictor homologue was identified from Apostichopus japonicus (designated as AjRictor) via the rapid amplification of cDNA ends (RACE). Spatial expression analysis indicated that AjRictor is ubiquitously expressed in all the examined tissues and has the highest transcript level in coelomocytes. Vibrio splendidus challenge in vivo and lipopolysaccharide (LPS) exposure in vitro could remarkably up-regulate the messenger RNA (mRNA) expression of AjRictor compared with the control group. AjRictor knockdown by 0.49- and 0.69-fold resulted in the significant decrease in endocytosis rate by 0.53- (P < 0.01) and 0.59-fold (P < 0.01) in vivo and in vitro compared with the control group, respectively. Similarly, the treatment of coelomocytes with rapamycin for 24 h and the destruction of the assembly of mTORC2 markedly decreased the endocytosis rate of the coelomocytes by 35.92% (P < 0.05). We detected the expression levels of endocytosis-related molecular markers after AjRictor knockdown and rapamycin treatment to further study the molecular mechanism between mTORC2 and endocytosis. Our results showed that AGC kinase family members (PKCα and Pan1) and the phosphorylation level of Akt^S473 were remarkably decreased after reducing mTORC2 activity; thus, mTORC2/Rictor plays a key role in the immune regulation of endocytosis in coelomocytes. Our current study indicates that mTORC2/Rictor is involved in the coelomocyte endocytosis of sea cucumber and plays an essential regulation role in defending pathogen invasion.

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.

Identification of molecular markers for superior quantitative traits in a novel sea cucumber strain by comparative microRNA-mRNA expression profiling

To investigate the adaptability of Apostichopus japonicus (A. japonicus) strain "Anyuan No. 1" in the South China Sea, field monitoring and microRNA-mRNA integrated analyses were conducted between "Anyuan No. 1" and a regular A. japonicus population from Wendeng (Shandong Province, as a control) in the Xiapu farming area in Fujian Province, China. The results showed that "Anyuan No. 1" exhibited greater body weight increase and a higher number of papillae compared to the control during two and a half months of field monitoring. Comparative microRNA (miRNA) and mRNA transcriptome analyses identified 12 differentially expressed miRNAs (DEMs) and 165 differentially expressed genes (DEGs) in "Anyuan No. 1" compared to the control. Long-chain specific acyl-CoA dehydrogenase (ACADL), transmembrane protein 251 (TMEM251), dehydrogenase/reductase SDR family protein 7-like (Dhrs7), insulin-like growth factor-binding protein 7 (IGFBP-7), CDK5 regulatory subunit-associated protein 1 (CDK5RAP1), visual pigment-like receptor peropsin, 39S ribosomal protein, miR-10, miR-153, miR-7, and miR-3529 were identified as gene and miRNA candidates correlated with superior economic traits in "Anyuan No. 1". Collectively, "Anyuan No. 1" is suitable for large-scale cultivation extension due to its better adaptability to the South China Sea area. Furthermore, we identified "miR10-ACADL" as a potential module for further molecular marker-assisted selective breeding of A. japonicus.

Transcriptomic analysis of sea cucumber (Holothuria leucospilota) coelomocytes revealed the echinoderm cytokine response during immune challenge

Background

The sea cucumber Holothuria leucospilota belongs to echinoderm, which is evolutionally the most primitive group of deuterostomes. Sea cucumber has a cavity between its digestive tract and the body wall that is filled with fluid and suspended coelomic cells similar to blood cells. The humoral immune response of the sea cucumber is based on the secretion of various immune factors from coelomocytes into the coelomic cavity. The aim of this study is to lay out a foundation for the immune mechanisms in echinoderms and their origins in chordates by using RNA-seq.

Results

Sea cucumber primary coelomocytes were isolated from healthy H. leucospilota and incubated with lipopolysaccharide (LPS, 10 μg/ml), polyinosinic-polycytidylic acid [Poly (I:C), 10 μg/ml] and heat-inactived Vibrio harveyi (10⁷ cell/ml) for 24 h, respectively. After high-throughput mRNA sequencing on an Illumina HiSeq2500, a de novo transcriptome was assembled and the Unigenes were annotated. Thirteen differentially expressed genes (DEGs) were selected randomly from our data and subsequently verified by using RT-qPCR. The results of RT-qPCR were consistent with those of the RNA-seq (R² = 0.61). The top 10 significantly enriched signaling pathways and immune-related pathways of the common and unique DEGs were screened from the transcriptome data. Twenty-one cytokine candidate DEGs were identified, which belong to 4 cytokine families, namely, BCL/CLL, EPRF1, IL-17 and TSP/TPO. Gene expression in response to LPS dose-increased treatment (0, 10, 20 and 50 μg/ml) showed that IL-17 family cytokines were significantly upregulated after 10 μg/ml LPS challenge for 24 h.

Conclusion

A de novo transcriptome was sequenced and assembled to generate the gene expression profiling across the sea cucumber coelomocytes treated with LPS, Poly (I:C) and V. harveyi. The cytokine genes identified in DEGs could be classified into 4 cytokine families, in which the expression of IL-17 family cytokines was most significantly induced after 10 μg/ml LPS challenge for 24 h. Our findings have laid the foundation not only for the research of molecular mechanisms related to the immune response in echinoderms but also for their origins in chordates, particularly in higher vertebrates.

Immune response of mollusk Onchidium struma to extremely low-frequency electromagnetic fields (ELF-EMF, 50 Hz) exposure based on immune-related enzyme activity and De novo transcriptome analysis

Along with rapid offshore and onshore wind power development in modern society, extremely low frequency electromagnetic fields (ELF-EMF) is produced extensively in the habits of aquatic organisms. However, the biological effects of ELF-EMF on aquatic organisms are almost sparse. In this study, Onchidium struma without shell was chosen to aim whether ELF-EMF can elicit immune response of mollusk based on immune-related enzyme activities and gene expression through high-throughput transcriptome sequencing. Three experimental groups, i.e. ELF-EMF unexposed control group (C), ELF-EMF (50 Hz, 100 μT) exposed E1 group, and ELF-EMF (50 Hz, 500 μT) exposed E2 group, were set, and coelomocytes were collected to analyze. The results showed that total coelomocyte and spherulocyte density in E1 group increased significantly compared to groups C and E2 (P < 0.05). There were no significant differences on amoebocyte and chromatocyte density among groups C, E1 and E2. ELF-EMF exposure could significantly increase immune-related enzyme activities in coelomic fluid of O. struma, including acidic phosphatase, alkaline phosphatase, antioxidative capacity, catalase, superoxide dismutase, and polyphenol oxidase (P < 0.05). A total of 54.32 Mb and 55.27 Mb raw reads with average length of 1520 bp were obtained from coelomocytes of O. struma in unexposed and exposed groups, respectively. There were 341 differentially expressed genes (DGEs) between unexposed and exposed groups, including 209 up-regulated and 132 down-regulated unigenes. All the DGEs were allocated to 14 Kyoto Encyclopedia of Genes and Genomes pathways, and five pathways were associated with immune response, including TLR/TNF/NOD-like receptor/MAPK/Fc epsilon RI signaling pathways. Altogether, short-term (to one week) exposure of O. struma to lower luxy density ELF-EMF (<500 μT) could elicit the immune response, and antioxidant system is recommended as indicators of immunological effects. Hopefully, this study will further provide insights into exploring biomarker for evaluation of the effect of ELF-EMF exposure on aquatic organisms regarding to field density, frequency and exposure duration, and provide good guidance for exploitation and utilization of renewable energy.

Genome-Wide Identification of Circular RNAs Revealed the Dominant Intergenic Region Circularization Model in Apostichopus japonicus

Circular RNAs (circRNAs) were recently recognized to act as competing endogenous RNAs and play roles in gene expression regulation. Previous studies in humans and silkworms have shown that circRNAs take part in immune regulation. Here, we conducted coelomocyte circRNA sequencing to explore its immune functions in healthy and skin ulceration syndrome (SUS)-diseased sea cucumbers. A total of 3,592 circRNAs were identified in libraries with diversified circularization patterns compared with animal models. The common intron-pairing-driven circularization models are not popular in sea cucumber genome, which was replaced with intergenic region circularization. The accuracy of these identified circRNAs was further validated by Sanger sequencing and RNase R-treated assays. Expression profile analysis indicated that 117 circRNAs were upregulated and 144 circRNAs were downregulated in SUS-diseased condition, of which 71.6% were intergenic-type circRNAs. The interaction network of differentially expressed circRNAs and microRNAs (miRNAs) was constructed and showed that miR-2008 and miR-31, detected with significantly differential expression in SUS-affected samples in a previous study, were predicted to be regulated by 10 and 11 differentially expressed circRNAs with more than 10 binding sites, respectively. Moreover, seven circRNAs were further validated by quantitative real-time PCR, whose variation trends were consistent with circRNA sequencing. All our results supported that intergenic-type circRNAs might have a dominant function in Apostichopus japonicas immune response by acting as miRNA regulators.

Divergent proteomics response of Apostichopus japonicus suffering from skin ulceration syndrome and pathogen infection

Skin ulceration syndrome (SUS) of sea cucumber is a common and serious disease that affects the stable development of Apostichopus japonicus in the culture industry. The part of sea cucumber that suffers from major injury and is directly observed is the body wall, in which protein variations should be the most direct evidence of the disease. To understand the response mechanisms of A. japonicus in SUS progression, we investigated protein changes in the body wall of diseased A. japonicus induced by Vibrio splendidus and individuals with natural diseases by isobaric tags for relative and absolute quantification (iTRAQ). About 119 proteins were identified in the two iTRAQ groups. A comparison of the protein expression profiles among two SUS conditions revealed that the mode of action induced by V. splendidus (Vs-SUS) was completely different from those in individuals with natural disease (ND-SUS). Most of the differentially expressed proteins (DEPs) (33 in 37 DEPs) were significantly depressed in the Vs-SUS group. Only 13 proteins in 27 DEPs showed similar trend to those in the ND-SUS group. Many important proteins involved in major intercellular signaling pathways associated with SUS disease were identified based on the KEGG and GO database search. Many proteins were located in the mitochondria and mainly involved in the oxidative stress pathway. Glutathione metabolism pathway was associated with reactive oxygen (ROS) production in the ND-SUS group. In the Vs-group, most of the proteins were concentrated in the cytoplasm and were related to immunity and extracellular matrix stability. In the ND-SUS group, the activity of key enzymes (CAT, GPx) that eliminate mitochondrial ROS production and structural stable protein (HSP60, HSP10) decreased, whereas those of complement proteins (C3, C3-2) that promoted ROS production was upregulated. This finding supported that oxidative damage caused by ROS might be the main effector for SUS in the ND-SUS group. The challenge with V. splendidus led to the breakdown of the defense capability of sea cucumber and suppressed the expression of immune-related proteins, such as C-type lectin, caspase, STAT, and cystatin. The downregulation of TIMP led to MMP1 overexpression. Members of the MMP family could directly degrade the extracellular matrix, which may be the main reason for the cell matrix degradation and induced SUS disease in the Vs-SUS group. Hence, ROS and extracellular matrix degradation enzymes could play important roles in the formation of SUS in sea cucumber. Results provide insights into the complex molecular mechanism of SUS in sea cucumber.

Integrative transcriptome analysis and discovery of genes involving in immune response of hypoxia/thermal challenges in the small abalone Haliotis diversicolor

In recent years, the abalone aquaculture industry has been threatened by the deteriorating environmental conditions, such as hypoxia and thermal stress in the hot summers. It is necessary to investigate the molecular mechanism in response to these environmental challenges, and subsequently understand the immune defense system. In this study, the transcriptome profiles by RNA-seq of hemocytes from the small abalone Haliotis diversicolor after exposure to hypoxia, thermal stress, and hypoxia plus thermal stress were established. A total of 103,703,074 clean reads were obtained and 99,774 unigenes were assembled. Of the 99,774 unigenes, 47,154 and 20,455 had homologous sequences in the Nr and Swiss-Prot protein databases, while 16,944 and 10,840 unigenes could be classified by COG or KEGG databases, respectively. RNAseq analysis revealed that the differentially expressed genes (DEGs) after challenges of hypoxia, thermal stress, or hypoxia plus thermal stress were 24,189, 29,165 and 23,665, among which more than 3000 genes involved in at least 230 pathways, including several classical immune-related pathways. The genes and pathways that were involved in immune response to hypoxia/thermal challenges were identified by transcriptome analysis and further validated by quantitative real-time PCR and RNAi technology. The findings in this study can provide information on H. diversicolor innate immunity to improve the abalone aquaculture industry, and the analysis of the potential immune-related genes in innate immunity signaling pathways and the obtained transcriptome data can provide an invaluable genetic resource for the study of the genome and functional genes.

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.

A novel fucolectin from Apostichopus japonicus with broad PAMP recognition pattern

F-type lectin (also known as fucolectin) is a newly identified family of fucose binding lectins with the sequence characters of a fucose binding motif and a unique lectin fold (the "F-type" fold). In the present study, a fucolectin was identified from sea cucumber Apostichopus japonicus (designated AjFL-1). The open reading frame (ORF) of AjFL-1 was of 546 bp, encoding a polypeptide of 181 amino acids with a predicted molecular mass of about 20 kDa. The deduced amino acid sequence of AjFL-1 shared 30%-40% similarity with the fucolectins from other animals. There were a typical F-type lectin domain (FLD) (residues 39-180) and a signal peptide (residues 1-24) in AjFL-1. The mRNA transcript of AjFL-1 could be detected by qRT-PCR in various tissues, such as intestinum, coelomocytes, respiratory tree, tentacle, and body wall, while undetectable in the gonads and longitudinal muscle. The mRNA expression level of AjFL-1 in coelomocytes was significantly up-regulated (47.06-fold to that in control group, p < 0.05) at 12 h after Vibrio splendidus challenge. Immunofluorescence assay showed that AjFL-1 protein was mainly distributed on the membrane, while few in cytoplasm of coelomocytes in sea cucumber. The recombinant AjFL-1 (rAjFL-1) could bind lipopolysaccharide (LPS), peptidoglycan (PGN), mannan (MAN) and fucose (FUC), and exhibited a broader binding activities towards Gram-negative bacterium Escherichia coli, Gram-positive bacterium Micrococcus luteus, as well fungus Pichia pastoris. In addition, rAjFL-1 could strongly promote the agglutination of fungus P. pastoris. These results indicated that AjFL-1 was a novel member of fucose-binding lectin family, which functioned as a pattern recognition receptor with broad spectrum of microbial recognition, and involved in innate immune response of sea cucumber.

Transcriptome analysis and discovery of genes involved in immune pathways from coelomocytes of Onchidium struma after bacterial challenge

Onchidium struma widely distributes in subtidal and low-tidal zones, which is considered to be an economical species with rich nutrition, a valuable biomonitor for heavy metal pollution and a representative species for evolution from ocean to land. However, there is limited genetic information available for O. struma development. This study compared transcriptomic profiles of coelomocytes from normal and bacteria infected O. struma by Illumina-based paired-end sequencing to explore the molecular immune mechanism of O. struma against bacterial infection. After assembly, a total of 92,450 unigenes with an average length of 1019 bp were obtained. Approximately 34,964 (37.82%) unigenes were annotated in the Nr NCBI database and 40.1% of unigenes were similar with that of Aplysia californica. Among them, 7609 unigenes were classified into three Gene Ontology (GO) categories: biological process (3250 unigenes, 42.7%), cellular component (2,281, 30.0%) and molecular function (2078 unigenes, 27.3%). A total of 22,776 unigenes were aligned to the Clusters of Orthologous Groups (COG) of proteins and classified into 25 functional categories. Following bacterial infection, 10,623 differently expressed unigenes (DEGs) were identified, including 7644 up-regulated and 2979 down-regulated unigenes. Further KEGG analysis annotated 11,681 DEGs to 42 pathways, and 11 pathways were identified to be related with diseases and immune system. To our knowledge, it was first time to analyze transcriptome profiles of O. struma. Results of the present study will provide valuable theoretical resources for future genetic and genomic research on O. struma. The research results will be helpful for improving the efficiency and quality of artificial breeding, establishing genetic linkage map, and enhancing health management for this species.

mcp, aer, cheB, and cheV contribute to the regulation of Vibrio alginolyticus (ND-01) adhesion under gradients of environmental factors

Adhesion is a key virulence factor of pathogens and can be affected by the environment. Our previously research with RNA-seq indicated that mcp, aer, cheB, and cheV might play roles in the regulation of adhesion in Vibrio alginolyticus (ND-01). In order to determine whether and how environmental factors affect adhesion through these genes, gene silencing was performed followed by quantitative real-time PCR (qRT-PCR), RNAi, transmission electron microscopy, and adhesion, capillary, and motility assays to verify how these genes influence adhesion. Silencing these genes led to deficiencies in adhesion, chemotaxis, flagellar assembly, and motility. The expression levels of cheA, cheW, and cheY, which are important genes closely related to the functions of mcp, aer, cheV, and cheB, were significantly downregulated in all of the RNAi groups. The expression of mcp, aer, cheV, and cheB under different gradients of temperature, pH, and salinity and after starvation for various durations was also detected, which showed that these genes were sensitive to certain environmental stresses, particularly pH and starvation. Our results indicated that mcp, aer, cheB, and cheV: (1) are necessary for ND-01 adhesion; (2) play key roles in the bacterial chemotaxis pathway by controlling the expression of downstream genes; (3) might affect adhesion by impacting motility, though motility is not the only route through which adhesion is affected; and (4) contribute to the regulation of ND-01 adhesion in natural environments with different temperatures, pH levels, and salinities as well as after various starvation periods.

De novo transcriptome assembly and developmental mode specific gene expression of Pygospio elegans

Species with multiple different larval developmental modes are interesting models for the study of mechanisms underlying developmental mode transitions and life history evolution. Pygospio elegans, a small, tube-dwelling polychaete worm commonly found in estuarine and marine habitats around the northern hemisphere, is one species with variable developmental modes. To provide new genomic resources for studying P. elegans and to address the differences in gene expression between individuals producing offspring with different larval developmental modes, we performed whole transcriptome Illumina RNA sequencing of adult worms from two populations and prepared a de novo assembly of the P. elegans transcriptome. The transcriptome comprises 66,233 unigenes, of which 33,807 contain predicted coding sequences, 26,448 have at least one functional annotation, and 3,076 are classified as putative long non-coding RNAs. We found more than 8,000 unigenes significantly differentially expressed between adult worms from populations producing either planktonic or benthic larvae. This comprehensive transcriptome resource for P. elegans adds to the available genomic data for annelids and can be used to uncover mechanisms allowing developmental variation in this and potentially other marine invertebrate species.

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.

Comparative transcriptome analysis of papilla and skin in the sea cucumber, Apostichopus japonicus

Papilla and skin are two important organs of the sea cucumber. Both tissues have ectodermic origin, but they are morphologically and functionally very different. In the present study, we performed comparative transcriptome analysis of the papilla and skin from the sea cucumber (Apostichopus japonicus) in order to identify and characterize gene expression profiles by using RNA-Seq technology. We generated 30.6 and 36.4 million clean reads from the papilla and skin and de novo assembled in 156,501 transcripts. The Gene Ontology (GO) analysis indicated that cell part, metabolic process and catalytic activity were the most abundant GO category in cell component, biological process and molecular funcation, respectively. Comparative transcriptome analysis between the papilla and skin allowed the identification of 1,059 differentially expressed genes, of which 739 genes were expressed at higher levels in papilla, while 320 were expressed at higher levels in skin. In addition, 236 differentially expressed unigenes were not annotated with any database, 160 of which were apparently expressed at higher levels in papilla, 76 were expressed at higher levels in skin. We identified a total of 288 papilla-specific genes, 171 skin-specific genes and 600 co-expressed genes. Also, 40 genes in papilla-specific were not annotated with any database, 2 in skin-specific. Development-related genes were also enriched, such as fibroblast growth factor, transforming growth factor-β, collagen-α2 and Integrin-α2, which may be related to the formation of the papilla and skin in sea cucumber. Further pathway analysis identified ten KEGG pathways that were differently enriched between the papilla and skin. The findings on expression profiles between two key organs of the sea cucumber should be valuable to reveal molecular mechanisms involved in the development of organs that are related but with morphological differences in the sea cucumber.

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

MicroRNAs (miRNAs) are small non-coding RNAs that are involved in many biological processes. To investigate the miRNAs related to skin ulceration syndrome (SUS) of Apostichopus japonicus, small RNA libraries of body wall, intestine, respiratory tree and coelomocytes from healthy and diseased A. japonicus were sequenced on Illumina Hiseq 2000 platform. A total of 247 conserved and 10 novel miRNAs were identified across all libraries. After pair-wise comparisons, 215 miRNAs in body wall, 36 in intestine, 2 in respiratory tree and 38 in coelomocytes showed significant expression differences. Further analyses were conducted on some tissue-specific differentially expressed miRNAs: miR-8 and miR-486-5p in body wall, miR-200-3p, let-7-5p and miR-125 in intestine, miR-278a-3p and bantam in respiratory, miR-10a and miR-184 in coelomocytes. Notably, these miRNAs in some species were reported to function in various physiological or pathological processes associated with immune regulations. Using stem-loop quantitative real time PCR, six representative miRNAs in four tissues were selected to validate the sequencing results. The Pearson's correlation coefficient (R) of the six miRNAs ranged from 0.777 to 0.948, which confirmed the consistency and accuracy between these two approaches. This study provides comprehensive expression and regulation patterns of functional miRNAs in different tissues and gives insights into the tissue-specific immune response mechanisms in SUS-infected A. japonicus.

Cathepsin L of the sea cucumber Apostichopus japonicus-molecular characterization and transcriptional response to Vibrio splendidus infection

Cathepsin L, a lysosomal endopeptidase, has been noted for its involvement in the innate immune response in invertebrates. Here, the cathepsin L cDNA of the sea cucumber Apostichopus japonicus (AjCatL) is identified from an EST library and then cloned by the rapid amplification of the cDNA ends (RACE) PCR. The full-length cDNA is 1678 bp long containing an open reading frame (ORF) of 1002 bp, an 80 bp 5' UTR and a 599 bp 3' UTR. The cDNA encodes 333 amino acid residues with a predicted molecular mass of 37.07 kDa and a theoretical isoelectric point (pI) of 5.01. The full-length AjCatL contains three active sites of eukaryotic thiol (cysteine) protease at positions 133-144, 278-288 and 295-314. Analysis of the predicted tertiary structure of prepro-CatL (17-333 aa) and mature-CatL (116-333 aa) reveals that the propeptide region (17-115 aa) blocks access to the substrate-binding cleft. Phylogenetic analysis shows that the AjCatL is clustered together with two other CatLs from Strongylocentrotus purpuratus. The enzymatic activity of AjCatL was verified using a substrate hydrolyzing assay with recombinant mAjCatL. Further analysis of real time-PCR demonstrates that the expression of AjCatL mRNA is significantly up-regulated in the coelomocytes in cases of infection with the common bacterial pathogen, Vibrio splendidus. This suggests that the AjCatL is likely to be involved in the immune response.