Molecular evolution of neuropeptide receptors with regard to maintaining high affinity to their authentic ligands

Emerging roles of APLN and APELA in the physiology and pathology of the female reproductive system

APLN, APELA and their common receptor APLNR (composing the apelinergic axis) have been described in various species with extensive body distribution and multiple physiological functions. Recent studies have witnessed emerging intracellular cascades triggered by APLN and APELA which play crucial roles in female reproductive organs, including hypothalamus-pituitary-gonadal axis, ovary, oviduct, uterus and placenta. However, a comprehensive summary of APLN and APELA roles in physiology and pathology of female reproductive system has not been reported to date. In this review, we aim to concentrate on the general characteristics of APLN and APELA, as well as their specific physiological roles in female reproductive system. Meanwhile, the pathological contexts of apelinergic axis dysregulation in the obstetrics and gynecology are also summarized here, suggesting its potential prospect as a diagnostic biomarker and/or therapeutic intervention in the polycystic ovary syndrome, ovarian cancer, preeclampsia and gestational diabetes mellitus.

The vasotocinergic system and its role in the regulation of stress in birds

The regulation of stress in birds includes a complex interaction of neural systems affecting the hypothalamic-pituitary-adrenal (HPA) axis. In addition to the hypothalamic paraventricular nucleus, a structure called the nucleus of the hippocampal commissure likewise affects the output of pituitary stress hormones and appears to be unique to avian species. Within the anterior pituitary, the avian V1a and V1b receptors were found in corticotropes. Based on our studies with central administration of hormones in the chicken, corticotropic releasing hormone (CRH) is a more potent ACTH secretagogue than arginine vasotocin (AVT). In contrast, when applied peripherally, AVT is more efficacious. Co-administration of AVT and CRH peripherally, resulted in a synergistic stimulation of corticosterone release. Data suggest receptor oligomerization as one possible mechanism. In birds, vasotocin receptors associated with stress responses include the V1a and V1b receptors. Three-dimensional, homology-based structural models of the avian V1aR were built to test agonists and antagonists for each receptor that were screened by molecular docking to map their binding sites on each receptor. Additionally, binding affinity values for each available peptide antagonist to the V1aR and V1bR were determined. An anterior pituitary primary culture system was developed to determine how effective each antagonist blocked the function of each receptor in culture when stimulated by a combination of AVT/CRH administration. Use of an antagonist in subsequent in vivo studies identified the V1aR in regulating food intake in birds. The V1aR was likewise found in circumventricular organs of the brain, suggesting a possible function in stress.

Identification of avian vasotocin receptor subtype-specific antagonists involved in the stress response of the chicken, Gallus gallus

Vasotocin 1a and 1b receptors (V1aR and V1bR) have been shown to play important roles in the neuroendocrine regulation of stress responses via the anterior pituitary (AP) of birds. To identify effective subtype-specific antagonists for the chicken V1aR (cV1aR) and cV1bR, potential antagonists to the mammalian V1R were screened against the cV1aR and cV1bR 3D structural models by molecular docking analysis with determination of binding pocket/amino acid residues involved in the interaction. The antagonistic effects of the selected ligands were examined by measuring pro-opiomelanocortin (POMC) heteronuclear RNA (hnPOMC) levels following the in vitro stress administration to primary chicken AP cells. Results of in silico analysis showed that the Manning compound and several other antagonists were bound to cV1bR with higher affinity than the natural agonist, arginine vasotocin (AVT). Similarities and differences in the antagonist-receptor binding interface with receptors were characterized for each ligand. Non-peptide mammalian V1bR antagonists, SSR-149415 and L-368899, were shown to be effective and had an additive effect in blocking POMC hnRNA expression in pituitary cell culture studies. SR-49059 antagonized the effect(s) of AVT/CRH on the downregulation of the cV1aR and the upregulation of the cCRH-R2 expression but not the cV1bR and cCRH-R1. The Manning compound antagonized the downregulation of cV1aR, cV1bR and cCRH-R1 and the upregulation of cCRH-R2 expression. The specificity of antagonists apparently resulted from unique differences in the interacting residues and their binding affinities. Collectively, these results provide valuable leads for future development of novel compounds capable of blocking or attenuating the AP stress response of avian species and perhaps other non-mammalian vertebrates as well.

The directed evolution of ligand specificity in a GPCR and the unequal contributions of efficacy and affinity

