Association of HLA-DRB1 locus with treatment response to abatacept or TNF inhibitors in patients with seropositive rheumatoid arthritis.

The strongest genetic risk factor for rheumatoid arthritis (RA) has been known as HLA-DRB1 based on amino acid positions 11, 71, and 74. This study analyzed the association between specific HLA-DRB1 locus and treatment response to abatacept or TNF inhibitors (TNFi) in patients with seropositive RA. A total of 374 Korean RA patients were treated with abatacept (n = 110) or TNFi (n = 264). Associations between HLA-DRB1 and treatment response after 6 months were analyzed using multivariable logistic regression. Seropositive RA patients with HLA-DRB1 shared epitope (SE) had a favorable response to abatacept (OR = 3.67, P = 0.067) and an inversely associated response to TNFi (OR 0.57, P = 0.058) based on EULAR response criteria, but the difference was not statistically significant in comparison to those without SE. In analyses using amino acid positions of HLA-DRB1, a significant association was found between valine at amino acid position 11 of SE and good response to abatacept (OR = 6.46, P = 5.4 × 10-3). The VRA haplotype also showed a good response to abatacept (OR = 4.56, P = 0.013), but not to TNFi. Our results suggest that treatment response to abatacept or TNFi may differ depending on HLA-DRB1 locus in seropositive RA, providing valuable insights for selecting optimal therapy.

Transcriptomic network analysis reveals key drivers of response to anti-TNF biologics in patients with rheumatoid arthritis.

OBJECTIVE: Anti-TNF biologics have been widely used to ameliorate disease activity in patients with rheumatoid arthritis (RA). However, a large fraction of patients show a poor response to these agents. Moreover, no clinically applicable predictive biomarkers have been established. This study aimed to identify response-associated biomarkers using longitudinal transcriptomic data in two independent RA cohorts. METHODS: RNA sequencing data from peripheral blood cell samples of Korean and Caucasian RA cohorts before and after initial treatment with anti-TNF biologics were analyzed to assess treatment-induced expression changes that differed between highly reliable excellent and null responders. Weighted correlation network, immune cell composition, and key driver analyses were performed to understand response-associated transcriptomic networks and cell types and their correlation with disease activity indices. RESULTS: In total, 305 response-associated genes showed significantly different treatment-induced expression changes between excellent and null responders. Co-expression network construction and subsequent key driver analysis revealed that 41 response-associated genes played a crucial role as key drivers of transcriptomic alteration in four response-associated networks involved in various immune pathways: type I interferon signalling, myeloid leucocyte activation, B cell activation, and NK cell/lymphocyte-mediated cytotoxicity. Transcriptomic response scores that we developed to estimate the individual-level degree of expression changes in the response-associated key driver genes were significantly correlated with the changes in clinical indices in independent patients with moderate or ambiguous response outcomes. CONCLUSIONS: This study provides response-specific treatment-induced transcriptomic signatures by comparing the transcriptomic landscape between patients with excellent and null responses to anti-TNF drugs at both gene and network levels.

Pregnancy Outcomes Associated With Biologic Agent Exposure in Patients With Several Rheumatic Diseases and Inflammatory Bowel Diseases.

