Sequential immunotherapy: towards cures for autoimmunity.

Despite major progress in the treatment of autoimmune diseases in the past two decades, most therapies do not cure disease and can be associated with increased risk of infection through broad suppression of the immune system. However, advances in understanding the causes of autoimmune disease and clinical data from novel therapeutic modalities such as chimeric antigen receptor T cell therapies provide evidence that it may be possible to re-establish immune homeostasis and, potentially, prolong remission or even cure autoimmune diseases. Here, we propose a 'sequential immunotherapy' framework for immune system modulation to help achieve this ambitious goal. This framework encompasses three steps: controlling inflammation; resetting the immune system through elimination of pathogenic immune memory cells; and promoting and maintaining immune homeostasis via immune regulatory agents and tissue repair. We discuss existing drugs and those in development for each of the three steps. We also highlight the importance of causal human biology in identifying and prioritizing novel immunotherapeutic strategies as well as informing their application in specific patient subsets, enabling precision medicine approaches that have the potential to transform clinical care.

High-throughput phenotyping with electronic medical record data using a common semi-supervised approach (PheCAP).

Phenotypes are the foundation for clinical and genetic studies of disease risk and outcomes. The growth of biobanks linked to electronic medical record (EMR) data has both facilitated and increased the demand for efficient, accurate, and robust approaches for phenotyping millions of patients. Challenges to phenotyping with EMR data include variation in the accuracy of codes, as well as the high level of manual input required to identify features for the algorithm and to obtain gold standard labels. To address these challenges, we developed PheCAP, a high-throughput semi-supervised phenotyping pipeline. PheCAP begins with data from the EMR, including structured data and information extracted from the narrative notes using natural language processing (NLP). The standardized steps integrate automated procedures, which reduce the level of manual input, and machine learning approaches for algorithm training. PheCAP itself can be executed in 1-2 d if all data are available; however, the timing is largely dependent on the chart review stage, which typically requires at least 2 weeks. The final products of PheCAP include a phenotype algorithm, the probability of the phenotype for all patients, and a phenotype classification (yes or no).

Correction: Epidemiology and treatment patterns of rheumatoid arthritis in a large cohort of Arab patients.

[This corrects the article DOI: 10.1371/journal.pone.0208240.].

Epidemiology and treatment patterns of rheumatoid arthritis in a large cohort of Arab patients.

OBJECTIVES: There is limited information on the epidemiology and treatment patterns of rheumatoid arthritis (RA) across the Arab region. We aim in this study to describe the demographic characteristics, clinical profile, and treatment patterns of patients of Arab ancestry with RA. METHODS: This is a cross sectional study of 895 patients with established rheumatoid arthritis enrolled from five sites (Jordan, Lebanon, Qatar, Kingdom of Saudi Arabia (KSA), and United Arab Emirates). Demographic characteristics, clinical profile, and treatment patterns are compared between the five countries. RESULTS: The majority of our patients are women, have an average disease duration of 10 years, are married and non-smokers, with completed secondary education. We report a high (>80%) ever-use of methotrexate (MTX) and steroids among our RA population, while the ever-use of disease modifying anti-rheumatic drugs (DMARDs) and TNF-inhibitors average around 67% and 33%, respectively. There are variations in RA treatment use between the five country sites. Highest utilization of steroids is identified in Jordan and KSA (p-value < 0.001), while the highest ever-use of TNF-inhibitors is reported in KSA (p-value < 0.001). CONCLUSION: Disparities in usage of RA treatments among Arab patients are noted across the five countries. National gross domestic product (GDP), as well as some other unique features in each country likely affect these. Developing treatment guidelines specific to this region could contribute in delivering standardized therapies to RA patients.

Phenome-wide association studies across large population cohorts support drug target validation.

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

Genomic atlas of the human plasma proteome.

