Mechanistic study of transcription factor Sox18 during heart development.

Heart development is a delicate and complex process regulated by coordination of various signaling pathways. In this study, we investigated the role of sox18 in heart development by modulating Wnt/ß-Catenin signaling pathways. Our spatiotemporal expression analysis revealed that sox18 is mainly expressed in the heart, branchial arch, pharyngeal arch, spinal cord, and intersegmental vessels at the tailbud stage of Xenopus tropicalis embryo. Overexpression of sox18 in the X. tropicalis embryos causes heart edema, while loss-of-function of sox18 can change the signal of developmental heart marker gata4 at different stages, suggesting that sox18 plays an essential role in the development of the heart. Knockdown of SOX18 in human umbilical vein endothelial cells suggests a link between Sox18 and ß-CATENIN, a key regulator of the Wnt signaling pathway. Sox18 negatively regulates islet1 and tbx3, the downstream factors of Wnt/ß-Catenin signaling, during the linear heart tube formation and the heart looping stage. Taken together, our findings highlight the crucial role of Sox18 in the development of the heart via inhibiting Wnt/ß-Catenin signaling.

SOX9-COL9A3-dependent regulation of choroid plexus epithelial polarity governs blood-cerebrospinal fluid barrier integrity.

The choroid plexus (CP) is an extensively vascularized neuroepithelial tissue that projects into the brain ventricles. The restriction of transepithelial transport across the CP establishes the blood-cerebrospinal fluid (CSF) barrier that is fundamental to the homeostatic regulation of the central nervous system microenvironment. However, the molecular mechanisms that control this process remain elusive. Here we show that the genetic ablation of Sox9 in the hindbrain CP results in a hyperpermeable blood-CSF barrier that ultimately upsets the CSF electrolyte balance and alters CSF protein composition. Mechanistically, SOX9 is required for the transcriptional up-regulation of Col9a3 in the CP epithelium. The reduction of Col9a3 expression dramatically recapitulates the blood-CSF barrier defects of Sox9 mutants. Loss of collagen IX severely disrupts the structural integrity of the epithelial basement membrane in the CP, leading to progressive loss of extracellular matrix components. Consequently, this perturbs the polarized microtubule dynamics required for correct orientation of apicobasal polarity and thereby impedes tight junction assembly in the CP epithelium. Our findings reveal a pivotal cascade of SOX9-dependent molecular events that is critical for construction of the blood-CSF barrier.

Temporal Quantitative Proteomic and Phosphoproteomic Profiling of SH-SY5Y and IMR-32 Neuroblastoma Cells during All-Trans-Retinoic Acid-Induced Neuronal Differentiation.

The human neuroblastoma cell lines SH-SY5Y and IMR-32 can be differentiated into neuron-like phenotypes through treatment with all-trans-retinoic acid (ATRA). After differentiation, these cell lines are extensively utilized as in vitro models to study various aspects of neuronal cell biology. However, temporal and quantitative profiling of the proteome and phosphoproteome of SH-SY5Y and IMR-32 cells throughout ATRA-induced differentiation has been limited. Here, we performed relative quantification of the proteomes and phosphoproteomes of SH-SY5Y and IMR-32 cells at multiple time points during ATRA-induced differentiation. Relative quantification of proteins and phosphopeptides with subsequent gene ontology analysis revealed that several biological processes, including cytoskeleton organization, cell division, chaperone function and protein folding, and one-carbon metabolism, were associated with ATRA-induced differentiation in both cell lines. Furthermore, kinase-substrate enrichment analysis predicted altered activities of several kinases during differentiation. Among these, CDK5 exhibited increased activity, while CDK2 displayed reduced activity. The data presented serve as a valuable resource for investigating temporal protein and phosphoprotein abundance changes in SH-SY5Y and IMR-32 cells during ATRA-induced differentiation.

Direct Interaction of Sox10 With Cadherin-19 Mediates Early Sacral Neural Crest Cell Migration: Implications for Enteric Nervous System Development Defects.

