The pelvic organs receive no parasympathetic innervation

The pelvic organs (bladder, rectum, and sex organs) have been represented for a century as receiving autonomic innervation from two pathways - lumbar sympathetic and sacral parasympathetic - by way of a shared relay, the pelvic ganglion, conceived as an assemblage of sympathetic and parasympathetic neurons. Using single-cell RNA sequencing, we find that the mouse pelvic ganglion is made of four classes of neurons, distinct from both sympathetic and parasympathetic ones, albeit with a kinship to the former, but not the latter, through a complex genetic signature. We also show that spinal lumbar preganglionic neurons synapse in the pelvic ganglion onto equal numbers of noradrenergic and cholinergic cells, both of which therefore serve as sympathetic relays. Thus, the pelvic viscera receive no innervation from parasympathetic or typical sympathetic neurons, but instead from a divergent tail end of the sympathetic chains, in charge of its idiosyncratic functions.

System-wide mapping of peptide-GPCR interactions in C. elegans

Neuropeptides and peptide hormones are ancient, widespread signaling molecules that underpin almost all brain functions. They constitute a broad ligand-receptor network, mainly by binding to G protein-coupled receptors (GPCRs). However, the organization of the peptidergic network and roles of many peptides remain elusive, as our insight into peptide-receptor interactions is limited and many peptide GPCRs are still orphan receptors. Here we report a genome-wide peptide-GPCR interaction map in Caenorhabditis elegans. By reverse pharmacology screening of over 55,384 possible interactions, we identify 461 cognate peptide-GPCR couples that uncover a broad signaling network with specific and complex combinatorial interactions encoded across and within single peptidergic genes. These interactions provide insights into peptide functions and evolution. Combining our dataset with phylogenetic analysis supports peptide-receptor co-evolution and conservation of at least 14 bilaterian peptidergic systems in C. elegans. This resource lays a foundation for system-wide analysis of the peptidergic network.

Receptor deorphanization in an echinoderm reveals kisspeptin evolution and relationship with SALMFamide neuropeptides

BACKGROUND

Kisspeptins are neuropeptides that regulate reproductive maturation in mammals via G-protein-coupled receptor-mediated stimulation of gonadotropin-releasing hormone secretion from the hypothalamus. Phylogenetic analysis of kisspeptin-type receptors indicates that this neuropeptide signaling system originated in a common ancestor of the Bilateria, but little is known about kisspeptin signaling in invertebrates.

RESULTS

Contrasting with the occurrence of a single kisspeptin receptor in mammalian species, here, we report the discovery of an expanded family of eleven kisspeptin-type receptors in a deuterostome invertebrate - the starfish Asterias rubens (phylum Echinodermata). Furthermore, neuropeptides derived from four precursor proteins were identified as ligands for six of these receptors. One or more kisspeptin-like neuropeptides derived from two precursor proteins (ArKPP1, ArKPP2) act as ligands for four A. rubens kisspeptin-type receptors (ArKPR1,3,8,9). Furthermore, a family of neuropeptides that act as muscle relaxants in echinoderms (SALMFamides) are ligands for two A. rubens kisspeptin-type receptors (ArKPR6,7). The SALMFamide neuropeptide S1 (or ArS1.4) and a 'cocktail' of the seven neuropeptides derived from the S1 precursor protein (ArS1.1-ArS1.7) act as ligands for ArKPR7. The SALMFamide neuropeptide S2 (or ArS2.3) and a 'cocktail' of the eight neuropeptides derived from the S2 precursor protein (ArS2.1-ArS2.8) act as ligands for ArKPR6.

CONCLUSIONS

Our findings reveal a remarkable diversity of neuropeptides that act as ligands for kisspeptin-type receptors in starfish and provide important new insights into the evolution of kisspeptin signaling. Furthermore, the discovery of the hitherto unknown relationship of kisspeptins with SALMFamides, neuropeptides that were discovered in starfish prior to the identification of kisspeptins in mammals, presents a radical change in perspective for research on kisspeptin signaling.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

ETS Proteins Bind with Glucocorticoid Receptors: Relevance for Treatment of Ewing Sarcoma