Although plasma proteins have important roles in biological processes and are the direct targets of many drugs, the genetic factors that control inter-individual variation in plasma protein levels are not well understood. Here we characterize the genetic architecture of the human plasma proteome in healthy blood donors from the INTERVAL study. We identify 1,927 genetic associations with 1,478 proteins, a fourfold increase on existing knowledge, including trans associations for 1,104 proteins. To understand the consequences of perturbations in plasma protein levels, we apply an integrated approach that links genetic variation with biological pathway, disease, and drug databases. We show that protein quantitative trait loci overlap with gene expression quantitative trait loci, as well as with disease-associated loci, and find evidence that protein biomarkers have causal roles in disease using Mendelian randomization analysis. By linking genetic factors to diseases via specific proteins, our analyses highlight potential therapeutic targets, opportunities for matching existing drugs with new disease indications, and potential safety concerns for drugs under development.

Large Scale Metabolic Profiling identifies Novel Steroids linked to Rheumatoid Arthritis.

Recent metabolomics studies of Rheumatoid Arthritis (RA) reported few metabolites that were associated with the disease, either due to small cohort sizes or limited coverage of metabolic pathways. Our objective is to identify metabolites associated with RA and its cofounders using a new untargeted metabolomics platform. Moreover, to investigate the pathomechanism of RA by identifying correlations between RA-associated metabolites. 132 RA patients and 104 controls were analyzed for 927 metabolites. Metabolites were tested for association with RA using linear regression. OPLS-DA was used to discriminate RA patients from controls. Gaussian Graphical Models (GGMs) were used to identify correlated metabolites. 32 metabolites are identified as significantly (Bonferroni) associated with RA, including the previously reported metabolites as DHEAS, cortisol and androstenedione and extending that to a larger set of metabolites in the steroid pathway. RA classification using metabolic profiles shows a sensitivity of 91% and specificity of 88%. Steroid levels show variation among the RA patients according to the corticosteroid treatment; lowest in those taking the treatment at the time of the study, higher in those who never took the treatment, and highest in those who took it in the past. Finally, the GGM reflects metabolite relations from the steroidogenesis pathway.

A Multinational Arab Genome-Wide Association Study Identifies New Genetic Associations for Rheumatoid Arthritis.

OBJECTIVE: Genetic factors underlying susceptibility to rheumatoid arthritis (RA) in Arab populations are largely unknown. This genome-wide association study (GWAS) was undertaken to explore the generalizability of previously reported RA loci to Arab subjects and to discover new Arab-specific genetic loci. METHODS: The Genetics of Rheumatoid Arthritis in Some Arab States Study was designed to examine the genetics and clinical features of RA patients from Jordan, the Kingdom of Saudi Arabia, Lebanon, Qatar, and the United Arab Emirates. In total, >7 million single-nucleotide polymorphisms (SNPs) were tested for association with RA overall and with seropositive or seronegative RA in 511 RA cases and 352 healthy controls. In addition, replication of 15 signals was attempted in 283 RA cases and 221 healthy controls. A genetic risk score of 68 known RA SNPs was also examined in this study population. RESULTS: Three loci (HLA region, intergenic 5q13, and 17p13 at SMTNL2/GGT6) reached genome-wide significance in the analyses of association with RA and with seropositive RA, and for all 3 loci, evidence of independent replication was demonstrated. Consistent with the findings in European and East Asian populations, the association of RA with HLA-DRB1 amino acid position 11 conferred the strongest effect (P = 4.8 × 10-16 ), and a weighted genetic risk score of previously associated RA loci was found to be associated with RA (P = 3.41 × 10-5 ) and with seropositive RA (P = 1.48 × 10-6 ) in this population. In addition, 2 novel associations specific to Arab populations were found at the 5q13 and 17p13 loci. CONCLUSION: This first RA GWAS in Arab populations confirms that established HLA-region and known RA risk alleles contribute strongly to the risk and severity of disease in some Arab groups, suggesting that the genetic architecture of RA is similar across ethnic groups. Moreover, this study identified 2 novel RA risk loci in Arabs, offering further population-specific insights into the pathophysiology of RA.

Brief Report: The Role of Rare Protein-Coding Variants in Anti-Tumor Necrosis Factor Treatment Response in Rheumatoid Arthritis.

