- Research article
- Open Access
The crowns have eyes: multiple opsins found in the eyes of the crown-of-thorns starfish Acanthaster planci
BMC Evolutionary Biology volume 18, Article number: 168 (2018)
Opsins are G protein-coupled receptors used for both visual and non-visual photoreception, and these proteins evolutionarily date back to the base of the bilaterians. In the current sequencing age, phylogenomic analysis has proven to be a powerful tool, facilitating the increase in knowledge about diversity within the opsin subclasses and, so far, at least nine types of opsins have been identified. Within echinoderms, opsins have been studied in Echinoidea and Ophiuroidea, which do not possess proper image forming eyes, but rather widely dispersed dermal photoreceptors. However, most species of Asteroidea, the starfish, possess true eyes and studying them will shed light on the diversity of opsin usage within echinoderms and help resolve the evolutionary history of opsins.
Using high-throughput RNA sequencing, we have sequenced and analyzed the transcriptomes of different Acanthaster planci tissue samples: eyes, radial nerve, tube feet and a mixture of tissues from other organs. At least ten opsins were identified, and eight of them were found significantly differentially expressed in both eyes and radial nerve, with R-opsin being the most highly expressed in the eye.
This study provides new important insight into the involvement of opsins in visual and nonvisual photoreception. Of relevance, we found the first indication of an r-opsin photopigment expressed in a well-developed visual eye in a deuterostome animal. Additionally, we provided tissue specific A. planci transcriptomes that will aid in future Evo Devo studies.
Light carries an immense amount of information about the surroundings. Direct light from the sun, the moon, or the stars is used by various animals to set diurnal or annual clocks and to set direction during navigational tasks. Light reflected from the surroundings guides innumerous different behaviors as it provides information about objects with unprecedented details and speed. Light reception is thus widespread in the animal kingdom but interestingly the molecular machinery behind light reception shares many common features across all phyla. In most cases examined so far the first step in the phototransduction, the absorption of the photons in metazoan, is mediated by a specific protein family called opsins . Opsins are seven transmembrane G protein-coupled receptors binding a chromophore, retinal, which undergoes a conformational change upon the absorption of light, thus triggering the rest of the transduction cascade. The coupling of the opsin and the G protein alpha subunit determines which phototransduction cascade, or particular function will take place [2,3,4,5]. During the last couple of decades, several molecular studies have examined the diversity of the opsin family, identifying at least three major clades (rhabdomeric, ciliary, and RGR/Go or group 4) . As the number of sequenced species grow, so does the number of opsin groups, enabling to better classify of previously considered lineage specific opsins (e.g. cnidopsins). Recently, a classification has been proposed which has nine clades of opsins in total [7, 8]. While Rhabdomeric (r-opsin) and, in some rare cases, go-opsins are the primary visual opsins in protostomes (such as insects, annelids and mollusks) and ciliary opsins (c-opsins) have been identified as the primary visual opsin in vertebrates, little is known about the opsin types used for vision in early branching deuterostomes (such as protochordates and ambulacrarians). What the work has also shown is that many animals have a surprisingly high number of opsin gene copies and that they can be expressed in almost any body region or organ . In many of these cases, the functions remain unknown and may well be outside light reception [10, 11].
Light reception is known from all major groups of echinoderms and is facilitated by different types of photoreceptors ranging from non-pigmented dermal photoreceptors to proper image forming eyes. One of the best studied examples of dermal photoreception is found in the brittlestar genus Ophiocoma [12, 13]. Another dermal photoreceptor system is found in sea urchins which has been suggested to support image forming vision [14, 15]. The genome has been sequenced for the sea urchin Strongylocentrotus purpuratus and eight opsin genes were found belonging to the opsin clades c-opsins, r-opsins, Go-opsins, peropsin, neuropsin and echinopsins A and B [7, 16]. The latter two groups were recently renamed as bathyopsins and chaopsins, respectively . The r-opsin Sp-opsin4 is expressed in cells at the base of the transparent tube feet and is putatively mediating the directional negative phototaxis described for a couple of species [15, 17, 18]. The brittle star Amphiura filiformis has even higher opsin diversity with at least 13 gene copies , but here little is known about expression patterns and behavioural roles.
Both dermal photoreception and proper eyes are known from several species of starfish (Asteroidea) [20,21,22] but it is unknown if it is opsin based and if so whether it is the same opsin in the two systems. The eyes are found in most non-burrowing starfish species at the tip of each arm sitting at the base of the unpaired terminal tube foot as a direct extension of the radial nerve. They are compound eyes and structurally they resemble the eyes of arch clams and fan worms [23, 24] with lens-less ommatidia typically 20–40 μm in diameter (Fig. 1). Depending on species, adult specimens have 50–300 ommatidia in each eye and recent studies have shown that this supports spatial resolution in the range of 8–17 degrees used for navigation [25,26,27]. These studies have also indicated that the ommatidia have a single population of photoreceptors which utilize an opsin with peak sensitivity in the deep blue part of the spectrum around 470 nm.
Recently, two draft genomes of the crown-of-thorns starfish (COTS) Acanthaster planci,Footnote 1 relative to animals collected from Okinawa, Japan and Great Barrier Reef (GBR), Australia, were released. These two genomes shared 98.8% nucleotide identity and were determined to be the same species , which is in agreement with previous classification of the Pacific ocean COTS being one species . Although A. planci is not the first Asteroidea with an assembled genome, it is the first species with well-defined eyes to have an assembled genome. The GBR genome was released along with annotations for ~ 24,500 protein coding genes , and this presents an opportunity to study the mechanisms behind image forming vision in an echinoderm. There are various theories on the origin of opsins and which opsins were present in the ancestral eye . Investigating a wide range of species with various types of opsin based photoreception will likely aid in the understanding of opsin origins. Here we have used tissue specific transcriptomics to investigate the differential expression of opsin genes in A. planci. We found ten different opsins belonging to seven different clades and by comparing expression levels in the eyes, in locomotory tube feet, in the radial nerve, and in a mixture of gonadal, stomach, and epidermal tissue, we are able to infer which opsins are most likely used in vision, in non-visual photoreception, and outside photoreception.
We obtained tissue specific transcriptomic data from the locomotive tube feet, the eyes found on the terminal tube feet, the radial nerve, and a mixture of gonad (unsexed), stomach, and epidermis from A. planci in biological triplicates in order to determine the differential expression of eye related genes, chiefly opsins. Each library produced 25–32 million reads, one exception (one eye library produced 18 million reads) and mapped the A. planci genome at a rate of ~ 63%. Initially we identified all possible opsin genes found in the available A. planci Great Barrier Reef genome  and classified them. We analyzed A. planci’s predicted proteins using a combination of Reciprocal Best Hits (RBH) BLAST against sequences collected from literature and hidden Markov model search with pantherSCORE2.0 , which returned thirteen putative opsin genes that share high sequence similarity. Closer examination using the genome revealed that two of these sequences were actually one fragmented opsin (gbr.65.47.t1 and gbr.65.48.t1), reducing the total putative opsins to twelve. We manually edited the sequence, which can be found in the opsin fasta file of Additional file 1. These sequences were included in a phylogenetic analysis totaling 169 sequences spanning 40 different species (see Additional file 1). Our phylogenetic analysis, from both Mr. Bayes  (Bayesian) and iqtree  (maximum-likelihood), revealed ten A. planci opsins belonging to 7 groups and two sequences clustering with the melatonin receptors outgroup (Fig. 2; Additional file 2: Figure S1): 1 rhabdomeric opsin (r-opsin), 4 ciliary opsins (c-opsins), 1 peropsin, 1 Go-opsin, 1 RGR opsin, 1 neuropsin, and 1 chaopsin. Chaopsins were first classified as echinoderm specific [7, 16] but were later found to group with several cnidarian opsins . A. planci possesses representatives of all previously identified echinoderm opsin groups with the exception of bathyopsin (former echinopsin A), which thus has yet to be identified in any starfish. In all cases, A. planci opsins grouped closest with those of Patiria miniata--an eyeless starfish--followed by Asteria rubens opsins.
Reads were quasi-mapped to the available A. planci transcriptome using Salmon (v0. 82) , and the differential expressed genes were identified using DESeq2  with Wald test and were identified as significantly differentially expressed if they had an adjusted p-value of 0.05. We compared the eyes, tube feet, and radial nerves to the mixture of gonad, stomach, and epidermis. Each tissue was sequenced in biological triplicates. Of 24,409 transcripts 2414 (9.9%), 2764 (11%) and 2719, (11%) were significantly higher expressed (adjusted p-value < 0.05) in the tube feet, eyes and radial nerve compared to the mixture. The number of transcripts that had lower expression than in the mixed tissue, was 2839 (12%), 3598 (15%) and 3490 (14%) in the tube feet, eyes and radial nerve respectively. In A. planci eyes, all opsins with the exception of ciliary opsin 1.1b and neuropsin showed higher expression in comparison to the mixed tissues (Fig. 3). This was also the case when comparing opsin expression in the radial nerve to the mixed sample but to a lesser degree (Additional file 3: Figure S2a). Expression of opsins in the tube feet was similar to the mixed tissue for 6 of the opsins but for c-opsin 1.1a, c-opsin 1.1b, c-opsin 1.3 and neuropsin. There was significantly lower expression (adjusted p-value < 0.05) in the tube feet (Additional file 3: Figure S2b). Additionally, we preformed qPCR on several opsins and the neural gene synaptotagmin 1 (syt1) to support out RNA-seq findings (Additional file 4: Figure S3).
In order to assess the putative functionality of the A. planci identified opsin sequences, an analysis of the key residues necessary for opsin function was performed (Table 1). In most opsins, the retinal binds to the Lysine K296 via a Schiff base bond, however the proton in the opsin protein is unstable and a counter ion is needed and often supplied by the highly conserved Glutamic acid (E113) residue. There are cases, however, where this residue is replaced by a Tyrosine (Y), Phenylalanine (F), Methionine (M), or Histidine (H) and the other highly conserved Glutamic acid residue, E181, serves as the counter ion . This is the case with the majority of the opsins in A. planci, where E113 are replaced with a Tyrosine, with the exception of Ap-Go-opsin which has an Isoleucine (I) in the position 113.
In addition to ten opsin sequences, we have also observed ten A. planci sequences that are potential G protein alpha subunits. Phylogenomic methods classified these sequences as 3 Gαs, 1 Gαo, 4 Gαi, 1 Gαq, and 1 Gα12 (Additional file 5: Figure S4). All identified G protein alpha subunits with the exception of 1 Gαs (gbr.231.19.t1), 1 Gαi (gbr.143.10.t1) and the Gα12 are higher expressed in the eyes of A. planci compared to the mixed tissue samples (Additional file 6: Fgure S5).
Ciliary and rhabdomeric opsins
There are four ciliary opsins identified in the A. planci genome, three of which are expressed in the eyes; c-opsin 1.1a, 1.2 and c-opsin 1.3. These three c-opsins were closely clustered on gbr_scaffold65 (~ 70 kb), and observed to be significantly differentially expressed in the animal, as neither are expressed in the mixed tissue (Fig. 4; Additional file 4: Figure S3). Of the Ap-c-opsins, 1.1a is the highest expressed, followed by Ap-c-opsin 1.2 and 1.3. C-opsin 1.2 was observed to not have the K296 residue, which is required for the formation of the Schiff base, but instead an Arginine (R) residue is found in this position. This is the only A. planci opsin missing this key residue. While the E113 counter ion is replaced with Tyrosine (Y) in all of the Ap-c-opsins, the E181 counter ion is present in Ap-c-opsin1.1a, completely missing in Ap-c-opsin 1.1b and 1.2 and replaced with an Alanine (A) in Ap-c-opsin 1.3. The Ap-c-opsin 1.2 and Ap-c-opsin 1.3 are missing other important motifs, the C187 and C110 disulfide bond motifs , respectively (Table 1).
A. planci’s r-opsin, on the other hand, is the most highly differentially expressed of the opsins and of any other genes when comparing eyes with the mixed tissue (Fig. 3). Further, the Ap-r-opsin was observed to be the most highly expressed in the starfish eye and its sequence features the Lysine residue (K296), critical for the Schiff base formation, and a putative counter ion (E181) (Table 1).
Chaopsin is a recently identified group of opsins. Ramirez et al.  found the formerly described groups of anthozoa I opsins  and the echinoderm echinopsin B  to cluster forming the chaopsin group. In agreement with D’Aniello et al. (2015) we found an A. planci’s chaopsin to cluster together with other ambulacrarian chaopsins between the r-opsin and c-opsin clades (Fig. 2). Ap-chaopsin is amongst the highest differentially expressed opsins in our A. planci transcriptomes, with ~ 9.7 log2 fold changes in the eye compared to the mixed tissue (Fig. 3). It is also expressed in the radial nerve, but to a far lesser degree, and are not significantly expressed in the mixed tissues or the tube feet (Fig. 4 and Additional file 3: Figure S2).
Peropsin RGR and go opsins
Peropsin and RGR opsin are the highest expressed of the opsins in the tube feet, mixed, and radial nerve, although both still have higher expression in the eyes (Fig. 4 and Additional file 3: Figure S2). The disulfide bond linkage C110/C187, counterion sites and the Lysine for the Schiff base formation are all present. However, both A. planci RGR-opsin and peropsin contain variations of the NPxxY motif, NAALQ and NPLMF, respectively (Table 1). Additionally, peropsin has a variation of the LxxxD motif, ASAGD. Ap-Go-opsin was expressed in both the eyes and radial nerve but was not found in the mixed sample or in the tube feet (Figs. 3 and 4).
Light sensing is an important aspect of life and much of it is mediated or initiated by the G-protein-coupled receptor proteins, opsins. The release of the annotated draft genome of A. planci has prompted us to investigate its opsin repertoire and expression in a tissue specific manner. This allowed us to classify the specific opsins and to infer possible function and further expand the knowledge of opsin evolution especially within deuterostomes. Ten opsins were identified spanning seven clades: r-opsin, c-opsin, Go-opsin, peropsin, neuropsin, RGR opsin and chaopsin. Opsins have also been sequenced from two other starfish species, Asterias rubens and Patiria miniata. Through a phylogenomic analysis, it was observed that A. planci opsins grouped closest to those in the eyeless P. miniata. This grouping is in accordance with the phylogenetic position of these starfish species , with A. planci as an Acanthasteridae more closely related to P. miniata, an Asterinidae, (both species belonging to Valvatida)- than to A. rubens belonging to Forcipulatacea. However, studies on tissue specific opsin expression are needed to reveal in which organs of the eyeless starfish P. miniata the respective opsin orthologs are expressed. To this point, it remains unclear if opsins potentially serving a visual function in the eye possessing A. planci might have switched functions in the eyeless representative. Alternatively, opsin expression in P. miniata might simply persist as a potential evolutionary remain, a finding known e.g. from blind cave salamanders. Their photoreceptors remain to express opsins despite the apparent loss of shading pigmentation and their highly degenerated morphology .
Electrophysiological recordings from the photoreceptors of A. planci have strongly suggested light absorption here to be utilizing a single opsin. Of the several opsins found to be expressed in the eyes of A. planci, r-opsin was the highest and most differentially expressed gene in the eye transcriptome when compared to mixed samples, thus suggesting this to be the opsin utilized for vision. Whereas rhabdomeric opsins have been described in many protostome species as the primary opsin for vision (reviewed in ) and as a non-visual opsin (melanopsin) in vertebrates (reviewed in ), no deuterostome eye has previously been found to express an r-opsin for primary vision. In the sea urchin S. purpuratus r-opsin in tube feet has been proposed to be involved in a visual context facilitating phototaxis but not in spatial vision , A. planci eyes have been previously demonstrated to perform proper spatial vision [25,26,27]. Our findings thus provide first indications for a deuterostome eye utilizing an r-opsin for spatial vision and might facilitate in identifying the split between r-opsin and c-opsin in complex eyes. However, those findings await final confirmation by future expression studies.
Ciliary opsins are well known for their role in vertebrate vision and brain function in some invertebrates . They are differentially expressed in A. planci eyes, with the exception of Ap-c-opsin 1.1b, which is not expressed in the eye and does not cluster on the same scaffold as the other 3 Ap-c-opsins. Interestingly, the three clustered opsins show a correspondence between their spatial ordering within the cluster (in the direction 3′ to 5′) and their expression level, with most 3′ Ap-c-opsin being the most highly expressed in both eye and radial nerve (Fig. 5). C-opsins appear clustered also in vertebrates, for example, human’s medium and long wavelength ciliary opsins (opsin1MW1, opsin1MW2, opsin1MW3 and opsin1LW) are clustered on chromosome X, while human’s short wavelength c-opsin is on chromosome 7 (Ensembl genome browser; ). Clustering of opsins has so far not been described outside deuterostomes, but it is also seen in the echinoderm, S. purpuratus where Go-opsins (opsin 3.1 and 3.2) clustered on scaffold26 (Echinobase; ). The significance of the occurrence of opsin clustering in echinoderms is not obvious and awaits more functional studies.
Chaopsin (Ap-chaopsin) is the second most differentially expressed opsin in A. planci eyes and it has many of the motifs necessary for phototransduction, including the NPxxY binding motif in the 7th transmembrane domain involved in coupling with the G protein. Little is known about the functions of chaopsins (or opsin5) but they have been identified in Echinoidea, Asteroidea, and Ophiuroidea [7, 19, 45,46,47]. Chaopsins have been found analyzing genomic data or transcriptomes of hypothesized photosensitive tissues, such as the tube feet in Strongylocentrotus droebachiensis  and Strongylocentrotus intermedius . A chaopsin was identified in the genome of the eyeless P. miniata but was not found in transcriptome data of the eye possessing starfish A. rubens. The lack of chaopsin in A. rubens is probably a methodological artifact, since general arm tissue, including radial nerve and tube feet, but not eyes, has been used to generate these transcriptomic data . In the Caribbean elkhorn coral, Acropora palmata, opsin3, which together with echinoderm opsin5 belongs to the chaopsin clade , has been demonstrated to couple with a Gq-protein in a light-dependent manner . This leads us to hypothesize that the Ap-chaopsin may be important for phototransduction in A. planci and potentially in all echinoderms. The exact role of this opsin remains elusive, though.
Peropsin and RGR opsin
Peropsin and RGR opsin were expressed in the tube feet, the mixed tissue, and the radial nerve of A. planci, but both have the highest expression in the eyes. They are considered as photoisomerase enzymes and not as photopigments, since they bind to all-trans retinaldehyde to regenerate 11-cis-retinoids for pigment regeneration. This has been observed for peropsin in the vertebrate retinal pigment epithelium , in cephalopod photoreceptors  and in the jumping spider . Knock-down mice , together with biochemical and spectroscopic studies in amphioxus , have demonstrated the same properties for RGR opsin in chordates. This family of opsins is thus important for visual pigment regeneration . While RGR opsins are known to not have the NPxxY binding motif in the 7th transmembrane domain involved in coupling with the G protein, this motif is often found in peropsins [51, 53]. This was not the case in A. planci, where both RGR opsin and peropsin have varying NPxxY binding motifs. This could alter peropsin’s ability to couple with G proteins, further supporting its function as photoisomerase. However, it is worth mentioning that, in chicken, the presence of both peropsin and RGR opsin is thought to serve in the visual cycle of the circadian clock .
In A. planci eye and radial nerve transcriptomes, we found expression of a Go-opsin along with putative Go alpha subunit proteins. Evidence for Go-opsins in animals is rare. These opsins interact with a specific G-protein that differs from those in the c-opsin as well as the r-opsin transduction cascade [6, 56]. In the retina of the bivalve Patinopecten yessonensis, Go-opsin is expressed in a layer of ciliary photoreceptor cells, which do not express a ciliary opsin . However, in the marine annelid Platynereis dumerilii, Go-opsin is co-expressed with two r-opsins in the photoreceptor cells of the larval eye . Knocking down Go-opsin did not lead to absence of phototaxis in P. dumerilii but did reduce the sensitivity to the blue-cyan part of the color spectrum (λmax = 488 nm). A similar absorbance spectrum (λmax = 483 nm) was observed in amphioxus, Branchinostoma belcheri, Go-opsin [53, 59]. While the morphological and expression data on Go-opsin in P. yessonensis point towards a functioning in a visual context, no other studies have found this. Our data does not reveal if A. planci Go-opsin is co-expressed with any other opsins inside the same cells, nonetheless, the fact that this opsin is expressed in the starfish eye along with the results from annelids and amphioxus, opens up the possibility that it is involved in spectral tuning of vision in A. planci. Such a tuning would indeed be consisted with the spectral sensitivity curves obtained by electrophysiology in this starfish species .
In conclusion, our findings demonstrate that the eye of the starfish A. planci expresses an entire set of ten different opsin proteins. This starfish has recently been demonstrated to possess spatial vision utilized when locating its preferred coral food. Ap-r-opsin is by far the most highly expressed photopigment in its eyes, and thus likely to be the photopigment facilitating these sophisticated visually guided behaviours. A. planci is thus one more echinoderm possessing a much more complex photoreceptor system than previously assumed and the variety of opsins found differentially expressed in various starfish tissues by our transcriptomic analyses sets the ground for comparative studies on evolutionary changes in photoreceptor function that occurred towards the vertebrate eye. In addition to investigating vision in A. planci this study thus provides tissue specific transcriptome data that will aid in future evolutionary studies.
Material and methods
The specimens of A. planci used in this study were hand collected on the Great Barrier Reef off the coast of Cairns, Australia. After collection the animals were kept in holding tanks with running seawater at 26 degrees for 2–3 days and then flown to Denmark. In Denmark they were kept under similar conditions at the Danish National Aquarium, The Blue Planet, where they were fed three times a week with a past of enriched squid and fish meat. Tissue samples were taken from four specimens with diameters of 15–23 cm. Three eyes which were dissected from the terminal tube feet, 3 eye-less locomotory tube feet, approx. 5 cm radial nerve, and pieces of the gonads, the stomachs, and the epidermis were sampled from each of the four specimens and stored in RNAlater at 4 °C. Even though microdissection was performed the eyes will likely contain some remains of the modified tube foot it is situated on. The samples of the radial nerve will also contain the epithelium covering the nerve ventrally. Two additional animals were collected at the coastal reefs of Guam and a total of 12 eyes and 10 cm radial nerve were taken from them directly after collection and stored in RNAlater at 4 °C.
RNA extraction and sequencing
The tissue samples were removed from the RNAlater, immediately frozen with liquid nitrogen and homogenized using a mortar and pestle. Powdered tissues were then dissolved in EuroGOLD RNAPure (EMR 506100) and processed using EUROzol RNA extraction protocol (EMR055100, euroclone), then subjected to LitCl (4 M) purification. Library preps and sequencing were done at Università degli Studi di Salerno in biological triplicates using SMART-Seq v4 Ultra Low Input RNA Kit to ensure enough RNA was available for each library. 10 ng of RNA was used for each library and we used 8 PCR cycles to amplify the samples.
Sequenced reads were examined using fastqc and then quality filtered and trimmed using trimmomatic (v0.33) . Quality controlled reads were quasi-mapped and quantified to v1 Great Barrier Reef Acanthaster planci transcriptome using salmon (v0. 82) . Transcripts per million (TPM), the normalized transcript counts . Differentially expressed genes were identified using DESeq2 . Wald test and transcripts were identified as significantly differentially expressed if they had an adjusted p-value ≤0.05, and − 1.5 ≥ log2FC ≥ 1.5. All scripts can be found at https://github.com/elijahlowe/Acanthaster_opsins.git. Additionally, reads were mapped onto the A. planci genome to determine mapping rates using transrate  default setting.
Identification of opsin sequences
Opsin protein sequences were collected from echinoderms [7, 19, 46, 47, 63,64,65,66], hemichordates , molluscs [56, 68, 69], arthropods [70, 71], vertebrates [55, 72,73,74,75,76,77] and annelids  covering 40 species and 159 opsin sequences. Additional 7 melatonin receptor sequences were used as outgroup. These sequences were retrieved from various databases including Echinobase , and NCBI , as well as from publications themselves, as described in Additional file 1. These sequences were used to perform Reciprocal Best Hits (RBH) BLAST against the A. planci (GBR) predicted proteins. Additionally, pantherSCORE2.0  was used against the A. planci predicted proteins to identify opsin sequences using hidden Markov models. The collected sequences were aligned with MAFFT (v7.215) [80, 81] using L-INS-i algorithm which is better designed for divergence sequences and performed well when benchmarked against other multiple sequence aligners . The aligned sequences were then trimmed using trimAl , removing gaps that occurred in 10% of the alignments while being sure to retain 60% of the total sequence length. Maximum likelihood phylogenetic trees were generated using the aligned sequences with iqtree  with 10,000 ultrafast bootstrap support  using the ‘a Bayesian-like transformation of aLRT’ (abayes) method , and GTR + G amino acid substitution model and visualised with figtree (v1.4.3) . Additional trees were then generated using bayesian tree was then generated using selected using MrBayes (v3.2.5)  50 million generations, with the GTR + G amino acid substitution model as well. Modifications such as additional labels and visual effects were done using inkscape.
First-strand cDNA was synthesized in a 20 μl reaction from 100 ng of total RNA using the SuperScript VILO cDNA synthesis kit (Invitrogen). The synthesis of the cDNA was performed with the following program: 25 °C for 10 min, 42 °C for 90 min, 85°C for 5 min. Specific primer sets for each gene were designed to amplify products 100–200 bp in length using the Primer3 program (Rozen and Skaletsky, 2000). Primers’ sequences were as follows: syt1forward, cgacccctacgtcaaagtgt; syt1 reverse, gacctcaccaatctggtcgt; Ef1a forward, ggtcattggccacgtagact; Ef1a reverse, gctctgccttcaacttgtcc; c-opsin forward, gctatcctggcgctgtactc; c-opsin reverse, tcagtgtgctcggtaggatg; r-opsin forward, cagatcgccaaagtgggtat; r-opsin reverse, gcggaactcttggctaacac; chaopsin forward, cttcctgtcagcctggactc; chaopsin reverse, aactgtgggctcatcaatcc.
Real time QPCR
Real time QPCR amplification was performed with diluted cDNA (1:4 dilution). Each 10 μl reaction contained 5 μl of SYBR Green reagent, 4 μl of forward and reverse primer mix (0,75 μM each) and 1 μl of 1:4 cDNA. Reactions were performed in technical triplicate with three different biological samples using the ViiA 7 REAL TIME PCR with SYBR Green chemistry (Applied Biosystems). The cycling condition was: 95 °C for 20 s, 40 cycles at 95 °C for 1 s and 60 °C for 20 s, 95 °C for 15 s, 60 °C 1 min, and a gradient from 60 °C to 95 °C for 15 min. All primer pairs were validated by QPCR against a positive (genomic DNA) and negative (water) control.
COTS is thought to be divided into four separate species distinguished by bodies of water --Red Sea, the Pacific, the Northern and the Southern Indian Ocean, reclassifying the Pacific ocean species as Acanthaster solaris . However, for consistency with the released genome we will refer to the GBR species as A. planci.
- A rubens:
Counts per million
Great Barrier Reef
- P. miniata :
Reciprocal Best Hits
Transcripts per million
Land MF, Nilsson D-E. Light and vision. In: Animal Eyes. New York: Oxford University Press; 2012. p. 23–45.
Strathmann M, Simon MI. G protein diversity: a distinct class of alpha subunits is present in vertebrates and invertebrates. Proc Natl Acad Sci U S A. 1990;87:9113–7.
Hildebrandt JD. Role of subunit diversity in signaling by heterotrimeric G proteins. Biochem Pharmacol. 1997;54:325–39.
Koyanagi M, Takano K, Tsukamoto H, Ohtsu K, Tokunaga F, Terakita A. Jellyfish vision starts with cAMP signaling mediated by opsin-G(s) cascade. Proc Natl Acad Sci U S A. 2008;105:15576–80.
Mason B, Schmale M, Gibbs P, Miller MW, Wang Q, Levay K, et al. Evidence for multiple phototransduction pathways in a reef-building coral. PLoS One. 2012;7:e50371.
Porter ML, Blasic JR, Bok MJ, Cameron EG, Pringle T, Cronin TW, et al. Shedding new light on opsin evolution. Proc Biol Sci. 2012;279:3–14.
D’Aniello S, Delroisse J, Valero-Gracia A, Lowe EK, Byrne M, Cannon JT, et al. Opsin evolution in the Ambulacraria. Mar Genomics. 2015;24(Pt 2):177–83.
Ramirez MD, Pairett AN, Pankey MS, Serb JM, Speiser DI, Swafford AJ, et al. The Last Common Ancestor of Most Bilaterian Animals Possessed at Least Nine Opsins. Genome Biol Evol. 2016;8:3640–52.
Suga H, Schmid V, Gehring WJ. Evolution and functional diversity of jellyfish opsins. Curr Biol. 2008;18:51–5.
Shen WL, Kwon Y, Adegbola AA, Luo J, Chess A, Montell C. Function of rhodopsin in temperature discrimination in Drosophila. Science. 2011;331:1333–6.
Leung NY, Montell C. Unconventional Roles of Opsins. Annu Rev Cell Dev Biol. 2017; Available from. https://doi.org/10.1146/annurev-cellbio-100616-060432.
Aizenberg J, Tkachenko A, Weiner S, Addadi L, Hendler G. Calcitic microlenses as part of the photoreceptor system in brittlestars. Nature. 2001;412:819–22.
Sumner-Rooney L, Rahman IA, Sigwart JD, Ullrich-Lüter E. Whole-body photoreceptor networks are independent of “lenses” in brittle stars. Proc Biol Sci. 2018;285. Available from. https://doi.org/10.1098/rspb.2017.2590.
Ullrich-Lüter EM, Dupont S, Arboleda E, Hausen H, Arnone MI. Unique system of photoreceptors in sea urchin tube feet. Proc Natl Acad Sci U S A. 2011;108:8367–72.
Ullrich-Lüter EM, D’Aniello S, Arnone MI. C-opsin expressing photoreceptors in echinoderms. Integr Comp Biol. 2013;53:27–38.
Raible F, Tessmar-Raible K, Arboleda E, Kaller T, Bork P, Arendt D, et al. Opsins and clusters of sensory G-protein-coupled receptors in the sea urchin genome. Dev Biol. 2006;300:461–75.
Blevins E, Johnsen S. Spatial vision in the echinoid genus Echinometra. J Exp Biol. 2004;207:4249–53.
Yerramilli D, Johnsen S. Spatial vision in the purple sea urchin Strongylocentrotus purpuratus (Echinoidea). J Exp Biol. 2009;213:249–55.
Delroisse J, Ullrich-Lüter E, Ortega-Martinez O, Dupont S, Arnone M-I, Mallefet J, et al. High opsin diversity in a non-visual infaunal brittle star. BMC Genomics. 2014;15:1035.
Yoshida M, Ohtsuki H. The phototactic behavior of the starfish, Asterias amurensis Lutken. Biol Bull. Marine Biological Laboratory. 1968;134:516–32.
Yoshida M. Extraocular Photoreception. Comparative Physiology and Evolution of Vision in Invertebrates. Berlin, Heidelberg: Springer; 1979. p. 581–640.
Yoshida M, Takasu N, Tamotsu S. Photoreception in Echinoderms. In: Photoreception and Vision in Invertebrates. Boston: Springer; 1984. p. 743–71.
Nilsson D-E. Eyes as Optical Alarm Systems in Fan Worms and Ark Clams. Philos Trans R Soc Lond B Biol Sci. 1994;346:195–212.
Bok MJ, Capa M, Nilsson D-E. Here, There and Everywhere: The Radiolar Eyes of Fan Worms (Annelida, Sabellidae). Integr Comp Biol. 2016;56:784–95.
Garm A, Nilsson D-E. Visual navigation in starfish: first evidence for the use of vision and eyes in starfish. Proc Biol Sci. 2014;281:20133011.
Petie R, Garm A, Hall MR. Crown-of-thorns starfish have true image forming vision. Front Zool. 2016;13:41.
Petie R, Hall MR, Hyldahl M, Garm A. Visual orientation by the crown-of-thorns starfish (Acanthaster planci). Coral Reefs. 2016;35:1139–50.
Vogler C, Benzie J, Lessios H, Barber P, Wörheide G. A threat to coral reefs multiplied? Four species of crown-of-thorns starfish. Biol Lett. 2008;4:696–9.
Haszprunar G, Vogler C, Wörheide G. Persistent Gaps of Knowledge for Naming and Distinguishing Multiple Species of Crown-of-Thorns-Seastar in the Acanthaster planci Species Complex. Diversity. Multidisciplinary Digital Publishing Institute. 2017;9:22.
Hall MR, Kocot KM, Baughman KW, Fernandez-Valverde SL, Gauthier MEA, Hatleberg WL, et al. The crown-of-thorns starfish genome as a guide for biocontrol of this coral reef pest. Nature. 2017;544:231–4.
Arendt D. Evolution of eyes and photoreceptor cell types. Int J Dev Biol. ijdb.ehu.es. 2003;47:563–71.
Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13:2129–41.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. academic.oup.com. 2012;61:539–42.
Minh BQ, Nguyen MAT, von Haeseler A. Ultrafast approximation for phylogenetic bootstrap. Mol Biol Evol. 2013;30:1188–95.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417–9.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Terakita A. The opsins. Genome Biol. 2005;6:213.
Karnik SS, Khorana HG. Assembly of functional rhodopsin requires a disulfide bond between cysteine residues 110 and 187. J Biol Chem. 1990;265:17520–4.
Hering L, Mayer G. Analysis of the opsin repertoire in the tardigrade Hypsibius dujardini provides insights into the evolution of opsin genes in panarthropoda. Genome Biol Evol. 2014;6:2380–91.
Mah C, Foltz D. Molecular phylogeny of the Valvatacea (Asteroidea: Echinodermata). Zool J Linn Soc. Blackwell Publishing Ltd. 2011;161:769–88.
Kos M, Bulog B, Szél A, Röhlich P. Immunocytochemical demonstration of visual pigments in the degenerate retinal and pineal photoreceptors of the blind cave salamander (Proteus anguinus). Cell Tissue Res. 2001;303:15–25.
Plachetzki DC, Serb JM, Oakley TH. New insights into the evolutionary history of photoreceptor cells. Trends Ecol Evol. 2005;20:465–7.
Arendt D, Tessmar-Raible K, Snyman H, Dorresteijn AW, Wittbrodt J. Ciliary photoreceptors with a vertebrate-type opsin in an invertebrate brain. Science. 2004;306:869–71.
Aken BL, Achuthan P, Akanni W, Amode MR, Bernsdorff F, Bhai J, et al. Ensembl 2017. Nucleic Acids Res. 2017;45:D635–42.
Cameron RA, Samanta M, Yuan A, He D, Davidson E. SpBase: the sea urchin genome database and web site. Nucleic Acids Res. 2009;37:D750–4.
Lesser MP, Carleton KL, Böttger SA, Barry TM, Walker CW. Sea urchin tube feet are photosensory organs that express a rhabdomeric-like opsin and PAX6. Proc Biol Sci. 2011;278:3371–9.
Zhao C, Ji N, Tian X, Feng W, Sun P, Wei J, et al. Opsin4, Opsin5, and Pax6 significantly increase their expression in recently settled juveniles of the sea urchinStrongylocentrotus intermedius(Echinodermata: Echinoidea). Invertebr Reprod Dev. 2015;59:119–23.
Delroisse J, Lanterbecq D, Eeckhaut I, Mallefet J, Flammang P. Opsin detection in the sea urchin Paracentrotus lividus and the sea star Asterias rubens. Cah Biol Mar. 2013;54:721–7.
Hao W, Fong HK. The endogenous chromophore of retinal G protein-coupled receptor opsin from the pigment epithelium. J Biol Chem. 1999;274:6085–90.
Hara T, Hara R. Regeneration of squid retinochrome. Nature. 1968;219:450–4.
Nagata T, Koyanagi M, Tsukamoto H, Terakita A. Identification and characterization of a protostome homologue of peropsin from a jumping spider. J Comp Physiol A Neuroethol Sens Neural Behav Physiol. Springer. 2010;196:51–9.
Chen P, Hao W, Rife L, Wang XP, Shen D, Chen J, et al. A photic visual cycle of rhodopsin regeneration is dependent on Rgr. Nat Genet. 2001;28:256–60.
Koyanagi M, Terakita A, Kubokawa K, Shichida Y. Amphioxus homologs of Go-coupled rhodopsin and peropsin having 11-cis- and all-trans-retinals as their chromophores. FEBS Lett. 2002;531:525–8.
Dowling JE. Chemistry of visual adaptation in the rat. Nature. 1960;188:114–8.
Bailey MJ, Cassone VM. Opsin photoisomerases in the chick retina and pineal gland: characterization, localization, and circadian regulation. Invest Ophthalmol Vis Sci. 2004;45:769–75.
Kojima D, Terakita A, Ishikawa T, Tsukahara Y, Maeda A, Shichida YA. Novel Go-mediated Phototransduction Cascade in Scallop Visual Cells. J Biol Chem. 1997;272:22979–82.
Arenas O, Osorno T, Malagón G, Pulido C, Gomez MDP, Nasi E. Molecular and functional identification of a novel photopigment in ciliary photoreceptors. J Gen Physiol. 2018; Available from. https://doi.org/10.1085/jgp.201711938.
Gühmann M, Jia H, Randel N, Verasztó C, Bezares-Calderón LA, Michiels NK, et al. Spectral Tuning of Phototaxis by a Go-Opsin in the Rhabdomeric Eyes of Platynereis. Curr Biol. 2015;25:2265–71.
Tsukamoto H, Terakita A, Shichida Y. A rhodopsin exhibiting binding ability to agonist all-trans-retinal. Proc Natl Acad Sci U S A. 2005;102:6303–8.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010;26:493–500.
Smith-Unna RD, Boursnell C, Patro R, Hibberd JM, Kelly S. TransRate: reference free quality assessment of de-novo transcriptome assemblies [Internet]. 2015. Available from: https://doi.org/10.1101/021626.
Ooka S, Katow T, Yaguchi S, Yaguchi J, Katow H. Spatiotemporal expression pattern of an encephalopsin orthologue of the sea urchin Hemicentrotus pulcherrimus during early development, and its potential role in larval vertical migration. Dev Growth Differ. 2010;52:195–207.
Cannon JT, Kocot KM, Waits DS, Weese DA, Swalla BJ, Santos SR, et al. Phylogenomic resolution of the hemichordate and echinoderm clade. Curr Biol. 2014;24:2827–32.
Mashanov VS, Zueva OR, García-Arrarás JE. Transcriptomic changes during regeneration of the central nervous system in an echinoderm. BMC Genomics. 2014;15:357.
Elphick MR, Semmens DC, Blowes LM, Levine J, Lowe CJ, Arnone MI, et al. Reconstructing SALMFamide Neuropeptide Precursor Evolution in the Phylum Echinodermata: Ophiuroid and Crinoid Sequence Data Provide New Insights. Front Endocrinol. 2015;6:2.
Freeman RM Jr, Wu M, Cordonnier-Pratt M-M, Pratt LH, Gruber CE, Smith M, et al. cDNA sequences for transcription factors and signaling proteins of the hemichordate Saccoglossus kowalevskii: efficacy of the expressed sequence tag (EST) approach for evolutionary and developmental studies of a new organism. Biol Bull. 2008;214:284–302.
Hall MD, Hoon MA, Ryba NJ, Pottinger JD, Keen JN, Saibil HR, et al. Molecular cloning and primary structure of squid (Loligo forbesi) rhodopsin, a phospholipase C-directed G-protein-linked receptor. Biochem J. 1991;274(Pt 1):35–40.
Hara-Nishimura I, Kondo M, Nishimura M, Hara R, Hara T. Cloning and nucleotide sequence of cDNA for rhodopsin of the squid Todarodes pacificus. FEBS Lett. 1993;317:5–11.
Carulli JP, Hartl DL. Variable rates of evolution among Drosophila opsin genes. Genetics. 1992;132:193–204.
Chou WH, Hall KJ, Wilson DB, Wideman CL, Townson SM, Chadwell LV, et al. Identification of a novel Drosophila opsin reveals specific patterning of the R7 and R8 photoreceptor cells. Neuron. 1996;17:1101–15.
Okano T, Kojima D, Fukada Y, Shichida Y, Yoshizawa T. Primary structures of chicken cone visual pigments: vertebrate rhodopsins have evolved out of cone visual pigments. Proc Natl Acad Sci U S A. 1992;89:5932–6.
Sun H, Gilbert DJ, Copeland NG, Jenkins NA, Nathans J. Peropsin, a novel visual pigment-like protein located in the apical microvilli of the retinal pigment epithelium. Proc Natl Acad Sci U S A. 1997;94:9893–8.
Tao L, Shen D, Pandey S, Hao W, Rich KA, Fong HK. Structure and developmental expression of the mouse RGR opsin gene. Mol Vis. 1998;4:25.
Provencio I, Rodriguez IR, Jiang G, Hayes WP, Moreira EF, Rollag MD. A novel human opsin in the inner retina. J Neurosci Soc Neuroscience. 2000;20:600–5.
Kasper G, Taudien S, Staub E, Mennerich D, Rieder M, Hinzmann B, et al. Different structural organization of the encephalopsin gene in man and mouse. Gene. 2002;295:27–32.
Sakami S, Kolesnikov AV, Kefalov VJ, Palczewski K. P23H opsin knock-in mice reveal a novel step in retinal rod disc morphogenesis. Hum Mol Genet. 2014;23:1723–41.
Arendt D, Tessmar K, de Campos-Baptista M-IM, Dorresteijn A, Wittbrodt J. Development of pigment-cup eyes in the polychaete Platynereis dumerilii and evolutionary conservation of larval eyes in Bilateria. Development. 2002;129:1143–54.
NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2016;44:D7–19.
Katoh K. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.
Katoh K, Standley DM. MAFFT: Iterative Refinement and Additional Methods. Methods in Molecular Biology. 2013:131–46.
Thompson JD, Linard B, Lecompte O, Poch O. A comprehensive benchmark study of multiple sequence alignment methods: current challenges and future perspectives. PLoS One. 2011;6:e18093.
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:1972–3.
Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74.
Anisimova M, Gil M, Dufayard J-F, Dessimoz C, Gascuel O. Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Syst Biol. 2011;60:685–99.
Morariu VI, Srinivasan BV, Raykar VC, Duraiswami R, Davis LS. Automatic online tuning for fast Gaussian summation. In: Koller D, Schuurmans D, Bengio Y, Bottou L, editors. Advances in Neural Information Processing Systems 21: Curran Associates, Inc; 2009. p. 1113–20.
Cai K, Klein-Seetharaman J, Hwa J, Hubbell WL, Khorana HG. Structure and function in rhodopsin: effects of disulfide cross-links in the cytoplasmic face of rhodopsin on transducin activation and phosphorylation by rhodopsin kinase. Biochemistry. 1999;38:12893–8.
Haskell-Luevano C, Cone RD, Monck EK, Wan YP. Structure activity studies of the melanocortin-4 receptor by in vitro mutagenesis: identification of agouti-related protein (AGRP), melanocortin agonist and synthetic peptide antagonist interaction determinants. Biochemistry. 2001;40:6164–79.
Fritze O, Filipek S, Kuksa V, Palczewski K, Hofmann KP, Ernst OP. Role of the conserved NPxxY(x)5,6F motif in the rhodopsin ground state and during activation. Proc Natl Acad Sci USA. 2003;100:2290–5.
We would like to thank Ronald Petie for help collecting the animals and Jérôme Delroisse for his expertise and assistance in reviewing the manuscripts. Additionally, we would like to thank Michigan State University for the use of their High Performance Computing Cluster (HPCC).
This work was partially supported by the Marie Curie ITN “Neptune” (grant number: 317172). AG was funded by the Danish Research Council # 4002–00284. EUL was supported by a grant of the German Research Foundation (DFT), grant no: UL 428/2–1. Funding bodies had no role in the design of the study and collection, analysis, nor interpretation of data nor in writing the manuscript.
Availability of data and materials
Raw sequencing data has been submitted to NCBI SRA under the accession number PRJNA422510. Additionally, scripts to reproduce the analysis and fasta files containing all sequences can be found in the data section of the github repository https://github.com/elijahlowe/Acanthaster_opsins and in supplementary material.
Ethics approval and consent to participate
No special permits are needed for the marine invertebrate species used in this work and that collection happened in accordance with local and international laws.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fasta file containing all 169 opsin and melatonin sequences used to generate phylogenetic tree. (FA 60 kb)
Figure S1. Bayesian tree of opsins sequences with support. Tree was generated using using MrBayes (v3.2.5)  50 million generations, with the GTR + G amino acid substitution model. (PNG 4279 kb)
Figure S2. Differential gene expression in all A. planci tissue samples with the opsins highlighted: c-opsins in red, go-opsins in green, chaopsin in black, neuropsin in purple, peropsin in yellow, r-opsin in blue and RGR opsin in orange. The y-axis in the log2 fold-change, as the distance from the y-axis increases the more differentially expressed a gene is in one tissue versus the other. The x-axis represents counts per million (CPM), an increase on this axis shows genes with more reads counts. (PNG 1039 kb)
Figure S3. Graphical representation of QPCR results. a) Relative expression to the elongation factor Ef1A gene of synaptotagmin 1 (syt1), r-opsin, chaopsin and c-opsin1a genes in mixed tissues, tube feet (TF), terminal tube feet (eye) and radial nerve (RN). Calculations from QPCR raw data used the formula 1.9-ΔCt, where 1.9 is the multiplier for amplification per PCR cycle, and ΔCt is the threshold cycle difference with Ef1a found for that sample. b) Fold enrichment of syt1, c-opsin1a, r-opsin and chaopsin genes in tube feet (TF), terminal tube feet (eye) and radial nerve (RN) compared to the mixed tissues. Calculations from QPCR raw data used the formula 1.9-ΔΔCt, where 1.9 is the multiplier for amplification per PCR cycle, and ΔΔCt is the ΔCt difference between mixed tissues and the other tissues. Data for each gene were normalized against the housekeeping Ef1a. All quantitative measurements were done in triplicate on the cDNA obtained from each biological replica tissue sample. Average calculations over the three technical replicas ± standard deviations are reported for each gene with error bars. (JPG 244 kb)
Figure S4. Phylogenomic tree of 92 G-protein alpha subunit sequences. A. planci sequences ID’s are highlighted in red. Of the 10 A. planci sequences 3 classified as Gαs, 1 as Gαo, 4 as Gαi, 1 as Gαq, and 1 as Gα12. (JPG 308 kb)
Figure S5. Differential gene expression in the A. planci eye samples compared to mixed tissues, with the G protein alpha subunits highlighted: 3 Gαs (red), 1 Gαo (green), 4 Gαi (black), 1 Gαq (purple), and 1 Gα12 (yellow). All identified g protein alpha subunits with the exception of 1 Gαs (gbr.231.19.t1), 1 Gαi (gbr.143.10.t1) and the Gα12 are show higher expression in the eyes of A. planci compared to the mixed tissue samples. (JPG 462 kb)