Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

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.

Cloning and functional characterization of the oxidative squalene cyclase gene in the deep-sea holothurian Chiridota sp

Saponins derived from holothurians have high potential medicinal value. However, the de novo synthesis of the derivatization of triterpenes is still unclear. Oxidative squalene cyclase (OSC) can catalyze 2,3-Oxidosqualene into diverse products that serve as important precursors for triterpene synthesis. However, the function of theOSCgene in Chiridotasp. hasnot been elucidated. In this study, an OSCgenederived from the deep-sea holothurianChiridota sp. was cloned and characterized functionally in a yeast system. The open reading frame of the OSC gene was 2086 bp, which encoded 695 amino acids. The Chiridota sp. OSC gene has a similarity of 66.89 % to the OSC of other holothurian species and 63.51 % to that of Acanthaster planci. The phylogenetic tree showed that the echinozoan OSCsclustered together, and then they formeda sister group to fungi and plant homologs. Chiridota sp. OSC catalyzed 2,3-Oxidosqualene into parkeol.Under high pressure, the relative enzymatic activity and stability of cyclase inChiridota sp. was higher than that in the shallow-sea holothurianStichopus horrens. The newly cloned OSC of Chiridota sp.provideskey information for the interpretation of the saponin synthesis pathway in deep-sea holothurians.

Circadian rhythm regulation in the sea cucumber Apostichopus japonicus: Insights into clock gene expression, photoperiod susceptibility, and neurohormone signaling

Sea cucumber Apostichopus japonicus displays the typical circadian rhythms. This present study investigated the molecular regulation of clock genes, as well as monoamines and melatonin, in multiple tissues of A. japonicus, responding to the photoperiod. In order to determine their pivotal role in circadian rhythms, the crucial clock genes, namely AjClock, AjArnt1, AjCry1, and AjTimeless, were identified and a comprehensive analysis of their expressions across various tissues in adult A. japonicus was conducted, revealing the potential existence of central and peripheral oscillators. Results demonstrated that the tissues of polian vesicle and nerve ring exhibited significant clock gene expression associated with the orchestration of circadian regulation, and that environmental light fluctuations exerted influence on the expression of these clock genes. However, a number of genes, such as AjArnt1 and AjCry1, maintained their circadian rhythmicity even under continuous light conditions. Moreover, we further investigated the circadian patterns of melatonin (MT), serotonin (5-HT), and dopamine (DA) secretion in A. japonicus, data that underscored the tissue-specific regulatory differences and the inherent adaptability to dynamic light environments. Collectively, these findings will provide the molecular mechanisms controlling the circadian rhythm in echinoderms and the candidate tissues playing the role of central oscillators in sea cucumbers.

Sulfated Triterpene Glycosides from the Far Eastern Sea Cucumber Cucumaria djakonovi: Djakonoviosides C1, D1, E1, and F1; Cytotoxicity against Human Breast Cancer Cell Lines; Quantitative Structure-Activity Relationships

Four new mono- and trisulfated triterpene penta- and tetraosides, djakonoviosides C1 (1), D1 (2), E1 (3), and F1 (4) were isolated from the Far Eastern sea cucumber Cucumaria djakonovi (Cucumariidae, Dendrochirotida), along with six known glycosides found earlier in other Cucumaria species. The structures of unreported compounds were established on the basis of extensive analysis of 1D and 2D NMR spectra as well as by HR-ESI-MS data. The set of compounds contains six different types of carbohydrate chains including two new ones. Thus, djakonovioside C1 (1) is characterized by xylose as the second residue, that was a branchpoint in the pentasaccharide chain. Meanwhile, only quinovose and rarely glucose have been found earlier in pentasaccharide chains branched at C-2 of the second sugar unit. Djakonovioside E1 (3) is characterized by a tetrasaccharide trisulfated chain, with glucose as the second residue. So, in the series of isolated glycosides, three types of sugars in the second position were presented: the most common, quinovose-in six compounds; glucose-in three substances; and the rare xylose-in one glycoside. The set of aglycones was composed of holostane- and non-holostane-type polycyclic systems; the latter comprised normal and reduced side chains. Noticeably, isokoreoside A (9), isolated from C. djakonovi, was a single glycoside having a 9(11)-double bond, indicating two oxidosqualenecyclases are operating in the process of the biosynthesis of aglycones. Some of the glycosides from C. djakonovi, which were characterized by pentasaccharide branched chains containing one to three sulfate groups, are chemotaxonomic features of the representatives of the genus Cucumaria. The assortment of sugar parts of Cucumaria's glycosides was broadened with previously undescribed penta- and tetrasaccharide moieties. The metabolic network of sugar parts and aglycones is constructed based on biogenetic relationships. The cytotoxic action of compounds 1-10, isolated from C. djakonovi, against human breast cancer cell lines was investigated along with the hemolytic activity. Erythrocytes were, as usual, more sensitive to the membranolytic action of the glycosides than cancer cells. The triple-negative breast cancer MDA-MB-231 cell line was more vulnerable to the action of glycosides in comparison with the other tested cancer cells, while the MCF-7 cell line was less susceptible to cytotoxic action. Djakonovioside E1 (3) demonstrated selective action against ER-positive MCF-7 and triple-negative MDA-MB-231 cell lines, while the toxic effect in relation to normal mammary epithelial cells (MCF-10A) was absent. Cucumarioside A2-5 (6) inhibited the formation and growth of colonies of cancer cells to 44% and tumor cell migration to 85% of the control. Quantitative structure-activity relationships (QSAR) were calculated on the basis of the correlational analysis of the physicochemical properties and structural features of the glycosidic molecules and their membranolytic activity. QSAR revealed the extremely complex nature of such relationships, but these calculations correlated well with the observed SAR.

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.

Chromosome-level genome assembly of the sea cucumber Apostichopus japonicus

Sea cucumber is a morphologically diverse and ecologically important clade of echinoderms. The sea cucumber Apostichopus japonicus is the most economically valuable species of sea cucumber. The initial assembly of the A. japonicus genome was released in 2017. However, this genome assembly is fragmented and lacks relative position information of genes on chromosomes. In this study, we produced a high-quality chromosome-level genome of A. japonicus using Pacbio HiFi long-reads and Hi-C sequencing data. The assembled A. japonicus genome spanned 671.60 Mb with a contig N50 size of 17.20 Mb and scaffold N50 size of 29.65 Mb. A total of 99.9% of the assembly was anchored to 23 chromosomes. In total, 19,828 genes were annotated, and 97.2% of BUSCO genes were fully represented. This high-quality genome of A. japonicus will not only aid in the development of sustainable aquaculture practices, but also lay a foundation for a deeper understanding of their genetic makeup, evolutionary history, and ecological adaptation.

