rs657075 (CSF2) Is Associated with the Disease Phenotype (BAS-G) of Ankylosing Spondylitis

Assessing the association of leukocyte telomere length with ankylosing spondylitis and rheumatoid arthritis: A bidirectional Mendelian randomization study

Background

Telomere length shortening can cause senescence and apoptosis in various immune cells, resulting in immune destabilization and ageing of the organism. In this study, we aimed to systematically assess the causal relationship of leukocyte telomere length (LTL) with ankylosing spondylitis (AS) and rheumatoid arthritis (RA) using a Mendelian randomization study.

Methods

LTL (n=472174) was obtained from the UK Biobank genome-wide association study pooled data. AS (n=229640), RA (n=212472) were obtained from FinnGen database. MR-Egger, inverse variance weighting, and weighted median methods were used to estimate the effects of causes. Cochran's Q test, MR Egger intercept test, MR-PRESSO, leave-one-out analysis, and funnel plots were used to look at sensitivity, heterogeneity, and multiple effects. Forward MR analysis considered LTL as the exposure and AS, RA as the outcome. Reverse MR analysis considered AS, RA as the exposure and LTL as the outcome.

Results

In the forward MR analysis, inverse variance-weighted and weighted median analysis results indicated that longer LTL might be associated with increased risk of AS (IVW: OR = 1.55, 95% CI: 1.14-2.11, p = 0.006). MR Egger regression analysis showed no pleiotropy between instrumental variables (IVs) (Egger intercept= 0.008, p = 0.294). The leave-one-out analysis showed that each single nucleotide polymorphism (SNP) of AS was robust to each outcome. No significant causal effects were found between AS, RA and LTL in the reverse MR analysis.

Conclusion

Longer LTL may be related with an increased risk of developing AS, and these findings provide a foundation for future clinical research on the causal association between LTL and AS.

Genetic variation in CSF2 (5q31.1) is associated with longitudinal susceptibility to pediatric malaria, severe malarial anemia, and all-cause mortality in a high-burden malaria and HIV region of Kenya

Plasmodium falciparum infections remain among the leading causes of morbidity and mortality in holoendemic transmission areas. Located within region 5q31.1, the colony-stimulating factor 2 gene (CSF2) encodes granulocyte–macrophage colony-stimulating factor (GM-CSF), a hematopoietic growth factor that mediates host immune responses. Since the effect of CSF2 variation on malaria pathogenesis remains unreported, we investigated the impact of two genetic variants in the 5q31.1 gene region flanking CSF2:g-7032 G > A (rs168681:G > A) and CSF2:g.64544T > C (rs246835:T > C) on the rate and timing of malaria and severe malarial anemia (SMA, Hb < 5.0 g/dL) episodes over 36 months of follow-up. Children (n = 1654, aged 2–70 months) were recruited from a holoendemic P. falciparum transmission area of western Kenya. Decreased incidence rate ratio (IRR) for malaria was conferred by inheritance of the CSF2:g.64544 TC genotype (P = 0.0277) and CSF2 AC/GC diplotype (P = 0.0015). Increased IRR for malaria was observed in carriers of the CSF2 AT/GC diplotype (P = 0.0237), while the inheritance of the CSF2 AT haplotype increased the IRR for SMA (P = 0.0166). A model estimating the longitudinal risk of malaria showed decreased hazard rates among CSF2 AC haplotype carriers (P = 0.0045). Investigation of all-cause mortality revealed that inheritance of the GA genotype at CSF2:g-7032 increased the risk of mortality (P = 0.0315). Higher risk of SMA and all-cause mortality were observed in younger children (P < 0.0001 and P = 0.0015), HIV-1(+) individuals (P < 0.0001 and P < 0.0001), and carriers of HbSS (P = 0.0342 and P = 0.0019). Results from this holoendemic P. falciparum area show that variation in gene region 5q31.1 influences susceptibility to malaria, SMA, and mortality, as does age, HIV-1 status, and inheritance of HbSS.

Insights into the pathogenesis of psoriatic arthritis from genetic studies

Psoriatic arthritis (PsA) is a relatively common inflammatory arthritis, a spondyloarthritis (SpA), that occurs most often in patients with psoriasis, a common immune-mediated inflammatory skin disease. Both psoriasis and PsA are highly heritable. Genetic and recent genomic studies have identified variants associated with psoriasis and PsA, but variants differentiating psoriasis from PsA are few. In this review, we describe recent developments in understanding the genetic burden of PsA, linkage, association and epigenetic studies. Using pathway analysis, we provide further insights into the similarities and differences between PsA and psoriasis, as well as between PsA and other immune-mediated inflammatory diseases, particularly ankylosing spondylitis, another SpA. Environmental factors that may trigger PsA in patients with psoriasis are also reviewed. To further understand the pathogenetic differences between PsA and psoriasis as well as other SpA, larger cohort studies of well-phenotyped subjects with integrated analysis of genomic, epigenomic, transcriptomic, proteomic and metabolomic data using interomic system biology approaches are required.

