The role of retinoic acid signaling in starfish metamorphosis
© The Author(s) 2018
Received: 6 February 2018
Accepted: 12 April 2018
Published: 21 April 2018
Although retinoic acid (RA) signaling plays a crucial role in the body patterning of chordates, its function in non-chordate invertebrates, other than its mediation of environmental cues triggering metamorphosis in cnidarians, is largely unknown. We investigated the role of RA signaling in the metamorphosis of starfish (Echinodermata).
We found that exogenous RA treatment induced metamorphosis in starfish larvae. In contrast, inhibitors of RA synthesis and RA receptors suppressed metamorphosis triggered by attachment to a substrate. Gene expressions of the RA signaling component were detected in competent larvae.
This study provides insight into the ancestral function of RA signaling, which is conserved in the metamorphosis of cnidarians and starfish.
Retinoic acid (RA) plays a critical role in the body patterning of chordates, such as the anterior–posterior patterning of the central nervous system and pharyngeal arches [1–3]. RA synthesized by retinal dehydrogenase (Raldh) regulates downstream gene expression through binding to the retinoic acid receptor (RAR) and the retinoid x receptor (RXR) heterodimer nuclear receptor. Although RXR was identified in various metazoan taxa including cnidarian, arthropod, and nematode , Raldh and RAR, as well as cytochrome P450 26 (CYP26), which degrades RA, had been described only in chordates, and thus RA signaling was thought to be specific to chordates, and the acquisition of the gene families for RA signaling was thought to be a key step in the evolution of the chordate body plan . However, Cañestro et al. reported that gene families of Raldh, RAR, and CYP26 are encoded in the genomes of non-chordate deuterostomes . Following additional genomic surveys of other invertebrates, the origin of RA signaling was pushed back to the common ancestor of metazoans [6, 7].
Despite these genomic surveys, the function of RA signaling in non-chordate deuterostomes remains largely unknown, other than the observation of pseudopodial cable growth on micromere-delivered cells from Hemicentrotus pulcherrimus after RA treatment . Sciarrino et al. reported that RA treatment on Paracentrotus lividus did not induce specific phenotype except for delaying development . Although Marlétaz et al. suggested that the promoter region of sea urchin Hox genes had an RA response element , no additional evidence of the effect of RA on Hox gene expression has been reported.
Insight into the ancestral function of RA signaling can be gained from the study of cnidarians, a basal lineage of animals. Fuchs et al. showed that RA signaling is involved in the life cycle transition of cnidarians . Medusozoan cnidarians have two distinct life cycle phases. Fertilized eggs first develop into planula larvae and then, through further development, become sessile polyps. When the polyps receive environmental stimuli, such as a temperature change, they metamorphose into medusae . Fuchs et al. reported that RA is an internal regulator of this process and suggested that RA works through binding with RXR . Because RXR is known to be involved in the metamorphosis of insects and frogs [11–16], RA signaling was suggested to also have a role in the life cycle regulation of metazoans.
Molecular mechanisms for metamorphosis have been investigated in echinoderms, primarily echinoids [17–28]. In the sea urchin Strongylocentrotus purpuratus, larvae acquired competence to metamorphose 4.5–6 weeks after fertilization  and commenced metamorphosis in response to substrates such as algal sheets . Previous studies suggested that thyroid hormone (TH) and histamine (HA) functioned as modulators of larval growth and the competent state [19–21, 24–27]. Although several works also revealed nitric oxide (NO) signaling functioned as a negative regulator upon settlement [17, 18], any role of RA signaling in metamorphosis is not reported in any taxa of echinoderm.
Here, we investigate the role of RA signaling in asteroid metamorphosis. We used the starfish Patiria pectinifera, as it acquires competence to metamorphose earlier than sea urchin species (about 1 week after fertilization). We show that exogenous RA treatment induced larval body absorption and juvenile rudiment development in competent larvae. Conversely, the inhibition of endogenous RA synthesis or binding of RA to RAR blocked metamorphosis. In addition, genes related to RA signaling were expressed in the juvenile and larval stages. These results suggest that RA signaling acts as a regulator of metamorphosis after settlement.
Sampling and culture
We collected adult specimens of P. pectinifera around the mainland of Japan: Manaduru (Kanagawa Prefecture), Asamushi (Aomori Prefecture), and Hiraiso (Ibaraki Prefecture). Fertilization and embryo rearing were performed from described previously . From 2 days post-fertilization (dpf), the larvae were fed Chaetoceros calcitrans (purchased as Sun-culture: Marine Tec, Aichi) by rearing in seawater with 50,000 cells/ml. We changed seawater with fresh algae every other day.
The RA signaling pathway was activated by exogenous all-trans RA (Sigma-Aldrich, St Louis, CAS number: 302-79-4). RA signaling was inhibited by the Raldh inhibitor N,N-diethylaminobenzaldehyde (DEAB, Tokyo Chemical Industry, Tokyo, CAS number: 120-21-8) or the RAR antagonist RO41-5253 (RO, Focus Biomolecules, Pennsylvania, CAS number: 144092-31-9). We prepared 100 mM, 1 M, and 50 mM stocks of RA, DEAB and RO, respectively, in dimethyl sulfoxide (DMSO). Larvae were incubated in 2 ml artificial seawater containing 2 µl of reagents or DMSO in 12-well plates at 22 °C. For the treatment without a substrate, 10 larvae were incubated in one well. In the experiments conducted to induce metamorphosis with a substrate, a single larva was cultured in one well to identify individuals. We observed metamorphosis proceed in the same manner in these densities of larvae. Coral sand from tanks of adult P. pectinifera was used as a substrate for P. pectinifera experiments . For cases in which drug treatment continued for more than 2 days, we changed the seawater with the same concentration of drugs every other day.
The metamorphosis ratio was calculated by dividing the number of metamorphosis commenced larvae by the number of treated larvae. In the case of reagent treatment with substrate, metamorphosis ratio was calculated by dividing the number of metamorphosed larvae by the number of settled larvae.
We used ANOVA to evaluate differences of reagent treatments on the metamorphosis or settlement, assuming above conditions and batches as main factors and “blocks,” respectively, because preliminary tests in some treatments revealed slight, but not essential, differences among batches. If we had two factors, we assessed an interaction of the two factors. Also, multiple comparison was performed by Tukey’s honestly significant difference test. The R statistical software was used for all statistical analyses .
As to the assumption of normal distribution, we use Kolmogorov–Smirnov test and confirmed no significant difference from the assumption of normality. Also, Bartlett’s test was used to test homogeneity of variances prior to ANOVA. In a few experiments, we found significant violations of the equal variance assumption and conducted Welch’s one-way ANOVA to test main effects without the assumption of equality. However, there were no differences in statistical conclusions between the main factor of ANOVA and the Welch’s method. We thus used ANOVA with assumption of equal variance in the following analyses because we needed multi-dimensional analyses simultaneously and because ANOVA is robust to deviations from equal variance, especially if numbers of sample are equal , and normality assumptions .
Identification of genes related to RA signaling in starfish
To gather sequence of RA signaling related genes, we performed transcriptome analysis on samples of blastula (13 h post-fertilization), brachiolaria larva (14 dpf) and adult epidermis in starfish P. pectinifera. Each sample was collected from three biological replicates: The samples of blastula and brachiolaria larva were collected from progeny of three distinct pairs of parents, and the adult epidermis was collected from three mother specimens. Total RNA was extracted using TRIzol reagent (Life Technologies, Carlsbad, CA, USA) and treated with DNase and cleaned up using the RNeasy kit (Qiagen, Hilden, Germany). Paired-end libraries (200 bp) were prepared using a TruSeq RNA Sample Preparation Kit (Illumina). Sequencing on the Hiseq 2000 platform (Illumina, San Diego, CA, USA) was performed at the National Institute of Basic Biology (NIBB) The raw reads were deposited in the DDBJ Sequence Reads Archives (DRA006662). The low-quality reads were filtered using the NGS QC Toolkit57 . De novo assembly was conducted by Trinity58 (version: trinityrnaseq_r20140717) .
We recovered nucleotide sequences coding for raldh, rar and rxr from gene models assembled from our de novo transcriptome data of larval stages for P. pectinifera. The orthologies of these genes were confirmed by molecular phylogenetic analysis (Additional file 1: Fig. S1, Additional file 2: Fig. S2). All sequences are shown in Additional file 3: Table S1 and aligned using MAFFT ver. 7 (https://mafft.cbrc.jp/alignment/server/) with the default parameters . Amino acid sites for tree construction were selected by trimAl using a gap threshold value of 0.8 . A best-fitting amino acid substitution model and a maximum likelihood tree were inferred using RAxML 8.2.0 . Confidence values were calculated after 1000 bootstrap runs.
Whole-mount in situ hybridization
cDNA was obtained by reverse transcription from larval RNA using PrimeScript 1st strand cDNA Synthesis kit (Takara, Shiga, Japan). Additional file 4: Table S2 shows the sequences of the primers used for the amplification of raldha, raldhb, raldhc, rar, and rxr. We used 40-bp reverse primers including a 20-bp T3 promoter sequence to synthesize Dig-labelled RNA probes for in situ hybridization. Dig-labelled RNA probes for raldha, raldhb, raldhc, rar, and rxr were transcribed from PCR products using T3 RNA polymerase (Roche, Mannheim, Germany). In situ hybridization was performed as described previously .
Competence for metamorphosis was acquired around 8 dpf in P. pectinifera
As in other echinoderm species, the transition from larval form to adult form proceeds via multiple steps in starfish. Development of the adult rudiment commences well before larvae acquire competence. Adult skeletogenesis was observed in penta-radial manner as early as 6 dpf, when development of brachiolar arms commences. Development of the brachiolar arms is tightly linked with the acquisition of competence. Murabe et al. indicated that adhesive papillae on the brachiolar arms sense environmental cues such as proper substrate for promoting metamorphosis . When competent larvae attach to proper substrate, they initiate metamorphosis. The metamorphosis is observed as absorption of larval body and development of pentaradial juvenile body with primary podia (the precursor structures of tube feet), which occur within 24 h after settlement (Fig. 1).
Although development of the brachiolar arms commences at 6 dpf, we do not know when the competence to receive cues for this metamorphosis is acquired. Thus, we first sought to identify the stage at which P. pectinifera larvae become competent to metamorphose upon attachment to a substrate.
Exogenous RA treatments with competent P. pectinifera larvae induced metamorphosis
These results indicate that all-trans RA has the potential to induce larval body absorption and juvenile rudiment development when added to competent larvae at 14 dpf.
The effect of RA treatment on metamorphosis was limited after larvae became competent
We investigated whether RA can also affect the timing of larval competence to respond to the metamorphosis cue. No larva treated with 1 µM all-trans RA at 5–7 dpf commenced metamorphosis, as judged by larval body absorption and juvenile rudiment development (Fig. 3f). RA treatment of larvae younger than 8 dpf had no apparent effect on morphology or behavior before 8 dpf. However, some of these larvae cultured in seawater with RA commenced metamorphosis at 8 dpf (Fig. 3f). Regardless of when the RA treatment started, larvae responded to RA and metamorphosed at 8 or 9 dpf, which is comparable with the stage at which larvae acquire competence to metamorphose during normal development (Figs. 2, 3f). Most larvae that were treated with 1 µM RA after 8 dpf commenced metamorphosis within 24 h (Fig. 3f). These timelines are similar to those induced by a substrate. Above results indicate that RA does not affect the development of competence for metamorphosis, but rather functions as a cue to commence metamorphosis when added to competent larvae.
RAR activation was required for the initiation of metamorphosis
These results suggest that exogenous RA drives metamorphosis through binding to RAR. Therefore, we investigated whether RAR plays any role when larvae respond to a substrate and commence metamorphosis. We treated competent larvae with 1 µM RO at 14 dpf. As shown above, more than 90% of 14-dpf larvae responded to the introduction of coral sand and commenced metamorphosis (Fig. 4k). However, when 1 µM RO was added together with the substrate to 14-dpf larvae, metamorphosis was significantly inhibited (P = 0.000 and 0.632 in RO and batch, respectively, ANOVA; Fig. 4k right, Additional file 6: Table S4). Interestingly, in contrast to metamorphosis, RO treatment did not affect the number of larvae attached to the substrate (P = 0.256 and 0.058 in RO and batch, respectively, ANOVA; Fig. 4k left, Additional file 7: Table S5). These results are consistent with the idea that RAR functions as an internal mediator of the signaling cue for metamorphosis received by the adhesive papillae .
Endogenous RA synthesis was essential for metamorphosis after settlement
To assess whether endogenous RA synthesis was required for metamorphosis, we treated 14-dpf larvae with DEAB, a known inhibitor of Raldh (Fig. 4l). We cultured larvae in seawater with 100 or 300 µM DEAB and coral sand as a substrate. DEAB-treated larvae crawled around the substrate surface, and most settled to the substrate with the brachiolar arms, as normal competent larvae do. The ratio of settled larvae in the DEAB treatments did not differ significantly (P = 0.728, ANOVA; Fig. 4m left, Additional file 8: Table S6), although significant difference was observed among batches (P = 0.010 in batch, ANOVA). However, DEAB significantly inhibited metamorphosis after settlement. Significance of DEAB was found in metamorphosis (P = 0.000 and 0.001 in DEAB and batch, respectively, ANOVA; Fig. 4m right, Additional file 9: Table S7). The metamorphosis rate was slightly lower in the 100 µM DEAB treatment than in the DMSO treatment. Inhibition of 300 µM DEAB was significantly effective with DMSO and 100 µM DEAB treatments (P = 0.000 and P = 0.000, Tukey’s test, respectively), although there was no significant difference between the latter two in metamorphosis (P = 0.624, Tukey’s test). These results indicate that endogenous RA synthesis was essential for the commencement of metamorphosis after receipt of the signaling cue by the adhesive papillae, whereas it was not essential for settlement on the substrate.
Temporal and spatial expression patterns of genes related to RA signaling in larvae
The results reported above suggest that RA signaling functions as an endogenous mediator of a metamorphosis cue upon settlement. To support this idea, we surveyed RA signaling genes in our transcriptome data for P. pectinifera. We investigated their orthologues by construction of phylogenetic trees (Additional file 1: Fig. S1, Additional file 2: Fig. S2). We identified three raldh genes in our transcriptome data: raldha, raldhb and raldhc (Additional file 1: Fig. S1). Single genes for rar and rxr were also identified from our transcriptome data (Additional file 2: Fig. S2).
RA signaling acts as a regulator of starfish metamorphosis
Of note, the expression of raldhs, rar, and rxr was observed well before competent larvae received environmental cues (Fig. 5a–o). These expression patterns suggest the presence of another regulatory mechanism that controls the onset of RA signaling upon attachment. Murabe et al. indicated that environmental cues were received by adhesive papillae on the brachiolar arms . Thus, the onset of RA signaling may be regulated by something in the adhesive papillae, likely through the nervous system. NO and HA are candidates for this regulator, as both have inhibitory effects on the onset of metamorphosis in sea urchins [11, 17, 24, 25]. Although no evidence of NO signaling or HA signaling has been reported for starfish, further studies on the crosstalk of these signals with RA signaling will facilitate a comprehensive understanding of the molecular regulatory mechanism of echinoderm metamorphosis.
The role of the nuclear receptor/RXR heterodimer in metamorphosis regulation is evolutionarily conserved
RA signaling mediated by a RAR/RXR heterodimer receptor plays a critical role in the development of the unique body plans of chordates, such as the anterior–posterior patterning of the central nervous system [1, 2, 7]. Therefore, the ancestral role of RA signaling has attracted the interest of several zoologists [5, 6]. Our findings suggest that RA signaling plays a role in the regulation of metamorphosis in ancestral deuterostomes. During the evolution of ancestral chordates, RA signaling may have been co-opted for the patterning of the chordate body plan via the recruitment of Hox genes as regulators. Further investigation of RA signaling in non-chordate deuterostomes may shed new light on the evolution of the vertebrate body plan.
The transition of distinct body plans through metamorphosis is widely observed in metazoans [11, 41], and several recent studies have revealed that conserved molecular regulatory mechanisms underlie metamorphosis [10, 11, 42]. Molecular mechanisms for metamorphosis in amphibians and insects have been investigated in detail [12–16]. In both taxa, metamorphosis is regulated by hormones received by receptors (TH for amphibians and an ecdysone receptor for insects), which make a heterodimer with RXR. A recent study by Fuchs indicated that RXR is also involved in cnidarian metamorphosis . Exogenous treatment with RA was shown to induce the metamorphic process of strobilation. Here, we present evidence that RA signaling is also involved in starfish metamorphosis, and provide additional evidence that components of RA signaling are conserved in this metamorphosis. However, amphibians and insects use different hormones for signaling, and RXR makes heterodimers with different counterparts accordingly. The suggested existence of a common ligand in cnidarians and starfish is thus intriguing and may reflect the ancestral molecular mechanism for metazoan metamorphosis. This hypothesis will be tested by further investigation of metamorphosis in various taxa, such as sea urchins, annelids and molluscs.
We demonstrated that RA signaling performs as the regulator of metamorphosis process in echinoderm (starfish). Considering RA is also core element of metamorphosis regulation in basal group of metazoan (cnidarian), our results suggest that ancestral function of RA signaling is involved in metamorphosis.
SY, YM, and HW contributed to the design of the experiments. SY and YM collected the samples. SY performed the experiments. SY and MH analyzed the data. All authors wrote the manuscript. All authors read and approved the final manuscript.
We thank Ei-ichi Okumura, Ritsu Kuraishi, and Akihito Wanishi for help in collecting P. pectinifera. We also thank the Ibaraki Prefectural Oarai Aquarium for providing the seawater for culturing adult starfish.
The authors declare that they have no competing interests.
Availability of data and materials
The sequences of starfish gene (raldha, raldhb, raldhc, rar, and rxr) were deposited to the DNA Data Bank of Japan (DDBJ; LC379260, LC379261 LC379262, LC379258 and LC379259, respectively). The dataset for phylogenetic analysis was supplied in Additional files 10 and 11. The raw reads of transcriptome are available from the DDBJ Sequence Reads Archives (DRA006662).
Consent for publication
Ethics approval and consent to participate
This work was supported by JSPS KAKENHI Grant No. 15KT0074.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gutierrez-Mazariegos J, Schubert M, Laudet V. Evolution of retinoic acid receptors and retinoic acid signaling. Subcell Biochem. 2014;70:55–73.View ArticlePubMedGoogle Scholar
- Rhinn M, Dolle P. Retinoic acid signalling during development. Development. 2012;139(5):843–58.View ArticlePubMedGoogle Scholar
- Marlétaz F, Holland LZ, Laudet V, Schubert M. Retinoic acid signaling and the evolution of chordates. Int J Biol Sci. 2006;2(2):38.View ArticlePubMedPubMed CentralGoogle Scholar
- Escriva H, Delaunay F, Laudet V. Ligand binding and nuclear receptor evolution. BioEssays. 2000;22(8):717–27.View ArticlePubMedGoogle Scholar
- Cañestro C, Postlethwait JH, Gonzàlez-Duarte R, Albalat R. Is retinoic acid genetic machinery a chordate innovation? Evol Dev. 2006;8(5):394–406.View ArticlePubMedGoogle Scholar
- Albalat R. The retinoic acid machinery in invertebrates: ancestral elements and vertebrate innovations. Mol Cell Endocrinol. 2009;313(1–2):23–35.View ArticlePubMedGoogle Scholar
- Campo-Paysaa F, Marletaz F, Laudet V, Schubert M. Retinoic acid signaling in development: tissue-specific functions and evolutionary origins. Genesis. 2008;46(11):640–56.View ArticlePubMedGoogle Scholar
- Kuno S, Kawamoto M, Okuyama M, Yasumasu I. Outgrowth of pseudopodial cables induced by all-trans retinoic acid in micromere-derived cells isolated from sea urchin embryos. Dev Growth Differ. 1999;41(2):193–9.View ArticlePubMedGoogle Scholar
- Sciarrino S, Matranga V. Effects of retinoic acid and dimethylsulfoxide on the morphogenesis of the sea urchin embryo. Cell Biol Int. 1995;19(8):675–80.View ArticlePubMedGoogle Scholar
- Fuchs B, Wang W, Graspeuntner S, Li Y, Insua S, Herbst EM, Dirksen P, Bohm AM, Hemmrich G, Sommer F, et al. Regulation of polyp-to-jellyfish transition in Aurelia aurita. Curr Biol. 2014;24(3):263–73.View ArticlePubMedGoogle Scholar
- Heyland A, Moroz LL. Signaling mechanisms underlying metamorphic transitions in animals. Integr Comp Biol. 2006;46(6):743–59.View ArticlePubMedGoogle Scholar
- Laudet V. The origins and evolution of vertebrate metamorphosis. Curr Biol. 2011;21(18):R726–37.View ArticlePubMedGoogle Scholar
- Buszczak M, Segraves WA. Drosophila metamorphosis: the only way is USP? Curr Biol. 1998;8(24):R879–82.View ArticlePubMedGoogle Scholar
- Thummel CS. From embryogenesis to metamorphosis: the regulation and function of Drosophila nuclear receptor superfamily members. Cell. 1995;83(6):871–7.View ArticlePubMedGoogle Scholar
- Buszczak M, Segraves WA. Insect metamorphosis: out with the old, in with the new. Curr Biol. 2000;10(22):R830–3.View ArticlePubMedGoogle Scholar
- Brown DD, Cai L. Amphibian metamorphosis. Dev Biol. 2007;306(1):20–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Bishop CD, Brandhorst BP. Development of nitric oxide synthase-defined neurons in the sea urchin larval ciliary band and evidence for a chemosensory function during metamorphosis. Dev Dyn. 2007;236(6):1535–46.View ArticlePubMedGoogle Scholar
- Bishop CD, Brandhorst BP. NO/cGMP signaling and HSP90 activity represses metamorphosis in the sea urchin Lytechinus pictus. Biol Bull. 2001;201(3):394–404.View ArticlePubMedGoogle Scholar
- Chino Y, Saito M, Yamasu K, Suyemitsu T, Ishihara K. Formation of the adult rudiment of sea urchins is influenced by thyroid hormones. Dev Biol. 1994;161(1):1–11.View ArticlePubMedGoogle Scholar
- Heyland A, Price DA, Bodnarova-Buganova M, Moroz LL. Thyroid hormone metabolism and peroxidase function in two non-chordate animals. J Exp Zool B Mol Dev Evol. 2006;306(6):551–66.View ArticlePubMedGoogle Scholar
- Johnson LG. Stage-dependent thyroxine effects on sea urchin development. N Z J Mar Freshw Res. 1998;32(4):531–6.View ArticleGoogle Scholar
- Miller AE, Heyland A. Iodine accumulation in sea urchin larvae is dependent on peroxide. J Exp Biol. 2013;216(Pt 5):915–26.View ArticlePubMedGoogle Scholar
- Naidenko TK. Induction of metamorphosis of two species of sea urchin from Sea of Japan. Mar Biol. 1996;126(4):685–92.View ArticleGoogle Scholar
- Swanson RL, Williamson JE, De Nys R, Kumar N, Bucknall MP, Steinberg PD. Induction of settlement of larvae of the sea urchin Holopneustes purpurascens by histamine from a host alga. Biol Bull. 2004;206(3):161–72.View ArticlePubMedGoogle Scholar
- Sutherby J, Giardini J-L, Nguyen J, Wessel G, Leguia M, Heyland A. Histamine is a modulator of metamorphic competence in Strongylocentrotus purpuratus (Echinodermata: Echinoidea). BMC Evol Biol. 2012;12(1):14.View ArticleGoogle Scholar
- Heyland A, Hodin J. Heterochronic developmental shift caused by thyroid hormone in larval sand dollars and its implications for phenotypic plasticity and the evolution of nonfeeding development. Evolution. 2004;58(3):524–38.View ArticlePubMedGoogle Scholar
- Heyland A, Reitzel AM, Hodin J. Thyroid hormones determine developmental mode in sand dollars (Echinodermata: Echinoidea). Evol Dev. 2004;6(6):382–92.View ArticlePubMedGoogle Scholar
- Johnson LG, Cartwright CM. Thyroxine-accelerated larval development in the crown-of-thorns starfish, Acanthaster planci. Biol Bull. 1996;190(3):299–301.View ArticlePubMedGoogle Scholar
- Rowley R. Settlement and recruitment of sea urchins (Strongylocentrotus spp.) in a sea-urchin barren ground and a kelp bed: are populations regulated by settlement or post-settlement processes? Mar Biol. 1989;100(4):485–94.View ArticleGoogle Scholar
- Koga H, Matsubara M, Fujitani H, Miyamoto N, Komatsu M, Kiyomoto M, Akasaka K, Wada H. Functional evolution of Ets in echinoderms with focus on the evolution of echinoderm larval skeletons. Dev Genes Evol. 2010;220(3–4):107–15.View ArticlePubMedGoogle Scholar
- Murabe N, Hatoyama H, Komatsu M, Kaneko H, Nakajima Y. Adhesive papillae on the brachiolar arms of brachiolaria larvae in two starfishes, Asterina pectinifera and Asterias amurensis, are sensors for metamorphic inducing factor(s). Dev Growth Differ. 2007;49(8):647–56.View ArticlePubMedGoogle Scholar
- Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017.Google Scholar
- Logan M. Biostatistical design and analysis using R: a practical guide. Chichester: Wiley; 2011.Google Scholar
- Zar JH. Biostatistical analysis. 4th ed. Englewood Cliffs: Prentice Hall; 1998.Google Scholar
- Patel RK, Jain M, Liu Z. NGS QC toolkit: a toolkit for quality control of next generation sequencing data. PLoS ONE. 2012;7(2):e30619.View ArticlePubMedPubMed CentralGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin YZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2017. https://doi.org/10.1093/bib/bbx108.PubMedGoogle Scholar
- Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25(15):1972–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Morino Y, Koga H, Tachibana K, Shoguchi E, Kiyomoto M, Wada H. Heterochronic activation of VEGF signaling and the evolution of the skeleton in echinoderm pluteus larvae. Evol Dev. 2012;14(5):428–36.View ArticlePubMedGoogle Scholar
- Degnan SM, Degnan BM. The initiation of metamorphosis as an ancient polyphenic trait and its role in metazoan life-cycle evolution. Philos Trans R Soc Lond B Biol Sci. 2010;365(1540):641–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Heyland A, Moroz LL. Cross-kingdom hormonal signaling: an insight from thyroid hormone functions in marine larvae. J Exp Biol. 2005;208(Pt 23):4355–61.View ArticlePubMedGoogle Scholar