Differential gene expression in the respiratory tree of the sea cucumber Apostichopus japonicus during aestivation

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.

Alterations in mitochondrial structure and function in response to environmental temperature changes in Apostichopus japonicus

Global temperatures have risen as a result of climate change, and the resulting warmer seawater will exert physiological stresses on many aquatic animals, including Apostichopus japonicus. It has been suggested that the sensitivity of aquatic poikilothermal animals to climate change is closely related to mitochondrial function. Therefore, understanding the interaction between elevated temperature and mitochondrial functioning is key to characterizing organisms' responses to heat stress. However, little is known about the mitochondrial response to heat stress in A. japonicus. In this work, we investigated the morphological and functional changes of A. japonicus mitochondria under three representative temperatures, control temperature (18 °C), aestivation temperature (25 °C) and heat stress temperature (32 °C) temperatures using transmission electron microscopy (TEM) observation of mitochondrial morphology combined with proteomics and metabolomics techniques. The results showed that the mitochondrial morphology of A. japonicus was altered, with decreases in the number of mitochondrial cristae at 25 °C and mitochondrial lysis, fracture, and vacuolization at 32 °C. Proteomic and metabolomic analyses revealed 103 differentially expressed proteins and 161 differential metabolites at 25 °C. At 32 °C, the levels of 214 proteins and 172 metabolites were significantly altered. These proteins and metabolites were involved in the tricarboxylic acid (TCA) cycle, substance transport, membrane potential homeostasis, anti-stress processes, mitochondrial autophagy, and apoptosis. Furthermore, a hypothetical network of proteins and metabolites in A. japonicus mitochondria in response to temperature changes was constructed based on proteomic and metabolomic data. These results suggest that the dynamic regulation of mitochondrial energy metabolism, resistance to oxidative stress, autophagy, apoptosis, and mitochondrial morphology in A. japonicus may play important roles in the response to elevated temperatures. In summary, this study describes the response of A. japonicus mitochondria to temperature changes from the perspectives of morphology, proteins, and metabolites, which provided a better understanding the mechanisms of mitochondrial regulation under environment stress in marine echinoderms.

Selection and validation of genes related to oxidative stress production and clearance in macrophages infected with Mycobacterium tuberculosis

Background

In the fight against tuberculosis, besides chemotherapy, the regulation of oxidative stress (OS) has also aroused people's interest in host-oriented therapy. However, there is limited research on the genes involved in reactive oxygen species (ROS) production and clearance in macrophages infected with Mycobacterium tuberculosis (MTB). This study analyzes and explores this to provide a basis for exploring new targets for antituberculosis treatments.

Methods

We established a macrophage model infected with MTB, counted intracellular bacteria, and determined the ROS produced using flow cytometry. We conducted ribonucleic acid sequencing, screened differentially expressed genes through transcriptomic methods, and validated the expression of them through reverse transcription-quantitative polymerase chain reaction.

Results

The ROS of macrophages increased with intracellular bacteria at 4 h after infection with MTB and reached its peak at 48 h, surpassing the uninfected macrophages (p < 0.05). A total of 1,613 differentially expressed genes were identified after infection with MTB, of which 458 were associated with ROS, with over 50% involved in the response of organelles and biological processes to stimuli. We analyzed and identified six genes. After macrophage infection with MTB, the expression of CAMK2B increased, whereas the expression of CYBB decreased (p < 0.05). The expression of GPX3 and SOD2 increased, whereas the expression of CAT decreased (p < 0.05).

Conclusion

The ROS-related differentially expressed genes between MTB infected and uninfected macrophages may be related to some organelles and involved in various biological processes, molecular functions, and signaling pathways. Among them, CAMK2B, GPX3, and SOD2 may be related to ROS.

Aestivation in Nature: Physiological Strategies and Evolutionary Adaptations in Hypometabolic States