OBJECTIVE: In many rheumatoid arthritis (RA) patients, disease is controlled with anti-tumor necrosis factor (anti-TNF) biologic therapies. However, in a significant number of patients, the disease fails to respond to anti-TNF therapy. We undertook the present study to examine the hypothesis that rare and low-frequency genetic variants might influence response to anti-TNF treatment. METHODS: We sequenced the coding region of 750 genes in 1,094 RA patients of European ancestry who were treated with anti-TNF. After quality control, 690 genes were included in the analysis. We applied single-variant association and gene-based association tests to identify variants associated with anti-TNF treatment response. In addition, given the key mechanistic role of TNF, we performed gene set analyses of 27 TNF pathway genes. RESULTS: We identified 14,420 functional variants, of which 6,934 were predicted as nonsynonymous 2,136 of which were further predicted to be "damaging." Despite the fact that the study was well powered, no single variant or gene showed study-wide significant association with change in the outcome measures disease activity or European League Against Rheumatism response. Intriguingly, we observed 3 genes, of 27 with nominal signals of association (P < 0.05), that were involved in the TNF signaling pathway. However, when we performed a rigorous gene set enrichment analysis based on association P value ranking, we observed no evidence of enrichment of association at genes involved in the TNF pathway (Penrichment = 0.15, based on phenotype permutations). CONCLUSION: Our findings suggest that rare and low-frequency protein-coding variants in TNF signaling pathway genes or other genes do not contribute substantially to anti-TNF treatment response in patients with RA.

The Rheumatoid Arthritis Risk Variant CCR6DNP Regulates CCR6 via PARP-1.

Understanding the implications of genome-wide association studies (GWAS) for disease biology requires both identification of causal variants and definition of how these variants alter gene function. The non-coding triallelic dinucleotide polymorphism CCR6DNP is associated with risk for rheumatoid arthritis, and is considered likely causal because allelic variation correlates with expression of the chemokine receptor CCR6. Using transcription activator-like effector nuclease (TALEN) gene editing, we confirmed that CCR6DNP regulates CCR6. To identify the associated transcription factor, we applied a novel assay, Flanking Restriction Enhanced Pulldown (FREP), to identify specific association of poly (ADP-ribose) polymerase 1 (PARP-1) with CCR6DNP consistent with the established allelic risk hierarchy. Correspondingly, manipulation of PARP-1 expression or activity impaired CCR6 expression in several lineages. These findings show that CCR6DNP is a causal variant through which PARP-1 regulates CCR6, and introduce a highly efficient approach to interrogate non-coding genetic polymorphisms associated with human disease.

Disciplined approach to drug discovery and early development.

Our modern health care system demands therapeutic interventions that improve the lives of patients. Unfortunately, decreased productivity in therapeutics research and development (R&D) has driven drug costs up while delivering insufficient value to patients. Here, I discuss a model of translational medicine that connects four components of the early R&D pipeline-causal human biology, therapeutic modality, biomarkers of target modulation, and proof-of-concept clinical trials. Whereas the individual components of this model are not new, technological advances and a disciplined approach to integrating all four areas offer hope for improving R&D productivity.

Loci associated with N-glycosylation of human IgG are not associated with rheumatoid arthritis: a Mendelian randomisation study.

OBJECTIVES: A recent study identified 16 genetic variants associated with N-glycosylation of human IgG. Several of the genomic regions where these single nucleotide polymorphisms (SNPs) reside have also been associated with autoimmune disease (AID) susceptibility, suggesting there may be pleiotropy (genetic sharing) between loci controlling both N-glycosylation and AIDs. We investigated this by testing variants associated with levels of IgG N-glycosylation for association with rheumatoid arthritis (RA) susceptibility using a Mendelian randomisation study, and testing a subset of these variants in a less well-powered study of treatment response and severity. METHODS: SNPs showing association with IgG N-glycosylation were analysed for association with RA susceptibility in 14 361 RA cases and 43 923 controls. Five SNPs were tested for association with response to anti-tumour necrosis factor (TNF) therapy in 1081 RA patient samples and for association with radiological disease severity in 342 patients. RESULTS: Only one SNP (rs9296009) associated with N-glycosylation showed an association (p=6.92×10(-266)) with RA susceptibility, although this was due to linkage disequilibrium with causal human leukocyte antigen (HLA) variants. Four regions of the genome harboured SNPs associated with both traits (shared loci); although statistical analysis indicated that the associations observed for the two traits are independent. No SNPs showed association with response to anti-TNF therapy. One SNP rs12342831 was modestly associated with Larsen score (p=0.05). CONCLUSIONS: In a large, well-powered cohort of RA patients, we show SNPs driving levels of N-glycosylation have no association with RA susceptibility, indicating colocalisation of associated SNPs are not necessarily indicative of a shared genetic background or a role for glycosylation in disease susceptibility.