The glucocorticoid receptor (GR) acts as a ubiquitous cortisol-dependent transcription factor (TF). To identify co-factors, we used protein-fragment complementation assays and found that GR recognizes FLI1 and additional ETS family proteins, TFs relaying proliferation and/or migration signals. Following steroid-dependent translocation of FLI1 and GR to the nucleus, the FLI1-specific domain (FLS) binds with GR and strongly enhances GR's transcriptional activity. This interaction has functional consequences in Ewing sarcoma (ES), childhood and adolescence bone malignancies driven by fusions between EWSR1 and FLI1. In vitro, GR knockdown inhibited the migration and proliferation of ES cells, and in animal models, antagonizing GR (or lowering cortisol) retarded both tumor growth and metastasis from bone to lung. Taken together, our findings offer mechanistic rationale for repurposing GR-targeting drugs for the treatment of patients with ES.

Low-frequency variation near common germline susceptibility loci are associated with risk of Ewing sarcoma

BACKGROUND

Ewing sarcoma (EwS) is a rare, aggressive solid tumor of childhood, adolescence and young adulthood associated with pathognomonic EWSR1-ETS fusion oncoproteins altering transcriptional regulation. Genome-wide association studies (GWAS) have identified 6 common germline susceptibility loci but have not investigated low-frequency inherited variants with minor allele frequencies below 5% due to limited genotyped cases of this rare tumor.

METHODS

We investigated the contribution of rare and low-frequency variation to EwS susceptibility in the largest EwS genome-wide association study to date (733 EwS cases and 1,346 unaffected controls of European ancestry).

RESULTS

We identified two low-frequency variants, rs112837127 and rs2296730, on chromosome 20 that were associated with EwS risk (OR = 0.186 and 2.038, respectively; P-value < 5×10-8) and located near previously reported common susceptibility loci. After adjusting for the most associated common variant at the locus, only rs112837127 remained a statistically significant independent signal (OR = 0.200, P-value = 5.84×10-8).

CONCLUSIONS

These findings suggest rare variation residing on common haplotypes are important contributors to EwS risk.

IMPACT

Motivate future targeted sequencing studies for a comprehensive evaluation of low-frequency and rare variation around common EwS susceptibility loci.

NPY/NPF-Related Neuropeptide FLP-34 Signals from Serotonergic Neurons to Modulate Aversive Olfactory Learning in Caenorhabditis elegans

Aversive learning is fundamental for animals to increase chances of survival. In addition to classical neurotransmitters, neuropeptides have emerged to modulate such complex behaviors. Among them, neuropeptide Y (NPY) is well known to promote aversive memory acquisition in mammals. Here we identify an NPY/neuropeptide F (NPF)-related neuropeptide system in Caenorhabditis elegans and show that this FLP-34/NPR-11 system is required for learning negative associations, a process that is reminiscent of NPY signaling in mammals. The Caenorhabditis elegans NPY/NPF ortholog FLP-34 displays conserved structural hallmarks of bilaterian-wide NPY/NPF neuropeptides. We show that it is required for aversive olfactory learning after pairing diacetyl with the absence of food, but not for appetitive olfactory learning in response to butanone. To mediate diacetyl learning and thus integrate the aversive food context with the diacetyl odor, FLP-34 is released from serotonergic neurons and signals through its evolutionarily conserved NPY/NPF GPCR, NPR-11, in downstream AIA interneurons. NPR-11 activation in the AIA integration center results in avoidance of a previously attractive stimulus. This study opens perspectives for a deeper understanding of stress conditions in which aversive learning results in excessive avoidance.SIGNIFICANCE STATEMENT Aversive learning evolved early in evolution to promote avoidance of dangerous and stressful situations. In addition to classical neurotransmitters, neuropeptides are emerging as modulators of complex behaviors, including learning and memory. Here, we identified the evolutionary ortholog of neuropeptide Y/neuropeptide F in the nematode Caenorhabditis elegans, and we discovered that it is required for olfactory aversive learning. In addition, we elucidated the neural circuit underlying this avoidance behavior, and we discovered a novel coordinated action of Caenorhabditis elegans neuropeptide Y/neuropeptide F and serotonin that could aid in our understanding of the molecular mechanisms underlying stress disorders in which excessive avoidance results in maladaptive behaviors.