G protein-coupled receptors (GPCRs) must discriminate between hundreds of related signal molecules. In order to better understand how GPCR specificity can arise from a common promiscuous ancestor, we used laboratory evolution to invert the specificity of the Saccharomyces cerevisiae mating receptor Ste2. This GPCR normally responds weakly to the pheromone of the related species Kluyveromyces lactis, though we previously showed that mutation N216S is sufficient to make this receptor promiscuous. Here, we found that three additional substitutions, A265T, Y266F and P290Q, can act together to confer a novel specificity for K. lactis pheromone. Unlike wild-type Ste2, this new variant does not rely on differences in binding affinity to discriminate against its non-preferred ligand. Instead, the mutation P290Q is critical for suppressing the efficacy of the native pheromone. These two alternative methods of ligand discrimination were mapped to specific amino acid positions on the peptide pheromones. Our work demonstrates that changes in ligand efficacy can drive changes in GPCR specificity, thus obviating the need for extensive binding pocket re-modeling.

Sensitive Periods, Vasotocin-Family Peptides, and the Evolution and Development of Social Behavior

Nonapeptides, by modulating the activity of neural circuits in specific social contexts, provide an important mechanism underlying the evolution of diverse behavioral phenotypes across vertebrate taxa. Vasotocin-family nonapeptides, in particular, have been found to be involved in behavioral plasticity and diversity in social behavior, including seasonal variation, sexual dimorphism, and species differences. Although nonapeptides have been the focus of a great deal of research over the last several decades, the vast majority of this work has focused on adults. However, behavioral diversity may also be explained by the ways in which these peptides shape neural circuits and influence social processes during development. In this review, I synthesize comparative work on vasotocin-family peptides during development and classic work on early forms of social learning in developmental psychobiology. I also summarize recent work demonstrating that early life manipulations of the nonapeptide system alter attachment, affiliation, and vocal learning in zebra finches. I thus hypothesize that vasotocin-family peptides are involved in the evolution of social behaviors through their influence on learning during sensitive periods in social development.

Evolutionary and Comparative Genomics to Drive Rational Drug Design, with Particular Focus on Neuropeptide Seven-Transmembrane Receptors

Seven transmembrane receptors (7TMRs), also known as G protein-coupled receptors, are popular targets of drug development, particularly 7TMR systems that are activated by peptide ligands. Although many pharmaceutical drugs have been discovered via conventional bulk analysis techniques the increasing availability of structural and evolutionary data are facilitating change to rational, targeted drug design. This article discusses the appeal of neuropeptide-7TMR systems as drug targets and provides an overview of concepts in the evolution of vertebrate genomes and gene families. Subsequently, methods that use evolutionary concepts and comparative analysis techniques to aid in gene discovery, gene function identification, and novel drug design are provided along with case study examples.

Molecular cloning, sequencing and phylogeny of vasotocin receptor genes in the air-breathing catfish Heteropneustes fossilis with sex dimorphic and seasonal variations in tissue expression

Vasotocin (VT) is the ortholog of vasopressin (VP) in non-mammalian vertebrates and is known for multiple functions. Teleost fishes have a complete repertoire of known VP/VT receptor subtypes (vasopressin type, VR): two V1A subtypes (V1Aa and V1Ab or V1a1 and V1a2) and five V2 subtypes (V2A1, V1A2, V2B1, V2B2 and V2C). Full-length cDNAs of v1a1, v1a2 and v2 (v2a1) with ORFs of 1,308, 1,137 and 1,527 bp, respectively, were cloned and characterized in the catfish Heteropneustes fossilis (Siluriformes, Ostariophysi). BLAST analysis revealed that the genes encoded three VT receptors, V1a1, V1a2 and V2 of 436, 379 and 509 amino acid residues, respectively. The predicted proteins showed typical features of the seven-transmembrane domain receptor core structure with hallmark triplets Asp-Arg-Tyr/Asp-Arg-His (DRY/DRH) and the variable intracellular loop III of vertebrate neurohypophysial hormone receptors. Phylogenetic analysis of the deduced protein sequences revealed that they clustered with the V1Aa, V1Ab and V2A1, respectively, of other teleosts. The V2R has a sequence identity of 70-76% with V2A1 than with the V2B type (sequence identity 43-49%). Semiquantitative PCR analysis showed that the receptor gene transcripts were expressed ubiquitously in the tissues examined (brain, pituitary, gonads, liver, muscle, kidney and gills) and displayed sex and seasonal fluctuations in a tissue-specific manner. The results form a basis for functional studies on the VT receptors in the catfish.

Ligand binding pocket formed by evolutionarily conserved residues in the glucagon-like peptide-1 (GLP-1) receptor core domain