BACKGROUND: This study aimed to analyze pregnancy outcomes based on biologic agents use in women using the nationwide population-based database. METHODS: The study used the claims database to identify women of childbearing age with several rheumatic (rheumatoid arthritis, juvenile idiopathic arthritis, ankylosing spondylitis, psoriatic arthritis) and inflammatory bowel diseases (Crohn's disease and ulcerative colitis) who had pregnancy-related codes between January 2010 and December 2019. We analyzed live births and adverse pregnancy outcomes based on the previous use of biologics. We also stratified the patients according to duration of biologic agent exposure before pregnancy and the use of biologics during pregnancy to analyze the pregnancy outcomes by subgroups. RESULTS: We identified 4,787 patients with pregnancy events. Among them, 1,034 (21.6%) used biologics before pregnancy. Live birth rate was not different between the biologics group and biologics naïve group (75.0% vs. 75.2%). Multivariate analyses showed that biologics use was associated with higher risk of intrauterine growth retardation (odds ratio [OR], 1.780) and lower risk of gestational diabetes mellitus (OR, 0.776) compared with biologics naïve. Biologics use during pregnancy was associated with higher risk of preterm delivery (OR, 1.859), preeclampsia/eclampsia (OR, 1.762), intrauterine growth retardation (OR, 3.487), and cesarean section (OR, 1.831), but lower risk of fetal loss (OR, 0.274) compared with biologics naïve. CONCLUSIONS: Although there was no difference in live birth rate between the biologics group and biologics naïve group, biologics use seems to be associated with several adverse pregnancy outcomes, especially in patients with biologics during pregnancy. Therefore, patients with biologics during pregnancy need to be carefully observed for adverse pregnancy outcomes.

Higher Genetic Risk Loads Confer More Diverse Manifestations and Higher Risk of Lupus Nephritis in Systemic Lupus Erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a highly heritable complex disorder with heterogeneous clinical manifestations. In this study, we aimed to identify the genetic risk load using clinical and serological manifestations in SLE patients. METHODS: We genotyped a total of 1,655 Korean patients with SLE (n = 1,243 as a discovery set and n = 412 as a replication set) using a customized genome-wide single-nucleotide polymorphism (SNP) array, KoreanChip. A weighted genetic risk score (wGRS) for an individual was calculated from 112 well-validated non-HLA SNPs and HLA haplotypes of SLE-risk loci. We analyzed associations between individual wGRS and clinical SLE subphenotypes and autoantibodies using multivariable linear or logistic regression adjusted by onset age, sex, and disease duration. RESULTS: Childhood-onset SLE (<16 years) conferred the highest genetic risk compared with adult-onset (16-50 years) or late-onset (>50 years) SLE (P = 6.8 × 10-6 ). High wGRS significantly increased associations with SLE manifestations, regardless of onset age, sex, and disease duration. Individual wGRS significantly correlated positively with more clinical American College of Rheumatology criteria (ß = 0.143, P = 1.8 × 10-6 ). Subphenotype analysis revealed significant associations between the highest and lowest wGRS quartile with risk of renal disorder (hazard ratio [HR] 1.74, P = 2.2 × 10-8 ) and anti-Sm antibody production (HR 1.85, P = 2.8 × 10-5 ). Higher wGRS markedly modulated the pathogenesis of proliferative and membranous lupus nephritis class III or IV (HR 1.98, P = 1.6 × 10-5 ) and class V (HR 2.79, P = 1.0 × 10-3 ), but especially lupus nephritis class V in anti-Sm-positive SLE (area under the curve 0.68, P = 1.8 × 10-4 ). CONCLUSION: Patients with SLE and high wGRS tended to have earlier age of SLE onset, higher anti-Sm antibody positivity, and more diverse clinical phenotypes. Genetic profiling may predict high risk for lupus nephritis and a diverse clinical course in SLE patients.

Clinical and genetic factors associated with radiographic damage in patients with ankylosing spondylitis.

OBJECTIVES: To identify clinical and genetic factors associated with severe radiographic damage in patients with ankylosing spondylitis (AS). METHODS: We newly generated genome-wide single nucleotide polymorphism data (833K) for 444 patients with AS. The severity of radiographic damage was assessed using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). To identify clinical and genetic factors associated with severe radiographic damage, multiple linear regression analyses were performed. Human AS-osteoprogenitor and control-osteoprogenitor cells were used for functional validation. RESULTS: The significant clinical factors of final mSASSS were baseline mSASSS (ß=0.796, p=3.22×10-75), peripheral joint arthritis (ß=-0.246, p=6.85×10-6), uveitis (ß=0.157, p=1.95×10-3), and smoking (ß=0.130, p=2.72×10-2) after adjusting for sex, age and disease duration. After adjusting significant clinical factors, the Ryanodine receptor 3 (RYR3) gene was associated with severe radiographic damage (p=1.00×10-6). For pathway analysis, the PI3K-Akt signalling pathway was associated with severe radiographic damage in AS (p=2.21×10-4, false discovery rate=0.040). Treatment with rhodamine B, a ligand of RYR3, dose-dependently induced matrix mineralisation of AS osteoprogenitors. However, the rhodamine B-induced accelerated matrix mineralisation was not definitive in control osteoprogenitors. Knockdown of RYR3 inhibited matrix mineralisation in SaOS2 cell lines. CONCLUSIONS: This study identified clinical and genetic factors that contributed to better understanding of the pathogenesis and biology associated with radiographic damage in AS.