Genetic influences on susceptibility to rheumatoid arthritis in African-Americans

Large meta-analyses of rheumatoid arthritis (RA) susceptibility in European (EUR) and East Asian (EAS) populations have identified >100 RA risk loci, but genome-wide studies of RA in African-Americans (AAs) are absent. To address this disparity, we performed an analysis of 916 AA RA patients and 1392 controls and aggregated our data with genotyping data from >100 000 EUR and Asian RA patients and controls. We identified two novel risk loci that appear to be specific to AAs: GPC5 and RBFOX1 (PAA < 5 × 10-9). Most RA risk loci are shared across different ethnicities, but among discordant loci, we observed strong enrichment of variants having large effect sizes. We found strong evidence of effect concordance for only 3 of the 21 largest effect index variants in EURs. We used the trans-ethnic fine-mapping algorithm PAINTOR3 to prioritize risk variants in >90 RA risk loci. Addition of AA data to those of EUR and EAS descent enabled identification of seven novel high-confidence candidate pathogenic variants (defined by posterior probability > 0.8). In summary, our trans-ethnic analyses are the first to include AAs, identified several new RA risk loci and point to candidate pathogenic variants that may underlie this common autoimmune disease. These findings may lead to better ways to diagnose or stratify treatment approaches in RA.

Identification of the key genes and long non‑coding RNAs in ankylosing spondylitis using RNA sequencing

Ankylosing spondylitis (AS) is an insidious and debilitating form of arthritis that involves the axial skeleton, and its etiology and pathogenesis remain unclear. In the present study, three patients with AS and three normal controls from our hospital were enrolled. RNA sequencing and bioinformatics analysis were performed in order to identify the differentially expressed (DE) mRNAs (DEmRNAs) and DE long non‑coding RNAs (DElncRNAs) between the patients with AS and normal controls. Construction of an AS‑specific protein‑protein interaction network, a weighted DElncRNA‑DEmRNA co‑expression network and functional annotation of the DEmRNAs co‑expressed with DElncRNAs was performed. Nearby cis‑targeted DEmRNAs or DElncRNAs were identified by searching for DEmRNAs that were transcribed within 100‑kb up‑ or downstream of DElncRNAs. Based on the Gene Expression Omnibus datasets GSE25101 and GSE73754, the expression of selected DEmRNAs and DElncRNAs were verified using published RNA sequencing data from blood samples, and receiver operating characteristic analysis of selected DEmRNAs was performed. Compared with the normal controls, 1,072 DEmRNAs and 372 DElncRNAs in the patients with AS were identified. Caspase recruitment domain family member 11 and DNA methyltransferase 1 have great diagnostic value for AS. MSTRG.8559 and LINC00987 were also identified as two hub DElncRNAs. The T‑cell receptor signaling pathway was a significantly enriched pathway of the DEmRNAs co‑expressed with DElncRNAs in patients with AS. In conclusion, the present study identified the key DEmRNAs and DElncRNAs in AS, which provides novel information for understanding the pathogenesis of AS and developing potential biomarkers for AS.

rs2841277 (PLD4) is associated with susceptibility and rs4672495 is associated with disease activity in rheumatoid arthritis

Rheumatoid arthritis (RA) is one of the most common autoimmune diseases, can lead to long-term joint damage, chronic pain, and loss of motor function in the hands, and may share some common genetic factors with other autoimmune disorders, such as ankylosing spondylitis (AS). Many single-nucleotide polymorphisms (SNPs) were reported by genome-wide association studies (GWASs) of RA, but some of them have not been examined in the Taiwanese population. In this study, for 15 SNPs reported in previous RA and AS GWASs, we investigated their association with RA in a Taiwanese population. Based on 334 RA patients recruited from the Taichung Veterans General Hospital and 16,036 healthy subjects from the Taiwan Biobank (TWB) project, we observed that subjects having minor allele C at rs2841277 (phospholipase D family, member 4 (PLD4)) have lower susceptibility of RA, compare to those having genotype TT (Odds ratio (OR) = 0.6, p = 3.0 × 10⁻⁶). Among the RA patients, we observed that subjects having GG at rs4672495 have a lower proportion of severe RA, compare to other subjects (OR = 0.09, p = 5.6 × 10⁻³). Results of a bioinformatics approach showed that rs2841277 is able to influence expression of LINC00638 and AHNAK2 and rs4672495 is able to influence the expression of B3GNT2. In summary, this study replicated an association of rs2841277 with RA susceptibility and showed an AS-associated SNP, rs4672495, is associated with RA activity in the Taiwanese population.