Glucagon-like peptide-1 (GLP-1) plays a pivotal role in glucose homeostasis through its receptor GLP1R. Due to its multiple beneficial effects, GLP-1 has gained great attention for treatment of type 2 diabetes and obesity. However, little is known about the molecular mechanism underlying the interaction of GLP-1 with the heptahelical core domain of GLP1R conferring high affinity ligand binding and ligand-induced receptor activation. Here, using chimeric and point-mutated GLP1R, we determined that the evolutionarily conserved amino acid residue Arg(380) flanked by hydrophobic Leu(379) and Phe(381) in extracellular loop 3 (ECL3) may have an interaction with Asp(9) and Gly(4) of the GLP-1 peptide. The molecular modeling study showed that Ile(196) at transmembrane helix 2, Met(233) at ECL1, and Asn(302) at ECL2 of GLP1R have contacts with His(1) and Thr(7) of GLP-1. This study may shed light on the mechanism underlying high affinity interaction between the ligand and the binding pocket that is formed by these conserved residues in the GLP1R core domain.

Molecular evolution of GPCRs: GLP1/GLP1 receptors

Glucagon-like peptide 1 (GLP1) is an intestinal incretin that regulates glucose homeostasis through stimulation of insulin secretion from pancreatic β-cells and inhibits appetite by acting on the brain. Thus, it is a promising therapeutic agent for the treatment of type 2 diabetes mellitus and obesity. Studies using synteny and reconstructed ancestral chromosomes suggest that families for GLP1 and its receptor (GLP1R) have emerged through two rounds (2R) of whole genome duplication and local gene duplications before and after 2R. Exon duplications have also contributed to the expansion of the peptide family members. Specific changes in the amino acid sequence following exon/gene/genome duplications have established distinct yet related peptide and receptor families. These specific changes also confer selective interactions between GLP1 and GLP1R. In this review, we present a possible macro (genome level)- and micro (gene/exon level)-evolution mechanisms of GLP1 and GLP1R, which allows them to acquire selective interactions between this ligand-receptor pair. This information may provide critical insight for the development of potent therapeutic agents targeting GLP1R.

Variations in the expression of vasotocin and isotocin receptor genes in the gilthead sea bream Sparus aurata during different osmotic challenges

The dynamic changes in mRNA expression levels for vasotocin (AVT) and isotocin (IT) receptor gene levels were assessed in a time-course response study in immature male specimens of the gilthead sea bream (Sparus aurata) submitted to hyper- (55‰ salinity) and hypo-osmotic (5‰ salinity) challenges. Two different cDNAs for the AVT receptor and one for the IT receptor (V1a2-type and V2-type AVTR, and ITR, respectively) were cloned by screening an S. aurata brain cDNA library. Genes for these receptors were expressed differentially and is nearly ubiquitously in 26 of the examined tissues. In the gills, both environmental salinity challenges up-regulated AVTR V1a2-type gene expression concomitantly with mRNA expression protein activity of Na(+), K(+)-ATPase gene expression and protein, whereas the AVTR V2-type and cystic fibrosis transmembrane conductance regulator (CFTR) mRNA levels were associated with mRNAs environmental salinity, indicating a possible connection between AVTRs and these transporters. In kidney, AVTR V1a2-type gene expression peaked rapidly and lasted only a short time (12-24h) in response to both osmotic challenges. In contrast, AVTR V2-type mRNA levels were enhanced in specimens exposed to hyperosmotic conditions, whereas they decreased under hypoosmotic environments, suggesting an antidiuretic role related to the vasoconstriction function. In the hypothalamus, only the expression of the AVTR V2-type gene was enhanced at 7 and 14 days under both experimental conditions. In the liver, both AVTRs had increased mRNA levels, with the upregulation of their AVTR V2-type gene increasing faster than the V1a2-type. The ITR gene was not sensitive to variations of external salinity in any of the analyzed tissues. Our results demonstrate the involvement of the vasotocinergic, but not the isotocinergic, pathway as well as the hypothalamic function, in the adjustments of both osmoregulatory and metabolic processes after osmotic challenges.

Does Kisspeptin Belong to the Proposed RF-Amide Peptide Family?

Kisspeptin (KISS) plays a key role in regulating reproduction by binding to its receptor, GPR54. Because of the Arg-Phe (RF) sequence at its carboxyl terminus, KISS has been proposed to be a member of the RF-amide peptide family consisting of neuropeptide FF (NPFF), neuropeptide VF (NPVF), pyroglutamylated RF-amide peptide (QRFP), and prolactin-releasing hormone (PRLH). Evolutionary relationships of protein families can be determined through phylogenetic analysis. However, phylogenetic analysis among related peptide families often fails to provide sufficient information because only short mature peptide sequences from full preprohormone sequences are conserved. Considering the concept of the coevolution of peptide ligands and their cognate receptors, evolutionary relationships among related receptor families provide clues to explore relationships between their peptides. Although receptors for NPFF, NPVF, and QRFP are phylogenetically clustered together, receptors for PRLH and KISS are on different branches of the phylogenetic tree. In particular, KISS has been proposed to be a member of the KISS/galanin/spexin family based on synteny analysis and the phylogenetic relationship between their receptors. This article discusses the evolutionary history of the receptors for the proposed RF-amide peptide family and proposes that, from an evolutionary aspect, KISS has emerged from an ancestor, which is distinct from those of the other RF-amide peptides, and so should be classed separately.

