Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation

Genome-Wide DNA Methylation Signatures of Sea Cucumber Apostichopus japonicus during Environmental Induced Aestivation

Organisms respond to severe environmental changes by entering into hypometabolic states, minimizing their metabolic rates, suspending development and reproduction, and surviving critical ecological changes. They come back to an active lifestyle once the environmental conditions are conducive. Marine invertebrates live in the aquatic environment and adapt to environmental changes in their whole life. Sea cucumbers and sponges are only two recently known types of marine organisms that aestivate in response to temperature change. Sea cucumber has become an excellent model organism for studies of environmentally-induced aestivation by marine invertebrates. DNA methylation, the most widely considered epigenetic marks, has been reported to contribute to phenotypic plasticity in response to environmental stress in aquatic organisms. Most of methylation-related enzymes, including DNA methyltransferases, Methyl-CpG binding domain proteins, and DNA demethylases, were up-regulated during aestivation. We conducted high-resolution whole-genome bisulfite sequencing of the intestine from sea cucumber at non-aestivation and deep-aestivation stages. Further DNA methylation profile analysis was also conducted across the distinct genomic features and entire transcriptional units. A different elevation in methylation level at internal exons was observed with clear demarcation of intron/exon boundaries during transcriptional unit scanning. The lowest methylation level occurs in the first exons, followed by the last exons and the internal exons. A significant increase in non-CpG methylation (CHG and CHH) was observed within the intron and mRNA regions in aestivation groups. A total of 1393 genes were annotated within hypermethylated DMRs (differentially methylated regions), and 749 genes were annotated within hypomethylated DMRs. Differentially methylated genes were enriched in the mRNA surveillance pathway, metabolic pathway, and RNA transport. Then, 24 hypermethylated genes and 15 hypomethylated genes were Retrovirus-related Pol polyprotein from transposon (RPPT) genes. This study provides further understanding of epigenetic control on environmental induced hypometabolism in aquatic organisms.

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.

Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae

The transient larva-bearing biphasic life cycle is the hallmark of many metazoan phyla, but how metazoan larvae originated remains a major enigma in animal evolution. There are two hypotheses for larval origin. The 'larva-first' hypothesis suggests that the first metazoans were similar to extant larvae, with later evolution of the adult-added biphasic life cycle; the 'adult-first' hypothesis suggests that the first metazoans were adult forms, with the biphasic life cycle arising later via larval intercalation. Here, we investigate the evolutionary origin of primary larvae by conducting ontogenetic transcriptome profiling for Mollusca-the largest marine phylum characterized by a trochophore larval stage and highly variable adult forms. We reveal that trochophore larvae exhibit rapid transcriptome evolution with extraordinary incorporation of novel genes (potentially contributing to adult shell evolution), and that cell signalling/communication genes (for example, caveolin and innexin) are probably crucial for larval evolution. Transcriptome age analysis of eight metazoan species reveals the wide presence of young larval transcriptomes in both trochozoans and other major metazoan lineages, therefore arguing against the prevailing larva-first hypothesis. Our findings support an adult-first evolutionary scenario with a single metazoan larval intercalation, and suggest that the first appearance of proto-larva probably occurred after the divergence of direct-developing Ctenophora from a metazoan ancestor.

Gender specific differences of immune competence in the sea cucumber Apostichopus japonicus before and after spawning

The gender differences of immunity have been elucidated in many vertebrates and invertebrates. However, the information of this difference was still not clear in the sea cucumber Apostichopus japonicus, which is one of the most valuable aquaculture species and susceptible to diseases caused by pathogen infection. In the present study, the transcriptome of coelomocytes from female and male A. japonicus before and after spawning was obtained by RNA-sequencing technology. A total of 4,538 and 8,248 differentially expressed genes were identified between female and male A. japonicus before and after spawning, respectively, indicating that the gender differences of gene expression profiles in A. japonicus were more remarkable after spawning. Further KEGG enrichment analyses were conducted for both male and female up-regulated genes before and after spawning. The results revealed that the capacity to kill pathogens in female A. japonicus might be more powerful than that in males no matter before and after spawning; the antioxidant ability in male A. japonicus was probably stronger than that in females after spawning; the complement system in male A. japonicus might be more effective than that in females after spawning; and the apoptosis was likely to be more serious in male A. japonicus before spawning. Moreover, we speculated that the fatty acid composition might be one of the inducements for gender specific immune differences of A. japonicus. Overall, the results of our study illustrated the global gender specific immune differences of A. japonicus and contributed to understanding of the molecular mechanisms underlying sea cucumber immune regulation.

