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

Sea cucumbers: an emerging system in evo-devo

A challenge for evolutionary developmental (evo-devo) biology is to expand the breadth of research organisms used to investigate how animal diversity has evolved through changes in embryonic development. New experimental systems should couple a relevant phylogenetic position with available molecular tools and genomic resources. As a phylum of the sister group to chordates, echinoderms extensively contributed to our knowledge of embryonic patterning, organ development and cell-type evolution. Echinoderms display a variety of larval forms with diverse shapes, making them a suitable group to compare the evolution of embryonic developmental strategies. However, because of the laboratory accessibility and the already available techniques, most studies focus on sea urchins and sea stars mainly. As a comparative approach, the field would benefit from including information on other members of this group, like the sea cucumbers (holothuroids), for which little is known on the molecular basis of their development. Here, we review the spawning and culture methods, the available morphological and molecular information, and the current state of genomic and transcriptomic resources on sea cucumbers. With the goal of making this system accessible to the broader community, we discuss how sea cucumber embryos and larvae can be a powerful system to address the open questions in evo-devo, including understanding the origins of bilaterian structures.

A novel short-type peptidoglycan recognition protein with unique polysaccharide recognition specificity in sea cucumber, Apostichopus japonicus

Pattern recognition receptors (PRRs) are the first line of immune defense in invertebrates against pathogen infection; they recognize pathogens and transmit signals to downstream immune pathways. Among these, peptidoglycan recognition proteins (PGRPs) are an important family in invertebrates that generally comprise of complicated isoforms. A comprehensive understanding of PGRPs in evolutionarily and economically important marine invertebrates, such as the sea cucumber, Apostichopus japonicus, is crucial. Previous studies have identified two PGRPs in sea cucumber, AjPGRP-S and AjPGRP-S1, and another novel short-type PGRP, AjPGRP-S3, was additionally identified here. The full-length cDNA sequence of AjPGRP-S3 was obtained here by PCR-RACE, followed by which showed its gene expression analyses by in situ hybridization that showed it to be relatively highly expressed in coelomocytes and tube feet. Based on an analysis of the recombinant protein, rAjPGRP-S3, a board-spectrum pathogen recognition ability was noted that covered diverse Gram-negative and -positive bacteria, and fungi. Moreover, according to the results of yeast two-hybridization, it was suggested that rAJPGRP-S3 interacted with multiple immune-related factors, including proteins involved in the complement system, extracellular matrix, vesicle trafficking, and antioxidant system. These findings prove the important functions of AjPGRP-S3 in the transduction of pathogen signals to downstream immune effectors and help explore the functional differences in the AjPGRP isoforms.

Transcriptomic Responses of a Lightly Calcified Echinoderm to Experimental Seawater Acidification and Warming during Early Development

Ocean acidification (OA) and ocean warming (OW) are potential obstacles to the survival and growth of marine organisms, particularly those that rely on calcification. This study investigated the single and joint effects of OA and OW on sea cucumber Apostichopus japonicus larvae raised under combinations of two temperatures (19 °C or 22 °C) and two pCO2 levels (400 or 1000 μatm) that reflect the current and end-of-21st-century projected ocean scenarios. The investigation focused on assessing larval development and identifying differences in gene expression patterns at four crucial embryo-larval stages (blastula, gastrula, auricularia, and doliolaria) of sea cucumbers, using RNA-seq. Results showed the detrimental effect of OA on the early development and body growth of A. japonicus larvae and a reduction in the expression of genes associated with biomineralization, skeletogenesis, and ion homeostasis. This effect was particularly pronounced during the doliolaria stage, indicating the presence of bottlenecks in larval development at this transition phase between the larval and megalopa stages in response to OA. OW accelerated the larval development across four stages of A. japonicus, especially at the blastula and doliolaria stages, but resulted in a widespread upregulation of genes related to heat shock proteins, antioxidant defense, and immune response. Significantly, the negative effects of elevated pCO2 on the developmental process of larvae appeared to be mitigated when accompanied by increased temperatures at the expense of reduced immune resilience and increased system fragility. These findings suggest that alterations in gene expression within the larvae of A. japonicus provide a mechanism to adapt to stressors arising from a rapidly changing oceanic environment.

