Expression of Foxp3 and Rorc Genes in Patients Undergoing Allogeneic Hematopoietic Cell Transplantation: Impact in the Development of Acute and Chronic Graft-Versus-Host Disease

Citation: Cabral CMV, Braga WMT, Arantes AM, de Oliveira JSR, Colleoni GWB, et al. (2018) Expression of Foxp3 and Rorc Genes in Patients Undergoing Allogeneic Hematopoietic Cell Transplantation: Impact in the Development of Acute and Chronic Graft-Versus-Host Disease. Int J Blood Res Disord 5:032. doi.org/10.23937/2469-5696/1410032 Accepted: September 05, 2018: Published: September 07, 2018 Copyright: © 2018 Cabral CMV, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Introduction
Graft versus host disease (GVHD) is a leading cause of morbidity and mortality following hematopoietic cell transplantation (HCT) [1]. Despite prophylactic immunosuppression, GVHD occurs in up to 50% of the patients, and innovative strategies to prevent and treat it are urgently needed. GVHD involves dysregulation of inflammatory cytokine cascades and donor cellular response against host antigens [2,3]. CD4 + T cells play a critical role in mediating GVHD and understanding the Abstract The aim of this study is to characterize Treg (FOXP3) and Th17 (RORC)-related genes in patients who underwent hematopoietic cell transplantation (HCT) and correlate them with graft-versus-host disease (GVHD) onset and patients' outcome. Fifty-one patients undergoing allogeneic HCT from matched related donors (MRD) were studied. Samples were collected before HCT, at acute (aGVHD) or chronic GVHD (cGVHD) onset, and 90 days post-HCT for those who did not have signs of aGVHD. FOXP3 and RORC gene expressions were assessed by real-time PCR (qPCR) as representatives of Treg and Th17 cells, respectively. We found increased expression of proinflammatory Th-17 related RORC in cGVHD compared with patient without GVHD; also, we found decrease in Treg-related FOXP3 expression in the samples at the time of GVHD when compared to its respective samples before HCT. Although these results could have therapeutic relevance, we could not identify genes' expression values to be used to predict GVHD onset, overall survival, non-relapse or relapse-related mortality.
role of each subset of these lymphocytes in the development of GVHD is critical to discover effective therapy and improve HCT outcome [4,5].
The challenge in allogeneic HCT is to find a balance between beneficial and harmful T-cell effects [6]. CD4 + CD25 + FOXP3 + Regulatory T (Treg) cells are key immunoregulatory cell population involved in the maintenance of immune tolerance and have been used to prevent and treat GVHD [4,7]. Murine GVHD models have demonstrated that depletion of CD25 + Tregs promotes GVHD [8,9] while, conversely, infusion of in vitro expanded CD4 + CD25 + Tregs can either prevent or treat GVHD [10,11]. Similar to mouse models, Tregs in the peripheral blood of human HCT recipients have been reported to be lower in the presence of acute GVHD (aGVHD) [4,12,13]. However, in chronic GVHD (cGVHD), this correlation between increased Tregs and decreased risk of GVHD has been absent [14,15] or even reversed [16].
In this context, it is clearly necessary to understand the role of these effector cells and target antigens for pathogenesis of graft-versus-tumor (GVT) effect [24]. Standardization of new techniques, such as the use of molecular biology as a tool to detect the behavior of the immune system in the post-HCT, particularly of T lymphocyte subpopulations, can be helpful to adopt cell therapy and is of fundamental importance to the development of new prophylactic and therapeutic strategies for this clinical syndrome. Therefore, the aim of this study was to characterize the expression of FOXP3 and RORC genes in peripheral blood of 51 patients who underwent allogeneic HCT from human leukocyte antigen (HLA) identical matched related donors (MRD) and correlate them with the development of acute and chronic GVHD. Furthermore, we aimed to determine the utility of both genes as GVHD biomarkers and their impact in clinical outcome after HCT.

Patients and transplantation procedures
It is a retrospective analysis of blood samples from a biorepository consisting of peripheral blood samples from 89 consecutive patients undergoing allogeneic HCT from MRD, from May 2004 to August 2007, in two Brazilians centers. All samples were collected within 24 hours of clinical signs of aGVHD or cGVHD onset, prior to initiation of corticosteroid therapy, after informed consent signature. In addition, samples were collected before HCT and at day +90 in patients who did not develop aGVHD. Among 89 patients, 51 were eligible for this study. Exclusion criteria were: Sepsis, sinusoidal obstruction syndrome, cytomegalovirus infection or others untreated infections during bone marrow recovery, benign hematologic diseases, relapsed leukemia, death, and failure to collect/technical problems. GVHD was histologically confirmed in 88% of patients and was staged using modified Gluksberg criteria for aGVHD [25] and NIH 2005 Criteria for cGVHD [26]. All patients underwent allogeneic HCT from MRD, with unmanipulated grafts. Thirty-five peripheral blood samples of blood donor volunteers were used as the control group.

