Multitrait Genome Association Analysis Identifies New Susceptibility Genes

Total Page:16

File Type:pdf, Size:1020Kb

Multitrait Genome Association Analysis Identifies New Susceptibility Genes Complex traits J Med Genet: first published as 10.1136/jmedgenet-2018-105437 on 30 August 2018. Downloaded from ORIGINAL ARTICLE Multitrait genome association analysis identifies new susceptibility genes for human anthropometric variation in the GCAT cohort Iván Galván-Femenía,1 Mireia Obón-Santacana,1,2 David Piñeyro,3 Marta Guindo-Martinez,4 Xavier Duran,1 Anna Carreras,1 Raquel Pluvinet,3 Juan Velasco,1 Laia Ramos,3 Susanna Aussó,3 J M Mercader,5,6 Lluis Puig,7 Manuel Perucho,8 David Torrents,4,9 Victor Moreno,2,10 Lauro Sumoy,3 Rafael de Cid1 ► Additional material is ABSTRact INTRODUCTION published online only. To view, Background Heritability estimates have revealed an Common disorders cause 85% of deaths in the please visit the journal online 1 (http:// dx. doi. org/ 10. 1136/ important contribution of SNP variants for most common European Union (EU). The increasing incidence jmedgenet- 2018- 105437). traits; however, SNP analysis by single-trait genome- and prevalence of cancer, cardiovascular diseases, wide association studies (GWAS) has failed to uncover chronic respiratory diseases, diabetes and mental For numbered affiliations see their impact. In this study, we applied a multitrait GWAS illness represent a challenge that leads to extra costs end of article. approach to discover additional factor of the missing for the healthcare system. Moreover, as European population is getting older, this scenario will be Correspondence to heritability of human anthropometric variation. Dr Rafael de Cid, GCAT lab Methods We analysed 205 traits, including diseases heightened in the next few years. Like complex Group, Program of Predictive identified at baseline in the GCAT cohort (Genomes For traits, many common diseases are complex inher- and Personalized Medicine of Life- Cohort study of the Genomes of Catalonia) (n=4988), ited conditions with genetic and environmental Cancer (PMPPC), Germans Trias determinants. Advancing in their understanding i Pujol Research Institute (IGTP), a Mediterranean adult population-based cohort study Crta. de Can Ruti, Badalona from the south of Europe. We estimated SNP heritability requires the use of multifaceted and long-term 08916, Spain; rdecid@ igtp. cat contribution and single-trait GWAS for all traits from prospective approaches. Cohort analyses provide 15 million SNP variants. Then, we applied a multitrait- an exceptional tool for dissecting the architecture Received 24 April 2018 of complex diseases by contributing knowledge Revised 19 July 2018 related approach to study genome-wide association to Accepted 21 July 2018 anthropometric measures in a two-stage meta-analysis for evidence-based prevention, as exemplified by with the UK Biobank cohort (n=336 107). the Framingham Heart Study2 or the European Results Heritability estimates (eg, skin colour, Prospective Investigation into Cancer and Nutrition alcohol consumption, smoking habit, body mass cohort study.3 index, educational level or height) revealed an In the last decades, high performance DNA geno- http://jmg.bmj.com/ important contribution of SNP variants, ranging from typing technology has fuelled genomic research 18% to 77%. Single-trait analysis identified 1785 in large cohorts, having been the most promising SNPs with genome-wide significance threshold. line in research on the aetiology of most common From these, several previously reported single-trait diseases. Genome-wide association studies (GWAS) hits were confirmed in our sample with LINC01432 have provided valuable information for many single (p=1.9×10−9) variants associated with male baldness, conditions.4 Despite the perception of the limitations on September 24, 2021 by guest. Protected copyright. LDLR variants with hyperlipidaemia (ICD-9:272) of the GWAS analyses, efforts combining massive (p=9.4×10−10) and variants in IRF4 (p=2.8×10−57), data deriving from whole-genome sequencing at SLC45A2 (p=2.2×10−130), HERC2 (p=2.8×10−176), population scale with novel conceptual and meth- OCA2 (p=2.4×10−121) and MC1R (p=7.7×10−22) odological analysis frameworks have been set forth associated with hair, eye and skin colour, freckling, to explore the last frontier of the missing herita- tanning capacity and sun burning sensitivity and the bility issue,5 driving the field of genomic research Fitzpatrick phototype score, all highly correlated cross- on complex diseases to a new age.6Pritchard and phenotypes. Multitrait meta-analysis of anthropometric colleagues recently proposed the breakthrough idea variation validated 27 loci in a two-stage meta- of the omnigenic character of genetic architecture © Author(s) (or their 7 employer(s)) 2018. Re-use analysis with a large British ancestry cohort, six of of diseases and complex traits. They suggested that permitted under CC BY-NC. No which are newly reported here (p value threshold beyond a handful of driver genes (ie, core genes) commercial re-use. See rights <5×10−9) at ZRANB2-AS2, PIK3R1, EPHA7, MAD1L1, directly connected to an illness, the missing herita- and permissions. Published bility could be accounted for by multiple genes (ie, by BMJ. CACUL1 and MAP3K9. Conclusion Considering multiple-related genetic peripheral genes) not clustered in functional path- To cite: Galván-Femenía I, phenotypes improve associated genome signal ways, but dispersed along the genome, explaining Obón-Santacana M, the pleiotropy frequently seen in most complex Piñeyro D, et al. J Med Genet detection. These results indicate the potential value Epub ahead of print: [please of data-driven multivariate phenotyping for genetic traits. Core genes have been already outlined by the include Day Month Year]. studies in large population-based cohorts to contribute GWAS approach, but most of the possible contrib- doi:10.1136/ to knowledge of complex traits. uting genes have been disregarded based on meth- jmedgenet-2018-105437 odological issues such as p value or lower minor Galván-Femenía I, et al. J Med Genet 2018;0:1–14. doi:10.1136/jmedgenet-2018-105437 1 Complex traits J Med Genet: first published as 10.1136/jmedgenet-2018-105437 on 30 August 2018. Downloaded from allele frequency (MAF). Pathway disturbances have also been a Phenome landmark in the search for genetic associations,8 but not always Baseline variables were obtained from a self-reported epidemi- appear to the root of the mechanism of inheritance of complex ological questionnaire and included biological traits, medical diseases, at least for peripheral genes.7 With this challenging diagnoses, drug use, lifestyle habits and sociodemographic vision, a multitrait genome association analysis of the whole and socioeconomic variables.18 Description of GCAT variables phenome9 becomes a more appropriate way to detect peripheral dataset is available at GCAT (see the URLs section). To keep gene variation effects and new network disturbances affecting as many as possible of the genotyped samples in the study, we core genes. Multitrait analysis approaches are developed for imputed anthropometric missing values (<1%) from the overall research of genetically complex conditions using raw or summa- distribution values using statistical approaches. Missing values ry-level data statistics from GWAS in order to explain the largest (<1%) for biological and anthropometric measures (height, possible amount of the covariation between SNPs and traits.10–15 weight, waist and hip circumference, systolic and diastolic blood The contribution of total genetic variation, known as heri- pressure and heart rate) were imputed by stratifying the whole tability (broad-sense heritability, h2), is estimated now from GCAT cohort by gender and age and using multiple imputa- genome-wide studies in large cohorts directly from SNP data tion by the fully conditional specification method, implemented (known as h2SNP). However, even if most disease conditions in the R mice package.19 For GWAS analysis, we retained all have a strong genetic basis, it is well known that our capacity to variables with at least five observations (n=205). For herita- find genetic effects depends on the overall genetic contribution of bility estimates, only variables with at least 500 individuals per the trait. Overall estimations differed depending on the ancestry, class were retained (n=96) for robustness. The description of sample ascertainment, gender and age of the population under the traits and measures included in this study is summarised in study. Recently, data from the UK Biobank determined genetic online supplementary table S1. contributions with a phenome-based approach16 and identi- fied a shared familial environment as a significant important Genotyping, relatedness and population structure factor besides genetic heritability values in 12 common diseases 17 Genotyping of the 5459 GCAT participants (GCATcore) was analysed. done using the Infinium Expanded Multi-Ethnic Genotyping In this study, we present new data on phenotype-wide estima- Array (MEGAEx) (ILLUMINA, San Diego, California, USA). A tion of the heritability of 205 complex traits (including diseases) customised cluster file was produced from the entire sample and new insights into the genetics of anthropometric traits dataset and used for joint calling. We applied PCA to detect any in a Mediterranean Caucasian population using a two-stage hidden substructure and the method of moments for the estima- meta-analysis approach with multiple-related phenotypes tion of identity by descent probabilities to exclude cases with (MRPs). cryptic relatedness. The extensive QC protocol used for
Recommended publications
  • Entrez Symbols Name Termid Termdesc 117553 Uba3,Ube1c
    Entrez Symbols Name TermID TermDesc 117553 Uba3,Ube1c ubiquitin-like modifier activating enzyme 3 GO:0016881 acid-amino acid ligase activity 299002 G2e3,RGD1310263 G2/M-phase specific E3 ubiquitin ligase GO:0016881 acid-amino acid ligase activity 303614 RGD1310067,Smurf2 SMAD specific E3 ubiquitin protein ligase 2 GO:0016881 acid-amino acid ligase activity 308669 Herc2 hect domain and RLD 2 GO:0016881 acid-amino acid ligase activity 309331 Uhrf2 ubiquitin-like with PHD and ring finger domains 2 GO:0016881 acid-amino acid ligase activity 316395 Hecw2 HECT, C2 and WW domain containing E3 ubiquitin protein ligase 2 GO:0016881 acid-amino acid ligase activity 361866 Hace1 HECT domain and ankyrin repeat containing, E3 ubiquitin protein ligase 1 GO:0016881 acid-amino acid ligase activity 117029 Ccr5,Ckr5,Cmkbr5 chemokine (C-C motif) receptor 5 GO:0003779 actin binding 117538 Waspip,Wip,Wipf1 WAS/WASL interacting protein family, member 1 GO:0003779 actin binding 117557 TM30nm,Tpm3,Tpm5 tropomyosin 3, gamma GO:0003779 actin binding 24779 MGC93554,Slc4a1 solute carrier family 4 (anion exchanger), member 1 GO:0003779 actin binding 24851 Alpha-tm,Tma2,Tmsa,Tpm1 tropomyosin 1, alpha GO:0003779 actin binding 25132 Myo5b,Myr6 myosin Vb GO:0003779 actin binding 25152 Map1a,Mtap1a microtubule-associated protein 1A GO:0003779 actin binding 25230 Add3 adducin 3 (gamma) GO:0003779 actin binding 25386 AQP-2,Aqp2,MGC156502,aquaporin-2aquaporin 2 (collecting duct) GO:0003779 actin binding 25484 MYR5,Myo1e,Myr3 myosin IE GO:0003779 actin binding 25576 14-3-3e1,MGC93547,Ywhah
    [Show full text]
  • Gene Regulation and Speciation in House Mice
    Downloaded from genome.cshlp.org on September 26, 2021 - Published by Cold Spring Harbor Laboratory Press Research Gene regulation and speciation in house mice Katya L. Mack,1 Polly Campbell,2 and Michael W. Nachman1 1Museum of Vertebrate Zoology and Department of Integrative Biology, University of California, Berkeley, California 94720-3160, USA; 2Department of Integrative Biology, Oklahoma State University, Stillwater, Oklahoma 74078, USA One approach to understanding the process of speciation is to characterize the genetic architecture of post-zygotic isolation. As gene regulation requires interactions between loci, negative epistatic interactions between divergent regulatory elements might underlie hybrid incompatibilities and contribute to reproductive isolation. Here, we take advantage of a cross between house mouse subspecies, where hybrid dysfunction is largely unidirectional, to test several key predictions about regulatory divergence and reproductive isolation. Regulatory divergence between Mus musculus musculus and M. m. domesticus was charac- terized by studying allele-specific expression in fertile hybrid males using mRNA-sequencing of whole testes. We found ex- tensive regulatory divergence between M. m. musculus and M. m. domesticus, largely attributable to cis-regulatory changes. When both cis and trans changes occurred, they were observed in opposition much more often than expected under a neutral model, providing strong evidence of widespread compensatory evolution. We also found evidence for lineage-specific positive se- lection on a subset of genes related to transcriptional regulation. Comparisons of fertile and sterile hybrid males identified a set of genes that were uniquely misexpressed in sterile individuals. Lastly, we discovered a nonrandom association between these genes and genes showing evidence of compensatory evolution, consistent with the idea that regulatory interactions might contribute to Dobzhansky-Muller incompatibilities and be important in speciation.
    [Show full text]
  • Molecular Mechanisms Underlying Noncoding Risk Variations in Psychiatric Genetic Studies
    OPEN Molecular Psychiatry (2017) 22, 497–511 www.nature.com/mp REVIEW Molecular mechanisms underlying noncoding risk variations in psychiatric genetic studies X Xiao1,2, H Chang1,2 and M Li1 Recent large-scale genetic approaches such as genome-wide association studies have allowed the identification of common genetic variations that contribute to risk architectures of psychiatric disorders. However, most of these susceptibility variants are located in noncoding genomic regions that usually span multiple genes. As a result, pinpointing the precise variant(s) and biological mechanisms accounting for the risk remains challenging. By reviewing recent progresses in genetics, functional genomics and neurobiology of psychiatric disorders, as well as gene expression analyses of brain tissues, here we propose a roadmap to characterize the roles of noncoding risk loci in the pathogenesis of psychiatric illnesses (that is, identifying the underlying molecular mechanisms explaining the genetic risk conferred by those genomic loci, and recognizing putative functional causative variants). This roadmap involves integration of transcriptomic data, epidemiological and bioinformatic methods, as well as in vitro and in vivo experimental approaches. These tools will promote the translation of genetic discoveries to physiological mechanisms, and ultimately guide the development of preventive, therapeutic and prognostic measures for psychiatric disorders. Molecular Psychiatry (2017) 22, 497–511; doi:10.1038/mp.2016.241; published online 3 January 2017 RECENT GENETIC ANALYSES OF NEUROPSYCHIATRIC neurodevelopment and brain function. For example, GRM3, DISORDERS GRIN2A, SRR and GRIA1 were known to involve in the neuro- Schizophrenia, bipolar disorder, major depressive disorder and transmission mediated by glutamate signaling and synaptic autism are highly prevalent complex neuropsychiatric diseases plasticity.
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • Proteome Analysis of Isolated Podocytes Reveals Stress Responses in Glomerular Sclerosis
    BASIC RESEARCH www.jasn.org Proteome Analysis of Isolated Podocytes Reveals Stress Responses in Glomerular Sclerosis Sybille Koehler,1,2 Alexander Kuczkowski,1 Lucas Kuehne,1 Christian Jüngst ,3 Martin Hoehne ,1,3 Florian Grahammer,4 Sean Eddy ,5 Matthias Kretzler ,5,6 Bodo B. Beck,7 Jörg Höhfeld,8 Bernhard Schermer,1,3 Thomas Benzing,1,3 Paul T. Brinkkoetter,1 and Markus M. Rinschen1,3,9 Due to the number of contributing authors, the affiliations are listed at the end of this article. ABSTRACT Background Understanding podocyte-specific responses to injury at a systems level is difficult because injury leads to podocyte loss or an increase of extracellular matrix, altering glomerular cellular composi- tion. Finding a window into early podocyte injury might help identify molecular pathways involved in the podocyte stress response. Methods We developed an approach to apply proteome analysis to very small samples of purified podo- cyte fractions. To examine podocytes in early disease states in FSGS mouse models, we used podocyte fractions isolated from individual mice after chemical induction of glomerular disease (with Doxorubicin or LPS). We also applied single-glomerular proteome analysis to tissue from patients with FSGS. Results Transcriptome and proteome analysis of glomeruli from patients with FSGS revealed an under- representation of podocyte-specific genes and proteins in late-stage disease. Proteome analysis of puri- fied podocyte fractions from FSGS mouse models showed an early stress response that includes perturbations of metabolic, mechanical, and proteostasis proteins. Additional analysis revealed a high correlation between the amount of proteinuria and expression levels of the mechanosensor protein Filamin-B.
    [Show full text]
  • Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights Into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle
    animals Article Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle Masoumeh Naserkheil 1 , Abolfazl Bahrami 1 , Deukhwan Lee 2,* and Hossein Mehrban 3 1 Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj 77871-31587, Iran; [email protected] (M.N.); [email protected] (A.B.) 2 Department of Animal Life and Environment Sciences, Hankyong National University, Jungang-ro 327, Anseong-si, Gyeonggi-do 17579, Korea 3 Department of Animal Science, Shahrekord University, Shahrekord 88186-34141, Iran; [email protected] * Correspondence: [email protected]; Tel.: +82-31-670-5091 Received: 25 August 2020; Accepted: 6 October 2020; Published: 9 October 2020 Simple Summary: Hanwoo is an indigenous cattle breed in Korea and popular for meat production owing to its rapid growth and high-quality meat. Its yearling weight and carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score) are economically important for the selection of young and proven bulls. In recent decades, the advent of high throughput genotyping technologies has made it possible to perform genome-wide association studies (GWAS) for the detection of genomic regions associated with traits of economic interest in different species. In this study, we conducted a weighted single-step genome-wide association study which combines all genotypes, phenotypes and pedigree data in one step (ssGBLUP). It allows for the use of all SNPs simultaneously along with all phenotypes from genotyped and ungenotyped animals. Our results revealed 33 relevant genomic regions related to the traits of interest.
    [Show full text]
  • Reconstructing Cell Cycle Pseudo Time-Series Via Single-Cell Transcriptome Data—Supplement
    School of Natural Sciences and Mathematics Reconstructing Cell Cycle Pseudo Time-Series Via Single-Cell Transcriptome Data—Supplement UT Dallas Author(s): Michael Q. Zhang Rights: CC BY 4.0 (Attribution) ©2017 The Authors Citation: Liu, Zehua, Huazhe Lou, Kaikun Xie, Hao Wang, et al. 2017. "Reconstructing cell cycle pseudo time-series via single-cell transcriptome data." Nature Communications 8, doi:10.1038/s41467-017-00039-z This document is being made freely available by the Eugene McDermott Library of the University of Texas at Dallas with permission of the copyright owner. All rights are reserved under United States copyright law unless specified otherwise. File name: Supplementary Information Description: Supplementary figures, supplementary tables, supplementary notes, supplementary methods and supplementary references. CCNE1 CCNE1 CCNE1 CCNE1 36 40 32 34 32 35 30 32 28 30 30 28 28 26 24 25 Normalized Expression Normalized Expression Normalized Expression Normalized Expression 26 G1 S G2/M G1 S G2/M G1 S G2/M G1 S G2/M Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage CCNE1 CCNE1 CCNE1 CCNE1 40 32 40 40 35 30 38 30 30 28 36 25 26 20 20 34 Normalized Expression Normalized Expression Normalized Expression 24 Normalized Expression G1 S G2/M G1 S G2/M G1 S G2/M G1 S G2/M Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage Cell Cycle Stage Supplementary Figure 1 | High stochasticity of single-cell gene expression means, as demonstrated by relative expression levels of gene Ccne1 using the mESC-SMARTer data. For every panel, 20 sample cells were randomly selected for each of the three stages, followed by plotting the mean expression levels at each stage.
    [Show full text]
  • WO 2019/079361 Al 25 April 2019 (25.04.2019) W 1P O PCT
    (12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) (19) World Intellectual Property Organization I International Bureau (10) International Publication Number (43) International Publication Date WO 2019/079361 Al 25 April 2019 (25.04.2019) W 1P O PCT (51) International Patent Classification: CA, CH, CL, CN, CO, CR, CU, CZ, DE, DJ, DK, DM, DO, C12Q 1/68 (2018.01) A61P 31/18 (2006.01) DZ, EC, EE, EG, ES, FI, GB, GD, GE, GH, GM, GT, HN, C12Q 1/70 (2006.01) HR, HU, ID, IL, IN, IR, IS, JO, JP, KE, KG, KH, KN, KP, KR, KW, KZ, LA, LC, LK, LR, LS, LU, LY, MA, MD, ME, (21) International Application Number: MG, MK, MN, MW, MX, MY, MZ, NA, NG, NI, NO, NZ, PCT/US2018/056167 OM, PA, PE, PG, PH, PL, PT, QA, RO, RS, RU, RW, SA, (22) International Filing Date: SC, SD, SE, SG, SK, SL, SM, ST, SV, SY, TH, TJ, TM, TN, 16 October 2018 (16. 10.2018) TR, TT, TZ, UA, UG, US, UZ, VC, VN, ZA, ZM, ZW. (25) Filing Language: English (84) Designated States (unless otherwise indicated, for every kind of regional protection available): ARIPO (BW, GH, (26) Publication Language: English GM, KE, LR, LS, MW, MZ, NA, RW, SD, SL, ST, SZ, TZ, (30) Priority Data: UG, ZM, ZW), Eurasian (AM, AZ, BY, KG, KZ, RU, TJ, 62/573,025 16 October 2017 (16. 10.2017) US TM), European (AL, AT, BE, BG, CH, CY, CZ, DE, DK, EE, ES, FI, FR, GB, GR, HR, HU, ΓΕ , IS, IT, LT, LU, LV, (71) Applicant: MASSACHUSETTS INSTITUTE OF MC, MK, MT, NL, NO, PL, PT, RO, RS, SE, SI, SK, SM, TECHNOLOGY [US/US]; 77 Massachusetts Avenue, TR), OAPI (BF, BJ, CF, CG, CI, CM, GA, GN, GQ, GW, Cambridge, Massachusetts 02139 (US).
    [Show full text]
  • Protein Interactions in the Cancer Proteome† Cite This: Mol
    Molecular BioSystems View Article Online PAPER View Journal | View Issue Small-molecule binding sites to explore protein– protein interactions in the cancer proteome† Cite this: Mol. BioSyst., 2016, 12,3067 David Xu,ab Shadia I. Jalal,c George W. Sledge Jr.d and Samy O. Meroueh*aef The Cancer Genome Atlas (TCGA) offers an unprecedented opportunity to identify small-molecule binding sites on proteins with overexpressed mRNA levels that correlate with poor survival. Here, we analyze RNA-seq and clinical data for 10 tumor types to identify genes that are both overexpressed and correlate with patient survival. Protein products of these genes were scanned for binding sites that possess shape and physicochemical properties that can accommodate small-molecule probes or therapeutic agents (druggable). These binding sites were classified as enzyme active sites (ENZ), protein–protein interaction sites (PPI), or other sites whose function is unknown (OTH). Interestingly, the overwhelming majority of binding sites were classified as OTH. We find that ENZ, PPI, and OTH binding sites often occurred on the same structure suggesting that many of these OTH cavities can be used for allosteric modulation of Creative Commons Attribution 3.0 Unported Licence. enzyme activity or protein–protein interactions with small molecules. We discovered several ENZ (PYCR1, QPRT,andHSPA6)andPPI(CASC5, ZBTB32,andCSAD) binding sites on proteins that have been seldom explored in cancer. We also found proteins that have been extensively studied in cancer that have not been previously explored with small molecules that harbor ENZ (PKMYT1, STEAP3,andNNMT) and PPI (HNF4A, MEF2B,andCBX2) binding sites. All binding sites were classified by the signaling pathways to Received 29th March 2016, which the protein that harbors them belongs using KEGG.
    [Show full text]
  • BRAF-MAD1L1 and BRAF-ZC3H7A: Novel Gene Fusion in Rhabdomyosarcoma and Lung Adenocarcinoma
    Case report: BRAF-MAD1L1 and BRAF-ZC3H7A: novel gene fusion in Rhabdomyosarcoma and Lung adenocarcinoma Da Jiang Fourth Hospital of Heibei Medical University Hui Jin Fourth Hospital of Heibei Medical University Xinliang Zhou Fourth Hospital of Heibei Medical University Shaoshuang Fan Fourth Hospital of Heibei Medical University Mengping Lei Origimed Mian Xu Origimed Xiaoyan Chen ( [email protected] ) The Hospital of Shunyi District Beijing Case Report Keywords: novel gene fusion, BRAF-MAD1L1, BRAF-ZC3H7A, Rhabdomyosarcoma, Lung adenocarcinoma Posted Date: August 26th, 2020 DOI: https://doi.org/10.21203/rs.3.rs-58002/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Page 1/13 Abstract Background Rhabdomyosarcoma (RMS) and lung adenocarcinoma (LADC) epitomizes the success of cancer prevention by the development of conventional therapy, but huge challenges remain in the therapy of advanced diseases. Case presentation We reported two cases of novel BRAF gene fusion. The rst case was a 34-year-old female with RMS harboring a BRAF-MAD1L1 fusion. She suffered tumor resection, recurrence and rapid progression. The second case was a 72-year-old female with LADC harboring a BRAF-ZC3H7A fusion, and she gained rapid progression after receiving a rst-line course of chemotherapy. Conclusions These two BRAF fusions retain the intact BRAF kinase domain (exon 11-18) and showed poor prognosis in RMS and LADC. Background Rhabdomyosarcoma (RMS) is an aggressive malignancy of soft-tissue in children and adolescents (1). There are 4.5 cases per 1 million in children per year and it accounts for 4–5% of all childhood malignancies (2).
    [Show full text]
  • Human Lectins, Their Carbohydrate Affinities and Where to Find Them
    biomolecules Review Human Lectins, Their Carbohydrate Affinities and Where to Review HumanFind Them Lectins, Their Carbohydrate Affinities and Where to FindCláudia ThemD. Raposo 1,*, André B. Canelas 2 and M. Teresa Barros 1 1, 2 1 Cláudia D. Raposo * , Andr1 é LAQVB. Canelas‐Requimte,and Department M. Teresa of Chemistry, Barros NOVA School of Science and Technology, Universidade NOVA de Lisboa, 2829‐516 Caparica, Portugal; [email protected] 12 GlanbiaLAQV-Requimte,‐AgriChemWhey, Department Lisheen of Chemistry, Mine, Killoran, NOVA Moyne, School E41 of ScienceR622 Co. and Tipperary, Technology, Ireland; canelas‐ [email protected] NOVA de Lisboa, 2829-516 Caparica, Portugal; [email protected] 2* Correspondence:Glanbia-AgriChemWhey, [email protected]; Lisheen Mine, Tel.: Killoran, +351‐212948550 Moyne, E41 R622 Tipperary, Ireland; [email protected] * Correspondence: [email protected]; Tel.: +351-212948550 Abstract: Lectins are a class of proteins responsible for several biological roles such as cell‐cell in‐ Abstract:teractions,Lectins signaling are pathways, a class of and proteins several responsible innate immune for several responses biological against roles pathogens. such as Since cell-cell lec‐ interactions,tins are able signalingto bind to pathways, carbohydrates, and several they can innate be a immuneviable target responses for targeted against drug pathogens. delivery Since sys‐ lectinstems. In are fact, able several to bind lectins to carbohydrates, were approved they by canFood be and a viable Drug targetAdministration for targeted for drugthat purpose. delivery systems.Information In fact, about several specific lectins carbohydrate were approved recognition by Food by andlectin Drug receptors Administration was gathered for that herein, purpose. plus Informationthe specific organs about specific where those carbohydrate lectins can recognition be found by within lectin the receptors human was body.
    [Show full text]
  • Identification of Potential Pathogenic Genes Associated with Osteoporosis
    610.BJBJR0010.1302/2046-3758.612.BJR-2017-0102 research-article2017 Freely available online OPEN ACCESS BJR RESEARCH Identification of potential pathogenic genes associated with osteoporosis Objectives B. Xia, Osteoporosis is a chronic disease. The aim of this study was to identify key genes in osteo- Y. Li, porosis. J. Zhou, Methods B. Tian, Microarray data sets GSE56815 and GSE56814, comprising 67 osteoporosis blood samples L. Feng and 62 control blood samples, were obtained from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified in osteoporosis using Limma pack- Jining No. 1 People’s age (3.2.1) and Meta-MA packages. Gene Ontology and Kyoto Encyclopedia of Genes and Hospital, Jining, Genomes enrichment analyses were performed to identify biological functions. Further- Shandong Province, more, the transcriptional regulatory network was established between the top 20 DEGs and China transcriptional factors using the UCSC ENCODE Genome Browser. Receiver operating char- acteristic (ROC) analysis was applied to investigate the diagnostic value of several DEGs. Results A total of 1320 DEGs were obtained, of which 855 were up-regulated and 465 were down- regulated. These differentially expressed genes were enriched in Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways, mainly associated with gene expres- sion and osteoclast differentiation. In the transcriptional regulatory network, there were 6038 interactions pairs involving 88 transcriptional factors. In addition, the quantitative reverse transcriptase-polymerase chain reaction result validated the expression of several genes (VPS35, FCGR2A, TBCA, HIRA, TYROBP, and JUND). Finally, ROC analyses showed that VPS35, HIRA, PHF20 and NFKB2 had a significant diagnostic value for osteoporosis.
    [Show full text]