Expression of Transcription Factors in CD4 + T Cells as Potential Biomarkers of Motor Complications in Parkinson's Disease.

BACKGROUND: Management of motor complications (MC) represents a major challenge in the long-term treatment of Parkinson's disease (PD) patients. In this context, the role of peripheral adaptive immunity may provide new insights, since neuroinflammatory mechanisms have been proved crucial in the disease.
OBJECTIVE: The aim of this study was to analyze the transcription factors genes involved in CD4 + T cells development to uncover specific molecular signatures in patients with (PMC) and without (WMC) motor complications.
METHODS: mRNA levels of CD4 + T lymphocytes transcription factor genes TBX21, STAT1, STAT3, STAT4, STAT6, RORC, GATA3, FOXP3, and NR4A2 were measured from 40 PD patients, divided into two groups according to motor complications. Also, 40 age- and sex-matched healthy controls were enrolled.
RESULTS: WMC patients had higher levels of STAT1 and NR4A2 (p = 0.004; p = 0.003), whereas in PMC we found higher levels of STAT6 (p = 0.04). Also, a ROC curve analysis confirmed STAT1 and NR4A2 as feasible biomarkers to discriminate WMC (AUC = 0.76, 95%CI 0.59-0.92, p = 0.005; AUC = 0.75, 95%CI 0.58-0.90, p = 0.007). Similarly, STAT6 detected PMC patients (AUC = 0.69, 95%CI 0.52-0.86, p = 0.037).
CONCLUSION: These results provide evidence of different molecular signatures in CD 4 + T cells of PD patients with and without MC, thus suggesting their potential as biomarkers of MC development.

INTRODUCTION

Levodopa represents the main symptomatic treatment for Parkinson’s disease (PD). Nonetheless, after a variable time, its administration can lead to the development of complications (MC) such as motor fluctuations (MF) and levodopa-induced dyskinesias (LIDs). MF, defined as the worsening or reappearance of motor symptoms (often in parallel with low levodopa plasma concentration), include wearing-off, delayed-ON, no-ON, random ON-OFF, and early morning or nocturnal akinesia [1]. LIDs comprise different hyperkinetic involuntary movements (most commonly chorea and dystonia), which can be peak dose-related or diphasic. Both MF and LIDs have a great impact and affect patients’ quality of life [1, 2]. Known predictors for MC development are the presence of nonmotor features as well as medications and demographic factors (e.g., younger age at PD onset and lower body mass index) [3]. The pathological background of MC is complex and not yet fully understood. Interestingly, peripheral adaptive immunity is known to be a key regulator of neuroinflammation both in the central nervous system and in the periphery: accordingly, on one side T cells can be found in the substantia nigra of PD brains and on the other CD4 + T lymphocytes subsets from peripheral blood may have a crucial role as well [4, 5].

We recently reported reduced CD4 + T cells in PD patients [6] with relatively increased Th1 resulting in a putative pro-inflammatory Th1 bias. Remarkably, the peripheral immune imbalance observed in PD patients was recapitulated by a peculiar transcription factor (TF) gene expression profile in CD4 + T cells [6]. More specifically, in comparison with control subjects, PD patients had lower levels of TBX21, STAT3, STAT4, RORC, and NR4A2, and higher levels of STAT6, GATA3, and FOXP3. A similar peculiar molecular signature in CD4 + T cells, strongly resembling cells from PD patients, was reported in idiopathic REM sleep behavior disorder (iRBD): in both iRBD subjects and PD patients lower levels of TBX21, STAT3, and STAT4, and higher levels of FOXP3 were reported. iRBD represents the strongest risk factor for prodromal PD, and these findings suggested early involvement of peripheral immunity in the disease [7].

Whether neuroinflammation may be involved in disease progression, and in particular in the development of MC, is however still unclear. Boi and coll-eagues showed in an animal model that thalidomide and 3,6’-dithiothalidomide significantly attenuate the severity of L-dopa induced dyskinesias by reducing tumor necrosis factor (TNF)-α levels in the striatum and substantia nigra pars reticulata and increasing the levels of interleukin (IL)-10 [8]. Moreover, amantadine (a noncompetitive antagonist of the N-methyl-D-aspartate glutamate receptor commonly used to treat LIDs) can act on microglia by inhibiting its inflammatory activation [9]. Based on available evidence, it has been recently suggested that in PD abnormally activated microglia and astrocytes lead to altered neuronal-glial communication, thus affecting synaptic activity and neuroplasticity and contributing to the development of LIDs [10]. Similarly, a study from Teema et al. [11] showed in a rat model that treatment with ibuprofen or piroxicam in combination with l-dopa ameliorated wearing-off at the end of week 10, delayed the development of dyskinesia, and decreased striatal cyclooxygenase-2 (COX-2) and vascular endothelial growth factor (VEGF) levels. Taken together, these results point to a potential role of the immune system in the development of MC by affecting both presynaptic and postsynaptic mechanisms.