Echinoderms provide missing link in the evolution of PrRP/sNPF-type neuropeptide signalling

Neuropeptide signalling systems comprising peptide ligands and cognate receptors are evolutionarily ancient regulators of physiology and behaviour. However, there are challenges associated with determination of orthology between neuropeptides in different taxa. Orthologs of vertebrate neuropeptide-Y (NPY) known as neuropeptide-F (NPF) have been identified in protostome invertebrates, whilst prolactin-releasing peptide (PrRP) and short neuropeptide-F (sNPF) have been identified as paralogs of NPY/NPF in vertebrates and protostomes, respectively. Here we investigated the occurrence of NPY/NPF/PrRP/sNPF-related signalling systems in a deuterostome invertebrate phylum - the Echinodermata. Analysis of transcriptome/genome sequence data revealed loss of NPY/NPF-type signalling, but orthologs of PrRP-type neuropeptides and sNPF/PrRP-type receptors were identified in echinoderms. Furthermore, experimental studies revealed that the PrRP-type neuropeptide pQDRSKAMQAERTGQLRRLNPRF-NH₂ is a potent ligand for a sNPF/PrRP-type receptor in the starfish Asterias rubens. Our findings indicate that PrRP-type and sNPF-type signalling systems are orthologous and originated as a paralog of NPY/NPF-type signalling in Urbilateria.

RPamide neuropeptides NLP-22 and NLP-2 act through GnRH-like receptors to promote sleep and wakefulness in C. elegans

Sleep and wakefulness are fundamental behavioral states of which the underlying molecular principles are becoming slowly elucidated. Transitions between these states require the coordination of multiple neurochemical and modulatory systems. In Caenorhabditis elegans sleep occurs during a larval transition stage called lethargus and is induced by somnogenic neuropeptides. Here, we identify two opposing neuropeptide/receptor signaling pathways: NLP-22 promotes behavioral quiescence, whereas NLP-2 promotes movement during lethargus, by signaling through gonadotropin-releasing hormone (GnRH) related receptors. Both NLP-2 and NLP-22 belong to the RPamide neuropeptide family and share sequence similarities with neuropeptides of the bilaterian GnRH, adipokinetic hormone (AKH) and corazonin family. RPamide neuropeptides dose-dependently activate the GnRH/AKH-like receptors GNRR-3 and GNRR-6 in a cellular receptor activation assay. In addition, nlp-22-induced locomotion quiescence requires the receptor gnrr-6. By contrast, wakefulness induced by nlp-2 overexpression is diminished by deletion of either gnrr-3 or gnrr-6. nlp-2 is expressed in a pair of olfactory AWA neurons and cycles with larval periodicity, as reported for nlp-22, which is expressed in RIA. Our data suggest that the somnogenic NLP-22 neuropeptide signals through GNRR-6, and that both GNRR-3 and GNRR-6 are required for the wake-promoting action of NLP-2 neuropeptides.

Cooperation of cancer drivers with regulatory germline variants shapes clinical outcomes

Pediatric malignancies including Ewing sarcoma (EwS) feature a paucity of somatic alterations except for pathognomonic driver-mutations that cannot explain overt variations in clinical outcome. Here, we demonstrate in EwS how cooperation of dominant oncogenes and regulatory germline variants determine tumor growth, patient survival and drug response. Binding of the oncogenic EWSR1-FLI1 fusion transcription factor to a polymorphic enhancer-like DNA element controls expression of the transcription factor MYBL2 mediating these phenotypes. Whole-genome and RNA sequencing reveals that variability at this locus is inherited via the germline and is associated with variable inter-tumoral MYBL2 expression. High MYBL2 levels sensitize EwS cells for inhibition of its upstream activating kinase CDK2 in vitro and in vivo, suggesting MYBL2 as a putative biomarker for anti-CDK2-therapy. Collectively, we establish cooperation of somatic mutations and regulatory germline variants as a major determinant of tumor progression and highlight the importance of integrating the regulatory genome in precision medicine.

