Phenome‐Wide Association Study of Autoantibodies to Citrullinated and Noncitrullinated Epitopes in Rheumatoid Arthritis

• Published in: Arthritis & rheumatology (Hoboken, N.J.) Vol. 69; no. 4; pp. 742 - 749
• Main Authors: Liao, Katherine P., Sparks, Jeffrey A., Hejblum, Boris P., Kuo, I‐Hsin, Cui, Jing, Lahey, Lauren J., Cagan, Andrew, Gainer, Vivian S., Liu, Weidong, Cai, T. Tony, Sokolove, Jeremy, Cai, Tianxi
• Format: Journal Article
• Language: English
• Published: United States Wiley Subscription Services, Inc 01.04.2017; Wiley; John Wiley and Sons Inc
• Subjects: Antigens; Applications; Arthritis; Arthritis, Rheumatoid - genetics; Arthritis, Rheumatoid - immunology; Autoantibodies; Autoantibodies - genetics; Bioinformatics; Bronchopneumonia; Citrulline; Codes; Electronic health records; Electronic medical records; Epitopes; Female; Fibrinogen; Fibronectin; Genome-Wide Association Study; Human health and pathology; Humans; Identification methods; Life Sciences; Male; Medical records; Middle Aged; Obliterative bronchiolitis; Patients; Peptides, Cyclic - immunology; Phenotype; Rheumatoid Arthritis; Rhumatology and musculoskeletal system; Statistics; Autoimmunity; Pulmonary fibrosis; Anti-CCP; Rheumatoid arthritis
• ISSN: 2326-5191; 2326-5205; 2326-5205; 2326-5191
• DOI: 10.1002/art.39974

Objective Patients with rheumatoid arthritis (RA) develop autoantibodies against a spectrum of antigens, but the clinical significance of these autoantibodies is unclear. Using a phenome‐wide association study (PheWAS) approach, we examined the association between autoantibodies and clinical subphenotypes of RA. Methods This study was conducted in a cohort of RA patients identified from the electronic medical records (EMRs) of 2 tertiary care centers. Using a published multiplex bead assay, we measured 36 autoantibodies targeting epitopes implicated in RA. We extracted all International Classification of Diseases, Ninth Revision (ICD‐9) codes for each subject and grouped them into disease categories (PheWAS codes), using a published method. We tested for the association of each autoantibody (grouped by the targeted protein) with PheWAS codes. To determine significant associations (at a false discovery rate [FDR] of ≤0.1), we reviewed the medical records of 50 patients with each PheWAS code to determine positive predictive values (PPVs). Results We studied 1,006 RA patients; the mean ± SD age of the patients was 61.0 ± 12.9 years, and 79.0% were female. A total of 3,568 unique ICD‐9 codes were grouped into 625 PheWAS codes; the 206 PheWAS codes with a prevalence of ≥3% were studied. Using the PheWAS method, we identified 24 significant associations of autoantibodies to epitopes at an FDR of ≤0.1. The associations that were strongest and had the highest PPV for the PheWAS code were autoantibodies against fibronectin and obesity (P = 6.1 × 10−4, PPV 100%), and that between fibrinogen and pneumonopathy (P = 2.7 × 10−4, PPV 96%). Pneumonopathy codes included diagnoses for cryptogenic organizing pneumonia and obliterative bronchiolitis. Conclusion We demonstrated application of a bioinformatics method, the PheWAS, to screen for the clinical significance of RA‐related autoantibodies. Using the PheWAS approach, we identified potentially significant links between variations in the levels of autoantibodies and comorbidities of interest in RA.