Aestivation is considered to be one of the "purest" hypometabolic states in nature, as it involves aerobic dormancy that can be induced and sustained without complex factors. Animals that undergo aestivation to protect themselves from environmental stressors such as high temperatures, droughts, and food shortages. However, this shift in body metabolism presents new challenges for survival, including oxidative stress upon awakening from aestivation, accumulation of toxic metabolites, changes in energy sources, adjustments to immune status, muscle atrophy due to prolonged immobility, and degeneration of internal organs due to prolonged food deprivation. In this review, we summarize the physiological and metabolic strategies, key regulatory factors, and networks utilized by aestivating animals to address the aforementioned components of aestivation. Furthermore, we present a comprehensive overview of the advancements made in aestivation research across major species, including amphibians, fish, reptiles, annelids, mollusks, and echinoderms, categorized according to their respective evolutionary positions. This approach offers a distinct perspective for comparative analysis, facilitating an understanding of the shared traits and unique features of aestivation across different groups of organisms.

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.

Mechanism underlying the toxicity of the microplastic fibre transfer in the sea cucumber Apostichopus japonicus

Microscopic plastic particles (0.1 µm-5 mm) are widespread hazardous pollutants, and microfibres (MFs) are their dominant shape in habitats. Previous field and laboratory studies have demonstrated that MFs enter the coelomic fluid of sea cucumbers from the water through the respiratory tree. However, the possible mechanism underlying the toxicity of this process is not well understood. Herein, RNA-Seq was used to examine the responses of the respiratory tree during the MF transfer process in the sea cucumber Apostichopus japonicus. Polyester synthetic MFs were used, and the number of transferred MFs was controlled to the amount reported from the field. The results showed that the MFs altered gene expression as the transfer process increased. The top genes regulated by MF transfer were mainly related to metabolic processes and signal transduction pathways, with upregulated genes following low MF transfer and downregulated genes following high MF transfer. Functional enrichment analysis revealed the pathways in which differentially expressed genes were enriched under different MF transfer scenarios. The transcriptomic findings were further supported by histological observations, which revealed injury and loss of cell components. This study contributes to understanding the effects of MFs in a valuable echinoderm species through transcriptomic and histological examinations.

The possible routes of microplastics uptake in sea cucumber Holothuria cinerascens (Brandt, 1835)

Investigating the routes of microplastic uptake in sea cucumber is crucial in this plastics pollution age considering their non-selective feeding process, nutritional, biomedical and ecological importance. The guts, respiratory trees and coelomic fluids of Holothuria cinerascens were sampled and examined for microplastic after exposure to fluorescent microplastic fragments and microfibres. The madreporite pore size was also determined. 90% of the animals sampled ingested microplastic fragments via their feeding tentacles. Microplastic ingested ranged from 0 to 24 fragments intestine-¹ (8.7 ± 2.11). All (100%) the animals sampled had microfibre in their coelomic fluid ranging from 32 to 227 microfibres coelom-¹ (79.58 ± 10.53). Microfibres were found attached to all undigested respiratory trees examined. Microfibres were also found in 57.8% of digested respiratory trees with a range of 0-12 microfibres respiratory tree^-1 (1.74 ± 0.66). Notably, there was no fluorescent microplastic fragments/microfibres found in the gut, coelomic fluid, and respiratory trees of animals in the control group. The madreporite pore size ranged from 0.59 to 2.90 μm (1.22 ± 0.03 μm). Microfibres found in the coelomic fluid were transferred from the respiratory trees because the size of microfibre used for this experiment was assumed larger for it to have passed through the gut into the coelom. Although the madreporite pore size is smaller than microfibre used in this study, madreporite can act as another channel to transport nanoplastics from the coelom into the animal tissue. This study showed that sea cucumber uptake microplastics from the environment using the feeding tentacles and the respiratory trees. Chemicals desorption from microplastics into H. cinerascens when ingested may pose a risk to the animals and their consumers. Further study is required to determine if microfibres transfer from the respiratory trees to the coelomic fluid is dependent on microfibre concentration and water temperature.