BACKGROUND AND AIMS: The enteric nervous system, which regulates many gastrointestinal functions, is derived from neural crest cells (NCCs). Defective NCC migration during embryonic development may lead to enteric neuropathies such as Hirschsprung's disease (hindgut aganglionosis). Sox10 is known to be essential for cell migration but downstream molecular events regulating early NCC migration have not been fully elucidated. This study aimed to determine how Sox10 regulates migration of sacral NCCs toward the hindgut using Dominant megacolon mice, an animal model of Hirschsprung's disease with a Sox10 mutation. METHODS: We used the following: time-lapse live cell imaging to determine the migration defects of mutant sacral NCCs; genome-wide microarrays, site-directed mutagenesis, and whole embryo culture to identify Sox10 targets; and liquid chromatography and tandem mass spectrometry to ascertain downstream effectors of Sox10. RESULTS: Sacral NCCs exhibited retarded migration to the distal hindgut in Sox10-null embryos with simultaneous down-regulated expression of cadherin-19 (Cdh19). Sox10 was found to bind directly to the Cdh19 promoter. Cdh19 knockdown resulted in retarded sacral NCC migration in vitro and ex vivo, whereas re-expression of Cdh19 partially rescued the retarded migration of mutant sacral NCCs in vitro. Cdh19 formed cadherin-catenin complexes, which then bound to filamentous actin of the cytoskeleton during cell migration. CONCLUSIONS: Cdh19 is a direct target of Sox10 during early sacral NCC migration toward the hindgut and forms cadherin-catenin complexes which interact with the cytoskeleton in migrating cells. Elucidation of this novel molecular pathway helps to provide insights into the pathogenesis of enteric nervous system developmental defects.

Proteomic analysis of zebrafish folliculogenesis identifies YB-1 (Ybx1/ybx1) as a potential gatekeeping molecule controlling early ovarian folliculogenesis.

As in mammals, ovarian folliculogenesis in teleosts also consists of two phases: the primary growth (PG) and secondary growth (SG) phases, which are analogous to the preantral and antral phases respectively in mammals. In this study, we performed a proteomic analysis on zebrafish follicles undergoing the PG-SG transition aiming to identify factors involved in the event. Numerous proteins showed significant changes, and the most prominent one was Y-box binding protein 1 (YB-1; Ybx1/ybx1), a transcription factor and mRNA-binding protein. YB-1 belongs to the Y-box binding protein family, which also includes the gonad-specific YB-2. Interestingly, phylogenetic analysis showed no YB-2 homolog in zebrafish. Although ybx1 mRNA was expressed in various tissues, its protein Ybx1 was primarily produced in the gonads, similar to YB-2 in other species. In the ovary, Ybx1 protein started to appear in early follicles newly emerged from the germ cell cysts, reached the highest level in late PG oocytes, but decreased precipitously when the follicles entered the SG phase. In PG follicles, Ybx1 might function as a key component of the messenger ribonucleoprotein particles (mRNPs) in association with other RNA-binding proteins. Similar to mammalian YB-1, zebrafish Ybx1 also contains functional signals that determine its intracellular localization. In conclusion, Ybx1 may play dual roles of YB-1 and YB-2 in zebrafish. In the ovary, Ybx1 binds mRNAs to stabilize them while preventing their translation. At PG-SG transition, Ybx1 is removed to release the masked mRNAs for translation into functional proteins, leading to follicle activation.

Millipede genomes reveal unique adaptations during myriapod evolution.