Bacillus licheniformis reduces mortality, improves symptoms, and controls residual pathogenicity in juvenile sea cucumbers Apostichopus japonicus infected with skin ulceration syndrome

Mass mortality caused by skin ulceration syndrome (SUS) is the bottle-neck for the sustainable aquaculture of the sea cucumber Apostichopus japonicus. In the present study, probiotic Bacillus licheniformis (0.25 × 109 CFU/g) was used as the treatment for A. japonicus infected with the SUS that caused by Vibrio harveyi. We found that B. licheniformis significantly reduced the number of infected sea cucumbers 5 days and 7 days after the treatment (group B), compared to those without B. licheniformis treatment (group C) (P < 0.001; P < 0.001). Further, the sea cucumbers fed B. licheniformis had significantly lower mortality at the end of the experiment (<10%) than that of those without the B. licheniformis treatment (>60%) (P < 0.001). These results suggest that the treatment of B. licheniformis is an effective method to reduce the mass mortality resulted from SUS in sea cucumber aquaculture. Further, 3-5 days of treatment significantly improved the adverse symptoms of SUS on the physiology and behavior of sea cucumbers, including the righting behavior, adhesion behavior, food consumption, fecal output and mobility. This indicates B. licheniformis treatment has the advantage in the recovery of sea cucumbers after SUS. Moreover, there was no significant difference observed in the physiology and behavior of sea cucumbers between the SUS infected sea cucumbers after the 7-day treatment of B. licheniformis and the healthy individuals. SUS infected sea cucumbers effectively returned to a stage of normalcy. Further, we found a significantly lower infected rate in sea cucumbers exposed to the culture water of group B (∼5%) than that of those in exposure to the culture water of group C (∼60%). This indicates that the treatment of B. licheniformis efficiently controls the residual pathogenicity of SUS in culture water. The present study demonstrated the effectiveness of B. licheniformis treatment as an environmentally friendly approach to reducing mortality, improving symptoms, and controlling residual pathogenicity in sea cucumber aquaculture.

Integrative mRNA-miRNA interaction analysis associated with the immune response of Strongylocentrotus intermedius to Vibrio harveyi infection

Strongylocentrotus intermedius is one of the most economically valuable sea urchin species in China and has experienced mass mortality owing to outbreaks of bacterial diseases such as black mouth disease. This has caused serious economic losses to the sea urchin farming industry. To investigate the immune response mechanism of S. intermedius with different tube feet colors in response to Vibrio harveyi infection, we examined the different tube feet-colored S. intermedius under V. harveyi challenge and compared their transcriptome and microRNA (miRNA) profiles using RNA-Seq. We obtained 1813 differentially expressed genes (DEGs), 28 DE miRNAs, and 303 DE miRNA-DEG pairs in different tube feet-colored S. intermedius under V. harveyi challenge. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the most significant DEGs were associated with the Notch signaling and phagosome pathways. The target genes of immune-related miRNAs (miR-71, miR-184, miR-193) and genes (CALM1, SPSB4, DMBT, CSRP1) in S. intermedius were predicted and validated. This study provides insight into the molecular mechanisms that regulate genes involved in the immune response of S. intermedius infected with V. harveyi.

Functions of lysin motif (LysM)-containing protein in antibacterial responses of sea cucumbers, Apostichopus japonicus