RNA extraction
Mononuclear cells from peripheral blood were separated by Ficoll-Hypaque®, according to the manufacturer's instructions, from 10 mL peripheral blood collected in EDTA. Immediately after this, and without previous sorting of specific cell types, total RNA extraction was performed using Trizol reagent (Invitrogen, Carslbad, CA, USA) following the manufacturer's instructions and stored at -80 °C until the time of use. The quality of samples was confirmed in 2% agarose gel stained with Sybr Safe. RNA concentrations were evaluated in spectrophotometer (NanoDrop, Wilmington, Delaware, USA).

Real time quantitative PCR (qPCR)
In this study qPCR technique was used due to the availability of an RNA Biobank from peripheral blood samples of patients undergoing allogeneic HCT at different times, before and after transplantation. PCR was performed using FOXP3 (Hs01085834_m1) and RORC (Hs01076112_m1) specific primers and probes using the 7500 Fast Real Time PCR System (Applied Biosystems, Foster City, CA). Briefly, experiments were performed using universal PCR Master Mix (Applied Biosystems, Foster City, CA), by Taqman method. To calculate the expression level of each gene in the samples, relative quantification was performed using the 2 -ΔΔCT method [27]. Results were normalized based on the amplification of three genes: ACTB, BGUS e GAPDH (TaqMan® Endogenous Control). All samples were analyzed in triplicate [28].

Statistical analysis
To compare FOXP3 and RORC expression in the different groups we used nonparametric Kruskal-Wallis test and, in cases of significant differences between the groups, we performed multiple comparisons by Mann-Whitney test with Bonferroni corrections. To evaluate the differences between the expression of FOXP3 and RORC in the same patient before transplant and at the time of GVHD, we used non-parametric Wilcoxon test for paired data. ROC curves were constructed to evaluate the sensitivity, specificity, and positive and negative predictive values of FOXP3 and RORC expression for acute and chronic GVHD. Overall survival (OS) was calculated from the date of transplant until death or last follow-up. Actuarial probabilities of OS were estimated according to the Kaplan-Meier method and the curves were compared using the log-rank test. Differences with a p < 0.05 were statistically significant. Statistical analysis was performed using SPSS Software version 14.0 (SPSS Inc., Chicago, IL) and Graphpad Prism 5 software (GraphPad, San Diego, CA).

Patients and transplant characteristics
Fifty-one adult allogeneic HCT recipients were enrolled in this study. Demographic and transplantation information are summarized in ( Table 1). The median age of the patients was 36 years, ranging from 16 to 69 years. 53% were male and 47% female. Acute leukemia/

B.
A.

FOXP3 and RORC expression and GVHD correlation
We found statistically significant differences when healthy controls, no GVHD, aGVHD and cGVHD where analyzed for FOXP3 (p = 0.006) and RORC expression (p = 0.0042) (Figure 1, Kruskal-Wallis test followed by Mann-Whitney test, with Bonferroni posttest). Among them, the clinically relevant information is: The increase of FOXP3 in cGVHD group compared with aGVHD and patients without GVHD ( Figure 1A) and the increase of RORC expression in cGVHD compared with patient without GVHD ( Figure 1B).
After this global analysis, samples of sixteen patients before HCT and their respective samples at aGVHD onset were evaluated in a paired analysis Figure 2. The expression of FOXP3 decreased at aGVHD onset (p = 0.0085, Figure 2C), whereas no increase was observed in RORC expression at aGVHD onset (p > 0.05, Figure  2B). Similarly, when performing paired analysis of eight patients who developed cGVHD, there was a decrease of FOXP3 expression (p = 0.0391, Figure 2B) at cGVHD onset but we did not find any statistically significant correlation with RORC expression at cGVHD onset (p > 0.05, Figure 3D).

FOXP3 and RORC expression as GVHD biomarkers
ROC curves were constructed to evaluate FOXP3 and RORC expression for aGVHD and cGVHD, using samples before GVHD as controls. The AUC value obtained for relative expression of FOXP3 in aGVHD was 0.554, with p > 0.05. Therefore, we do not find a marker capable of differentiating transplanted patients and pre-GVHD controls. The same occurred for both genes, FOXP3 and RORC, in aGVHD or cGVHD ( Table 2). In this scenario, the ideal value of FOXP3 and RORC relative expression to express high sensitivity and specificity was not found ( Table 2).