ALK mutation dynamics and clonal evolution in a neuroblastoma model exhibiting two ALK mutations

The ALK gene is a major oncogene of neuroblastoma cases exhibiting ALK activating mutations. Here, we characterized two neuroblastoma cell lines established from a stage 4 patient at diagnosis either from the primary tumor (PT) or from the bone marrow (BM). Both cell lines exhibited similar genomic profiles. All cells in the BM-derived cell line exhibited an ALK F1174L mutation, whereas this mutation was present in only 5% of the cells in the earliest passages of the PT-derived cell line. The BM-derived cell line presented with a higher proliferation rate in vitro and injections in Nude mice resulted in tumor formation only for the BM-derived cell line. Next, we observed that the F1174L mutation frequency in the PT-derived cell line increased with successive passages. Further Whole Exome Sequencing revealed a second ALK mutation, L1196M, in this cell line. Digital droplet PCR documented that the allele fractions of both mutations changed upon passages, and that the F1174L mutation reached 50% in late passages, indicating clonal evolution. In vitro treatment of the PT-derived cell line exhibiting the F1174L and L1196M mutations with the alectinib inhibitor resulted in an enrichment of the L1196M mutation. Using xenografts, we documented a better efficacy of alectinib compared to crizotinib on tumor growth and an enrichment of the L1196M mutation at the end of both treatments. Finally, single-cell RNA-seq analysis was consistent with both mutations resulting in ALK activation. Altogether, this study provides novel insights into ALK mutation dynamics in a neuroblastoma model harbouring two ALK mutations.

Abstract A13: Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility

Ewing sarcoma (EWS), a pediatric tumor predominantly occurring in children of European ancestry, is characterized by the EWSR1-FLI1 fusion oncogene. We performed a genome-wide association study (GWAS) of 749 EWS cases and 1,378 unaffected individuals of European ancestry. Our study replicated previously reported susceptibility loci at 1p36.22, 10q21.3, and 15q15.1 and identified new loci at 6p25.1, 8q24.23, 20p11.22, and 20p11.23 (P-values &lt;5x10-8). Effect estimates exhibited odds ratios (ORs) in excess of 1.7, which is high for cancer GWAS, and striking in light of the rarity of EWS cases in familial cancer syndromes. In expression quantitative trait locus (eQTL) analyses, we observed altered expression patterns for plausible candidate genes at 6p25.1 with RREB1, a RAS-responsive element, and at 20p11.23 with KIZ, a centrosomal stabilization protein. The 20p11.22 locus is also near NKX2-2, a highly overexpressed gene in EWS. Interestingly, most loci reside near GGAA repeat sequences where binding of the EWSR1-FLI1 fusion protein occurs. The seven EWS susceptibility loci discovered in only 749 cases make EWS one of the most productive GWAS studied cancers when considering a locus-to-case discovery ratio. The substantive estimated ORs suggest a distinctive underlying genetic architecture for EWS in which moderate-risk SNPs likely constitute a significant fraction and may interact with EWSR1-FLI1 binding.

Citation Format: Mitchell J. Machiela, Thomas G.P. Grünewald, Didier Surdez, Stephanie Reynaud, Olivier Mirabeau, Eric Karlins, Rebeca Alba Rubio, Sakina Zaidi, Sandrine Grossetete-Lalami, Stelly Ballet, Eve Lapouble, Valérie Laurence, Jean Michon, Gaelle Pierron, Heinrich Kovar, Nathalie Gaspar, Udo Kontny, Anna González-Neira, Piero Picci, Javier Alonso, Ana Patino-Garcia, Nadège Corradini, Neal D. Freedman, Nathaniel Rothman, Casey L. Dagnall, Laurie Burdett, Kristine Jones, Michelle Manning, Kathleen Wyatt, Weiyin Zhou, Meredith Yeager, David G. Cox, Robert N. Hoover, Javed Khan, Gregory T. Armstrong, Wendy M. Leisenring, Smita Bhatia, Leslie L. Robison, Uta Dirksen, Markus Metzler, Wolfgang Hartmann, Konstantin Strauch, Thomas Kirchner, Andreas E. Kulozik, Lindsay M. Morton, Lisa Mirabello, Margaret A. Tucker, Franck Tirode, Stephen Chanock, Olivier Delattre. Genome-wide association study identifies multiple new loci associated with Ewing sarcoma susceptibility [abstract]. In: Proceedings of the AACR Special Conference: Pediatric Cancer Research: From Basic Science to the Clinic; 2017 Dec 3-6; Atlanta, Georgia. Philadelphia (PA): AACR; Cancer Res 2018;78(19 Suppl):Abstract nr A13.