Although neuroinflammation is increasingly reg-arded as a key mechanism in MC occurring during PD progression, no information exists so far on the possible involvement of peripheral immunity. In the present study we, therefore, decided to explore in CD4 + T cells of PD patients with and without MC, the TF gene expression profile which we previously reported as dysregulated in both established PD as well as in iRBD [6, 7]. The identification of suitable molecular signatures could indeed be crucial for the early identification of subjects at risk of MC development, as well as to find new therapeutic targets and eventually delay MC onset in PD.

MATERIALS AND METHODS

Subjects

This is a retrospective case-control study: patients with an established diagnosis of idiopathic PD [12] were selected from an electronic database established in the context of a research program aimed at the study of peripheral immunity in PD [6, 7, 13]. The database currently includes 205 PD patients and 94 healthy controls (HC). At the time of enrollment, all patients underwent neurological assessment, performed by experienced neurologists in Movement Disorders. Motor symptoms and disease staging were assessed in “ON” condition using respectively the Unified Parkinson’s Disease Rating Scale (UPDRS) part III and the Hoehn and Yahr scale [14, 15]. Other relevant data were collected, including duration and type of antiparkinsonian treatment, age at PD onset, disease duration, interval between PD diagnosis and MC onset, interval between MC onset and blood withdrawal, and levodopa equivalent daily dose or LEDD [16]. For the aforementioned research program protocol, established exclusion criteria were: a history of chronic autoimmune diseases or cancer and administration of immunomodulatory treatment. A peripheral blood venous sample was obtained at enrollment and used to assess complete blood count and peripheral CD4 + T lymphocytes TF gene expression.

For the present study, we carefully revised the records of the 205 PD patients: drug-naïve patients (n = 69) and subjects with a later diagnosis of atypical parkinsonism (n = 4) were excluded from the analysis. Only drug-treated patients with comprehensive clinical and CD4 + T cells gene expression data available were included (n = 132). We subsequently selected those with motor complications (PMC), defined as those subjects who reported one or more of the following (wearing-off, delayed-ON, no-ON, random ON-OFF, early morning/nocturnal akinesia, peak-dose or diphasic dyskinesias). Of the 25 PMC identified, 5 patients were excluded because of moderate/severe dementia (Mini-Mental State Examination < 19 [17]). For each PMC, we selected a PD patient without motor complications (WMC), nearly perfectly matched with respect to sex, UPDRS part III score (±2 points), and duration of antiparkinsonian and levodopa treatment (±2 years). Finally, HC were also selected and matched to each PMC and WMC patient, according to sex and age (±2 years). The institutional Ethics Committee approved the study protocol (CE 65/16) and all subjects gave their informed written consent. The study was conducted according to the Declaration of Helsinki and following the international research ethical principles involving human subjects.

Isolation of CD4 + T cells and real-time PCR

Detailed procedures for isolation of peripheral blood mononuclear cells (PBMC) from whole blood, immunomagnetic isolation of CD4 + T cells, and total RNA extraction and reverse transcription, are described elsewhere [6, 7]. Expression of the TF genes TBX21, STAT1, STAT3, STAT4, STAT6, RORC, GATA3, FOXP3, and NR4A2 was measured by Real-Time PCR. To start Real-Time PCR reactions 2μl aliquots were obtained (cDNA final concentration: 1μg/μl) and the following thermal protocol was used: 20 s at 95°C (x 1, hot start); 2-step cycles as follows: 1 s at 95°C, 20 s at 60°C (x 40). Assays were performed in triplicate for each sample, and levels of mRNA were finally expressed as 2-DCt where DCt = [Ct (sample) - Ct (housekeeping gene)]. Relative expression was determined by normalization to the expression of RPS18, which is the gene for 18S cDNA.

Statistical analysis