FOXP3 and RORC expression and HCT outcomes
Overall survival (OS) at 1 year was 60.8%. During the follow up, there were 29 deaths: Eight (28%) related to disease relapse and 21 (78%) attributable to transplantation-related events. Five (24%) patients died from aGVHD, nine (43%) from cGVHD complications and seven (33%) from infection. OS curves were built using the not develop any GVHD, 10 (20%) showed only aGVHD, 18 (35%) cGVHD and 13 (25%) patients developed both   The higher expression of FOXP3 at cGVHD onset compared with controls, no GVHD and aGVHD contradicts several studies, which showed a decreased frequency of Treg from peripheral blood of patients with cGVHD [29,30]. However, the most clinically relevant information in this study is the increase of RORC expression in cGVHD compared to patient without GVHD. There are few studies in the literature investigating the role of Th17 lymphocytes in cGVHD. Dander, et al. [31] also observed a significant increase in the number of Th17 cells in the blood of patients with active cGVHD. However, other studies in humans have failed to demonstrate the increase expression of Th17 cells in peripheral blood after HCT at the time of GVHD [32,33].
We performed a paired analysis of sixteen patients with their respective samples before HCT and at the time of aGVHD and we observed a statistically significant decrease in FOXP3 expression at aGVHD onset, confirming previous clinical studies [12,15,34,35]. Similarly, we conducted a matched analysis among patients who developed cGVHD, and also decreased expression of FOXP3 was found, consistent with some published work [29,[36][37][38]. There was no difference in the expression of RORC at the time of aGVHD and cGVHD. Dlubek, et al. [39] investigated the presence of IL-17 producing cells in peripheral blood of patients undergoing HCT. Although all patients had shown increase of cells producing IL-17, a minority of patients developed aGVHD and these had lower proportions of CD4 + cells producing IL-17 on the median FOXP3 (1.67) and RORC (1.10) expression values, measured at the time of aGVHD. The small number of patients within cGVHD group did not allow such analysis. FOXP3 expression ( Figure 3A, Figure 3B and Figure  3C) and RORC ( Figure 3D, Figure 3E and Figure 3F) did not present impact on OS and relapse-related or not related mortality.

Discussion
Our study confirmed the role of Th17 and Treg related genes in the pathogenesis of GVHD. We found increased expression of RORC in cGVHD compared to patient with no GVHD. We also found a decreased FOXP3 expression at the time of GVHD onset when compared to its respective samples before HCT, in patients undergoing allogeneic HCT from MRD. Although these results could have therapeutic relevance, we could not identify genes' expression values to be used to predict GVHD onset, overall survival, non-relapse or relapse-related mortality.
Initially, performed an analysis of all patients assigned in this study (N = 51) to assess the expression of FOXP3 and RORC in different times proposed: At aGVHD and cGVHD onset and at D + 90 for patients who did not develop any GVHD. This last moment was chosen with the aim of having more homogeneous groups and before the possibility of occurring cGVHD. Our results demonstrated a similar behavior of both genes between healthy controls and patients who did not develop GVHD.  ulation of CD CD4 + FOXP3 + CD25 + by flow cytometry [4]. qPCR has also the advantage of being highly sensitive and needs a small amount of material, relevant features when it is intended to study a numerically small cell population or when searching for a biomarker [46]. aGVHD onset. One hypothesis suggested by the author was that these cells could have migrated to the affected tissues during the clinical manifestations of GVHD, justifying its absence in the peripheral blood. Other study evaluated not only the peripheral blood as well as the skin of patients with cutaneous aGVHD through biopsy and also failed to show an increase in IL-17 producing cells in both sites [32]. Fulton, et al. [40] demonstrated in a murine model that the use of RORC-deficient donor T cells results in attenuated aGVHD severity and lethality. Thus, although there is strong evidence of involvement of Th17 cells in the pathophysiology of acute and chronic GVHD, further elucidation of the mechanisms of action of these cells in GVHD is needed [41].
GVHD remains one of the major determinants of survival, post-transplant morbidity and quality of life of patients submitted to allogeneic HCT, especially when immunosuppression is required to control this complication. The presence of moderate to severe aGVHD itself confers adverse prognosis for patients, since half of them do not respond to traditional therapy based on steroids. Therefore, it is essential to develop tools to predict and better understand the biological mechanisms responsible for this type of immune response [37].
An effective immune response is required for successful tumor eradication and numerous strategies are being tested to enhance specificity, maximize GVT, and minimize GVHD after adoptive immunotherapy. In this scenario, Treg immunotherapy is a promising strategy for the therapeutic control of GVHD [7,24]. Several early-phase clinical trials were performed where Tregs were given to HCT recipients to prevent GVHD [42][43][44]. Although a limited number of patients were treated in these trials, the initial data suggested that Tregs diminish GVHD without increasing the risk of disease recurrence or infection. Interestingly, despite their clinical use, it is still unclear how Tregs function to prevent GVHD. Regulatory T cells may have adverse effect: Tregs are a major concern in cancer immunology where they have documented inhibitory activity on antitumor immunity. Just as T-cell depletion has reduced the graft-versus-malignancy effects of allogeneic hematopoietic transplantation, infusion of Tregs may potentially increase the risk of recurrent malignancy. Therefore, this can only be assessed in larger scale human trials and the main challenge in hematopoietic transplantation remains the prevention of GVHD without impairing graft-versus malignancy effects [45].
One limitation of this study is the small number of cases. Other limitation is that most published studies have used flow cytometry to characterize Treg and Th17 lymphocytes populations, not qPCR. It is noteworthy that FOXP3 transcripts analysis in CD4 + population by qPCR showed comparable results on the same pop-