De novo assembly transcriptome analysis reveals the preliminary molecular mechanism of pigmentation in juveniles of the hard clam Mercenaria mercenaria

Color plays a vital function in camouflage, sexual selection, immunity, and evolution. Mollusca possess vivid shell colors and pigmentation starts at the juvenile stage. The hard clam Mercenaria mercenaria is a widely cultivated bivalve of high economic value. To explore the molecular mechanism of pigmentation in juvenile clams, here, we performed RNA-Seq analysis on non-pigmented, white, and red M. mercenaria specimens. Clean reads were assembled into 358,285 transcripts and 149,234 unigenes, whose N50 lengths were 2107 bp and 1567 bp, respectively. Differentially expressed genes were identified and analyzed for KEGG enrichment. "Melanoma/Melanogenesis", "ABC transporters", and "Porphyrin and chlorophyll metabolism" pathways appeared to be associated with pigmentation. Pathways related to carotenoid metabolism seemed to also play a vital role in pigmentation in juveniles. Our results provide new insights into the formation of shell color in juvenile hard clams.

Plasticity of Respiratory Function Accommodates High Oxygen Demand in Breeding Sea Cucumbers

Physiological plasticity allows animals to adjust their physiological function to abiotic and biotic variations. It has been mostly studied in the context of response to external factors and not much is known on how animals adjust their physiology to cope with variations in internal conditions. The process of reproduction implies gonadal maturation and other internal changes, bringing various challenges to the animal such as an increased demand for energy and oxygen. Here, the capacity of the sea cucumber, Apostichopus japonicus to adjust its respiratory function and physiological mechanisms during reproduction was studied using a time-lapse videography and metabolomics approach. The results showed that reproduction caused a significant increase in oxygen consumption in A. japonicus. Interestingly, breeding sea cucumbers can accommodate the high oxygen demand by accelerating respiratory rate. However, to maintain a necessary high level of respiratory activity during reproduction, sea cucumbers need consume large amounts of adenosine triphosphate (ATP). In addition, the metabolomic data suggests that oxidative stress and hormone regulation are the physiological mechanisms linking reproduction and respiratory function. Altogether, these findings suggest that plasticity of respiratory function is an effective tactic to cope with high oxygen demand during reproduction.

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.

Estivation-responsive microRNAs in a hypometabolic terrestrial snail

When faced with extreme environmental conditions, the milk snail (Otala lactea) enters a state of dormancy known as estivation. This is characterized by a strong reduction in metabolic rate to <30% of normal resting rate that is facilitated by various behavioural, physiological, and molecular mechanisms. Herein, we investigated the regulation of microRNA in the induction of estivation. Changes in the expression levels of 75 highly conserved microRNAs were analysed in snail foot muscle, of which 26 were significantly upregulated during estivation compared with controls. These estivation-responsive microRNAs were linked to cell functions that are crucial for long-term survival in a hypometabolic state including anti-apoptosis, cell-cycle arrest, and maintenance of muscle functionality. Several of the microRNA responses by snail foot muscle also characterize hypometabolism in other species and support the existence of a conserved suite of miRNA responses that regulate environmental stress responsive metabolic rate depression across phylogeny.

De Novo assembly and comparative transcriptome analyses of purple and green morphs of Apostichopus japonicus during body wall pigmentation process

Pigmentation processes provide a traceable and relevant trait for understanding key issues in evolutionary biology such as adaptation, speciation and the maintenance of balanced polymorphisms. The sea cucumber Apostichopus japonicus, which has nutritive and medical properties, is considered the most valuable commercial species in many parts of Asia. Compared with the green morph, the purple morph is rare and has great appeal to consumers. However, little is currently known about the molecular mechanism of body color formation in A. japonicus, even in echinoderm. Here, we employ illumina sequencing to examine expression patterns of the gene network underlying body wall development in purple and green morphs of A. japonicus. Overall, the number of down-regulated genes in the green morph was significantly more than in the purple morph during the pigmentation stage. We observed dynamic expression patterns of a large number of pigment, regulation and growth genes from the "Melanogenesis", "Melanoma", "Wnt signaling pathway", "Notch signaling pathway", "epithelium development", "epidermal growth factor receptor binding","growth factor activity" and "growth", including contrasting expression patterns of these genes in green and purple morph. This study provides comprehensive lists of differentially expressed genes during body wall development in the green and purple morphs, revealing potential candidate genes that may be involved in regulating body color formation and polymorphism. These data will provide valuable information for future genetic studies on sea cucumbers elucidating the molecular mechanisms underlying pigmentation, and may support the culture of desirable color morphs.

Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation

