Triterpenoid saponin biosynthesis genes and their expression patterns during the development of sea cucumber Apostichopus japonicus*

Yujia YANG , Xiaoni LI , Lina SUN ,**

1 Laboratory for Evolution and Development, Institute of Evolution and Marine Biodiversity, Ocean University of China, Qingdao 266003, China

2 CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences,Qingdao 266071, China

Abstract Sea cucumbers are echinoderm species with a leathery skin and an elongated body. Natural compounds, including saponins, have been previously isolated from sea cucumber. The majority of saponins derived from this organism typically belong to the triterpene glycosides, which exhibit various biological eff ects (e.g., antifungal, cytotoxic, hemolytic, and cytostatic immunomodulatory activities). In the present study, 30 saponin biosynthesis enzymes, namely 2 AACTs, 2 FPSs, 2 HMGSs, 2 OSCs, 1 SS, 2 SEs, and 19 UGTs, were found in the genome of sea cucumber Apostichopus japonicus. Compared with upstream saponin biosynthesis enzymes (i.e., AACTs, FPSs, HMGSs, OSCs, and SSs), the downstream saponin biosynthesis enzymes (SEs and UGTs) shared lower amino acid sequence similarity with the corresponding genes in other echinoderms. It was proposed that the enzyme structures comprised both α-helices and β-sheets. Notably, only a low percentage of β-sheets were present in FPSs, OSCs, and SS. The saponin biosynthesis enzymes showed dynamic expression patterns during f ive critical developmental stages of A. japonicus (fertilized oocytes, blastula, gastrula, doliolaria, and penractula). The present study involved elucidation of putative saponin biosynthesis pathways in sea cucumber and provides a valuable platform for further investigation of saponin biosynthesis in echinoderms.

Keyword: triterpene saponin biosynthesis; sea cucumber; enzyme; expression; development

1 INTRODUCTION

Saponins are a class of chemical compounds,which are abundant in various plant species, including the soapwort plant, maple, and ginseng. They are characterized by a hydrophobic triterpene or sterol backbone and a hydrophilic carbohydrate chain(Mugford and Osbourn, 2012). Saponins are broadly classif ied into steroid glycosides, triterpenoids,spirostanols, and furostanol saponins (Vincken et al.,2007). It is known that saponins and their biosynthetic intermediates exhibit diff erent biological activities and play particularly important roles in the pharmaceutical, cosmetic, and food industries (Moses et al., 2014; Xiao et al., 2019). Saponins are the main bioactive compounds in numerous plants and folk medicines, especially traditional Chinese remedies(Bahrami and Franco, 2015). As bioactive molecules,saponins have been used as vaccine adjuvants (Kensil,1996), dietary supplements, and nutritional food ingredients (Zhao et al., 2018). While at certain concentrations some saponins are toxic to coldblooded organisms, they can enhance nutrient absorption (e.g., oat and spinach) and aid digestion(Forester and Hartmut, 2006). These natural products are understood to be plant derived. Only a few marine organisms, such as echinoderms, marine sponges, and octocorals, have been shown to produce saponins(Stonik, 2001; Francis et al., 2002; Skropeta, 2008;Mitu et al., 2017). Sea cucumberApostichopusjaponicus, also called marine ginseng in China, is a species of echinoderm. It has been demonstrated that sea cucumber can generate numerous agricultural,nutraceutical, pharmaceutical, and cosmeceutical products (Bahrami et al., 2014). Sea cucumbers are characterized by slow motion and the absence of prominent structural defenses (Moghimipour and Handali, 2015). The abundance of saponins in the body wall and viscera was considered as a physical barrier of sea cucumber (Van Dyck et al., 2010).Highly diverse saponins, including 59 triterpene glycosides, have been previously detected in the organism (Caulier et al., 2011). The structural architecture of the saponins is associated with intraindividual variability, which leads to specif ic functions(Caulier et al., 2011). The majority of saponins isolated from sea cucumber are triterpene glycosides,which exhibit various biological eff ects (e.g.,antifungal, cytotoxic, hemolytic, and cytostatic immunomodulatory activities) (Maier et al., 2001). In addition to playing important roles in the defense mechanisms of sea cucumber, it has been speculated that saponins also inhibit oocyte maturation during spawning and control gametogenesis (Bahrami et al.,2014).

