Early human migration determines the risk of being attacked by wolves: ethnic gene diversity on the development of systemic lupus erythematosus.

The prevalence of systemic lupus erythematosus (SLE) varies significantly based on ethnicity rather than geographic distribution; thus, the prevalence is higher in Asian, Hispanic, and Black African populations than in European populations. The risk of developing lupus nephritis (LN) is the highest among Asian populations. Therefore, we hypothesize that human genetic diversity between races has occurred through the early human migration and human genetic adaptation to various environments, with a particular focus on pathogens. Additionally, we compile the currently available evidence on the ethnic gene diversity of SLE and how it relates to disease severity. The human leukocyte antigen (HLA) locus is well established as associated with susceptibility to SLE; specific allele distributions have been observed across diverse populations. Notably, specific amino acid residues within these HLA loci demonstrate significant associations with SLE risk. The non-HLA genetic loci associated with SLE risk also varies across diverse ancestries, implicating distinct immunological pathways, such as the type-I interferon and janus kinase-signal transducers and activators of transcription (JAK-STAT) pathways in Asians, the type-II interferon signaling pathway in Europeans, and B cell activation pathway in Africans. Furthermore, assessing individual genetic susceptibility using genetic risk scores (GRS) for SLE helps to reveal the diverse prevalence, age of onset, and clinical phenotypes across different ethnicities. A higher GRS increases the risk of LN and the severity of SLE. Therefore, understanding ethnic gene diversity is crucial for elucidating disease mechanisms and SLE severity, which could enable the development of novel drugs specific to each race.

Correction: Associations of upper respiratory mucosa microbiota with Rheumatoid arthritis, autoantibodies, and disease activity.

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

Associations of upper respiratory mucosa microbiota with Rheumatoid arthritis, autoantibodies, and disease activity.

The lung is recognized as a site for initiating the formation of self-antigen and autoimmune responses in rheumatoid arthritis (RA). We aimed to investigate the association of upper respiratory microbiota with RA, autoantibody production, and disease activity. Forty-six patients with RA and 17 controls were examined. Nasopharyngeal swab samples were sequenced for microbiome profiling using the V3-V4 region of the 16S rRNA gene. The microbial diversity and relative abundance were compared between RA patients and controls. Correlation analyses were conducted to evaluate the relationship between microbial abundance and clinical markers such as autoantibodies and disease activity. Microbial diversity analysis revealed no major differences between RA patients and healthy controls. However, beta diversity analysis indicated a subtle distinction in microbial composition (unweighted UniFrac distance) between the two groups (P = 0.03), hinting at a minor subset of microbiota associated with disease status. Differential abundance analysis uncovered specific taxa at various taxonomic levels, including Saccharibacteria (TM7) [O-1] (PFDR = 2.53 × 10-2), TM7 [F-1] (PFDR = 5.20 × 10-3), Microbacterium (PFDR = 3.37 × 10-4), and Stenotrophomonas (PFDR = 2.57 × 10-3). The relative abundance of ten genera correlated significantly with anti-cyclic citrullinated peptide (anti-CCP) antibody levels (PFDR < 0.05) and 11 genera were significantly associated with disease activity markers, including ESR, CRP, DAS28-ESR, and DAS-CRP (PFDR < 0.05). In particular, Saccharibacteria TM7 [G-3] and Peptostreptococcaceae [XI] [G-1] were correlated with all disease activity biomarkers. Dysbiosis in the upper respiratory mucosa is associated with RA, anti-CCP antibody levels, and disease activity.

Unveiling the dynamics of B lymphocytes in systemic lupus erythematosus patients treated with belimumab through longitudinal single-cell RNA sequencing.

OBJECTIVES: Unraveling the mechanisms underlying treatment response for targeted therapeutics in systemic lupus erythematosus (SLE) patients is challenging due to the limited understanding of diverse responses of circulating immune cells, particularly B cells. We investigated B lymphocyte dynamics during anti-BAFF treatment, utilizing longitudinal single-cell transcriptome data. METHODS: We conducted single-cell RNA sequencing on PBMCs in four Korean SLE patients before and after belimumab treatment at the following time points: 2 weeks, 1, 3, 6, and 12 months. RESULTS: Analyzing over 73 000 PBMCs, we identified 8 distinct subsets of B cells and plasmablasts and analyzed dynamic changes within these cell subsets: initial declines in naive and transitional B cells followed by an increase at three months, contrasted by an initial increase and subsequent decrease in memory B cells by the third month. Meanwhile, plasmablasts exhibited a consistent decline throughout treatment. B cell activation pathways, specifically in naive and memory B cells, were downregulated during the third and sixth months. These findings were validated at the protein level throughout the first four weeks of treatment using flow cytometry. Comparative analysis with bulk transcriptome data from 22 Japanese SLE patients showed increased NR4A1 expression six months post-belimumab treatment, indicating its role in restricting self-reactive B cells, thereby contributing to the biological responses of anti-BAFF treatment. CONCLUSION: The observed B cell dynamics provided insights into the immunological mechanisms underlying the therapeutic effects of anti-BAFF in SLE patients. Furthermore, it underscores the need for research in predicting drug responses based on immune profiling.

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.

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.

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.

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.