Echinoderms exhibit several fascinating evolutionary innovations that are rarely seen in the animal kingdom, but how these animals attained such features is not well understood. Here we report the sequencing and analysis of the genome and extensive transcriptomes of the sea cucumber Apostichopus japonicus, a species from a special echinoderm group with extraordinary potential for saponin synthesis, aestivation and organ regeneration. The sea cucumber does not possess a reorganized Hox cluster as previously assumed for all echinoderms, and the spatial expression of Hox7 and Hox11/13b potentially guides the embryo-to-larva axial transformation. Contrary to the typical production of lanosterol in animal cholesterol synthesis, the oxidosqualene cyclase of sea cucumber produces parkeol for saponin synthesis and has "plant-like" motifs suggestive of convergent evolution. The transcriptional factors Klf2 and Egr1 are identified as key regulators of aestivation, probably exerting their effects through a clock gene-controlled process. Intestinal hypometabolism during aestivation is driven by the DNA hypermethylation of various metabolic gene pathways, whereas the transcriptional network of intestine regeneration involves diverse signaling pathways, including Wnt, Hippo and FGF. Decoding the sea cucumber genome provides a new avenue for an in-depth understanding of the extraordinary features of sea cucumbers and other echinoderms.

Differential gene expression in the intestine of sea cucumber (Apostichopus japonicus) under low and high salinity conditions

Sea cucumber, Apostichopus japonicus is an important species for aquaculture, and its behavior and physiology can change in response to changing salinity conditions. For this reason, it is important to understand the molecular responses of A. japonicus when exposed to ambient changes in salinity. In this study, RNA-Seq provided a general overview of the gene expression profiles in the intestine of A. japonicus exposed to high salinity (SD40), normal salinity (SD30) and low salinity (SD20) environments. Screening for differentially expressed genes (DEGs) using the NOISeq method identified 109, 100, and 89 DEGs based on a fold change of ≥2 and divergence probability ≥0.8 according to the comparisons of SD20 vs. SD30, SD20 vs.SD40, and SD30 vs. SD40, respectively. Gene ontology analysis showed that the terms "metabolic process" and "catalytic activity" comprised the most enriched DEGs. These fell into the categories of "biological process" and "molecular function". While "cell" and "cell part" had the most enriched DEGs in the category of "cellular component". With these DEGs mapping to 2119, 159, and 160 pathways in the Kyoto Encyclopedia of Genes and Genomes database. Of these 51, 2, and 57 pathways were significantly enriched, respectively. The osmosis-specific DEGs identified in this study of A. japonicus will be important targets for further studies to understand the biochemical mechanisms involved with the adaption of sea cucumbers to changes in salinity.

Identification and functional characterisation of 5-HT4 receptor in sea cucumber Apostichopus japonicus (Selenka)

