Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis

Osteoarthritis is a common complex disease imposing a large public-health burden. Here, we performed a genome-wide association study for osteoarthritis, using data across 16.5 million variants from the UK Biobank resource. After performing replication and meta-analysis in up to 30,727 cases and 297,191 controls, we identified nine new osteoarthritis loci, in all of which the most likely causal variant was noncoding. For three loci, we detected association with biologically relevant radiographic endophenotypes, and in five signals we identified genes that were differentially expressed in degraded compared with intact articular cartilage from patients with osteoarthritis. We established causal effects on osteoarthritis for higher body mass index but not for triglyceride levels or genetic predisposition to type 2 diabetes.Genome-wide association study for osteoarthritis using data from UK Biobank identifies loci for knee- and hip-specific disease. Functional analyses of chondrocytes provide further insight into candidate causal genes.

[1]  E. Zeggini,et al.  Integrative epigenomics, transcriptomics and proteomics of patient chondrocytes reveal genes and pathways involved in osteoarthritis , 2016, bioRxiv.

[2]  E. Zeggini,et al.  Pathways to understanding the genomic aetiology of osteoarthritis , 2017, Human molecular genetics.

[3]  Qian Chen,et al.  CD95-induced osteoarthritic chondrocyte apoptosis and necrosis: dependency on p38 mitogen-activated protein kinase , 2006, Arthritis research & therapy.

[4]  J. Buxbaum,et al.  A SPECTRAL APPROACH INTEGRATING FUNCTIONAL GENOMIC ANNOTATIONS FOR CODING AND NONCODING VARIANTS , 2015, Nature Genetics.

[5]  Francisco S. Roque,et al.  A large-scale analysis of tissue-specific pathology and gene expression of human disease genes and complexes , 2008, Proceedings of the National Academy of Sciences.

[6]  Gunwoo Kim,et al.  Transforming growth factor alpha controls the transition from hypertrophic cartilage to bone during endochondral bone growth. , 2012, Bone.

[7]  Jie Huang,et al.  Discovery and refinement of genetic loci associated with cardiometabolic risk using dense imputation maps , 2016, Nature Genetics.

[8]  Yin Xiao,et al.  ERK-1/2 and p38 in the regulation of hypertrophic changes of normal articular cartilage chondrocytes induced by osteoarthritic subchondral osteoblasts. , 2010, Arthritis and rheumatism.

[9]  T. Aigner,et al.  Transforming growth factor alpha suppression of articular chondrocyte phenotype and Sox9 expression in a rat model of osteoarthritis. , 2007, Arthritis and rheumatism.

[10]  Ying Zhang,et al.  A review of crosstalk between MAPK and Wnt signals and its impact on cartilage regeneration , 2014, Cell and Tissue Research.

[11]  S. Ebrahim,et al.  'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? , 2003, International journal of epidemiology.

[12]  Tom R. Gaunt,et al.  LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis , 2016 .

[13]  A. Palmieri,et al.  Comparison between osteoblasts derived from human dental pulp stem cells and osteosarcoma cell lines , 2008, Cell biology international.

[14]  J. Shendure,et al.  A general framework for estimating the relative pathogenicity of human genetic variants , 2014, Nature Genetics.

[15]  J. Mort,et al.  Rho/ROCK and MEK/ERK activation by transforming growth factor-α induces articular cartilage degradation , 2010, Laboratory Investigation.

[16]  C. Hiller,et al.  Self-reported knee pain and disability among healthy individuals: reference data and factors associated with the Knee injury and Osteoarthritis Outcome Score (KOOS) and KOOS-Child. , 2017, Osteoarthritis and Cartilage.

[17]  Jonathan A C Sterne,et al.  The impact of residual and unmeasured confounding in epidemiologic studies: a simulation study. , 2007, American journal of epidemiology.

[18]  Rob Patro,et al.  Salmon provides fast and bias-aware quantification of transcript expression , 2017, Nature Methods.

[19]  David T. Broome,et al.  Mitogen-activated protein kinase phosphatase-1: function and regulation in bone and related tissues , 2016, Connective tissue research.

[20]  Yin Xiao,et al.  Osteoarthritic cartilage chondrocytes alter subchondral bone osteoblast differentiation via MAPK signalling pathway involving ERK1/2. , 2010, Bone.

[21]  Trairak Pisitkun,et al.  Large-scale proteomics and phosphoproteomics of urinary exosomes. , 2009, Journal of the American Society of Nephrology : JASN.

[22]  Gabor T. Marth,et al.  A global reference for human genetic variation , 2015, Nature.

[23]  P. Casey,et al.  Protein prenylation: unique fats make their mark on biology , 2016, Nature Reviews Molecular Cell Biology.

[24]  Tom R. Gaunt,et al.  The UK10K project identifies rare variants in health and disease , 2016 .

[25]  M. Daly,et al.  Proteins Encoded in Genomic Regions Associated with Immune-Mediated Disease Physically Interact and Suggest Underlying Biology , 2011, PLoS genetics.

[26]  J. Jordan,et al.  Associations of educational attainment, occupation and community poverty with knee osteoarthritis in the Johnston County (North Carolina) osteoarthritis project , 2011, Arthritis research & therapy.

[27]  M. Daly,et al.  LD Score regression distinguishes confounding from polygenicity in genome-wide association studies , 2014, Nature Genetics.

[28]  M. Viergever,et al.  CHECK (Cohort Hip and Cohort Knee): similarities and differences with the Osteoarthritis Initiative , 2008, Annals of the rheumatic diseases.

[29]  R. Derynck,et al.  Epidermal growth factor receptor function is necessary for normal craniofacial development and palate closure , 1999, Nature Genetics.

[30]  Joris M. Mooij,et al.  MAGMA: Generalized Gene-Set Analysis of GWAS Data , 2015, PLoS Comput. Biol..