Association analysis of copy numbers of FC-gamma receptor genes for rheumatoid arthritis and other immune-mediated phenotypes.

Segmental duplications (SDs) comprise about 5% of the human genome and are enriched for immune genes. SD loci often show copy numbers variations (CNV), which are difficult to tag with genotyping methods. CNV in the Fcγ receptor region (FCGR) has been suggested to be associated with rheumatic diseases. The objective of this study was to delineate association of FCGR-CNV with rheumatoid arthritis (RA), coeliac disease and Inflammatory bowel disease incidence. We developed a method to accurately quantify CNV in SD loci based on the intensity values from the Immunochip platform and applied it to the FCGR locus. We determined the method's validity using three independent assays: segregation analysis in families, arrayCGH, and whole genome sequencing. Our data showed the presence of two separate CNVs in the FCGR locus. The first region encodes FCGR2A, FCGR3A and part of FCGR2C gene, the second encodes another part of FCGR2C, FCGR3B and FCGR2B. Analysis of CNV status in 4578 individuals with RA and 5457 controls indicated association of duplications in the FCGR3B gene in antibody-negative RA (P=0.002, OR=1.43). Deletion in FCGR3B was associated with increased risk of antibody-positive RA, consistently with previous reports (P=0.023, OR=1.23). A clear genotype-phenotype relationship was observed: CNV polymorphisms of the FCGR3A gene correlated to CD16A expression (encoded by FCGR3A) on CD8 T-cells. In conclusion, our method allows determining the CNV status of the FCGR locus, we identified association of CNV in FCGR3B to RA and showed a functional relationship between CNV in the FCGR3A gene and CD16A expression.

Somatic Variation of T-Cell Receptor Genes Strongly Associate with HLA Class Restriction.

Every person carries a vast repertoire of CD4+ T-helper cells and CD8+ cytotoxic T cells for a healthy immune system. Somatic VDJ recombination at genomic loci that encode the T-cell receptor (TCR) is a key step during T-cell development, but how a single T cell commits to become either CD4+ or CD8+ is poorly understood. To evaluate the influence of TCR sequence variation on CD4+/CD8+ lineage commitment, we sequenced rearranged TCRs for both α and ß chains in naïve T cells isolated from healthy donors and investigated gene segment usage and recombination patterns in CD4+ and CD8+ T-cell subsets. Our data demonstrate that most V and J gene segments are strongly biased in the naïve CD4+ and CD8+ subsets with some segments increasing the odds of being CD4+ (or CD8+) up to five-fold. These V and J gene associations are highly reproducible across individuals and independent of classical HLA genotype, explaining ~11% of the observed variance in the CD4+ vs. CD8+ propensity. In addition, we identified a strong independent association of the electrostatic charge of the complementarity determining region 3 (CDR3) in both α and ß chains, where a positively charged CDR3 is associated with CD4+ lineage and a negatively charged CDR3 with CD8+ lineage. Our findings suggest that somatic variation in different parts of the TCR influences T-cell lineage commitment in a predominantly additive fashion. This notion can help delineate how certain structural features of the TCR-peptide-HLA complex influence thymic selection.

Methods to Develop an Electronic Medical Record Phenotype Algorithm to Compare the Risk of Coronary Artery Disease across 3 Chronic Disease Cohorts.

