The evolution of neurosensation provides opportunities and constraints for phenotypic plasticity

Phenotypic plasticity is widely regarded as important for enabling species resilience to environmental change and for species evolution. However, insight into the complex mechanisms by which phenotypic plasticity evolves in nature is limited by our ability to reconstruct evolutionary histories of plasticity. By using part of the molecular mechanism, we were able to trace the evolution of pre-feeding phenotypic plasticity across the class Echinoidea and identify the origin of plasticity at the base of the regular urchins. The neurosensory foundation for plasticity was ancestral within the echinoids. However, coincident development of the plastic trait and the neurosensory system was not achieved until the regular urchins, likely due to pleiotropic effects and linkages between the two colocalized systems. Plasticity continues to evolve within the urchins with numerous instances of losses associated with loss of sensory abilities and neurons, consistent with a cost of maintaining these capabilities. Thus, evidence was found for the neurosensory system providing opportunities and constraints to the evolution of phenotypic plasticity.


Scientific Reports
| (2022) 12:11883 | https://doi.org/10.1038/s41598-022-15583-y www.nature.com/scientificreports/ to remain due to diminished selection pressure. For example, if predator-induced plasticity is lost in a species or line of Daphnia, an artificially-induced expression of juvenile hormones 21,22 may still produce a phenotype that mimics the predator-induced form. We take advantage of knowledge of part of the mechanism for phenotypic plasticity in sea urchin larvae. The feeding structure of many species of sea urchin vary with food concentration throughout larval development e.g. [23][24][25][26][27] , including during the pre-feeding stage. When food is abundant, post-oral arm length is shorter. When food is scarce, post-oral arm length is longer. Plasticity during the pre-feeding stage must be sensory driven, since food is not yet ingested 24,28 . Although the sensory receptor remains unknown, it has been established that sensation of food initiates a dopamine signal which is received by a dopamine type-2 receptor to inhibit post-oral arm elongation 28 . This optimizes arm development and associated feeding potential relative to maternal lipid expenditure 28 . There are distinct phylogenetic limits to when this phenotypic plasticity in arm elongation could have first evolved in echinoderms. While both Echinoidea (urchins and sand dollars) and Ophiuroidea (brittle stars) have a pluteus larval form with skeletal supports, morphological and molecular phylogenies support these as convergent forms that evolved independently [29][30][31] . Thus, it is likely that the plasticity of the pluteus feeding arms (including the skeletal elements), also evolved independently.
Here, we trace the gains and losses of pre-feeding phenotypic plasticity across the echinoids using not only the phenotypic outcome but also part of the underlying developmental signaling mechanism to identify the origin of pre-feeding plasticity. To do this, echinoids will first be surveyed for the pre-feeding response to food to determine when phenotypic plasticity evolved (Genesis of Phenotypic Plasticity). Then, to distinguish between evolutionary losses and convergent gains, species lacking a phenotypic response to food will be used to test whether they still retained a phenotypic response to activation of dopamine type D 2 receptors (DRD2) (Remnant Signaling Mechanism). Finally, to test if neural development constrained the evolution of phenotypic plasticity, we characterized the temporal and spatial development of putative dopaminergic neurons (TH-positive) throughout Echinoidea (Dopaminergic Neural Development).

