Identification of a new Kir6 potassium channel and comparison of properties of Kir6 subtypes by structural modelling and molecular dynamics.

ATP-sensitive potassium ion channels (KATP) are transmembrane proteins that modulate insulin release and muscle contraction. KATP channels are composed of two types of subunit, Kir6 and SUR, which exist in two and three isoforms respectively with different tissue distribution. In this work, we identify a previously undescribed ancestral vertebrate gene encoding a Kir6-related protein that we have named Kir6.3, which may not have a SUR binding partner, unlike the other two Kir6 proteins. Whereas Kir6.3 was lost in amniotes including mammals, it is still present in several early-diverging vertebrate lineages such as frogs, coelacanth, and rayfinned fishes. Molecular dynamics (MD) simulations using homology models of Kir6.1, Kir6.2, and Kir6.3 from the coelacanth Latimeria chalumnae showed that the three proteins exhibit subtle differences in their dynamics. Steered MD simulations of Kir6-SUR pairs suggest that Kir6.3 has a lower binding affinity for the SUR proteins than either Kir6.1 or Kir6.2. As we found no additional SUR gene in the genomes of the species that have Kir6.3, it most likely forms a lone tetramer. These findings invite studies of the tissue distribution of Kir6.3 in relation to the other Kir6 as well as SUR proteins to determine the functional roles of Kir6.3.

Ancient multiplicity in cyclic nucleotide-gated (CNG) cation channel repertoire was reduced in the ancestor of Olfactores before re-expansion by whole genome duplications in vertebrates.

Cyclic nucleotide-gated (CNG) cation channels are important heterotetrameric proteins in the retina, with different subunit composition in cone and rod photoreceptor cells: three CNGA3 and one CNGB3 in cones and three CNGA1 and one CNGB1 in rods. CNGA and CNGB subunits form separate subfamilies. We have analyzed the evolution of the CNG gene family in metazoans, with special focus on vertebrates by using sequence-based phylogeny and conservation of chromosomal synteny to deduce paralogons resulting from the early vertebrate whole genome duplications (WGDs). Our analyses show, unexpectedly, that the CNGA subfamily had four sister subfamilies in the ancestor of bilaterians and cnidarians that we named CNGC, CNGD, CNGE and CNGF. Of these, CNGC, CNGE and CNGF were lost in the ancestor of Olfactores while CNGD was lost in the vertebrate ancestor. The remaining CNGA and CNGB genes were expanded by a local duplication of CNGA and the subsequent chromosome duplications in the basal vertebrate WGD events. Upon some losses, this resulted in the gnathostome ancestor having three members in the visual CNGA subfamily (CNGA1-3), a single CNGA4 gene, and two members in the CNGB subfamily (CNGB1 and CNGB3). The nature of chromosomal rearrangements in the vertebrate CNGA paralogon was resolved by including the genomes of a non-teleost actinopterygian and an elasmobranch. After the teleost-specific WGD, additional duplicates were generated and retained for CNGA1, CNGA2, CNGA3 and CNGB1. Furthermore, teleosts retain a local duplicate of CNGB3. The retention of duplicated CNG genes is explained by their subfunctionalisation and photoreceptor-specific expression. In conclusion, this study provides evidence for four previously unknown CNG subfamilies in metazoans and further evidence that the early vertebrate WGD events were instrumental in the evolution of the vertebrate visual and central nervous systems.

Dynamic evolution of transient receptor potential vanilloid (TRPV) ion channel family with numerous gene duplications and losses.

