Application of multiple statistical tests to enhance mass spectrometry-based biomarker discovery

<p>Abstract</p> <p>Background</p> <p>Mass spectrometry-based biomarker discovery has long been hampered by the difficulty in reconciling lists of discriminatory peaks identified by different laboratories for the same diseases studied. We describe a multi-statistical analysis procedure that combines several independent computational methods. This approach capitalizes on the strengths of each to analyze the same high-resolution mass spectral data set to discover consensus differential mass peaks that should be robust biomarkers for distinguishing between disease states.</p> <p>Results</p> <p>The proposed methodology was applied to a pilot narcolepsy study using logistic regression, hierarchical clustering, t-test, and CART. Consensus, differential mass peaks with high predictive power were identified across three of the four statistical platforms. Based on the diagnostic accuracy measures investigated, the performance of the consensus-peak model was a compromise between logistic regression and CART, which produced better models than hierarchical clustering and t-test. However, consensus peaks confer a higher level of confidence in their ability to distinguish between disease states since they do not represent peaks that are a result of biases to a particular statistical algorithm. Instead, they were selected as differential across differing data distribution assumptions, demonstrating their true discriminatory potential.</p> <p>Conclusion</p> <p>The methodology described here is applicable to any high-resolution MALDI mass spectrometry-derived data set with minimal mass drift which is essential for peak-to-peak comparison studies. Four statistical approaches with differing data distribution assumptions were applied to the same raw data set to obtain consensus peaks that were found to be statistically differential between the two groups compared. These consensus peaks demonstrated high diagnostic accuracy when used to form a predictive model as evaluated by receiver operating characteristics curve analysis. They should demonstrate a higher discriminatory ability as they are not biased to a particular algorithm. Thus, they are prime candidates for downstream identification and validation efforts.</p

Genomic correlates of glatiramer acetate adverse cardiovascular effects lead to a novel locus mediating coronary risk

Glatiramer acetate is used therapeutically in multiple sclerosis but also known for adverse effects including elevated coronary artery disease (CAD) risk. The mechanisms underlying the cardiovascular side effects of the medication are unclear. Here, we made use of the chromosomal variation in the genes that are known to be affected by glatiramer treatment. Focusing on genes and gene products reported by drug-gene interaction database to interact with glatiramer acetate we explored a large meta-analysis on CAD genome-wide association studies aiming firstly, to investigate whether variants in these genes also affect cardiovascular risk and secondly, to identify new CAD risk genes. We traced association signals in a 200-kb region around genomic positions of genes interacting with glatiramer in up to 60 801 CAD cases and 123 504 controls. We validated the identified association in additional 21 934 CAD cases and 76 087 controls. We identified three new CAD risk alleles within the TGFB1 region on chromosome 19 that independently affect CAD risk. The lead SNP rs12459996 was genome-wide significantly associated with CAD in the extended meta-analysis (odds ratio 1.09, p = 1.58×10-12). The other two SNPs at the locus were not in linkage disequilibrium with the lead SNP and by a conditional analysis showed p-values of 4.05 × 10-10 and 2.21 × 10-6. Thus, studying genes reported to interact with glatiramer acetate we identified genetic variants that concordantly with the drug increase the risk of CAD. Of these, TGFB1 displayed signal for association. Indeed, the gene has been associated with CAD previously in both in vivo and in vitro studies. Here we establish genome-wide significant association with CAD in large human samples.This work was supported by grants from the Fondation Leducq (CADgenomics: Understanding CAD Genes, 12CVD02), the German Federal Ministry of Education and Research (BMBF) within the framework of the e:Med research and funding concept (e:AtheroSysMed, grant 01ZX1313A-2014 and SysInflame, grant 01ZX1306A), and the European Union Seventh Framework Programme FP7/2007-2013 under grant agreement no HEALTH-F2-2013-601456 (CVgenes-at-target). Further grants were received from the DFG as part of the Sonderforschungsbereich CRC 1123 (B2). T.K. was supported by a DZHK Rotation Grant. I.B. was supported by the Deutsche Forschungsgemeinschaft (DFG) cluster of excellence ‘Inflammation at Interfaces’. F.W.A. is supported by a Dekker scholarship-Junior Staff Member 2014T001 - Netherlands Heart Foundation and UCL Hospitals NIHR Biomedical Research Centre