The neuropeptide complement of the marine annelid Platynereis dumerilii

Background

The marine annelid Platynereis dumerilii is emerging as a powerful lophotrochozoan experimental model for evolutionary developmental biology (evo-devo) and neurobiology. Recent studies revealed the presence of conserved neuropeptidergic signaling in Platynereis, including vasotocin/neurophysin, myoinhibitory peptide and opioid peptidergic systems. Despite these advances, comprehensive peptidome resources have yet to be reported.

Results

The present work describes the neuropeptidome of Platynereis. We established a large transcriptome resource, consisting of stage-specific next-generation sequencing datasets and 77,419 expressed sequence tags. Using this information and a combination of bioinformatic searches and mass spectrometry analyses, we increased the known proneuropeptide (pNP) complement of Platynereis to 98. Based on sequence homology to metazoan pNPs, Platynereis pNPs were grouped into ancient eumetazoan, bilaterian, protostome, lophotrochozoan, and annelid families, and pNPs only found in Platynereis. Compared to the planarian Schmidtea mediterranea, the only other lophotrochozoan with a large-scale pNP resource, Platynereis has a remarkably full complement of conserved pNPs, with 53 pNPs belonging to ancient eumetazoan or bilaterian families. Our comprehensive search strategy, combined with analyses of sequence conservation, also allowed us to define several novel lophotrochozoan and annelid pNP families. The stage-specific transcriptome datasets also allowed us to map changes in pNP expression throughout the Platynereis life cycle.

Conclusion

The large repertoire of conserved pNPs in Platynereis highlights the usefulness of annelids in comparative neuroendocrinology. This work establishes a reference dataset for comparative peptidomics in lophotrochozoans and provides the basis for future studies of Platynereis peptidergic signaling.

Insights into the molecular evolution of oxytocin receptor ligand binding

The design and development of selective ligands for the human OT (oxytocin) and AVP (arginine vasopressin) receptors is a big challenge since the different receptor subtypes and their native peptide ligands display great similarity. Detailed understanding of the mechanism of OT's interaction with its receptor is important and may assist in the ligand- or structure-based design of selective and potent ligands. In the present article, we compared 69 OT- and OT-like receptor sequences with regards to their molecular evolution and diversity, utilized an in silico approach to map the common ligand interaction sites of recently published G-protein-coupled receptor structures to a model of the human OTR (OT receptor) and compared these interacting residues within a selection of different OTR sequences. Our analysis suggests the existence of a binding site for OT peptides within the common transmembrane core region of the receptor, but it appears extremely difficult to identify receptor or ligand residues that could explain the selectivity of OT to its receptors. We remain confident that the presented evolutionary overview and modelling approach will aid interpretation of forthcoming OTR crystal structures.

Regulatory decisions on endocrine disrupting chemicals should be based on the principles of endocrinology

For years, scientists from various disciplines have studied the effects of endocrine disrupting chemicals (EDCs) on the health and wellbeing of humans and wildlife. Some studies have specifically focused on the effects of low doses, i.e. those in the range that are thought to be safe for humans and/or animals. Others have focused on the existence of non-monotonic dose-response curves. These concepts challenge the way that chemical risk assessment is performed for EDCs. Continued discussions have clarified exactly what controversies and challenges remain. We address several of these issues, including why the study and regulation of EDCs should incorporate endocrine principles; what level of consensus there is for low dose effects; challenges to our understanding of non-monotonicity; and whether EDCs have been demonstrated to produce adverse effects. This discussion should result in a better understanding of these issues, and allow for additional dialog on their impact on risk assessment.

Identification of antagonists to the vasotocin receptor sub-type 4 (VT4R) involved in stress by molecular modelling and verification using anterior pituitary cells

The vasotocin receptor family is homologous to the mammalian vasopressin G-protein coupled receptor (GPCR) family. The vasotocin receptor 2 (VT2R) and 4 (VT4R) have recently been shown to play important role(s) in the neuroendocrine regulation of stress in birds. A homology-based structural model of VT4R of the domestic chicken, Gallus gallus, was built using the sophisticated SYBYL-X suite. The structure of VT4R built with and without extra- and intracellular unstructured loops showed a seven-helix transmembrane domain, which is a characteristic feature of GPCRs. Several agonists and antagonists were screened by molecular docking to map their potential binding sites on the structure of VT4R. Interestingly, the presence of the N-terminal, intracellular and extracellular loops and C-terminal amino acid sequences emerging from the transmembrane domains during molecular docking appeared to influence the binding interface of the peptide agonists and peptide/non-peptide antagonists on the VT4R. The presence of unstructured loops, however, did not affect the relative binding affinity ranking of the peptide antagonists to VT4R. In general, the natural ligand, arginine vasotocin and the peptide/non-peptide antagonists were observed to be more deeply buried in the receptor. Results of in vitro inhibition experiments, using cultured anterior pituitary cells, showed excellent agreement with the binding affinity of the antagonists predicted by molecular docking. The results of this study provide valuable clues for the rational design of novel pharmaceutical compounds capable of blocking or attenuating the stress response.