Variables were expressed as counts and percentages when categorical and as mean±standard deviation when continuous. Normal distribution of data was assessed using the Shapiro-Wilk test. Differences between groups were analyzed through ANOVA analysis of variance, after testing for the ind-ependence of observations, normality of residuals, homogeneity of variances (Levene statistics) and checking for outliers, or with non-parametric Mann-Whitney U test if these assumptions were not met. Comparisons between categorical variables were ass-essed using Fisher’s exact test. When adjusting for covariates in the two-way ANCOVA model, ass-umptions for homogeneity of regression slopes and absence of interaction between each covariate and factors were tested. A receiver operating characteristic (ROC) curve analysis was carried out to esta-blish the biomarkers’ discriminatory power. Area under the curve (AUC) and significance values were obtained and AUC values interpretation was determined according to Mandrekar [18]. Optimal cut-offs were chosen by coordinate tracing of the ROC curve according to Youden’s index analysis. Sensitivity, specificity, positive, and negative likelihood ratios (LR+, LR-), positive and negative predictive values (PPV, NPV) were computed. Biomarkers’ combination was explored by estimating predicted probabilities from a logistic regression model and using the values obtained as the test variable in the ROC analysis. The significance level was set to p < 0.05. All analyses were performed using SPSS Version 25 (IBM Corporation, Armonk, USA) and Graphpad Prism version 8 (GraphPad Software Inc., San Diego, USA).

RESULTS

For this cross-sectional study we identified 20 PMC, 20 matched WMC, and 40 age- and sex-matched HC. Complete clinical and demographic data are shown in Table 1. PMC reported wearing off (100%), morning and nocturnal akinesia (15%), delayed or no-ON/random ON-OFF (20%), and peak-dose dyskinesias (45%). Fisher’s exact test and Mann-Whitney U test didn’t reveal any statistically significant difference between groups for any variable. Complete blood count showed slightly increased values of mean corpuscular hemoglobin (MCH) and mean corpuscular hemoglobin concentration (MCHC) in PMC compared to WMC (data not shown). Detailed data of TF mRNA levels in PMC, WMC and HC are included in Supplementary Table 1. Compared to HC, PD patients displayed significantly higher levels of STAT1, GATA3, STAT6, and FOXP3. When considering TF mRNA levels between PD subgroups, we found significantly higher levels of STAT1 and NR4A2 in WMC (respectively p = 0.004, p = 0.003), whereas PMC had higher levels of STAT6 (p = 0.04). Complete Real-Time PCR data, including TBX21/STAT4 (pro-Th1), STAT3/RORC (pro-Th17), and GATA3 (pro-Th2), can be seen in Fig. 1. When controlling in a two-way ANCOVA model for remaining potential confounders such as age, sex [19, 20], and disease duration, between-group differences were statistically significant for STAT1 and NR4A2 (see Table 2). We, therefore, evaluated the ability of TF mRNA levels in CD4 + T cells to discriminate between PD subgroups using ROC curve analysis: STAT1 and NR4A2 showed good AUC values (respectively 0.76, 95%CI 0.59–0.92, p = 0.005; 0.75, 95%CI 0.58–0.90, p = 0.007) thus efficiently identifying WMC. The combination of STAT1 and NR4A2 did not increase significantly the discrimination of WMC (AUC = 0.77, p = 0.003). On the other hand, STAT6 allowed MC identification (AUC = 0.69, 95%CI 0.52–0.86, p = 0.037), see Fig. 2. Choosing for STAT1 an optimal cut-off value of 1.99E–004, the LR+and LR- were respectively 12 and 0.4 (sensitivity = 60%, 95%CI 38.6%–78.1%; specificity = 95%, 95%CI 76.4%–99.7%; PPV = 92.3%, NPV = 70.4%). Considering for NR4A2 a cut-off value of 1.09E–004, LR+was 5.5 and the LR–0.5 (sensitivity = 55%, 95%CI 34.2%–74.2%; specificity = 90%, 95%CI 69.9%–98.2%; PPV = 84.6%, NPV = 66.7%). The optimal cut-off value for STAT6 was found at 9.157E–006, with a LR+value of 2.5 and a LR- value of 1.6 (sensitivity = 50%, 95%CI 29.93%–70.07%; specificity = 80%, 95%CI 58.40%–91.93%; PPV = 71%, NPV = 61%). Cut-off values for each TF with sensitivity, specificity, LR+, LR-, PPV and NPV are summarized in Table 3.

Table 1