Abstract 2970: Multiple new susceptibility loci identified in genome-wide association study of Ewing sarcoma

Ewing sarcoma (EWS) is a rare pediatric tumor predominantly occurring in children of European ancestry and is characterized by the pathognomonic EWSR1-FLI1 fusion oncogene. To identify germline susceptibility loci associated with EWS risk, we performed a genome-wide association study (GWAS) meta-analysis of 749 EWS cases and 1,378 unaffected individuals of European ancestry from sample collections within the Institut Curie, National Cancer Institute and the Childhood Cancer Survivor Study. Our study replicated previously reported susceptibility loci at 1p36.22, 10q21.3 and 15q15.1 as well as identified novel loci at 6p25.1, 8q24.23, 20p11.22 and 20p11.23 (P-values&lt;5×10-8). These seven EWS susceptibility loci discovered in only 749 cases make EWS one of the most productive GWAS studied cancers when considering a locus to case discovery ratio. All estimated effect estimates were high for cancer GWAS with odds ratios in excess of 1.7 observed. These high per allele effects among relatively common germline variants are striking in light of the rarity of EWS cases and lack of evidence of EWS as part of a familial cancer syndrome and therefore suggest a distinctive genetic architecture for EWS. Interestingly, in silico bioinformatics analysis identified that most EWS susceptibility loci reside near GGAA nucleotide repeat sequences where binding of the EWSR1-FLI1 fusion protein occurs. ChIP-seq analyses confirmed in vivo binding of EWSR1-FLI1, suggesting germline variation in these regions could alter EWSR1-FLI1 binding and potentially deregulate neighboring genes. To identify genes with allele specific expression differences, we carried out expression quantitative trait locus (eQTL) analyses at each identified EWS susceptibility locus. We identified eQTLs for plausible candidate genes at 6p25.1 with RREB1, a RAS-responsive element, and at 20p11.23 with KIZ, a centrosomal stabilization protein. We also noted the 20p11.22 locus is near NKX2-2, a highly overexpressed gene in EWS, although no eQTL was observed in our expression data. Furthermore, knockdown of EWSR1-FLI1 in EWS cell lines indicated a more than 2-fold difference in expression of RREB1 and NKX2-2, further supporting the role of specific regulation of these genes by EWSR1-FLI1 and suggesting RREB1 and NKX2-2 may be transcription factors involved in core regulatory circuitries of EWS. Overall, our study suggests a distinctive underlying genetic architecture for EWS in which moderate risk common germline variants interact with EWSR1-FLI1 binding to alter expression of nearby target genes.

Citation Format: Mitchell J. Machiela, Thomas G. Grünewald, Didier Surdez, Stephanie Reynaud, Olivier Mirabeau, Eric Karlins, Rebeca Alba Rubio, Sakina Zaidi, Sandrine Grossetete-Lalami, Stelly Ballet, Eve Lapouble, Valérie Laurence, Jean Michon, Gaelle Pierron, Heinrich Kovar, Nathalie Gaspar, Udo Kontny, Anna González-Neira, Piero Picci, Javier Alonso, Ana Patino-Garcia, Nadège Corradini, Neal D. Freedman, Nathaniel Rothman, Casey L. Dagnall, Laurie Burdette, Kristine Jones, Michelle Manning, Kathleen Wyatt, Weiyin Zhou, Meredith Yeager, David G. Cox, Robert N. Hoover, Javed Khan, Gregory T. Armstrong, Wendy M. Leisenring, Smita Bhatia, Leslie L. Robison, Uta Dirksen, Markus Metzler, Wolfgang Hartmann, Konstantin Strauch, Thomas Kirchner, Andreas E. Kulozik, Lindsay M. Morton, Lisa Mirabello, Margaret A. Tucker, Franck Tirode, Stephen J. Chanock, Olivier Delattre. Multiple new susceptibility loci identified in genome-wide association study of Ewing sarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2970.

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.