The transient receptor potential vanilloid (TRPV) ion channel family is involved in multiple sensory and physiological functions including thermosensing and temperature-dependent neuroendocrine regulation. The objective of the present study was to investigate the number, origin and evolution of TRPV genes in metazoans, with special focus on the impact of the vertebrate whole-genome duplications (WGD). Gene searches followed by phylogenetic and synteny analyses revealed multiple previously undescribed TRPV genes. The common ancestor of Cnidaria and Bilateria had three TRPV genes that became four in the deuterostome ancestor. Two of these were lost in the vertebrate ancestor. The remaining two genes gave rise to two TRPV subfamilies in vertebrates, consisting of subtypes 1, 2, 3, 4, 9 and 5, 6, 7, 8, respectively. This gene expansion resulted from the two basal vertebrate WGD events (1R and 2R) and three local duplications before the radiation of gnathostomes. TRPV1, 4 and 5 have been retained in all gnathostomes investigated, presumably reflecting important functions. TRPV7 and 8 have been lost independently in various lineages but are still retained in cyclostomes, actinistians (coelacanth), amphibians, prototherians and basal actinopterygians (Polypteridae). TRPV3 and 9 are present in extant elasmobranchs, while TRPV9 was lost in the osteichthyan ancestor and TRPV3 in the actinopterygian ancestor. The coelacanth has retained the ancestral osteichthyan repertoire of TRPV1, 3, 4, 5, 7 and 8. TRPV2 arose in the tetrapod ancestor. Duplications of TRPV5 occurred independently in various lineages, such as cyclostomes, chondrichthyans, anuran amphibians, sauropsids, mammals (where the duplicate is called TRPV6), and actinopterygians (Polypteridae and Esocidae). After the teleost-specific WGD (3R) only TRPV1 retained its duplicate, whereas TRPV4 and 5 remained as single genes. Both 3R-paralogs of TRPV1 were kept in some teleost species, while one paralog was lost in others. The salmonid-specific WGD (4R) duplicated TRPV1, 4, and 5 leading to six TRPV genes. The largest number was found in Xenopus tropicalis with no less than 15 TRPV genes. This study provides a comprehensive evolutionary scenario for the vertebrate TRPV family, revealing additional TRPV types and proposing a phylogeny-based classification of TRPV across metazoans.

Correction: Copy number of pancreatic polypeptide receptor gene NPY4R correlates with body mass index and waist circumference.

[This corrects the article DOI: 10.1371/journal.pone.0194668.].

Protein kinase C family evolution in jawed vertebrates.

Protein kinase C (PKC) was one of the first kinases identified in human cells. It is now known to constitute a family of kinases that respond to diacylglycerol, phosphatidylserine and for some family members, Ca2+. They have a plethora of different functions, such as cell cycle regulation, immune response and memory formation. In mammals, 12 PKC family members have been described, usually divided into 4 different subfamilies. We present here a comprehensive evolutionary analysis of the PKC genes in jawed vertebrates with special focus on the impact of the two tetraploidizations (1R and 2R) before the radiation of jawed vertebrates and the teleost tetraploidization (3R), as illuminated by synteny and paralogon analysis including many neighboring gene families. We conclude that the vertebrate predecessor had five PKC genes, as tunicates and lancelets still do, and that the PKC family should therefore ideally be organized into five subfamilies. The 1R and 2R events led to a total of 12 genes distributed among these five subfamilies. All 12 genes are still present in some of the major lineages of jawed vertebrates, including mammals, whereas birds and cartilaginous fishes have lost one member. The 3R event added another nine genes in teleosts, bringing the total to 21 genes. The zebrafish, a common experimental model animal, has retained 19. We have found no independent gene duplications. Thus, the genome doublings completely account for the complexity of this gene family in jawed vertebrates and have thereby had a huge impact on their evolution.

The Evolution of Oxytocin and Vasotocin Receptor Genes in Jawed Vertebrates: A Clear Case for Gene Duplications Through Ancestral Whole-Genome Duplications.