Adverse effects of polystyrene nanoplastics on sea cucumber Apostichopus japonicus and their association with gut microbiota dysbiosis

The mariculture environment is a sink of microplastics (MPs) due to its enclosed nature and mass use of plastics. Nanoplastics (NPs) are MPs with a diameter <1 μm that have a more toxic effect on aquatic organisms than other MPs. However, little is known about the underlying mechanisms of NP toxicity on mariculture species. Here, we performed a multi-omics investigation to explore gut microbiota dysbiosis and associated health problems induced by NPs in juvenile sea cucumber Apostichopus japonicus, a commercially and ecologically important marine invertebrate. We observed significant differences in gut microbiota composition after 21 days of NP exposure. Ingestion of NPs significantly increased core gut microbes, especially Rhodobacteraceae and Flavobacteriaceae families. Additionally, gut gene expression profiles were altered by NPs, especially those related to neurological diseases and movement disorders. Correlation and network analyses indicated close relationships between transcriptome changes and gut microbiota variation. Furthermore, NPs induced oxidative stress in sea cucumber intestines, which may be associated with intraspecies variation in Rhodobacteraceae in the gut microbiota. The results suggested that NPs were harmful to the health of sea cucumbers, and they highlighted the importance of the gut microbiota in the responses to NP toxicity in marine invertebrates.

Discovering a uniform functional trade-off of the CBC-type 2,3-oxidosqualene cyclases and deciphering its chemical logic

Many functionally promiscuous plant 2,3-oxidosqualene cyclases (OSCs) have been found, but complete functional reshaping is rarely reported. In this study, we have identified two new plant OSCs: a unique protostadienol synthase (AoPDS) and a common cycloartenol synthase (AoCAS) from Alisma orientale (Sam.) Juzep. Multiscale simulations and mutagenesis experiments revealed that threonine-727 is an essential residue responsible for protosta-13 (17),24-dienol biosynthesis in AoPDS and that the F726T mutant completely reshapes the native function of AoCAS into a PDS function to yield almost exclusively protosta-13 (17),24-dienol. Unexpectedly, various native functions were uniformly reshaped into a PDS function by introducing the phenylalanine → threonine substitution at this conserved position in other plant and non-plant chair-boat-chair-type OSCs. Further computational modeling elaborated the trade-off mechanisms of the phenylalanine → threonine substitution that leads to the PDS activity. This study demonstrates a general strategy for functional reshaping by using a plastic residue based on the decipherment of the catalytic mechanism.

Bisphenol A (BPA) Directly Activates the G Protein-Coupled Estrogen Receptor 1 and Triggers the Metabolic Disruption in the Gonadal Tissue of Apostichopus japonicus

