Integrative analysis highlighted susceptibility genes for rheumatoid arthritis (2023)

Section snippets

Overview of the study design

First, we performed SMR analysis to detect putative causal genes for RA by integrating data from RA GWAS with data from each of four eQTL studies. Second, we tested the probability of the identified genes clustering into a specific biological process by functional annotation enrichment analysis. Third, we evaluated the associations between mRNA levels of the identified genes and RA in blood cells and Synovial tissue. Fourth, the connections between novel genes and known RA genes were evaluated

Inference of RA causal genes

The SMR analysis identified a long list of genes as potential causal genes underlying RA GWAS associations (PSMR<5×10–6) (Supplementary Tables 1–3). Among these genes, 119 and 53 (140 unique) were found to be associated with RA in Europeans (Supplementary Table 2) and Asians (Supplementary Table 3), respectively. Of these genes, 52 passed the HEIDI tests (PHEIDI>0.05), suggesting that there was no heterogeneity, and the expression of these genes and RA were affected by same variants.


The current study represented an effort to identify RA-associated genes by integrating GWAS data with eQTL data using the SMR method. We identified a list of potential causal genes with high confidence for RA. The identified genes were highly connected and enriched in specific biological pathways which were related to RA. Moreover, the findings have highlighted a potential functional role of IFNAR2 in RA pathology.

Previous GWASs have confirmed hundreds of genetic variants for RA susceptibility.

CRediT authorship contribution statement

Xing-Bo Mo: Conceptualization, Methodology, Software, Data curation, Writing - original draft, Visualization, Investigation. Yang-Hua Sun: . Yong-Hong Zhang: Writing - review & editing. Shu-Feng Lei: Supervision, Writing - review & editing.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.


The study was supported by Natural Science Foundation of China (81872681, 31401079 and 81473046), the Startup Fund from Soochow University (Q413900313), Natural Science Fund for Colleges and Universities in Jiangsu Province (18KJB330005), and a Project of the Priority Academic Program Development of Jiangsu Higher Education Institutions.