Variation in gene transcript profiles of two V1a-type arginine vasotocin receptors among sexual phases of bluehead wrasse (Thalassoma bifasciatum)

The neurohypophyseal hormone arginine vasotocin (AVT) mediates behavioral and reproductive plasticity in vertebrates, and has been linked to the behavioral changes associated with protogyny in the bluehead wrasse (Thalassoma bifasciatum). In this study, we sequenced full-length cDNAs encoding two distinct V1a-type AVT receptors (v1a1 and v1a2) from the bluehead wrasse, and examined variation in brain and gonadal abundance of these receptor transcripts among sexual phases. End point RT-PCR revealed that v1a1 and v1a2 transcripts varied in tissue distribution, with v1a1 receptor mRNAs at greatest levels in the telencephalon, hypothalamus, optic tectum, cerebellum and testis, and v1a2 receptor transcripts most abundant in the hypothalamus, cerebellum and gills. In the brain, v1a1 and v1a2 mRNAs both localized by in situ hybridization to the dorsal and ventral telencephalon, the preoptic area of the hypothalamus, the ventral hypothalamus and lateral recess of the third ventricle. Quantitative real-time RT-PCR revealed that relative abundance of these two receptor mRNAs varied significantly in brain and gonad with sexual phase. Relative levels of v1a2 mRNAs were greater in whole brain and isolated hypothalamus of terminal phase (TP) male wrasse compared to initial phase (IP) males or females. In the gonad, v1a1 mRNAs were at levels 2.5-fold greater in the testes of IP males - and 4-5-fold greater in the testes of TP males - compared to the ovaries of females. These results provide evidence that V1a-type AVT receptor transcript abundance in the hypothalamus and gonads of bluehead wrasse varies in patterns linked to sexual phase, and bestow a foundation for future studies investigating how differential expression of v1a1 and v1a2 teleost AVT receptors links to behavioral status and gonadal function in fish more broadly.

Molecular Coevolution of Neuropeptides Gonadotropin-Releasing Hormone and Kisspeptin with their Cognate G Protein-Coupled Receptors

The neuropeptides gonadotropin-releasing hormone (GnRH) and kisspeptin (KiSS), and their receptors gonadotropin-releasing hormone receptor (GnRHR) and kisspeptin receptor (KiSSR) play key roles in vertebrate reproduction. Multiple paralogous isoforms of these genes have been identified in various vertebrate species. Two rounds of genome duplication in early vertebrates likely contributed to the generation of these paralogous genes. Genome synteny and phylogenetic analyses in a variety of vertebrate species have provided insights into the evolutionary origin of and relationship between paralogous genes. The paralogous forms of these neuropeptides and their receptors have coevolved to retain high selectivity of the ligand–receptor interaction. These paralogous forms have become subfunctionalized, neofunctionalized, or dysfunctionalized during evolution. This article reviews the evolutionary mechanism of GnRH/GnRHR and KiSS/KiSSR, and the fate of the duplicated paralogs in vertebrates.

Structural and molecular conservation of glucagon-like Peptide-1 and its receptor confers selective ligand-receptor interaction

Glucagon-like peptide-1 (GLP-1) is a major player in the regulation of glucose homeostasis. It acts on pancreatic beta cells to stimulate insulin secretion and on the brain to inhibit appetite. Thus, it may be a promising therapeutic agent for the treatment of type 2 diabetes mellitus and obesity. Despite the physiological and clinical importance of GLP-1, molecular interaction with the GLP-1 receptor (GLP1R) is not well understood. Particularly, the specific amino acid residues within the transmembrane helices and extracellular loops of the receptor that may confer ligand-induced receptor activation have been poorly investigated. Amino acid sequence comparisons of GLP-1 and GLP1R with their orthologs and paralogs in vertebrates, combined with biochemical approaches, are useful to determine which amino acid residues in the peptide and the receptor confer selective ligand-receptor interaction. This article reviews how the molecular evolution of GLP-1 and GLP1R contributes to the selective interaction between this ligand-receptor pair, providing critical clues for the development of potent agonists for the treatment of diabetes mellitus and obesity.

The pyrokinin/ pheromone biosynthesis-activating neuropeptide (PBAN) family of peptides and their receptors in Insecta: evolutionary trace indicates potential receptor ligand-binding domains

