Abstract
Background. Oral chronic graft-versus-host disease (cGVHD) impacts quality of life of patients after allogeneic hematopoietic stem cell transplantation (allo-HSCT). However, its precise pathogenesis remains unknown, with potential associations with differential microRNA (miRNA) expression and the TGF-β/Smad signaling pathway.
Objectives. This study aims to explore miRNA expression profiles in the peripheral blood of oral cGVHD patients, focusing on miRNA-769-5p and its relationship with Smad2.
Materials and methods. Peripheral venous blood samples were collected for RNA extraction from 8 patients with oral cGVHD, 8 patients without cGVHD and 8 participants from the healthy control group. The miRNA library was constructed using the Illumina Hiseq 2500 platform. We focused on identifying miRNAs associated with the TGF-β/Smad signaling pathway and subsequently conducted validation experiments. The oral cGVHD and without cGVHD groups were each expanded to include 15 individuals. Peripheral blood samples were subjected to polymerase chain reaction (PCR) analysis to assess miRNA levels and to evaluate Smad2 mRNA levels in peripheral blood mononuclear cells (PBMC). Additionally, enzyme-linked immunosorbent assay (ELISA) was conducted to determine the Smad2 protein levels in peripheral blood.
Results. The most significantly differentially expressed miRNAs among the 3 groups were miRNA-505-5p and miRNA-769-5p. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated an enrichment of the target genes of miRNA-769-5p in the TGF-β signaling pathway. It was observed that miRNA-769-5p expression was higher in patients without oral cGVHD in comparison to those with oral cGVHD. Receiver operating characteristic (ROC) analysis demonstrated that miRNA-769-5p holds diagnostic value for oral cGVHD. As a target of miRNA-769-5p, Smad2 mRNA exhibited a negative correlation with it. Moreover, both Smad2 mRNA and protein levels were higher in patients with oral cGVHD as opposed to those without cGVHD.
Conclusions. Differential expression of miRNAs, particularly the downregulation of miRNA-769-5p, may influence the development of oral cGVHD by diminishing its inhibitory effect on the TGF-β/Smad signaling pathway through its interaction with Smad2.
Key words: oral mucosa, allogeneic hematopoietic stem cell transplantation, miRNAs, chronic graft-versus-host-disease
Background
Chronic graft-versus-host disease (cGVHD) is the leading cause of death following allogeneic hematopoietic stem cell transplantation (allo-HSCT).1 As an alloimmune and autoimmune disease, it is characterized by lichenoid changes and fibrosis affecting various tissues, thereby compromising organ function.2 Its clinical signs manifest first in the mucosal tissues.3 The oral cavity is the primary organ affected by cGVHD, with oral manifestations occurring in 45–83% of cGVHD patients.4 Encompassing lichen planus-like changes, erythema, ulceration, mucoceles, etc., the clinical presentation of oral cGVHD is diverse and can affect any sites within the oral cavity, including the oral mucosa, salivary glands and periodontium.5, 6 Lichen planus-like changes, marked by hyperkeratotic white lines and lacy-appearing lesions on the oral mucosa, serve as the diagnostic feature of oral cGVHD among its oral manifestations.7 Our previous study showed that these lichen planus-like changes persist even when the systemic symptoms of cGVHD in patients are well controlled.8 The repercussions of oral cGVHD significantly impact patient’s quality of life, particularly in terms of eating and nutrient absorption.9 Furthermore, oral cGVHD poses a risk for the development of oral squamous cell carcinomas (OSCCs).10, 11 Given the adverse effects of oral cGVHD on both oral and systemic health, raising awareness is crucial for understanding its occurrence and progression.
The pathogenesis of cGVHD, especially in the oral context, is intricate and remains controversial. Recent studies have shed light on the role of microRNAs (miRNAs) in cGVHD pathogenesis.12 Reikvam et al. were the first to report miRNA serum profiles in cGVHD,13 and miRNA profiles from plasma extracellular vesicles (EVs) have been considered markers of cGVHD onset.14 However, owing to the distinct microenvironments of various organs and tissues, the miRNAs associated with the disease also exhibit variations.15, 16 Consequently, it prompts the question of whether miRNA profiles also change in oral cGVHD patients and how these changes impact the disease.
In our previous research, we uncovered an imbalance of cytokines associated with the TGF-β/Smad pathway in oral cGVHD patients.17 TGF-β/Smad represents a major subfamily within the TGF-β signaling pathway. Phosphorylation of drosophila mothers against decapentaplegic protein 2 (Smad2) serves to activate this pathway, subsequently exerting its effects.18 Research supports the involvement of the TGF-β/Smad signaling pathway in the initiation and progression of cGVHD,19 and it has been shown that therapeutic intervention targeting this pathway can effectively manage cGVHD.20 Furthermore, it is worth noting that Smads can be subject to regulation by miRNAs.21, 22 This leads us to question whether patients with oral cGVHD may exhibit aberrant miRNA expression profiles that influence the TGF-β/Smad pathway and contribute to the development of oral cGVHD.
Objectives
This study aimed to investigate miRNA expression profiles in the peripheral blood of patients affected by oral cGVHD, with a specific focus on miRNA-769-5p and its relationship with Smad2. We aim to provide an initial exploration of the correlation between miRNA-769-5p and Smad2 and their potential influences on oral cGVHD.
Materials and methods
Study groups and approval
The oral cGVHD group comprised patients whose clinical manifestations of cGVHD exclusively affected the oral cavity. The group without cGVHD comprised patients who underwent transplantation more than a year prior and did not develop cGVHD. The healthy control group comprised healthy volunteers.
All study participants were recruited from Guangxi Medical University (Nanning, China). Informed consent, consistent with the principles of the Declaration of Helsinki, was obtained from all study participants. The research protocol received approval from the Research Ethics Committee of Guangxi Medical University (approval No. IRB: 20170301-4). All procedures adhered to the relevant guidelines and regulations.
Study design
Eight individuals were recruited for sampling in each group for RNA sequencing analysis. To further validate the RNA sequencing results in both the oral cGVHD and without cGVHD groups, an additional 7 individuals were included in each group, resulting in a total of 15 individuals in each group for polymerase chain reaction (PCR)-based testing and analysis.
Clinical examination
The oral cGVHD diagnosis and grading were determined based on established diagnostic and scoring criteria.7 Oral discomfort was assessed using a Visual Analogue Scale (VAS). Demographic information was collected through a structured questionnaire, and clinical evaluations and examinations were recorded.
RNA sequencing
Peripheral venous blood was collected from 8 participants in each of the 3 groups. Total RNA was isolated from the whole blood using TRIZOL (Invitrogen, Waltham, USA) to create the miRNA library. Sequencing and the subsequent analysis of miRNA differential expression were performed on the Illumina Hiseq 2500 platform by LC-Bio (Hangzhou, China).
qPCR
Peripheral venous blood samples were collected, and total RNA was extracted using TRIzol (Invitrogen) to measure miRNA levels. Concurrently, peripheral blood mononuclear cells (PBMCs) were isolated and subjected to RNA extraction to determine the Smad2 mRNA levels. The RNA was subsequently reverse-transcribed into cDNA utilizing TruScript™ First Strand cDNA Synthesis Kit. Quantitative real-time polymerase chain reaction (qPCR) was conducted using qTOWER 2.2 system (Analytik Jena AG, Jena, Germany). Cel-miR-39-5p served as the control for miRNA detection, while GAPDH acted as the control for Smads mRNA. Primer sequences are provided in Table 1. Data analysis was performed using the following formula: R = 2−(ΔCt sample − ΔCt control).
ELISA
Plasma samples from participants in the oral cGVHD and without cGVHD groups were collected and subjected to analysis using a human Smad2 immunoblot assay enzyme-linked immunosorbent assay (ELISA) kit (Cusabio, Wuhan, China). Spectrophotometric measurements were taken at 450 nm using the Wellscan Mk3 instrument (Thermo Fisher Scientific, Waltham, USA).
Luciferase reporter assay
The 3’-UTR of Smad2, containing either wild-type (WT) or mutated (MUT) binding sites for miRNA-769-5p, was amplified and inserted into the pGL3 vector to generate the plasmids pGL3-WT-Smad2-3’-UTR or pGL3-MUT-Smad2-3’-UTR, respectively. In the luciferase reporter assay, HEK-293 cells were co-transfected with the luciferase reporter vectors and miRNA-769-5p mimics, or their respective negative controls, using Lipofectamine LTX (Invitrogen). The pRL-TK plasmid (Promega, Madison, USA) served as a normalization control. After 24 h of incubation, luciferase activity was measured using a Luciferase Reporter Assay kit (Promega) following the manufacturer’s instructions.
Statistical analyses
RNA sequencing results were expressed as the normalized fragments per kilo base per million mapped reads (FPKM) (log base per million reads) and subjected to analysis of variance (ANOVA) testing. The miRNAs with differential expression were defined based on a Benjamini–Hochberg adjusted p-value <5% and |log2 (fold-change FC) | ≥1. Significant differential miRNAs were then analyzed for target gene prediction using TargetScan (https://www.targetscan.org/vert_80/)23 and miRWalk (http://mirwalk.umm.uni-heidelberg.de).24 The resulting overlapping datasets were subsequently utilized for the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis with the aid of OmicStudio (LC-Bio). In the following experiments, we employed the Shapiro–Wilk normality test to evaluate the distribution of the data (Supplementary Table 1). The miRNA expression displayed skewness and underwent transformation using the Box–Cox method to conform to a normal distribution. Differences in miRNA and Smad2 expression between groups were assessed using t-test. Diagnostic analysis was conducted through receiver operating characteristic (ROC) curve analysis. Furthermore, after assigning ranks to miRNA-769 based on its expression values, Spearman’s correlation method was utilized to assess its association with Smad2. A p-value of less than 0.05 was considered statistically significant. The statistical software IBM SPSS v. 22.0 (IBM Corp., Armonk, USA) was utilized for data analysis.
Results
miRNAs differential expression among allo-HSCT patients and healthy controls
The demographic and clinical characteristics of the study participants are presented in Table 2. Among all 3 groups, a total of 518 miRNAs were expressed. After conducting an ANOVA analysis (details not presented), it was found that 191 of these miRNAs exhibited differential expression (Figure 1). The most significantly differentially expressed miRNAs were identified as miRNA-505-5p and miRNA-769-5p (Table 3), sorted in descending order of significance based on p-values. Enrichment analysis of these 2 miRNAs was conducted using KEGG pathway functional annotation (Figure 2, Figure 3). The KEGG results indicated that target genes of miRNA-769-5p were enriched in the TGF-β signaling pathway (Figure 3). While high-throughput sequencing indicated a modest reduction in miRNA-769-5p abundance in the oral cGVHD group compared to the without cGVHD group, the observed difference was not statistically significant.
Lower expressions of miRNA-769-5p in patients with oral cGVHD
To further validate the differential expression of miRNA-769-5p, we expanded the sample size to include 15 patients in both the oral cGVHD and without cGVHD groups (Table 4), and conducted verification through PCR analysis. It was observed that the expression of miRNA-769-5p was higher in patients without oral cGVHD compared to those with oral cGVHD (Table 5, Figure 4). Among 30 patients after allo-HSCT (both the oral cGVHD and without cGVHD groups), ROC analysis indicated that miRNA-769-5p holds a certain diagnostic value for oral cGVHD, with an area under the curve (AUC) of 0.809, signifying a moderate diagnostic efficacy (95% confidence interval (95% CI): 0.657, 0.961; p = 0.004, Figure 5).
Higher expressions of Smad2 in patients with oral cGVHD
Following predictions made by TargetScan and miRWalk, Smad2 was validated as a potential target gene of miRNA-769-5p through a dual-luciferase test (Figure 6). Among 30 patients after allo-HSCT, a negative correlation was observed between miRNA-769-5p and Smad2 mRNA (r = –0.635, 95% CI: –0.8136, –0.3460; p < 0.001, Figure 7). Additionally, both Smad2 mRNA and protein levels were higher in the patients with oral cGVHD compared to those without cGVHD (Table 5, Figure 8).
Discussion
The study successfully established an expression profile of circulating miRNAs in oral cGVHD and made initial observations indicating a phenomenon of a low expression of miRNA-769-5p and high expression of Smad2 in oral cGVHD patients.
Through the sequencing of miRNAs in peripheral blood samples of patients, we identified a total of 518 miRNAs expressed across all 3 groups. Subsequent ANOVA analysis revealed differential expression in 191 of these miRNAs. Guided by the ranking of p-values, we selected the top 2 miRNAs for KEGG pathway analysis, which unveiled an enrichment of target genes for miRNA-769-5p within the TGF-β signaling pathway. Based on the sequencing results, we expanded the sample size and conducted PCR validation, which revealed the lower expression of miRNA-769-5p in oral cGVHD. Importantly, it was ascertained that miRNA-769-5p holds diagnostic significance in the context of oral cGVHD.
Research on miRNA-769-5p has primarily centered on its involvement in cancer. Studies have revealed its function as an oncogene in gastric cancer,25 hepatocellular carcinoma26 and prostate cancer.27 Conversely, in oral squamous cell carcinoma (OSCC), elevated miRNA-769-5p expression has been found to inhibit tumor cell proliferation, migration and invasion while promoting apoptosis.28 Studies have shown that patients after allo-HSCT are 2–3 times more likely to develop solid tumors than healthy individuals.29 About 1/3 of solid tumors are secondary tumors in the skin and mucosa of recipients of transplants, with OSCC accounting for 50% of these patients.30, 31 Consequently, the question arises as to whether the low expression of miRNA-769-5p in oral cGVHD patients is one of the factors predisposing them to the development of oral cancer, a matter that warrants further consideration.
Furthermore, the study unveiled that the overexpression of miR-769-5p has the potential to alleviate tissue fibrosis by modulating the TGF-β/Smad signaling pathway.32 It is well-established that the TGF-β/Smad signaling pathway can stimulate tissue fibrosis.33 Both mRNA and protein levels were higher in the oral cGVHD group than those in patients without cGVHD. The high expression of Smad2 hinted at the activation of TGF-β/Smads, which leads to oral cGVHD. In this study, after being predicted using TargetScan and miRWalk, we verified that Smad2 was the target gene of miRNA-769-5p using the dual-luciferase test. A negative relationship between miRNA-769-5p and Smad2 mRNA was also observed. Fibrosis is a feature of oral cGVHD, and the effective treatment of cGVHD often involves inhibiting the TGF-β/Smad signaling pathway.34 Consequently, the high expression of miRNA-769-5p in patients without oral cGVHD may also act to suppress the expression of SMADs, subsequently inhibiting the activation of the TGF-β/Smad signaling pathway, ultimately leading to a reduction in fibrosis and symptom relief.
Additionally, it is noteworthy that the TGF-β/Smad signaling pathway is involved in Th17/Treg immunoreactions.35 In patients suffering from skin cGVHD, the study identified a dysregulation in cytokine secretion, such as TGF-β, alongside an imbalance in Th17/Treg cell ratios.36 In our earlier research, we also observed imbalances in cytokines related to Th17 and Treg cells in oral cGVHD.17 Previous studies have demonstrated that downregulating Smad2 can mitigate the severity of cGVHD20 by not only reducing fibrosis37 but also by inhibiting the generation of Th17.38 The investigations further emphasized that by inhibiting the TGF-β/Smad pathway, Th17-mediated immune responses could be effectively controlled, thus ameliorating murine cardiac transplant rejection.39 Therefore, the high expression of miRNA-769-5p in patients without oral cGVHD may serve to modulate the immune response by suppressing the TGF-β/Smad pathway, ultimately reducing the manifestation of oral cGVHD symptoms.
Limitations
In this study, we have examined the distinctive expression patterns of miRNAs in patients affected by oral cGVHD. We have also undertaken an initial exploration of the potential relationship between miRNA-769-5p and Smad2 differential expression and oral cGVHD. Nevertheless, the precise mechanistic interactions of miRNAs influencing the development of diseases require further in-depth exploration and validation through in vitro experiments. Additionally, during the sequencing phase, we observed lower expression levels of miRNA-769-5p in the peripheral blood of healthy volunteers compared to individuals with oral cGVHD and those without it. However, this study did not proceed to validate these findings using PCR. Subsequent research could involve further validation and exploration of the implications of elevated miRNA-769-5p expression in post-transplant patients.
Conclusions
Differential expression of miRNAs, particularly the downregulation of miRNA-769-5p, may influence the development of oral cGVHD by diminishing its inhibitory effect on the TGF-β/Smad signaling pathway through its interaction with Smad2. This interaction potentially plays a significant role in the pathogenesis of oral cGVHD and offers a promising avenue for further research and therapeutic exploration.
Supplementary data
The Supplementary materials are available at https://doi.org/10.5281/zenodo.10551750. The package includes the following file:
Supplementary Table 1. Normality test for the data distribution.
Data availability
The datasets generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.
Consent for publication
Not applicable.