BACKGROUND: Typically, algorithms to classify phenotypes using electronic medical record (EMR) data were developed to perform well in a specific patient population. There is increasing interest in analyses which can allow study of a specific outcome across different diseases. Such a study in the EMR would require an algorithm that can be applied across different patient populations. Our objectives were: (1) to develop an algorithm that would enable the study of coronary artery disease (CAD) across diverse patient populations; (2) to study the impact of adding narrative data extracted using natural language processing (NLP) in the algorithm. Additionally, we demonstrate how to implement CAD algorithm to compare risk across 3 chronic diseases in a preliminary study. METHODS AND RESULTS: We studied 3 established EMR based patient cohorts: diabetes mellitus (DM, n = 65,099), inflammatory bowel disease (IBD, n = 10,974), and rheumatoid arthritis (RA, n = 4,453) from two large academic centers. We developed a CAD algorithm using NLP in addition to structured data (e.g. ICD9 codes) in the RA cohort and validated it in the DM and IBD cohorts. The CAD algorithm using NLP in addition to structured data achieved specificity >95% with a positive predictive value (PPV) 90% in the training (RA) and validation sets (IBD and DM). The addition of NLP data improved the sensitivity for all cohorts, classifying an additional 17% of CAD subjects in IBD and 10% in DM while maintaining PPV of 90%. The algorithm classified 16,488 DM (26.1%), 457 IBD (4.2%), and 245 RA (5.0%) with CAD. In a cross-sectional analysis, CAD risk was 63% lower in RA and 68% lower in IBD compared to DM (p<0.0001) after adjusting for traditional cardiovascular risk factors. CONCLUSIONS: We developed and validated a CAD algorithm that performed well across diverse patient populations. The addition of NLP into the CAD algorithm improved the sensitivity of the algorithm, particularly in cohorts where the prevalence of CAD was low. Preliminary data suggest that CAD risk was significantly lower in RA and IBD compared to DM.

TRAF1/C5 but not PTPRC variants are potential predictors of rheumatoid arthritis response to anti-tumor necrosis factor therapy.

BACKGROUND: The aim of our work was to replicate, in a Southern European population, the association reported in Northern populations between PTPRC locus and response to anti-tumor necrosis factor (anti-TNF) treatment in rheumatoid arthritis (RA). We also looked at associations between five RA risk alleles and treatment response. METHODS: We evaluated associations between anti-TNF treatment responses assessed by DAS28 change and by EULAR response at six months in 383 Portuguese patients. Univariate and multivariate linear and logistic regression analyses were performed. In a second step to confirm our findings, we pooled our population with 265 Spanish patients. RESULTS: No association was found between PTPRC rs10919563 allele and anti-TNF treatment response, neither in Portuguese modeling for several clinical variables nor in the overall population combining Portuguese and Spanish patients. The minor allele for RA susceptibility, rs3761847 SNP in TRAF1/C5 region, was associated with a poor response in linear and logistic univariate and multivariate regression analyses. No association was observed with the other allellic variants. Results were confirmed in the pooled analysis. CONCLUSION: This study did not replicate the association between PTPRC and the response to anti-TNF treatment in our Southern European population. We found that TRAF1/C5 risk RA variants potentially influence anti-TNF treatment response.

Accuracy of the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) as a research tool for identification of patients with uveitis and scleritis.

PURPOSE: To report on the accuracy of the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes for identifying patients with polymyalgia rheumatica (PMR) and concurrent noninfectious inflammatory ocular conditions in a large healthcare organization database. METHODS: Queries for patients with PMR and uveitis or scleritis were executed in two general teaching hospitals' databases. Patients with ocular infections or other rheumatologic conditions were excluded. Patients with PMR and ocular inflammation were identified, and medical records were reviewed to confirm accuracy. RESULTS: The query identified 10,697 patients with the ICD-9-CM code for PMR and 4154 patients with the codes for noninfectious inflammatory ocular conditions. The number of patients with both PMR and noninfectious uveitis or scleritis by ICD-9-CM codes was 66. On detailed review of the charts of these 66 patients, 31 (47%) had a clinical diagnosis of PMR, 43 (65%) had noninfectious uveitis or scleritis, and only 20 (30%) had PMR with concurrent noninfectious uveitis or scleritis confirmed based on clinical notes. CONCLUSIONS: While the use of ICD-9-CM codes has been validated for medical research of common diseases, our results suggest that ICD-9-CM codes may be of limited value for epidemiological investigations of diseases which can be more difficult to diagnose. The ICD-9-CM codes for rarer diseases (PMR, uveitis and scleritis) did not reflect the true clinical problem in a large proportion of our patients. This is particularly true when coding is performed by physicians outside the area of specialty of the diagnosis.