The neuronal and neuroendocrine peptides oxytocin (OT) and vasotocin (VT), including vasopressins, have six cognate receptors encoded by six receptor subtype genes in jawed vertebrates. The peptides elicit a broad range of responses that are specifically mediated by the receptor subtypes including neuronal functions regulating behavior and hormonal actions on reproduction and water/electrolyte balance. Previously, we have demonstrated that these six receptor subtype genes, which we designated VTR1A, VTR1B, OTR, VTR2A, VTR2B and VTR2C, arose from a syntenic ancestral gene pair, one VTR1/OTR ancestor and one VTR2 ancestor, through the early vertebrate whole-genome duplications (WGD) called 1R and 2R. This was supported by both phylogenetic and chromosomal conserved synteny data. More recently, other studies have focused on confounding factors, such as the OTR/VTR orthologs in cyclostomes, to question this scenario for the origin of the OTR/VTR gene family; proposing instead less parsimonious interpretations involving only one WGD followed by complex series of chromosomal or segmental duplications. Here, we have updated the phylogeny of the OTR/VTR gene family, including a larger number of vertebrate species, and revisited seven representative neighboring gene families from our previous conserved synteny analyses, adding chromosomal information from newer high-coverage genome assemblies from species that occupy key phylogenetic positions: the polypteriform fish reedfish (Erpetoichthys calabaricus), the cartilaginous fish thorny skate (Amblyraja radiata) and a more recent high-quality assembly of the Western clawed frog (Xenopus tropicalis) genome. Our analyses once again add strong support for four-fold symmetry, i.e., chromosome quadruplication in the same time window as the WGD events early in vertebrate evolution, prior to the jawed vertebrate radiation. Thus, the evolution of the OTR/VTR gene family can be most parsimoniously explained by two WGD events giving rise to the six ancestral genes, followed by differential gene losses of VTR2 genes in different lineages. We also argue for more coherence and clarity in the nomenclature of OT/VT receptors, based on the most parsimonious scenario.

Corticotropin-Releasing Hormone (CRH) Gene Family Duplications in Lampreys Correlate With Two Early Vertebrate Genome Doublings.

The ancestor of gnathostomes (jawed vertebrates) is generally considered to have undergone two rounds of whole genome duplication (WGD). The timing of these WGD events relative to the divergence of the closest relatives of the gnathostomes, the cyclostomes, has remained contentious. Lampreys and hagfishes are extant cyclostomes whose gene families can shed light on the relationship between the WGDs and the cyclostome-gnathostome divergence. Previously, we have characterized in detail the evolution of the gnathostome corticotropin-releasing hormone (CRH) family and found that its five members arose from two ancestral genes that existed before the WGDs. The two WGDs resulted, after secondary losses, in one triplet consisting of CRH1, CRH2, and UCN1, and one pair consisting of UCN2 and UCN3. All five genes exist in representatives for cartilaginous fishes, ray-finned fishes, and lobe-finned fishes. Differential losses have occurred in some lineages. We present here analyses of CRH-family members in lamprey and hagfish by comparing sequences and gene synteny with gnathostomes. We found five CRH-family genes in each of two lamprey species (Petromyzon marinus and Lethenteron camtschaticum) and two genes in a hagfish (Eptatretus burgeri). Synteny analyses show that all five lamprey CRH-family genes have similar chromosomal neighbors as the gnathostome genes. The most parsimonious explanation is that the lamprey CRH-family genes are orthologs of the five gnathostome genes and thus arose in the same chromosome duplications. This suggests that lampreys and gnathostomes share the same two WGD events and that these took place before the lamprey-gnathostome divergence.

Oxytocin Receptors Regulate Social Preference in Zebrafish.

With a strong tendency to socialise, the zebrafish is a useful model to study social behaviour, with implications for better treatments of social impairments, for instance in autism spectrum disorders. Although oxytocin is crucial for social behaviour in mammals, the importance of the fish orthologue - isotocin or zebrafish oxytocin (zOT) - for social behaviour in zebrafish is unclear. The aims of this study were firstly, to elucidate the receptor specificity of zOT and the related vasotocin or zebrafish vasopressin (zVP; the orthologue of mammalian vasopressin) and the nonpeptidergic oxytocin receptor antagonist L-368,899, and secondly to investigate if L-368,899 inhibits social preference in zebrafish. The potencies of ligands were evaluated for zOT/zVP family receptors in HEK293 cells. Adult and larval zebrafish were treated with L-368,899 or vehicle and subsequently assessed for social behaviour and anxiety (adults only). The antagonist L-368,899 specifically inhibited the two zOT receptors, but not the two zVP-1 receptors. The antagonist decreased social preference in adult and larval zebrafish. It did not affect anxiety in adults. These results indicate that endogenous zOT, and possibly zVP, is involved in social behaviour in zebrafish via either or both of the two zOT receptors, and show promise for future explorations of the anatomy and evolution of networks underlying social behaviour.