The pyrokinin/pheromone biosynthesis-activating neuropeptide (PBAN) family of G-protein-coupled receptors and their ligands have been identified in various insects. Physiological functions of pyrokinin peptides include muscle contraction, whereas PBAN regulates, among other functions, pheromone production in moths which indicates the pleiotropic nature of these peptides. Based on the alignment of annotated genomic sequences, the pyrokinin/PBAN family of receptors have similarity with the corresponding structures of the capa or periviscerokinin receptors of insects and the neuromedin U receptors of vertebrates. In our study, evolutionary trace (ET) analysis on the insect receptor sequences was conducted to predict the putative ligand recognition and binding sites. The ET analysis of four class-specific receptors indicated several amino acid residues that are conserved in the transmembrane domains. The receptor extracellular domains exhibit several class-specific amino acid residues, which could indicate putative domains for activation of these receptors by ligand recognition and binding.

Revisiting the evolution of gonadotropin-releasing hormones and their receptors in vertebrates: secrets hidden in genomes

Gonadotropin-releasing hormone (GnRH) and its G protein-coupled receptor, GnRHR, play a pivotal role in the control of reproduction in vertebrates. To date, many GnRH and GnRHR genes have been identified in a large variety of vertebrate species using conventional biochemical and molecular biological tools in combination with bioinformatic tools. Phylogenetic approaches, primarily based on amino acid sequence identity, make it possible to classify these multiple GnRHs and GnRHRs into several lineages. Four vertebrate GnRH lineages GnRH1, GnRH2, GnRH3, and GnRH4 (for lamprey) are well established. Four vertebrate GnRHR lineages have also been proposed-three for nonmammalian GnRHRs and mammalian GnRHR2 as well as one for mammalian GnRHR1. However, these phylogenetic analyses cannot fully explain the evolutionary origins of each lineage and the relationships among the lineages. Rapid and vast accumulation of genome sequence information for many vertebrate species, together with advances in bioinformatic tools, has allowed large-scale genome comparison to explore the origin and relationship of gene families of interest. The present review discusses the evolutionary mechanism of vertebrate GnRHs and GnRHRs based on extensive genome comparison. In this article, we focus only on vertebrate genomes because of the difficulty in comparing invertebrate and vertebrate genomes due to their marked divergence.

Molecular coevolution of kisspeptins and their receptors from fish to mammals

Kisspeptin and its receptor, GPR54, play a pivotal role in vertebrate reproduction. Recent advances in bioinformatic tools combined with comparative genomics have led to the identification of a large number of kisspeptin and GPR54 genes in a variety of vertebrate species. Genome duplications may have produced at least two isoforms of both ligand (KiSS1 and KiSS2) and receptor (GPR54-1 and GPR54-2). Additional isoforms of kisspeptin (KiSS1b) and GPR54 (GPR54-1b) have been found in an amphibian species, Xenopus (Silurana) tropicalis. Here, we describe the evolutionary lineages of these kisspeptin and GPR54 isoforms using genome synteny and phylogenetic analyses, and possible molecular interactions between kisspeptin and GPR54 subtypes based on ligand-receptor selectivity. Together, kisspeptin and GPR54 provide an excellent model for understanding molecular coevolution of the peptide ligand and GPCR pairs.

Tyr1 and Ile7 of glucose-dependent insulinotropic polypeptide (GIP) confer differential ligand selectivity toward GIP and glucagon-like peptide-1 receptors

Glucagon like peptide-1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP) are incretin hormones released in response to food intake and potentiate insulin secretion from pancreatic beta cells through their distinct yet related G protein-coupled receptors, GLP1R and GIPR. While GLP-1 and GIP exhibit similarity in their N-terminal sequence and overall alpha-helical structure, GLP-1 does not bind to GIPR and vice versa. To determine which amino acid residues of these peptide ligands are responsible for specific interaction with their respective receptors, we generated mutant GIP in which several GLP-1-specific amino acid residues were substituted for the original amino acids. The potency of the mutant ligands was examined using HEK293 cells transfected with GLP1R or GIPR expression plasmids together with a cAMP-responsive element-driven luciferase (CRE-luc) reporter plasmid. A mutated GIP peptide in which Tyr(1), Ile(7), Asp(15), and His(18) were replaced by His, Thr, Glu, and Ala, respectively, was able to activate both GLP1R and GIPR with moderate potency. Replacing the original Tyr(1) and/or Ile(7) in the N-terminal moiety of this mutant peptide allowed full activation of GIPR but not of GLP1R. However, reintroducing Asp(15) and/or His(18) in the central alpha-helical region did not significantly alter the ligand potency. These results suggest that Tyr/His(1) and Ile/Thr(7) of GIP/GLP-1 peptides confer differential ligand selectivity toward GIPR and GLP1R.