The lysin motif (LysM)-containing protein is one of widespread pattern-recognition receptors in prokaryotes and eukaryotes. Numerous LysM-containing gene sequences are present in gene databases; however, few have been well characterized, especially in echinoderms. In this study, the full-length cDNA of a novel LysM-containing gene was obtained from the sea cucumber Apostichopus japonicus, named AjLysM-1, using polymerase chain reaction (PCR) combined with rapid amplification of cDNA ends. We prepared and expressed recombinant AjLysM-1 protein (rAjLysM-1) and determined its pathogen-recognition ability by enzyme-linked immunosorbent and immunofluorescence assays. We also analyzed the tissue expression pattern and response to immune challenges of AjLysM-1 using quantitative real-time reverse transcription-PCR and in situ hybridization. The AjLysM-1 protein was predicted to be an intracellular non-secreted LysM-containing protein, highly homologous to the same protein in other marine echinoderms. AjLysM-1 transcripts were highest expressed in coelomocytes and were strikingly induced by challenge with representative bacterial and fungal polysaccharides. rAjLysM-1 showed weak binding to mannan, Pseudoalteromonas nigrifaciens, and Shewanella baltica, implying that AjLysM-1 might provide inadequate defense against Gram-negative bacteria and fungi. Notably, rAjLysM-1 also interacted with tyrosine protein kinase and filamin-B, indicating that it could be involved in focal adhesion in A. japonicus. These findings improve our understanding of the functions of LysM-containing proteins in marine echinoderms.

Genome-wide identification m6A modified circRNAs revealed their key roles in skin ulceration syndrome disease development in Apostichopus japonicus

Circular RNAs (circRNAs) are novel endogenous non-coding RNAs (ncRNAs) and can be acted as competing endogenous RNAs (ceRNAs) to regulate microRNA (miRNA) and downstream gene expression. Recently, m⁶A modification has been found in circRNA, and m⁶A circRNAs also play important roles in various biological processes and a variety of diseases. Our previous study had been demonstrated that circRNAs were differentially expressed in skin ulceration syndrome (SUS) diseased sea cucumber Apostichopus japonicus. However, whether the function of circRNAs are dependent on m⁶A levels are largely unknown. Here, we firstly investigated the genome-wide map of m⁶A circRNAs in sea cucumbers with different stages of Vibrio splendidus challenge, that's Control group, SUS-diseased group, and SUS-resistant group. MeRIP-seq revealed that m⁶A abundances were enriched in circRNAs in all three groups, especially for SUS-resistant group. Among them, more than 62% of modified circRNAs harbor only a single m⁶A peak and about 55% of m⁶A sites in circRNAs were derived from sense overlapping in each group. After V. splendidus infection, we found that most of m⁶A peaks in circRNAs were upregulated and less were downregulated in both SUS-diseased and SUS-resistant groups when compared with Control. Furthermore, GO analysis indicated that the host genes of circRNAs with dysregulated m⁶A peaks in SUS-diseased and SUS-resistant groups were both mainly enriched in the adhesion pathway. More importantly, we discovered that more than 50% m⁶A circRNAs showed a positive correlation between the circRNAs expression and m⁶A methylation levels both in SUS-diseased and SUS-resistant groups. Therefore, a core circRNA-miRNA-mRNA (ceRNA) network whether influenced by m⁶A modification was constructed based on conjoint analysis. Our results indicated that several selected m⁶A circRNAs bind with miRNAs were mainly targeting to ubiquitylation system and adhesion pathway. What's more, three candidate m⁶A circRNAs and three target genes were validated by MeRIP-qPCR and qPCR, whose m⁶A levels in circRNA and mRNA expressions were consistent with disease occurrence or disease resistance. All of our current findings suggested that m⁶A circRNAs could play important roles during pathogen infection and might be served as a new molecular biomarker in SUS disease diagnose of A. japonicus.

Sea Cucumber Body Vesicular Syndrome Is Driven by the Pond Water Microbiome via an Altered Gut Microbiota