Serotonin (5-HT) is an important neurotransmitter and neuromodulator that controls a variety of sensory and motor functions through 5-HT receptors (5-HTRs). The 5-HT4R subfamily is linked to Gs proteins, which activate adenylyl cyclases (ACs), and is involved in many responses in peripheral organs. In this study, the 5-HT4R from Apostichopus japonicus (Aj5-HT4R) was identified and characterised. The cloned full-length Aj5-HT4R cDNA is 1,544 bp long and contains an open reading frame 1,011 bp in length encoding 336 amino acid proteins. Bioinformatics analysis of the Aj5-HT4R protein indicated this receptor was a member of class A G protein coupled receptor (GPCR) family. Further experiments using Aj5-HT4R-transfected HEK293 cells demonstrated that treatment with 5-HT triggered a significant increase in intracellular cAMP level in a dose-dependent manner and induced a rapid internalisation of Aj5-HT4R fused with enhanced green fluorescent protein (Aj5-HT4R-EGFP) from the cell surface into the cytoplasm. In addition, the transcriptional profiles of Aj5-HT4R in aestivating A. japonicas and phosphofructokinase (AjPFK) in 5-HT administrated A. japonicus have been analysed by real-time PCR assays. Results have led to a basic understanding of Aj5-HT4R in A. japonicus, and provide a foundation for further exploration of the cell signaling and regulatory functions of this receptor.

Effect of acute salinity stress on ion homeostasis, Na+/K+-ATPase and histological structure in sea cucumber Apostichopus japonicus

Background

Sea cucumbers (Apostichopus japonicus) are an imperiled fauna exposed to a variety of environmental condition such as salinity and studies are urgently needed to assess their effects to guide aquaculture efforts. The effects of acute salinity stress on coelomic fluid osmotic pressure, ion concentrations, the activity of Na⁺/K⁺-ATPase in respiratory trees and the histological variations were measured to evaluate the salinity tolerance of sea cucumbers.

Results

Significant correlations in osmotic pressure were observed between coelomic fluid and ambient environmental salinity. In coelomic fluid, Na⁺ concentration was observed fluctuated during salinity 18 psu and the inflection point presented at the 6 h. The Na⁺/K⁺-ATPase activity in respiratory trees indicated the “U-shaped” fluctuant change and the change trend was opposite with the Na⁺ concentration. The ions (K⁺, Cl⁻) concentration decreased and showed the same tendency at salinity 40 psu with salinity 18 psu. The total coelomocytes counts and phagocytosis of coelomic fluid Na⁺/K⁺-ATPase activity indicated fluctuating changes under different salinity stress. Histological variation revealed a negative relation between decreasing salt concentration and tissue integrity. Tissue damages were significantly observed in intestines, muscles and tube feet under low salinity environment (18, 23 and 27 psu). The connective tissue in intestines of A. japonicus exposed to 18 and 23 psu damaged and partly separated from the mucosal epithelium. The significant variations occurred in tube feet, which presented the swelling in connective tissue and a fracture in longitudinal muscles under low salinity (18 psu). The morphological change of tube feet showed the shrinkage of connective tissue under high salinity (40 psu). The amount of infusoria in the respiratory trees decreased or even disappeared in salinity treatment groups (18 and 23 psu).

Conclusion

The results inferred that osmoconformity and ionoregulation were seen in sea cucumbers, which contributed to understand the salinity regulatory mechanisms of A. japonicus under acute salinity stress.

Transcriptomic Analysis of Differentially Expressed Genes During Larval Development of Rapana venosa by Digital Gene Expression Profiling

During the life cycle of shellfish, larval development, especially metamorphosis, has a vital influence on the dynamics, distribution, and recruitment of natural populations, as well as seed breeding. Rapana venosa, a carnivorous gastropod, is an important commercial shellfish in China, and is an ecological invader in the United States, Argentina, and France. However, information about the mechanism of its early development is still limited, because research in this area has long suffered from a lack of genomic resources. In this study, 15 digital gene expression (DGE) libraries from five developmental stages of R. venosa were constructed and sequenced on the IIIumina Hi-Sequation 2500 platform. Bioinformaticsanalysis identified numerous differentially and specifically expressed genes, which revealed that genes associated with growth, nervous system, digestive system, immune system, and apoptosis participate in important developmental processes. The functional analysis of differentially expressed genes was further implemented by gene ontology, and Kyoto encyclopedia of genes and genomes enrichment. DGE profiling provided a general picture of the transcriptomic activities during the early development of R. venosa, which may provide interesting hints for further study. Our data represent the first comparative transcriptomic information available for the early development of R. venosa, which is a prerequisite for a better understanding of the physiological traits controlling development.

