Publications by Year: 2020

Huizinga TWJ, Holers MV, Anolik J, Brenner MB, Buckley CD, Bykerk V, Connolly SE, Deane KD, Guo J, Hodge M, Hoffmann S, Nestle F, Pitzalis C, Raychaudhuri S, Yamamoto K, Li Z, Klareskog L. Disruptive innovation in rheumatology: new networks of global public–private partnerships are needed to take advantage of scientific progress [Internet]. Annals of the Rheumatic Diseases 2020;79(5):553-555. Publisher's Version
Amariuta T, Luo Y, Knevel R, Okada Y, Raychaudhuri S. Advances in genetics toward identifying pathogenic cell states of rheumatoid arthritis [Internet]. Immunological Reviews 2020;294(1):188-204. Publisher's VersionAbstract
Rheumatoid arthritis (RA) risk has a large genetic component (~60%) that is still not fully understood. This has hampered the design of effective treatments that could promise lifelong remission. RA is a polygenic disease with 106 known genome-wide significant associated loci and thousands of small effect causal variants. Our current understanding of RA risk has suggested cell-type-specific contexts for causal variants, implicating CD4 + effector memory T cells, as well as monocytes, B cells and stromal fibroblasts. While these cellular states and categories are still mechanistically broad, future studies may identify causal cell subpopulations. These efforts are propelled by advances in single cell profiling. Identification of causal cell subpopulations may accelerate therapeutic intervention to achieve lifelong remission.
Svensson MND, Zoccheddu M, Yang S, Nygaard G, Secchi C, Doody KM, Slowikowski K, Mizoguchi F, Humby F, Hands R, Santelli E, Sacchetti C, Wakabayashi K, Wu DJ, Barback C, Ai R, Wang W, Sims GP, Mydel P, Kasama T, Boyle DL, Galimi F, Vera D, Tremblay ML, Raychaudhuri S, Brenner MB, Firestein GS, Pitzalis C, Ekwall A-KH, Stanford SM, Bottini N. Synoviocyte-targeted therapy synergizes with TNF inhibition in arthritis reversal [Internet]. Science Advances 2020;6(26):eaba4353. Publisher's VersionAbstract
Fibroblast-like synoviocytes (FLS) are joint-lining cells that promote rheumatoid arthritis (RA) pathology. Current disease-modifying antirheumatic agents (DMARDs) operate through systemic immunosuppression. FLS-targeted approaches could potentially be combined with DMARDs to improve control of RA without increasing immunosuppression. Here, we assessed the potential of immunoglobulin-like domains 1 and 2 (Ig1&2), a decoy protein that activates the receptor tyrosine phosphatase sigma (PTPRS) on FLS, for RA therapy. We report that PTPRS expression is enriched in synovial lining RA FLS and that Ig1&2 reduces migration of RA but not osteoarthritis FLS. Administration of an Fc-fusion Ig1&2 attenuated arthritis in mice without affecting innate or adaptive immunity. Furthermore, PTPRS was down-regulated in FLS by tumor necrosis factor (TNF) via a phosphatidylinositol 3-kinase–mediated pathway, and TNF inhibition enhanced PTPRS expression in arthritic joints. Combination of ineffective doses of TNF inhibitor and Fc-Ig1&2 reversed arthritis in mice, providing an example of synergy between FLS-targeted and immunosuppressive DMARD therapies.
Krebs K, Bovjin J, Zheng N, Lepamets M, Censin JC, Jurgenson T, Sarg D, Abner E, Laisk T, Luo Y, Skotte L, Geller F, Feenstra B, Wang W, Auton A, Auton A, Raychaudhuri S, Esko T, Metspalu A, Laur S, Roden DM, Wei W-Q, Holmes MV, Lindgren CM, Phillips EJ, Magi R, Milani L, Fadista J. Genome-wide Study Identifies Association between HLA-B 55:01 and Self-Reported Penicillin Allergy [Internet]. American Journal of Human Genetics 2020;107(4):612-621. Publisher's VersionAbstract
Hypersensitivity reactions to drugs are often unpredictable and can be life threatening, underscoring a need for understanding their underlying mechanisms and risk factors. The extent to which germline genetic variation influences the risk of commonly reported drug allergies such as penicillin allergy remains largely unknown. We extracted data from the electronic health records of more than 600,000 participants from the UK, Estonian, and Vanderbilt University Medical Center's BioVU biobanks to study the role of genetic variation in the occurrence of self-reported penicillin hypersensitivity reactions. We used imputed SNP to HLA typing data from these cohorts to further fine map the human leukocyte antigen (HLA) association and replicated our results in 23andMe's research cohort involving a total of 1.12 million individuals. Genome-wide meta-analysis of penicillin allergy revealed two loci, including one located in the HLA region on chromosome 6. This signal was further fine-mapped to the HLA-B55:01 allele (OR 1.41 95% CI 1.33-1.49, p value 2.04 × 10-31) and confirmed by independent replication in 23andMe's research cohort (OR 1.30 95% CI 1.25-1.34, p value 1.00 × 10-47). The lead SNP was also associated with lower lymphocyte counts and in silico follow-up suggests a potential effect on T-lymphocytes at HLA-B55:01. We also observed a significant hit in PTPN22 and the GWAS results correlated with the genetics of rheumatoid arthritis and psoriasis. We present robust evidence for the role of an allele of the major histocompatibility complex (MHC) I gene HLA-B in the occurrence of penicillin allergy.
Bonfa E, Gossec L, Isenberg DA, Li Z, Raychaudhuri S. How COVID-19 is changing rheumatology clinical practice [Internet]. Nature Reviews. Rheumatology 2020; Publisher's VersionAbstract
The emergence of COVID-19 in early 2020 led to unprecedented changes to rheumatology clinical practice worldwide, including the closure of research laboratories, the restructuring of hospitals and the rapid transition to virtual care. As governments sought to slow and contain the spread of the disease, rheumatologists were presented with the difficult task of managing risks, to their patients as well as to themselves, while learning and implementing new systems for remote health care. Consequently, the COVID-19 pandemic led to a transformation in health infrastructures and telemedicine that could become powerful tools for rheumatologists, despite having some limitations. In this Viewpoint, five experts from different regions discuss their experiences of the pandemic, including the most challenging aspects of this unexpected transition, the advantages and limitations of virtual visits, and potential opportunities going forward.
Choi W, Luo Y, Raychaudhuri S, Han B. HATK: HLA analysis toolkit [Internet]. Bioinformatics 2020; Publisher's VersionAbstract
Fine-mapping human leukocyte antigen (HLA) genes involved in disease susceptibility to individual alleles or amino acid residues has been challenging. Using information regarding HLA alleles obtained from HLA typing, HLA imputation or HLA inference, our software expands the alleles to amino acid sequences using the most recent IMGT/HLA database and prepares a dataset suitable for fine-mapping analysis. Our software also provides useful functionalities, such as various association tests, visualization tools and nomenclature conversion.
Orange DE, Yao V, Sawicka K, Fak J, Frank MO, Parveen S, Blachere NE, Hale C, Zhang F, Raychaudhuri S, Troyanskaya OG, Darnell RB. RNA Identification of PRIME Cells Predicting Rheumatoid Arthritis Flares [Internet]. New England Journal of Medicine 2020;383(3):218-228. Publisher's VersionAbstract