Recommended articles (6)

  • Research article

    Retention of endothelial progenitor cells in bone marrow in a murine model of endogenous tissue plasminogen activator (tPA) deficiency in response to critical limb ischemia

    International Journal of Cardiology, Volume 170, Issue 3, 2014, pp. 394-405

    This study tested the hypothesis that tissue plasminogen activator (tPA) is crucial for regulating endothelial progenitor cell (EPC) mobilization from bone marrow to circulation in murine critical limb ischemia (CLI) by ligating the left femoral artery.

    Wild-type (C57BL/6) (n=40) mice were equally divided into group 1A (sham control), group 2A (CLI), group 3A [control-tPA (4.0mg/kg)] and group 4A [CLI-tPA (intravenously at 3h after CLI)]. Similarly, tPA knock-out (tPA−/−) mice (n=40) were equally divided into group 1B (sham control), group 2B (CLI), group 3B [control-tPA (4.0mg/kg)], and group 4B (CLI-tPA).

    The circulating levels of EPCs (C-kit/CD31+, Sca-1/KDR+, CXCR4/CD34+) were lower in groups 1B and 2B than in groups 1A and 2A, respectively (all p<0.01), and were reversed after tPA treatment (3B vs. 3A or 4B vs. 4A, p>0.05) at 6h and 18h post-CLI. Levels of these biomarkers decreased again 14days after CLI in tPA−/− mice compared to those in wild-type between the respective groups (all p<0.01). Laser Doppler flowmetry showed a higher ratio of ischemic-to-normal blood flow in 2A than in 2B and in 4A than in 4B by day 14 after CLI (all p<0.05). Angiogenesis at protein (CXCR4, SDF-1α, VEGF) and cellular (CXCR4+, SDF-1α+, and CD31+ cells) levels was highest in animals with CLI-tPA, significantly higher in mice with CLI only than in sham controls for both wild-type and tPA−/− mice (p<0.01).

    tPA played an essential role in augmenting circulating EPCs, angiogenesis, and blood flow in the ischemic limb in a murine model.

  • Research article

    Associations with Corneal Hysteresis in a Population Cohort: Results from 96 010 UK Biobank Participants

    Ophthalmology, Volume 126, Issue 11, 2019, pp. 1500-1510

    To describe the distribution of corneal hysteresis (CH) in a large cohort and explore its associated factors and possible clinical applications.

    Cross-sectional study within the UK Biobank, a large cohort study in the United Kingdom.

    We analyzed CH data from 93 345 eligible participants in the UK Biobank cohort, aged 40 to 69 years.

    All analyses were performed using left eye data. Linear regression models were used to evaluate associations between CH and demographic, lifestyle, ocular, and systemic variables. Piecewise logistic regression models were used to explore the relationship between self-reported glaucoma and CH.

    Corneal hysteresis (mmHg).

    The mean CH was 10.6 mmHg (10.4 mmHg in male and 10.8 mmHg in female participants). After adjusting for covariables, CH was significantly negatively associated with male sex, age, black ethnicity, self-reported glaucoma, diastolic blood pressure, and height. Corneal hysteresis was significantly positively associated with smoking, hyperopia, diabetes, systemic lupus erythematosus (SLE), greater deprivation (Townsend index), and Goldmann-correlated intraocular pressure (IOPg). Self-reported glaucoma and CH were significantly associated when CH was less than 10.1 mmHg (odds ratio, 0.86; 95% confidence interval, 0.79–0.94 per mmHg CH increase) after adjusting for covariables. When CH exceeded 10.1 mmHg, there was no significant association between CH and self-reported glaucoma.

    In our analyses, CH was significantly associated with factors including age, sex, and ethnicity, which should be taken into account when interpreting CH values. In our cohort, lower CH was significantly associated with a higher prevalence of self-reported glaucoma when CH was less than 10.1 mmHg. Corneal hysteresis may serve as a biomarker aiding glaucoma case detection.

  • Research article

    The plasminogen receptor, Plg‐RKT, is essential for mammary lobuloalveolar development and lactation

    Journal of Thrombosis and Haemostasis, Volume 16, Issue 5, 2018, pp. 919-932

    Lactational competence requires plasminogen, the zymogen of the serine protease, plasmin. Plg‐RKT is a unique transmembrane plasminogen receptor that promotes plasminogen activation to plasmin on cell surfaces. Plg‐RKT−/− mice are viable, but no offspring of Plg‐RKT−/− female mice survive to weaning.

    We investigated potential lactational failure in Plg‐RKT−/− mice and addressed causal mechanisms.

    Fibrin accumulation, macrophage infiltration, processing of extracellular matrix components, effects of genetic deletion of fibrinogen, expression of fibrosis genes, and proliferation and apoptosis of epithelial cells were examined in lactating mammary glands of Plg‐RKT−/− and Plg‐RKT+/+ mice.

    Milk was not present in the stomachs of offspring of Plg‐RKT−/− female mice and the pups were rescued by foster mothers. Although the mammary ductal tree developed normally in Plg‐RKT−/− glands, lobuloalveolar development was blocked by a hypertrophic fibrotic stroma and infiltrating macrophages were present. A massive accumulation of fibrin was also present in Plg‐RKT−/− alveoli and ducts. Although this accumulation was decreased when Plg‐RKT−/− mice were made genetically heterozygous for fibrinogen, defects in lobuloalveolar development were not rescued by fibrinogen heterozygosity. Transcriptional profiling revealed that EGF was downregulated 12‐fold in Plg‐RKT−/−glands. Furthermore, proliferation of epithelial cells was not detectable. In addition, the pro‐survival protein, Mcl‐1, was markedly downregulated and apoptosis was observed in Plg‐RKT−/− but not Plg‐RKT+/+glands.

    Plg‐RKT is essential for lactogenesis and functions to maintain the appropriate stromal extracellular matrix environment, regulate epithelial cell proliferation and apoptosis, and, by regulating fibrinolysis, preserve alveolar and ductal patency.

  • Research article

    Enclosings of λ-fold 5-cycle systems for u=2

    Discrete Mathematics, Volume 338, Issue 5, 2015, pp. 743-765

    A k-cycle system of a multigraph G is an ordered pair (V,C) where V is the vertex set of G and C is a set of k-cycles, the edges of which partition the edges of G. A k-cycle system of λKv is known as a λ-fold k-cycle system of order v. A k-cycle system of λKv(V,C) is said to be enclosed in a k-cycle system of (λ+m)Kv+u(VU,P) if CP and u,m1. In this paper the enclosing problem for 5-cycle systems is settled in the general situation where the three parameters λ, m, and v are allowed to vary freely and u is constrained to the difficult case of adding two vertices. New graph theoretic approaches are introduced to handle this situation developing an avenue of research that is of interest in its own right.

  • Research article

    Interleukin-2 maintains the survival of interleukin-17+ gamma/delta T cells in inflammation and autoimmune diseases

    International Immunopharmacology, Volume 86, 2020, Article 106721

    There is increasing appreciation of the critical pathogenic role of IL-17 in inflammation and autoimmune diseases, which could be produced from both adaptive Th17 cells and innate γδ T cells. Existing evidences suggest that IL-2 is important for in vivo accumulation of IL-17+ γδ T cells, leaving the mechanisms still elusive. Herein, using lupus-prone MRL/lpr mice, we demonstrated that splenic γδ T cells were potent IL-17 producers at the onset of lupus, which could be diminished by in vivo IL-2 neutralization. Additional in vivo results showed that neutralization of IL-2 also significantly deleted the IL-17-producing γδ T cells in ovalbumin (OVA) /CFA-immunized B6 mice. Using splenic γδ T cells from OVA/CFA-immunized B6 mice, we further demonstrated that IL-2 could induce IL-17 production alone or together with IL-1β or IL-23 or anti-TCRγδ. Mechanism studies demonstrated that IL-2 could support the survival of γδ T cells, rather than induce the proliferation. Through specific pharmacologic inhibitor, we demonstrated that IL-2 could maintain that RORγt expression of γδ T cells in a STAT5-dependent manner. Collectively, this study suggested that the interplay between IL and 2 and other pro-inflammatory cytokines could trigger the rapid IL-17 production from innate γδ T cells, thus to orchestrate an inflammatory response before the development of adaptive Th17 cells.

  • Research article

    Molecular mechanisms of autophagic memory in pathogenic T cells in human arthritis

    Journal of Autoimmunity, Volume 94, 2018, pp. 90-98

    T-cell resilience is critical to the immune pathogenesis of human autoimmune arthritis. Autophagy is essential for memory T cell generation and associated with pathogenesis in rheumatoid arthritis (RA). Our aim here was to delineate the role and molecular mechanism of autophagy in resilience and persistence of pathogenic T cells from autoimmune arthritis.

    We demonstrated “Autophagic memory” as elevated autophagy levels in CD4+ memory T cells compared to CD4+ naive T cells and in Jurkat Human T cell line trained with starvation stress. We then showed increased levels of autophagy in pathogenic CD4+ T cells subsets from autoimmune arthritis patients. Using RNA-sequencing, transcription factor gene regulatory network and methylation analyses we identified MYC as a key regulator of autophagic memory. We validated MYC levels using qPCR and further demonstrated that inhibiting MYC increased autophagy. The present study proposes the novel concept of autophagic memory and suggests that autophagic memory confers metabolic advantage to pathogenic T cells from arthritis and supports its resilience and long term survival. Particularly, suppression of MYC imparted the heightened autophagy levels in pathogenic T cells. These studies have a direct translational valency as they identify autophagy and its metabolic controllers as a novel therapeutic target.

© 2020 Elsevier B.V. All rights reserved.

Top Articles
Latest Posts
Article information

Author: Duane Harber

Last Updated: 21/04/2023

Views: 6389

Rating: 4 / 5 (71 voted)

Reviews: 94% of readers found this page helpful

Author information

Name: Duane Harber

Birthday: 1999-10-17

Address: Apt. 404 9899 Magnolia Roads, Port Royceville, ID 78186

Phone: +186911129794335

Job: Human Hospitality Planner

Hobby: Listening to music, Orienteering, Knapping, Dance, Mountain biking, Fishing, Pottery

Introduction: My name is Duane Harber, I am a modern, clever, handsome, fair, agreeable, inexpensive, beautiful person who loves writing and wants to share my knowledge and understanding with you.