Circadian regulation of phosphodiesterase 6 genes in zebrafish differs between cones and rods: Implications for photopic and scotopic vision.

A correlation is known to exist between visual sensitivity and oscillations in red opsinand rhodopsin gene expression in zebrafish, both regulated by the clock gene. This indicates that an endogenous circadian clock regulates behavioural visual sensitivity, apart from the regulation exerted by the pineal organ. However, the specific mechanisms for cones (photopic vision) and rods (scotopic vision) are poorly understood. In this work, we performed gene expression, cosinor and immunohistochemical analyses to investigate other key genes involved in light perception, encoding the different subunits of phosphodiesterase pde6 and transducin GαT, in constant lighting conditions and compared to normal light-dark conditions. We found that cones display prominent circadian oscillations in mRNA levels for the inhibitory subunit gene pde6ha that could contribute to the regulation of photopic sensitivity by preventing overstimulation in photopic conditions. In rods, the mRNA levels of the inhibitory subunit gene pde6ga oscillate under normal conditions and dampen down in constant light but continue oscillating in constant darkness. There is an increase in total relative expression for pde6gb in constant conditions. These observations, together with previous data, suggest a complex regulation of the scotopic sensitivity involving endogenous and non-endogenous components, possibly present also in other teleost species. The GαT genes do not display mRNA oscillations and therefore may not be essential for the circadian regulation of photosensitivity. In summary, our results support different regulation for the zebrafish photopic and scotopic sensitivities and suggest circadian regulation of pde6ha as a key factor regulating photopic sensitivity, while the regulatory mechanisms in rods appear to be more complex.

The Neuropeptide Y Y2 Receptor Is Coexpressed with Nppb in Primary Afferent Neurons and Y2 Activation Reduces Histaminergic and IL-31-Induced Itch.

Itch stimuli are detected by specialized primary afferents that convey the signal to the spinal cord, but how itch transmission is regulated is still not completely known. Here, we investigated the roles of the neuropeptide Y (NPY)/Y2 receptor system on scratch behavior. The inhibitory Y2 receptor is expressed on mouse primary afferents, and intrathecal administration of the Y2 agonist peptide YY (PYY)3-36 reduced scratch episode frequency and duration induced by compound 48/80, an effect that could be reversed by intrathecal preadministration of the Y2 antagonist BIIE0246. Also, scratch episode duration induced by histamine could be reduced by PYY3-36 In contrast, scratch behavior induced by α-methyl-5HT, protease-activated receptor-2-activating peptide SLIGRL, chloroquine, topical dust mite extract, or mechanical itch induced by von Frey filaments was unaffected by stimulation of Y2 Primary afferent neurons expressing the Npy2r gene were found to coexpress itch-associated markers such as natriuretic peptide precursor b, oncostatin M receptor, and interleukin (IL) 31 receptor A. Accordingly, intrathecal PYY3-36 reduced the scratch behavior induced by IL-31. Our findings imply that the NPY/Y2 system reduces histaminergic and IL-31-associated itch through presynaptic inhibition of a subpopulation of itch-associated primary afferents. SIGNIFICANCE STATEMENT: The spinal neuropeptide Y system dampens scratching behavior induced by histaminergic compounds and interleukin 31, a cytokine involved in atopic dermatitis, through interactions with the Y2 receptor. The Y2 receptor is expressed by primary afferent neurons that are rich in itch-associated neurotransmitters and receptors such as somatostatin, natriuretic peptide precursor b, and interleukin 31 receptors.

