fetching data ...

OP0017 (2021)
VALIDATION OF GWAS-IDENTIFIED VARIANTS FOR ANTI-TNF DRUG RESPONSE IN RHEUMATOID ARTHRITIS: A META-ANALYSIS OF THREE LARGE COHORTS
J. M. Sánchez-Maldonado1, R. Cáliz Cáliz2, M. Á. López-Nevot3, A. Moñiz-Díez1, A. J. Cabrera-Serrano1, H. Canhão4, R. Ter Horst5, A. Escudero Contreras6, S. Sorensen7, M. L. Hetland8, M. Á. Ferrer González2, B. Glintborg8,9, I. Filipescu10, E. Perez-Pampín11, P. Conesa-Zamora12, J. Swierkot13, A. Den Broeder14, S. De Vita15, S. E. Hove Jaconsen16, L. Quartuccio15, Y. LI5, M. Netea5, M. J. H. Coenen14, V. Andersen8,9, J. E. Fonseca17, K. Bogunia-Kubik18, E. Collantes Estevez6, J. Sainz1,19, on behalf of REPAIR consortium
1GENYO. Center for Genomics and Oncological Research, Genomic Oncology, Granada, Spain
2Virgen de las Nieves University Hospital, Rheumatology Deparment, Granada, Spain
3Virgen de las Nieves University Hospital, Immunology Department, Granada, Spain
4Universidade Nova de Lisboa, CEDOC, EpiDoC Unit, NOVA Medical School and National School of Public Health, Lisboa, Portugal
5Radboud University Nijmegen Medical Center, Department of Internal Medicine and Radboud Center for Infectious Diseases, Nijmegen, Netherlands
6Reina Sofía Hospital/IMIBIC/University of Córdoba, Rheumatology Deparment, Córdoba, Spain
7Hospital of Southern Jutland, Focused Research Unit for Molecular Diagnostic and Clinical Research, Aabenraa, Denmark
8Rigshospitalet, The DANBIO Registry, The Danish Rheumatologic Biobank and Copenhagen Center for Arthritis Research (COPECARE), Center for Rheumatology ad Spine Diseases, Centre of Head and Orthopaedics, Glostrup, Denmark
9University of Copenhagen, Department of Clinical, Faculty of Health and Medical Sciences, Copenhagen, Denmark
10University of Medicine and Pharmacy “Iuliu Hatieganu” Cluj-Napoca, Rheumatology Department, Cluj-Napoca, Romania
11University Hospital of Santiago de Compostela, Rheumatology Unit, Santiago de Compostela, Spain
12Santa Lucía University Hospital, Clinical Analysis department, Cartagena, Spain
13Wroclaw Medical University, Department and Clinic of Rheumatology and Internal Medicine, Warsaw, Poland
14Radboud University Medical Center, Radboud Institute for Health Sciences, Department of Human Genetics, Nijmegen, Netherlands
15University of Udine, Department of Medical Area, Udine, Italy
16University Hospital of Southern Jutland, Department of Biochemistry and Immunology, Aabenraa, Denmark
17Hospital de Santa Maria, CHLN, Rheumatology and Metabolic Bone Diseases Department, Lisbon, Portugal
18Polish Academy of Sciences Wrocław, Hirszfeld Institute of Immunology and Experimental Therapy, Warsaw, Poland
19University of Granada, Department of Medicine, Granada, Spain

Background: The interplay between genetics and drug response in rheumatoid arthritis (RA) has shown that response to biologics varies between individuals and that a large proportion of patients show no clinical improvement (Plenge and Bridges, 2011). Despite the disappointing scenario, to date, only a few genetic markers have been consistently identified and we are far from being able to optimize drug dosing or prioritize drug combinations based on genetic findings.


Objectives: With this background, we sought to validate the association of GWAS-identified variants for response to TNF inhibitors (TNFi) in a two-stage case control association study and to shed some light into the functional role of the most interesting markers.


