Institute of Oceanology, Chinese Academy of Sciences
Article Information
- JI Jinlin, WANG Zhenhui, ZHU Wei, LI Qiang, WANG Yinan
- Proteomic analysis provides insights into the function of Polian vesicles in the sea cucumber Apostichopus japonicus post-evisceration
- Journal of Oceanology and Limnology, 40(5): 2028-2041
- http://dx.doi.org/10.1007/s00343-021-1198-3
Article History
- Received Jun. 28, 2021
- accepted in principle Sep. 23, 2021
- accepted for publication Dec. 4, 2021
2 Dalian Ocean University, Dalian 116023, China
The sea cucumber Apostichopus japonicus (phylum Echinodermata; class Holothuroidea; family Stichopodidae) is an important marine resource. Fisheries of this species support the livelihoods of coastal dwellers along the northwest coast of the Pacific Ocean from 35°N to 44°N, including in Russia, China, Japan, and Korea. Sea cucumber evisceration, which occurs when animals are exposed to intense stimuli or adverse environmental conditions, is well known (Wang et al., 2015). A consequence of this evisceration is that sea cucumbers undergo major tissue loss, metabolic rate depression, and immune system modification (Zang et al., 2012). Visceral regeneration following natural or induced evisceration is spatially and temporally separated from other developmental events. Evisceration provides an excellent model to study organogenesis and tissue regeneration, and therefore it has been extensively studied in, for example, research on regeneration of the intestine (Sun et al., 2013, 2017b), coelomocytes (Li et al., 2018a; Shi et al., 2020), and respiratory trees (Dolmatov and Ginanova, 2009).
Like other echinoderms, sea cucumbers possess a water vascular system that provides hydraulic pressure to the tentacles and tube feet, allowing them to move and collect food items (Liao, 1997). As the main accessory structures of the water vascular system in sea cucumber, Polian vesicles are bag-like accessory structures that consist of mesothelium, muscular coat, connective tissue, and ciliated epithelium, and a lumen that is filled with coelomic fluid (Baccetti and Rosati, 1968). Polian vesicles occur in different numbers in sea cucumbers, and they can hold water vascular fluid under slight pressure, until required. In addition to their function in the water vascular system, Polian vesicles are considered to play important roles in the inflammatory response under antigenic provocation, such as erythrocyte injection, and a phylogenetic relationship between sea cucumber Polian vesicles and the vertebrate lymphoreticular system has been proposed (Smith, 1978; D'Ancona et al., 1989). More recently, Polian vesicles were proved to be one of the hematopoietic tissues and to play a key role in the generation of coelomocytes in sea cucumber (Li et al., 2019). Despite these important functions, studies into the role of Polian vesicles in sea cucumber post-evisceration are lacking, which is an issue considering that the Polian vesicle is the only organ to remain within the coelom.
Proteomic analysis involves the systematic identification and quantification of the complete complement of proteins of a biological system. The quantitative measurements enable comparisons to be made between different samples to detect differentially abundant proteins (DAPs) associated with, for example, cellular responses to a range of events. Isobaric stable isotope labeling is being increasingly applied in proteomic analysis. Comparative proteomics has been widely used in studies on A. japonicus to investigate pathogen-host interactions and identify the proteins involved (Zhang et al., 2014; Lv et al., 2019), to reveal sex differences in humoral immunity and physiological characteristics associated with reproduction (Jiang et al., 2020), to detect regulatory mechanisms associated with intestine regeneration (Sun et al., 2017b), and to identify key proteins involved in stress resistance (Xu et al., 2016; Huo et al., 2019).
Previously, we analyzed the transcriptomes of A. japonicus Polian vesicles pre-evisceration and at 6-h post-evisceration to detect related gene expression patterns. Of 2 752 identified differentially expressed genes (DEGs), 1 453 were up-regulated, and 1 299 were down-regulated. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified various DEGs and pathways involved in cell biosynthesis and immunity (Shi et al., 2020). Here, we perform a Tandem Mass Tag (TMT)-based proteomics analysis of Polian vesicles by comparing protein abundances in samples obtained pre-evisceration, and at 6-h and 3-d post-evisceration, and combine this with published sea cucumber genomic data. The results greatly expands the knowledge of the functions of Polian vesicles in A. japonicus and provides new insights into the mechanisms of evisceration and regeneration.
2 MATERIAL AND METHOD 2.1 Sample preparationHealthy sea cucumbers (A. japonicus) weighing 60.1±10.2 g were purchased from an aquatic farm in Rizhao, Shandong Province, China, and acclimated to the laboratory environment in an aquarium with aerated seawater at 18±1 ℃ for 10 d. After acclimation, sea cucumbers were randomly selected and divided into three replicates (10 cucumbers in each) in each of three groups: pre-evisceration (PV0h), 6-h post-evisceration (PV6h), and 3-d post-evisceration (PV3d). Evisceration in groups PV6h and PV3d was induced by intracoelomic injection of 1.2 mL of 0.35-mol/L KCl. In all three groups the ventral surface of A. japonicus was opened longitudinally with an aseptic scalpel, and the Polian vesicle was excised (Fig. 1a). The vesicle wall was washed three times with 1×phosphate buffered saline (PBS). The Polian vesicles from the 10 sea cucumbers in each replicate for each group were combined into a single sample. In total, nine samples (3 replicates × 3 groups) were prepared and placed in liquid nitrogen before TMT sequencing.
2.2 Protein extraction and digestionEach sample was milled in liquid nitrogen, then lysed with lysis buffer (50-mmol/L Tris-HCL, pH 8.0, 8-mol/L urea, 0.2% sodium dodecyl sulfate (SDS)) and ultrasonicated for 5 min on ice. The supernatant of the lysate was collected after centrifugation (12 000×g, 15 min, 4 ℃) and the total protein content was measured using a Bradford Protein Assay Kit (Sangon, China). Sample extracts were reduced with 2-mmol/L dithiothreitol (DTT) (56 ℃) and alkylated in sufficient iodoacetamide (room temperature in darkness) for 1 h each. Four times the volume of precooled acetone was added to the samples and mixed by vortexing. After incubation at -20 ℃ for 2 h or more, the mixture was centrifuged (12 000×g, 15 min, 4 ℃), and the precipitated pellet was collected and washed twice using cold acetone, then dissolved in buffer containing 0.1-mol/L tetraethylammonium bromide (TEAB) (pH 8.5) and 8-mol/L urea. The protein concentration was measured using a Bradford Protein Assay Kit. Precisely 0.1-mg protein from each sample was digested for 16 h at 37 ℃ with 3-μg Trypsin Gold (Promega). The peptides were desalted using a C18 cartridge solid-phase extraction (SPE) from Sep-Pak (Waters, MA, USA) and then vacuum dried.
2.3 TMT labeling of peptides and high-performance liquid chromatography (HPLC) fractionationThe peptides were labeled with TMT 10-plex reagent using a TMT10plexTM Isobaric Label Reagent Set (Thermo Fisher) in accordance with the manufacturer's protocols. One unit of labeling reagent was used for 100-μg peptides. The peptides and labeling reagent were dissolved in 100 L of 0.1-mol/L TEAB and 41-μL acetonitrile, respectively, then co-incubated for 1 h. The reaction was terminated using ammonium hydroxide. Different peptide samples labeled by TMT were equally mixed and then desalted using a desalting spin column (Thermo Fisher, 89852). The labeled peptides were separated using a C18 column (Waters BEH C18 4.6 mm × 250 mm, 5 μm) connected to a Rigol L3000 HPLC system operating at 1 mL/min. The temperature of the column oven was set at 50 ℃. Mobile phases A (2% acetonitrile) and B (98% acetonitrile) were adjusted to pH 10.0 using ammonium hydroxide and used for gradient elution. The solvent gradient was set at: 3% B for 0 min; 3%– 8% B for 6 s; 8%–18% B for 11 min 54 s; 18%–32% B for 11 min; 32%–45% B for 7 min; 45%–80% B for 3 min; 80% B for 5 min; 80%–5% B for 6 s, and 5% B for 6 min 54 s. The tryptic peptides were monitored at 214 nm and the eluent was collected every minute, then merged into 15 fractions. All of the fractions were dried under vacuum and resuspended in 0.1% formic acid for subsequent assays.
2.4 Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysisShotgun proteomics analyses were performed using an EASY-nLC 1200 system (Thermo Fisher) coupled with an Orbitrap Q Exactive HF-X mass spectrometer (Thermo Fisher) operating in the data-dependent acquisition (DDA) mode. A total of 2-μg peptides reconstituted in 0.1% formic acid was injected onto an Acclaim PepMap100 C18 Nano-Trap column (2 cm×100 μm, 5 μm). Peptides were separated on a ReproSil-Pur 120 C18-AQ analytical column (15 cm ×150 μm, 1.9 μm) using a linear gradient of eluent B (0.1% formic acid in 80% acetonitrile) in eluent A (0.1% formic acid in H2O). The solvent gradient was set at: 5%–10% B for 2 min; 10%–40% B for 105 min; 40%–50% B for 5 min; 50%–90% B for 3 min; and 90%–100% B for 5 min. The separated peptides were analyzed using a Q-Exactive HF-X mass spectrometer with spray voltage 2.3 kV and capillary temperature 320 ℃. Full MS scans from 350–1 500 m/z were conducted at 60 000 resolution (at 200 m/z), with an automatic gain control (AGC) target value of 3×106 and maximum ion injection time of 20 ms. A maximum of 40 precursor ions with the highest abundance from the full MS scan was selected for high energy collision dissociation (HCD) fragment analysis at 45 000 resolution (at 200 m/z), with an AGC target value of 1×105, maximum ion injection time of 45 ms, normalized collision energy 32%, intensity threshold 8.3×105, and dynamic exclusion parameter 60 s. LC-MS/MS testing was performed by Novogene Co., Ltd. (Beijing, China).
2.5 Protein identification and quantificationThe resulting spectra were searched individually against the A. japonicus reference genome (National Center for Biotechnology Information GCA_002754855.1) using Proteome Discoverer 2.2 (Thermo Fisher). Proteins that contained at least one unique peptide were identified with a false discovery rate of < 0.01 for the peptide and protein. Reporter Quantification (TMT 10-plex) was used for TMT quantification. The significance of differences in the protein composition among different groups was determined by the Mann-Whitney U Test. P-value < 0.05 and fold change > 1.20 or < 0.83 were set as the thresholds for screening DAPs. The mean values of the abundances of the DAPs in each group were calculated and fuzzy C-means clustering of DAPs was analyzed using the "TCseq" package in R. The raw proteomics data was submitted to the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the iProX partner repository (Ma et al., 2019).
2.6 Protein annotation and functional analysis of the DAPsGene ontology (GO) and InterPro (IPR) analyses of the DAPs were performed using InterProScan 5 (Jones et al., 2014). The Clusters of Orthologous Groups (COG; http://www.geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/) databases were used to analyze protein families and pathways. GO and KEGG pathway enrichment analyses of the DAPs were performed and terms or pathways with P-values < 0.05 were considered significantly enriched. The fold change of annotated DAPs was visualized using the "heatmap" package in R. SignalP 4.1 (Petersen et al., 2011) and TMHMM (Krogh et al., 2001) were used to predict the presence of signal peptides and transmembrane domains, respectively, in the DAPs.
3 RESULT 3.1 Overview of the proteomics dataA total of 444 064 spectra were obtained by TMT-LC-MS/MS proteomic analysis from the nine samples (comprising 10 pooled sea cucumber Polian vesicles in each of three replicates in three groups). After filtering to eliminate low-scoring spectra, 49 863 peptides were obtained, with 8 453 proteins identified from the pre- and post-evisceration samples. Of the identified proteins, 8 120, 5 391, 3 555, and 7 368 proteins were assigned terms from the KEGG, GO, COG, and IPR databases, respectively (Fig. 1b). Most of the proteins (96%; 8 120/8 453) were assigned to KEGG pathways, including cellular processes, environmental information processing, genetic information processing, human diseases, and metabolism and organismal systems (Fig. 1c). The identification, annotation, and quantitative analysis details for these proteins are available in Supplementary Table S1.
3.2 DAPs identified post-eviscerationPairwise comparisons among the different groups detected 222 DAPs in the post-evisceration Polian vesicle samples, with P-values < 0.05 and fold change > 1.20 or < 0.83 (Fig. 2; Supplementary Table S2). Only 35 proteins were significantly differentially abundant in the PV6h vs. PV0h comparison, whereas 161 proteins were significantly differentially abundant in the PV3d vs. PV0h comparison. Another 71 proteins were significantly differentially abundant in the PV3d vs. PV6h comparison (Table 1). Many of the DAPs were common in at least two of the comparisons, and a number of DAPs were specific in the PV3d vs. PV6h comparison (Fig. 2). The fuzzy C-means clustering divided the 222 DAPs into six clusters, with 69, 40, 32, 13, 54, and 14 proteins in clusters 1–6, respectively (Fig. 3).
3.3 Functions of DAPs post-eviscerationAmong the 222 DAPs, 35 were unannotated, including 15 DAPs with limited information about their domains (Supplementary Table S3). Considering the protein descriptions from genome annotations and the annotations assigned by searches against the various databases, we manually classified the other 187 DAPs. Most of the annotated DAPs were associated with cell growth and proliferation, immune response and wound healing, substance transport and metabolism, cytoskeleton/cilia/flagella, extracellular matrix (ECM), energy production and conversion, protein synthesis and modification, and signal recognition and transduction (Fig. 4; Supplementary Table S4). The DAPs from the PV3d vs. PV0h comparison were associated with all these biology processes or functions. While the DAPs from the PV6h vs. PV0h comparison were mainly associated with substance transport and metabolism, immune response and wound healing, cell growth and proliferation, and cytoskeleton/cilia/flagella.
The SignalP 4.1 and TMHMM analyses detected 24 secretory proteins and three transmembrane proteins among the 222 DAPs (Supplementary Table S5). Among the secretory proteins, six were associated with immune response and wound healing, and five were associated with ECM. Seven of the secretory DAPs were unannotated proteins, and two of them had a Plasminogen-Apple-Nematode (PAN)/Apple domain. The annotated secretory and transmembrane DAPs are marked in Fig. 4.
3.4 GO and KEGG pathway enrichment analysis of DAPs after eviscerationThe GO functional enrichment analysis identified 19 DAPs in the PV6h group (compared with the PV0h group) that were assigned significantly enriched GO terms (P < 0.05); three were under the biological process (BP) category, seven were under the molecular function (MF) category, and one was under the cellular component (CC) category (Fig. 5a). Most enriched MF terms were associated with substance transport, particularly lipid transport. The MF term insulin-like growth factor binding (1 DAP) and BP term regulation of cell growth (1 DAP) were associated with cell growth. One KEGG pathway, fat digestion and absorption, was significantly enriched, and the down-regulated DAP annotated as ATP-binding cassette sub-family A member 2-like was among the DAPs annotated with this pathway.
In the PV3d vs. PV0h comparison, 95 DAPs were assigned 28 significantly enriched GO terms (P < 0.05) under the BP (7), MF (18), and CC (3) categories (Fig. 5c). The top five enriched GO terms were transporter activity (MF, 11 DAPs), receptor activity (MF, 8 DAPs), single-organism transport (BP, 8 DAPs), substrate-specific transporter activity (MF, 8 DAPs), and DNA binding (MF, 7 DAPs). Additional to lipid transport, multiple substance transport-related terms were enriched. The significantly enriched KEGG pathways (P < 0.05) were cholinergic synapse, metabolism of materials (glycerolipid metabolism, fat digestion and absorption, vitamin digestion and absorption), signaling pathways (estrogen, peroxisome proliferator-activated receptor (PPAR), erythroblastic leukemia viral oncogene B (ErbB), and nucleotide-binding oligomerization domain (NOD)-like receptor), and RNA polymerase.
In the PV3d vs. PV6h comparison, 38 DAPs were assigned 26 significantly enriched GO terms (P < 0.05) under the BP (6), CC (4), and MF (16) categories (Fig. 5b). The enriched GO terms included membrane (CC, 14 DAPs), integral component of membrane (CC, 9 DAPs), transport (BP, 7 DAPs), transporter activity (MF, 6 DAPs), and three terms, receptor activity (MF), single-organism transport (BP), and substrate-specific transporter activity (MF), had five DAPs each. Many DAPs were enriched in CC terms related to membrane. ECM-related terms were also enriched, as were multiple substance transport-related terms, including iron ion transport (BP, 1 DAP). The significantly enriched KEGG pathways (P < 0.05) were ferroptosis, glutamatergic synapse, signaling pathways (PPAR, gonadotropin-releasing hormone (GnRH), estrogen, and relaxin), cell adhesion (cell adhesion molecules, focal adhesion), proteoglycans in cancer, and glycerolipid metabolism.
4 DISCUSSIONThe Polian vesicle may have an important role in the survival and regeneration of sea cucumber because it is the only organ that remains in the coelom after evisceration. While complete visceral regeneration in sea cucumber requires 2 to 3 weeks, coelomocytes can recover to pre-evisceration levels within 6-h post-evisceration, and some intracavity tissues gradually regenerate from 3-d post-evisceration, indicating that notable changes occur at the early stage of regeneration (Zhang et al., 2017; Li et al., 2018a). In this study, we performed TMT-based proteomics of Polian vesicles to investigate the role of this vesicle in the early regeneration process.
Compared with pre-evisceration levels, some proteins present in the Polian vesicle showed differential abundance at 6-h and 3-d post-evisceration, indicating an active response of the Polian vesicle to the pressure of evisceration. Most of the annotated DAPs were associated with cell growth and proliferation, immune response and wound healing, substance transport and metabolism, cytoskeleton/cilia/flagella, ECM, energy production and conversion, protein synthesis and modification, and signal recognition and transduction. Compared with pre-evisceration levels, more DAPs were found at 3-d post-evisceration than at 6-h post-evisceration, and these DAPs were involved in multiple biological processes. Our results indicated that a wide range of biological processes was induced in Polian vesicles in response to evisceration and that the response was gradually strengthened.
During evisceration of A. japonicus, almost all the coelomocytes in the coelom cavity are lost along with the coelomic fluid and internal organs. However, the number of coelomocytes in the coelom cavity was found to recover rapidly at 6-h post-evisceration (Li et al., 2018a). Considering that the Polian vesicle is full of coelomic fluid as well as coelomocytes and connected to the coelom cavity through the water vascular system, we presumed that the Polian vesicle had an important role in this replenishment by pumping out the trapped coelomocytes. Furthermore, both inner and outer epithelia of the Polian vesicle have microvilli and cilia on their surfaces, which might facilitate the migration and flowing of coelomocytes. Consistent with this speculation, we found that several cilia/flagella-related proteins were significantly up-regulated at PV6h in this study, including serine/threonine-protein kinase Unc51-like kinase 4 (ULK4), which may be involved in cerebrospinal fluid circulation and ciliogenesis (Liu et al., 2016), and hydrocephalus-inducing protein-like (HYDIN), which is a ciliary protein with a unique role in ciliary motility (Smith, 2007). The migration of coelomocytes replenishes the coelomocytes in the coelom cavity as an emergency response; however, continual and enhanced hematopoiesis processes are necessary to make up for the huge loss of coelomocytes.
We previously demonstrated that the Polian vesicle, rete mirabile, and respiratory tree were fundamental tissues involved in hematopoiesis in A. japonicus under normal conditions (Li et al., 2019). Hemocyte proliferation is known to take place in specialized hematopoietic organs that provide an environment where undifferentiated blood stem cells are able to self-renew and, at the same time, generate offspring that differentiate into different blood cell types (Grigorian and Hartenstein, 2013). In a previous study, a number of DEGs related to cell proliferation and differentiation were found in Polian vesicles at 6-h post-evisceration compared with pre-evisceration levels, and these DEGs were enriched in several signaling pathways involved in cell proliferation and differentiation (Shi et al., 2020). Here, we found that many DAPs related to cell growth and proliferation were consistently and remarkably regulated post-evisceration, including up-regulated growth factor-related proteins, such as the insulin-like growth factor binding protein (IGFBP) (Wirthgen et al., 2016), transforming growth factor ig-h3 (beta ig-h3-like) (Oh et al., 2005), and growth factor receptor-bound protein 2 (Grb2L) (Ge et al., 2017), as well as a oncogene partner called fibroblast growth factor receptor 1 and oncogene partner 2 (FOP2) (Jin et al., 2011). In contrast, several proteins reported to suppress tumors or cell proliferation were down-regulated, including cub and sushi multiple domains 2 (Zhang and Song, 2014), tripartite motif-containing protein 56 (TRIM56) (Chen et al., 2018), phosphatase and tensin homolog (PTEN) (Gupta and Leslie, 2016), caspase recruitment domain-containing protein 9 (Card9) (Pan et al., 2020), and zinc finger CCHC domain-containing protein 10 (ZCCHC10) (Ohira et al., 2019). Additionally, GO terms and KEGG pathways associated with cell proliferation and differentiation were significantly enriched, such as ErbB receptor (Masunaga et al., 2017) and PPAR signaling pathway (Lee et al., 2015; Li et al., 2018b; Magadum and Engel, 2018). These results consistently suggested that an active hematopoiesis process was initiated in the Polian vesicle in response to evisceration. Notably, some DAPs, such as up-regulated serine/threonine-protein kinase 10 (Stk10), also known as lymphocyte-oriented kinase (LOK), (Walter et al., 2003) and apoptosis-stimulating of p53 protein 1 (ASPP1) (Yamashita et al., 2015) were reported to be present in high abundance in hematopoietic tissue or hematopoietic stem cells, further confirming the hematopoietic function of the Polian vesicle. Taken together, these results suggest that the Polian vesicle participates in coelomocyte restoration post-evisceration as both a storage organ and hematopoietic tissue.
The Polian vesicle was reported to be involved in inflammatory responses to some adverse stresses (Smith, 1978; D'Ancona and Rizzuto, 1989). When almost all the coelomocytes are excreted during evisceration, sea cucumber immunity is weakened, and therefore the Polian vesicle might play an important role in increasing immunity post-evisceration. We previously demonstrated that several immune-related genes were differentially expressed and that immune-related pathways were enriched in the Polian vesicle post-evisceration, indicating a role for the Polian vesicle in immunity in A. japonicus (Shi et al., 2020). Here, we reported various DAPs associated with immune response and wound healing in the Polian vesicle after evisceration, as well as the enrichment of immune-related NOD-like receptor signaling pathway (Platnich and Muruve, 2019). These results indicated a strengthened immune state in Polian vesicles responding to evisceration. In addition, most of the up-regulated DAPs predicted as secretory proteins, such as acute-phase protein α-macroglobulin (α2MG) (Varma et al., 2017), and mucosa protein Fcgbp and agglutinin DMBT1 (Houben et al., 2019), may also function in the coelom cavity, especially for wound healing after evisceration.
The functions of the Polian vesicle in coelomocytes migration, proliferation, immune response, and wound healing are based on essential supplies of substances and energy. Energy and nutrients are also extremely important for survival and regeneration of sea cumber because food intake is totally ceased due to intestine lose after evisceration. Various metabolic and physiological parameters were found to change dramatically post-evisceration, with lipids and carbohydrates being consumed as energy sources (Tan et al., 2008; Sun et al., 2017a). Here, we found that many proteins associated with substance transport and metabolism, and energy production and conversion, were significantly increased in the Polian vesicle, especially 3-d post-evisceration. These DAPs were predicted to be involved in transport and metabolism of multiple kinds of substances (amino acid, carbohydrate, coenzyme, inorganic ion, lipid, and secondary metabolites). Corresponding GO terms and KEGG pathways were also enriched, including those associated with lipid transport and metabolism. These results indicate an augment metabolic progression in the Polian vesicle response to evisceration, which would contribute to the urgent needs of survival and regeneration. The sea cucumber Polian vesicle was also considered to be a primitive "excretory" organ that may have adapted to facilitate the passage of waste-filled amoebocytes through its tissues into the coelom (Baccetti and Rosati, 1968). Hence, we speculated that the Polian vesicle might facilitate excretion of dysfunctional metabolites and disintegrating coelomocytes in response to an adverse stimulus such as evisceration.
The ECM is recognized to play key regulatory roles because it orchestrates cell signaling, and various functions, properties, and morphology. Among them, proteins/glycoproteins such as collagens and tenascins, and their cell receptors such as integrins, are involved in cell adhesion. Some enzymes such as matrix metalloproteinases (MMPs), a family of endopeptidases involved in the degradation of various proteins in ECM, and their natural inhibitors, tissue inhibitors of metalloproteinases (TIMPs), play important roles in cell proliferation, migration, differentiation angiogenesis, and immunity regulation (Cui et al., 2017; Dolmatov et al., 2021; Karamanos et al., 2021). We found that ECM-related proteins were significantly regulated in Polian vesicle post-evisceration, indicating an active ECM remodeling response to evisceration, which might be related to coelomocyte proliferation and migration, and immunity regulation in Polian vesicles. Beside these, several DAPs might be secreted and then play roles in regeneration. Dedifferentiation of specialized cells in the remnant of organ and the migration of dedifferentiated cells are fundamental for regeneration in sea cucumber (Wang et al., 2015; Dolmatov, 2021). In these processes, ECM transformation occurred to promote the migration of dedifferentiated cells and the formation of new tissues. Two MMP genes, ajMMP-2 and ajMMP-16, have been cloned in A. japonicus, and were thought to form a road for migrating cells and accelerate cell migration (Miao et al., 2017). Because these genes were expressed only at the initial stage of regeneration, it was suggested that they may be involved in degradation of the esophageal ECM and dedifferentiation of coelomic epithelial cells and enterocytes (Miao et al., 2017; Dolmatov et al., 2021). In the present study, several members of MMP and TIMP with varied regulation trends were found, indicating a complicated process of ECM remodeling. However, the up-regulation and down-regulation of secreted proteins MMP16-like and TIMP-like in PV3d respectively suggested that they might contribute to visceral regeneration through promoting dedifferentiation and migration.
It should be mentioned that a number of unannotated DAPs were found in the present study, especially seven that were identified as secretory proteins. These DAPs might be specific to echinoderms and may play important roles in the recovery of sea cucumber post-evisceration. Further estimation of their functions using more methods is needed.
5 CONCLUSIONThe Polian vesicle exhibited a wide range of biological functions in response to evisceration in A.japonicus, and the response gradually strengthened. The Polian vesicle played an indispensable role in the restoration of coelomocytes after evisceration both as a storage and hematopoietic tissue. Immune processes were elevated in the Polian vesicle, and some DAPs may be secreted into the coelom cavity and contribute to wound healing after evisceration. To support these functions, augment metabolic progression was induced in Polian vesicle and was supposed to meet the urgent needs of nutrition and energy for survival and regeneration. In addition, several ECM-related DAPs may be secreted and then play roles in regeneration, especially by promoting dedifferentiation and migration. Collectively, these results provide new insights into the functions of Polian vesicles in A. japonicus post-evisceration.
6 DATA AVAILABILITY STATEMENTThe raw proteomic data have been deposited in the ProteomeXchange Consortium via the iProX partner repository, with the dataset identifier PXD025117.
Electronic supplementary material
Supplementary material (Supplementary Tables S1–S5) is available in the online version of this article at https://doi.org/10.1007/s00343-021-1198-3.
Baccetti B, Rosati F. 1968. The fine structure of the Polian vesicles of Holothurians. Zeitschrift für Zellforschung und Mikroskopische Anatomie, 90(1): 148-160.
DOI:10.1007/BF00496708 |
Chen Y, Zhao J, Li D Z, et al. 2018. TRIM56 suppresses multiple myeloma progression by activating TLR3/TRIF signaling. Yonsei Medical Journal, 59(1): 43-50.
DOI:10.3349/ymj.2018.59.1.43 |
Cui N, Hu M, Khalil R A. 2017. Biochemical and biological attributes of matrix metalloproteinases. Progress in Molecular Biology and Translational Science, 147: 1-73.
DOI:10.1016/bs.pmbts.2017.02.005 |
D'Ancona G, Rizzuto L, Canicattí C. 1989. Engagement of Polian vesicles during Holothuria Polii response to erythrocyte injection. Journal of Invertebrate Pathology, 54(3): 300-305.
DOI:10.1016/0022-2011(89)90113-4 |
Dolmatov I Y, Ginanova T T. 2009. Post-autotomy regeneration of respiratory trees in the holothurian Apostichopus japonicus (Holothuroidea, Aspidochirotida). Cell and Tissue Research, 336(1): 41-58.
DOI:10.1007/s00441-009-0761-6 |
Dolmatov I Y, Nizhnichenko V A, Dolmatova L S. 2021. Matrix metalloproteinases and tissue inhibitors of metalloproteinases in echinoderms: structure and possible functions. Cells, 10(9): 2331.
DOI:10.3390/cells10092331 |
Dolmatov I Y. 2021. Molecular aspects of regeneration mechanisms in holothurians. Genes, 12(2): 250.
DOI:10.3390/genes12020250 |
Ge S F, Xiong Y, Wu X P, et al. 2017. Role of growth factor receptor-bound 2 in CCl4-induced hepatic fibrosis. Biomedicine & Pharmacotherapy, 92: 942-951.
DOI:10.1016/j.biopha.2017.05.142 |
Grigorian M, Hartenstein V. 2013. Hematopoiesis and hematopoietic organs in arthropods. Development Genes and Evolution, 223(1-2): 103-115.
DOI:10.1007/s00427-012-0428-2 |
Gupta A, Leslie N R. 2016. Controlling PTEN (phosphatase and tensin homolog) stability: a dominant role for lysine66. Journal of Biological Chemistry, 291(35): 18465-18473.
DOI:10.1074/jbc.M116.727750 |
Houben T, Harder S, Schlüter H, et al. 2019. Different forms of TFF3 in the human saliva: heterodimerization with IgG Fc Binding Protein (FCGBP). International Journal of Molecular Science, 20(20): 5000.
DOI:10.3390/ijms20205000 |
Huo D, Sun L N, Zhang L B, et al. 2019. Global-warming-caused changes of temperature and oxygen alter the proteomic profile of sea cucumber Apostichopus japonicus. Journal of Proteomics, 193: 27-43.
DOI:10.1016/j.jprot.2018.12.020 |
Jiang J W, Zhao Z L, Pan Y J, et al. 2020. Proteomics reveals the gender differences in humoral immunity and physiological characteristics associated with reproduction in the sea cucumber Apostichopus japonicus. Journal of Proteomics, 217: 103687.
DOI:10.1016/j.jprot.2020.103687 |
Jin Y X, Zhen Y, Haugsten E M, et al. 2011. The driver of malignancy in KG-1a leukemic cells, FGFR1OP2-FGFR1, encodes an HSP90 addicted oncoprotein. Cellular Signalling, 23(11): 1758-1766.
DOI:10.1016/j.cellsig.2011.06.010 |
Jones P, Binns D, Chang H Y, et al. 2014. InterProScan 5: genome-scale protein function classification. Bioinformatics, 30(9): 236-1240.
DOI:10.1093/bioinformatics/btu031 |
Karamanos N K, Theocharis A D, Piperigkou Z, et al. 2021. A guide to the composition and functions of the extracellular matrix. FEBS Journal, 288(24): 6850-6912.
DOI:10.1111/febs.15776 |
Krogh A, Larsson B, Von Heijne G, et al. 2001. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. Journal of Molecular Biology, 305(3): 567-580.
DOI:10.1006/jmbi.2000.4315 |
Lee H Y, Gao X F, Barrasa M I, et al. 2015. PPAR-α and glucocorticoid receptor synergize to promote erythroid progenitor self-renewal. Nature, 522(7557): 474-477.
DOI:10.1038/nature14326 |
Li Q, Ren Y, Liang C L, et al. 2018a. Regeneration of coelomocytes after evisceration in the sea cucumber, Apostichopus japonicus. Fish & Shellfish Immunology, 76: 266-271.
DOI:10.1016/j.fsi.2018.03.013 |
Li Q, Ren Y, Luan L L, et al. 2019. Localization and characterization of hematopoietic tissues in adult sea cucumber, Apostichopus japonicus. Fish & Shellfish Immunology, 84: 1-7.
DOI:10.1016/j.fsi.2018.09.058 |
Li Y, Jin D X, Xie W X, et al. 2018b. PPAR-γ and Wnt regulate the differentiation of MSCs into adipocytes and osteoblasts respectively. Current Stem Cell Research & Therapy, 13(3): 185-192.
DOI:10.2174/1574888X12666171012141908 |
Liao Y L. 1997. Fauna Sinica, Echinodermata, Holothurioidea. Science Press, Beijing. 334p. (in Chinese)
|
Liu M, Guan Z L, Shen Q, et al. 2016. Ulk4 is essential for ciliogenesis and CSF flow. Journal of Neuroscience, 36(29): 7589-7600.
DOI:10.1523/JNEUROSCI.0621-16.2016 |
Lv Z M, Guo M, Li C H, et al. 2019. Divergent proteomics response of Apostichopus japonicus suffering from skin ulceration syndrome and pathogen infection. Comparative Biochemistry and Physiology Part D: Genomics and Proteomics, 30: 196-205.
DOI:10.1016/j.cbd.2019.03.003 |
Ma J, Chen T, Wu S F, et al. 2019. iProX: an integrated proteome resource. Nucleic Acids Research, 47(D1): D1211-D1217.
DOI:10.1093/nar/gky869 |
Magadum A, Engel F B. 2018. PPARβ/δ: linking metabolism to regeneration. International Journal of Molecular Sciences, 19(7): 2013.
DOI:10.3390/ijms19072013 |
Masunaga H, Sugimoto Y, Magi S, et al. 2017. Robustness analysis of the detailed kinetic model of an ErbB signaling network by using dynamic sensitivity. PLoS One, 12(5): e0178250,.
DOI:10.1371/journal.pone.0178250 |
Miao T, Wan Z X, Sun L N, et al. 2017. Extracellular matrix remodeling and matrix metalloproteinases (ajMMP-2 like and ajMMP-16 like) characterization during intestine regeneration of sea cucumber Apostichopus japonicus. Comparative Biochemistry and Physiology Part B: Biochemistry and Molecular Biology, 212: 12-23.
DOI:10.1016/j.cbpb.2017.06.011 |
Oh J E, Kook J K, Min B M. 2005. βig-h3 induces keratinocyte differentiation via modulation of involucrin and transglutaminase expression through the integrin α3β1 and the phosphatidylinositol 3-kinase/Akt signaling pathway. Journal of Biological Chemistry, 280(22): 21629-21637.
DOI:10.1074/jbc.M412293200 |
Ohira T, Kojima H, Kuroda Y, et al. 2019. PITX1 protein interacts with ZCCHC10 to regulate hTERT mRNA transcription. PLoS One, 14(8): e0217605.
DOI:10.1371/journal.pone.0217605 |
Pan L Y, Tan Y T, Wang B, et al. 2020. Caspase recruitment domain containing protein 9 suppresses non-small cell lung cancer proliferation and invasion via inhibiting MAPK/p38 pathway. Cancer Research and Treatment, 52(3): 867-885.
DOI:10.4143/crt.2019.606 |
Petersen T N, Brunak S, Von Heijne G, et al. 2011. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nature Methods, 8(10): 785-786.
DOI:10.1038/nmeth.1701 |
Platnich J M, Muruve D A. 2019. NOD-like receptors and inflammasomes: a review of their canonical and non-canonical signaling pathways. Archives of Biochemistry and Biophysics, 670: 4-14.
DOI:10.1016/j.abb.2019.02.008 |
Shi W B, Zhang J L, Wang Y N, et al. 2020. Transcriptome analysis of sea cucumber (Apostichopus japonicus) Polian vesicles in response to evisceration. Fish & Shellfish Immunology, 97: 108-113.
DOI:10.1016/j.fsi.2019.12.016 |
Smith A C. 1978. A proposed phylogenetic relationship between sea cucumber Polian vesicles and the vertebrate lymphoreticular system. Journal of Invertebrate Pathology, 31(3): 353-357.
DOI:10.1016/0022-2011(78)90228-8 |
Smith E F. 2007. Hydin seek: finding a function in ciliarymotility. Journal of Cell Biology, 176(4): 403-404.
DOI:10.1083/jcb.200701113 |
Sun L N, Sun J C, Xu Q Z, et al. 2017a. Metabolic responses to intestine regeneration in sea cucumbers Apostichopus japonicus. Comparative Biochemistry and Physiology Part D: Genomics and Proteomics, 22: 32-38.
DOI:10.1016/j.cbd.2017.02.003 |
Sun L N, Xu D X, Xu Q Z, et al. 2017b. iTRAQ reveals proteomic changes during intestine regeneration in the sea cucumber Apostichopus japonicus. Comparative Biochemistry and Physiology Part D: Genomics and Proteomics, 22: 39-49.
DOI:10.1016/j.cbd.2017.02.004 |
Sun L N, Yang H S, Chen M Y, et al. 2013. RNA-seq reveals dynamic changes of gene expression in key stages of intestine regeneration in the sea cucumber Apostichopus japonicas. PLoS One, 8(8): e69441.
DOI:10.1371/journal.pone.0069441 |
Tan Y K, Li X, Duan J J. 2008. Energy metabolism and biochemical composition of Apostichopus japonicus during the regeneration of viscera. Journal of Fishery Sciences of China, 15(4): 683-688.
DOI:10.3321/j.issn:1005-8737.2008.04.021.(inChinesewithEnglishabstract) |
Varma V R, Varma S, An Y, et al. 2017. Alpha-2 macroglobulin in Alzheimer's disease: a marker of neuronal injury through the RCAN1 pathway. Molecular Psychiatry, 22(1): 13-23.
DOI:10.1038/mp.2016.206 |
Walter S A, Cutler Jr R E, Martinez R, et al. 2003. Stk10, a new member of the polo-like kinase kinase family highly expressed in hematopoietic tissue. Journal of Biological Chemistry, 278(20): 18221-18228.
DOI:10.1074/jbc.M212556200 |
Wang T M, Sun L N, Chen M Y. 2015. Chapter 11 - Aestivation and regeneration. Developments in Aquaculture and Fisheries Science, 39: 177-209.
DOI:10.1016/B978-0-12-799953-1.00011-8 |
Wirthgen E, Höflich C, Spitschak M, et al. 2016. Quantitative Western ligand blotting reveals common patterns and differential features of IGFBP-fingerprints in domestic ruminant breeds and species. Growth Hormone & IGF Research, 26: 42-49.
DOI:10.1016/j.ghir.2015.11.001 |
Xu D X, Sun L N, Liu S L, et al. 2016. Understanding the heat shock response in the sea cucumber Apostichopus japonicus, using iTRAQ-based proteomics. International Journal of Molecular Sciences, 17(2): 150.
DOI:10.3390/ijms17020150 |
Yamashita M, Nitta E, Suda T. 2015. Maintenance of hematopoietic stem cell integrity and regulation of leukemogenesis by p53 and its coactivator Aspp1. The Japanese Journal of Clinical Hematology, 56(12): 2426-2433.
DOI:10.11406/rinketsu.56.2426 |
Zang Y Q, Tian X L, Dong S L, et al. 2012. Growth, metabolism and immune responses to evisceration and the regeneration of viscera in sea cucumber, Apostichopus japonicus. Aquaculture, 358-359: 50-60.
DOI:10.1016/j.aquaculture.2012.06.007 |
Zhang P, Li C H, Li Y, et al. 2014. Proteomic identification of differentially expressed proteins in sea cucumber Apostichopus japonicus coelomocytes after Vibrio splendidus infection. Developmental & Comparative Immunology, 44(2): 370-377.
DOI:10.1016/j.dci.2014.01.013 |
Zhang R, Song C. 2014. Loss of CSMD1 or 2 may contribute to the poor prognosis of colorectal cancer patients. Tumor Biology, 35(5): 4419-4423.
DOI:10.1007/s13277-013-1581-6 |
Zhang X J, Sun L N, Yuan J B, et al. 2017. The sea cucumber genome provides insights into morphological evolution and visceral regeneration. PLoS Biology, 15(10): e2003790.
DOI:10.1371/journal.pbio.2003790 |