1 Fairfax Centre, OX5 2PB Gosford and Water Eaton, UK
2 Department of Psychiatry and Psychotherapy, Section for Dementia Research, University of Tübingen, 72076 Tübingen, Germany
3 Research Center on Aging, Graduate Program in Immunology, Faculty of Medicine and Health Sciences, University of Sherbrooke, Sherbrooke, QC H3A 0B8, Canada
4 Hertie Institute for Clinical Brain Research, Department of Neurodegeneration, University of Tübingen and German Center for Neurodegenerative Diseases, 72076 Tübingen, Germany
5 Institute of Medical Virology and Epidemiology of Viral Diseases, University of Tübingen, 72076 Tübingen, Germany
6 Department of Neurology, Kiel University, 24118 Kiel, Germany
7 Department of Immunology, University of Tübingen, 72076 Tübingen, Germany
8 Cancer Solutions Program, Health Sciences North Research Institute, Sudbury, ON P3E 3B7, Canada
Abstract
Parkinson’s Disease (PD) is associated with dysregulated/chronic inflammation. The immune system has multiple roles including beneficial effects such as clearing alpha synuclein aggregates. However, peripheral immune cells entering the brain may also contribute to inflammation and neurodegeneration. To identify which cells might have a negative impact and could be potential therapeutic targets, we compared immune signatures of patients and healthy controls.
Multicolor flow cytometry was used to determine the frequencies of major immune cell subsets in peripheral blood mononuclear cells (PBMCs) of PD patients and controls. Because of the major impact of Cytomegalovirus (CMV) infection on the distribution of immune cell subsets, particularly cluster of differentiation (CD)8+ T-cells, all participants were tested for CMV seropositivity.
Although the cohort of 35 PD patients exhibited the well-established T-cell differentiation signature driven by CMV infection, there were no differences in the frequencies of differentiated or pro-inflammatory T-cells, B-cells or natural killer cells (NK-cells) attributable to the disease. However, percentages of myeloid-derived suppressor cells (MDSCs) were higher in PD patients than controls. Moreover, percentages of CD14+CD16+ (intermediate) monocytes expressing the C-C chemokine receptor type 5 (CCR5) correlated with disease severity assessed by the Movement Disorder Society’s revised version of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS) score and disease duration.
A comprehensive evaluation of the major subsets of circulating immune cells in PD patients revealed differences in myeloid cells between PD and healthy controls and some correlation of monocyte abundance with disease severity.
Keywords
- Parkinson’s disease
- CMV
- T-cells
- B-cells
- NK-cells
- MDSC
- CCR5
- monocytes
- inflammation
- neurodegeneration
- peripheral immune cells
Parkinson’s Disease (PD) is characterised by an age-associated progressive
degeneration of dopaminergic neurons in the substantia nigra [1] and the
accumulation of alpha-synuclein (
To be able to address the most promising candidates for immunotherapy, a
comprehensive overview of immune signatures in PD compared to healthy controls is
needed, especially as there is clear evidence in the literature that subsets of
inflammatory cells are associated with PD pathology. For example, the levels of
peripheral dendritic cells, especially those secreting pro-inflammatory
cytokines, were found to be altered in PD patients compared to controls and an
impact on disease severity was proposed [10, 11]. A neuroprotective role of
natural killer cells (NK-cells) (by supporting
Infection with the persistent herpesvirus cytomegalovirus (CMV) and high anti-CMV Immunoglobulin G (IgG) antibody titers are related to an inverted cluster of differentiation (CD)4:CD8 T-cell ratio andmay be associated with age-independent morbidity, mortality and cognitive and functional decline [16]. Another point to consider: T-helper 17 (Th17) cells affect BBB integrity [17] which can promote enhanced migration of peripheral cells into the brain and mediation of dopaminergic neurodegeneration as shown in a mouse model [18]. Therefore we investigated the capacity of T-cells to produce cytokines and determined T-helper cell profiles in the context of CMV-serostatus and PD. Taken together, functional and immunophentyping analysis leads to a better understanding of the complexity of PD immunity. Our study revealed an example of the immune network in PD with similar T-cell profiles between PD patients and controls and on the other hand a correlation of the frequencies of C-C chemokine receptor type 5 (CCR5)-expressing CD14+CD16+ monocytes with disease severity and disease duration.
PD patients without signs of dementia and their partners as controls free from any neurological disease and living in the same environment were recruited from the outpatient clinic at the Neurodegenerative Department of the University of Tübingen, Germany. Clinical data were collected (Table 1) with patients’ written consent and approval by the Ethics Committee of the Medical Faculty of the University of Tübingen (project number 480/2015BO2) [19, 20]. PD patients fulfilled the United Kingdom Parkinson’s Disease Society (UKPDS) Brain Bank Diagnostic criteria [21]. The severity of motor malfunctions and disease progression was assessed with the Hoehn and Yahr [22] score and the Movement Disorders Society-revised version of the Unified Parkinson’s Disease Rating scale (MDS-UPDRS) [23].
| Controls | Controls | PD patients | PD patients | |
| CMV serostatus | CMV- | CMV+ | CMV- | CMV+ |
| Number | 9 | 13 | 17 | 18 |
| Median age | 68 | 60 | 68 | 68 |
| Age range | 58–74 | 53–77 | 55–74 | 43–73 |
| Gender | 3F/6M | 6F/7M | 7F/10M | 6F/12M |
| Median MDS-UPDRS I | 2 | 0 | 7 | 6 |
| Median MDS-UPDRS II | 0 | 0 | 10 | 6 |
| Median MDS-UPDRS III | 1 | 1 | 29 | 28 |
| Median MDS-UPDRS IV | 0 | 0 | 0 | 0 |
| Median total MDS-UPDRS | 2 | 2 | 45 | 42 |
| Total MDS-UPDRS range | 0–13 | 0–11 | 24–112 | 14–68 |
| H&Y score range | 0–2 | 0–1 | 1–4 | 1–4 |
| Median disease duration | 9 years | 5 years | ||
| Disease duration range | 2–13 years | 1–10 years |
CMV, Cytomegalovirus; PD, Parkinson’s Disease; MDS-UPDRS, Movement Disorder Society-revised version of the Unified Parkinson’s Disease Rating scale; H&Y, Hoehn and Yahr Scale; M, male; F, female.
Plasma was collected and CMV serostatus determined via a recombinant CMV IgG
immunoblot (Mikrogen, Neuried, Germany) using six different targets (IE1, p150,
CM2, p65, gB1 and gB2). Peripheral blood mononuclear cells (PBMCs) were isolated from whole ethylenediaminetetraacetic acid (EDTA) blood using
FicoLite FicoLite-H(human) (Cat.No.: GTF1511YK, Linaris, Dossenheim, Germany) gradient centrifugation. PBMCs in 10% dimethylsulfoxide (DMSO) (research grade, Cat.No.: 20385.01,
Serva, Heidelberg, Germany) and
20% fetalbovineSerum (FBS) (Cat.No.: F0804, Sigma, Darmstadt, Germany) were frozen at –80 °C and subsequently stored in
liquid nitrogen. Cells were thawed and after blocking non-specific antibody
binding with 1% Gamunex (human IgG from Talecris, Research Triangle Park, NC, USA) dead cells were stained with
ethidium monoazide bromide (Cat.No.: 40015, Biotium, Fremont, CA, USA). Different surface marker stainings in phosphatebufferedsaline
(PBS) (Cat.No.: H15-002, PAA, Paching, Austria) with 2% (v/v) FBS (Cat.No.: F0804, Sigma, Darmstadt, Germany), 2 mM EDTA (Cat.No.: 11278, Serva, Heidelberg, Germany) and 0.01% (w/v) sodiumazide (NaN3) (Ca.No.: S-8032, Sigma, Darmstadt, Germany) were performed using the antibodies shown in Supplementary Table
1. To assess T-cell function 0.5
To investigate possible interactions between the systemic inflammation in PD described in the literature [5] and T-cell differentiation, the expression of different surface markers was analysed. To exclude the confounding effects of CMV infection, patients and controls were dichotomized into seropositives and seronegatives. First, the overall frequencies of CD4+ and CD8+ T-cells were quantified. Higher CD8+ T-cell frequencies in CMV+ individuals were found while patients and controls did not display different patterns except a slightly higher frequency of CD4+CD8+ T-cells (T-cells associated with chronic inflammatory disorders [24]) in PD patients. All differences between groups are summarized in Table 2. Also regarding the specific T-cell phenotypes the differences were limited to CMV- compared to CMV+ study participants, but not affected by disease. T-cells from patients and controls showed a similar differentiation profile. Naïve T-cells were characterised as CD45RO-CD45RAhiCD27+CD28+ and late-differentiated T-cells as CD45RO+CD45RA-CD27-CD28- or CD95+ or CD57+. Applying the main published models [25] we found a shift from early- to late-differentiated T-cells in CMV+ compared to CMV- individuals, although not all marker combinations led to statistically significant results and CD8+CD45RO+ T-cells showed contrary results with slightly lower frequencies in CMV+ than CMV- PD patients. The differences between CMV- and CMV+ individuals were more pronounced for CD4+ T-cells than CD8+ T-cells and in patients than controls.
| Population | Parental | C– vs. C+ | PD– vs. PD+ | C– vs. PD– | C+ vs. PD+ |
| N | 7 vs. 11 | 14 vs. 14 | 7 vs. 14 | 11 vs. 14 | |
| CD4+ | CD3+ | − | ↓ | − | − |
| CD8+ | CD3+ | ↑↑ | ↑↑ | − | − |
| CD4+:CD8+ | CD3+ | − | ↓↓ | − | − |
| CD4–CD8– | CD3+ | − | − | − | − |
| CD4+CD8+ | CD3+ | ↑ | − | ↑ | − |
| CD27+CD28+ | CD4+ | ↓↓ | ↓↓ | − | − |
| CD27–CD28– | CD4+ | ↑↑ | ↑↑↑ | − | − |
| CD45RO–CD45RAhi | CD4+ | − | − | − | − |
| CD45RO+CD45RA– | CD4+ | − | − | − | − |
| CD45RO–CD45RAhiCD27+CD28+ | CD4+ | − | − | − | − |
| CD45RO+CD45RA–CD27–CD28– | CD4+ | ↑↑ | ↑↑↑ | − | − |
| CD57+ | CD4+ | ↑↑ | ↑↑ | − | − |
| CD95+ | CD4+ | − | − | − | − |
| CD27+CD28+ | CD8+ | − | ↓↓ | − | − |
| CD27–CD28– | CD8+ | ↑ | ↑↑↑ | − | − |
| CD27+CD28– | CD8+ | − | − | − | − |
| CD27–CD28+ | CD8+ | − | − | − | − |
| CD45RO–CD45RAhi | CD8+ | − | − | − | − |
| CD45RO+CD45RA– | CD8+ | − | ↓ | − | − |
| CD45RO+CD45RA+ | CD8+ | − | ↓ | − | − |
| CD45RO–CD45RAhiCD27+CD28+ | CD8+ | − | − | − | − |
| CD45RO+CD45RA–CD27–CD28– | CD8+ | − | − | − | − |
| CD57+ | CD8+ | − | ↑ | − | − |
| CD57+ | CD4–CD8– | − | − | − | − |
Footnotes: − p
To detect key players possibly contributing to systemic inflammation in PD,
PBMCs were mitogenically stimulated and intracellular cytokines analysed by
flow-cytometry for a panel of pro- and anti-inflammatory cytokines. CMV drives
T-cell differentiation, as demonstrated by phenotypic analysis and also by
functionality. However, the only statistically significant differences between
groups in this analysis were limited to a higher frequency of unstimulated (basal
ex-vivo) CD4+ T-cells producing pro-inflammatory interferon gamma (IFN
First, the frequencies of all B-cells (CD19+) and plasmablasts (CD19+CD20–) within the total leukocyte population (CD45+) were analysed. Then CD27+CD43+ B1 cell frequencies were investigated within CD19+ and CD20+ cells. Next, CD24 and CD38 expression was plotted against each other and double-positive cells gated to define transitional B-cells (CD38hiCD24+). Further, CD38–CD24–, CD38hiCD27+ and CD40hi B-cells were gated. To determine differentiation states, IgD expression was plotted against CD27. IgD+CD27– were considered early- and double-negative B-cells late-differentiated. No significant differences between any of the populations analysed were found. These results are summarized in Supplementary Table 3.
We next analysed the frequencies of the different major innate leukocytes but again found no differences regarding immature NK-cells (CD56hiCD3–), mature NK-cells (CD56loCD3–) or NKT-like cells (CD14–CD3+CD16+). NKT-cell (CD56+CD3+) frequencies were higher in CMV+ than in CMV– controls and also in PD patients but decreased over the disease duration. No differences were observed when comparing patients and controls. Monocyte subsets (CD3–CD20–CD56–CD14+CD16–, CD3–CD20–CD56–CD14+CD16+, CD3–CD20–CD56–CD14loCD16+) within CD45+ PBMCs also revealed no differences between groups. However percentages of CD14–CD3–CD16–CD20+ cells negatively correlated with total MDS-UPDRS values and disease duration in CMV- PD patients. CCR5 expression was not different on B-cells, but frequencies of CCR5+CD14+CD16+ monocytes positively correlated with disease duration and total MDS-UPDRS values. Several different myeloid-derived suppressor cells (MDSCs) were also different between the groups as follows: CD14+CD15–CD11blo human leucocyte antigen(HLA)-DR- higher in CMV+ vs. CMV– PD patients, CD14–HLA-DR–CD33–CD11blo lower in CMV+ vs. CMV– controls and higher in CMV– PD patients vs. controls. Notably, frequencies of CD14+CD15+HLA-DR–CD11b+CD33+ MDSCs correlated positively with total MDS-UPDRS values in CMV+ PD patients (Table 3).
| Population | C– vs. C+ | PD– vs. PD+ | C– vs. PD– | C+ vs. PD+ | PD– dis.duration | PD+ dis.duration. | PD– MDS-UPDRS | PD+ MDS-UPDRS |
| n | 9 vs. 13 | 17 vs. 17 | 9 vs. 17 | 13 vs. 17 | 17 | 17 | 17 | 17 |
| CD56+CD3hi | − | − | − | − | − | − | − | − |
| CD56+CD3+ | ↑↑ | ↑ | − | − | − | ↓ | − | − |
| CD56hiCD3– | − | − | − | − | − | − | − | − |
| CD56loCD3– | − | − | − | − | − | − | − | − |
| CD14–CD3–CD16+ | − | − | − | − | − | − | − | − |
| CD14–CD3–CD16+HLA-DR– | − | − | − | − | − | − | − | − |
| CD14–CD3–CD16+HLA-DR+ | − | − | − | − | − | − | − | − |
| CD14–CD3+CD16+ | − | − | − | − | − | − | − | − |
| CD14–CD3–CD16– | − | − | − | − | − | − | ↓ | − |
| CD14–CD3+CD16– | − | − | − | − | − | − | − | − |
| CD14–CD3+CD16–CCR5+ | − | − | − | − | − | − | − | ↓ |
| CD14–CD3–CD16–CD20+ | − | − | − | − | ↓ | − | ↓ | − |
| CD20+CCR5+ | − | − | − | − | − | − | − | − |
| lin–CD14+CD16+ | − | − | − | − | − | − | − | − |
| lin–CD14+CD16– | − | − | − | − | − | − | − | − |
| lin–CD14loCD16+ | − | − | − | − | − | − | − | − |
| lin–CD14–CD16– | − | − | − | − | − | − | − | ↓ |
| CD14+CD16–CCR5+ | − | − | − | − | − | − | − | − |
| CD14+CD16+CCR5+ | − | − | − | − | ↑↑ | − | ↑↑ | − |
| CD14–HLA-DR– | − | − | − | − | − | − | − | − |
| CD14–HLA-DR–CD33–CD11blo | − | ↓↓ | ↑ | − | − | − | − | − |
| CD15hi | − | − | − | − | − | − | − | − |
| CD14–CD15+ | − | − | − | − | − | − | − | − |
| CD14–CD15+CD11blo | − | − | − | − | − | − | − | − |
| CD14+CD15+ | − | − | − | − | − | − | − | − |
| CD14+CD15+HLA-DR– | − | − | − | − | − | − | − | − |
| CD14+CD15+HLA-DR–CD11b+CD33+ | − | − | − | − | − | − | − | ↑ |
| CD14–CD15– | − | − | − | − | − | − | − | − |
| CD14+CD15– | − | − | − | ↑ | − | − | − | − |
| CD14+CD15-CD11bloHLA-DR– | − | − | − | ↑ | − | ↓ | − | − |
| CD11b+CD33– | − | − | − | − | − | − | − | − |
| CD11b+CD33+ | − | − | − | − | − | − | − | ↑ |
| CD11b–CD33– | − | − | − | − | − | − | − | − |
| CD11b–CD33+ | − | − | − | − | − | − | − | − |
Footnotes: – p
Systemic inflammation has been documented in several neurodegenerative diseases. Peripheral immune profiles differ between patients and controls in several published studies, mostly in Alzheimer’s Disease (AD) where a shift from early- to late-differentiated T-cells is observed [25, 26, 27]. In this disease also higher frequencies of activated T-cells have been reported [25] as well as T-cells expressing chemokine receptors, especially CCR6 [27], which is a key player for brain-homing [28]. To the best of our knowledge, for PD the situation is less clear, especially as previous studies generally failed to account for the effects of CMV infection [15, 29, 30, 31]. Therefore we performed comprehensive phenotyping of the main T-cell subsets as well as other immune cell types in peripheral blood of PD patients and a balanced control cohort. Because CMV has been shown to have a great impact on T-cell differentiation profiles [32], we analysed CMV– and CMV+ study participants separately. The previously reported decrease in the frequency of naïve T-cells and the accumulation of late-stage differentiated T-cells in CMV+ compared to CMV- individuals [32] could be confirmed in the present study. In PD patients the difference between CMV-seronegative and seropositive individuals was slighter greater, but this does not necessarily mean that CMV can have a greater impact due to an exhausted immune system in PD. While CMV– and CMV+ PD patients were age-matched, in the controls the CMV– participants were slightly older than the CMV-seropositives. It is known that chronological age is also independently associated with greater T-cell differentiation status [33, 34] and it could be that this counteracted the effect of CMV in our study, such that fewer differences between CMV– and CMV+ individuals were seen. That the disease does not play a main role can be definitely concluded, because comparing patients and controls revealed no significant differences regarding the phenotypes of either CD4+ or CD8+ T-cell subsets. Other studies have also reported no relevant differences in CD4+ phenotypes, but a shift to more early differentiated CD8+ T-cells, even in the context of CMV [35, 36, 37].
In PD, a pro-inflammatory environment was reported, so we asked whether this might be associated with higher frequencies of cytokine-producing T-cells as the source in patients. However, the frequencies of anti-inflammatory cytokine-positive T-cells were not lower in PD patients and neither were the levels of pro-inflammatory cytokine-producing T-cells higher. A possible reason why a shift from early- to late-differentiated T-cells was observed in AD [25, 26, 27], but not here in PD patients might be the different stimulatory environment; cerebrospinal fluid (CSF) of PD patients was reported to differ from CSF of AD patients and healthy controls in respect to selected proteins [38], which could act as immune stimulatory antigens or induce cell death. This could lead to immune dysfunction as reported for Tregs in PD patients [39]. Calopa et al. [40] detected an increase in apoptosis of peripheral CD4+ T-cells in PD patients. As a consequence the level of CD4+ T-cells was lower in patients than controls. This was mainly related to cells which were considered as naïve (CD4+CD45RA+), while the frequencies of CD25+ activated CD4+ T-cells were higher in PD patients [40]. In the present study, lower frequencies of CD4+ T-cells within the total CD3+ T-cells were observed only in CMV+ patients relative to CMV+ controls. A reason for the difference between our results and those of Calopa et al. [40] could be that they referred to the lymphocyte count (decreased in PD) and we only analysed frequencies. Further they revealed higher CD95 (FAS) expression on CD4+ T-cells of PD patients than controls [40] and another study found an association between CD45RO+CD95+CD4+ effector memory T-cells and the motor part of the MDS-UPDRS III score, independent of age and disease duration [39]. We also failed to find correlations with disease duration, but also with MDS-UPDRS values, although we used total MDS-UPDRS values and not only MDS-UPDRS III.
In parallel to T-cell phenotypes, we determined frequencies of Th17 cells in vitro as well as other CD4+ and CD8+ T-cells producing pro- or anti-inflammatory cytokines, but again found no significant differences between patients and controls or between CMV- and CMV+ individuals as reported for Th17 cells in other studies [41]. Tregs did counteract this effect in mice [18], so possibly regulatory cells play an important role in PD. In our study frequencies of CD14+CD15–and CD14+CD15–CD11bloHLA-DR– cells were higher in PD patients compared to controls, but only in CMV+ patients. In addition, it is noteworthy that other MDSC subsets, CD14+CD15+HLA-DR–CD11b+CD33+ and CD11b+CD33+ cells, correlated positively with total MDS-UPDRS values, but only in CMV+ donors. This could be a sign that CMV activates the immune system which as a consequence can react more readily to other targets, such as PD-related molecules. Other leucocyte subsets with suppressive function such as CD4+ CD25+ regulatory T-cells can reduce neuroinflammation and production of reactive oxygen species via IL-10 and neurotrophic factors (reviewed in [5]). Frequencies of Forkhead-Box-Protein P3 (FoxP3)+CD4+ regulatory T-cells were reported to be higher in older than in younger people, but not different in AD or PD patients. On the other hand, the suppressive activity of these regulatory T-cells was stronger in a previous patient cohort [42].
In PD patients and animal models increased neurodegeneration has been associated
with systemic inflammation originating in the periphery [5, 43]. In this context,
peripheral immune cells have been shown to be activated, infiltrate the brain and
contribute to neuroinflammation [6, 44] and neurodegeneration, as shown in a
mouse model of PD [45]. The hypothesis in this latter model was that peripheral
inflammation starts the cascade and promotes activation of microglia, which
induces oxidative stress in dopaminergic neurons. As a consequence,
Monocytic cells might reduce disease burden by phagocytosis of
To comprehensively assess immune signatures, we also investigated the frequencies and phenotypes of other leukocytes including B-cells and NK-cells, but found no significant differences between either CMV- and CMV+ individuals or between PD patients and controls. However, this does not necessarily imply that they play no role in PD pathology. For example, higher NK-cell activation as reflected by natural killer cell 2D (NKG2D) expression was reported in early PD [31], but we did not test this marker. For the B-cell phenotypes too few samples for a correlation with disease severity were available, but the frequencies of total CD20+ cells within the CD45+ correlated negatively with disease duration and total MDS-UPDRS values. This was also only seen in CMV-seronegative patients. In another study age-related decrease of B-cell numbers was reported, but only for the very old over 85 years of age [49]. As reported for T-cells, age is also a driver of B-cell differentiation and higher frequencies of late-differentiated IgG+IgD–CD27– B-cells, with a reduced capacity to respond to new antigens (maybe also PD-related antigens) were found in older adults [50]. This finding was related to CD19+CD38–CD24– B-cells, which show low expression of immunosenescence-associated marker CD180, produce pro-inflammatory TNF, and whose frequencies were higher in older individuals [51]. Because PD pathology might cause B-cell stress, we investigated different subsets and phenotypes, but found no impact of the disease.
From a comprehensive analysis of immune signatures in CMV-seronegative and CMV-seropositive PD patients, significantly higher frequencies of certain MDSCs were found in PD patients than controls, and CCR5 expression correlated with total MDS-UPDRS values hinting at higher cell migration with increased disease severity. No significant differences in frequencies of other immune cells emerged from this study, overall implying that the interaction between the (circulating) immune system and brain pathology/function in PD is very specific and potentially best visible in CMV-negative PD patients. This study was considered a pilot study to identify PD immune signatures covering a broad overview of immune cells in a limited number of patients. Further studies are needed to confirm these findings in larger cohorts. Thereby a focus on MDSC populations, chemokine receptors and regulatory T-cells (previously found relevant) could be considered.
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
DG, TF, DB, WM and GP designed the research study. DG, LO performed the research. LO, CS and WM managed the biobank. CS supported the study as biobank manager and recruiter of study participants and organising clinical data. KH performed CMV-serostatus testing. DG and LO analyzed the data. DG and GP wrote the manuscript. All authors contributed to editorial changes in the manuscript. All authors read and approved the final manuscript. All authors have participated sufficiently in the work and agreed to be accountable for all aspects of the work.
Study participants were recruited by the outpatient clinic at the Neurodegenerative Department of the University of Tübingen, Germany and clinical data and blood samples collected and biobanked with patients’ written consent. The study was carried out in accordance with the guidelines of the Declaration of Helsinki. The Ethics Committee of the Medical Faculty of the University of Tübingen approved this study under project number 480/2015BO2.
We thank the patients and their families for agreeing to participate in this study.
This research received no external funding.
The authors declare no conflict of interest.
Supplementary material associated with this article can be found, in the online version, at https://doi.org/10.31083/JIN26393.
References
Publisher’s Note: IMR Press stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