Neuropeptide precursors and neuropeptides in the sea cucumber Apostichopus japonicus: a genomic, transcriptomic and proteomic analysis

The sea cucumber Apostichopus japonicus is a foodstuff with very high economic value in China, Japan and other countries in south-east Asia. It is at the heart of a multibillion-dollar industry and to meet demand for this product, aquaculture methods and facilities have been established. However, there are challenges associated with optimization of reproduction, feeding and growth in non-natural environments. Therefore, we need to learn more about the biology of A. japonicus, including processes such as aestivation, evisceration, regeneration and albinism. One of the major classes of molecules that regulate physiology and behaviour in animals are neuropeptides, and a few bioactive peptides have already been identified in A. japonicus. To facilitate more comprehensive investigations of neuropeptide function in A. japonicus, here we have analysed genomic and transcriptomic sequence data and proteomic data to identify neuropeptide precursors and neuropeptides in this species. We identified 44 transcripts encoding neuropeptide precursors or putative neuropeptide precursors, and in some instances neuropeptides derived from these precursors were confirmed by mass spectrometry. Furthermore, analysis of genomic sequence data enabled identification of the location of neuropeptide precursor genes on genomic scaffolds and linkage groups (chromosomes) and determination of gene structure. Many of the precursors identified contain homologs of neuropeptides that have been identified in other bilaterian animals. Precursors of neuropeptides that have thus far only been identified in echinoderms were identified, including L- and F-type SALMFamides, AN peptides and others. Precursors of several peptides that act as modulators of neuromuscular activity in A. japonicus were also identified. The discovery of a large repertoire of neuropeptide precursors and neuropeptides provides a basis for experimental studies that investigate the physiological roles of neuropeptide signaling systems in A. japonicus. Looking ahead, some of these neuropeptides may have effects that could be harnessed to enable improvements in the aquaculture of this economically important species.

Triterpenoids in Echinoderms: Fundamental Differences in Diversity and Biosynthetic Pathways

Echinoderms form a remarkable phylum of marine invertebrates that present specific chemical signatures unique in the animal kingdom. It is particularly the case for essential triterpenoids that evolved separately in each of the five echinoderm classes. Indeed, while most animals have Δ5-sterols, sea cucumbers (Holothuroidea) and sea stars (Asteroidea) also possess Δ7 and Δ9(11)-sterols, a characteristic not shared with brittle stars (Ophiuroidea), sea urchins (Echinoidea), and crinoids (Crinoidea). These particular Δ7 and Δ9(11) sterols emerged as a self-protection against membranolytic saponins that only sea cucumbers and sea stars produce as a defense mechanism. The diversity of saponins is large; several hundred molecules have been described in the two classes of these saponins (i.e., triterpenoid or steroid saponins). This review aims to highlight the diversity of triterpenoids in echinoderms by focusing on sterols and triterpenoid glycosides, but more importantly to provide an updated view of the biosynthesis of these molecules in echinoderms.

Predicting Growth Traits with Genomic Selection Methods in Zhikong Scallop (Chlamys farreri)

Selective breeding is a common and effective approach for genetic improvement of aquaculture stocks with parental selection as the key factor. Genomic selection (GS) has been proposed as a promising tool to facilitate selective breeding. Here, we evaluated the predictability of four GS methods in Zhikong scallop (Chlamys farreri) through real dataset analyses of four economical traits (e.g., shell length, shell height, shell width, and whole weight). Our analysis revealed that different GS models exhibited variable performance in prediction accuracy depending on genetic and statistical factors, but non-parametric method, including reproducing kernel Hilbert spaces regression (RKHS) and sparse neural networks (SNN), generally outperformed parametric linear method, such as genomic best linear unbiased prediction (GBLUP) and BayesB. Furthermore, we demonstrated that the predictability relied mainly on the heritability regardless of GS methods. The size of training population and marker density also had considerable effects on the predictive performance. In practice, increasing the training population size could better improve the genomic prediction than raising the marker density. This study is the first to apply non-linear model and neural networks for GS in scallop and should be valuable to help develop strategies for aquaculture breeding programs.