MUC5B promoter variant rs35705950, rare but significant susceptibility locus in rheumatoid arthritis-interstitial lung disease with usual interstitial pneumonia in Asian populations.

BACKGROUND: MUC5B variant rs35705950 is the common and most significant risk variant for rheumatoid arthritis-interstitial lung disease (RA-ILD) in Western populations. However, little is known about its significant association with RA-ILD in Asian populations. We here investigate the association of rs35705950 with Korean patients with RA-ILD. METHODS: In this cross-sectional study, we genotyped rs35705950 in 2444 patients with RA. Among them, 683 patients with RA who have chest CT were divided into RA-ILD and RA-noILD. RA-ILD was classified as usual interstitial pneumonia (UIP) and other than UIP. The associations of rs35705950 with RA-ILD and its subtype were analysed using multivariable regression adjusted for age at RA diagnosis. Meta-analysis of a previously reported Japanese dataset and Korean dataset obtained for this study was conducted. RESULTS: The minor allele (T) frequency of rs35705950 was 0.37%, 1.43% and 2.38% in 2444 patients with RA, 105 patients with RA-ILD and 63 patients with UIP, respectively. Genotypic association of rs35705950 with RA-ILD was insignificant (OR 2.49, 95% CI 0.64 to 9.69, p=0.187), but showed significant association with UIP (OR 4.90, 95% CI 1.23 to 19.59, p=0.024) compared with RA-noILD. In meta-analysis (123 UIP and 878 RA-noILD) combining our data with previously reported Japanese data, this variant was found to be significantly associated with UIP (OR 3.51, 95% CI 1.19 to 10.37, p=0.023). CONCLUSION: MUC5B variant rs35705950 is a rare but significant risk factor for Asian patients with RA-ILD with UIP, suggesting a sharing of the genetic background between Asian and Western populations.

Genome-wide association study of the response of patients with diabetic macular edema to intravitreal Anti-VEGF injection.

Diabetic macular edema (DME), a complication of diabetes mellitus, is a leading cause of adult-onset blindness worldwide. Recently, intravitreal anti-VEGF injection has been used as a first-line treatment. This study analyzed the association between the genetic profile of patients with DME and their response to treatment. Intravitreal anti-VEGF injections were administered monthly for three months to Korean patients diagnosed with DME, who were classified into two groups depending on whether they responded to anti-VEGF therapy or showed recurrence within six months. Peripheral blood samples were used for genetic analyses. Genome-wide association analysis results sowed that the genes DIRC3 on chromosome 2 (rs16857280, p = 1.2 × 10-6), SLCO3A1 on chromosome 15 (rs12899055, p = 2.5 × 10-6), and RAB2A on chromosome 8 (rs2272620, p = 4.6 × 10-6) were associated with treatment response to intravitreal anti-VEGF injection. SLC35F1, TMEM132D, KIAA0368, HPCAL1, IGF2BP3, SPN2S, COL23A1, and CREB5 were also related to treatment response (p < 5.0 × 10-5). Using the KEGG pathway analysis, RAB2A and CREB5 were found to be associated with AMPK signaling related to VEGF (p = 0.018). The identified genetic biomarkers can elucidate the factors affecting patient response to intravitreal anti-VEGF injection and help select appropriate therapeutic strategy.

Distinct transcriptome architectures underlying lupus establishment and exacerbation.

Systemic lupus erythematosus (SLE) is a complex autoimmune disease involving multiple immune cells. To elucidate SLE pathogenesis, it is essential to understand the dysregulated gene expression pattern linked to various clinical statuses with a high cellular resolution. Here, we conducted a large-scale transcriptome study with 6,386 RNA sequencing data covering 27 immune cell types from 136 SLE and 89 healthy donors. We profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively. We then identified candidate biological processes unique to each signature. This study suggested the clinical value of disease-activity signatures, which were associated with organ involvement and therapeutic responses. However, disease-activity signatures were less enriched around SLE risk variants than disease-state signatures, suggesting that current genetic studies may not well capture clinically vital biology. Together, we identified comprehensive gene signatures of SLE, which will provide essential foundations for future genomic and genetic studies.

Biological insights into systemic lupus erythematosus through an immune cell-specific transcriptome-wide association study.

OBJECTIVE: Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis. METHODS: We built gene expression predictive models in blood B cells, CD4+ and CD8+ T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches. RESULTS: TWAS identified 171 genes for SLE (p<1.0×10-5); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10-8). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10-9) around CD83. For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1, and that presence of the SLE risk allele decreased ACAP1 expression. CONCLUSIONS: Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations.