Apostichopus japonicus (sea cucumber) is one of the most valuable aquaculture species in China; however, different diseases can limit its economic development. Recently, a novel disease, body vesicular syndrome (BVS), was observed in A. japonicus aquaculture. Diseased animals displayed no obvious phenotypic characteristics; however, after boiling at the postharvest stage, blisters, lysis, and body ruptures appeared. In this study, a multiomics strategy incorporating analysis of the gut microbiota, the pond microbiome, and A. japonicus genotype was established to investigate BVS. Detailed analyses of differentially expressed proteins (DEPs) and metabolites suggested that changes in cell adhesion structures, caused by disordered fatty acid β-oxidation mediated by vitamin B5 deficiency, could be a putative BVS mechanism. Furthermore, intestinal dysbacteriosis due to microbiome variations in pond water was considered a potential reason for vitamin B5 deficiency. Our BVS index, based on biomarkers identified from the A. japonicus gut microbiota, was a useful tool for BVS diagnosis. Finally, vitamin B5 supplementation was successfully used to treat BVS, suggesting an association with BVS etiology.

IMPORTANCE Body vesicular syndrome (BVS) is a novel disease in sea cucumber aquaculture. As no phenotypic features are visible, BVS is difficult to confirm during aquaculture and postharvest activities, until animals are boiled. Therefore, BVS could lead to severe economic losses compared with other diseases in sea cucumber aquaculture. In this study, for the first time, we systematically investigated BVS pathogenesis and proposed an effective treatment for the condition. Moreover, based on the gut microbiota, we established a noninvasive diagnostic method for BVS in sea cucumber.

RNA-Seq reveals divergent gene expression between larvae with contrasting trophic modes in the poecilogonous polychaete Boccardia wellingtonensis

The polychaete Boccardia wellingtonensis is a poecilogonous species that produces different larval types. Females may lay Type I capsules, in which only planktotrophic larvae are present, or Type III capsules that contain planktotrophic and adelphophagic larvae as well as nurse eggs. While planktotrophic larvae do not feed during encapsulation, adelphophagic larvae develop by feeding on nurse eggs and on other larvae inside the capsules and hatch at the juvenile stage. Previous works have not found differences in the morphology between the two larval types; thus, the factors explaining contrasting feeding abilities in larvae of this species are still unknown. In this paper, we use a transcriptomic approach to study the cellular and genetic mechanisms underlying the different larval trophic modes of B. wellingtonensis. By using approximately 624 million high-quality reads, we assemble the de novo transcriptome with 133,314 contigs, coding 32,390 putative proteins. We identify 5221 genes that are up-regulated in larval stages compared to their expression in adult individuals. The genetic expression profile differed between larval trophic modes, with genes involved in lipid metabolism and chaetogenesis over expressed in planktotrophic larvae. In contrast, up-regulated genes in adelphophagic larvae were associated with DNA replication and mRNA synthesis.

Comparison of microbial communities and the antibiotic resistome between prawn mono- and poly-culture systems

Antibiotic resistance genes (ARGs) in mariculture sediments pose a potential risk to public health due to their ability to transfer from environmental bacteria to human pathogens. Long term, this may reduce pathogen susceptibility to antibiotics in medical settings. In recent years, the poly-culture of multiple species has become a popular mariculture approach in China, thanks to its environmental and economic benefits. However, differences in microbial communities and antibiotic resistome between mono- and poly-culture systems are still unclear. In this study, microbial community composition and profiles of entire (microbial DNA) and mobile (plasmid and phage) ARGs in prawn mono- and poly-culture systems were investigated using metagenomics. The abundance of several viruses and human pathogens were enhanced in prawn poly-culture ponds, when compared to monoculture systems. In contrast, sediments from poly-culture systems had a lower diversity and ARG abundance when compared to mono-culture approaches. These ARG variations were predominantly related to mobile genetic elements. Prawn mariculture activities exerted a unique selectivity for ARGs in plasmids, and this selectivity was not influenced by culture methods. The findings of this study have important implications for the selection of mariculture systems in preventing pollution with ARGs.