Background: Rheumatoid arthritis, like many inflammatory diseases, is characterized by episodes of quiescence and exacerbation (flares). The molecular events leading to flares are unknown.

Methods: We established a clinical and technical protocol for repeated home collection of blood in patients with rheumatoid arthritis to allow for longitudinal RNA sequencing (RNA-seq). Specimens were obtained from 364 time points during eight flares over a period of 4 years in our index patient, as well as from 235 time points during flares in three additional patients. We identified transcripts that were differentially expressed before flares and compared these with data from synovial single-cell RNA-seq. Flow cytometry and sorted-blood-cell RNA-seq in additional patients were used to validate the findings.

Results: Consistent changes were observed in blood transcriptional profiles 1 to 2 weeks before a rheumatoid arthritis flare. B-cell activation was followed by expansion of circulating CD45-CD31-PDPN+ preinflammatory mesenchymal, or PRIME, cells in the blood from patients with rheumatoid arthritis; these cells shared features of inflammatory synovial fibroblasts. Levels of circulating PRIME cells decreased during flares in all 4 patients, and flow cytometry and sorted-cell RNA-seq confirmed the presence of PRIME cells in 19 additional patients with rheumatoid arthritis.

Conclusions: Longitudinal genomic analysis of rheumatoid arthritis flares revealed PRIME cells in the blood during the period before a flare and suggested a model in which these cells become activated by B cells in the weeks before a flare and subsequently migrate out of the blood into the synovium. (Funded by the National Institutes of Health and others.).