Functional characterization in vitro of twelve naturally occurring variants of the human pancreatic polypeptide receptor NPY4R.

Obesity has become a global health problem and therefore understanding of the mechanisms regulating hunger and satiety is of utmost importance for the development of new treatment strategies. The Y4 receptor, encoded by the NPY4R gene, and its ligand pancreatic polypeptide (PP) have been reported to mediate a satiety signal. Multiple genetic studies have reported an association between NPY4R copy number and body weight. The gene also displays several SNP variants, many of which lead to amino acid differences, making it interesting to study. We have investigated the functional properties of 12 naturally occurring amino acid sequence variants of the Y4 and interpret the results in relation to sequence conservation and our structural model of the human Y4 receptor protein. Three receptor variants, Cys201ECL2Tyr, Val2716.41Leu and Asn3187.49Asp, were found to completely lose functional response, measured as inositol phosphate turnover, while retaining membrane expression. They display high sequence conservation and have important roles in the receptor structure. For two receptor variants the potency of PP was significantly decreased, Cys34NTSer (EC50 = 2.9 nM, p < .001) and Val1353.46Met (EC50 = 3.0 nM, p < .01), compared to wild-type Y4 (EC50 = 0.68 nM). Cys34 forms a disulphide bond with Cys298, linking the N-terminal part to ECL3. The Val1353.46Met variant has an amino acid replacement located in the TM3 helix, one helix turn above the highly conserved ERH motif. This position has influence on the network of residues involved in receptor activation and subsequent inactivation. Sequence conservation and the structural model are consistent with these results. The remaining seven positions had no significant effect on the receptor's functional response compared to wild-type Y4. These positions display more variation during evolution. Understanding of the interactions between the Y4 receptor and its native PP agonist and the effects of amino acid variation on its functional response will hopefully lead to future therapeutic possibilities.

Copy number determination of the gene for the human pancreatic polypeptide receptor NPY4R using read depth analysis and droplet digital PCR.

BACKGROUND: Copy number variation (CNV) plays an important role in human genetic diversity and has been associated with multiple complex disorders. Here we investigate a CNV on chromosome 10q11.22 that spans NPY4R, the gene for the appetite-regulating pancreatic polypeptide receptor Y4. This genomic region has been challenging to map due to multiple repeated elements and its precise organization has not yet been resolved. Previous studies using microarrays were interpreted to show that the most common copy number was 2 per genome. RESULTS: We have investigated 18 individuals from the 1000 Genomes project using the well-established method of read depth analysis and the new droplet digital PCR (ddPCR) method. We find that the most common copy number for NPY4R is 4. The estimated number of copies ranged from three to seven based on read depth analyses with Control-FREEC and CNVnator, and from four to seven based on ddPCR. We suggest that the difference between our results and those published previously can be explained by methodological differences such as reference gene choice, data normalization and method reliability. Three high-quality archaic human genomes (two Neanderthal and one Denisova) display four copies of the NPY4R gene indicating that a duplication occurred prior to the human-Neanderthal/Denisova split. CONCLUSIONS: We conclude that ddPCR is a sensitive and reliable method for CNV determination, that it can be used for read depth calibration in CNV studies based on already available whole-genome sequencing data, and that further investigation of NPY4R copy number variation and its consequences are necessary due to the role of Y4 receptor in food intake regulation.

Evolution of vertebrate nicotinic acetylcholine receptors.

