Escudero Castelán N , Semmens DC , Guerra LAY , Zandawala M , Dos Reis M , Slade SE , Scrivens JH , Zampronio CG , Jones AM , Mirabeau O , Elphick MR .

	Fig. 1. Phylogenetic analysis of bilaterian kisspeptin-type receptors, including A. rubens receptors ArKPR1-11. The phylogenetic tree, which was constructed using the maximum-likelihood method and rooted with galanin/allatostatin-A-type receptors as an outgroup, comprises three distinct clades, with bootstrap support > 90. Clade 1 includes ArKPR2-4 (ArubKPR2-4) and chordate kisspeptin-type receptors. Clade 2 includes ArKPR1 and protostome (annelid, mollusk) kisspeptin-type receptors. Clade 3 comprises three sub-clades: clade 3.1 includes ArKPR8-9 (ArubKPR8-9) and related receptors in other echinoderms and the hemichordate S. kowalevskii, clade 3.2 includes ArKPR5-7 (ArubKPR5-7) and related receptors in other echinoderms and the hemichordate S. kowalevskii, and clade 3.3 includes ArKPR10-11 (ArubKPR10-11) and related receptors in other echinoderms. The stars represent branch support (bootstrap 1000 replicates) and the pastel-colored backgrounds represent taxonomic groups (see key). The arrowheads label the eleven A. rubens kisspeptin-type receptors ArKPR1-11. The scale bar indicates amino acid substitutions per site. Receptors for which ligands have been identified experimentally in this study or other published studies are colored in blue. Species names are as follows: Apla, Acanthaster planci; Ajap, Apostichopus japonicus; Arub, Asterias rubens, Anjap Annessia japonica; Bbel, Branchiostoma belcheri; Ctel, Capitella teleta; Cgig, Crassostrea gigas; Hsap, Homo sapiens; Lcha, Latimeria chalumnae; Locu, Lepisosteus oculatus; Lgig, Lottia gigantea; Mmus, Mus musculus; Pdum, Platynereis dumerilii; Pbiv, Python bivittatus; Skow, Saccoglossus kowalevskii; Spur, Strongylocentrotus purpuratus; Xtro, Xenopus tropicalis. Accession numbers for the sequences of the receptors included in this tree are listed in Additional file 3
	Fig. 2. Comparative analysis of the structure and chromosomal location of genes encoding kisspeptin-type receptors. A Comparison of the exon/intron structure of genes encoding kisspeptin-type receptors in A. rubens (Ar), Homo sapiens (Hs), Xenopus tropicalis (Xt), Lepisosteus oculatus (Lo), and Crassostrea gigas (Cg). Exons are shown as rectangles, with non-coding regions white and protein-coding regions black or colored (regions encoding predicted transmembrane domains 1–7 are shown in red, orange, yellow, green, blue, dark purple, and light purple, respectively). Introns are shown as lines, with length underneath. The two introns highlighted with a pink bracket are a conserved feature (see also Additional file 9). Accession numbers for the transcripts/genes represented in this figure are in Additional file 2. B Diagram showing the chromosomal locations of genes encoding ArKPR1-11 in A. rubens. Genes encoding (i) ArKPR1-4, (ii) ArKPR5-9; ArKPR11, and (iii) ArKPR10 are located on chromosomes 15, 10, and 4, respectively. Curved lines linking genes correspond to clades in Fig. 1: light purple is clade 1, green is clade 3.1, orange is clade 3.2, and black is clade 3.3. C Diagram showing the relative locations and orientations of genes encoding (i) ArKPR1-4 (clade 1 and 2 of Fig. 1) in a region of chromosome 15, (ii) ArKPR5-7 (clade 3.2 of Fig. 1) in a region of chromosome 10, (iii) ArKPR8-9 (clade 3.1 and 3.3 of Fig. 1) in a region of chromosome 10, and (iv) ArKPR10 (clade 3.3 of Fig. 1). The length (including exons and introns) and orientation of the genes is indicated by the purple arrows and the distance between genes (number of bases) is stated underneath the intervening black lines. The absence (0 genes) or presence (number of genes) of other genes in between those encoding the kisspeptin-type receptors is also stated (see Additional file 11 for data)
	Fig. 3. Neuropeptides identified as candidate ligands for kisspeptin-type receptors in echinoderms aligned with chordate kisspeptins. A Sequences of the A. rubens kisspeptin-type precursors ArKPP1 and ArKPP2 and SALMFamide-type precursors ArL-SALMFaP and ArF-SALMFaP. N-terminal signal peptide is shown in blue, neuropeptides predicted or shown to be derived from these proteins are shown in red (with C-terminal glycine that is a substrate for amidation shown in orange) and monobasic or dibasic cleavage sites are shown in green. Neuropeptides (red) are named in accordance with the precursor they are derived from and their relative position in the precursor, which explains the nomenclature used in B–D. Likewise, annotated precursor sequences for other neuropeptides included in B and C are shown in Additional file 13. B Alignment of ArKP1.1, ArKP1.2, and KP1-type neuropeptides from other echinoderms (blue) with chordate kisspeptin-type neuropeptides (yellow). C Alignment ArKP2.1, ArKP2.2, and KP2-type neuropeptides from other echinoderms (blue) with chordate kisspeptin-type neuropeptides (yellow). D Alignment of L-type SALMFamide precursor-derived peptides from A. rubens (ArS1.1, ArS1.3, ArS1.4, and S1.7) and other echinoderms (blue) with chordate kisspeptin-type neuropeptides (yellow). E Alignment of F-type SALMFamide precursor-derived peptides from A. rubens (ArS2.2, ArS2.3, and S2.8) and other echinoderms (blue) with chordate kisspeptin-type neuropeptides (yellow). Conserved residues are highlighted in black or gray. Abbreviations: Ar, Asterias rubens; Sp, Strongylocentrotus purpuratus; Aj, Apostichopus japonicus, Hs, Homo sapiens, Lo, Lepisosteus oculatus, Bf, Branchiostoma floridae. The accession numbers for the precursor sequences used for this figure are in Additional file 13
	Fig. 4. Structure and chromosomal location of genes encoding candidate ligands for A. rubens kisspeptin-type receptors. A Comparison of the exon/intron structure of genes encoding precursors of candidate ligands for kisspeptin-type receptors in A. rubens (Ar) and genes encoding kisspeptin-type precursors in Homo sapiens (Hs), Xenopus tropicalis (Xt), and Lepisosteus oculatus (Lo). Exons are shown as rectangles, with non-coding regions white and protein-coding regions black or colored (regions encoding the N-terminal signal peptide, neuropeptides and predicted monobasic/dibasic cleavage sites are shown in blue, red, and green, respectively). Introns are shown as lines, with length underneath. The presence of a phase 1 intron interrupting the coding sequence between exons encoding the N-terminal signal peptide and one or more neuropeptides is a feature that is conserved between the A. rubens KPP1 gene and vertebrate kisspeptin precursor genes, providing evidence of orthology, whereas ArKPP2 is encoded by a single exon. The presence of an intron that interrupts the coding sequence between an exon encoding the N-terminal signal peptide and an exon encoding multiple neuropeptides is also a feature of the two A. rubens SALMFamide precursor genes, but this is a phase 0 intron. The presence and position of introns shown here in A. rubens neuropeptide precursor genes are conserved in orthologs in other echinoderm species (see Additional file 17). The accession numbers for the sequences of the precursors shown in this figure are listed in Additional file 17. B Diagram showing the chromosomal locations of four genes encoding precursors of candidate ligands for kisspeptin-type receptors in A. rubens: ArKPP1, ArKPP2, ArL-SALMFaP, and ArF-SALMFaP. Genes encoding ArKPP1 and ArKPP2 are located on chromosome 22 and genes encoding ArL-SALMFaP and ArF-SALMFaP are located on chromosome 13, indicating that these are two pairs of paralogous genes that originated by gene duplication and intrachromosomal translocation
	Fig. 5. Identification of neuropeptides that act as ligands for A. rubens kisspeptin-type receptors. A ArKPR1 is only activated by ArKP1.2 (EC50 = 1.16 × 10−8 M). B ArKPR3 is only activated by ArKP2.2 (EC50 = 5.80 × 10−10 M). C ArKPR6 is activated by multiple neuropeptides at high concentrations (> 1 µM) but is only activated by a ‘cocktail’ of neuropeptides derived from the F-type SALMFamide precursor at lower concentrations (EC50 = 1.33 × 10−9 M). D ArKPR7 is activated by multiple neuropeptides at high concentrations (> 1 µM) but is only activated by a ‘cocktail’ of neuropeptides derived from the L-type SALMFamide precursor at lower concentrations (EC50 = 2.43 × 10−9 M). E ArKPR8 is activated by both ArKP1.1 and ArKP1.2 at high concentrations (> 1 µM) but is only activated by ArKP1.1 at lower concentrations (EC50 = 3.88 × 10−7 M). F ArKPR9 is activated by ArKP1.1, ArKP1.2, and a ‘cocktail’ of neuropeptides derived from the L-type SALMFamide precursor, but ArKP1.2 is the most potent ligand for this receptor (EC50 = 3.07 × 10−.8 M). Key: green = ArKP1.1; pink = ArKP1.2; purple = ArKP2.2; blue = ‘cocktail’ of neuropeptides derived from the L-type SALMFamide precursor (ArS1.1–7); red = ‘cocktail’ of neuropeptides derived from the F-type SALMFamide precursor (ArS2.1–8). Each point represents mean values (± S.E.M.) from at least four independent experiments, with each experiment performed in triplicate. Luminescence is expressed as a percentage of the maximal response observed in each experiment. The source data for these experiments is provided in Additional file 15

