Remission is a primary end point of in clinical practice and trials of treatments for rheumatoid arthritis (RA). The 2011 American College of Rheumatology/European League Against Rheumatism (ACR/EULAR) remission criteria were developed to provide a consensus definition of remission. This study aimed to assess the concordance between the new remission criteria and the physician’s clinical judgment of remission and also to identify factors that affect the discordance between these two approaches. A total of 3,209 patients with RA were included from the KORean Observational Study Network for Arthritis (KORONA) database. The frequency of remission was evaluated based on each approach. The agreement between the results was estimated by Cohen’s kappa (κ). Patients with remission according to the 2011 ACR/EULAR criteria (i.e. the Boolean criteria) and/or physician judgment (n = 855) were divided into three groups: concordant remission, the Boolean criteria only, and physician judgment only. Multinomial logistic regression analysis was used to identify factors responsible for the assignment of patients with remission to one of the discordant groups rather than the concordant group. The remission rates using the Boolean criteria and physician judgment were 10.5% and 19.9%, respectively. The agreement between two approaches for remission was low (κ = 0.226) and the concordant remission rate was only 5.5% (n = 177). Pain affected classification in both discordant groups, whereas fatigue was associated with remission only by physician clinical judgment. The Boolean criteria were more stringent than clinical judgment. Patient subjective symptoms such as pain and fatigue were associated with discordance between the two approaches.
Arthritis, Rheumatoid; Remission; Discordance
This study aimed to refine the interaction between cigarette smoking and human leukocyte antigen (HLA) polymorphisms in seropositive rheumatoid arthritis (RA), in the context of a recent amino-acid based HLA model for RA susceptibility.
We imputed HLA amino acids and classical alleles from case-control Immunochip array data of 3,588 Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA; case/control=1654/1934), 589 Nurses’ Health Study (NHS; 229/360) and 2,125 Korean (1390/735) subjects. We examined interaction effects between heavy smoking (>10 pack-years) and genetic risk score (GRS) from RA-associated amino-acid positions (11, 13, 71 and 74 in HLA-DRβ1; 9 in HLA-B; and 9 in HLA-DPβ1) and from HLA-DRβ1 four-amino-acid haplotypes with an attributable proportion due to interaction (AP) using the additive interaction model.
Heavy smoking and all investigated HLA amino-acid positions and haplotypes were associated with RA susceptibility in all three populations. In the interaction analysis, we found a significant deviation from the expected additive joint effect between heavy smoking and the HLA-DRβ1 amino-acid haplotype in all three studies (0.416≤AP≤0.796). We further identified the key interacting variants as being located at amino-acid positions 11 and 13 of HLA-DRβ1 but not the other RA-risk amino-acid positions in all populations. At the positions 11 and 13, there were similar patterns between RA-risk effects and interaction effects of residues.
Our findings of significant gene-environment interaction effects implicate that a physical interaction between citrullinated auto-antigens produced by smoking and HLA-DR molecules is characterized by the HLA-DRβ1 four-amino-acid haplotype, primarily by the positions 11 and 13.
Objective. Although two recent randomized placebo-controlled trials of rituximab (RTX) failed to demonstrate efficacy in systemic lupus erythematosus (SLE), clinicians continue to use off-label RTX for cases refractory to current treatments. We evaluated the effectiveness and safety of rituximab for patients with refractory SLE in Korea. Methods. We retrospectively analyzed multicenter patients treated with RTX in Korea. Results. 39 SLE patients treated with RTX were included in the following manner: lupus nephritis 43.6%, hematologic 33.3%, arthritis 7.8%, myositis 7.8%, and others 7.7%. All patients had responded poorly to at least one conventional immunosuppressive agent (mean 2.5 ± 1.1, cyclophosphamide 43.6%, mycophenolate mofetil 48.7%, and other drugs) before RTX. Clinical improvements (complete or partial remission) occurred in patients with renal disease, hematologic disease, arthritis, myositis, and other manifestations at 6 months after RTX. The SLEDAI score was significantly decreased from 10.8 ± 7.1 at baseline to 6.7 ± 4.0 at 6 months, 6.2 ± 4.1 at 12 months, and 5.5 ± 3.6 at 24 months after RTX (P < 0.05). Among 28 clinical responders, 4 patients experienced a relapse of disease at 25 ± 4 months. Infections were noted in 3 patients (7.7%).
Conclusion. RTX could be an effective and relatively safe therapeutic option in patients with severe refractory SLE until novel B-cell depletion therapy is available.
Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. We sought to identify the genetic loci associated with SLE in a Korean population by performing an unbiased genome-wide association scan.
A total of 1,174 Korean SLE cases and 4,248 population controls were genotyped with strict quality control measures and analyzed for association. For select variants, replication was tested in an independent set of 1,412 SLE cases and 1,163 population controls of Korean and Chinese ancestries.
Eleven regions outside the HLA exceeded genome-wide significance (P<5×10−8). A novel SNP-SLE association was identified between FCHSD2 and P2RY2 peaking at rs11235667 (P = 1.0×10−8, odds ratio (OR) = 0.59) on a 33kb haplotype upstream to ATG16L2. Replication for rs11235667 resulted in Pmeta-rep=0.001 and Pmeta-overall=6.67×10−11 (OR=0.63). Within the HLA region, association peaked in the Class II region at rs116727542 with multiple independent effects. Classical HLA allele imputation identified HLA-DRB1*1501 and HLA-DQB1*0602, both highly correlated, as most strongly associated with SLE. We replicated ten previously established SLE risk loci: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1 and IRAK1-MECP2. Of these loci, we identified previously unreported independent second effects in TNFAIP3 and TNFSF4 as well as differences in the association for a putative causal variant in the WDFY4 region.
Further studies are needed to identify true SLE risk effects in other suggestive loci and to identify the causal variant(s) in the regions of ATG16L2, FCHSD2, and P2RY2.
Systemic lupus erythematosus (SLE) has a strong but incompletely understood genetic architecture. We conducted an association study with replication in 4,492 SLE cases and 12,675 controls from six East-Asian cohorts, to identify novel and better localize known SLE susceptibility loci. We identified 10 novel loci as well as 20 known loci with genome-wide significance. Among the novel loci, the most significant was GTF2IRD1-GTF2I at 7q11.23 (rs73366469, Pmeta=3.75×10−117, OR=2.38), followed by DEF6, IL12B, TCF7, TERT, CD226, PCNXL3, RASGRP1, SYNGR1 and SIGLEC6. We localized the most likely functional variants for each locus by analyzing epigenetic marks and gene regulation data. Ten putative variants are known to alter cis- or trans-gene expression. Enrichment analysis highlights the importance of these loci in B- and T-cell biology. Together with previously known loci, the explained heritability of SLE increases to 24%. Novel loci share functional and ontological characteristics with previously reported loci, and are possible drug targets for SLE therapeutics.
Considerable sharing of disease alleles among populations is well-characterized in autoimmune disorders (e.g., rheumatoid arthritis), but there are some exceptional loci showing heterogenic association among populations. Here we investigated genetic variants with distinct effects on the development of rheumatoid arthritis in Asian and European populations. Ancestry-related association heterogeneity was examined using Cochran’s homogeneity tests for the disease association data from large Asian (n = 14,465; 9,299 discovery subjects and 5,166 validation subjects; 4 collections) and European (n = 45,790; 11 collections) rheumatoid arthritis case-control cohorts with Immunochip and genome-wide SNP array data. We identified significant heterogeneity between the two ancestries for the common variants in the GTF2I locus (PHeterogeneity = 9.6 × 10−9 at rs73366469) and showed that this heterogeneity was due to an Asian-specific association effect (ORMeta = 1.37 and PMeta = 4.2 × 10−13 in Asians; ORMeta = 1.00 and PMeta = 1.00 in Europeans). Trans-ancestral comparison and bioinfomatics analysis revealed a plausibly causal or disease-variant-tagging SNP (rs117026326; in linkage disequilibrium with rs73366469), whose minor allele is common in Asians but rare in Europeans. In conclusion, we identified largest-ever effect on Asian rheumatoid arthritis across human non-HLA regions at GTF2I by heterogeneity mapping followed by replication studies, and pinpointed a possible causal variant.
Systemic lupus erythematosus (SLE) predominantly affects women in their reproductive years and has a significant impact on childbearing. We investigated the influence of personal decision on family size among Korean women with SLE and factors that affect the decisions. A case-control study comparing childbearing history and decisions of 112 SLE patients and 135 controls was performed. Women with SLE participating in the Network for Lupus Clinical Research in South Korea and matching controls between ages of 18-45, who are/were married or living with a partner were included. Data regarding socio-demographics, reproductive history, and childbearing decisions were collected through a survey using a standardized questionnaire and medical record review. More women with SLE reported at least one pregnancy (85.7% vs. 71.9%, P = 0.009) or at least one live birth (85.7% vs. 71.9%, P = 0.003) compared with controls. Mean number of pregnancies was significantly higher (2.4 ± 1.6 vs. 1.4 ± 1.3, P < 0.001), and mean number of live births was significantly lower in women with SLE (1.2 ± 0.8 vs. 1.6 ± 0.8, P < 0.001). Significantly more women with SLE made the decision not to have children compared with controls (54.5% vs. 40.7%, P = 0.031), and health-related concerns were the major cause of the decision. Other socio-demographic factors did not influence the decision to limit childbearing in SLE women. The disease-related concerns had significant impact on family size and childbearing decisions among Korean women with SLE.
Systemic Lupus Erythematosus; Family Size; Childbearing Decisions; Korean
Patients who develop an active tuberculosis infection during tumor necrosis factor (TNF) inhibitor treatment typically discontinue TNF inhibitor and receive standard anti-tuberculosis treatment. However, there is currently insufficient information on patient outcomes following resumption of TNF inhibitor treatment during ongoing anti- tuberculosis treatment. Our study was designed to investigate the safety of resuming TNF inhibitors in ankylosing spondylitis (AS) patients who developed tuberculosis as a complication of the use of TNF inhibitors.
Through the nationwide registry of the Korean Society of Spondyloarthritis Research, 3929 AS patients who were prescribed TNF inhibitors were recruited between June 2003 and June 2014 at fourteen referral hospitals. Clinical information was analyzed about the patients who experienced tuberculosis after exposure to TNF inhibitors. The clinical features of resumers and non-resumers of TNF inhibitors were compared and the outcomes of tuberculosis were surveyed individually.
Fifty-six AS patients were treated for tuberculosis associated with TNF inhibitors. Among them, 23 patients resumed TNF inhibitors, and these patients were found to be exposed to TNF inhibitors for a longer period of time and experienced more frequent disease flare-up after discontinuation of TNF inhibitors compared with those who did not resume. Fifteen patients resumed TNF inhibitors during anti-tuberculosis treatment (early resumers) and 8 after completion of anti-tuberculosis treatment (late resumers). Median time to resuming TNF inhibitor from tuberculosis was 3.3 and 9.0 months in the early and late resumers, respectively. Tuberculosis was treated successfully in all resumers and did not relapse in any of them during follow-up (median 33.8 [IQR; 20.8–66.7] months).
Instances of tuberculosis were treated successfully in our AS patients, even when given concomitantly with TNF inhibitors. We suggest that early resumption of TNF inhibitors in AS patients could be safe under effective coverage of tuberculosis.
The genetic association of HLA-DRB1 with rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) is well documented, but association with other HLA-DR beta genes (HLA-DRB3, HLA-DRB4 and HLA-DRB5) has not been thoroughly studied, despite their similar functions and chromosomal positions. We examined variants in all functional HLA-DR beta genes in RA and SLE patients and controls, down to the amino-acid level, to better understand disease association with the HLA-DR locus. To this end, we improved an existing HLA reference panel to impute variants in all protein-coding HLA-DR beta genes. Using the reference panel, HLA variants were inferred from high-density SNP data of 9,271 RA-control subjects and 5,342 SLE-control subjects. Disease association tests were performed by logistic regression and log-likelihood ratio tests. After imputation using the newly constructed HLA reference panel and statistical analysis, we observed that HLA-DRB1 variants better accounted for the association between MHC and susceptibility to RA and SLE than did the other three HLA-DRB variants. Moreover, there were no secondary effects in HLA-DRB3, HLA-DRB4, or HLA-DRB5 in RA or SLE. Of all the HLA-DR beta chain paralogs, those encoded by HLA-DRB1 solely or dominantly influence susceptibility to RA and SLE.
The aim of this study was to investigate the prevalence of osteoporosis in rheumatoid arthritis (RA) patients and to analyze the risk factors in these patients using the KORean Observational study Network for Arthritis (KORONA) database.
Among the RA patients in the KORONA who were recruited between July 2009 and December 2011, postmenopausal women with bone mineral density (BMD) results within one year from the time of KORONA enrollment were included in this study. The baseline characteristics of patients in three groups, defined by BMD results, were compared. The BMD measurement rates and prevalence of osteoporosis in the study patients were calculated in accordance with age and gender subgroups. Multivariable logistic regression analysis was used to explore the association between osteoporosis and demographics and disease-related risk factors.
Of 1322 postmenopausal woman patients with RA in whom BMD was measured within one year of study enrollment, 619 patients (46.8 %) were in the osteoporosis group (T-score ≤ −2.5 SD). RA patients with osteoporosis had a higher frequency of previous fractures than those in other groups, especially fractures of the femur (p = 0.004) and wrist (p = 0.042). Advanced age (≥70 years; OR = 2.28, 95 % CI: 1.40–3.58), lower body mass index (<25; OR = 2.14, 95 % CI:1.52–3.02), longer disease duration (≥10 years; OR = 1.46, 95 % CI: 1.07–2.00), higher cumulative glucocorticoid dose (OR = 1.03, 95 % CI: 1.01–1.05), and higher Health Assessment Questionnaire score (OR = 1.37, 95 % CI:1.11–1.69) were independent risk factors for osteoporosis.
A large percentage (90.8 %) of RA patients enrolled in the KORONA cohort had osteoporosis and osteopenia. Nevertheless, BMD measurement rates in this population remained low, despite high risk groups of fractures.
Osteoporosis; Rheumatoid arthritis; Frequency; Risk factors
Systemic lupus erythematosus (SLE; OMIM 152700) is characterised by the production of antibodies to nuclear antigens. We previously identified variants in complement receptor 2 (CR2/CD21) that were associated with decreased risk of SLE. This study aimed to identify the causal variant for this association.
Genotyped and imputed genetic variants spanning CR2 were assessed for association with SLE in 15 750 case-control subjects from four ancestral groups. Allele-specific functional effects of associated variants were determined using quantitative real-time PCR, quantitative flow cytometry, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP)-PCR.
The strongest association signal was detected at rs1876453 in intron 1 of CR2 (pmeta=4.2×10−4, OR 0.85), specifically when subjects were stratified based on the presence of dsDNA autoantibodies (case-control pmeta=7.6×10−7, OR 0.71; case-only pmeta=1.9×10−4, OR 0.75). Although allele-specific effects on B cell CR2 mRNA or protein levels were not identified, levels of complement receptor 1 (CR1/CD35) mRNA and protein were significantly higher on B cells of subjects harbouring the minor allele (p=0.0248 and p=0.0006, respectively). The minor allele altered the formation of several DNA protein complexes by EMSA, including one containing CCCTC-binding factor (CTCF), an effect that was confirmed by ChIP-PCR.
These data suggest that rs1876453 in CR2 has long-range effects on gene regulation that decrease susceptibility to lupus. Since the minor allele at rs1876453 is preferentially associated with reduced risk of the highly specific dsDNA autoantibodies that are present in preclinical, active and severe lupus, understanding its mechanisms will have important therapeutic implications.
Systemic Lupus Erythematosus; Autoantibodies; Gene Polymorphism; B cells
Rheumatoid arthritis (RA) has been known to be associated with increased risk of cardiovascular disease (CVD). The aim of this study was to investigate the effects of Tai Chi exercise on CVD risk in elderly women with RA.
In total, 56 female patients with RA were assigned to either a Tai Chi exercise group (29 patients) receiving a 3-month exercise intervention once a week or a control group (27 patients) receiving general information about the benefits of exercise. All participants were assessed at baseline and at 3 months for RA disease activity (Disease Activity Score 28 and Routine Assessment of Patient Index Data 3), functional disability (Health Assessment Questionnaire), CVD risk factors (blood pressure, lipids profile, body composition, and smoking), and three atherosclerotic measurements: carotid intima-media thickness, flow-mediated dilatation (FMD), and brachial-ankle pulse wave velocity (baPWV).
FMD, representative of endothelial function, significantly increased in the Tai Chi exercise group (initial 5.85 ± 2.05 versus 3 months 7.75 ± 2.53 %) compared with the control group (initial 6.31 ± 2.12 versus 3 months 5.78 ± 2.13 %) (P = 1.76 × 10−3). Moreover, baPWV, representative of arterial stiffness, significantly decreased in the Tai Chi exercise group (initial 1693.7 ± 348.3 versus 3 months 1600.1 ± 291.0 cm/s) compared with the control group (initial 1740.3 ± 185.3 versus 3 months 1792.8 ± 326.1 cm/s) (P = 1.57 × 10−2). In addition, total cholesterol decreased significantly in the Tai Chi exercise group compared with the control group (−7.8 ± 15.5 versus 2.9 ± 12.2 mg/dl, P = 2.72 × 10−2); other changes in RA-related characteristics were not significantly different between the two groups. Tai Chi exercise remained significantly associated with improved endothelial function (FMD; P = 4.32 × 10−3) and arterial stiffness (baPWV; P = 2.22 × 10−2) after adjustment for improvement in total cholesterol level.
Tai Chi exercise improved endothelial dysfunction and arterial stiffness in elderly women with RA, suggesting that it can be a useful behavioral strategy for CVD prevention in patients with RA.
Electronic supplementary material
The online version of this article (doi:10.1186/s13075-015-0893-x) contains supplementary material, which is available to authorized users.
Rheumatoid arthritis; Tai Chi; Cardiovascular risk
Previous studies have emphasized ethnically heterogeneous human leukocyte antigen (HLA) classical allele associations to rheumatoid arthritis (RA) risk. We fine-mapped RA risk alleles within the major histocompatibility complex (MHC) in 2782 seropositive RA cases and 4315 controls of Asian descent. We applied imputation to determine genotypes for eight class I and II HLA genes to Asian populations for the first time using a newly constructed pan-Asian reference panel. First, we empirically measured high imputation accuracy in Asian samples. Then we observed the most significant association in HLA-DRβ1 at amino acid position 13, located outside the classical shared epitope (Pomnibus = 6.9 × 10−135). The individual residues at position 13 have relative effects that are consistent with published effects in European populations (His > Phe > Arg > Tyr ≅ Gly > Ser)—but the observed effects in Asians are generally smaller. Applying stepwise conditional analysis, we identified additional independent associations at positions 57 (conditional Pomnibus = 2.2 × 10−33) and 74 (conditional Pomnibus = 1.1 × 10−8). Outside of HLA-DRβ1, we observed independent effects for amino acid polymorphisms within HLA-B (Asp9, conditional P = 3.8 × 10−6) and HLA-DPβ1 (Phe9, conditional P = 3.0 × 10−5) concordant with European populations. Our trans-ethnic HLA fine-mapping study reveals that (i) a common set of amino acid residues confer shared effects in European and Asian populations and (ii) these same effects can explain ethnically heterogeneous classical allelic associations (e.g. HLA-DRB1*09:01) due to allele frequency differences between populations. Our study illustrates the value of high-resolution imputation for fine-mapping causal variants in the MHC.
The presence of the HLA-B*27 allele is a major risk factor for the development of ankylosing spondylitis (AS), which causes chronic inflammation of the spine and other sites. We investigated residual effects outside HLA-B within the major histocompatibility complex (MHC) region in the Korean population.
Using the Korean HLA reference panel, we inferred the classic HLA alleles and amino-acid residues of the six HLA genes (HLA-A, -B, -C,-DPB1, -DQB1, and -DRB1) and MHC single-nucleotide polymorphisms in 3820 Korean subjects, including 654 Korean cases of AS and 3166 controls, who were genotyped by using Immunochip. Logistic regression and log-likelihood ratio tests were used in AS association tests for imputed markers.
The most significant associations were identified at amino-acid positions in the epitope-binding site of HLA-B (P = 1.71 × 10−481 at position 70, P = 7.20 × 10−479 at position 97, and P = 2.54 × 10−484 at positions 114), highlighting the risk effect of the HLA-B*27 allele and the protective effects of other classic alleles. A secondary effect was located at the leucine at amino-acid position 116 in the epitope-binding site of HLA-C (P = 1.69 × 10−14), completely tagging the HLA-C*15:02 allele. This residue had a large effect in HLA-B*27-negative patients (odds ratio = 6.6, 95 % confidence interval = 3.8 to 11.4).
The four amino-acid positions of HLA-B and -C account for most of the associations between AS and MHC in the Korean population. This finding updates the list of AS susceptibility loci and provides new insight into AS pathogenesis mediated by MHC class I molecules.
Electronic supplementary material
The online version of this article (doi:10.1186/s13075-015-0855-3) contains supplementary material, which is available to authorized users.
We have previously demonstrated that semimature dendritic cell- (smDC-) based immunotherapy is effective for the treatment of collagen-induced arthritis (CIA) prior to disease onset. In the present study, we examined the efficacy of combination therapy with smDCs and methotrexate (MTX) in advanced CIA with a score of 2-3. Combination therapy with low-dose MTX and type II collagen- (CII-) pulsed smDCs (CII-smDCs) was more effective in inhibiting disease progression than high or low-dose MTX alone or a combination of high dose MTX and CII-smDCs. The effect of CII-smDCs alone was also comparable to the combination therapy. CD4+Foxp3+ Treg populations and IL-10 secretion markedly increased, and CII-specific autoreactive T cells decreased in mice treated with CII-smDCs alone or in combination with MTX. Combination therapy reduced the secretion of interferon-γ (IFN-γ) and IL-17 with little influence on the IL-4 secretion in the mixed leukocyte reaction. These results imply that the combination therapy with low-dose MTX and smDCs is effective in controlling advanced CIA by enhancing Treg population and suppresses antigen-specific Th1/Th17 immunity, rather than initiating Th1 to Th2 immune deviation. Our findings provide a better understanding of the DC therapy in combination with MTX for the treatment of patients with rheumatoid arthritis (RA).
A highly polygenic etiology and high degree of allele-sharing between ancestries have been well-elucidated in genetic studies of rheumatoid arthritis. Recently, the high-density genotyping array Immunochip for immune disease loci identified 14 new rheumatoid arthritis risk loci among individuals of European ancestry. Here, we aimed to identify new rheumatoid arthritis risk loci using Korean-specific Immunochip data.
We analyzed Korean rheumatoid arthritis case-control samples using the Immunochip and GWAS array to search for new risk alleles of rheumatoid arthritis with anti-citrullinated peptide antibodies. To increase power, we performed a meta-analysis of Korean data with previously published European Immunochip and GWAS data, for a total sample size of 9,299 Korean and 45,790 European case-control samples.
We identified 8 new rheumatoid arthritis susceptibility loci (TNFSF4, LBH, EOMES, ETS1–FLI1, COG6, RAD51B, UBASH3A and SYNGR1) that passed a genome-wide significance threshold (p<5×10−8), with evidence for three independent risk alleles at 1q25/TNFSF4. The risk alleles from the 7 new loci except for the TNFSF4 locus (monomorphic in Koreans), together with risk alleles from previously established RA risk loci, exhibited a high correlation of effect sizes between ancestries. Further, we refined the number of SNPs that represent potentially causal variants through a trans-ethnic comparison of densely genotyped SNPs.
This study demonstrates the advantage of dense-mapping and trans-ancestral analysis for identification of potentially causal SNPs. In addition, our findings support the importance of T cells in the pathogenesis and the fact of frequent overlap of risk loci among diverse autoimmune diseases.
Rheumatoid arthritis; Gene polymorphism; Anti-CCP
Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (PEA = 1.01 × 10−54, PHS = 3.68 × 10−10, PAA = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10−9), and rs13306575 in HS and KR (PHS = 7.04 × 10−7, PKR = 3.30 × 10−3). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10−7), implying that SLE predisposing variants were tagged. Significant SNP–SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance (‘missing heritability’) of complex diseases like SLE.
Genetic variations of human leukocyte antigen (HLA) genes within the major histocompatibility complex (MHC) locus are strongly associated with disease susceptibility and prognosis for many diseases, including many autoimmune diseases. In this study, we developed a Korean HLA reference panel for imputing classical alleles and amino acid residues of several HLA genes. An HLA reference panel has potential for use in identifying and fine-mapping disease associations with the MHC locus in East Asian populations, including Koreans. A total of 413 unrelated Korean subjects were analyzed for single nucleotide polymorphisms (SNPs) at the MHC locus and six HLA genes, including HLA-A, -B, -C, -DRB1, -DPB1, and -DQB1. The HLA reference panel was constructed by phasing the 5,858 MHC SNPs, 233 classical HLA alleles, and 1,387 amino acid residue markers from 1,025 amino acid positions as binary variables. The imputation accuracy of the HLA reference panel was assessed by measuring concordance rates between imputed and genotyped alleles of the HLA genes from a subset of the study subjects and East Asian HapMap individuals. Average concordance rates were 95.6% and 91.1% at 2-digit and 4-digit allele resolutions, respectively. The imputation accuracy was minimally affected by SNP density of a test dataset for imputation. In conclusion, the Korean HLA reference panel we developed was highly suitable for imputing HLA alleles and amino acids from MHC SNPs in East Asians, including Koreans.
Although it has been suggested that rare coding variants could explain the substantial missing heritability, very few sequencing studies have been performed in rheumatoid arthritis (RA). We aimed to identify novel functional variants with rare to low frequency using targeted exon sequencing of RA in Korea.
We analyzed targeted exon sequencing data of 398 genes selected from a multifaceted approach in Korean RA patients (n = 1,217) and controls (n = 717). We conducted a single-marker association test and a gene-based analysis of rare variants. For meta-analysis or enrichment tests, we also used ethnically matched independent samples of Korean genome-wide association studies (GWAS) (n = 4,799) or immunochip data (n = 4,722).
After stringent quality control, we analyzed 10,588 variants of 398 genes from 1,934 Korean RA case controls. We identified 13 nonsynonymous variants with nominal association in single-variant association tests. In a meta-analysis, we did not find any novel variant with genome-wide significance for RA risk. Using a gene-based approach, we identified 17 genes with nominal burden signals. Among them, VSTM1 showed the greatest association with RA (P = 7.80 × 10−4). In the enrichment test using Korean GWAS, although the significant signal appeared to be driven by total genic variants, we found no evidence for enriched association of coding variants only with RA.
We were unable to identify rare coding variants with large effect to explain the missing heritability for RA in the current targeted resequencing study. Our study raises skepticism about exon sequencing of targeted genes for complex diseases like RA.
Electronic supplementary material
The online version of this article (doi:10.1186/s13075-014-0447-7) contains supplementary material, which is available to authorized users.
A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological datasets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here, we performed a genome-wide association study (GWAS) meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single nucleotide polymorphisms (SNPs). We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 1012–4. We devised an in-silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci (cis-eQTL)6, and pathway analyses7–9 – as well as novel methods based on genetic overlap with human primary immunodeficiency (PID), hematological cancer somatic mutations and knock-out mouse phenotypes – to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.
Systemic lupus erythematosus (SLE; OMIM 152700) is a chronic autoimmune disease for which the aetiology includes genetic and environmental factors. ITGAM, integrin αΜ (complement component 3 receptor 3 subunit) encoding a ligand for intracellular adhesion molecule (ICAM) proteins, is an established SLE susceptibility locus. This study aimed to evaluate the independent and joint effects of genetic variations in the genes that encode ITGAM and ICAM.
The authors examined several markers in the ICAM1–ICAM4–ICAM5 locus on chromosome 19p13 and the single ITGAM polymorphism (rs1143679) using a large-scale case–control study of 17 481 unrelated participants from four ancestry populations. The single marker association and gene–gene interaction were analysed for each ancestry, and a meta-analysis across the four ancestries was performed.
The A-allele of ICAM1–ICAM4–ICAM5 rs3093030, associated with elevated plasma levels of soluble ICAM1, and the A-allele of ITGAM rs1143679 showed the strongest association with increased SLE susceptibility in each of the ancestry populations and the trans-ancestry meta-analysis (ORmeta=1.16, 95% CI 1.11 to 1.22; p=4.88×10−10 and ORmeta=1.67, 95% CI 1.55 to 1.79; p=3.32×10−46, respectively). The effect of the ICAM single-nucleotide polymorphisms (SNPs) was independent of the effect of the ITGAM SNP rs1143679, and carriers of both ICAM rs3093030-AA and ITGAM rs1143679-AA had an OR of 4.08 compared with those with no risk allele in either SNP (95% CI 2.09 to 7.98; p=3.91×10−5).
These findings are the first to suggest that an ICAM–integrin-mediated pathway contributes to susceptibility to SLE.
Recently, interferon gamma releasing assay has been recommended to compensate the tuberculin skin test (TST) for screening for latent tuberculosis infection (LTBI). Although it improved the detection of LTBI before treatment with tumor necrosis factor blocker, its application to immune suppressed patients is limited. We report a case of peritoneal tuberculosis (TB) developed in a patient who tested positive for TST and QuantiFERON-TB Gold (QFT-G) before infliximab therapy, to emphasize the importance of monitoring during treatment. A 52-year-old woman presented with abdominal distension. She had been diagnosed with seropositive rheumatoid arthritis six years ago. She had started taking infliximab six months ago. All screening tests for TB were performed and the results of all were negative. At admission, the results of repeated TST and QFT-G tests were positive. Histopathological examination confirmed peritoneal TB. The patient started anti-TB therapy and the symptoms were relieved.
Peritonitis, Tuberculous; Infliximab
The genetic association of the major histocompatibility complex (MHC) to rheumatoid arthritis risk has commonly been attributed to HLA-DRB1 alleles. Yet controversy persists about the causal variants in HLA-DRB1 and the presence of independent effects elsewhere in the MHC. Using existing genome-wide SNP data in 5,018 seropositive cases and 14,974 controls, we imputed and tested classical alleles and amino acid polymorphisms for HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, and DRB1 along with 3,117 SNPs across the MHC. Conditional and haplotype analyses reveal that three amino acid positions (11, 71 and 74) in HLA-DRβ1, and single amino acid polymorphisms in HLA-B (position 9) and HLA-DPβ1 (position 9), all located in the peptide-binding grooves, almost completely explain the MHC association to disease risk. This study illustrates how imputation of functional variation from large reference panels can help fine-map association signals in the MHC.
Hepatitis B virus (HBV) may be reactivated after chemotherapy or immunosuppressive therapy, and therefore administration of antiviral agents before such treatment is recommended. Most reported cases of reactivation are associated with high doses of immunosuppressive agents or combination therapy. We present a case of a previously inactive HBV carrier with an acute severe flare-up during a long-term, very-low-dose (2.5 mg/day) steroid treatment for rheumatoid arthritis. We suggest that even a minimal dose of single-regimen oral steroid can cause reactivation of indolent, inactive HBV.
Hepatitis B virus; Liver failure; Steroids; Virus activation
Recent studies have identified a number of novel rheumatoid arthritis (RA) loci in Caucasian populations. In this study, we sought to determine whether the genetic variants at 4q27, 6q23, CCL21, TRAF1/C5, and CD40 identified in Caucasians are also associated with RA in a Korean case-control collection. We also comprehensively evaluated the genetic variation within PTPN22, a well established autoimmune disease-associated gene.
We designed a Sequenom iPlex experiment to thoroughly evaluate the PTPN22 linkage disequilibrium region using tag SNPs and disease-associated SNPs at 5 other previously reported Caucasian RA-associated loci in 1123 RA Korean RA patients and 1008 ethnically matched controls. We also re-sequenced the PTPN22 gene to look for novel coding variants that might be contributing to disease in this population.
None of the Caucasian RA susceptibility loci contributed significantly to disease in Koreans. Tag SNPs covering the PTPN22 linkage disequilibrium block, while polymorphic, did not reveal any disease association and re-sequencing did not identify any new common coding region variants in this population. The 6q23 and 4q27 SNPs assayed were non-polymorphic in this population and the TRAF1/C5, CD40, and CCL21 SNPs did not show any evidence for association.
Caucasian and Korean rheumatoid arthritis have different genetic risk factors. While patients of different ethnic groups share the HLA region as a major genetic risk locus, most other genes shown to be significantly associated with disease in Caucasians appear not to play a role in Korean RA.