Based on the characteristic reactions, the saponin biosynthesis enzymes could be categorized into four types, i.e., oxidosqualene cyclases (OSCs),cytochrome P450 monooxygenases (P450s), UDPglycosyltransferases (UGTs), and other tailoring enzymes (Moses et al., 2014). In recent years, saponin biosynthesis enzymes have been studied in several plant species, including oat (Mugford et al., 2013),Lotusjaponicus(Krokida et al., 2013), tomato, and potato (Adzhubei et al., 2013; Itkin et al., 2013).Identif ication of biosynthetic genes involved in new pathways is essential to synthesize specialized saponin metabolites in echinoderms such as sea cucumber. A recent study found that the OSC1 and OSC2 genes inA.japonicusexhibited high sequence divergence and contained plant-like motifs (Li et al.,2018). Moreover, a transition from lanosterol production to parkeol production was determined,which indicated convergent evolution in sea cucumber(Li et al., 2018). Another transcriptome study revealed the secondary metabolite biosynthesis pathways,including saponin biosynthesis, in the body wall of sea cucumberHolothuriascabra(Mitu et al., 2017).Core saponin biosynthesis genes (e.g.,hydroxymethylglutaryl-CoA synthase, farnesyl diphosphate synthase, squalene synthase, and squalene epoxidase) in the body wall ofH.scabrashared homology with those of plants, suggesting that the saponin biosynthesis pathways were highly conserved (Mitu et al., 2017). Nevertheless, even though the saponin biosynthesis genes were previously identif ied in another sea cucumber species, i.e.,H.scabra, the complete set of saponin biosynthesis enzymes and their expression patterns during critical developmental stages have never been investigated in the genome ofA.japonicus. In this study, we comprehensively studied the biosynthesis enzymes responsible for the generation of diverse saponin metabolites, including acetyl-CoA acetyltransferase(AACT), farnesyl diphosphate synthase (FPS),hydroxymethylglutaryl-CoA synthase (HMGS),oxidosqualene cyclase (OSC), squalene synthase(SS), squalene epoxidase (SE), and UDPglycosyltransferase (UGT). Additionally, their putative mechanisms inA.japonicuswere investigated. The f indings of the current study will facilitate further research into the saponin biosynthesis in echinoderms.

2 MATERIAL AND METHOD

2.1 Gene identif ication and sequence analysis