Epidemiology of a SKin Ulceration Disease (SKUD) in the sea cucumber Holothuria scabra with a review on the SKUDs in Holothuroidea (Echinodermata)

Aquacultivated sea cucumbers often suffer from SKin Ulceration Diseases (SKUDs). SKUDs have been observed in six holothuroid species from nine countries. All SKUDs present a similar symptom-the skin ulceration-and can be induced by bacteria, viruses, or abiotic factors. We here provide an update on SKUDs in holothuroids and analyse the case of the SKUD observed in Holothuria scabra in Madagascar. Field observations revealed a seasonality of the disease (i.e. wintertime maximum peak). Morphological analyses of integument ulcers showed that sea cucumbers react by forming a collagen fibre plug. Metagenomic analyses revealed a higher proportion of Vibrionaceae (Gammaproteobacteria) in ulcers in comparison to the healthy integument of the same individuals. Experimental infection assays were performed with ulcer crude extracts and bacteria isolated from these extracts (e.g. Vibrio parahaemolyticus) but did not significantly induce skin ulceration. Our results suggest that the disease is not induced by a pathogen or, at the very least, that the pathogen is not found within the ulcers as the disease is not transmissible by contact. An initial cause of the SKUD in Madagascar might be the repeated and prolonged exposures to cold temperatures. Opportunistic bacteria could settle in the dermis of ulcerated individuals and promote the ulcer extension. We propose a general nomenclature for SKUDs based on the acronym of the disease, the affected sea cucumber species (e.g. Hs for Holothuria scabra), the concerned region using an ISO code 3166-2 (e.g. MG for Madagascar), the description date (e.g. 20 for the year 2020), and, when known, the inducing agent (first letter of the general taxon, b for bacteria, v for virus in currently known cases; a a if it is an abiotic inducing parameter; nothing if the inducing cause has not been precisely identified). The disease described in this work will be designated under the name SKUD Hs-MG-20.

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.

A real-time isothermal amplification based portable microfluidic system for simple and reliable detection of Vibrio splendidus

The spread of infectious diseases among aquaculture species has a serious impact on the aquaculture industry. Simple, specific and low-cost detection methods are urgently needed for early diagnosis and timely treatment, particularly for on-site identifying and tracking of pathogens. Vibrio splendidus (V. splendidus) is regarded as one of the main pathogenic bacteria causing skin ulcerative syndrome in cultured sea cucumbers, leading to massive mortality and severe economic losses. We herein present a microfluidic-based real-time fluorogenic loop-mediated isothermal amplification (LAMP) system for simple and reliable detection of V. splendidus. A LAMP primer set with six primers (arsB1) specifically targeting the arsB gene of V. splendidus was successfully designed and tested on the portable microfluidic system for the first time. Only a single step of sample loading using a pipette is required to fill an array of reaction wells (with 10 or 18 wells) in a disposable chip for multiplex detection. A dedicated plastic shell is then utilized to tightly seal the openings of the chip by buckling to prevent contamination and evaporation. Up to four chips (one sample per chip) can be held in the stand-alone and inexpensive microdevice simultaneously, enabling on-demand detection of multiple samples in a single run. Reproducible (relatively low intra- and inter-chip variability) and sensitive (as few as ∼20 CFU, Colony-Forming Units, per reaction well) on-chip arsB1-LAMP assay was demonstrated by using diluted lysate of V. splendidus. A linear standard curve (R2 > 0.98) was attained over the template concentration range of 5 × 103 to 5 × 106 CFU mL-1. V. splendidus can be detected in samples containing different bacteria, indicating the feasibility of the portable microfluidic LAMP system for parallel detection of multiple bacterial pathogens. The proposed on-chip LAMP assay is simple to operate, reliable for amplification, flexible in detection and cost-effective in instrumentation and testing, holding great potential for on-site rapid detection and routine monitoring of aquaculture pathogens.