The sea cucumber, Apostichopus japonicus, is a marine benthic organism that feeds on small benthic particulate matter and is easily affected by pollutants. Bisphenol A (BPA, 4,4'-isopropylidenediphenol) has been identified as an endocrine disruptor. It is ubiquitously detectable in oceans and affects a variety of marine animals. It functions as an estrogen analog and typically causes reproductive toxicity by interfering with the endocrine system. To comparatively analyze the reproductive effects of estradiol (E₂) and BPA on sea cucumbers, we identified a G protein-coupled estrogen receptor 1 (GPER1) in A. japonicus and investigated its effects on reproduction. The results showed that BPA and E₂ exposure activated A. japonicus AjGPER1, thereby mediating the mitogen-activated protein kinase signaling pathways. High-level expression of AjGPER1 in the ovarian tissue was confirmed by qPCR. Furthermore, metabolic changes were induced by 100 nM (22.83 μg/L) BPA exposure in the ovarian tissue, leading to a notable increase in the activities of trehalase and phosphofructokinase. Overall, our findings suggest that AjGPER1 is directly activated by BPA and affects sea cucumber reproduction by disrupting ovarian tissue metabolism, suggesting that marine pollutants pose a threat to the conservation of sea cucumber resources.

The Holothuria leucospilota genome elucidates sacrificial organ expulsion and bioadhesive trap enriched with amyloid-patterned proteins

Some tropical sea cucumbers of the family Holothuriidae can efficiently repel or even fatally ensnare predators by sacrificially ejecting a bioadhesive matrix termed the Cuvierian organ (CO), so named by the French zoologist Georges Cuvier who first described it in 1831. Still, the precise mechanisms for how adhesiveness genetically arose in CO and how sea cucumbers perceive and transduce danger signals for CO expulsion during defense have remained unclear. Here, we report the first high-quality, chromosome-level genome assembly of Holothuria leucospilota, an ecologically significant sea cucumber with prototypical CO. The H. leucospilota genome reveals characteristic long-repeat signatures in CO-specific outer-layer proteins, analogous to fibrous proteins of disparate species origins, including spider spidroin and silkworm fibroin. Intriguingly, several CO-specific proteins occur with amyloid-like patterns featuring extensive intramolecular cross-β structures readily stainable by amyloid indicator dyes. Distinct proteins within the CO connective tissue and outer surface cooperate to give the expelled matrix its apparent tenacity and adhesiveness, respectively. Genomic evidence offers further hints that H. leucospilota directly transduces predator-induced mechanical pressure onto the CO surface through mediation by transient receptor potential channels, which culminates in acetylcholine-triggered CO expulsion in part or in entirety. Evolutionarily, innovative events in two distinct regions of the H. leucospilota genome have apparently spurred CO's differentiation from the respiratory tree to a lethal defensive organ against predators.

Amantadine Toxicity in Apostichopus japonicus Revealed by Proteomics

Amantadine exposure can alter biological processes in sea cucumbers, which are an economically important seafood in China. In this study, amantadine toxicity in Apostichopus japonicus was analyzed by oxidative stress and histopathological methods. Quantitative tandem mass tag labeling was used to examine changes in protein contents and metabolic pathways in A. japonicus intestinal tissues after exposure to 100 µg/L amantadine for 96 h. Catalase activity significantly increased from days 1 to 3 of exposure, but it decreased on day 4. Superoxide dismutase and glutathione activities were inhibited throughout the exposure period. Malondialdehyde contents increased on days 1 and 4 but decreased on days 2 and 3. Proteomics analysis revealed 111 differentially expressed proteins in the intestines of A. japonicus after amantadine exposure compared with the control group. An analysis of the involved metabolic pathways showed that the glycolytic and glycogenic pathways may have increased energy production and conversion in A. japonicus after amantadine exposure. The NF-κB, TNF, and IL-17 pathways were likely induced by amantadine exposure, thereby activating NF-κB and triggering intestinal inflammation and apoptosis. Amino acid metabolism analysis showed that the leucine and isoleucine degradation pathways and the phenylalanine metabolic pathway inhibited protein synthesis and growth in A. japonicus. This study investigated the regulatory response mechanisms in A. japonicus intestinal tissues after exposure to amantadine, providing a theoretical basis for further research on amantadine toxicity.

Gut Microbiota and Metabolites May Play a Crucial Role in Sea Cucumber Apostichopus japonicus Aestivation

The constant increase in temperatures under global warming has led to a prolonged aestivation period for Apostichopus japonicus, resulting in considerable losses in production and economic benefits. However, the specific mechanism of aestivation has not been fully elucidated. In this study, we first tried to illustrate the biological mechanisms of aestivation from the perspective of the gut microbiota and metabolites. Significant differences were found in the gut microbiota of aestivating adult A. japonicus (AAJSD group) compared with nonaestivating adult A. japonicus (AAJRT group) and young A. japonicus (YAJRT and YAJSD groups) based on 16S rRNA gene high-throughput sequencing analysis. The abundances of Desulfobacterota, Myxococcota, Bdellovibrionota, and Firmicutes (4 phyla) in the AAJSD group significantly increased. Moreover, the levels of Pseudoalteromonas, Fusibacter, Labilibacter, Litorilituus, Flammeovirga, Polaribacter, Ferrimonas, PB19, and Blfdi19 genera were significantly higher in the AAJSD group than in the other three groups. Further analysis of the LDA effect size showed that species with significant variation in abundance in the AAJSD group, including the phylum Firmicutes and the genera Litorilituus, Fusibacter, and Abilibacter, might be important biomarkers for aestivating adult A. japonicus. In addition, the results of metabolomics analysis showed that there were three distinct metabolic pathways, namely biosynthesis of secondary metabolites, tryptophan metabolism, and sesquiterpenoid and triterpenoid biosynthesis in the AAJSD group compared with the other three groups. Notably, 5-hydroxytryptophan was significantly upregulated in the AAJSD group in the tryptophan metabolism pathway. Moreover, the genera Labilibacter, Litorilituus, Ferrimonas, Flammeovirga, Blfdi19, Fusibacter, Pseudoalteromonas, and PB19 with high abundance in the gut of aestivating adult A. japonicus were positively correlated with the metabolite 5-HTP. These findings suggest that there may be potential biological associations among the gut microbiota, metabolites, and aestivation in A. japonicus. This work may provide a new perspective for further understanding the aestivation mechanism of A. japonicus.

Identification and functional analysis of Dmrt1 gene and the SoxE gene in the sexual development of sea cucumber, Apostichopus japonicus

Members of the Doublesex and Mab-3-related transcription factor (Dmrt) gene family handle various vital functions in several biological processes, including sex determination/differentiation and gonad development. Dmrt1 and Sox9 (SoxE in invertebrates) exhibit a very conserved interaction function during testis formation in vertebrates. However, the dynamic expression pattern and functional roles of the Dmrt gene family and SoxE have not yet been identified in any echinoderm species. Herein, five members of the Dmrt gene family (Dmrt1, 2, 3a, 3b and 5) and the ancestor SoxE gene were identified from the genome of Apostichopus japonicus. Expression studies of Dmrt family genes and SoxE in different tissues of adult males and females revealed different expression patterns of each gene. Transcription of Dmrt2, Dmrt3a and Dmrt3b was higher expressed in the tube feet and coelomocytes instead of in gonadal tissues. The expression of Dmrt1 was found to be sustained throughout spermatogenesis. Knocking-down of Dmrt1 by means of RNA interference (RNAi) led to the downregulation of SoxE and upregulation of the ovarian regulator foxl2 in the testes. This indicates that Dmrt1 may be a positive regulator of SoxE and may play a role in the development of the testes in the sea cucumber. The expression level of SoxE was higher in the ovaries than in the testes, and knocking down of SoxE by RNAi reduced SoxE and Dmrt1 expression but conversely increased the expression of foxl2 in the testes. In summary, this study indicates that Dmrt1 and SoxE are indispensable for testicular differentiation, and SoxE might play a functional role during ovary differentiation in the sea cucumber.

Cadmium nitrate and DNA methylation in gastropods: comparison between ovotestis and hepatopancreas

Dietary ingestion is the main route of exposure to hazardous contaminants in land animals. Cadmium, a high-profile toxic metal, affects living systems at different organismal levels, including major storage organs (liver, kidneys), key organs for species survival (gonads), and epigenetic networks regulating gene expression. 5-methylcytosine (5mC) is the most common and best-characterized epigenetic mark among different modified nucleosides in DNA. This important player in methylation-driven gene expression is impacted by cadmium in sentinel terrestrial vertebrates. However, limited information exists regarding its impact on macroinvertebrates, especially land snails commonly used as (eco)toxicological models. We first investigate the methylomic effects of dietary cadmium given as cadmium nitrate on terrestrial mollusks. Mature specimens of the common brown garden snail, Cornu aspersum, were continuously exposed for four weeks to environmentally-relevant cadmium levels. We determined global genomic DNA methylation in hepatopancreas and ovotestis, as well as changes in the methylation status of CG pairs at the 5' region close to the transcription site of gene encoding the Cd-selective metallothionein (Cd-MT). Weight gain/loss, hypometabolism tendency, and survival rates were also assessed. Although this exposure event did not adversely affect survival, gastropods exposed to the highest Cd dose revealed a significant reduction in body weight and a significant increase in hypometabolic behavior. The hepatopancreas, but not the ovotestis, displayed significant hypermethylation, but only for the aforementioned specimens. We also found that the 5' end of the Cd-MT gene was unmethylated in both organs and its methylation status was insensitive to cadmium exposure. Our results are important since they provide scientists, for the first time, with quantitative data on DNA methylation in gastropod ovotestis and refine our understanding of Cd epigenetic effects on terrestrial mollusks.

The Isolation, Structure Elucidation and Bioactivity Study of Chilensosides A, A1, B, C, and D, Holostane Triterpene Di-, Tri- and Tetrasulfated Pentaosides from the Sea Cucumber Paracaudina chilensis (Caudinidae, Molpadida)

Five new triterpene (4,4,14-trimethylsterol) di-, tri- and tetrasulfated pentaosides, chilensosides A (1), A₁ (2), B (3), C (4), and D (5) were isolated from the Far-Eastern sea cucumber Paracaudina chilensis. The structures were established on the basis of extensive analysis of 1D and 2D NMR spectra and confirmed by HR-ESI-MS data. The structural variability of the glycosides concerned the pentasaccharide chains. Their architecture was characterized by the upper semi-chain consisting of three sugar units and the bottom semi-chain of two sugars. Carbohydrate chains of compounds 2–5 differed in the quantity and positions of sulfate groups. The interesting structural features of the glycosides were: the presence of two sulfate groups at C-4 and C-6 of the same glucose residue in the upper semi-chain of 1, 2, 4, and 5 and the sulfation at C-3 of terminal glucose residue in the bottom semi-chain of 4 that makes its further elongation impossible. Chilensoside D (5) was the sixth tetrasulfated glycoside found in sea cucumbers. The architecture of the sugar chains of chilensosides A–D (1–5), the positions of sulfation, the quantity of sulfate groups, as well as the aglycone structures, demonstrate their similarity to the glycosides of the representatives of the order Dendrochirotida, confirming the phylogenetic closeness of the orders Molpadida and Dendrochirotida. The cytotoxic activities of the compounds 1–5 against human erythrocytes and some cancer cell lines are presented. Disulfated chilensosides A₁ (2) and B (3) and trisulfated chilensoside C (4) showed significant cytotoxic activity against human cancer cells.

Cloning and characterization of a phosphomevalonate kinase gene that is involved in saponin biosynthesis in the sea cucumber Apostichopus japonicus

The sea cucumber Apostichopus japonicus is one of the most dominant and economically important aquaculture species in China. Saponin, which possesses notable biological and pharmacological properties, is a key determinant of the nutritional and health value of A. japonicus. In the present study, we amplified the full-length cDNA of a phosphomevalonate kinase (PMK) gene (named AjPMK) using rapid amplification of cDNA ends (RACE). Subsequently, we engineered a recombinant AjPMK (rAjPMK) protein and assessed its enzymatic activity by enzyme-linked immunosorbent assay (ELISA). Proteins that interact with rAjPMK were screened and identified via pull-down assay combined with liquid chromatography with tandem mass spectrometry (LC-MS/MS). We found that the full-length cDNA of AjPMK contained 1354 bp and an open reading frame (ORF) of 612 bp. The AjPMK protein was predicted not to contain a signal peptide but to contain a phosphonolate kinase domain seen in higher eukaryotes and a P-loop with a relatively conserved nucleoside triphosphate hydrolase domain. The molecular weight of the AjPMK protein was estimated to be 23.81 kDa, and its isoelectric point was predicted to be 8.72. Phylogenetic analysis showed that AjPMK had a closer evolutionary relationship with genes from starfish than with those of other selected species. Besides, we found that rAjPMK synthesized mevalonate-5-diphosphate, interacted either directly or indirectly with crucial pattern recognition receptors (PRRs) and was regulated by immune-related processes, including antioxidative reactions, stress resistance responses and enzyme hydrolysis. Moreover, AjPMK also interacted with farnesyl pyrophosphate synthase, an enzyme reported to be involved in saponin biosynthesis. Together, our findings implied that AjPMK may be directly involved in saponin biosynthesis and the regulation of various innate immune processes.

Characterization of oxidosqualene cyclases from Trichosanthes cucumerina L. reveals key amino acids responsible for substrate specificity of isomultiflorenol synthase

Two key amino acids of isomultiflorenol synthase, Y125 and M254, were first proposed. They could be associated with the production of isomultiflorenol. Oxidosqualene cyclases (OSCs) are the first committed enzymes in the triterpenoid biosynthesis by converting 2,3-oxidosqualene to specific triterpenoid backbones. Thus, these enzymes are potential targets for developing plant-active compounds through the study of triterpenoid biosynthesis. We applied transcriptome information and metabolite profiling from Trichosanthes cucumerina L. to define the diversity of triterpenoids in this plant through OSCs. Isomultiflorenol synthase and cucurbitadienol synthase were previously identified in this plant. Here, three new OSCs, TcBAS, TcLAS, and TcCAS, were cloned and functionally characterized as β-amyrin synthase, lanosterol synthase, and cycloartenol synthase activities, respectively. We also took advantage of the multiple sequence alignment and molecular docking of OSCs exhibiting in this plant and other plant OSCs to identify key residues associated with isomultiflorenol synthase specificity. Two novel key amino acids, referred to the Y125 and M254, were first discovered. These results provide information on a possible catalytic mechanism for plant OSCs that produce specific products.

Development of a high-throughput SNP array for sea cucumber (Apostichopus japonicus) and its application in genomic selection with MCP regularized deep neural networks

High-throughput single nucleotide polymorphism (SNP) genotyping assays are powerful tools for genetic studies and genomic breeding applications for many species. Though large numbers of SNPs have been identified in sea cucumber (Apostichopus japonicus), but, as yet, no high-throughput genotyping platform is available for this species. In this study, we designed and developed a high-throughput 24 K SNP genotyping array named HaishenSNP24K for A. japonicus, based on the multi-objective-local optimization (MOLO) algorithm and HD-Marker genotyping method. The SNP array exhibited a relatively high genotyping call rate (> 96%), genotyping accuracy (>95%) and exhibited highly polymorphic in sea cucumber populations. In addition, we also assessed its application in genomic selection (GS). Deep neural networks (DNN) that can capture the complicated interactions of genes have been proposed as a promising tool in GS for SNP-based genomic prediction of complex traits in animal breeding. To overcome the problem of over-fitting when using the HaishenSNP24K array as high-dimensional DNN input, we developed minmax concave penalty (MCP) regularization for sparse deep neural networks (DNN-MCP) that finds an optimal sparse structure of a DNN by minimizing the square error subject to the non-convex penalty MCP on the parameters (weights and biases). Compared to two linear models, namely RR-GBLUP and Bayes B, and the nonlinear model DNN, DNN-MCP has greatly improved the genomic prediction ability for three quantitative traits (e.g., wet weight, dry weight and survival time) in the sea cucumber population. To the best of our knowledge, this is the first work to develop a high-throughput SNP array for A. japonicus and a new model DNN-MCP for genomic prediction of complex traits in GS. The present results provide evidence that supports the HaishenSNP24K array with DNN-MCP will be valuable for genetic studies and molecular breeding in A. japonicus.

Biosynthesis of saponin defensive compounds in sea cucumbers

Soft-bodied slow-moving sea creatures such as sea stars and sea cucumbers lack an adaptive immune system and have instead evolved the ability to make specialized protective chemicals (glycosylated steroids and triterpenes) as part of their innate immune system. This raises the intriguing question of how these biosynthetic pathways have evolved. Sea star saponins are steroidal, while those of the sea cucumber are triterpenoid. Sterol biosynthesis in animals involves cyclization of 2,3-oxidosqualene to lanosterol by the oxidosqualene cyclase (OSC) enzyme lanosterol synthase (LSS). Here we show that sea cucumbers lack LSS and instead have two divergent OSCs that produce triterpene saponins and that are likely to have evolved from an ancestral LSS by gene duplication and neofunctionalization. We further show that sea cucumbers make alternate sterols that confer protection against self-poisoning by their own saponins. Collectively, these events have enabled sea cucumbers to evolve the ability to produce saponins and saponin-resistant sterols concomitantly.

Sea stars and sea cucumbers biosynthesize protective glycosylated steroids and triterpenes via divergent oxidosqualene cyclases (OSCs) that produce these distinct saponins in different species as well as in different tissues of a single species.

Structures and Biologic Activity of Chitonoidosides I, J, K, K1 and L-Triterpene Di-, Tri- and Tetrasulfated Hexaosides from the Sea Cucumber Psolus chitonoides

Five new triterpene di-, tri- and tetrasulfated hexaosides (chitonoidosides I (1), J (2), K (3), K₁ (4) and L (5)) were isolated from the Far-Eastern sea cucumber Psolus chitonoides, collected near Bering Island (Commander Islands) from a depth of 100–150 m. The structural variability of the glycosides concerned both the aglycones (with 7(8)- or 9(11)-double bonds) and carbohydrate chains differing from each other by the third sugar residue (Xyl or sulfated by C-6 Glc) and/or by the fourth—terminal in the bottom semi-chain—residue (Glc or sulfated by C-6 MeGlc) as well as by the positions of a sulfate group at C-4 or C-6 in the sixth—terminal in the upper semi-chain—residue (MeGlc). Hemolytic activities of these compounds 1–5 against human erythrocytes as well as cytotoxicity against human cancer cell lines, HeLa, DLD-1 and HL-60, were studied. The hexaosides, chitonoidosides K (3) and L (5) with four sulfate groups, were the most active against tumor cells in all the tests. Noticeably, the sulfate group at C-4 of MeGlc6 did not decrease the membranolytic effect of 5 as compared with 3, having the sulfate group at C-6 of MeGlc6. Erythrocytes were, as usual, more sensitive to the action of the studied glycosides than cancer cells, although the sensitivity of leukemia promyeloblast HL-60 cells was higher than that of other tumor cells. The glycosides 1 and 2 demonstrated some weaker action in relation to DLD-1 cells than against other tumor cell lines. Chitonoidoside K₁ (4) with a hydroxyl at C 25 of the aglycone was not active in all the tests. The metabolic network formed by the carbohydrate chains of all the glycosides isolated from P. chitonoides as well as the aglycones biosynthetic transformations during their biosynthesis are discussed and illustrated with schemes.

Identification of sex determination locus in sea cucumber Apostichopus japonicus using genome-wide association study

Background

Sex determination mechanisms are complicated and diverse across taxonomic categories. Sea cucumber Apostichopus japonicus is a benthic echinoderm, which is the closest group of invertebrates to chordate, and important economic and ecologically aquaculture species in China. A. japonicus is dioecious, and no phenotypic differences between males and females can be detected before sexual maturation. Identification of sex determination locus will broaden knowledge about sex-determination mechanism in echinoderms, which allows for the identification of sex-linked markers and increases the efficiency of sea cucumber breeding industry.

Results

Here, we integrated assembly of a novel chromosome-level genome and resequencing of female and male populations to investigate the sex determination mechanisms of A. japonicus. We built a chromosome-level genome assembly AJH1.0 using Hi-C technology. The assembly AJH1.0 consists of 23 chromosomes ranging from 22.4 to 60.4 Mb. To identify the sex-determination locus of A. japonicus, we conducted genome-wide association study (GWAS) and analyses of distribution characteristics of sex-specific SNPs and fixation index F_ST. The GWAS analysis showed that multiple sex-associated loci were located on several chromosomes, including chromosome 4 (24.8%), followed by chromosome 9 (10.7%), chromosome 17 (10.4%), and chromosome 18 (14.1%). Furthermore, analyzing the homozygous and heterozygous genotypes of plenty of sex-specific SNPs in females and males confirmed that A. japonicus might have a XX/XY sex determination system. As a physical region of 10 Mb on chromosome 4 included the highest number of sex-specific SNPs and higher F_ST values, this region was considered as the candidate sex determination region (SDR) in A. japonicus.

Conclusions

In the present study, we integrated genome-wide association study and analyses of sex-specific variations to investigate sex determination mechanisms. This will bring novel insights into gene regulation during primitive gonadogenesis and differentiation and identification of master sex determination gene in sea cucumber. In the sea cucumber industry, investigation of molecular mechanisms of sex determination will be helpful for artificial fertilization and precise breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08632-3.

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.

Exploiting the gut microbiota to predict the origins and quality traits of cultured sea cucumbers

Nowadays, the true economic and nutritional value of food is underpinned by both origin and quality traits, more often expressed as increased quality benefits derived from the origin source. Gut microbiota contribute to food metabolism and host health, therefore, it may be suitable as a qualifying indicator of origin and quality of economic species. Here, we investigated relationships between the gut microbiota of the sea cucumber (Apostichopus japonicus), a valuable aquaculture species in Asia, with their origins and quality metrics. Based on data from 287 intestinal samples, we generated the first biogeographical patterns for A. japonicus gut microbiota from origins across China. Importantly, A. japonicus origins were predicted using the random forest model that was constructed using 20 key gut bacterial genera, with 97.6% accuracy. Furthermore, quality traits such as saponin, fat and taurine were also successfully predicted by random forest models based on gut microbiota, with approximately 80% consistency between predicted and true values. We showed that substantial variations existed in the gut microbiota and quality variables in A. japonicus across different origins, and we also demonstrated the great potential of gut microbiota to track A. japonicus origins and predict their quality traits.

The genome of an apodid holothuroid (Chiridota heheva) provides insights into its adaptation to a deep-sea reducing environment

Cold seeps and hydrothermal vents are deep-sea reducing environments that are characterized by lacking oxygen and photosynthesis-derived nutrients. Most animals acquire nutrition in cold seeps or hydrothermal vents by maintaining epi- or endosymbiotic relationship with chemoautotrophic microorganisms. Although several seep- and vent-dwelling animals hosting symbiotic microbes have been well-studied, the genomic basis of adaptation to deep-sea reducing environment in nonsymbiotic animals is still lacking. Here, we report a high-quality genome of Chiridota heheva Pawson & Vance, 2004, which thrives by extracting organic components from sediment detritus and suspended material, as a reference for nonsymbiotic animal's adaptation to deep-sea reducing environments. The expansion of the aerolysin-like protein family in C. heheva compared with other echinoderms might be involved in the disintegration of microbes during digestion. Moreover, several hypoxia-related genes (Pyruvate Kinase M2, PKM2; Phospholysine Phosphohistidine Inorganic Pyrophosphate Phosphatase, LHPP; Poly(A)-specific Ribonuclease Subunit PAN2, PAN2; and Ribosomal RNA Processing 9, RRP9) were subject to positive selection in the genome of C. heheva, which contributes to their adaptation to hypoxic environments.

Xenophagy of invasive bacteria is differentially activated and modulated via a TLR-TRAF6-Beclin1 axis in echinoderms

In marine environments, organisms are confronted with numerous microbial challenges, although the differential regulation of xenophagy in response to different pathogenic bacterial species remains relatively unknown. Here, we addressed this issue using Apostichopus japonicus as a model. We identified 39 conserved autophagy-related genes by genome-wide screening, which provided a molecular basis for autophagy regulation in sea cucumbers. Furthermore, xenophagy of two Gram-negative bacteria, Vibrio splendidus and Escherichia coli, but not a Gram-positive bacteria, Micrococcus luteus, was observed in different autophagy assays. Surprisingly, a significantly higher autophagy capacity was found in the E. coli–challenged group than in the V. splendidus–challenged group. To confirm these findings, two different lipopolysaccharides, LPS^{V. splendidus} and LPS^E. coli, were isolated; we found that these LPS species differentially activated coelomocyte xenophagy. To explore the molecular mechanism mediating differential levels of xenophagy, we used an siRNA knockdown assay and confirmed that LPS^{V. splendidus}-mediated xenophagy was dependent on an AjTLR3-mediated pathway, whereas LPS^E. coli-mediated xenophagy was dependent on AjToll. Moreover, the activation of different AjTLRs resulted in AjTRAF6 ubiquitination and subsequent activation of K63-linked ubiquitination of AjBeclin1. Inversely, the LPS^{V. splendidus}-induced AjTLR3 pathway simultaneously activated the expression of AjA20, which reduced the extent of K63-linked ubiquitination of AjBeclin1 and impaired the induction of autophagy; however, this finding was no t evident with LPS^E. coli. Our present results provide the first evidence showing that xenophagy could be differentially induced by different bacterial species to yield differential autophagy levels in echinoderms.

Proteomics reveals gender differences in physiological characteristics of the gonads and tube feet of the sea cucumber, Apostichopus japonicus

The sea cucumber Apostichopus japonicus is an important aquaculture species in China because of its high nutritional and medicinal values. Gender, as a factor affecting the physiology of organisms, is always considered when improving the breeding efficiency of economically important animals. In the present study, protein expression profiles of the gonads and tube feet of male and female A. japonicus were investigated using a comparative proteomics approach. A set of 7499 proteins were identified, which covered a broad range of functions based on function annotations. A significant difference in protein expression profiles was observed between the gonads and tube feet of A. japonicus; gonads showed more apparent gender differences than tube feet. Moreover, the findings revealed that male A. japonicus had more specific functions and most of these functions were associated with energy consumption. Further analyses suggested that the regulation of ERK activity and the capacity of tyrosine production and virus immunity might be more powerful in male and female A. japonicus, respectively. Some candidate proteins were also recognized as potential targets for gender identification of A. japonicus. Overall, our study provides new insights into the understanding of molecular mechanisms underlying gender-based physiological differences in A. japonicus. SIGNIFICANCE: The current study aimed to reveal gender differences in the physiological characteristics of gonads and tube feet of the sea cucumber A. japonicus. To the best of our knowledge, this is the first proteomics study to analyze the differences in the protein expression profiles of external organs between male and female A. japonicus. The analysis revealed gender differences in the protein expression profiles of both gonads and tube feet of A. japonicus, and the gender differences in gonads were quite apparent. Moreover, according to the recognition of differentially expressed proteins and the enrichment analyses based on Kyoto Encyclopedia of Genes and Genomes, a draft view of how the physiological functions of A. japonicus were affected by gender was obtained. Male A. japonicus could have more specific functions related to energy consumption than females. The regulation of ERK activity and virus immunity might be more robust in male and female A. japonicus, respectively. Some candidate proteins were also recognized as potential targets for gender identification of A. japonicus. The findings presented here will improve the understanding of researchers about the molecular mechanisms underlying gender-based differences in A. japonicus and contribute to the meticulous breeding of A. japonicus.

Effects of Cadmium Sulfate on the Brown Garden Snail Cornu aspersum: Implications for DNA Methylation

An extensive literature exists regarding the cellular, physiological, and genetic effects of cadmium (Cd)—A highly toxic, but commonly used trace metal in modern industry. However, limited data are available on its epigenetic effects, especially for terrestrial sentinel invertebrates. We determined Cd retention, total DNA methylation, and the methylation status of 5′ end of the Cd-MT gene in the hepatopancreas of the brown garden snail, Cornu aspersum, fed Cd sulfate for four weeks. Bodyweight changes and survival were also measured. Hepatopancreas cadmium increased in a dose-dependent manner from the third-lowest dose onward, with very large amounts being found for the highest treatment group. However, no mortalities occurred, irrespective of dietary Cd dose. We identified significant genome-wide hypermethylation in specimens given the highest dose, which overlapped with a significant bodyweight decrease. The Cd-MT gene showed an unmethylated 5′ end of the Cd-MT gene and this status was not affected by cadmium exposure. Hepatopancreas DNA methylation is as sensitive as bodyweight to non-lethal concentrations of dietary Cd given as cadmium sulfate but less responsive than tissue accumulation. Such an exposure event, by contrast, does not affect the methylation status of the Cd-MT gene 5′ end.

tRNA copy number and codon usage in the sea cucumber genome provide insights into adaptive translation for saponin biosynthesis

Genomic tRNA copy numbers determine cytoplasmic tRNA abundances, which in turn influence translation efficiency, but the underlying mechanism is not well understood. Using the sea cucumber Apostichopus japonicus as a model, we combined genomic sequence, transcriptome expression and ecological food resource data to study its codon usage adaptation. The results showed that, unlike intragenic non-coding RNAs, transfer RNAs (tRNAs) tended to be transcribed independently. This may be attributed to their specific Pol III promoters that lack transcriptional regulation, which may underlie the correlation between genomic copy number and cytoplasmic abundance of tRNAs. Moreover, codon usage optimization was mostly restrained by a gene's amino acid sequence, which might be a compromise between functionality and translation efficiency for stress responses were highly optimized for most echinoderms, while enzymes for saponin biosynthesis (LAS, CYPs and UGTs) were especially optimized in sea cucumbers, which might promote saponin synthesis as a defence strategy. The genomic tRNA content of A. japonicus was positively correlated with amino acid content in its natural food particles, which should promote its efficiency in protein synthesis. We propose that coevolution between genomic tRNA content and codon usage of sea cucumbers facilitates their saponin synthesis and survival using food resources with low nutrient content.

Progress in the Studies of Triterpene Glycosides From Sea Cucumbers (Holothuroidea, Echinodermata) Between 2017 and 2021

Structural diversity of triterpene glycosides produced by sea cucumbers or holothurians (Holothuroidea, Echinodermata) is extremely high, although all of them are either lanostane derivatives or, rarely, products of their molecular rearrangements. The majority of them are holostane derivatives possessing an 18(20)-lanostane lactone as aglycone. They contain carbohydrate chains consisting of one to six monosaccharide units including sulfated ones. The glycosides demonstrate interesting biological activities, mainly caused by membranolytic action, namely cytotoxic, ichthyotoxic, antifungal, and hemolytic properties, as well as a series of additional effects at sub-toxic doses, including immunomodulatory, and cancer preventive. This review summarizes the literature data concerning structures and biological activities of all the new triterpene glycosides isolated from sea cucumbers during 2017 to 2021.

A pan-metazoan concept for adult stem cells: the wobbling Penrose landscape

Adult stem cells (ASCs) in vertebrates and model invertebrates (e.g. Drosophila melanogaster) are typically long‐lived, lineage‐restricted, clonogenic and quiescent cells with somatic descendants and tissue/organ‐restricted activities. Such ASCs are mostly rare, morphologically undifferentiated, and undergo asymmetric cell division. Characterized by ‘stemness’ gene expression, they can regulate tissue/organ homeostasis, repair and regeneration. By contrast, analysis of other animal phyla shows that ASCs emerge at different life stages, present both differentiated and undifferentiated phenotypes, and may possess amoeboid movement. Usually pluri/totipotent, they may express germ‐cell markers, but often lack germ‐line sequestering, and typically do not reside in discrete niches. ASCs may constitute up to 40% of animal cells, and participate in a range of biological phenomena, from whole‐body regeneration, dormancy, and agametic asexual reproduction, to indeterminate growth. They are considered legitimate units of selection. Conceptualizing this divergence, we present an alternative stemness metaphor to the Waddington landscape: the ‘wobbling Penrose’ landscape. Here, totipotent ASCs adopt ascending/descending courses of an ‘Escherian stairwell’, in a lifelong totipotency pathway. ASCs may also travel along lower stemness echelons to reach fully differentiated states. However, from any starting state, cells can change their stemness status, underscoring their dynamic cellular potencies. Thus, vertebrate ASCs may reflect just one metazoan ASC archetype.

Diverse triterpene skeletons are derived from the expansion and divergent evolution of 2,3-oxidosqualene cyclases in plants

Triterpenoids are one of the largest groups of secondary metabolites and exhibit diverse structures, which are derived from C30 skeletons that are biosynthesized via the isoprenoid pathway by cyclization of 2,3-oxidosqualene. Triterpenoids have a wide range of biological activities, and are used in functional foods, drugs, and as industrial materials. Due to the low content levels in their native plants and limited feasibility and efficiency of chemical synthesis, heterologous biosynthesis of triterpenoids is the most promising strategy. Herein, we classified 121 triterpene alcohols/ketones according to their conformation and ring numbers, among which 51 skeletons have been experimentally characterized as the products of oxidosqualene cyclases (OSCs). Interestingly, 24 skeletons that have not been reported from nature source were generated by OSCs in heterologous expression. Comprehensive evolutionary analysis of the identified 152 OSCs from 75 species in 25 plant orders show that several pentacyclic triterpene synthases repeatedly originated in multiple plant lineages. Comparative analysis of OSC catalytic reaction revealed that stabilization of intermediate cations, steric hindrance, and conformation of active center amino acid residues are primary factors affecting triterpene formation. Optimization of OSC could be achieved by changing of side-chain orientations of key residues. Recently, methods, such as rationally design of pathways, regulation of metabolic flow, compartmentalization engineering, etc., were introduced in improving chassis for the biosynthesis of triterpenoids. We expect that extensive study of natural variation of large number of OSCs and catalytical mechanism will provide basis for production of high level of triterpenoids by application of synthetic biology strategies.

Clustering genomic organization of sea cucumber miRNAs impacts their evolution and expression

Echinoderms are marine deuterostomes with fascinating adaptation features such as aestivation and organ regeneration. However, post-transcriptional gene regulation by microRNAs (miRNAs) underlying these features are largely unexplored. Here, using homology-based and de novo approaches supported by expression data, we provided a comprehensive annotation of miRNA genes in the sea cucumber Apostichopus japonicus. By linkage and phylogenic analyses, we characterized miRNA genomic organization, evolutionary history and expression regulation. The results showed that sea cucumbers evolved a large number of new miRNAs, which tended to form polycistronic clusters via tandem duplication that had been especially active in the echinoderms. Most new miRNAs were weakly expressed, but miRNA clustering increased the expression level of clustered new miRNAs. The most abundantly expressed new miRNAs were organized in a single tandem cluster (cluster n2), which was activated during aestivation and intestine regeneration. Overall, our analyses suggest that clustering of miRNAs is important for their evolutionary origin, expression control, and functional cooperation.

Comparative transcriptome analysis reveals changes in gene expression in sea cucumber (Holothuria leucospilota) in response to acute temperature stress

Ambient temperature is an important abiotic factor that influences growth performance and physiological functions in sea cucumbers. To understand the molecular responses of the sea cucumber Holothuria leucospilota to acute temperature stress, we performed a de novo transcriptome analysis of body wall tissue from H. leucospilota exposed to 2 hoursh of acute heat (35 ± 1 °C) and cold stress (15 ± 1 °C). A total of 99,015 unigenes were obtained after assembly of the sequenced reads. Compared with a control group maintained at 25.0 ± 1 °C, 1169 differentially expressed unigenes (DEGs) were identified after heat stress, 781 were up-regulated and 388 were down-regulated. After cold stress, 1464 DEGs were identified; 900 were up-regulated and 564 were down-regulated. The annotation of DEGs revealed that heat shock proteins play important roles in protecting H. leucospilota from high temperature stress. Furthermore, KEGG pathway enrichment analysis showed that the categories: "Ribosome" and "Protein processing in endoplasmic reticulum" were strongly affected by heat stress. These two pathways are associated with biosynthesis and processing of proteins, and refolding of misfolded proteins. The lipid metabolism pathways "Sphingolipid metabolism" and "Ether lipid metabolism", were affected by cold stress. The RNA-Seq results for eight selected DEGs were verified the expression by quantitative real-time PCR analysis. Our results will improve the understanding of the molecular response mechanisms of H. leucospilota to ambient temperature stress.

Triterpene Glycosides from the Far Eastern Sea Cucumber Thyonidium (=Duasmodactyla) kurilensis (Levin): The Structures, Cytotoxicities, and Biogenesis of Kurilosides A3, D1, G, H, I, I1, J, K, and K1

Nine new mono-, di-, and trisulfated triterpene penta- and hexaosides, kurilosides A₃ (1), D₁ (2), G (3), H (4), I (5), I₁ (6), J (7), K (8), and K₁ (9) and two desulfated derivatives, DS-kuriloside L (10), having a trisaccharide branched chain, and DS-kuriloside M (11), having hexa-nor-lanostane aglycone with a 7(8)-double bond, have been isolated from the Far-Eastern deep-water sea cucumber Thyonidium (=Duasmodactyla) kurilensis (Levin) and their structures were elucidated based on 2D NMR spectroscopy and HR-ESI mass-spectrometry. Five earlier unknown carbohydrate chains and two aglycones (having a 16β,(20S)-dihydroxy-fragment and a 16β-acetoxy,(20S)-hydroxy fragment) were found in these glycosides. All the glycosides 1–9 have a sulfate group at C-6 Glc, attached to C-4 Xyl1, while the positions of the other sulfate groups vary in different groups of kurilosides. The analysis of the structural features of the aglycones and the carbohydrate chains of all the glycosides of T. kurilensis showed their biogenetic relationships. Cytotoxic activities of the compounds 1–9 against mouse neuroblastoma Neuro 2a, normal epithelial JB-6 cells, and erythrocytes were studied. The highest cytotoxicity in the series was demonstrated by trisulfated hexaoside kuriloside H (4), having acetoxy-groups at C(16) and C(20), the latter one obviously compensated the absence of a side chain, essential for the membranolytic action of the glycosides. Kuriloside I₁ (6), differing from 4 in the lacking of a terminal glucose residue in the bottom semi-chain, was slightly less active. The compounds 1–3, 5, and 8 did not demonstrate cytotoxic activity due to the presence of hydroxyl groups in their aglycones.

Nutritional Aspects of Ecologically Relevant Phytochemicals in Ruminant Production

This review provides an update of ecologically relevant phytochemicals for ruminant production, focusing on their contribution to advancing nutrition. Phytochemicals embody a broad spectrum of chemical components that influence resource competence and biological advantage in determining plant species' distribution and density in different ecosystems. These natural compounds also often act as plant defensive chemicals against predatorial microbes, insects, and herbivores. They may modulate or exacerbate microbial transactions in the gastrointestinal tract and physiological responses in ruminant microbiomes. To harness their production-enhancing characteristics, phytochemicals have been actively researched as feed additives to manipulate ruminal fermentation and establish other phytochemoprophylactic (prevent animal diseases) and phytochemotherapeutic (treat animal diseases) roles. However, phytochemical-host interactions, the exact mechanism of action, and their effects require more profound elucidation to provide definitive recommendations for ruminant production. The majority of phytochemicals of nutritional and pharmacological interest are typically classified as flavonoids (9%), terpenoids (55%), and alkaloids (36%). Within flavonoids, polyphenolics (e.g., hydrolyzable and condensed tannins) have many benefits to ruminants, including reducing methane (CH4) emission, gastrointestinal nematode parasitism, and ruminal proteolysis. Within terpenoids, saponins and essential oils also mitigate CH4 emission, but triterpenoid saponins have rich biochemical structures with many clinical benefits in humans. The anti-methanogenic property in ruminants is variable because of the simultaneous targeting of several physiological pathways. This may explain saponin-containing forages' relative safety for long-term use and describe associated molecular interactions on all ruminant metabolism phases. Alkaloids are N-containing compounds with vast pharmacological properties currently used to treat humans, but their phytochemical usage as feed additives in ruminants has yet to be exploited as they may act as ghost compounds alongside other phytochemicals of known importance. We discussed strategic recommendations for phytochemicals to support sustainable ruminant production, such as replacements for antibiotics and anthelmintics. Topics that merit further examination are discussed and include the role of fresh forages vis-à-vis processed feeds in confined ruminant operations. Applications and benefits of phytochemicals to humankind are yet to be fully understood or utilized. Scientific explorations have provided promising results, pending thorough vetting before primetime use, such that academic and commercial interests in the technology are fully adopted.

Marine natural products

This review covers the literature published between January and December in 2018 for marine natural products (MNPs), with 717 citations (706 for the period January to December 2018) referring to compounds isolated from marine microorganisms and phytoplankton, green, brown and red algae, sponges, cnidarians, bryozoans, molluscs, tunicates, echinoderms, mangroves and other intertidal plants and microorganisms. The emphasis is on new compounds (1554 in 469 papers for 2018), together with the relevant biological activities, source organisms and country of origin. Reviews, biosynthetic studies, first syntheses, and syntheses that led to the revision of structures or stereochemistries, have been included. The proportion of MNPs assigned absolute configuration over the last decade is also surveyed.

Comparative metabolomic analysis of the body wall from four varieties of the sea cucumber Apostichopus japonicus

The sea cucumber, Apostichopus japonicas, is an important economic species with high nutritive value. In recent years, driven by significant market demand, the sea cucumber breeding industry has developed rapidly. Body color and number of papillae are important factors that determine the value and price of sea cucumbers. In this study, metabolite profiling of four sea cucumber varieties (green, white, purple and spiny) was performed using ultra-performance liquid chromatography-quadrupole-time-of-flight-mass spectrometry (UPLC-QTOF-MS) combined with multivariate analysis. Orthogonal partial least-squares discriminant analysis (OPLS-DA) clearly discriminated the body wall metabolites of four sea cucumber varieties. Differential metabolites included fatty acids, phospholipids, and sugars. KEGG metabolic pathway analysis revealed that carbohydrate, protein, and phospholipid metabolism were highly conserved among the varieties. These results provide a comprehensive insight into differences in the metabolite profile of four A. japonicus varieties and a deeper understanding of sea cucumber varieties breeding.

Identification and functional characterization of natural resistance-associated macrophage protein 2 from sea cucumber Apostichopus japonicus

As a member of natural resistance-associated macrophage protein (Nramp) family, Nramp2 conservatively exists in the cell membrane across species and is essential for normal iron homeostasis in an H⁺-dependent manner. Withholding available iron represents an important host defense strategy. However, the function of Nramp2 in response to invading pathogens is largely unknown in invertebrates. In this study, a unique echinoderm Nramp2 was identified from sea cucumber Apostichopus japonicus (designated as AjNramp2). The cDNA sequence of AjNramp2 was 2360 bp, with a putative open reading frame of 1713 bp, encoding a typical Nramp domain containing protein with 570 amino acid residues. Structural analysis revealed that AjNramp2 consisted of highly conserved helix regions similar with the human Nramp2. Spatial expression analysis revealed that AjNramp2 was ubiquitously expressed in all examined tissues, with the highest level found in the intestine. Immunohistochemistry assay showed that AjNramp2 was mainly located in the cellular membrane in coelomocytes. Vibrio splendidus challenge and lipopolysaccharide (LPS) stimulation could significantly promote the expression of AjNramp2, which was consistent with the cellular iron level in coelomocytes. Moreover, when the expression of AjNramp2 was knocked down by siRNA-AjNramp2, the cellular iron level was coordinately decreased in coelomocytes under LPS stimulation. Taken together, results indicated that AjNramp2 serves as an iron transport receptor to withhold available iron and may contribute to the nutritional immunity defense system of sea cucumber.

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.