Predictive Factors for Renal Response in Lupus Nephritis: A Single-center Prospective Cohort Study.

Objective: To identify the predictive factors for renal response in patients with lupus nephritis (LN). Methods: Patients and data were extracted from a prospective systemic lupus erythematosus cohort in Korea, in which clinical data were collected at 0, 3, 6, and 12 months after induction therapy. Treatment response of LN were evaluated as a complete response (CR), partial response (PR), or non-response (NR) at 3, 6, and 12 months, respectively. Predictive factors for CR at 6 months were evaluated using multivariable Poisson regression analysis. Results: A total of 75 patients with LN who underwent biopsy was enrolled. The mean age at diagnosis of LN was 28.9±9.7 years, and 68 (90.7%) were female. The frequencies of classes III, IV, III+V, IV+V, and V were 20.0%, 44.0%, 16.0%, 12.0%, and 8.0%, respectively. Compared to relapsed LN, new-onset LN showed a lower percentage of glomerulosclerosis (45.5% vs. 76.2%, p=0.013). The overall proportions of CR, PR, and NR at 6 and 12 months were 52.0%, 26.7%, 21.3% and 50.7%, 24.0%, 25.3%, respectively. In multivariate analysis, age at enrollment (odds ratio [OR]=1.02, p=0.022), relapsed LN (OR=0.71, p=0.037), anti-Ro antibody (OR=0.67, p=0.014), and class III LN (OR=1.48, p=0.001) were associated with CR at 6 months. Conclusion: In our prospective cohort, class III LN was a good predictive factor for CR at 6 months in patients with LN, whereas younger age, relapsed LN, and anti-Ro antibody were poor predictive factors.

Novel susceptibility loci for steroid-associated osteonecrosis of the femoral head in systemic lupus erythematosus.

Osteonecrosis of the femoral head (ONFH) involves necrosis of bone and bone marrow of the femoral head caused by ischemia with unknown etiology. Previous genetic studies on ONFH failed to produce consistent results, presumably because ONFH has various causes with different genetic backgrounds and the underlying diseases confounded the associations. Steroid-associated ONFH (S-ONFH) accounts for one-half of all ONFH, and systemic lupus erythematosus (SLE) is a representative disease underlying S-ONFH. We performed a genome-wide association study (GWAS) to identify genetic risk factors for S-ONFH in patients with SLE. We conducted a two-staged GWAS on 636 SLE patients with S-ONFH and 95 588 non-SLE controls. Among the novel loci identified, we determined S-ONFH-specific loci by comparing allele frequencies between SLE patients without S-ONFH and non-SLE controls. We also used Korean datasets comprising 148 S-ONFH cases and 37 015 controls to assess overall significance. We evaluated the functional annotations of significant variants by in silico analyses. The Japanese GWAS identified 4 significant loci together with 12 known SLE susceptibility loci. The four significant variants showed comparable effect sizes on S-ONFH compared with SLE controls and non-SLE controls. Three of the four loci, MIR4293/MIR1265 [odds ratio (OR) = 1.99, P-value = 1.1 × 10-9)], TRIM49/NAALAD2 (OR = 1.65, P-value = 4.8 × 10-8) and MYO16 (OR = 3.91, P-value = 4.9 × 10-10), showed significant associations in the meta-analysis with Korean datasets. Bioinformatics analyses identified MIR4293, NAALAD2 and MYO16 as candidate causal genes. MIR4293 regulates a PPARG-related adipogenesis pathway relevant to S-ONFH. We identified three novel susceptibility loci for S-ONFH in SLE.

Polygenic Risk Scores have high diagnostic capacity in ankylosing spondylitis.