Insights into high-pressure acclimation: comparative transcriptome analysis of sea cucumber Apostichopus japonicus at different hydrostatic pressure exposures

Background

Global climate change is predicted to force the bathymetric migrations of shallow-water marine invertebrates. Hydrostatic pressure is proposed to be one of the major environmental factors limiting the vertical distribution of extant marine invertebrates. However, the high-pressure acclimation mechanisms are not yet fully understood.

Results

In this study, the shallow-water sea cucumber Apostichopus japonicus was incubated at 15 and 25 MPa at 15 °C for 24 h, and subjected to comparative transcriptome analysis. Nine samples were sequenced and assembled into 553,507 unigenes with a N50 length of 1204 bp. Three groups of differentially expressed genes (DEGs) were identified according to their gene expression patterns, including 38 linearly related DEGs whose expression patterns were linearly correlated with hydrostatic pressure, 244 pressure-sensitive DEGs which were up-regulated at both 15 and 25 MPa, and 257 high-pressure-induced DEGs which were up-regulated at 25 MPa but not up-regulated at 15 MPa.

Conclusions

Our results indicated that the genes and biological processes involving high-pressure acclimation are similar to those related to deep-sea adaptation. In addition to representative biological processes involving deep-sea adaptation (such as antioxidation, immune response, genetic information processing, and DNA repair), two biological processes, namely, ubiquitination and endocytosis, which can collaborate with each other and regulate the elimination of misfolded proteins, also responded to high-pressure exposure in our study. The up-regulation of these two processes suggested that high hydrostatic pressure would lead to the increase of misfolded protein synthesis, and this may result in the death of shallow-water sea cucumber under high-pressure exposure.

Nonspecific expression of fertilization genes in the crown-of-thorns Acanthaster cf. solaris: Unexpected evidence of hermaphroditism in a coral reef predator

The characterization of gene expression in gametes has advanced our understanding of the molecular basis for ecological variation in reproductive success and the evolution of reproductive isolation. These advances are especially significant for ecologically important keystone predators such as the coral-eating crown-of-thorns sea stars (COTS, Acanthaster) which are the most influential predator species in Indo-Pacific coral reef ecosystems and the focus of intensive management efforts. We used RNA-seq and transcriptome assemblies to characterize the expression of genes in mature COTS gonads. We described the sequence and domain organization of eight genes with sex-specific expression and well known functions in fertilization in other echinoderms. We found unexpected expression of genes in one ovary transcriptome that are characteristic of males and sperm, including genes that encode the sperm-specific guanylate cyclase receptor for an egg pheromone, and the sperm acrosomal protein bindin. In a reassembly of previously published RNA-seq data from COTS testes, we found a complementary pattern: strong expression of four genes that are otherwise well known to encode egg-specific fertilization proteins, including the egg receptor for bindin (EBR1) and the acrosome reaction-inducing substance in the egg coat (ARIS1, ARIS2, ARIS3). We also found histological evidence of both eggs and sperm developing in the same gonad in several COTS individuals from a parallel study. These results suggest the occurrence of hermaphrodites, and the potential for reproductive assurance via self-fertilization. Our findings have implications for management of COTS populations, especially in consideration of the large size and massive fecundity of these sea stars.

Transcriptome sequencing reveals phagocytosis as the main immune response in the pathogen-challenged sea urchin Strongylocentrotus intermedius

