Dataset Information

Integrative blood-derived epigenetic and transcriptomic analysis reveals the potential regulatory role of DNA methylation in ankylosing spondylitis.

ABSTRACT:

Background

The currently known risk loci could explain a small proportion of the heritability of ankylosing spondylitis (AS). Epigenetics might account for the missing heritability. We aimed to seek more novel AS-associated DNA methylation alterations and delineate the regulatory effect of DNA methylation and gene expression with integrated analysis of methylome and transcriptome.

Methods

Epigenome-wide DNA methylation and mRNA expression were profiled in peripheral blood mononuclear cells (PBMCs) from 45 individuals (AS: health controls (HCs) = 30:15) with high-throughput array. The methylome was validated in an independent cohort (AS: HCs = 12:12). Pearson correlation analysis and causal inference tests (CIT) were conducted to determine potentially causative regulatory effects of methylation on mRNA expression.

Results

A total of 4794 differentially methylated positions (DMPs) were identified associated with AS, 2526 DMPs of which were validated in an independent cohort. Both cohorts highlighted T cell receptor (TCR) signaling and Th17 differentiation pathways. Besides, AS patients manifested increased DNA methylation variability. The methylation levels of 158 DMPs were correlated with the mRNA expression levels of 112 genes, which formed interconnected network concentrated on Th17 cell differentiation and TCR signaling pathway (LCK, FYN, CD3G, TCF7, ZAP70, CXCL12, and PLCG1). We also identified several cis-acting DNA methylation and gene expression changes associated with AS risk, which might regulate the cellular mechanisms underlying AS.

Conclusions

Our studies outlined the landscapes of epi-signatures of AS and several methylation-gene expression-AS regulatory axis and highlighted the Th17 cell differentiation and TCR signaling pathway, which might provide innovative molecular targets for therapeutic interventions for AS.

SUBMITTER: Xiao M

PROVIDER: S-EPMC8728943 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundCircular RNA (circRNA) is a type of closed circular noncoding RNA (ncRNA), mostly formed by back-splicing or alternative splicing of pre-messenger RNA (mRNA). The aim of this study was to explore the expression profile of circRNA in peripheral blood mononuclear cells (PBMCs) of patients with ankylosing spondylitis (AS) and discover potential molecular markers of AS.MethodsThe circRNA microarray technology was used to detect the expression of circRNAs in the peripheral blood of 6 patients with AS and 6 healthy controls (HC). To screen the differentially expressed circRNAs by fold change (FC) and P value, these differentially expressed circRNAs were analyzed by bioinformatics. In 60 cases of AS and 30 cases of HC, 4 circRNAs were subjected to real-time fluorescence quantitative polymerase chain reaction (RT-qPCR), and their correlation with various clinical indicators was analyzed. Finally, the receiver operating characteristic (ROC) curve was used to analyze their potential as AS diagnostic markers.ResultsThe microarray results showed that there were 1369 significantly differently expressed (P < 0.05, FC > 1.5) circRNAs between the AS and HC groups (675 upregulated and 694 downregulated). The results of bioinformatics analysis suggested that they were mainly involved in "enzyme binding," "adenosine ribonucleotide binding," "MAPK signaling pathway", etc. The RT-qPCR results showed that the expressions of hsa_circRNA_001544 (U = 486.5, P < 0.05) and hsa_circRNA_102532 (U = 645, P < 0.05) were significantly different between the AS group and the HC group. The AS group was further divided into two subgroups: active AS (ASA) and stable AS (ASS). After analysis, it was found that compared with the HC group, hsa_circRNA_001544 was significantly increased in both ASA (U = 214, P < 0.05) and ASS groups (U = 273, P < 0.05), while hsa_circRNA_008961 (U = 250, P < 0.05) and hsa_circRNA_102532 (U = 295, P < 0.05) were only significantly increased in the ASA group. Furthermore, hsa_circRNA_012732 was significantly different between the ASA and ASS groups (U = 194, P < 0.05), and there was no statistical significance among the remaining groups. Correlation analysis results showed that hsa_circRNA_012732 was negatively correlated with Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), high-sensitivity C-reactive protein (hsCRP), and globulin (GLOB) and positively correlated with lymphocyte count (LY), mean corpusular volume, and albumin (ALB), and hsa_circRNA_008961 was negatively correlated with platelet (PLT) count. ROC curve analysis showed that hsa_circRNA_001544 (95% CI = 0.610-0.831, P < 0.05) and hsa_circRNA_102532 (95% CI = 0.521-0.762, P < 0.05) were statistically significant, and their area under curve (AUC) values were 0.720 and 0.642, respectively.ConclusionsThere are differentially expressed circRNAs in PBMCs of AS patients, and they may be involved in the occurrence and development of AS. Among these differentially expressed circRNAs, hsa_circRNA_012732 has the potential to become an indicator of disease activity, and hsa_circRNA_001544 has the potential to become a molecular marker for AS diagnosis.

Project description:Introduction: A number of genetic-association studies have identified genes contributing to AS susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a “snapshot” of the sampled cells activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: 18 active AS patients, classified according to the New York criteria. and 18 gender-and age-matched controls were profiled using Illumina HT-12 Whole-Genome Expression BeadChips which carry cDNAs for 48000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan Low Density Arrays (TLDAs). Results: 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a p-value <0.0005 (80% confidence level of false discovery rate). Forty seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 down-regulated 1.3-2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300 which modulate cartilage and bone metabolism. Conclusion: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease. RNA was extracted from whole blood using PAXGene tubes. 16 AS patients with active disease and 16 gender- and age-matched controls were analysed.

Project description:We have compared synovial biopsies from ankylosing spondylitis and undifferentiated spondylitis patients with healthy controls and osteoarthritis patients Objective: In spondylarthropies, whole-genome gene expression profiling studies have been limited to peripheral blood to date. By undertaking a study in knee synovial biopsies from spondylarthropy (SpA) and ankylosing spondylitis (AS) patients we aimed to identified joint-specific candidate genes and pathways. These pathways may mediate systemic inflammation driven joint damaging processes and more specifically, the osteoproliferation that is characteristic of these conditions. Methods: RNA was extracted from six seronegative SpA, two AS, three osteoarthritis (OA) and four normal control knee synovial biopsies. Whole genome expression profiling was undertaken using the Illumina DASL system, which assays 24000 cDNA probes. Differentially expressed candidate genes were then validated using quantitative PCR and immunohistochemistry. Results: 416 differentially expressed genes were identified that clearly delineated between AS/SpA and control groups. Pathway analysis showed altered gene-expression in oxidoreductase activity, osteoblast activity, B-cell associated, matrix catabolic, and metabolic pathways. The inflammatory mediator, MMP3, was strongly upregulated in AS/SpA samples and the Wnt pathway inhibitors DKK3 and Kremen1 were downregulated. Conclusion: Pathways mediating both systemic inflammation as well as local tissue changes were identified. This suggests initial systemic inflammation in spondylarthropies transfers to and persists in the local joint environment, subsequently mediating changes in genes directly involved in the destructive tissue remodelling. Fifteen knee synovial biopsy tissue samples consisting of six seronegative spondyloarthropy (SpA), two ankylosing spondylitis (AS), three osteoarthritis (OA) and four normal control biopsies were obtained from the Synovial Tissue Bank at the Repatriation General Hospital in Adelaide, South Australia with the appropriate ethical approval (Supplementary Table 1). All patients provided informed written consent.