Global view of the evolution and diversity of metazoan neuropeptide signaling

Abstract

Neuropeptides are signaling molecules that commonly act via G protein-coupled receptors (GPCRs) and are generated in neurons by proneuropeptide (pNP) cleavage. Present in both cnidarians and bilaterians, neuropeptides represent an ancient and widespread mode of neuronal communication. Due to the inherent difficulties of analyzing highly diverse and repetitive pNPs, the relationships among different families are often elusive. Using similarity-based clustering and sensitive similarity searches, I obtained a global view of metazoan pNP diversity and evolution. Clustering revealed a large and diffuse network of sequences connected by significant sequence similarity encompassing one-quarter of all families. pNPs belonging to this cluster were also identified in the early-branching neuronless animal Trichoplax adhaerens. Clustering of neuropeptide GPCRs identified several orthology groups and allowed the reconstruction of the phyletic distribution of receptor families. GPCR phyletic distribution closely paralleled that of pNPs, indicating extensive conservation and long-term coevolution of receptor-ligand pairs. Receptor orthology and intermediate sequences also revealed the homology of pNPs so far considered unrelated, including allatotropin and orexin. These findings, together with the identification of deuterostome achatin and luqin and protostome opioid pNPs, extended the neuropeptide complement of the urbilaterian. Several pNPs were also identified from the hemichordate Saccoglossus kowalevskii and the cephalochordate Branchiostoma floridae, elucidating pNP evolution in deuterostomes. Receptor-ligand conservation also allowed ligand predictions for many uncharacterized GPCRs from nonmodel species. The reconstruction of the neuropeptide-signaling repertoire at deep nodes of the animal phylogeny allowed the formulation of a testable scenario of the evolution of animal neuroendocrine systems.

Keywords: Petromyzon; Platynereis; proenkephalin.

Fig. 1.

BLOSUM62 cluster map of metazoan pNP families. Nodes correspond to pNPs and are colored based on taxonomy. Edges correspond to BLAST connections of P value >1e-5.

Fig. 2.

PAM30 map of the CC. The largest cluster in the PAM30 map was defined using linkage clustering and optimized separately. Nodes correspond to pNPs and are colored based on taxonomy. Edges represent BLAST connections of P value > 1e-5. Color code is the same as in Fig. 1.

Fig. 3.

BLOSUM62 cluster map of class-A neuropeptide GPCRs. Nodes correspond to class-A GPCRs and are colored based on taxonomy. Edges represent BLAST connections of P value > 1e-30. Color code is the same as in Fig. 1.