OBJECTIVE: We sought to test the hypothesis that Polygenic Risk Scores (PRSs) have strong capacity to discriminate cases of ankylosing spondylitis (AS) from healthy controls and individuals in the community with chronic back pain. METHODS: PRSs were developed and validated in individuals of European and East Asian ethnicity, using data from genome-wide association studies in 15 585 AS cases and 20 452 controls. The discriminatory values of PRSs in these populations were compared with other widely used diagnostic tests, including C-reactive protein (CRP), HLA-B27 and sacroiliac MRI. RESULTS: In people of European descent, PRS had high discriminatory capacity with area under the curve (AUC) in receiver operator characteristic analysis of 0.924. This was significantly better than for HLA-B27 testing alone (AUC=0.869), MRI (AUC=0.885) or C-reactive protein (AUC=0.700). PRS developed and validated in individuals of East Asian descent performed similarly (AUC=0.948). Assuming a prior probability of AS of 10% such as in patients with chronic back pain under 45 years of age, compared with HLA-B27 testing alone, PRS provides higher positive values for 35% of patients and negative predictive values for 67.5% of patients. For PRS, in people of European descent, the maximum positive predictive value was 78.2% and negative predictive value was 100%, whereas for HLA-B27, these values were 51.9% and 97.9%, respectively. CONCLUSIONS: PRS have higher discriminatory capacity for AS than CRP, sacroiliac MRI or HLA-B27 status alone. For optimal performance, PRS should be developed for use in the specific ethnic groups to which they are to be applied.

Genetic variants shape rheumatoid arthritis-specific transcriptomic features in CD4+ T cells through differential DNA methylation, explaining a substantial proportion of heritability.

OBJECTIVE: CD4+ T cells have been suggested as the most disease-relevant cell type in rheumatoid arthritis (RA) in which RA-risk non-coding variants exhibit allele-specific effects on regulation of RA-driving genes. This study aimed to understand RA-specific signatures in CD4+ T cells using multi-omics data, interpreting inter-omics relationships in shaping the RA transcriptomic landscape. METHODS: We profiled genome-wide variants, gene expression and DNA methylation in CD4+ T cells from 82 patients with RA and 40 healthy controls using high-throughput technologies. We investigated differentially expressed genes (DEGs) and differential methylated regions (DMRs) in RA and localised quantitative trait loci (QTLs) for expression and methylation. We then integrated these based on individual-level correlations to inspect DEG-regulating sources and investigated the potential regulatory roles of RA-risk variants by a partitioned-heritability enrichment analysis with RA genome-wide association summary statistics. RESULTS: A large number of RA-specific DEGs were identified (n=2575), highlighting T cell differentiation and activation pathways. RA-specific DMRs, preferentially located in T cell regulatory regions, were correlated with the expression levels of 548 DEGs mostly in the same topologically associating domains. In addition, expressional variances in 771 and 83 DEGs were partially explained by expression QTLs for DEGs and methylation QTLs (meQTLs) for DEG-correlated DMRs, respectively. A large number of RA variants were moderately to strongly correlated with meQTLs. DEG-correlated DMRs, enriched with meQTLs, had strongly enriched heritability of RA. CONCLUSION: Our findings revealed that the methylomic changes, driven by RA heritability-explaining variants, shape the differential expression of a substantial fraction of DEGs in CD4+ T cells in patients with RA, reinforcing the importance of a multidimensional approach in disease-relevant tissues.

Clinical and Genetic Risk Factors Associated With the Presence of Lupus Nephritis.

Objective: To elucidate whether clinical features and the weighted genetic risk score (wGRS) were associated with the presence of lupus nephritis (LN). Methods: We retrospectively divided patients with systemic lupus erythematosus (SLE, n=1,078) into biopsy-proven LN (n=507) and non-LN groups (non-LN, n=571) Baseline clinical features, serologic markers, and the wGRS were collected The wGRS was calculated from 112 non-human leukocyte antigen (non-HLA) loci and HLA-DRß1 amino acid haplotypes for SLE Associations among clinical features, wGRS, and the presence of LN were identified. Results: In the multivariate analysis, patients with LN were younger at diagnosis (odds ratio [OR]=0.97, p<0.001), had more pleuritis (OR=2.44, p<0.001) and pericarditis (OR=1.62, p=0.029), had a higher detection rate of anti-double stranded deoxyribonucleic acid (anti-dsDNA antibodies, OR=2.22, p<0.001), anti-Smith antibodies (anti-Sm antibodies, OR=1.70, p=0.002), low level of complement (OR=1.37, p=0.043) and absence of antiphospholipid antibodies (aPL antibodies, OR=1.60, p=0.002), and had higher wGRS (OR=1.16, p=0.012) Mediation analysis suggested that anti-Sm antibodies and low complement could be mediators in the relationship between high wGRS and the presence of LN. Conclusion: Onset age, pleuritis, pericarditis, several serologic markers, and wGRS were associated with the presence of LN Anti-Sm antibodies and low complement appeared to mediate the indirect relationship between wGRS and the presence of LN.