Clinical and demographic characteristics of PD patients and healthy controls

Number of patients (N = 40)	PMC (N = 20)	WMC (N = 20)	HC (N = 40)	P PMC vs WMC	P PD vs HC
Age (years), mean (SD)	65.75 (10.3)	70 (8.6)	69.2 (8.6)	0.13	0.519
Sex (M/F)	12/8	14/6	26/14	0.74	1.0
Age at PD onset (years), mean (SD)	57.15 (9.5)	62.65 (7.8)		0.07
Disease duration, mean (SD)	8.60 (4.9)	7.60 (3.6)		0.62
UPDRS part III, mean (SD)	15 (7.5)	16.3 (6)		0.36
Hoehn and Yahr stage
stage 1, n (%)	8 (40%)	13 (65%)		0.20
stage 2, n (%)	10 (50%)	6 (30%)		0.33
stage 3-4, n (%)	2 (10%)	1 (5%)		0.50
Antiparkinsonian therapy
Levodopa, n (%)	19 (95%)	17 (85%)		0.60
Levodopa+COMTI, n (%)	6 (30%)	2 (10%)		0.23
DA, n (%)	12 (60%)	14 (70%)		0.74
MAOIs, n (%)	12 (60%)	12 (60%)		1.0
Duration of antiparkinsonian treatment (years), mean (SD)	7.25 (4.5)	6.2 (3.99)		0.47
Duration of levodopa treatment (years), mean (SD)	4.28 (2.7)	4.85 (2.8)		0.51
LEDD (mg/die), mean (SD)	745.5 (353)	625.5 (218.8)		0.18
Interval PD diagnosis-MC onset (years), mean (SD)	5.1 (2.7)
Interval MC onset-blood withdrawal (years), mean (SD)	2.31 (1)
Motor complications, n (%)	Wearing off, 20 (100%); morning/nocturnal akinesia, 3 (15%); delayed-ON/no- ON/random ON-OFF, 4 (20%); peak dose dyskinesias, 9 (45%)

COMTI, Catechol-O-methyltransferase inhibitors; DA, dopamine agonists; HC, healthy controls; LEDD, levodopa equivalent daily dose; MAOIs, monoamine oxidase inhibitors; MC, motor complications; PD, Parkinson’s disease; PMC, patients with motor complications; WMC, patients without motor complications; UPDRS, Unified Parkinson’s Disease Rating Scale.

Fig. 1

Transcription factors mRNA levels in CD4 + T lymphocytes in PD patients with (PMC) and without (WMC) motor complications. Data are plotted as medians with 25°–75° percentiles (boxes) and min–max values (whiskers).

Table 2

Two-way ANCOVA for STAT1, NR4A2, and STAT6 adding age, disease duration, and sex as covariates in the model

	F	p	Partial eta squared
STAT1
Age	F(1,34) = 0.051	0.823	0.001
Disease duration	F(1,34) = 1.683	0.203	0.047
Sex	F(1,34) = 0.037	0.848	0.001
Groups	F(1,34) = 5.844	0.021	0.147
NR4A2
Age	F(1,34) = 0.616	0.438	0.018
Disease duration	F(1,34) = 4.948	0.033	0.127
Sex	F(1,34) = 0.856	0.361	0.025
Groups	F(1,34) = 6.243	0.017	0.155
STAT6
Age	F(1,34) = 0.182	0.672	0.050
Disease duration	F(1,34) = 2.174	0.150	0.060
Sex	F(1,34) = 0.234	0.631	0.007
Groups	F(1,34) = 3.091	0.088	0.083

Fig. 2

ROC curves of transcription factors mRNA levels as candidate biomarkers to discriminate between patients with and without motor complications. Also AUC and p values are shown.

Table 3

Sensitivity, specificity, positive likelihood ratio (LR+), negative likelihood ratio (LR-), positive predictive value (PPV), and negative predictive value (NPV) for established cut-offs

Genes	Cut-off value	Sensitivity (95%CI)	Specificity (95%CI)	LR+	LR-	PPV	NPV
STAT1	1.99E-004	60%(38.6–78.1)	95%(76.4–99.7)	12	0.4	92.3%	70.4%
NR4A2	1.09E-004	55%(34.2–74.2)	90%(69.9–98.2)	5.5	0.5	84.6%	66.7%
STAT6	9.157E-006	50%(29.93–70)	80%(58.40–91.9)	2.5	1.6	71%	61%