Identification of multiple vasotocin receptor cDNAs in teleost fish: sequences, phylogenetic analysis, sites of expression, and regulation in the hypothalamus and gill in response to hyperosmotic challenge

Vasopressin and its homolog vasotocin regulate hydromineral balance, stress responses, and social behaviors in vertebrates. In mammals, the functions of vasopressin are mediated via three classes of membrane-bound receptors: V1a-type, V1b-type and V2-type. To date, however, only a single class of vasotocin receptor has been identified in teleost fish. Here, cDNAs encoding three putative vasotocin receptors - two distinct V1a-type receptor paralogs (V1a1 and V1a2) and a previously undescribed V2-type receptor (V2) - and a single isotocin receptor were isolated and sequenced from the Amargosa pupfish (Cyprinodon nevadensis amargosae). RT-PCR revealed that mRNAs for these receptors differed in expression patterns with V1a1 mRNAs abundant in the brain, pituitary and testis, V1a2 transcripts at greatest levels in brain, heart and muscle, V2 transcripts most common in the gills, heart and kidney, and isotocin receptor mRNAs abundant in the midbrain, pituitary and gonads. In response to an acute hyperosmotic challenge, pro-vasotocin and V2 mRNA levels in the hypothalamus decreased, while transcripts of V1a1 in the hypothalamus and V1a2 in the gills increased. Partial transcripts for structurally related V2-type, as well as multiple V1a-type, receptors were also identified in other teleosts, suggesting that multiple vasotocin receptors may be present in many Actinopterygii fishes.

Spatial autocorrelation of amino Acid replacement rates in the vasopressin receptor family

Evolutionary rates of sites can be independent of one another or correlated in some fashion. Significant spatial autocorrelation was observed for site amino acid replacement rates in vasopressin receptor family proteins (VPRs). Spatial autocorrelation of rates is the propensity of residues to lie near other residues of similar rate in the folded protein structure. Optimal correlation occurred at a distance suggesting that residues in contact had correlated rates. As another way to study the same phenomenon, VPR was partitioned into >40 x 10 A(3) contiguous spatial clusters for amino acid replacement rate estimation. Partitioning was done without preconception of functional regions of the protein and with a random partition control. Cluster rates exhibited an overdispersed distribution suggesting that rates were not randomly distributed in the spatial partitions. In tests, cluster partitioning improved maximum likelihood and Bayesian likelihood models for VPR evolution. Spatial clusters with outlier rates, or lineage-specific clusters differing in rate, proved to contain VPR features likely to be under selection. Thus the spatial autocorrelation observed is probably not just a statistical finding, but likely has an evolutionary basis in protein function.

A novel biological role of tachykinins as an up-regulator of oocyte growth: identification of an evolutionary origin of tachykininergic functions in the ovary of the ascidian, Ciona intestinalis

Tachykinins (TKs) and their receptors have been shown to be expressed in the mammalian ovary. However, the biological roles of ovarian TKs have yet to be verified. Ci-TK-I and Ci-TK-R, characterized from the protochordate (ascidian), Ciona intestinalis, are prototypes of vertebrate TKs and their receptors. In the present study, we show a novel biological function of TKs as an inducible factor for oocyte growth using C. intestinalis as a model organism. Immunostaining demonstrated the specific expression of Ci-TK-R in test cells residing in oocytes at the vitellogenic stage. DNA microarray and real-time PCR revealed that Ci-TK-I induced gene expression of several proteases, including cathepsin D, chymotrypsin, and carboxy-peptidase B1, in the ovary. The enzymatic activities of these proteases in the ovary were also shown to be enhanced by Ci-TK-I. Of particular significance is that the treatment of Ciona oocytes with Ci-TK-I resulted in progression of growth from the vitellogenic stage to the post-vitellogenic stage. The Ci-TK-I-induced oocyte growth was blocked by a TK antagonist or by protease inhibitors. These results led to the conclusion that Ci-TK-I enhances growth of the vitellogenic oocytes via up-regulation of gene expression and enzymatic activities of the proteases. This is the first clarification of the biological roles of TKs in the ovary and the underlying essential molecular mechanism. Furthermore, considering the phylogenetic position of ascidians as basal chordates, we suggest that the novel TK-regulated oocyte growth is an "evolutionary origin" of the tachykininergic functions in the ovary.

Molecular evolution of neuropeptides in the genus Drosophila

The first genomic and chemical characterization of fruit fly neuropeptides outside Drosophila melanogaster provides insights into the evolution of the neuropeptidome in this genus.

Background