BACKGROUND: Many physiological processes are influenced by nicotinic acetylcholine receptors (nAChR), ranging from neuromuscular and parasympathetic signaling to modulation of the reward system and long-term memory. Due to the complexity of the nAChR family and variable evolutionary rates among its members, their evolution in vertebrates has been difficult to resolve. In order to understand how and when the nAChR genes arose, we have used a broad approach of analyses combining sequence-based phylogeny, chromosomal synteny and intron positions. RESULTS: Our analyses suggest that there were ten subunit genes present in the vertebrate predecessor. The two basal vertebrate tetraploidizations (1R and 2R) then expanded this set to 19 genes. Three of these have been lost in mammals, resulting in 16 members today. None of the ten ancestral genes have kept all four copies after 2R. Following 2R, two of the ancestral genes became triplicates, five of them became pairs, and three seem to have remained single genes. One triplet consists of CHRNA7, CHRNA8 and the previously undescribed CHRNA11, of which the two latter have been lost in mammals but are still present in lizards and ray-finned fishes. The other triplet consists of CHRNB2, CHRNB4 and CHRNB5, the latter of which has also been lost in mammals. In ray-finned fish the neuromuscular subunit gene CHRNB1 underwent a local gene duplication generating CHRNB1.2. The third tetraploidization in the predecessor of teleosts (3R) expanded the repertoire to a total of 31 genes, of which 27 remain in zebrafish. These evolutionary relationships are supported by the exon-intron organization of the genes. CONCLUSION: The tetraploidizations explain all gene duplication events in vertebrates except two. This indicates that the genome doublings have had a substantial impact on the complexity of this gene family leading to a very large number of members that have existed for hundreds of millions of years.

Evolution of the Muscarinic Acetylcholine Receptors in Vertebrates.

The family of muscarinic acetylcholine receptors (mAChRs) consists of five members in mammals, encoded by the CHRM1-5 genes. The mAChRs are G-protein-coupled receptors, which can be divided into the following two subfamilies: M2 and M4 receptors coupling to Gi/o; and M1, M3, and M5 receptors coupling to Gq/11. However, despite the fundamental roles played by these receptors, their evolution in vertebrates has not yet been fully described. We have combined sequence-based phylogenetic analyses with comparisons of exon-intron organization and conserved synteny in order to deduce the evolution of the mAChR receptors. Our analyses verify the existence of two ancestral genes prior to the two vertebrate tetraploidizations (1R and 2R). After these events, one gene had duplicated, resulting in CHRM2 and CHRM4; and the other had triplicated, forming the CHRM1, CHRM3, and CHRM5 subfamily. All five genes are still present in all vertebrate groups investigated except the CHRM1 gene, which has not been identified in some of the teleosts or in chicken or any other birds. Interestingly, the third tetraploidization (3R) that took place in the teleost predecessor resulted in duplicates of all five mAChR genes of which all 10 are present in zebrafish. One of the copies of the CHRM2 and CHRM3 genes and both CHRM4 copies have gained introns in teleosts. Not a single separate (nontetraploidization) duplicate has been identified in any vertebrate species. These results clarify the evolution of the vertebrate mAChR family and reveal a doubled repertoire in zebrafish, inviting studies of gene neofunctionalization and subfunctionalization.

The Neuropeptide Y System Regulates Both Mechanical and Histaminergic Itch.

Itch is a somatosensory modality that serves to alert an organism to harmful elements removable by scratching, such as parasites and chemical irritants. Recently, ablation or silencing of neuropeptide Y (NPY)-expressing spinal interneurons was reported to selectively enhance mechanical itch, whereas chemical itch was unaffected. We examined the effect of activating the NPY/Y1 receptor system on scratch behavior in mice. We found that intrathecal administration of the Y1 agonist [Leu31,Pro34]-NPY (LP-NPY) attenuated itch behavior induced by application of 0.07 g von Frey filament in the nape of the neck compared with saline treatment, indicating that activation of the spinal NPY/Y1 system dampens mechanical itch. However, intrathecal administration of LP-NPY also attenuated chemically induced scratching provoked by intradermal application of histamine or the mast cell degranulator 48/80 (histaminergic itch), and the latter effect could be reversed by administration of the Y1 antagonist BIBO3304. Intrathecal application of the native nonselective agonist NPY also attenuated histamine or 48/80-induced scratching. Our analyses emphasize the importance of including additional quantitative parameters to characterize the full spectrum of itch behavior and show that the NPY/Y1 system dampens both mechanically and chemically induced scratching and hence is shared by the two submodalities of itch.

Copy number of pancreatic polypeptide receptor gene NPY4R correlates with body mass index and waist circumference.