Summary: A review of the revolutionary advances in our knowledge of the evolution of neuropeptide signalling systems that have been enabled by comparative genomics and neuropeptide receptor deorphanisation.

Synaptic and peptidergic connectome of a neurosecretory center in the annelid brain

Neurosecretory centers in animal brains use peptidergic signaling to influence physiology and behavior. Understanding neurosecretory center function requires mapping cell types, synapses, and peptidergic networks. Here we use transmission electron microscopy and gene expression mapping to analyze the synaptic and peptidergic connectome of an entire neurosecretory center. We reconstructed 78 neurosecretory neurons and mapped their synaptic connectivity in the brain of larval Platynereis dumerilii, a marine annelid. These neurons form an anterior neurosecretory center expressing many neuropeptides, including hypothalamic peptide orthologs and their receptors. Analysis of peptide-receptor pairs in spatially mapped single-cell transcriptome data revealed sparsely connected networks linking specific neuronal subsets. We experimentally analyzed one peptide-receptor pair and found that a neuropeptide can couple neurosecretory and synaptic brain signaling. Our study uncovered extensive networks of peptidergic signaling within a neurosecretory center and its connection to the synaptic brain.

Two genomes of highly polyphagous lepidopteran pests (Spodoptera frugiperda, Noctuidae) with different host-plant ranges

Emergence of polyphagous herbivorous insects entails significant adaptation to recognize, detoxify and digest a variety of host-plants. Despite of its biological and practical importance - since insects eat 20% of crops - no exhaustive analysis of gene repertoires required for adaptations in generalist insect herbivores has previously been performed. The noctuid moth Spodoptera frugiperda ranks as one of the world’s worst agricultural pests. This insect is polyphagous while the majority of other lepidopteran herbivores are specialist. It consists of two morphologically indistinguishable strains (“C” and “R”) that have different host plant ranges. To describe the evolutionary mechanisms that both enable the emergence of polyphagous herbivory and lead to the shift in the host preference, we analyzed whole genome sequences from laboratory and natural populations of both strains. We observed huge expansions of genes associated with chemosensation and detoxification compared with specialist Lepidoptera. These expansions are largely due to tandem duplication, a possible adaptation mechanism enabling polyphagy. Individuals from natural C and R populations show significant genomic differentiation. We found signatures of positive selection in genes involved in chemoreception, detoxification and digestion, and copy number variation in the two latter gene families, suggesting an adaptive role for structural variation.

Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study

IMPORTANCE

The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation.

OBJECTIVE

To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases.

DATA SOURCES

Genomewide association studies (GWAS) published up to January 15, 2015.

STUDY SELECTION

GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available.

DATA EXTRACTION AND SYNTHESIS

Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population.

MAIN OUTCOMES AND MEASURES

Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation.

RESULTS

Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]).

CONCLUSIONS AND RELEVANCE

It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

Transcriptomic identification of starfish neuropeptide precursors yields new insights into neuropeptide evolution