Neuropeptides comprise the most diverse group of neuronal signaling molecules. They often occur as multiple sequence-related copies within single precursors (the prepropeptides). These multiple sequence-related copies have not arisen by gene duplication, and it is debated whether they are mutually redundant or serve specific functions. The fully sequenced genomes of 12 Drosophila species provide a unique opportunity to study the molecular evolution of neuropeptides.

Results

We data-mined the 12 Drosophila genomes for homologs of neuropeptide genes identified in Drosophila melanogaster. We then predicted peptide precursors and the neuropeptidome, and biochemically identified about half of the predicted peptides by direct mass spectrometric profiling of neuroendocrine tissue in four species covering main phylogenetic lines of Drosophila. We found that all species have an identical neuropeptidome and peptide hormone complement. Calculation of amino acid distances showed that ortholog peptide copies are highly sequence-conserved between species, whereas the observed sequence variability between peptide copies within single precursors must have occurred prior to the divergence of the Drosophila species.

Conclusion

We provide a first genomic and chemical characterization of fruit fly neuropeptides outside D. melanogaster. Our results suggest that neuropeptides including multiple peptide copies are under stabilizing selection, which suggests that multiple peptide copies are functionally important and not dispensable. The last common ancestor of Drosophila obviously had a set of neuropeptides and peptide hormones identical to that of modern fruit flies. This is remarkable, since drosophilid flies have adapted to very different environments.

Molecular identification and functional characterization of the kisspeptin/kisspeptin receptor system in lower vertebrates

The KISS1 gene encodes the kisspeptin neuropeptide, which activates the KISS1 receptor (KISS1R; G protein-coupled receptor 54; GPR54) and participates in neuroendocrine regulation of GnRH secretion. To study the physiological function(s) and evolutionary conservation of KISS1, we cloned opossum, Xenopus, and zebrafish kiss1 cDNAs. Processing zebrafish, Xenopus, or opossum KISS proteins would liberate a carboxy-terminal amidated peptide with 52, 54, or 53 amino acid residues, respectively. Phylogenetic analysis of all known vertebrate KISS1 peptides showed clear clustering of the sequences according to canonical vertebrate classes. The zebrafish kiss1 gene consists of two exons and one intron. Real-time PCR analysis of two kiss1R cloned from zebrafish brain found expression of kiss1, kiss1ra, and kiss1rb, with kiss1ra-more similar to other piscine Kiss1 receptors-highly expressed in the gonads and kiss1rb in other nonbrain tissues. In females kiss1 mRNA levels gradually increased during the first few weeks of life to peak in fish with ovaries containing mature oocytes, while in males kiss1 mRNA levels peaked after 6 wk postfertilization when the testes exhibited initial stages of spermatogenesis and decreased after puberty. Zebrafish kiss1ra and kiss1rb were expressed differentially with similar patterns in both genders. These results indicate that the Kiss1/Kiss1r system may participate in puberty initiation in fish as well. Like human KISS1R, Kiss1ra transduces its activity via the PKC pathway, whereas Kiss1rb does so via both PKC and PKA pathways. The human KISS1R was highly activated by both huKISS10amide and zfKISS10amide, whereas both zebrafish Kiss1 receptor types were less sensitive to amidation.

Extracellular loop 3 (ECL3) and ECL3-proximal transmembrane domains VI and VII of the mesotocin and vasotocin receptors confer differential ligand selectivity and signaling activity

Mesotocin (MT) and vasotocin (VT) are the nonmammalian orthologs of mammalian oxytocin (OT) and arginine vasopressin (AVP), respectively. The OT/AVP family of peptides has arisen from gene duplication but has evolved to possess high selectivity toward their cognate receptors. The process of molecular evolution of receptors to confer high selectivity to their cognate ligands, however, is poorly understood. We constructed a series of reciprocal chimeras using a pair of bullfrog MT receptor (MTR) and VT1 receptor (VT1R) DNA fragments. Among the MTR/VT1R chimeras, the MTR chimera containing a region from transmembrane domain (TMD) VI to the carboxyl-terminal tail (C-tail) of VT1R showed an increased sensitivity to VT, while a chimeric VT1R containing TMD VI to C-tail of MTR showed an increased sensitivity to MT. Further dissection of domains using additional chimeras demonstrated that the receptor with the fragment containing extracellular loop 3 (ECL3) and ECL3-proximal TMDs VI and VII of MTR increased MT selectivity. This fragment is also important for receptor conformation that permits the signaling ability of the receptor. Particularly, the amino acids Val/Ile(6.54) in TMD VI and Pro/Glu(7.29) in ECL3 appear to be involved in this activity, since double mutation of these amino acids completely blocked signaling activity while maintaining ligand binding activity. Mutations at these residues in human OT and AVP 1a receptors markedly decreased receptor signaling activity. This study provides clues for understanding molecular coevolution of the OT/AVP peptides and their receptors with regard to receptor-ligand binding and receptor signaling activity.