The Myriapoda, composed of millipedes and centipedes, is a fascinating but poorly understood branch of life, including species with a highly unusual body plan and a range of unique adaptations to their environment. Here, we sequenced and assembled 2 chromosomal-level genomes of the millipedes Helicorthomorpha holstii (assembly size = 182 Mb; shortest scaffold/contig length needed to cover 50% of the genome [N50] = 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (assembly size = 449 Mb, N50 = 26.78 Mb mainly on 17 pseudomolecules). Unique genomic features, patterns of gene regulation, and defence systems in millipedes, not observed in other arthropods, are revealed. Both repeat content and intron size are major contributors to the observed differences in millipede genome size. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute (AGO) proteins for loading small RNAs are duplicated in both millipedes, but unlike in insects, an AGO duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs-including species-specific microRNA arm switching-providing differential gene regulation is also obtained. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules found in centipedes. We identify sets of genes associated with the ozadene that play roles in chemical defence as well as antimicrobial activity. Macro-synteny analyses revealed highly conserved genomic blocks between the 2 millipedes and deuterostomes. Collectively, our analyses of millipede genomes reveal that a series of unique adaptations have occurred in this major lineage of arthropod diversity. The 2 high-quality millipede genomes provided here shed new light on the conserved and lineage-specific features of millipedes and centipedes. These findings demonstrate the importance of the consideration of both centipede and millipede genomes-and in particular the reconstruction of the myriapod ancestral situation-for future research to improve understanding of arthropod evolution, and animal evolutionary genomics more widely.

Combination of untargeted and targeted proteomics for secretome analysis of L-WRN cells.

Organoid culture is a promising biomedical technology that requires specialized growth factors. Recently, a recombinant L-WRN cell line has been extensively used to generate conditioned medium (L-CM) for organoid culture. Nevertheless, methods for evaluating the stability of the L-WRN cells have been limited. In this study, a novel proteomics-based approach was developed to analyze the secretome of the cells. Serum-free L-CM was lyophilized, precipitated by trichloroacetic acid, and desalted prior to analysis by liquid chromatography-tandem mass spectrometry. Data-dependent acquisition (DDA) was conducted for the untargeted secretome profiling of the cells, and parallel reaction monitoring (PRM) was applied for the targeted quantification of the Wnt3A, R-spondin3, and noggin proteins (WRNs). This study also compared the performance of two types of PRM methods, namely MS1-independent PRM and MS1-dependent PRM, that can be executed on an Orbitrap instrument. The results showed that the growth of mouse intestinal organoids was closely related to the use of L-CM. The composition of L-CM could be markedly affected by the medium collection scheme. A total of 1725, 2302, and 2681 proteins were identified from the L-CM collected on day 5, day 9, and day 13, respectively. The MS1-independent PRM outperformed the MS1-dependent PRM and effectively quantified the WRNs with high repeatability and specificity. In conclusion, by integrating untargeted and targeted proteomics, this study develops a mass spectrometry-based method for the secretome analysis and quality control of the L-WRN cells. The methodology and findings of the present work will benefit future studies on organoids and secretomes.

Comprehending the allergen repertoire of shrimp for precision molecular diagnosis of shrimp allergy.

BACKGROUND: Clinical management of shrimp allergy is hampered by the lack of accurate tests. Molecular diagnosis has been shown to more accurately reflect the clinical reactivity but the full spectrum of shrimp allergens and their clinical relevance are yet to be established. We therefore sought to comprehend the allergen repertoire of shrimp, investigate and compare the sensitization pattern and diagnostic value of the allergens in allergic subjects of two distinct populations. METHODS: Sera were collected from 85 subjects with challenge-proven or doctor-diagnosed shrimp allergy in Hong Kong and Thailand. The IgE-binding proteins of Penaeus monodon were probed by Western blotting and identified by mass spectrometry. Recombinant shrimp allergens were synthesized and analyzed for IgE sensitization by ELISA. RESULTS: Ten IgE-binding proteins were identified, and a comprehensive panel of 11 recombinant shrimp allergens was generated. The major shrimp allergens among Hong Kong subjects were troponin C (Pen m 6) and glycogen phosphorylase (Pen m 14, 47.1%), tropomyosin (Pen m 1, 41.2%) and sarcoplasmic-calcium binding protein (Pen m 4, 35.3%), while those among Thai subjects were Pen m 1 (68.8%), Pen m 6 (50.0%) and fatty acid-binding protein (Pen m 13, 37.5%). Component-based tests yielded significantly higher area under curve values (0.77-0.96) than shrimp extract-IgE test (0.70-0.75). Yet the best component test differed between populations; Pen m 1-IgE test added diagnostic value only in the Thai cohort, whereas sensitizations to other components were better predictors of shrimp allergy in Hong Kong patients. CONCLUSION: Pen m 14 was identified as a novel shrimp allergen predictive of challenge outcome. Molecular diagnosis better predicts shrimp allergy than conventional tests, but the relevant component is population dependent.

Identification of Diverse Stress-Responsive Xylem Sap Peptides in Soybean.

Increasing evidence has revealed that plant secretory peptides are involved in the long-distance signaling pathways that help to regulate plant development and signal stress responses. In this study, we purified small peptides from soybean (Glycine max) xylem sap via o-chlorophenol extraction and conducted an in-depth peptidomic analysis using a mass spectrometry (MS) and bioinformatics approach. We successfully identified 14 post-translationally modified peptide groups belonging to the peptide families CEP (C-terminally encoded peptides), CLE (CLAVATA3/embryo surrounding region-related), PSY (plant peptides containing tyrosine sulfation), and XAP (xylem sap-associated peptides). Quantitative PCR (qPCR) analysis showed unique tissue expression patterns among the peptide-encoding genes. Further qPCR analysis of some of the peptide-encoding genes showed differential stress-response profiles toward various abiotic stress factors. Targeted MS-based quantification of the nitrogen deficiency-responsive peptides, GmXAP6a and GmCEP-XSP1, demonstrated upregulation of peptide translocation in xylem sap under nitrogen-deficiency stress. Quantitative proteomic analysis of GmCEP-XSP1 overexpression in hairy soybean roots revealed that GmCEP-XSP1 significantly impacts stress response-related proteins. This study provides new insights that root-to-shoot peptide signaling plays important roles in regulating plant stress-response mechanisms.

PrSM-Level Side-by-Side Comparison of Online LC-MS Methods with Intact Histone H3 and H4 Proteoforms.

The heterogeneity of histone H3 proteoforms makes histone H3 top-down analysis challenging. To enhance the detection coverage of the proteoforms, performing liquid chromatography (LC) front-end to mass spectrometry (MS) detection is recommended. Here, using optimized electron-transfer/high-energy collision dissociation (EThcD) parameters, we have conducted a proteoform-spectrum match (PrSM)-level side-by-side comparison of reversed-phase LC-MS (RPLC-MS), "dual-gradient" weak cation-exchange/hydrophilic interaction LC-MS (dual-gradient WCX/HILIC-MS), and "organic-rich" WCX/HILIC-MS on the top-down analyses of H3.1, H3.2, and H4 proteins extracted from a HeLa cell culture. While both dual-gradient WCX/HILIC and organic-rich WCX/HILIC could resolve intact H3 and H4 proteoforms by the number of acetylations, the organic-rich method could enhance the separations of different trimethyl/acetyl near-isobaric H3 proteoforms. In comparison with RPLC-MS, both of the WCX/HILIC-MS methods enhanced the qualities of the H3 PrSMs and remarkably improved the range, reproducibility, and confidence in the identifications of H3 proteoforms.

Analysis of Soybean Long Non-Coding RNAs Reveals a Subset of Small Peptide-Coding Transcripts.

Long non-coding RNAs (lncRNAs) are defined as non-protein-coding transcripts that are at least 200 nucleotides long. They are known to play pivotal roles in regulating gene expression, especially during stress responses in plants. We used a large collection of in-house transcriptome data from various soybean (Glycine max and Glycine soja) tissues treated under different conditions to perform a comprehensive identification of soybean lncRNAs. We also retrieved publicly available soybean transcriptome data that were of sufficient quality and sequencing depth to enrich our analysis. In total, RNA-sequencing data of 332 samples were used for this analysis. An integrated reference-based, de novo transcript assembly was developed that identified ∼69,000 lncRNA gene loci. We showed that lncRNAs are distinct from both protein-coding transcripts and genomic background noise in terms of length, number of exons, transposable element composition, and sequence conservation level across legume species. The tissue-specific and time-specific transcriptional responses of the lncRNA genes under some stress conditions may suggest their biological relevance. The transcription start sites of lncRNA gene loci tend to be close to their nearest protein-coding genes, and they may be transcriptionally related to the protein-coding genes, particularly for antisense and intronic lncRNAs. A previously unreported subset of small peptide-coding transcripts was identified from these lncRNA loci via tandem mass spectrometry, which paved the way for investigating their functional roles. Our results also highlight the present inadequacy of the bioinformatic definition of lncRNA, which excludes those lncRNA gene loci with small open reading frames from being regarded as protein-coding.

Proteomic Analysis of the Venom of Jellyfishes Rhopilema esculentum and Sanderia malayensis.

Venomics, the study of biological venoms, could potentially provide a new source of therapeutic compounds, yet information on the venoms from marine organisms, including cnidarians (sea anemones, corals, and jellyfish), is limited. This study identified the putative toxins of two species of jellyfish-edible jellyfish Rhopilema esculentum Kishinouye, 1891, also known as flame jellyfish, and Amuska jellyfish Sanderia malayensis Goette, 1886. Utilizing nano-flow liquid chromatography tandem mass spectrometry (nLC-MS/MS), 3000 proteins were identified from the nematocysts in each of the above two jellyfish species. Forty and fifty-one putative toxins were identified in R. esculentum and S. malayensis, respectively, which were further classified into eight toxin families according to their predicted functions. Amongst the identified putative toxins, hemostasis-impairing toxins and proteases were found to be the most dominant members (>60%). The present study demonstrates the first proteomes of nematocysts from two jellyfish species with economic and environmental importance, and expands the foundation and understanding of cnidarian toxins.

Hepatocellular carcinoma-derived exosomes promote motility of immortalized hepatocyte through transfer of oncogenic proteins and RNAs.

Exosomes are increasingly recognized as important mediators of cell-cell communication in cancer progression through the horizontal transfer of RNAs and proteins to neighboring or distant cells. Hepatocellular carcinoma (HCC) is a highly malignant cancer, whose metastasis is largely influenced by the tumor microenvironment. The possible role of exosomes in the interactions between HCC tumor cell and its surrounding hepatic milieu are however largely unknown. In this study, we comprehensively characterized the exosomal RNA and proteome contents derived from three HCC cell lines (HKCI-C3, HKCI-8 and MHCC97L) and an immortalized hepatocyte line (MIHA) using Ion Torrent sequencing and mass spectrometry, respectively. RNA deep sequencing and proteomic analysis revealed exosomes derived from metastatic HCC cell lines carried a large number of protumorigenic RNAs and proteins, such as MET protooncogene, S100 family members and the caveolins. Of interest, we found that exosomes from motile HCC cell lines could significantly enhance the migratory and invasive abilities of non-motile MIHA cell. We further demonstrated that uptake of these shuttled molecules could trigger PI3K/AKT and MAPK signaling pathways in MIHA with increased secretion of active MMP-2 and MMP-9. Our study showed for the first time that HCC-derived exosomes could mobilize normal hepatocyte, which may have implication in facilitating the protrusive activity of HCC cells through liver parenchyma during the process of metastasis.

Synergistic induction of apoptosis by methylseleninic acid and cisplatin, the role of ROS-ERK/AKT-p53 pathway.

Cisplatin-based therapy is one of the most important chemotherapy treatments for cancers. However, its efficacy is greatly limited by drug resistance and undesirable side effects. Therefore, it is of great importance to develop chemosensitizing agents to cisplatin. In the present study, we demonstrated the strategy to use methylseleninic acid (MeSe) as a synergistic agent of cisplatin and elucidated their action mechanisms. The combination of MeSe and cisplatin exhibited synergistic anticancer efficacy and achieved greater selectivity between cancer cell and normal cell. By inducing intracellular oxidative stress, MeSe potentiated cisplatin-induced DNA damage and led to enhanced p53 phosphorylation, followed by increased activation of both mitochondrial and death receptor pathway. Down-regulation of phosphorylated AKT and ERK also played important roles in the synergistic effects of MeSe and cisplatin. Our results suggested that the strategy to apply MeSe as a synergistic agent to cisplatin could be a highly efficient way to achieve anticancer synergism by targeting the intracellular redox system. MeSe might be a candidate for clinical application as a chemosensitizer to cisplatin-based therapy for cancer treatments, especially for hepatocellular carcinoma.

A novel selenadiazole derivative induces apoptosis in human glioma cells by dephosphorylation of AKT.

Selenadiazole derivatives are synthetic organoselenium compounds with improved anticancer activity and greater selectivity than inorganic selenium. In this study, 4-(benzo[c][1,2,5]selenadiazol-6-yl)-benzene-1,2-diamine (BSBD) was shown to induce time- and dose-dependent apoptosis in SWO-38 human glioma cells by accumulation of a sub-G1 cell population, DNA fragmentation, nuclear condensation, caspase activation and poly(ADP-ribose) polymerase (PARP) cleavage. Further mechanistic investigation showed that BSBD treatment induced dephosphorylation of AKT and DNA damage-mediated activation of p53, leading to extensive apoptosis through the mitochondrial pathway. Our findings suggest that BSBD represents a potential human glioma therapeutic.

A reference-grade genome of the xerophyte Ammopiptanthus mongolicus sheds light on its evolution history in legumes and drought-tolerance mechanisms.

Plants that grow in extreme environments represent unique sources of stress-resistance genes and mechanisms. Ammopiptanthus mongolicus (Leguminosae) is a xerophytic evergreen broadleaf shrub native to semi-arid and desert regions; however, its drought-tolerance mechanisms remain poorly understood. Here, we report the assembly of a reference-grade genome for A. mongolicus, describe its evolutionary history within the legume family, and examine its drought-tolerance mechanisms. The assembled genome is 843.07 Mb in length, with 98.7% of the sequences successfully anchored to the nine chromosomes of A. mongolicus. The genome is predicted to contain 47 611 protein-coding genes, and 70.71% of the genome is composed of repetitive sequences; these are dominated by transposable elements, particularly long-terminal-repeat retrotransposons. Evolutionary analyses revealed two whole-genome duplication (WGD) events at 130 and 58 million years ago (mya) that are shared by the genus Ammopiptanthus and other legumes, but no species-specific WGDs were found within this genus. Ancestral genome reconstruction revealed that the A. mongolicus genome has undergone fewer rearrangements than other genomes in the legume family, confirming its status as a "relict plant". Transcriptomic analyses demonstrated that genes involved in cuticular wax biosynthesis and transport are highly expressed, both under normal conditions and in response to polyethylene glycol-induced dehydration. Significant induction of genes related to ethylene biosynthesis and signaling was also observed in leaves under dehydration stress, suggesting that enhanced ethylene response and formation of thick waxy cuticles are two major mechanisms of drought tolerance in A. mongolicus. Ectopic expression of AmERF2, an ethylene response factor unique to A. mongolicus, can markedly increase the drought tolerance of transgenic Arabidopsis thaliana plants, demonstrating the potential for application of A. mongolicus genes in crop improvement.

Comparative proteomics analysis of selenium responses in selenium-enriched rice grains.

By foliar fortification with selenite, selenium (Se)-enriched rice with a higher Se content and grain yield has been generated. However, the regulatory mechanisms of Se response in rice grains remain unknown; therefore, we carried out a comparative proteomics study in Se-enriched rice grains by using two approaches including two-dimensional gel electrophoresis (2-DE)-coupled MALDI-TOF/TOF MS and 1-DE/LC-FTICR-MS-coupled label-free quantification. By comparison between Se treatment and control, 62 and 250 abundance changed proteins were identified from 2-DE and 1-DE, respectively. By functional classification, proteins involved in metabolism, cell redox regulation, and seed nutritional storage were the most highly affected by Se accumulation. The up-regulation of late embryogenesis abundant proteins as well as proteins involved in sucrose synthesis and other metabolism pathways may contribute to the earlier maturation and higher yield of the Se-enriched rice. In addition, there have been six proteins identified to contain selenoamino acid modification, which is the first identification of selenoproteins in higher plants. In conclusion, our study provided novel insights into Se response in rice grains at the proteome level, which are expected to be highly useful for dissecting the Se response pathways in rice and for the production of Se-enriched rice in the future.

Genome of the sea anemone Exaiptasia pallida and transcriptome profiles during tentacle regeneration.

Cnidarians including sea anemones, corals, hydra, and jellyfishes are a group of animals well known for their regeneration capacity. However, how non-coding RNAs such as microRNAs (also known as miRNAs) contribute to cnidarian tissue regeneration is poorly understood. Here, we sequenced and assembled the genome of the sea anemone Exaiptasia pallida collected in Hong Kong waters. The assembled genome size of E. pallida is 229.21 Mb with a scaffold N50 of 10.58 Mb and BUSCO completeness of 91.1%, representing a significantly improved genome assembly of this species. The organization of ANTP-class homeobox genes in this anthozoan further supported the previous findings in jellyfishes, where most of these genes are mainly located on three scaffolds. Tentacles of E. pallida were excised, and both mRNA and miRNA were sequenced at 9 time points (0 h, 6 h, 12 h, 18 h, 1 day, 2, 3, 6, and 8 days) from regenerating tentacles. In addition to the Wnt signaling pathway and homeobox genes that are shown to be likely involved in tissue regeneration as in other cnidarians, we have shown that GLWamide neuropeptides, and for the first time sesquiterpenoid pathway genes could potentially be involved in the late phase of cnidarian tissue regeneration. The established sea anemone model will be useful for further investigation of biology and evolution in, and the effect of climate change on this important group of animals.