Neuropeptides are evolutionarily ancient mediators of neuronal signalling in nervous systems. With recent advances in genomics/transcriptomics, an increasingly wide range of species has become accessible for molecular analysis. The deuterostomian invertebrates are of particular interest in this regard because they occupy an ‘intermediate' position in animal phylogeny, bridging the gap between the well-studied model protostomian invertebrates (e.g. Drosophila melanogaster, Caenorhabditis elegans) and the vertebrates. Here we have identified 40 neuropeptide precursors in the starfish Asterias rubens, a deuterostomian invertebrate from the phylum Echinodermata. Importantly, these include kisspeptin-type and melanin-concentrating hormone-type precursors, which are the first to be discovered in a non-chordate species. Starfish tachykinin-type, somatostatin-type, pigment-dispersing factor-type and corticotropin-releasing hormone-type precursors are the first to be discovered in the echinoderm/ambulacrarian clade of the animal kingdom. Other precursors identified include vasopressin/oxytocin-type, gonadotropin-releasing hormone-type, thyrotropin-releasing hormone-type, calcitonin-type, cholecystokinin/gastrin-type, orexin-type, luqin-type, pedal peptide/orcokinin-type, glycoprotein hormone-type, bursicon-type, relaxin-type and insulin-like growth factor-type precursors. This is the most comprehensive identification of neuropeptide precursor proteins in an echinoderm to date, yielding new insights into the evolution of neuropeptide signalling systems. Furthermore, these data provide a basis for experimental analysis of neuropeptide function in the unique context of the decentralized, pentaradial echinoderm bauplan.

Chimeric EWSR1-FLI1 regulates the Ewing sarcoma susceptibility gene EGR2 via a GGAA microsatellite

Deciphering the ways in which somatic mutations and germline susceptibility variants cooperate to promote cancer is challenging. Ewing sarcoma is characterized by fusions between EWSR1 and members of the ETS gene family, usually EWSR1-FLI1, leading to the generation of oncogenic transcription factors that bind DNA at GGAA motifs. A recent genome-wide association study identified susceptibility variants near EGR2. Here we found that EGR2 knockdown inhibited proliferation, clonogenicity and spheroidal growth in vitro and induced regression of Ewing sarcoma xenografts. Targeted germline deep sequencing of the EGR2 locus in affected subjects and controls identified 291 Ewing-associated SNPs. At rs79965208, the A risk allele connected adjacent GGAA repeats by converting an interspaced GGAT motif into a GGAA motif, thereby increasing the number of consecutive GGAA motifs and thus the EWSR1-FLI1-dependent enhancer activity of this sequence, with epigenetic characteristics of an active regulatory element. EWSR1-FLI1 preferentially bound to the A risk allele, which increased global and allele-specific EGR2 expression. Collectively, our findings establish cooperation between a dominant oncogene and a susceptibility variant that regulates a major driver of Ewing sarcomagenesis.

The Evolution and Variety of RFamide-Type Neuropeptides: Insights from Deuterostomian Invertebrates

Five families of neuropeptides that have a C-terminal RFamide motif have been identified in vertebrates: (1) gonadotropin-inhibitory hormone (GnIH), (2) neuropeptide FF (NPFF), (3) pyroglutamylated RFamide peptide (QRFP), (4) prolactin-releasing peptide (PrRP), and (5) Kisspeptin. Experimental demonstration of neuropeptide-receptor pairings combined with comprehensive analysis of genomic and/or transcriptomic sequence data indicate that, with the exception of the deuterostomian PrRP system, the evolutionary origins of these neuropeptides can be traced back to the common ancestor of bilaterians. Here, we review the occurrence of homologs of vertebrate RFamide-type neuropeptides and their receptors in deuterostomian invertebrates - urochordates, cephalochordates, hemichordates, and echinoderms. Extending analysis of the occurrence of the RFamide motif in other bilaterian neuropeptide families reveals RFamide-type peptides that have acquired modified C-terminal characteristics in the vertebrate lineage (e.g., NPY/NPF), neuropeptide families where the RFamide motif is unique to protostomian members (e.g., CCK/sulfakinins), and RFamide-type peptides that have been lost in the vertebrate lineage (e.g., luqins). Furthermore, the RFamide motif is also a feature of neuropeptide families with a more restricted phylogenetic distribution (e.g., the prototypical FMRFamide-related neuropeptides in protostomes). Thus, the RFamide motif is both an ancient and a convergent feature of neuropeptides, with conservation, acquisition, or loss of this motif occurring in different branches of the animal kingdom.

Molecular evolution of peptidergic signaling systems in bilaterians