Knevel R, le Cessie S, Terao CC, Slowikowski K, Cui J, Huizinga TWJ, Costenbader KH, Liao KP, Karlson EW, Raychaudhuri S. Using genetics to prioritize diagnoses for rheumatology outpatients with inflammatory arthritis [Internet]. Science Translational Medicine 2020;12(545) Full TextAbstract
Multiple slowly progressing diseases initially present with inflammatory arthritis, and it can be difficult to clinically differentiate these conditions. Knevel et al. show that genetic data could be used to triage inflammatory arthritis–causing diagnoses at a patient’s first visit, improving the likelihood of a correct initial diagnosis and potentially expediting appropriate treatment. Their genetic diagnostic tool, here optimized for rheumatic disease diagnosis, could, in principle, be calibrated for other phenotypically similar diseases with different underlying genetics.It is challenging to quickly diagnose slowly progressing diseases. To prioritize multiple related diagnoses, we developed G-PROB (Genetic Probability tool) to calculate the probability of different diseases for a patient using genetic risk scores. We tested G-PROB for inflammatory arthritis–causing diseases (rheumatoid arthritis, systemic lupus erythematosus, spondyloarthropathy, psoriatic arthritis, and gout). After validating on simulated data, we tested G-PROB in three cohorts: 1211 patients identified by International Classification of Diseases (ICD) codes within the eMERGE database, 245 patients identified through ICD codes and medical record review within the Partners Biobank, and 243 patients first presenting with unexplained inflammatory arthritis and with final diagnoses by record review within the Partners Biobank. Calibration of G-probabilities with disease status was high, with regression coefficients from 0.90 to 1.08 (1.00 is ideal). G-probabilities discriminated true diagnoses across the three cohorts with pooled areas under the curve (95% CI) of 0.69 (0.67 to 0.71), 0.81 (0.76 to 0.84), and 0.84 (0.81 to 0.86), respectively. For all patients, at least one disease could be ruled out, and in 45% of patients, a likely diagnosis was identified with a 64% positive predictive value. In 35% of cases, the clinician’s initial diagnosis was incorrect. Initial clinical diagnosis explained 39% of the variance in final disease, which improved to 51% (P < 0.0001) after adding G-probabilities. Converting genotype information before a clinical visit into an interpretable probability value for five different inflammatory arthritides could potentially be used to improve the diagnostic efficiency of rheumatic diseases in clinical practice.
Asgari S, Luo Y, Belbin GM, Bartell E, Calderon R, Slowikowski K, Contreras C, Yataco R, Galea JT, Jimenez J, Coit JM, Farroñay C, Nazarian RM, O’Connor TD, Dietz HC, Hirschhorn J, Guio H, Lecca L, Kenny EE, Freeman E, Murray MB, Raychaudhuri S. A positively selected, common, missense variant in FBN1 confers a 2.2 centimeter reduction of height in the Peruvian population [Internet]. Nature 2020;582(7811):234-239. NCBI LinkAbstract
Peruvians are among the shortest people in the world. To understand the genetic basis of short stature in Peru, we examined an ethnically diverse group of Peruvians and identified a novel, population-specific, missense variant in FBN1 (E1297G) that is significantly associated with lower height in the Peruvian population. Each copy of the minor allele (frequency = 4.7%) reduces height by 2.2 cm (4.4 cm in homozygous individuals). This is the largest effect size known for a common height-associated variant. This variant shows strong evidence of positive selection within the Peruvian population and is significantly more frequent in Native American populations from coastal regions of Peru compared to populations from the Andes or the Amazon, suggesting that short stature in Peruvians is the result of adaptation to the coastal environment.One Sentence Summary A mutation found in Peruvians has the largest known effect on height for a common variant. This variant is specific to Native American ancestry.
Ishigaki K, Akiyama M, Kanai M, Takahashi A, Kawakami E, Sugishita H, Sakaue S, Matoba N, Low S-K, Okada Y, Terao C, Amariuta T, Gazal S, Kochi Y, Horikoshi M, Suzuki K, Ito K, Momozawa Y, Hirata M, Matsuda K, Ikeda M, Iwata N, Ikegawa S, Kou I, Tanaka T, Nakagawa H, Suzuki A, Hirota T, Tamari M, Chayama K, Miki D, Mori M, Nagayama S, Daigo Y, Miki Y, Katagiri T, Ogawa O, Obara W, Ito H, Yoshida T, Imoto I, Takahashi T, Tanikawa C, Suzuki T, Sinozaki N, Minami S, Yamaguchi H, Asai S, Takahashi Y, Yamaji K, Takahashi K, Fujioka T, Takata R, Yanai H, Masumoto A, Koretsune Y, Kutsumi H, Higashiyama M, Murayama S, Minegishi N, Suzuki K, Tanno K, Shimizu A, Yamaji T, Iwasaki M, Sawada N, Uemura H, Tanaka K, Naito M, Sasaki M, Wakai K, Tsugane S, Yamamoto M, Yamamoto K, Murakami Y, Nakamura Y, Raychaudhuri S, Inazawa J, Yamauchi T, Kadowaki T, Kubo M, Kamatani Y. Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases [Internet]. Nature Genetics 2020;52(7):669-679. Publisher's VersionAbstract
The overwhelming majority of participants in current genetic studies are of European ancestry1–3, limiting our genetic understanding of complex disease in non-European populations. To address this, we aimed to elucidate polygenic disease biology in the East Asian population by conducting a genome-wide association study (GWAS) with 212,453 Japanese individuals across 42 diseases. We detected 383 independent signals in 331 loci for 30 diseases, among which 45 loci were novel (P < 5 × 10-8). Compared with known variants, novel variants have lower frequency in European populations but comparable frequency in East Asian populations, suggesting the advantage of this study in discovering these novel variants. Three novel signals were in linkage disequilibrium (r2 > 0.6) with missense variants which are monomorphic in European populations (1000 Genomes Project) including rs11235604(p.R220W of ATG16L2, a autophagy-related gene) associated with coronary artery disease. We further investigated enrichment of heritability within 2,868 annotations of genome-wide transcription factor occupancy, andidentified 378 significant enrichments across nine diseases (FDR < 0.05) (e.g. NF-κB for immune-related diseases). This large-scale GWAS in a Japanese population provides insights into the etiology of common complex diseases and highlights the importance of performing GWAS in non-European populations.
Gutierrez-Arcelus M#, Baglaenko Y#, Arora J, Hannes S, Luo Y, Amariuta T, Teslovich N, Rao DA, Ermann J, Jonsson AH, for Consortium NHLBIT-OPM (TOPM), Navarrete C, Rich SS, Taylor KD, Rotter JI, Gregersen PK, Esko T, Brenner MB, Raychaudhuri S. Allele-specific expression changes dynamically during T cell activation in HLA and other autoimmune loci [Internet]. Nature Genetics 2020;52:247-253. Publisher's VersionAbstract
Genetic studies have revealed that autoimmune susceptibility variants are over-represented in memory CD4+ T cell regulatory elements1-3. Understanding how genetic variation affects gene expression in different T cell physiological states is essential for deciphering genetic mechanisms of autoimmunity4,5. Here, we characterized the dynamics of genetic regulatory effects at eight time points during memory CD4+ T cell activation with high-depth RNA-seq in healthy individuals. We discovered widespread, dynamic allele-specific expression across the genome, where the balance of alleles changes over time. These genes were enriched fourfold within autoimmune loci. We found pervasive dynamic regulatory effects within six HLA genes. HLA-DQB1 alleles had one of three distinct transcriptional regulatory programs. Using CRISPR-Cas9 genomic editing we demonstrated that a promoter variant is causal for T cell-specific control of HLA-DQB1 expression. Our study shows that genetic variation in cis-regulatory elements affects gene expression in a manner dependent on lymphocyte activation status, contributing to the interindividual complexity of immune responses.
Wei K#, Korsunsky I#, Marshall JL, Gao A, Watts GFM, Major T, Croft AP, Watts J, Blazar P, Lange J, Thornhill T, Filer A, Raza K, Donlin LT, Accelerating Medicines Partnership-Rheumatoid arthritis/Systemic Lupus Erythematosus (AMP RA/SLE), Siebel CW, Buckley CD, Raychaudhuri S*, Brenner* MB. Notch signalling drives synovial fibroblast identity and arthritis pathology [Internet]. Nature 2020;582(7811):259-264. Publisher's Version
Slowikowski K#, Nguyen# HN, Noss EH, Simmons DP, Mizoguchi F, Watts GFM, Gurish MF, Brenner MB*, Raychaudhuri S*. CUX1 and IκBζ (NFKBIZ) mediate the synergistic inflammatory response to TNF and IL-17A in stromal fibroblasts. [Internet]. Proc Natl Acad Sci. 2020;117(10):5532-5541. Publisher's VersionAbstract
The role of stromal fibroblasts in chronic inflammation is unfolding. In rheumatoid arthritis, leukocyte-derived cytokines TNF and IL-17A work together, activating fibroblasts to become a dominant source of the hallmark cytokine IL-6. However, IL-17A alone has minimal effect on fibroblasts. To identify key mediators of the synergistic response to TNF and IL-17A in human synovial fibroblasts, we performed time series, dose-response, and gene-silencing transcriptomics experiments. Here we show that in combination with TNF, IL-17A selectively induces a specific set of genes mediated by factors including cut-like homeobox 1 (CUX1) and IκBζ (NFKBIZ). In the promoters of CXCL1, CXCL2, and CXCL3, we found a putative CUX1-NF-κB binding motif not found elsewhere in the genome. CUX1 and NF-κB p65 mediate transcription of these genes independent of LIFR, STAT3, STAT4, and ELF3. Transcription of NFKBIZ, encoding the atypical IκB factor IκBζ, is IL-17A dose-dependent, and IκBζ only mediates the transcriptional response to TNF and IL-17A, but not to TNF alone. In fibroblasts, IL-17A response depends on CUX1 and IκBζ to engage the NF-κB complex to produce chemoattractants for neutrophil and monocyte recruitment.