ABSTRACT
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.
Subject(s)
Lupus Erythematosus, Systemic , Transcriptome , Humans , Lupus Erythematosus, Systemic/genetics , Sequence Analysis, RNAABSTRACT
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.
Subject(s)
Femur Head Necrosis , Lupus Erythematosus, Systemic , Steroids , Carboxypeptidases/genetics , Carrier Proteins/genetics , Femur Head , Femur Head Necrosis/chemically induced , Femur Head Necrosis/complications , Femur Head Necrosis/genetics , Genetic Predisposition to Disease , Genome-Wide Association Study , Humans , Lupus Erythematosus, Systemic/genetics , MicroRNAs/genetics , Myosin Heavy Chains/genetics , Polymorphism, Single Nucleotide , Steroids/adverse effectsABSTRACT
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.
ABSTRACT
S6K1 has been implicated in a number of key metabolic responses, which contribute to obesity. Critical among these is the control of a transcriptional program required for the commitment of mesenchymal stem cells to the adipocytic lineage. However, in contrast to its role in the cytosol, the functions and targets of nuclear S6K1 are unknown. Here, we show that adipogenic stimuli trigger nuclear translocation of S6K1, leading to H2BS36 phosphorylation and recruitment of EZH2 to H3, which mediates H3K27 trimethylation. This blocks Wnt gene expression, inducing the upregulation of PPARγ and Cebpa and driving increased adipogenesis. Consistent with this finding, white adipose tissue from S6K1-deficient mice exhibits no detectable H2BS36 phosphorylation or H3K27 trimethylation, whereas both responses are highly elevated in obese humans or in mice fed a high-fat diet. These findings define an S6K1-dependent mechanism in early adipogenesis, contributing to the promotion of obesity.
Subject(s)
Adipocytes/enzymology , Adipogenesis , Adipose Tissue/enzymology , Enhancer of Zeste Homolog 2 Protein/metabolism , Histones/metabolism , Obesity/enzymology , Protein Processing, Post-Translational , Ribosomal Protein S6 Kinases, 70-kDa/metabolism , Ribosomal Protein S6 Kinases, 90-kDa/metabolism , Adipose Tissue/pathology , Adiposity , Animals , CCAAT-Enhancer-Binding Proteins/metabolism , Disease Models, Animal , Enhancer of Zeste Homolog 2 Protein/genetics , Epigenesis, Genetic , HeLa Cells , Histones/genetics , Humans , Male , Methylation , Mice , Obesity/genetics , Obesity/pathology , PPAR gamma/genetics , PPAR gamma/metabolism , Phosphorylation , Ribosomal Protein S6 Kinases, 70-kDa/genetics , Ribosomal Protein S6 Kinases, 90-kDa/genetics , Transcription, Genetic , Transfection , Wnt Proteins/genetics , Wnt Proteins/metabolism , Wnt Signaling PathwayABSTRACT
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.
Subject(s)
Spondylitis, Ankylosing , Humans , Spondylitis, Ankylosing/diagnostic imaging , Spondylitis, Ankylosing/genetics , Spondylitis, Ankylosing/drug therapy , Phosphatidylinositol 3-Kinases , Ryanodine Receptor Calcium Release Channel , Radiography , Spine/pathology , Disease Progression , Severity of Illness IndexABSTRACT
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.
ABSTRACT
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.
Subject(s)
Arthritis, Rheumatoid , Biological Products , Inflammatory Bowel Diseases , Infant, Newborn , Humans , Female , Pregnancy , Pregnancy Outcome , Biological Factors , Fetal Growth Retardation , Cesarean Section , Inflammatory Bowel Diseases/drug therapy , Biological Products/adverse effectsABSTRACT
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.
ABSTRACT
Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10-75) and its proxy position 11 (P = 1.1x10-67), followed by HLA-DRB1-37 (P = 4.5x10-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10-8), HLA-DPB1-35 (P = 9.0x10-16), HLA-DQB1-37 (P = 2.7x10-14), and HLA-B-9 (P = 6.5x10-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10-14; DRB1-13, P = 2.9x10-13; DRB1-30, P = 3.9x10-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.
Subject(s)
Amino Acid Sequence , Autoantibodies/immunology , Disease Susceptibility , Histocompatibility Antigens Class II/chemistry , Histocompatibility Antigens Class II/immunology , Histocompatibility Antigens Class I/chemistry , Histocompatibility Antigens Class I/immunology , Lupus Erythematosus, Systemic/etiology , Alleles , Amino Acid Substitution , Asian People , Female , Genetic Predisposition to Disease , Genetic Variation , Histocompatibility Antigens Class I/genetics , Histocompatibility Antigens Class II/genetics , Humans , Male , Odds Ratio , Polymorphism, Single NucleotideABSTRACT
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.
Subject(s)
Arthritis, Rheumatoid/genetics , Arthritis, Rheumatoid/immunology , CD4-Positive T-Lymphocytes/immunology , DNA Methylation/genetics , DNA Methylation/immunology , Adult , Aged , Female , Genetic Variation , Genome-Wide Association Study , Humans , Male , Middle Aged , TranscriptomeABSTRACT
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.
Subject(s)
Back Pain/diagnosis , Chronic Pain/diagnosis , Multifactorial Inheritance , Sacroiliac Joint/diagnostic imaging , Spondylitis, Ankylosing/diagnosis , Adult , Asian People , Back Pain/genetics , Back Pain/metabolism , C-Reactive Protein/metabolism , Case-Control Studies , Chronic Pain/genetics , Chronic Pain/metabolism , Female , HLA-B27 Antigen/genetics , Humans , Magnetic Resonance Imaging , Male , Reproducibility of Results , Risk Factors , Spondylitis, Ankylosing/genetics , Spondylitis, Ankylosing/metabolism , White PeopleABSTRACT
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.
Subject(s)
Asian People/genetics , Genetic Loci/genetics , Genetic Predisposition to Disease/ethnology , Lupus Erythematosus, Systemic/ethnology , Lupus Erythematosus, Systemic/genetics , Adult , Bayes Theorem , Case-Control Studies , China/epidemiology , China/ethnology , Asia, Eastern/ethnology , Female , Genetic Predisposition to Disease/epidemiology , Genetic Variation , Genome-Wide Association Study , Genotype , Humans , Japan/epidemiology , Japan/ethnology , Lupus Erythematosus, Systemic/epidemiology , Male , Middle Aged , Prevalence , Republic of Korea/epidemiology , Republic of Korea/ethnologyABSTRACT
Crohn's disease (CD) and ulcerative colitis (UC) are the major types of chronic inflammatory bowel disease (IBD) characterized by recurring episodes of inflammation of the gastrointestinal tract. Although it is well established that human leukocyte antigen (HLA) is a major risk factor for IBD, it is yet to be determined which HLA alleles or amino acids drive the risks of CD and UC in Asians. To define the roles of HLA for IBD in Asians, we fine-mapped HLA in 12 568 individuals from Korea and Japan (3294 patients with CD, 1522 patients with UC and 7752 controls). We identified that the amino acid position 37 of HLA-DRß1 plays a key role in the susceptibility to CD (presence of serine being protective, P = 3.6 × 10-67, OR = 0.48 [0.45-0.52]). For UC, we confirmed the known association of the haplotype spanning HLA-C*12:02, HLA-B*52:01 and HLA-DRB1*1502 (P = 1.2 × 10-28, OR = 4.01 [3.14-5.12]).
Subject(s)
Colitis, Ulcerative/genetics , Crohn Disease/genetics , Genetic Predisposition to Disease , HLA-DRB1 Chains/genetics , Inflammatory Bowel Diseases/genetics , Alleles , Amino Acid Substitution/genetics , Amino Acids/chemistry , Amino Acids/genetics , Asian People/genetics , Colitis, Ulcerative/pathology , Crohn Disease/pathology , Female , Genetic Association Studies , Genotype , HLA-DRB1 Chains/chemistry , Haplotypes/genetics , Humans , Inflammatory Bowel Diseases/pathology , Japan , Male , Protein Conformation , Republic of KoreaABSTRACT
Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.
Subject(s)
Alleles , Lupus Erythematosus, Systemic/genetics , Polymorphism, Genetic , Quantitative Trait Loci , STAT1 Transcription Factor/genetics , STAT4 Transcription Factor/genetics , Female , Humans , Lupus Erythematosus, Systemic/epidemiology , Male , Risk FactorsABSTRACT
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.
Subject(s)
Arthritis, Rheumatoid/genetics , Genetic Predisposition to Disease/genetics , Case-Control Studies , Genome-Wide Association Study , Humans , Polymorphism, Single Nucleotide , Republic of KoreaABSTRACT
We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.
Subject(s)
Genetic Predisposition to Disease , Genome-Wide Association Study , Lupus Erythematosus, Systemic/genetics , Adaptor Proteins, Signal Transducing/genetics , Asian People , Autophagy-Related Proteins/genetics , Chemokine CCL22/genetics , DNA-Binding Proteins , Gene Expression Regulation/genetics , Genotype , Humans , Kruppel-Like Transcription Factors/genetics , Lupus Erythematosus, Systemic/epidemiology , Lupus Erythematosus, Systemic/pathology , Polymorphism, Single Nucleotide/genetics , Ribonuclease H/genetics , Transcription FactorsABSTRACT
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.
Subject(s)
Arthritis, Rheumatoid/complications , Arthroscopy/methods , Knee Joint/surgery , Osteoarthritis/complications , Popliteal Cyst/surgery , Aged , Female , Humans , Middle Aged , Popliteal Cyst/complications , Treatment OutcomeABSTRACT
Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.
Subject(s)
Lupus Erythematosus, Systemic/genetics , Promoter Regions, Genetic , Transcription, Genetic , src-Family Kinases/genetics , Alleles , Chromosomes, Human, Pair 8 , Electrophoretic Mobility Shift Assay , Female , Genetic Predisposition to Disease , Haplotypes , Humans , Male , Polymorphism, Single NucleotideABSTRACT
To compare the characteristics of rheumatoid arthritis (RA) patients receiving either biosimilar or originator infliximab and to identify the effectiveness and safety of biosimilar infliximab in RA patients in real-world practice. RA patients who started either biosimilar or originator infliximab were selected using the prospective biologic disease-modifying anti-rheumatic drugs (DMARDs) registry: BIOlogics Pharmacoepidemiologic StudY (BIOPSY). Baseline characteristics of the two groups were compared, and short-term treatment outcomes, including DAS28-ESR and HAQ-DI scores, were compared after initiation of biosimilar or originator infliximab. The drug retention rates of the two groups were also compared. A total of 100 RA patients, 55 biosimilar, and 45 originator infliximab users were included in this analysis. Baseline characteristics of age, disease duration, and previous or current medications were similar in the two groups. Baseline DAS28-ESR was higher in the originator infliximab group (6.3 ± 1.1 vs. 5.8 ± 1.1, p = 0.02). The early DAS28-ESR remission rates observed 7.9 ± 1.8 months after starting biosimilar and originator infliximab were 15.0 and 25.0%, respectively (p = 0.47). The change in HAQ-DI did not differ between the two groups (0.4 ± 0.7 vs. 0.4 ± 0.8, p = 0.94). Patients treated with biosimilar infliximab in clinical practice had lower disease activity at the start of treatment than those receiving originator infliximab. Biosimilar infliximab was well-tolerated, safe, and of similar clinical effectiveness to originator infliximab. Larger number of patient and longer follow-up data will be needed to confirm the effectiveness and safety of biosimilar infliximab in clinical practice.
Subject(s)
Antirheumatic Agents/therapeutic use , Arthritis, Rheumatoid/drug therapy , Biological Products/therapeutic use , Biosimilar Pharmaceuticals/therapeutic use , Infliximab/therapeutic use , Tumor Necrosis Factor-alpha/antagonists & inhibitors , Adult , Antirheumatic Agents/adverse effects , Arthritis, Rheumatoid/diagnosis , Arthritis, Rheumatoid/immunology , Biological Products/adverse effects , Biosimilar Pharmaceuticals/adverse effects , Female , Humans , Infliximab/adverse effects , Male , Middle Aged , Patient Reported Outcome Measures , Prospective Studies , Registries , Remission Induction , Republic of Korea , Time Factors , Treatment Outcome , Tumor Necrosis Factor-alpha/immunologyABSTRACT
OBJECTIVE: To compare the clinical effectiveness of two treatment strategies for active rheumatoid arthritis (RA) refractory to conventional synthetic disease-modifying anti-rheumatic drugs (csDMARDs): starting TNF inhibitors (TNFIs) or changing csDMARDs. METHODS: We used two nationwide Korean RA registries for patient selection. TNFI users were selected from the BIOPSY, which is an inception cohort of RA patients starting biologic DMARDs. As a control group, we selected RA patients with moderate or high disease activity from the KORONA database whose treatment was changed to other csDMARDs. After comparing baseline characteristics between the two groups in either unmatched or propensity score matched cohorts, we compared potential differences in the 1-year remission rate as a primary outcome and changes in HAQ-DI and EQ-5D scores as secondary outcomes. RESULTS: A total of 356 TNFI starters and 586 csDMARD changers were identified from each registry as unmatched cohorts, and 294 patients were included in the propensity score matched cohort. In the intention-to-treat analysis, TNFI starters had higher 1-year remission rates than csDMARD changers in both unmatched (19.1 vs. 18.4%, p < 0.01) and matched cohorts (19.7 vs. 15.0%, p < 0.01). In per protocol analysis, TNFI starters had much higher remission rates in unmatched (37.2 vs. 28.0%, p = 0.04) and matched cohorts (35.4 vs. 19.1%, p = 0.04). However, in matched cohorts, no significant differences were observed between two groups in HAQ-DI and EQ-5D scores. CONCLUSIONS: We compared the clinical effectiveness of the two treatment strategies for active RA refractory to csDMARDs. TNFI starters showed higher 1-year remission rates than csDMARD changers.