Fine-mapping of prostate cancer susceptibility loci in a large meta-analysis identifies candidate causal variants

Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling. © 2018 The Author(s).Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling. © 2018 The Author(s).Peer reviewe

Germline variation at 8q24 and prostate cancer risk in men of European ancestry

Chromosome 8q24 is a susceptibility locus for multiple cancers, including prostate cancer. Here we combine genetic data across the 8q24 susceptibility region from 71,535 prostate cancer cases and 52,935 controls of European ancestry to define the overall contribution of germline variation at 8q24 to prostate cancer risk. We identify 12 independent risk signals for prostate cancer (p < 4.28 × 10−15), including three risk variants that have yet to be reported. From a polygenic risk score (PRS) model, derived to assess the cumulative effect of risk variants at 8q24, men in the top 1% of the PRS have a 4-fold (95%CI = 3.62–4.40) greater risk compared to the population average. These 12 variants account for ~25% of what can be currently explained of the familial risk of prostate cancer by known genetic risk factors. These findings highlight the overwhelming contribution of germline variation at 8q24 on prostate cancer risk which has implications for population risk stratification

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer.

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead

NIST Interlaboratory Study on Glycosylation Analysis of Monoclonal Antibodies: Comparison of Results from Diverse Analytical Methods

Glycosylation is a topic of intense current interest in the development of biopharmaceuticals because it is related to drug safety and efficacy. This work describes results of an interlaboratory study on the glycosylation of the Primary Sample (PS) of NISTmAb, a monoclonal antibody reference material. Seventy-six laboratories from industry, university, research, government, and hospital sectors in Europe, North America, Asia, and Australia submit- Avenue, Silver Spring, Maryland 20993; 22Glycoscience Research Laboratory, Genos, Borongajska cesta 83h, 10 000 Zagreb, Croatia; 23Faculty of Pharmacy and Biochemistry, University of Zagreb, A. Kovacˇ ic´ a 1, 10 000 Zagreb, Croatia; 24Department of Chemistry, Georgia State University, 100 Piedmont Avenue, Atlanta, Georgia 30303; 25glyXera GmbH, Brenneckestrasse 20 * ZENIT / 39120 Magdeburg, Germany; 26Health Products and Foods Branch, Health Canada, AL 2201E, 251 Sir Frederick Banting Driveway, Ottawa, Ontario, K1A 0K9 Canada; 27Graduate School of Advanced Sciences of Matter, Hiroshima University, 1-3-1 Kagamiyama Higashi-Hiroshima 739–8530 Japan; 28ImmunoGen, 830 Winter Street, Waltham, Massachusetts 02451; 29Department of Medical Physiology, Jagiellonian University Medical College, ul. Michalowskiego 12, 31–126 Krakow, Poland; 30Department of Pathology, Johns Hopkins University, 400 N. Broadway Street Baltimore, Maryland 21287; 31Mass Spec Core Facility, KBI Biopharma, 1101 Hamlin Road Durham, North Carolina 27704; 32Division of Mass Spectrometry, Korea Basic Science Institute, 162 YeonGuDanji-Ro, Ochang-eup, Cheongwon-gu, Cheongju Chungbuk, 363–883 Korea (South); 33Advanced Therapy Products Research Division, Korea National Institute of Food and Drug Safety, 187 Osongsaengmyeong 2-ro Osong-eup, Heungdeok-gu, Cheongju-si, Chungcheongbuk-do, 363–700, Korea (South); 34Center for Proteomics and Metabolomics, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, The Netherlands; 35Ludger Limited, Culham Science Centre, Abingdon, Oxfordshire, OX14 3EB, United Kingdom; 36Biomolecular Discovery and Design Research Centre and ARC Centre of Excellence for Nanoscale BioPhotonics (CNBP), Macquarie University, North Ryde, Australia; 37Proteomics, Central European Institute for Technology, Masaryk University, Kamenice 5, A26, 625 00 BRNO, Czech Republic; 38Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany; 39Department of Biomolecular Sciences, Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany; 40AstraZeneca, Granta Park, Cambridgeshire, CB21 6GH United Kingdom; 41Merck, 2015 Galloping Hill Rd, Kenilworth, New Jersey 07033; 42Analytical R&D, MilliporeSigma, 2909 Laclede Ave. St. Louis, Missouri 63103; 43MS Bioworks, LLC, 3950 Varsity Drive Ann Arbor, Michigan 48108; 44MSD, Molenstraat 110, 5342 CC Oss, The Netherlands; 45Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, 5–1 Higashiyama, Myodaiji, Okazaki 444–8787 Japan; 46Graduate School of Pharmaceutical Sciences, Nagoya City University, 3–1 Tanabe-dori, Mizuhoku, Nagoya 467–8603 Japan; 47Medical & Biological Laboratories Co., Ltd, 2-22-8 Chikusa, Chikusa-ku, Nagoya 464–0858 Japan; 48National Institute for Biological Standards and Control, Blanche Lane, South Mimms, Potters Bar, Hertfordshire EN6 3QG United Kingdom; 49Division of Biological Chemistry & Biologicals, National Institute of Health Sciences, 1-18-1 Kamiyoga, Setagaya-ku, Tokyo 158–8501 Japan; 50New England Biolabs, Inc., 240 County Road, Ipswich, Massachusetts 01938; 51New York University, 100 Washington Square East New York City, New York 10003; 52Target Discovery Institute, Nuffield Department of Medicine, University of Oxford, Roosevelt Drive, Oxford, OX3 7FZ, United Kingdom; 53GlycoScience Group, The National Institute for Bioprocessing Research and Training, Fosters Avenue, Mount Merrion, Blackrock, Co. Dublin, Ireland; 54Department of Chemistry, North Carolina State University, 2620 Yarborough Drive Raleigh, North Carolina 27695; 55Pantheon, 201 College Road East Princeton, New Jersey 08540; 56Pfizer Inc., 1 Burtt Road Andover, Massachusetts 01810; 57Proteodynamics, ZI La Varenne 20–22 rue Henri et Gilberte Goudier 63200 RIOM, France; 58ProZyme, Inc., 3832 Bay Center Place Hayward, California 94545; 59Koichi Tanaka Mass Spectrometry Research Laboratory, Shimadzu Corporation, 1 Nishinokyo Kuwabara-cho Nakagyo-ku, Kyoto, 604 8511 Japan; 60Children’s GMP LLC, St. Jude Children’s Research Hospital, 262 Danny Thomas Place Memphis, Tennessee 38105; 61Sumitomo Bakelite Co., Ltd., 1–5 Muromati 1-Chome, Nishiku, Kobe, 651–2241 Japan; 62Synthon Biopharmaceuticals, Microweg 22 P.O. Box 7071, 6503 GN Nijmegen, The Netherlands; 63Takeda Pharmaceuticals International Co., 40 Landsdowne Street Cambridge, Massachusetts 02139; 64Department of Chemistry and Biochemistry, Texas Tech University, 2500 Broadway, Lubbock, Texas 79409; 65Thermo Fisher Scientific, 1214 Oakmead Parkway Sunnyvale, California 94085; 66United States Pharmacopeia India Pvt. Ltd. IKP Knowledge Park, Genome Valley, Shamirpet, Turkapally Village, Medchal District, Hyderabad 500 101 Telangana, India; 67Alberta Glycomics Centre, University of Alberta, Edmonton, Alberta T6G 2G2 Canada; 68Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2 Canada; 69Department of Chemistry, University of California, One Shields Ave, Davis, California 95616; 70Horva´ th Csaba Memorial Laboratory for Bioseparation Sciences, Research Center for Molecular Medicine, Doctoral School of Molecular Medicine, Faculty of Medicine, University of Debrecen, Debrecen, Egyetem ter 1, Hungary; 71Translational Glycomics Research Group, Research Institute of Biomolecular and Chemical Engineering, University of Pannonia, Veszprem, Egyetem ut 10, Hungary; 72Delaware Biotechnology Institute, University of Delaware, 15 Innovation Way Newark, Delaware 19711; 73Proteomics Core Facility, University of Gothenburg, Medicinaregatan 1G SE 41390 Gothenburg, Sweden; 74Department of Medical Biochemistry and Cell Biology, University of Gothenburg, Institute of Biomedicine, Sahlgrenska Academy, Medicinaregatan 9A, Box 440, 405 30, Gothenburg, Sweden; 75Department of Clinical Chemistry and Transfusion Medicine, Sahlgrenska Academy at the University of Gothenburg, Bruna Straket 16, 41345 Gothenburg, Sweden; 76Department of Chemistry, University of Hamburg, Martin Luther King Pl. 6 20146 Hamburg, Germany; 77Department of Chemistry, University of Manitoba, 144 Dysart Road, Winnipeg, Manitoba, Canada R3T 2N2; 78Laboratory of Mass Spectrometry of Interactions and Systems, University of Strasbourg, UMR Unistra-CNRS 7140, France; 79Natural and Medical Sciences Institute, University of Tu¨ bingen, Markwiesenstrae 55, 72770 Reutlingen, Germany; 80Bijvoet Center for Biomolecular Research and Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands; 81Division of Bioanalytical Chemistry, Amsterdam Institute for Molecules, Medicines and Systems, Vrije Universiteit Amsterdam, de Boelelaan 1085, 1081 HV Amsterdam, The Netherlands; 82Department of Chemistry, Waters Corporation, 34 Maple Street Milford, Massachusetts 01757; 83Zoetis, 333 Portage St. Kalamazoo, Michigan 49007 Author’s Choice—Final version open access under the terms of the Creative Commons CC-BY license. Received July 24, 2019, and in revised form, August 26, 2019 Published, MCP Papers in Press, October 7, 2019, DOI 10.1074/mcp.RA119.001677 ER: NISTmAb Glycosylation Interlaboratory Study 12 Molecular & Cellular Proteomics 19.1 Downloaded from https://www.mcponline.org by guest on January 20, 2020 ted a total of 103 reports on glycan distributions. The principal objective of this study was to report and compare results for the full range of analytical methods presently used in the glycosylation analysis of mAbs. Therefore, participation was unrestricted, with laboratories choosing their own measurement techniques. Protein glycosylation was determined in various ways, including at the level of intact mAb, protein fragments, glycopeptides, or released glycans, using a wide variety of methods for derivatization, separation, identification, and quantification. Consequently, the diversity of results was enormous, with the number of glycan compositions identified by each laboratory ranging from 4 to 48. In total, one hundred sixteen glycan compositions were reported, of which 57 compositions could be assigned consensus abundance values. These consensus medians provide communityderived values for NISTmAb PS. Agreement with the consensus medians did not depend on the specific method or laboratory type. The study provides a view of the current state-of-the-art for biologic glycosylation measurement and suggests a clear need for harmonization of glycosylation analysis methods. Molecular & Cellular Proteomics 19: 11–30, 2020. DOI: 10.1074/mcp.RA119.001677.L