The saponin biosynthesis enzymes were identif ied based on theA.japonicusgenome assembly (http://www.genedatabase.cn/aja_genome_20161129.html)from a previous publication (Zhang et al., 2017).Sequence alignments were carried out using the ClustalW algorithm in DNAMAN software (Lynnon Biosoft, USA). The identif ication of domains was carried out by SMART (Schultz et al., 2000) and the Batch CD-search tool in the Conserved Domain Database (Marchler-Bauer et al., 2010) of NCBI. The copy numbers of saponin biosynthesis enzymes from diff erent species were gained from the NCBI database.

2.2 Phylogenetic analysis and structural and function prediction analysis

Phylogenetic analysis was conducted using amino acid sequences of enzymes responsible for saponin biosynthesis fromA.japonicusand selected species retrieved from NCBI, Uniprot, and Ensembl,including sea urchin, sea star, and ginseng. Multiple alignments of protein sequences were conducted using the MUSCLE method (Edgar, 2004). The phylogenetic trees of the protein sequences were built using the maximum likelihood method and JTT+I+G model in MEGA 6 (Darriba et al., 2011; Tamura et al.,2013).

The 3D structural models were built based on the templates matching the target sequences in Swiss-Model (Guex and Peitsch, 1997). The models were visualized using the UCSF Chimera program(Pettersen et al., 2004). Sequence analysis and protein structure prediction were analyzed by PredictProtein(http://www.predictprotein.org) (Rost et al., 2004).The Gene Ontology (GO) terms were also predicted using PredictProtein (Rost et al., 2004).

2.3 Gene expression analysis

The expression patterns of crucial saponin biosynthesis enzymes during the developmental stages ofA.japonicuswere analyzed with our unpublished data. The RNA-seq datasets were obtained from f ive developmental stages ofA.japonicus(fertilized oocytes, blastula, gastrula,doliolaria, and penractula). Accurate transcript quantif ication was analyzed with RSEM (Li and Dewey, 2011) using default parameters (http://deweylab.github.io/RSEM/). A Heatmap of expression patterns was built using GraphPad Prism 8. The putative triterpene saponins biosynthetic pathway was plotted. It was based on literature searches on the biosynthesis of triterpenoid saponins in other species of sea cucumber, includingStichopushorrens,Eupentactafraudatrix, andHolothuriascabra(Makarieva et al., 1993; Silchenko et al., 2012;Mitu et al., 2017; Liu et al., 2018; Claereboudt et al.,2019).

3 RESULT

3.1 Phylogenetic and sequence analysis of saponin biosynthesis enzymes in A. japonicus

In theA.japonicusgenome, 30 genes responsible for saponin biosynthesis were found (Table 1):

Upstream enzymes [2 genes for AACTs (ACAT1 and ACAT2), 2 genes for HMGS (HMGS1 and HMGS2), and 2 genes for FPSs (FPS1 and FPS2)].

Midstream enzymes [2 genes for OSCs (OSC1 and OSC2), 1 gene for SS, 2 genes for SEs (SE1 and SE2)].

Downstream enzyme [19 genes for UGTs, 3 enzymes for UGT1 subfamily, 3 enzymes for UGT2A subfamily, 4 for UGT2B subfamily, 9 for UGT2C subfamily].

In the overall phylogenetic analysis (Figs.1 & 2),all of the saponin biosynthesis enzymes (i.e., AACTs,FPSs, HMGSs, OSCs, SS, SEs, and UGTs) inA.japonicuswere f irst clustered with their counterparts in other echinoderm species, namely sea star and sea urchin. These genes were then grouped with the orthologous genes in ginseng. Two types of AACT enzymes were identif ied inA.japonicusand sea star,i.e., cytosolic and mitochondrial-like biocatalysts.Detailed phylogenetic analysis revealed that these two types of AACT enzymes were clustered into two separate clades. InA.japonicus, most UGTs were grouped with homologous copies in the genome and clustered with the corresponding sea urchin genes.

The two types of AACTs, ACAT1 and ACAT2,both included a conserved domain of the Thiolase_N superfamily (Fig.3). Furthermore, FPS1, FPS2, and SS contained a conserved domain of the Isoprenoid_Biosyn_C1 superfamily. A conserved domain of the HMG_CoA_synt_C superfamily was identif ied in HMGS1 and HMGS2. In addition, OSC1 and OSC2 exhibited a conserved domain of the ISOPREN_C2_like superfamily, while SE1 and SE2 contained a conserved domain of squalene monooxygenase.Lastly, UGTs shared a conserved domain of the UDPglucuronosyltransferase (UDPGT) superfamily.

Analysis of the saponin biosynthesis enzymes revealed conserved tertiary structures comparable to those found in plants. The secondary structures of AACTs, HMGSs, SEs, and UGTs were composed of both α-helices and β-sheets. Notably, the secondary structures of FPSs, OSCs, and SSs displayed a relatively high level of α-helices and a low percentage of β-sheets (Fig.4 & Table 2). Specif ically, the secondary structure composition of ACAT2 was 34.13% helices and 11.74% strands. Moreover, the secondary structure of ACAT1 was made up of 46.32% helices and 8.59% strands. The FPSs were composed of 49.21%-55.06% helices and 2.86%-3.16% strands, whereas the secondary structures of the HMGSs were comprised of 33.60%-36.17%helices and 10.22%-12.77% strands. Additionally,the OSCs consisted of 36.53%-42.08% helices and 3.70%-5.94% strands. The SEs were made up of 27.76%-37.55% helices and 12.85%-21.49% strands.The SS enzyme was composed of 67.23% helices and 0.72% strands. The UGTs were comprised of 29.86%-48.50% helices and 9.01%-16.99% strands.

Fig.1 Phylogenetic analysis of enzymes responsible for saponin biosynthesis using sequences from A. japonicus, sea star, sea urchin, and ginseng

Fig.2 Phylogenetic analysis of UDP-glucuronosyltransferase using sequences from A. japonicus, sea star, sea urchin, and ginseng

Most of the enzymes responsible for the saponin biosynthesis were categorized into the corresponding molecular function ontology via homology to the annotated proteins (Table 2). ACAT1 and ACAT2 belonged to the acetyl-CoA C-acyltransferase activity(GO: 0003988). Furthermore, HMGSs, OSCs, SS,and UGTs were categorized into the hydroxymethylglutaryl-CoA synthase (GO:0004421), oxidosqualene cyclase (GO: 0031559),squalene synthase (GO: 0051996), and UDPglycosyltransferase (GO: 0008194) activity,respectively. Lastly, FPSs and SEs were classif ied into the transferase (GO: 0016740) and oxidoreductase(GO: 0016491) activity, correspondingly.

Table 1 Identif ication of enzymes responsible for saponin biosynthesis in A. japonicus

3.2 Copy numbers of saponin biosynthesis enzymes

Overall, the copy numbers of the enzymes responsible for the saponin biosynthesis signif icantly varied between organisms, particularly those of UGTs(Table 3). It is noteworthy that the copy numbers of SS in ginseng were higher than those in the echinoderm species. Some duplicated saponin biosynthesis enzymes in the genome ofA.japonicusmight be a result of independent tandem duplication events. Two copies of the lanosterol synthase genes (FPS1 and FPS2) were determined as tandemly arrayed genes on scaff old879. Moreover, multiple copies of UGT2C1 were identif ied in the genomeA.japonicus. Two copies (UGT2A1-1 and UGT2A1-2) and two additional paralogous copies (UGT2C1-8 and UGT2C1-9) were tandemly arrayed genes on scaff old28 and scaff old1550, respectively. Lastly,UGT2B37 and UGT2C1-4 were tandemly arrayed genes on scaff old285.

InA.japonicus, ACAT2 exhibited 57.55%sequence identity with the corresponding gene in the sea star. Notably, ACAT1 shared a 45.67% and 44.06% sequence identity with the homologous genes in the sea star and sea urchin. The HMGSs ofA.japonicusshared amino acid sequence similarities of 44.32%-64.35% and 43.16%-64.55% with the sea star and sea urchin. Additionally, the OSCs inA.japonicusdisplayed 29.40%-36.30% and 29.62%-36.30% sequence identity with the corresponding genes in the sea star and sea urchin. The SS enzyme inA.japonicusshowed sequence identity of 57.51% and 53.43% with the relevant genes in the sea star and sea urchin. The SEs inA.japonicusshared an amino acid identity of 26.83%-39.22% with the corresponding protein in the sea star. However, the SEs inA.japonicusshared a relatively low amino acid sequence similarity of approximately 1.36%-15.62%with the corresponding genes in the sea urchin. TheUGTs inA.japonicusexhibited low amino acid sequence similarities of approximately 9.69% and 11.36% with the sea star and sea urchin, respectively.Thus, the obtained results suggested that overall, the SEs and UGTs in A. japonicus shared low amino acid sequence similarities to the corresponding genes in other echinoderms.

Table 2 Summary of secondary structure component in saponin biosynthesis enzymes and their Gene Ontology term

Table 3 Summary of enzymes (copy numbers) for saponin biosynthesis in A. japonicus, sea urchin, sea star,and ginseng

Fig.3 Conserved domain analysis of enzymes responsible for saponin biosynthesis in A. japonicus

Similar toA.japonicus, we also determined seven types of enzymes responsible for saponin biosynthesis in ginseng, i.e., AACTs, FPSs, HMGSs, OSCs, SS,SEs, and UGTs. The ginseng ACAT2 shared a sequence similarity of 43.10% with the relevant gene inA.japonicus. Moreover, FPSs shared 48.23%-48.58% sequence similarity with the corresponding genes inA.japonicus. The ginseng HMGSs shared an amino acid similarity of 39.54%-52.06% with the homologous genes inA.japonicus. The OSCs in ginseng exhibited 33.33%-35.80% identities with the related genes inA.japonicus. The ginseng SEs shared a sequence identity of 37.28%-44.00% with the homologous genes inA.japonicus. Furthermore, the ginseng SS showed a 44.56%-48.89% amino acid sequence identity with the corresponding proteins inA.japonicus. It is recognized that sequence identities of most of the saponin biosynthesis enzymes have been evolutionarily conserved across the ginseng and marine organisms. Nonetheless, 19 and 225 UGTs were identif ied inA.japonicusand ginseng,respectively. Intriguingly, it was one of the largest gene families in ginseng (Xu et al., 2017); however,not in other plant species. The sequence similarities of UGTs betweenA.japonicusand ginseng ranged from 5.76% to 23.50%.

3.3 Expression of saponin biosynthesis enzymes during developmental stages

AACTs were expressed at the highest levels at the doliolaria stage (Fig.5a). Moreover, ACAT2 was expressed at a higher level than ACAT1 at all f ive developmental stages. The expression levels of FPSs only marginally changed during the f ive developmental stages. It is noteworthy that the expression of HMGS1 was the highest at the gastrula stage. Similarly, OSC1 also showed the highest expression at the gastrula stage, whereas OSC2 was highly expressed at the remaining four stages (fertilized oocytes, blastula,doliolaria, and penractula). SE1 and SE2 were highly expressed at the two early stages (fertilized oocytes and blastula). SS exhibited higher expression levels at the gastrula, doliolaria, and penractula stages. In contrast, UGT1-B7, UGT1-1, UGT1-7C, UGT2A1-1,and UGT2A1-3 showed relatively high expression levels at the fertilized oocytes, blastula, and gastrula stages (Fig.5b). The expression levels of UGTs were relatively low at all developmental stages. Among them, UGT2A1-1 showed the highest expression levels, while UGT2A1-2 was not detected in the analyzed RNA-seq datasets.

3.4 Putative mechanisms of saponin biosynthesis in A. japonicus

Fig.4 Tertiary structures of saponin biosynthesis enzymes in A. japonicus

In the present study, the putative mechanisms of saponin biosynthesis inA.japonicuswere predicted based on previous studies concerning triterpenoid saponins biosynthesis. As demonstrated in Fig.6, the saponin biosynthesis involves several important biological reactions. Initially, AACTs generate acetoacetyl-CoA from acetyl-CoA (Haapalainen et al., 2007). HMGSs subsequently catalyze the reaction between acetyl-CoA and acetoacetyl-CoA to generate HMG-CoA (Rokosz et al., 1994). HMG-CoA reductase then facilitates the synthesis of mevalonate from HMG-CoA (Luskey and Stevens, 1985).Furthermore, mevalonate-5-diphosphate is biosynthesized from mevalonate by the action of mevalonate kinase and phosphomevalonate kinase(Miziorko, 2011). Subsequently, isopentenyl-PP and dimethylallyl-PP are generated via the mevalonic acid(MVA) pathway (Chang et al., 2013). FPSs catalyze the reaction of dimethylallyl pyrophosphate to form farnesyl-PP (Qian et al., 2017), while SS facilitates the conversion of farnesyl-PP to squalene (Kim et al.,2011). SEs then catalyze the transformation of squalene to (S)-2,3-oxidosqualene (Zhao et al., 2010).InA.japonicus, OSCs exhibiting a plant-like motif produces parkeol instead of lanosterol (Li et al.,2018). Lastly, UGTs modify and produce diverse triterpenoid saponins (Seki et al., 2015).

4 DISCUSSION

Fig.5 RSEM-calculate expression levels of saponin biosynthesis enzymes in A. japonicus during developmental stages

Fig.6 Putative triterpene saponins biosynthetic pathway in A. japonicus

Saponins are a large family of amphiphilic glycosides of steroids and triterpenes (Yang et al.,2014). Although saponins have been widely identif ied and isolated from plants, their determination in aquaculture species has been limited. Many pharmacological activity studies have focused on the sea cucumber saponins (Chen et al., 2018; Ding et al.,2018; Kwon et al., 2018; Meng et al., 2018; Han et al.,2019). It has been found that enzymes, such as OSCs,P450s, UGTs, and several specialized transferases,are essential in the saponin biosynthesis (Adzhubei et al., 2013; Moses et al., 2014). A recent study on the evolution of OSC genes inA.japonicus(Li et al.,2018) as well as other work on the body wall transcriptome ofH.scabra(Mitu et al., 2017)indicated that sea cucumber exhibited a plant-like saponin biosynthesis pathway.

In this work, we identif ied and characterized 30 enzymes responsible for the saponin biosynthesis by analyzing the genome ofA.japonicus. The key saponin biosynthesis enzymes were classif ied into seven types, i.e., two AACTs, two FPSs, two HMGSs,two OSCs, one SS, two SEs, and 19 UGTs. The annotation of the enzymes responsible for saponin biosynthesis was based on the phylogenetic and sequence analyses. Sea cucumber, also called marine ginseng, displays similar appearance and medicinal properties to the plant ginseng (Bordbar et al., 2011).Both sea cucumber and ginseng contain a high abundance of potential anticancer compounds,including triterpene glycosides (Wang et al., 2007;Aminin et al., 2015). We found that the sequences of the upstream triterpenoid saponin biosynthesis genes in echinoderms and ginseng exhibited high sequence similarities. The comparison of animal and plant OSCs revealed that the OSCs inA.japonicuscontained plant-like motifs lacking in sea urchin and starf ish (Li et al., 2018). In addition, UGTs inA.japonicusshared relatively low amino acid sequence similarities with other echinoderms (9.69%and 11.36%), while those betweenA.japonicusand ginseng were 5.76%-23.50%. It is recognized that UGT enzymes participate in the generation of structurally diverse triterpenoid saponins via glycosylation reactions using a UDP-based sugar donor (Seki et al., 2015). It was previously reported that the tandem array of UGTs inBarbareavulgarisled to higher order bisdesmosidic glycosylation and the formation of various saponins (Erthmann et al.,2018). It was hypothesized that inA.japonicus, the low homology of UGTs with the corresponding genes in other echinoderms might result in the production of structurally and functionally diverse saponins. The conserved domains (Fig.3) as well as secondary structures of the key enzymes in the saponin biosynthesis pathway were evaluated based on their predicted protein sequences. AACTs, HMGSs, SEs,and UGTs were composed of both α-helices and β-sheets, whereas the secondary structures of FPSs,OSCs, and SS contained a relatively high level of α-helices and a low percentage of β-sheets.

Saponin biosynthesis is an important biological process occurring during the growth and development of numerous plants (Haralampidis et al., 2002).Nevertheless, research concerning the biosynthesis of saponins during the developmental stages of sea cucumber is limited. In this work, we determined that the enzymes involved in the saponin biosynthesis pathway exhibited dynamic expression patterns during the f ive critical stages of development ofA.japonicus(fertilized oocytes, blastula, gastrula,doliolaria, and penractula). Mitochondrial ACAT1 and cytosolic ACAT2 showed diff erential expression in various tissues of non-human primates (Lee et al.,2000). It was proposed that the two enzymes had diff erent functions, as ACAT1 was ubiquitous, while ACAT2 was specif ically expressed in the liver and intestine (Lee et al., 2000). InA.japonicus, ACAT2 was expressed at a higher level than ACAT1 at all f ive developmental stages. OSC1 showed higher expression at the gastrula stage, whereas OSC2 was highly expressed at the remaining four stages (i.e.,fertilized oocytes, blastula, doliolaria, and penractula).Conversely, inStichopushorrens, during the triterpenoid biosynthesis, OSC1 and OSC2 were highly expressed in the intestine (Liu et al., 2018).Both inA.japonicus(Li et al., 2018) andHolothuriaarenicola(Cordeiro et al., 1988), OSC1 and OSC2 were modif ied to synthesize parkeol instead of lanosterol. Furthermore, inA.japonicus, most UGTs were expressed at low levels during the critical development period, with the exception of the UGT1 family and UGT2A subfamily. UGT2A1-1 showed the highest expression levels. It is known that UGTs play essential roles in the structural diversity of saponins (Seki et al., 2015). Hence, higher expression levels of the UGT1 and UGT2A subfamilies might contribute to the formation of various saponins during the early development ofA.japonicus.

The MVA and non-mevalonate pathway (MEP)pathways are critical in saponin biosynthesis (Zhao et al., 2014). The MEP pathway is a plant-specialized process, while the MVA pathway plays an essential role in the saponin biosynthesis as well as the production of saponin precursors inA.japonicus. InA.japonicus, parkeol is biosynthesized by OSCs (Li et al., 2018), and then modif ied by various UGTs to produce diverse structural triterpene saponins. In the present study, we presented putative mechanisms for the biosynthesis of saponin inA.japonicus(Fig.6);however, further verif ication of the proposed mechanisms is required.

5 CONCLUSION

Thirty genes responsible for biosynthesis of saponin were identif ied inA.japonicus. These included upstream (AACTs, FPSs, HMGSs),midstream (OSCs, SEs, SS), and downstream (UGTs)enzymes. The genomic organization and evolution analysis of the entire set of saponin biosynthesis enzymes inA.japonicusprovided further understanding of crucial saponin biosynthesis components. It was determined that UGTs shared lower amino acid sequence similarity with the corresponding genes in other echinoderms than upstream and midstream enzymes, which could contribute to the production of more diverse types of saponins. The dynamic expression patterns of the saponin biosynthesis enzymes during developmental stages are important for understanding the role of saponin biosynthesis in the growth and development of sea cucumber. The putative mechanisms proposed herein provide a complete view of the saponin biosynthesis inA.japonicusand will be valuable for the evaluation of saponin biosynthesis in other echinoderms.

6 DATA AVAILABILITY STATEMENT

The datasets analyzed during the current study are available at NCBI under the accession number PRJNA354676.