De Novo Adult Transcriptomes of Two European Brittle Stars: Spotlight on Opsin-Based Photoreception

Next generation sequencing (NGS) technology allows to obtain a deeper and more complete view of transcriptomes. For non-model or emerging model marine organisms, NGS technologies offer a great opportunity for rapid access to genetic information. In this study, paired-end Illumina HiSeqTM technology has been employed to analyse transcriptomes from the arm tissues of two European brittle star species, Amphiura filiformis and Ophiopsila aranea. About 48 million Illumina reads were generated and 136,387 total unigenes were predicted from A. filiformis arm tissues. For O. aranea arm tissues, about 47 million reads were generated and 123,324 total unigenes were obtained. Twenty-four percent of the total unigenes from A. filiformis show significant matches with sequences present in reference online databases, whereas, for O. aranea, this percentage amounts to 23%. In both species, around 50% of the predicted annotated unigenes were significantly similar to transcripts from the purple sea urchin, the closest species to date that has undergone complete genome sequencing and annotation. GO, COG and KEGG analyses were performed on predicted brittle star unigenes. We focused our analyses on the phototransduction actors involved in light perception. Firstly, two new echinoderm opsins were identified in O. aranea: one rhabdomeric opsin (homologous to vertebrate melanopsin) and one RGR opsin. The RGR-opsin is supposed to be involved in retinal regeneration while the r-opsin is suspected to play a role in visual-like behaviour. Secondly, potential phototransduction actors were identified in both transcriptomes using the fly (rhabdomeric) and mammal (ciliary) classical phototransduction pathways as references. Finally, the sensitivity of O.aranea to monochromatic light was investigated to complement data available for A. filiformis. The presence of microlens-like structures at the surface of dorsal arm plate of O. aranea could potentially explain phototactic behaviour differences between the two species. The results confirm (i) the ability of these brittle stars to perceive light using opsin-based photoreception, (ii) suggest the co-occurrence of both rhabdomeric and ciliary photoreceptors, and (iii) emphasise the complexity of light perception in this echinoderm class.

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.

Identification and Characterization of MicroRNAs from Longitudinal Muscle and Respiratory Tree in Sea Cucumber (Apostichopus japonicus) Using High-Throughput Sequencing

MicroRNAs (miRNAs), as a family of non-coding small RNAs, play important roles in the post-transcriptional regulation of gene expression. Sea cucumber (Apostichopus japonicus) is an important economic species which is widely cultured in East Asia. The longitudinal muscle (LTM) and respiratory tree (RPT) are two important tissues in sea cucumber, playing important roles such as respiration and movement. In this study, we identified and characterized miRNAs in the LTM and RPT of sea cucumber (Apostichopus japonicus) using Illumina HiSeq 2000 platform. A total of 314 and 221 conserved miRNAs were identified in LTM and RPT, respectively. In addition, 27 and 34 novel miRNAs were identified in the LTM and RPT, respectively. A set of 58 miRNAs were identified to be differentially expressed between LTM and RPT. Among them, 9 miRNAs (miR-31a-3p, miR-738, miR-1692, let-7a, miR-72a, miR-100b-5p, miR-31b-5p, miR-429-3p, and miR-2008) in RPT and 7 miRNAs (miR-127, miR-340, miR-381, miR-3543, miR-434-5p, miR-136-3p, and miR-300-3p) in LTM were differentially expressed with foldchange value being greater than 10. A total of 14,207 and 12,174 target genes of these miRNAs were predicted, respectively. Functional analysis of these target genes of miRNAs were performed by GO analysis and pathway analysis. This result provided in this work will be useful for understanding biological characteristics of the LTM and RPT of sea cucumber and assisting molecular breeding of sea cucumber for aquaculture.