[31]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[32]  J. Cheverud,et al.  A simple correction for multiple comparisons in interval mapping genome scans , 2001, Heredity.

[33]  Frank P Luyten,et al.  GDF5 deficiency in mice is associated with instability-driven joint damage, gait and subchondral bone changes , 2010, Annals of the rheumatic diseases.

[34]  Shankung Lin,et al.  Zinc finger protein ZFP36L1 promotes osteoblastic differentiation but represses adipogenic differentiation of mouse multipotent cells , 2017, Oncotarget.

[35]  T. Capellini,et al.  Ancient selection for derived alleles at a GDF5 enhancer influencing human growth and osteoarthritis risk , 2017, Nature Genetics.

[36]  Brian V. Jenkins,et al.  A mouse forward genetics screen identifies LISTERIN as an E3 ubiquitin ligase involved in neurodegeneration , 2009, Proceedings of the National Academy of Sciences.

[37]  G. Walko,et al.  Plectin–intermediate filament partnership in skin, skeletal muscle, and peripheral nerve , 2013, Histochemistry and Cell Biology.

[38]  S. Murakami,et al.  Constitutive activation of MKK6 in chondrocytes of transgenic mice inhibits proliferation and delays endochondral bone formation. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[39]  S. Ochiai,et al.  Oxidative stress reaction in the meniscus of Bach 1 deficient mice: Potential prevention of meniscal degeneration , 2008, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[40]  Gregory A. Poland,et al.  Fine Mapping Causal Variants with an Approximate Bayesian Method Using Marginal Test Statistics , 2015, Genetics.

[41]  M. Lotz,et al.  Bach1 deficiency reduces severity of osteoarthritis through upregulation of heme oxygenase-1 , 2015, Arthritis Research & Therapy.

[42]  P. Donnelly,et al.  A new multipoint method for genome-wide association studies by imputation of genotypes , 2007, Nature Genetics.

[43]  Susumu Goto,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 2000, Nucleic Acids Res..

[44]  Charity W. Law,et al.  voom: precision weights unlock linear model analysis tools for RNA-seq read counts , 2014, Genome Biology.

[45]  D. Moore,et al.  Reduced limb length and worsened osteoarthritis in adult mice after genetic inhibition of p38 MAP kinase activity in cartilage. , 2008, Arthritis and rheumatism.

[46]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[47]  Beatriz Gámez,et al.  p38 MAPK Signaling in Osteoblast Differentiation , 2016, Front. Cell Dev. Biol..

[48]  H. Stefánsson,et al.  Whole-genome sequencing identifies rare genotypes in COMP and CHADL associated with high risk of hip osteoarthritis , 2017, Nature Genetics.

[49]  C. Cooper,et al.  An increased rate of falling leads to a rise in fracture risk in postmenopausal women with self-reported osteoarthritis: a prospective multinational cohort study (GLOW) , 2012, Annals of the rheumatic diseases.

[50]  V. Moskvina,et al.  On multiple‐testing correction in genome‐wide association studies , 2008, Genetic epidemiology.

[51]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[52]  Hashem A. Shihab,et al.  MR-Base: a platform for systematic causal inference across the phenome using billions of genetic associations , 2016, bioRxiv.

[53]  S. Goldring,et al.  Osteoarthritis: a disease of the joint as an organ. , 2012, Arthritis and rheumatism.

[54]  Hollis G. Potter,et al.  Author Manuscript , 2013 .

[55]  Timothy E. McAlindon,et al.  Osteoarthritis: New Insights , 2002 .

[56]  Yusuke Nakamura,et al.  A functional polymorphism in the 5′ UTR of GDF5 is associated with susceptibility to osteoarthritis , 2007, Nature Genetics.

[57]  T. Spector,et al.  Osteoarthritis: New Insights. Part 1: The Disease and Its Risk Factors , 2000, Annals of Internal Medicine.

[58]  N. Eriksson,et al.  Genome-wide analysis identifies 12 loci influencing human reproductive behavior , 2016 .

[59]  Yun Li,et al.  METAL: fast and efficient meta-analysis of genomewide association scans , 2010, Bioinform..

[60]  I. Welch,et al.  Context-specific protection of TGFα null mice from osteoarthritis , 2016, Scientific Reports.

[61]  A. Hofman,et al.  Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis , 2016, PLoS genetics.

[62]  Yu Sun,et al.  Bioinformatics analysis of gene expression profiles of osteoarthritis. , 2015, Acta histochemica.

[63]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[64]  Michael Boehnke,et al.  Recommended Joint and Meta‐Analysis Strategies for Case‐Control Association Testing of Single Low‐Count Variants , 2013, Genetic epidemiology.

[65]  German Tischler,et al.  biobambam: tools for read pair collation based algorithms on BAM files , 2013, Source Code for Biology and Medicine.

[66]  Yin Xiao,et al.  Inhibition of p38 pathway leads to OA-like changes in a rat animal model. , 2012, Rheumatology.

[67]  John Loughlin,et al.  An SNP in the 5'-UTR of GDF5 is associated with osteoarthritis susceptibility in Europeans and with in vivo differences in allelic expression in articular cartilage. , 2007, Human molecular genetics.

[68]  R. Pearson,et al.  Histopathology grading systems for characterisation of human knee osteoarthritis--reproducibility, variability, reliability, correlation, and validity. , 2011, Osteoarthritis and cartilage.

[69]  M. Englund,et al.  Association between occupation and knee and hip replacement due to osteoarthritis: a case-control study , 2010, Arthritis research & therapy.

[70]  F. Cicuttini,et al.  How Are Obesity and Body Composition Related to Patellar Cartilage? A Systematic Review , 2017, The Journal of Rheumatology.