Multiple genetic studies have linked copy number variation (CNV) in different genes to body mass index (BMI) and obesity. A CNV on chromosome 10q11.22 has been associated with body weight. This CNV region spans NPY4R, the gene encoding the pancreatic polypeptide receptor Y4, which has been described as a satiety-stimulating receptor. We have investigated CNV of the NPY4R gene and analysed its relationship to BMI, waist circumference and self-reported dietary intake from 558 individuals (216 men and 342 women) representing a wide BMI range. The copy number for NPY4R ranged from 2 to 8 copies (average 4.6±0.8). Rather than the expected negative correlation, we observed a positive correlation between NPY4R copy number and BMI as well as waist circumference in women (Pearson's r = 0.267, p = 2.65×10−7 and r = 0.256, p = 8×10−7, respectively). Each additional copy of NPY4R correlated with 2.6 kg/m2 increase in BMI and 5.67 cm increase in waist circumference (p = 2.8×10−5 and p = 6.2×10−5, respectively) for women. For men, there was no statistically significant correlation between CNV and BMI. Our results suggest that NPY4R genetic variation influences body weight in women, but the exact role of this receptor appears to be more complex than previously proposed.

Structural basis of ligand binding modes at the neuropeptide Y Y1 receptor.

Neuropeptide Y (NPY) receptors belong to the G-protein-coupled receptor superfamily and have important roles in food intake, anxiety and cancer biology 1,2 . The NPY-Y receptor system has emerged as one of the most complex networks with three peptide ligands (NPY, peptide YY and pancreatic polypeptide) binding to four receptors in most mammals, namely the Y1, Y2, Y4 and Y5 receptors, with different affinity and selectivity 3 . NPY is the most powerful stimulant of food intake and this effect is primarily mediated by the Y1 receptor (Y1R) 4 . A number of peptides and small-molecule compounds have been characterized as Y1R antagonists and have shown clinical potential in the treatment of obesity 4 , tumour 1 and bone loss 5 . However, their clinical usage has been hampered by low potency and selectivity, poor brain penetration ability or lack of oral bioavailability 6 . Here we report crystal structures of the human Y1R bound to the two selective antagonists UR-MK299 and BMS-193885 at 2.7 and 3.0 Å resolution, respectively. The structures combined with mutagenesis studies reveal the binding modes of Y1R to several structurally diverse antagonists and the determinants of ligand selectivity. The Y1R structure and molecular docking of the endogenous agonist NPY, together with nuclear magnetic resonance, photo-crosslinking and functional studies, provide insights into the binding behaviour of the agonist and for the first time, to our knowledge, determine the interaction of its N terminus with the receptor. These insights into Y1R can enable structure-based drug discovery that targets NPY receptors.

Evolution of neuropeptide signalling systems.

Neuropeptides are a diverse class of neuronal signalling molecules that regulate physiological processes and behaviour in animals. However, determining the relationships and evolutionary origins of the heterogeneous assemblage of neuropeptides identified in a range of phyla has presented a huge challenge for comparative physiologists. Here, we review revolutionary insights into the evolution of neuropeptide signalling that have been obtained recently through comparative analysis of genome/transcriptome sequence data and by 'deorphanisation' of neuropeptide receptors. The evolutionary origins of at least 30 neuropeptide signalling systems have been traced to the common ancestor of protostomes and deuterostomes. Furthermore, two rounds of genome duplication gave rise to an expanded repertoire of neuropeptide signalling systems in the vertebrate lineage, enabling neofunctionalisation and/or subfunctionalisation, but with lineage-specific gene loss and/or additional gene or genome duplications generating complex patterns in the phylogenetic distribution of paralogous neuropeptide signalling systems. We are entering a new era in neuropeptide research where it has become feasible to compare the physiological roles of orthologous and paralogous neuropeptides in a wide range of phyla. Moreover, the ambitious mission to reconstruct the evolution of neuropeptide function in the animal kingdom now represents a tangible challenge for the future.