Results
Genesis of phenotypic plasticity. There was no evidence of shortened post-oral arms in the presence of food in the cidaroid Eucidaris tribuloides (Fig. 1, Table 1). E. tribuloides begins feeding before the post-oral arms have substantially elongated. The canonical plastic response was not observed in any of the irregular species tested, Echinarachnius parma, Dendraster excentricus, Encope michelini, and Leodia sexiesperforata (Fig. 1). Postoral arm lengths were significantly different in the presence of food for the keyhole sand dollars E. michelini (Student's t test, p < 0.001) and L. sexiesperforata (F 2,125 = 15.697, p < 0.001) (Fig. 1, Table 1). However, the response was in the opposite direction of the previously described canonical response. Larvae exposed to food concentrations had significantly longer post-oral arms than those without food. This elongation may be a specific response of the keyhole sand dollars (Mellitidae), although changes in post-oral arm length in response to food concentration followed a similar but non-significant trend for E. parma (F 2,213 = 1.434, p = 0.241, Table 1). Pre-feeding plasticity was not detected for the common sand dollar D. excentricus in the experiments (Fig. 1, F 1,115 = 1.010, p = 0.316). The lack of the canonical plasticity in D. excentricus is consistent with the results for the other irregular urchins. Significant pre-feeding phenotypic responses to food abundance were only detected within the regular urchins. Arbacia punctulata, Lytechinus variegatus variegatus, and Strongylocentrotus purpuratus 28 all had significantly shorter post-oral arm lengths in the presence of high food (Fig. 2, Table 2). Three other taxa tested, Echinometra lucunter, Lytechinus pictus, and Lytechinus variegatus carolinus, did not significantly respond to changes in food concentration (Fig. 2, Table 2). Figure 3 shows that consistent with an ancestral origin within the regular urchins, activation of dopamine type-D 2 receptors with the selective agonist, quinpirole, inhibited post-oral arm elongation in all of the regular urchins tested, including those without the response to food (E. lucunter, L. pictus, and L. variegatus carolinus). With each increasing quinpirole concentration treatment, from 0 to 25 µm to 50 µm (37.5 µm for L. pictus), post-oral arm length decreased.

Remnant signaling mechanism.
Dopaminergic neural development. The member of the most basal group, E. tribuloides, developed TH-positive lateral ganglia near the future post-oral arms before feeding and before arm elongation (Fig. 4A). After feeding starts, TH-positive neurons are also detected in the oral ganglia around the mouth and associated with the stomach. The irregular urchins investigated have altered dopaminergic development ( Fig. 4B-D). Tyrosine hydroxylase first appears after feeding begins in D. excentricus, E. parma, and E. michelini. In both D. excentricus and E. parma, TH-positive neurons are detected in the mouth and gut, but not as lateral ganglia near the post-oral arms. Only within regular urchins does the development of the post-oral arms and TH-positive neurons coincide during the pre-feeding stage. Lateral dopaminergic neurons developed during the prism stage, at approximately the time of arm elongation in A. punctulata, L. pictus, L. variegatus variegatus, L. variegatus carolinus, and S. purpuratus 28 (Fig. 5). At the onset of feeding, TH-positive neurons appear around the mouth as oral ganglia and begin to appear in the stomach. The number of TH-positive neurons associated with lateral ganglia near the post oral arms vary between species. Both A. punctulata and L. variegatus subspp. develop multiple TH-positive neurons along the post-oral arms. Fewer TH-positive neurons develop in L. pictus and the shorter S. purpuratus arms. Echinometra lucunter is the exception-this species does not develop TH-positive neurons until post-feeding and even then, the lateral ganglia appear to be absent. www.nature.com/scientificreports/ Figure 1. Canonical pre-feeding plasticity is absent in the basal Cidaroids and irregular urchins. Change in post-oral arm length at initiation of feeding averaged across families with food concentration in the Cidaroids and Irregularia. Phylogenetic tree is not scaled to divergence. Error bars, ± standard error of the mean. Letters denote a significant difference between food treatments at p < 0.05 (Table 1).

Discussion
Our data suggest that pre-feeding phenotypic plasticity of the post-oral arms arose in the regular urchins and has continued to evolve within the clade. The cidaroid E. tribuloides was selected for the pre-feeding response experiment because molecular and morphological data place cidaroids as the most basal extant taxa within Echinoidea [31][32][33] . The lack of pre-feeding plasticity in E. tribuloides was not unexpected, due to a temporal mismatch between the timing of arm elongation and the onset of feeding. These results support a more recent origin of plasticity within the Echinoids. Irregular urchins have elongated arms during the pre-feeding stage and some are known to alter the length of their post-oral arms in response to food after feeding starts e.g. 23,34 D. excentricus has demonstrable phenotypic plasticity after feeding starts 23,25 and has been previously reported to have the canonical pre-feeding response 24 . The differences between our observations and those of Miner 24 could be due to the different populations tested (Goleta, CA vs Orcas Island, WA) or our ability to detect the small magnitude of change (~ 5% reduction 24 ). However, D. excentricus and S. purpuratus responded morphologically to different cues (soluble vs algal bound, respectively), which is consistent with convergent evolution 24 .
An evolutionary origin of pre-feeding plasticity at the base of the regular urchins is in contrast to phenotypic plasticity that occurs after feeding starts, when additional and more reliable cues, such as metabolic byproducts, could be used to assess food availability. Feeding plasticity has not been reported for any of the basal echinoids tested to date (2 of 2 cidaroids, 3 of 3 diademids). However, both the irregular (3 of 5 species) and regular (9 of 12 species) urchins have taxa that exhibit phenotypic responses to food after feeding starts 23,[25][26][27][35][36][37][38][39][40][41][42][43] . Interestingly, the two species of irregular urchins reported to lack feeding plasticity are the mellitid keyhole sand dollars, E. michelini and L. sexiesperforata 34 , which also lacked canonical pre-feeding plasticity here. Similarly, the three species of regular urchins lacking post-feeding plasticity were species in the genus Echinometra 35 , including E. lucunter, which also lacked canonical pre-feeding plasticity here. This may be a recent loss isolated to the genus Echinometra, as feeding plasticity was reported in the Echinometrid Heliocidaris tuberculata 36 . Knowledge of the mechanism(s) underlying plasticity during the planktonic feeding stage would again provide the ability to discriminate between evolutionary losses and multiple convergent gains. The basal position of A. punctulata within the regular urchins supports the interpretation that plasticity is ancestral within the clade and that there have been multiple losses (Fig. 2). However, the alternative of multiple convergent evolutionary events within the regular urchins is also a possibility based on these data.
Dopamine signaling through DRD2 is required for the presence of food to inhibit arm elongation in the regular urchin S. purpuratus 28 . If plasticity is ancestral within the regular urchins, we would expect that all of the regular urchins would use this same neural signaling mechanism and that even those that lost the plastic response might still retain this signaling remnant. Alternatively, if plasticity evolved convergently multiple times, we might expect differences in the neural signaling mechanism and no response to dopamine signaling in those species without plasticity. Our results suggest that development of the dopaminergic neurons in the lateral ganglia was already in place within the ancestral echinoids, the cidaroids. This provides a foundational component that could have later facilitated the evolution of pre-feeding phenotypic plasticity. However, in the cidaroids, post-oral arm elongation does not occur until days after feeding starts. So, the timing of skeletal elongation may have been an ancestral constraint on the evolution of pre-feeding plasticity.
We propose that pleiotropic effects or gene linkage associated with the temporal shift in arm elongation altered development of dopaminergic neurons in the lateral ganglia (Fig. 6). This would explain the loss or temporal shift in the development of dopaminergic lateral ganglia in the irregular urchins investigated. The lateral ganglia develops within the lateral/boundary ectoderm, where epithelial-mesenchymal signaling is known to coordinate skeletal elongation 30,44-47 . Thus, it is possible that changes in the signaling milieu to advance skeletal elongation could have suppressed dopaminergic development. In support of this, many of the genes within the skeletogenic gene regulatory network 48 , including FGF, Pax 2/5/8, Wnt5, and Otp 44,49,50 , also have roles in dopaminergic development in other systems e.g. [51][52][53] . A decoupling of gene expression or function in the regular urchins would be necessary to allow for the coincident development of dopaminergic lateral ganglia and post-oral arms during the pre-feeding stage. The loss of plasticity in the regular urchin E. lucunter could represent a reversion to the irregular-like state with early skeletal elongation and delayed neural development. Thus, dynamic changes in the development of the lateral ganglia throughout Echinoidea are likely to have both constrained and provided opportunities for plasticity.
Changes in neural development cannot account for the differences in plasticity within all of the regular urchins. L. variegatus carolinus and L. pictus develop TH-positive cells at the appropriate time and place (Fig. 5B,C) and respond to activation of DRD2 (Fig. 3B,C). However, they lack the pre-feeding phenotypic  Pre-feeding plasticity has dynamically evolved within the regular urchins. Change in post-oral arm length at initiation of feeding averaged across families with food concentration in Echinacea, the regular urchins. *Data from a single family of A. punctulata is presented for clarity, though food treatment was significant across all families tested (Table 2). Phylogenetic tree is not scaled to divergence. Green dotted lines denote taxa with canonical pre-feeding plasticity. Error bars, ± standard error of the mean. Letters denote a significant difference between food treatments at p < 0.05 (Table 2). www.nature.com/scientificreports/ response to food concentration (Fig. 2). Since we can detect TH-positive cells in the mouth and gut, we do not believe that the absence of TH-positive lateral ganglia is due to a detection issue. The lack of early dopaminergic lateral ganglia is consistent with the lack of feeding arm plasticity detected within the irregular urchins ( Fig. 1) and suggests that neural development may have constrained the evolution of pre-feeding plasticity within this clade. This indicates that the change responsible for the loss of plasticity occurred upstream of the dopamine receptor-during the neurosensory process. Given that there is evidence for rapid evolution of putative sensory receptors in sea urchins 18 that could affect developmental plasticity, we hypothesize that changes in sensory receptor expression or sequence caused the loss of arm plasticity. Sea urchins have a large repertoire of GPCRs 18 and immune receptors 54,55 that could act in the sensation of food. Both immunity receptors and GPCRs are often found in large tandem arrays of genes and pseudogenes suggestive of gene duplications. In the purple sea urchin, 979 GPCRs have been identified-comprising nearly 3% of the predicted proteins. Two groups of these GPCRs have rapidly expanded and are most similar to vertebrate olfactory receptors 18 . Although innate immunity receptors are generally believed to have ancient origins and minimal subsequent evolution, there is genomic evidence from sea urchins for extensive radiations in this group as well, with 10-20 fold more genes in the purple sea urchin than in humans 54,55 . Rapid evolution of sensory receptors is also consistent with the recent evolutionary loss of the response between relatively close sister species (~ 3 million years) and subspecies (less than a million years) in the genus Lytechinus 56 . However, the identity of the sensory receptor and its evolution remains to be determined. Thus, multiple distinct changes in neural development-timing (heterochrony), number (heterometry), and location (heterotopy)-could have contributed to the evolutionary constraints and opportunities for pre-feeding phenotypic plasticity in sea urchin larvae. To conclude, our data demonstrate the power of a comparative approach to understand the evolutionary dynamics of phenotypic plasticity when part of the molecular mechanism is known. Once within an evolutionary context, we were able to assess the role of neural development in constraining the evolution of plasticity. In this case, ancestral neural development provided a foundational opportunity rather than a constraint. Instead, we propose that interactions between neural development and development of the plastic trait constrained the  www.nature.com/scientificreports/ rise of phenotypic plasticity, and decoupling was necessary to allow for the advent of plasticity. Once established, phenotypic plasticity has continued to evolve dynamically both through changes in neural development and potential evolution of sensory receptors.

Materials and methods
Embryo and larval culture. Adult echinoids were obtained from the following vendors for broodstock: Embryos were cultured using standard methods at densities of 1-5 embryos ml −1 in artificial seawater (ASW) at 21 °C for tropical species or 15 °C for temperate species. Larvae were treated with 5000, 7500 or 10,000 cells ml −1 of the algae Dunaliella sp., dependent on the size of larval spp., to assay for the developmentalresponse to food. Algal concentration was determined using a hemocytometer. Larvae of the regular echinoids were also treated with the specific type-D 2 receptor agonist 57 , quinpirole, at late gastrula stage or prism stage to test for conservation of the dopamine-signaling mechanism. Doses of 0, 25, and 50 µM were used. The highest dose was decreased to 37.5 µM for L. pictus due to sickness in this species at 50 µM.
Quantification of skeletal lengths. Post-oral arm and body rod lengths were assayed just before feeding begins as in Adams et al. 28 . The time post fertilization varied with each species and was experimentally determined by observing algal particles within the gut. All collections were done when algae were observed in less than 50% of the larvae's guts. Larvae were randomly sampled from each treatment, such that sample sizes varied but all were n ≥ 20 individuals each. Larvae were squash mounted on microscope slides to position the skeletal elements in the same plane, then imaged on a Zeiss Axiovert 200 M or Zeiss Axiovert A1 inverted microscope at 20 × under differential inference contrast (DIC) which readily identifies the birefringent skeletal elements. The skeletal lengths were quantified from the digital images using Zen Lite software (Carl Zeiss MicroImaging).

Statistical analyses.
The response of post-oral arm length to algal and quinpirole treatments was assessed using a two-way ANCOVA, where perturbation treatment (food or quinpirole) and biological replicate (malefemale cross) were fixed effects. Body rod length was included as a covariate. We used post-hoc Bonferronicorrected pair-wise comparisons when effects were significant at p < 0.05. Experiments were replicated with two or more sets of non-related full siblings (male-female crosses) for all species except E. michelini and L. sexiesperformata, due to limitations in obtaining ripe broodstock. For these species, only one male and female were available yielding one set of full siblings; thus, an one-way ANOVA (E. michelini) or Student's two-tailed t-tests (L. sexiesperformata) were used. All datasets were determined to be normal based on probability distribution plots. All statistical analyses were done in SYSTAT v10 with output to three decimal places, thus exact p values are given if p > 0.001.   28 on two stages of larvae, (1) just after the initial elongation of the post-oral arms and (2) after feeding started. When tyrosine hydroxylase was not detected at these developmental stages, later stage larvae were also assayed to ensure that the antibody worked in all species tested. Specificity of the antibody in echinoids was established in S. purpuratus by morpholino knock down of tyrosine hydroxylase 28 . Larvae were imaged using a Zeiss Axiovert 200 M epifluorescent inverted microscope with an optically sectioning Zeiss LSM 710 Confocal microscope at 20×. Stacked images were prepared using Imaris (Bitplane Inc., St. Paul, MN).

Data availability
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.