ObjectiveRecently, a genome-wide association study identified an association between RA-associated interstitial lung disease (ILD) and RPA3-UMAD1 rs12702634 in the Japanese population, especially... Show moreObjectiveRecently, a genome-wide association study identified an association between RA-associated interstitial lung disease (ILD) and RPA3-UMAD1 rs12702634 in the Japanese population, especially for patients with a usual interstitial pneumonia (UIP) pattern. We aimed to replicate this association in a European population and test for interaction with MUC5B rs35705950.MethodsIn this genetic case–control association study, patients with RA and ILD and controls with RA and no ILD were included from France, the USA and the Netherlands. Only cases and controls from European genetic ancestries determined by principal components analysis were included in the analyses. RA was defined by the 1987 ACR or 2010 ACR/EULAR criteria and ILD by chest high-resolution CT scan, except in the control dataset from the Netherlands, where the absence of ILD was determined by chart review. Patients were genotyped for RPA3-UMAD1 rs12702634 and MUC5B rs35705950. Associations were tested using logistic regression adjusted for sex, age at RA onset, age at ILD onset or at certified absence of ILD, tobacco smoking status and country of origin.ResultsAmong the 883 patients included, 322 were RA-ILD cases (36.5%). MUC5B rs35705950 was strongly associated with RA-ILD in all datasets {combined adjusted odds ratio [OR] 2.9 [95% CI 2.1, 3.9], P = 1.1 × 10−11. No association between RPA3-UMAD1 rs12702634 and RA-ILD was observed [combined OR 1.2 (95% CI 0.8, 1.6), P = 0.31. No interaction was found between RPA3-UMAD1 rs12702634 and MUC5B rs35705950 (P = 0.70).ConclusionOur findings did not support a contribution of RPA3-UMAD1 rs12702634 to the overall RA-ILD susceptibility in the European population. Show less
Multi-ancestry genome-wide association analyses identify 124 risk loci for rheumatoid arthritis, of which 34 are novel. A polygenic risk score based on multi-ancestry data showed comparable... Show moreMulti-ancestry genome-wide association analyses identify 124 risk loci for rheumatoid arthritis, of which 34 are novel. A polygenic risk score based on multi-ancestry data showed comparable performance between populations of European and East Asian ancestries.Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 x 10(-8)), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA. Show less
Author summaryArteriviruses are a rapidly expanding family of positive-stranded RNA viruses, which includes economically important veterinary pathogens like equine arteritis virus (EAV) and two... Show moreAuthor summaryArteriviruses are a rapidly expanding family of positive-stranded RNA viruses, which includes economically important veterinary pathogens like equine arteritis virus (EAV) and two species of porcine reproductive and respiratory syndrome virus (PRRSV-1 and PRRSV-2). In our previous studies, we uncovered an unprecedented arterivirus gene expression mechanism: a highly efficient -2 programmed ribosomal frameshift (PRF) that is controlled by an interaction of viral protein nsp1ss with specific RNA sequences and host poly(C) binding proteins. It is used by PRRSVs, and most other arteriviruses, to efficiently produce a previously unknown nonstructural protein variant, nsp2TF. In this study, we demonstrate that PRRSV nsp2TF interacts with the two major arteriviral envelope proteins, GP5 and M, whose heterodimerization in the secretory pathway is a key step in envelope protein trafficking and virus assembly. Our findings suggest that nsp2TF promotes arterivirus assembly by antagonizing the ubiquitination-dependent proteasomal degradation of GP5 and M proteins. This mechanism is based on the DUB activity of the PLP2 protease domain located within the N-terminal region of nsp2TF. To our knowledge, this is the first study to demonstrate that viruses can express a DUB that functions specifically to counteract the ubiquitination and degradation of key viral structural proteins.Arteriviruses are enveloped positive-strand RNA viruses that assemble and egress using the host cell's exocytic pathway. In previous studies, we demonstrated that most arteriviruses use a unique -2 ribosomal frameshifting mechanism to produce a C-terminally modified variant of their nonstructural protein 2 (nsp2). Like full-length nsp2, the N-terminal domain of this frameshift product, nsp2TF, contains a papain-like protease (PLP2) that has deubiquitinating (DUB) activity, in addition to its role in proteolytic processing of replicase polyproteins. In cells infected with porcine reproductive and respiratory syndrome virus (PRRSV), nsp2TF localizes to compartments of the exocytic pathway, specifically endoplasmic reticulum-Golgi intermediate compartment (ERGIC) and Golgi complex. Here, we show that nsp2TF interacts with the two major viral envelope proteins, the GP5 glycoprotein and membrane (M) protein, which drive the key process of arterivirus assembly and budding. The PRRSV GP5 and M proteins were found to be poly-ubiquitinated, both in an expression system and in cells infected with an nsp2TF-deficient mutant virus. In contrast, ubiquitinated GP5 and M proteins did not accumulate in cells infected with the wild-type, nsp2TF-expressing virus. Further analysis implicated the DUB activity of the nsp2TF PLP2 domain in deconjugation of ubiquitin from GP5/M proteins, thus antagonizing proteasomal degradation of these key viral structural proteins. Our findings suggest that nsp2TF is targeted to the exocytic pathway to reduce proteasome-driven turnover of GP5/M proteins, thus promoting the formation of GP5-M dimers that are critical for arterivirus assembly. Show less
Objective To identify interactions between genetic factors and current or recent smoking in relation to risk of developing systemic lupus erythematosus (SLE). Methods For the study, 673 patients... Show moreObjective To identify interactions between genetic factors and current or recent smoking in relation to risk of developing systemic lupus erythematosus (SLE). Methods For the study, 673 patients with SLE (diagnosed according to the American College of Rheumatology 1997 updated classification criteria) were matched by age, sex, and race (first 3 genetic principal components) to 3,272 control subjects without a history of connective tissue disease. Smoking status was classified as current smoking/having recently quit smoking within 4 years before diagnosis (or matched index date for controls) versus distant past/never smoking. In total, 86 single-nucleotide polymorphisms and 10 classicHLAalleles previously associated with SLE were included in a weighted genetic risk score (wGRS), with scores dichotomized as either low or high based on the median value in control subjects (low wGRS being defined as less than or equal to the control median; high wGRS being defined as greater than the control median). Conditional logistic regression models were used to estimate both the risk of SLE and risk of anti-double-stranded DNA autoantibody-positive (dsDNA+) SLE. Additive interactions were assessed using the attributable proportion (AP) due to interaction, and multiplicative interactions were assessed using a chi-square test (with 1 degree of freedom) for the wGRS and for individual risk alleles. Separate repeated analyses were carried out among subjects of European ancestry only. Results The mean +/- SD age of the SLE patients at the time of diagnosis was 36.4 +/- 15.3 years. Among the 673 SLE patients included, 92.3% were female and 59.3% were dsDNA+. Ethnic distributions were as follows: 75.6% of European ancestry, 4.5% of Asian ancestry, 11.7% of African ancestry, and 8.2% classified as other ancestry. A high wGRS (odds ratio [OR] 2.0,P= 1.0 x 10(-51)versus low wGRS) and a status of current/recent smoking (OR 1.5,P= 0.0003 versus distant past/never smoking) were strongly associated with SLE risk, with significant additive interaction (AP 0.33,P= 0.0012), and associations with the risk of anti-dsDNA+ SLE were even stronger. No significant multiplicative interactions with the total wGRS (P= 0.58) or with theHLA-only wGRS (P= 0.06) were found. Findings were similar in analyses restricted to only subjects of European ancestry. Conclusion The strong additive interaction between an updated SLE genetic risk score and current/recent smoking suggests that smoking may influence specific genes in the pathogenesis of SLE. Show less
Knevel, R.; Cessie, S. le; Terao, C.C.; Slowikowski, K.; Cui, J.; Huizinga, T.W.J.; ... ; Raychaudhuri, S. 2020
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... Show moreIt 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. Show less
BackgroundPrevious studies of radiological damage in rheumatoid arthritis (RA) have used candidate-gene approaches, or evaluated single genome-wide association studies (GWAS). We undertook the... Show moreBackgroundPrevious studies of radiological damage in rheumatoid arthritis (RA) have used candidate-gene approaches, or evaluated single genome-wide association studies (GWAS). We undertook the first meta-analysis of GWAS of RA radiological damage to: (1) identify novel genetic loci for this trait; and (2) test previously validated variants.MethodsSeven GWAS (2,775 RA cases, of a range of ancestries) were combined in a meta-analysis. Radiological damage was assessed using modified Larsen scores, Sharp van Der Heijde scores, and erosive status. Single nucleotide polymophsim (SNP) associations with radiological damage were tested at a single time-point using regression models. Primary analyses included age and disease duration as covariates. Secondary analyses also included rheumatoid factor (RF). Meta-analyses were undertaken in trans-ethnic and European-only cases.ResultsIn the trans-ethnic primary meta-analysis, one SNP (rs112112734) in close proximity to HLA-DRB1, and strong linkage disequilibrium with the shared-epitope, attained genome-wide significance (P = 4.2x10(-8)). In the secondary analysis (adjusting for RF) the association was less significant (P = 1.7x10(-6)). In both trans-ethnic primary and secondary meta-analyses 14 regions contained SNPs with associations reaching P<5x10(-6); in the European primary and secondary analyses 13 and 10 regions contained SNPs reaching P<5x10(-6), respectively. Of the previously validated SNPs for radiological progression, only rs660895 (tagging HLA-DRB1*04:01) attained significance (P = 1.6x10(-5)) and had a consistent direction of effect across GWAS.ConclusionsOur meta-analysis confirms the known association between the HLA-DRB1 shared epitope and RA radiological damage. The lack of replication of previously validated non-HLA markers highlights a requirement for further research to deliver clinically-useful prognostic genetic markers. Show less
Spiliopoulou, A.; Colombo, M.; Plant, D.; Nair, N.; Cui, J.; Coenen, M.J.H.; ... ; McKeigue, P. 2019
Objective. Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors... Show moreObjective. Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy. Methods. A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (Delta DAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models. Results. Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the Delta DAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99). Conclusion. Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections. Show less
To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases... Show moreTo identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 x 10(-8)) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles. Show less