Ayres, BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics. 2012, Pubmed

Ayres, BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics. 2012, Pubmed
Baubet, Chimeric green fluorescent protein-aequorin as bioluminescent Ca2+ reporters at the single-cell level. 2000, Pubmed
Biran, Molecular identification and functional characterization of the kisspeptin/kisspeptin receptor system in lower vertebrates. 2008, Pubmed
Burke, A genomic view of the sea urchin nervous system. 2006, Pubmed , Echinobase
Cai, Biochemical, Anatomical, and Pharmacological Characterization of Calcitonin-Type Neuropeptides in Starfish: Discovery of an Ancient Role as Muscle Relaxants. 2018, Pubmed , Echinobase
Capella-Gutiérrez, trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. 2009, Pubmed
Cardoso, Allatostatin-type A, kisspeptin and galanin GPCRs and putative ligands as candidate regulatory factors of mantle function. 2016, Pubmed
Chen, Neuropeptide precursors and neuropeptides in the sea cucumber Apostichopus japonicus: a genomic, transcriptomic and proteomic analysis. 2019, Pubmed , Echinobase
Chen, TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. 2020, Pubmed
Chieu, Identification of neuropeptides in the sea cucumber Holothuria leucospilota. 2019, Pubmed , Echinobase
Comninos, Emerging Roles of Kisspeptin in Sexual and Emotional Brain Processing. 2018, Pubmed
Dong, Intraperitoneal Treatment of Kisspeptin Suppresses Appetite and Energy Expenditure and Alters Gastrointestinal Hormones in Mice. 2020, Pubmed
Dudek, Kisspeptin and Metabolism: The Brain and Beyond. 2018, Pubmed
Dufour, Origin and Evolution of the Neuroendocrine Control of Reproduction in Vertebrates, With Special Focus on Genome and Gene Duplications. 2020, Pubmed
Díaz-Miranda, Pharmacological action of the heptapeptide GFSKLYFamide in the muscle of the sea cucumber Holothuria glaberrima (Echinodermata). 1995, Pubmed , Echinobase
Elphick, The Evolution and Variety of RFamide-Type Neuropeptides: Insights from Deuterostomian Invertebrates. 2014, Pubmed , Echinobase
Elphick, The SALMFamides: a new family of neuropeptides isolated from an echinoderm. 1991, Pubmed , Echinobase
Elphick, SALMFamide salmagundi: the biology of a neuropeptide family in echinoderms. 2014, Pubmed , Echinobase
Elphick, Isolation of the neuropeptide SALMFamide-1 from starfish using a new antiserum. 1991, Pubmed , Echinobase
Elphick, Reconstructing SALMFamide Neuropeptide Precursor Evolution in the Phylum Echinodermata: Ophiuroid and Crinoid Sequence Data Provide New Insights. 2015, Pubmed , Echinobase
Elphick, Distribution and action of SALMFamide neuropeptides in the starfish Asterias rubens. 1995, Pubmed , Echinobase
Franssen, The kisspeptin receptor: A key G-protein-coupled receptor in the control of the reproductive axis. 2018, Pubmed
Frickey, CLANS: a Java application for visualizing protein families based on pairwise similarity. 2004, Pubmed
Guindon, New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. 2010, Pubmed
Guindon, A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. 2003, Pubmed
Hoang, UFBoot2: Improving the Ultrafast Bootstrap Approximation. 2018, Pubmed
Jékely, Global view of the evolution and diversity of metazoan neuropeptide signaling. 2013, Pubmed
Kim, Coevolution of the spexin/galanin/kisspeptin family: Spexin activates galanin receptor type II and III. 2014, Pubmed
Kotani, The metastasis suppressor gene KiSS-1 encodes kisspeptins, the natural ligands of the orphan G protein-coupled receptor GPR54. 2001, Pubmed
Koziol, Precursors of neuropeptides and peptide hormones in the genomes of tardigrades. 2018, Pubmed
Krzywinski, Circos: an information aesthetic for comparative genomics. 2009, Pubmed
Lapatto, Kiss1-/- mice exhibit more variable hypogonadism than Gpr54-/- mice. 2007, Pubmed
Li, Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. 2006, Pubmed
Lin, Pedal peptide/orcokinin-type neuropeptide signaling in a deuterostome: The anatomy and pharmacology of starfish myorelaxant peptide in Asterias rubens. 2017, Pubmed , Echinobase
Lin, Cellular localization of relaxin-like gonad-stimulating peptide expression in Asterias rubens: New insights into neurohormonal control of spawning in starfish. 2017, Pubmed , Echinobase
Liu, IBS: an illustrator for the presentation and visualization of biological sequences. 2015, Pubmed
Mayorova, Localization of Neuropeptide Gene Expression in Larvae of an Echinoderm, the Starfish Asterias rubens. 2016, Pubmed , Echinobase
Menschaert, A hybrid, de novo based, genome-wide database search approach applied to the sea urchin neuropeptidome. 2010, Pubmed , Echinobase
Mills, Kisspeptin and the control of emotions, mood and reproductive behaviour. 2018, Pubmed
Minh, IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. 2020, Pubmed
Mirabeau, Molecular evolution of peptidergic signaling systems in bilaterians. 2013, Pubmed
Mita, Inhibitory effect of a SALMFamide neuropeptide on secretion of gonad-stimulating substance from radial nerves in the starfish Asterina pectinifera. 2004, Pubmed , Echinobase
Muir, AXOR12, a novel human G protein-coupled receptor, activated by the peptide KiSS-1. 2001, Pubmed
Nash, The KISS1 metastasis suppressor: mechanistic insights and clinical utility. 2006, Pubmed
Newman, Tissue distribution of the SALMFamide neuropeptides S1 and S2 in the starfish Asterias rubens using novel monoclonal and polyclonal antibodies. I. Nervous and locomotory systems. 1995, Pubmed , Echinobase
Nguyen, Insights Into Sexual Maturation and Reproduction in the Norway Lobster (Nephrops norvegicus) via in silico Prediction and Characterization of Neuropeptides and G Protein-coupled Receptors. 2018, Pubmed
Nguyen, IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. 2015, Pubmed
Odekunle, Ancient role of vasopressin/oxytocin-type neuropeptides as regulators of feeding revealed in an echinoderm. 2019, Pubmed , Echinobase
Ohtaki, Metastasis suppressor gene KiSS-1 encodes peptide ligand of a G-protein-coupled receptor. 2001, Pubmed
Pasquier, Molecular evolution of GPCRs: Kisspeptin/kisspeptin receptors. 2014, Pubmed
Pasquier, Comparative evolutionary histories of kisspeptins and kisspeptin receptors in vertebrates reveal both parallel and divergent features. 2012, Pubmed
Pasquier, Looking for the bird Kiss: evolutionary scenario in sauropsids. 2014, Pubmed
Pinilla, Kisspeptins and reproduction: physiological roles and regulatory mechanisms. 2012, Pubmed
Ranwez, Improvement of distance-based phylogenetic methods by a local maximum likelihood approach using triplets. 2002, Pubmed
Ronquist, MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. 2012, Pubmed
Rowe, Discovery of a second SALMFamide gene in the sea urchin Strongylocentrotus purpuratus reveals that L-type and F-type SALMFamide neuropeptides coexist in an echinoderm species. 2010, Pubmed , Echinobase
Seminara, The GPR54 gene as a regulator of puberty. 2003, Pubmed
Semmens, Transcriptomic identification of starfish neuropeptide precursors yields new insights into neuropeptide evolution. 2016, Pubmed , Echinobase
Semmens, Discovery of a novel neurophysin-associated neuropeptide that triggers cardiac stomach contraction and retraction in starfish. 2013, Pubmed , Echinobase
Smith, The neuropeptidome of the Crown-of-Thorns Starfish, Acanthaster planci. 2017, Pubmed , Echinobase
Suwansa-Ard, Transcriptomic discovery and comparative analysis of neuropeptide precursors in sea cucumbers (Holothuroidea). 2018, Pubmed , Echinobase
Tian, Urbilaterian origin of paralogous GnRH and corazonin neuropeptide signalling pathways. 2016, Pubmed , Echinobase
Tinoco, Ancient role of sulfakinin/cholecystokinin-type signalling in inhibitory regulation of feeding processes revealed in an echinoderm. 2021, Pubmed , Echinobase
Tinoco, Characterization of NGFFYamide Signaling in Starfish Reveals Roles in Regulation of Feeding Behavior and Locomotory Systems. 2018, Pubmed , Echinobase
Tolson, Impaired kisspeptin signaling decreases metabolism and promotes glucose intolerance and obesity. 2014, Pubmed
Veenstra, Similarities between decapod and insect neuropeptidomes. 2016, Pubmed
Veenstra, Coleoptera genome and transcriptome sequences reveal numerous differences in neuropeptide signaling between species. 2019, Pubmed
Wang, Demonstration of a Functional Kisspeptin/Kisspeptin Receptor System in Amphioxus With Implications for Origin of Neuroendocrine Regulation. 2017, Pubmed
Wang, Existence and functions of a kisspeptin neuropeptide signaling system in a non-chordate deuterostome species. 2020, Pubmed , Echinobase
Yañez-Guerra, Discovery and functional characterisation of a luqin-type neuropeptide signalling system in a deuterostome. 2018, Pubmed , Echinobase
Yañez-Guerra, Echinoderms provide missing link in the evolution of PrRP/sNPF-type neuropeptide signalling. 2020, Pubmed , Echinobase
Zandawala, Discovery of novel representatives of bilaterian neuropeptide families and reconstruction of neuropeptide precursor evolution in ophiuroid echinoderms. 2017, Pubmed , Echinobase
Zhang, Molecular and functional characterization of somatostatin-type signalling in a deuterostome invertebrate. 2020, Pubmed , Echinobase
Zhang, Somatostatin-type and allatostatin-C-type neuropeptides are paralogous and have opposing myoregulatory roles in an echinoderm. 2022, Pubmed , Echinobase
d'Anglemont de Tassigny, Hypogonadotropic hypogonadism in mice lacking a functional Kiss1 gene. 2007, Pubmed
de Roux, Hypogonadotropic hypogonadism due to loss of function of the KiSS1-derived peptide receptor GPR54. 2003, Pubmed