Allele-Specific Quantification of HLA-DRB1 Transcripts Reveals Imbalanced Allelic Expression That Modifies the Amino Acid Effects in HLA-DRß1.

OBJECTIVE: HLA association fine-mapping studies have shown the effects of missense variants in HLA-DRB1 on rheumatoid arthritis (RA) susceptibility, prognosis, and autoantibody production. However, the phenotypic effects of expression changes in HLA-DRB1 remain poorly understood. Therefore, we investigated the allele-specific expression of HLA-DRB1 and its effect on an HLA-DRß1 structure-associated trait in RA. METHODS: We quantified the allele-specific expression of each HLA-DRB1 3-field classic allele in 48 Korean RA patients with anti-citrullinated protein antibodies (ACPAs) and 319 healthy European subjects by using both RNA sequencing and HLA-DRB1 genotype data to calculate the relative expression strength of multiple HLA-DRB1 alleles (n = 14 in Koreans and n = 25 in Europeans) in each population. The known association between ACPA level and alanine at position 74 of HLA-DRß1 in ACPA-positive RA was revisited to understand the phenotypic effect of allele-specific expression of HLA-DRB1 by modeling multivariate logistic regression with the genomic dosage or relative expression dosage of Ala-74 in 2 independent sets of 1,723 Korean RA patients with ACPA. RESULTS: The relative expression strength was highly allele-specific, causing imbalanced allelic expression in HLA-DRB1 heterozygotes. The association between HLA-DRß1 Ala-74 and ACPA level in RA was better explained by relative expression dosage of Ala-74 than by the genomic dosage (change in Akaike's information criterion = -6.98). Moreover, the expression variance of Ala-74 in Ala-74 heterozygotes with no genomic variance of Ala-74 was significantly associated with ACPA level (P = 2.26 × 10-3 ). CONCLUSION: Our findings illustrate the advantage of integrating quantitative and qualitative changes in HLA-DRB1 into a single model for understanding HLA-DRB1 associations.

Meta-analysis of 208370 East Asians identifies 113 susceptibility loci for systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations. METHODS: We newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations. RESULTS: We identified 113 genetic regions including 46 novel loci at genome-wide significance (p<5×10-8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (rg=-0.242) and non-albumin protein (rg=0.238). CONCLUSION: This study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

Genome-wide association study in a Korean population identifies six novel susceptibility loci for rheumatoid arthritis.

OBJECTIVE: Genome-wide association studies (GWAS) in rheumatoid arthritis (RA) have discovered over 100 RA loci, explaining patient-relevant RA pathogenesis but showing a large fraction of missing heritability. As a continuous effort, we conducted GWAS in a large Korean RA case-control population. METHODS: We newly generated genome-wide variant data in two independent Korean cohorts comprising 4068 RA cases and 36 487 controls, followed by a whole-genome imputation and a meta-analysis of the disease association results in the two cohorts. By integrating publicly available omics data with the GWAS results, a series of bioinformatic analyses were conducted to prioritise the RA-risk genes in RA loci and to dissect biological mechanisms underlying disease associations. RESULTS: We identified six new RA-risk loci (SLAMF6, CXCL13, SWAP70, NFKBIA, ZFP36L1 and LINC00158) with pmeta<5×10-8 and consistent disease effect sizes in the two cohorts. A total of 122 genes were prioritised from the 6 novel and 13 replicated RA loci based on physical distance, regulatory variants and chromatin interaction. Bioinformatics analyses highlighted potentially RA-relevant tissues (including immune tissues, lung and small intestine) with tissue-specific expression of RA-associated genes and suggested the immune-related gene sets (such as CD40 pathway, IL-21-mediated pathway and citrullination) and the risk-allele sharing with other diseases. CONCLUSION: This study identified six new RA-associated loci that contributed to better understanding of the genetic aetiology and biology in RA.

MHC associations of ankylosing spondylitis in East Asians are complex and involve non-HLA-B27 HLA contributions.

BACKGROUND: The association of HLA-B*27 with AS is amongst the strongest of any known association of a common variant with any human disease. Nonetheless, there is strong evidence indicating that other HLA-B alleles are involved in the disease. European ethnicity studies have demonstrated risk associations with HLA-B*40 and multiple other HLA-B, HLA-A, and HLA class II alleles, and demonstrated that in that ethnic group, the amino acid sequence at position 97 in HLA-B is the key determinant of HLA associations with AS. A recent study in Korean AS cases and controls additionally identified association at HLA-C*15:02. In the current study, we examined the MHC associations of AS in an expanded East Asian cohort. METHODS: A total of 1637 Chinese, Taiwanese, and Korean AS cases meeting the modified New York Criteria for AS, and 1589 ethnically matched controls, were genotyped with the Illumina Immunochip, including a dense coverage of the MHC region. HLA genotypes and amino acid composition were imputed using the SNP2HLA programme using the Han-MHC reference panel based on the data of Han Chinese subjects (n = 9689), and association tested using logistic regression controlling for population stratification effects. RESULTS: A strong association was seen with HLA-B*27 (odds ratio (OR) = 205.3, P = 5.76 × 10-244). Controlling for this association, the strongest risk association is seen with HLA-C*15 at genome-wide significant level (OR = 7.62, P = 9.30 × 10-19), and confirmed association is also seen with HLA-B*40 at suggestive level (OR = 1.65, P = 2.54 × 10-4). At amino acid level, the strongest association seen in uncontrolled analysis was with histidine at position 114 in HLA-B (P = 7.24 × 10-241), but conditional analyses suggest that the primary amino acid associations are with lysine at position 70 and asparagine at position 97. Restriction of the ERAP1 association with HLA-B27-positive AS, previously reported in European subjects, was confirmed in East Asians. CONCLUSIONS: This study confirms in East Asians that the HLA associations of AS are multiple, including previously reported associations at HLA-B*27, HLA-B*40, and HLA-C*15, as well as novel association with HLA-DQB1*04. The HLA-B associations are driven by the amino acids at positions 70 and 97, in the B pocket of HLA-B.

Deletion at 2q14.3 is associated with worse response to TNF-α blockers in patients with rheumatoid arthritis.

BACKGROUND: Structural variations such as copy number variations (CNVs) have a functional impact on various human traits. This study profiled genome-wide CNVs in Korean patients with rheumatoid arthritis (RA) to investigate the efficacy of treatment with TNF-α blockers. METHODS: A total of 357 Korean patients with RA were examined for the efficacy of TNF-α blocker treatment. Disease activity indexes were measured at baseline and 6 months after the treatment. The patients were classified as responders and non-responders based on the change in disease activity indexes according to the EULAR response criteria. CNVs in the same patients were profiled using fluorescence signal intensity data generated by a genome-wide SNP array. The association of CNVs with response to TNF-α blockers was analyzed by multivariate logistic regression accounting for genetic background and clinical factors including body mass index, gender, baseline disease activity, TNF-α blocker used, and methotrexate treatment. RESULTS: The study subjects varied in their responses to TNF-α blockers and had 286 common CNVs in autosomes. We identified that the 3.8-kb deletion at 2q14.3 in 5% of the subjects was associated with response to TNF-α blockers (1.37 × 10- 5 ≤ P ≤ 4.07 × 10- 4) at a false discovery rate threshold of 5%. The deletion in the identified CNV was significantly more frequent in the non-responders than in the responders, indicating worse response to TNF-α blockers in the deletion carriers. The 3.8-kb deletion at 2q14.3 is located in an intergenic region with the binding sites of two transcription factors, MAFF and MAFK. CONCLUSIONS: This study obtained the CNV landscape of Korean patients with RA and identified the common regional deletion associated with poor response to treatment with TNF-α blockers.

Successful arthroscopic treatment of refractory and complicated popliteal cyst associated with rheumatoid arthritis in combination with osteoarthritis: case series and literature review.

Although popliteal cysts are most frequently identified in patients with osteoarthritis (OA), they may occur in patients with rheumatoid arthritis (RA), in which serious complicated cases such as cyst rupture can be developed. The objective of this study was to report four patients with RA (six knees) in combination with OA with a brief review of literature of previous similar published cases. This is a retrospective review of case records of patients with refractory and/or complicated popliteal cysts, who have successfully treated with arthroscopic intervention. We suggest that arthroscopic interventions such as radical debridement, synovectomy, biomechanical valve excision, and/or cystectomy should be considered in patients with refractory and complicated popliteal cysts associated with RA or RA in combination with OA.