Polygenic Risk Modelling for Prediction of Epithelial Ovarian Cancer Risk

Funder: Funding details are provided in the Supplementary MaterialAbstractPolygenic risk scores (PRS) for epithelial ovarian cancer (EOC) have the potential to improve risk stratification. Joint estimation of Single Nucleotide Polymorphism (SNP) effects in models could improve predictive performance over standard approaches of PRS construction. Here, we implemented computationally-efficient, penalized, logistic regression models (lasso, elastic net, stepwise) to individual level genotype data and a Bayesian framework with continuous shrinkage, “select and shrink for summary statistics” (S4), to summary level data for epithelial non-mucinous ovarian cancer risk prediction. We developed the models in a dataset consisting of 23,564 non-mucinous EOC cases and 40,138 controls participating in the Ovarian Cancer Association Consortium (OCAC) and validated the best models in three populations of different ancestries: prospective data from 198,101 women of European ancestry; 7,669 women of East Asian ancestry; 1,072 women of African ancestry, and in 18,915 BRCA1 and 12,337 BRCA2 pathogenic variant carriers of European ancestry. In the external validation data, the model with the strongest association for non-mucinous EOC risk derived from the OCAC model development data was the S4 model (27,240 SNPs) with odds ratios (OR) of 1.38(95%CI:1.28–1.48,AUC:0.588) per unit standard deviation, in women of European ancestry; 1.14(95%CI:1.08–1.19,AUC:0.538) in women of East Asian ancestry; 1.38(95%CI:1.21-1.58,AUC:0.593) in women of African ancestry; hazard ratios of 1.37(95%CI:1.30–1.44,AUC:0.592) in BRCA1 pathogenic variant carriers and 1.51(95%CI:1.36-1.67,AUC:0.624) in BRCA2 pathogenic variant carriers. Incorporation of the S4 PRS in risk prediction models for ovarian cancer may have clinical utility in ovarian cancer prevention programs.</jats:p