DISCUSSION

This study provides evidence that PMC and WMC exhibit striking differences in lymphocytes TF. Distinct molecular patterns have previously been observed in drug-naïve/drug-treated PD patients and healthy controls using the same TF gene panel [6] but no significant data were found concerning disease progression. It is well recognized that STAT1, together with TBX21 and STAT4, drives Th1 differentiation whereas STAT3/RORC are involved in Th17, STAT6/GATA3 in Th2, and FOXP3/NR4A2 in Treg development [21-26]. More in detail, STAT1 and STAT3, in concert with Jmjd3 gene, activate microglia thus driving the production of neurotoxic molecules such as proinflammatory cytokines, chemokines, and nitric oxide [27]. The role of NR4A2 is fairly complex since it is crucial for the development and specification of midbrain dopamine neurons [28] and Le et al. found that its reduced expression increased the vulnerability of mesencephalic dopaminergic neurons to N-methyl–4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced injury [29]. Moreover, Montarolo et al. observed that in peripheral blood mononuclear cells (PBMC) from PD patients there is a marked down-regulated expression of the whole NR4A family subsets (NR4A1, NRA4A2, NRA4A3) [30].

The present study, on one side confirms the role of CD4 + T lymphocytes in discriminating PD patients from HC, while on the other points out novel and intriguing findings. We demonstrated that remarkable differences in STAT1, NR4A2, and STAT6 expression can be detected in PD patients with and without MC. In our previous study [6] these TF genes were similar in drug-treated and drug-naïve patients, thus opening the question of whether therapy or disease progression may influence their expression. Accordingly, the confounding effect of longer disease duration in our PD group (8.10±4.3 years) has been carefully addressed in the ANCOVA model, confirming the prominent role of STAT1 and NR4A2 in identifying WMC subgroup.

The role of neuroinflammation in MC has been explored in various preclinical studies: Barnum et al. tested the anti-dyskinetic effect of corticosterone in 6-hydroxydopamine (6-OHDA)-lesioned rats (made dyskinetic by l-DOPA chronic treatment) and observed that l-DOPA-treated animals showed increased striatal expression of pro-inflammatory cytokines, particularly interleukin-1beta (IL-1β), which was prevented by corticosterone administration [31]. Another study found that there are strong correlations between LIDs occurrence and expression of a pro-inflammatory microglia phenotype: in 6-OHDA-lesioned rats, pulsatile l-DOPA administration but not continuous subcutaneous l-DOPA infusion determined progressive development of abnormal involuntary movements [32].

Nevertheless, evidence from human studies is lacking and the association between peripheral immunity and long-term treatment complications is still elusive. In this scenario, the present study provides for the first time distinctive molecular patterns that discriminate PMC and WMC. Among limitations, the retrospective case-control design, the small sample size which limited subgroup analysis in the PMC group, and missing information related to circulating CD4 + T lymphocytes should be mentioned. However, the limited number of patients recruited was due to the strict matching procedure performed to control for epidemiological and clinical variables: for each PMC we selected a paired patient WMC with similar features. In addition, even though gene expression does not necessarily mirror gene transcript activity and should be interpreted just as a general indication of CD4 + T cells involvement, in previous studies we observed that lower TBX21, STAT3, and STAT4 and higher FOXP3 were associated with reduced Th2, Th17, and Treg, with a relative increase of Th1 and increased production of interferon-γ and tumor necrosis factor (TNF)-α [6]. Based on this scenario, we therefore suggest that the present findings about TF mRNA levels in PD-associated MC likely underlie specific phenotypic and functional immune profiles possibly peculiar of MC, which deserve careful assessment in future investigations.

The present study addressed only molecular patterns in CD4 + T cells from PD patients, providing evidence that STAT1 and NR4A2 show respectively excellent and good positive likelihood ratios and better specificity than sensitivity in identifying the lack of motor complications: a thorough analysis of circulating CD4 + T cells subsets and longitudinal prospective studies are now warranted to elucidate the potential immunological background of long-term treatment complications in PD.

CONCLUSIONS

To the best of our knowledge, this is the first study providing distinctive molecular signatures of PD patients with and without motor complications. Though exploratory, these results shed light on the suitable involvement of the peripheral immune system, thus opening new landscapes in the therapeutic management of PD patients.

CONFLICT OF INTEREST

The authors declare that they have no relevant conflict of interest.

Click here for additional data file.