The clarification of host immune responses to causative bacteria of spotting disease in the sea urchin Strongylocentrotus intermedius is vital to preventing and controlling this disease, especially to selective breeding for disease resistance. For this purpose, sea urchins were challenged with the causative bacterium Vibrio sp. to obtain spotting diseased and undiseased samples. We conducted next-generation sequencing to assess the key genes/pathways in control (CG), diseased (DG), and undiseased (UG) groups. A total of 454.1 million clean reads were obtained and assembled into 23,899 UniGenes with an N50 of 1359 bp, with 86.11% of them matching the genome sequence of the sea urchin S. purpuratus. A total of 8415 UniGenes were mapped to the non-redundant database. Salmon expression analysis revealed 725 significantly differentially expressed genes (DEGs) among CG, DG, and UG. These DEGs were enriched into 72 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including a core set of immune correlated pathways notably in the phagosome, vitamin digestion and absorption, Wnt signaling, and Notch signaling pathways. DG was evidently upregulated in these immune pathways and could enhance phagocytosis directly or indirectly. Thus, phagocytosis was the main coelomic cellular immune response in S. intermedius challenged by spotting disease causative bacterium. The expression patterns of 10 DEGs were confirmed via RT-qPCR, and the expression levels were consistent with the results of RNA-seq. Furthermore, 9899 SSRs were identified, and 123,692, 151,827, and 114,368 candidate SNPs were identified from CG, DG, and UG, respectively. These results provide basic information for our understanding of sea urchin antibacterial immunity.

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.

RNA Sequencing Analysis to Capture the Transcriptome Landscape during Tenderization in Sea Cucumber Apostichopus japonicus

Sea cucumber (Apostichopus japonicus) is an economically significant species in China having great commercial value. It is challenging to maintain the textural properties during thermal processing due to the distinctive physiochemical structure of the A. japonicus body wall (AJBW). In this study, the gene expression profiles associated with tenderization in AJBW were determined at 0 h (CON), 1 h (T_1h), and 3 h (T_3h) after treatment at 37 °C using Illumina HiSeq™ 4000 platform. Seven-hundred-and-twenty-one and 806 differentially expressed genes (DEGs) were identified in comparisons of T_1h vs. CON and T_3h vs. CON, respectively. Among these DEGs, we found that two endogenous proteases-72 kDa type IV collagenase and matrix metalloproteinase 16 precursor-were significantly upregulated that could directly affect the tenderness of AJBW. In addition, 92 genes controlled four types of physiological and biochemical processes such as oxidative stress response (3), immune system process (55), apoptosis (4), and reorganization of the cytoskeleton and extracellular matrix (30). Further, the RT-qPCR results confirmed the accuracy of RNA-sequencing analysis. Our results showed the dynamic changes in global gene expression during tenderization and provided a series of candidate genes that contributed to tenderization in AJBW. This can help further studies on the genetics/molecular mechanisms associated with tenderization.

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.

Comparative Transcriptome Analysis Reveals Growth-Related Genes in Juvenile Chinese Sea Cucumber, Russian Sea Cucumber, and Their Hybrids

Heterosis is important for sea cucumber breeding, but its molecular mechanism remains largely unexplored. In this study, parental lines of Apostichopus japonicus from Russia (R) and China (C) were used to construct hybrids (CR and RC) by reciprocal crossing. We examined the transcriptional profiles of the hybrids (CR and RC) and the purebreds (CC and RR) at different developmental times. A total of 60.27 Gb of clean data was obtained, and 176,649 unigenes were identified, of which 50,312 unigenes were annotated. A total of 414,536 SNPs were identified. A total of 7011 differentially expressed genes (DEGs) were obtained between the purebreds and hybrids at 45 days after fertilization (DAF), and a total of 8218 DEGs were obtained between the purebreds and hybrids at 75 DAF. In addition, a total of 7652 DEGs were obtained between 45 DAF and 75 DAF. The significant DEGs were mainly involved in the MAPK and FOXO signaling pathways, especially in the Ras-Raf-MEK1/2-ERK module, which may be a key regulator of development and growth in juvenile A. japonicus. In addition, we also identified key growth-related genes, such as fgfs, igfs, megfs and hgfs, which were upregulated in the hybrids (RC and CR); these genes may play important roles in heterosis in A. japonicus. Our study provides fundamental information on the molecular mechanisms underlying heterosis in sea cucumber and might suggest strategies for the selection of rapidly growing strains of sea cucumber in aquaculture.