Peptide hormones and their receptors are widespread in metazoans, but the knowledge we have of their evolutionary relationships remains unclear. Recently, accumulating genome sequences from many different species have offered the opportunity to reassess the relationships between protostomian and deuterostomian peptidergic systems (PSs). Here we used sequences of all human rhodopsin and secretin-type G protein-coupled receptors as bait to retrieve potential homologs in the genomes of 15 bilaterian species, including nonchordate deuterostomian and lophotrochozoan species. Our phylogenetic analysis of these receptors revealed 29 well-supported subtrees containing mixed sets of protostomian and deuterostomian sequences. This indicated that many vertebrate and arthropod PSs that were previously thought to be phyla specific are in fact of bilaterian origin. By screening sequence databases for potential peptides, we then reconstructed entire bilaterian peptide families and showed that protostomian and deuterostomian peptides that are ligands of orthologous receptors displayed some similarity at the level of their primary sequence, suggesting an ancient coevolution between peptide and receptor genes. In addition to shedding light on the function of human G protein-coupled receptor PSs, this work presents orthology markers to study ancestral neuron types that were probably present in the last common bilaterian ancestor.

Evolution of dopamine receptor genes of the D1 class in vertebrates

The receptors of the dopamine neurotransmitter belong to two unrelated classes named D₁ and D₂. For the D₁ receptor class, only two subtypes are found in mammals, the D_1A and D_1B, receptors, whereas additional subtypes, named D_1C, D_1D, and D_1X, have been found in other vertebrate species. Here, we analyzed molecular phylogeny, gene synteny, and gene expression pattern of the D₁ receptor subtypes in a large range of vertebrate species, which leads us to propose a new view of the evolution of D₁ dopamine receptor genes. First, we show that D_1C and D_1D receptor sequences are encoded by orthologous genes. Second, the previously identified Cypriniform D_1X sequence is a teleost-specific paralog of the D_1B sequences found in all groups of jawed vertebrates. Third, zebrafish and several sauropsid species possess an additional D₁-like gene, which is likely to form another orthology group of vertebrate ancestral genes, which we propose to name D_1E. Ancestral jawed vertebrates are thus likely to have possessed four classes of D₁ receptor genes—D_1A, D_1B(X), D_1C(D), and D_1E—which arose from large-scale gene duplications. The D_1C receptor gene would have been secondarily lost in the mammalian lineage, whereas the D_1E receptor gene would have been lost independently in several lineages of modern vertebrates. The D_1A receptors are well conserved throughout jawed vertebrates, whereas sauropsid D_1C receptors have rapidly diverged, to the point that they were misidentified as D_1D. The functional significance of the D_1C receptor loss is not known. It is possible that the function may have been substituted with D_1A or D_1B receptors in mammals, following the disappearance of D_1C receptors in these species.

Identification of novel peptide hormones in the human proteome by hidden Markov model screening

Peptide hormones are small, processed, and secreted peptides that signal via membrane receptors and play critical roles in normal and pathological physiology. The search for novel peptide hormones has been hampered by their small size, low or restricted expression, and lack of sequence similarity. To overcome these difficulties, we developed a bioinformatics search tool based on the hidden Markov model formalism that uses several peptide hormone sequence features to estimate the likelihood that a protein contains a processed and secreted peptide of this class. Application of this tool to an alignment of mammalian proteomes ranked 90% of known peptide hormones among the top 300 proteins. An analysis of the top scoring hypothetical and poorly annotated human proteins identified two novel candidate peptide hormones. Biochemical analysis of the two candidates, which we called spexin and augurin, showed that both were localized to secretory granules in a transfected pancreatic cell line and were recovered from the cell supernatant. Spexin was expressed in the submucosal layer of the mouse esophagus and stomach, and a predicted peptide from the spexin precursor induced muscle contraction in a rat stomach explant assay. Augurin was specifically expressed in mouse endocrine tissues, including pituitary and adrenal gland, choroid plexus, and the atrio-ventricular node of the heart. Our findings demonstrate the utility of a bioinformatics approach to identify novel biologically active peptides. Peptide hormones and their receptors are important diagnostic and therapeutic targets, and our results suggest that spexin and augurin are novel peptide hormones likely to be involved in physiological homeostasis.