Methods: The discovery population consisted of 1361 RA patients ascertained through the REPAIR consortium and the DANBIO registry. RA patients fulfilled the 1987 revised American College of Rheumatology (ACR) and the ACR/EULAR 2010 classification criteria. The validation cohort included 706 Dutch RA patients from the DREAM registry. The study followed the Declaration of Helsinki and study participants gave their written informed consent to participate in the study, which was approved by the ethical review committee of participant institutions. Twenty-seven single-nucleotide polymorphisms (SNPs) were selected through a literature search of relevant GWAS. Linear regression analysis adjusted for age, sex and country of origin was used to determine the association between GWAS-identified SNPs and changes in DAS28 (ΔDAS28) after 3 or 6 months of treatment. The meta-analysis of both populations was performed using a fixed effect model. Correction for multiple testing was performed using the Bonferroni method but also considering the number of inheritance models tested ( P =0.0009). To assess the role of the most interesting markers in modulating immune responses, stimulation experiments in whole blood, peripheral mononuclear cells (PBMCs) and monocyte-derived macrophages using a large number of pathogens and microbiome bacteria were performed in 408 subjects from the Human Functional Genomic Project cohort. We also evaluated the correlation of these SNPs with plasmatic levels of 108 inflammatory proteins, 7 serum steroid hormones and counts of 91 blood-derived immune cell populations.


Results: The meta-analysis of the discovery cohort and DREAM registry including 2067 RA patients treated with TNFi revealed an overall association of the LINC02549 rs7767069 SNP with a decreased drop in DAS28 that remained significant after correction for multiple testing (per-allele OR Meta =0.83, P Meta =0.000077; P Het =0.61). In addition, the meta-analysis of these large cohorts showed that each copy of the LARRC55 rs717117G allele significantly decreased the drop in DAS28 in RF-positive patients (per-allele OR Meta =0.67, P =0.00058; P Het =0.06) whereas an opposite but not significant effect was found in RF-negative subjects (per-allele OR Meta =1.38, P =0.10; P Interaction =0.00028; P Het =0.45). Interestingly, the meta-analysis also showed potentially interesting but not statistically significant overall and RF-specific associations for the MAFB rs6071980 and CNTN5 rs1813443 SNPs with ΔDAS28 (per-allele OR Meta_rs6071980 =0.84, P =0.0059; P Het =0.63 and OR Meta_rs1813443_RF+ =0.81, P =0.0059; P Het =0.69 and OR Meta_rs1813443_RF- =1.00, P =0.99; P Het =0.12; P Interaction =0.032). Although analysis of functional data is ongoing, so far, we found that carriers of the LARRC55 rs717117G allele showed decreased levels of IL6 after stimulation of PBMCs with Borrelia burgdorferi and Escherichia Coli bacteria ( P =0.00046 and 0.00044), which suggested a reduced IL6-mediated anti-inflammatory effect of this marker to worsen the response to TNFi.


Conclusion: This study confirmed the influence of the LINC02549 and LARRC55 loci to determine the response to TNFi in RA patients and a weak effect of the MAFB and CNTN5 loci that needs to be further investigated.


REFERENCES:

[1]Plenge RM et al 2011. Arthritis Rheum 63, 590-3.


ACKNOWLEDGEMENTS: We thank all participants who have agreed to participate in this study. Authors also thank María Dolores Casares, Ángeles Molina, Carmen Oloriz for the collection of Spanish samples and Hans Jurgen Hoffmann, Marianne Thomsen, Vibeke Østergaard Thomsen, Malene Rohr Andersen, Lise Lotte B. Laursen, Helle Jørgensen, Ram Benny Christian Dessau, Niels Steen Krogh, Ulla Vogel, Paal Skytt Andersen, Ivan Brandslund, Steffen Bank, Frederik Trier Møller, Nikolai Toft and Niels Møller Andersen for the participation in collection and purification of Danish samples. We also thank the Danish Departments of Rheumatology for their implication in the collection of clinical data from RA patients included in the DANBIO cohort and the Danish Rheumatologic Biobank. Likewise, we would like to thank Teun van Herwaarden for steroid hormone measurements in serum samples from subjects ascertained through the HFGP initiative.


Disclosure of Interests: None declared


Citation: Ann Rheum Dis, volume 80, supplement 1, year 2021, page 9
Session: Genomics, genetic basis of disease and functional genomics (Oral Presentations)