- 1Department of Human Neuroscience, Sapienza University of Rome, Rome, Italy
- 2Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS) Neuromed, Isernia, Italy
- 3Department of Clinical and Experimental Sciences Section, Brescia University, Brescia, Italy
Tourette syndrome (TS) and early-onset obsessive-compulsive disorder (OCD) are frequently associated and conceptualized as distinct phenotypes of a common disease spectrum. However, the nature of their relationship is still largely unknown on a pathophysiological level. In this study, early structural white matter (WM) changes investigated through diffusion tensor imaging (DTI) were compared across four groups of drug-naïve children: TS-pure (n = 16), TS+OCD (n = 14), OCD (n = 10), and 11 age-matched controls. We analyzed five WM tracts of interest, i.e., cortico-spinal tract (CST), anterior thalamic radiations (ATR), inferior longitudinal fasciculus (ILF), corpus callosum (CC), and cingulum and evaluated correlations of DTI changes to symptom severity. Compared to controls, TS-pure and TS+OCD showed a comparable pattern of increased fractional anisotropy (FA) in CST, ATR, ILF and CC, with FA changes displaying negative correlation to tic severity. Conversely, in OCD, FA decreased in all WM tracts (except for the cingulum) compared to controls and negatively correlated to symptoms. We demonstrate different early WM microstructural alterations in children with TS-pure/TS+OCD as opposed to OCD. Our findings support the conceptualization of TS+OCD as a subtype of TS while suggesting that OCD is characterized by independent pathophysiological mechanisms affecting WM development.
Introduction
The clinical hallmark of Tourette syndrome (TS) is the presence of multiple motor and phonic tics lasting at least 1 year (1). Although tics are the defining feature, TS is a multifaceted neurodevelopmental disorder in which several other neuropsychiatric conditions may occur (2). Obsessive-compulsive disorder (OCD) is present in up to 61.5% of individuals with TS (3) and is characterized by higher genetic heritability and disease severity as compared to the adult form (4, 5). Evidence is accumulating in support for the common genetic liability of TS and OCD (6, 7), which has raised new interest for their comorbidity (TS+OCD), in terms of both clinical course and treatment. However, although TS and OCD are both thought to reflect brain network disorders (8–10), their clinical association and overlap have been poorly understood on a pathophysiological level, possibly because of the lack of studies comparing the comorbid condition with the “pure” phenotypes.
To fill this knowledge gap, we have recently used diffusion tensor imaging (DTI) and resting-state functional MRI (rsfMRI) to investigate structural as well as functional brain changes in drug-naïve children with TS, TS+OCD, and OCD (11, 12). Those studies demonstrated that TS+OCD shares common neural correlates with TS, but not with OCD, in terms of structural WM changes and functional connectivity abnormalities in fronto-cerebellar and fronto-parietal networks (11, 12). However, our previous DTI analysis (12) focused on cerebellar connections by studying WM microstructural integrity of the three cerebellar peduncles, thus not clarifying the possible pathophysiological contribution of early changes in supratentorial WM bundles in TS and OCD.
In the present study, we investigated microstructural changes in those brain WM bundles that previous studies found to be altered in TS and/or OCD (13–23). We studied pediatric drug-naïve TS patients without comorbidities (TS-pure), TS patients with OCD (TS+OCD) and pure OCD patients in comparison to a control group, to provide evidence on primary WM abnormalities and their relationship to clinical severity. Based on our prior works on the same patient cohort (11, 12), we hypothesize that: (i) TS-pure and TS+OCD share common microstructural WM alterations and may be regarded as a unitary group, and (ii) TS-pure/TS+OCD and OCD are underpinned by different WM changes. We also hypothesize that the degree of WM structural integrity in patients with TS and OCD may be predicted by severity measures of tics and obsessive-compulsive symptoms. The study has been designed to control for several confounding factors such as disease duration, comorbidities, medication, and age. In particular, pediatric drug-naïve populations with definite comorbidity profiles are crucial to identify primary neuroimaging correlates in neurodevelopment disorders such as TS and OCD.
Methods
Participants
Seventy drug-naïve children were recruited from a specialized outpatient clinic for TS and related disorders (Child and Adolescent Neuropsychiatry Unit, Department of Human Neurosciences, Sapienza University of Rome, Italy). Patients were diagnosed according to DSM-5 criteria by a child neuropsychiatrist experienced in TS, OCD, and related disorders. Nineteen participants received a diagnosis of TS-pure, 19 of TS+OCD, 17 of OCD. As controls, we selected 15 age-matched participants with episodic tension headache who were headache-free during the MRI scans. Symptom severity was assessed using the Yale Global Tic Severity Scale [YGTSS, total tic severity scale: max. 50, without impairment score, (24)] and the Children's Yale-Brown Obsessive-Compulsive Scale [CY-BOCS, (25)]. Children diagnosed with TS and CY-BOCS scores above 14 were identified as patients with the comorbid condition and included in the TS+OCD group. All subjects had normal cognitive profile (IQ ≥ 70) according to Wechsler Intelligence Scale for Children III (WISC-III) full scale and were right-handed, as assessed by the Edinburgh handedness inventory (26). Other developmental disorders (e.g., Attention-Deficit/Hyperactivity Disorder, ADHD) or psychiatric disorders were excluded using the Schedule for Affective Disorders and Schizophrenia [K-SADS-PL, (27)] administered to both parents. Parents or guardians provided written informed consent. The study was approved by the Institutional Review Board and conformed to the Declaration of Helsinki.
MRI parameters
After clinical evaluation, all participants underwent an MRI scan performed with 3.0 T scanner (Magnetic Verio; Siemens, Erlangen, Germany) with a 12-channel head coil designed for parallel imaging (GRAPPA, Generalized Autocalibrating Partial Parallel Acquisition) using a standardized protocol. Noise reduction headphones were used to prune the scanner noise. Head positioning was standardized using canthomeatal landmarks, and the head motion was minimized by inserting foam pads between the participant's head and the head coil. Moreover, participants were instructed to lie as still as possible during and between scans. MRI included: (i) Diffusion tensor imaging (DTI, single-shot echo-planar spin-echo sequence with 30 directions, TR = 12.200 ms, TE = 94 ms, FOV = 192 mm2, matrix = 96 × 96, b = 0 and 1,000 s/mm2, axial 2-mm-thick slices, no gap); (ii) 3D T1-weighted MPRAGE (TR = 1,900 ms, TE = 2.93 ms, 176 sagittal 1-mm-thick sections, without a gap, flip angle = 9°, FOV = 260 mm2, matrix = 256 × 256); (iii) dual turbo spin-echo proton density and T2- weighted images (DP-T2) (TR = 3,320 ms, TE = 10/103 ms, FOV = 220 mm2, matrix = 384 x 384, 25 axial 4-mm thick slices, 30% gap) acquired to rule out possible concomitant brain focal lesions.
Data analysis
Preprocessing
MRI preprocessing was performed via FMRIB's software library (FSL), version 5.0.9 (http://fsl.fmrib.ox.ac.uk). Considering that the prevalence of motion in the child population is high in neuroimaging data (28), the DTI scan of each subject was first visually inspected and entirely removed in presence of artifacts even in one single volume alone. Secondly, using the FSL tools “eddy” and “topup,” correction of the susceptibility-by-movement interactions was performed. Eddy current corrected DTI files were further subjected to DTI Fit using FMRIB's diffusion toolbox to generate fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) maps.
Tract–based spatial statistics
Firstly, FA maps were subjected to the Tract-Based Spatial Statistical (TBSS) tool for the voxel-wise statistical analysis (28). To deal with age-dependent skull characteristics, we created a pediatric skeleton template by averaging the FA maps, thresholded at mean FA = 0.2 and binarized, of all participants in the same space of a reference subject. Mean FA derived from the FA of the patients was therefore used as a pediatric template. In order to run the TBSS, linear registration was automatically performed to register FA images of all participants onto the pediatric template. Similarly, MD, AD, and RD maps were registered to the pediatric template to obtain comparable images in the same standardized space.
We selected five WM tracts chosen to be of interest (TOIs) (29) according to previous findings showing consistent WM bundles alterations in TS or OCD patients (14, 17, 20, 22, 30, 31). The five TOIs, i.e., anterior thalamic radiations (ATR), corpus callosum (CC), cortico-spinal tract (CST), inferior longitudinal fasciculus (ILF), and cingulum, were first selected in MNI standard space using the JHU ICBM-DTI-81 White-Matter atlas and then were linearly registered onto the pediatric template space. The TOIs in the pediatric template space were further used to create binarized masks for further voxel-wise inter-group analysis.
Statistical analysis
We performed the Shapiro-Wilk normality test to check for normal distribution of the demographic and clinical data. One-way ANOVA and post-hoc unpaired t-tests (two-tailed, alpha = 0.05, unequal variance) were performed to investigate differences among groups with respect to age. The Chi-square test was also exploited to check for sex distribution among the groups. Differences in clinical scores among pure TS, TS+OCD, and OCD were analyzed with the Mann-Whitney U test. Analyses were performed with SPSS (Statistical Package for the Social Sciences, https://www.ibm.com/analytics/spss-statistics-software).
Using the FSL Randomize (n = 5,000 permutations), non-parametric statistics were performed to investigate voxel-wise differences within the five pre-selected microstructural WM tracts between pairs of the study cohorts and to further compute correlations with clinical measures (i.e., YGTSS and CYBOCS). Statistical significance was set at p < 0.05, corrected for false discovery rate (FDR). Age and sex were included as covariates of no interest.
Results
Demographic and clinical characteristics
The clinical details of the study cohorts are summarized in Table 1. From the original population of 70 children, 51 were included in the study. Overall, 19 children were excluded due to excessive head motion (n = 10) or inability to complete the MRI scan (n = 9). Children who excessively moved their head were mainly diagnosed with TS or TS+OCD (respectively 3 and 4 cases), while two children were diagnosed with OCD, and one was an healthy control. On the contrary, 5 children with OCD, one child with TS+OCD and 3 controls could not complete the scan. Sixteen subjects were included in the TS-pure group (15 males, mean age ± standard deviation: 9.7 ± 2.1 years), 14 in TS+OCD (10 males, 10.2 ± 2.1 years old), 10 in OCD (7 males, 10.9±2.5 years old) and 11 age-matched children in the control group (2 males, 9.9 ± 1.2 years old).
TS-pure, TS+OCD, OCD and controls were not statistically different for age [F (3, 46) = 0.77, p = 0.51, partial η2 = 0.048). Conversely, chi-square test revealed uneven sex distribution between TS-pure and controls [χ2 (1, N = 27) = 15.90, p < 0.001], and between OCD and controls [χ2 (1, N = 21) = 48.90, p < 0.001]. Mann-Whitney U test revealed no significant difference in YGTSS scores between TS-pure and TS+OCD (U = 108.5, p = 0.88), as well as in CY-BOCS scores between OCD and TS+OCD (U = 65, p = 0.53).
WM microstructural alterations in TS, TS+OCD and OCD
TS-pure patients showed higher FA and lower MD, AD, RD than controls within the right ATR, the genu, body, splenium of CC, the CST, and the ILF bilaterally (Figure 1A). Since the voxel-wise results obtained by analyzing MD, AD, and RD perfectly matched those obtained by FA analysis, the relative maps are shown as Supplementary material.
Figure 1. Fractional anisotropy (FA) differences between (A) TS-pure and controls, (B) TS+OCD and controls, (C) OCD and controls, (D) TS and OCD, at anterior thalamic radiation (ATR), corpus callosum (CC), corticospinal tract (CST), inferior longitudinal fasciculus (ILF). (A): higher FA in TS-pure than in controls, (B): higher FA in TS+OCD than in controls, (C): lower FA in OCD than in controls, (D): higher FA in TS than in OCD. Results were obtained within the mask of ATR, CC, CST, and ILF. Results are presented in the whole brain FA skeleton mask derived from the complete set of participants. FA results were corrected for multiple comparisons at the false discovery rate (FDR) of p < 0.05. Red color: Higher FA differences, Blue color: Lower FA differences, TS-pure: participants with pure Tourette syndrome (TS), OCD: participants with obsessive compulsive disorder, TS+OCD: TS participants with comorbid condition, TS: participants with TS-pure and TS+OCD.
TS+OCD patients exhibited overlapping DTI features, e.g., FA, MD, AD, RD, as TS-pure in respect to controls (Figure 1B), i.e., DTI parameters were comparable in TS-pure and TS+OCD within any of the investigated WM tracts. Based on such findings, TS-pure and TS+OCD were merged for the following analyses into a joint group named “TS.”
OCD patients showed opposite DTI alterations with respect to controls, i.e., lower FA and higher MD, AD, RD within the ATR bilaterally, the genu, body, splenium of CC, the CST and ILF bilaterally (Figure 1C).
The direct comparison between TS and OCD revealed that TS had higher FA and lower MD, AD, RD than OCD within the four aforementioned WM bundles (Figure 1D).
No changes in DTI parameters were found in the cingulum in any group comparison.
White matter regions exhibiting differences in FA between groups are detailed in Table 2. MD, AD, RD alterations in all the four group comparisons are available as Supplementary Figures 1–3.
Clinical-neuroradiological correlations
In Figures 2A,B, the correlations between FA alterations with clinical scores are displayed for TS and OCD groups. For both TS and OCD, the maps of MD, AD, RD alterations and correlation with clinical severity are shown as Supplementary Figure 4.
Figure 2. Correlations of fractional anisotropy (FA) abnormalities in (A) TS with YGTSS, (B) OCD with CYBOCS at the anterior thalamic radiation (ATR), corpus callosum (CC), corticospinal tract (CST), inferior longitudinal fasciculus (ILF). Results were obtained within the mask of ATR, CC, CST, and ILF. Results are presented in the whole brain FA skeleton mask derived from the complete set of participants. Results were corrected for multiple comparisons at the false discovery rate (FDR) of p < 0.05. Blue color: Negative correlation, TS: participants with pure Tourette and comorbid condition [TS+(TS+OCD)], OCD: participants with obsessive compulsive disorder, YGTSS: Yale Global Tic Severity Scale, CYBOCS: Children's Yale Brown Obsessive-Compulsive Scale.
Concerning the severity of tics scored by the YGTSS, in the TS group, there was a negative correlation between FA and YGTSS within ATR (i.e., anterior limb of internal capsule), CC (i.e., body, splenium, genu), CST (i.e., cerebral peduncles) and ILF (i.e., posterior thalamic radiation, and inferior fronto-occipital fasciculus) (Figure 2A). Accordingly, within the same tracts, a positive correlation was found between the other diffusivity parameters, i.e., MD, AD, RD, and YGTSS. No significant correlation was observed between DTI parameters and YGTSS in any of the five selected WM tracts when TS-pure and TS+OCD were analyzed separately.
Concerning the severity of OCD symptoms as tested by the CY-BOCS, in the TS+OCD, there was no significant correlation between DTI parameters and CY-BOCS scores. In OCD, CY-BOCS negatively correlated with FA and positively with MD, AD, and RD within ATR (i.e., anterior limb of internal capsule), CC (i.e., body, splenium, genu), CST (i.e., cerebral peduncles), and ILF (i.e., posterior thalamic radiation, and inferior fronto-occipital fasciculus) (Figure 2B).
Correlational analyses of FA with clinical scales are reported in Table 3.
Table 3. Brain areas within white matter tracts showing significant correlations between fractional anisotropy (FA) and clinical data.
Discussion
The present DTI study extends our prior observations (11, 12) demonstrating abnormalities in several WM tracts in drug-naïve TS-pure, TS+OCD, and OCD. Such WM microstructural changes represent early-stage correlates, which are not affected by long disease duration, medication, or any other comorbidity.
WM microstructural alterations in TS
TS-pure and TS+OCD showed common DTI changes with respect to controls, i.e., increased FA in the corpus callosum, anterior thalamic radiations, corticospinal tract, and inferior longitudinal fasciculus. Previous studies in children with TS reported variable DTI findings, i.e., increased FA (32), reduced FA (15, 33), or no variation at all (22, 34). Despite such inconsistencies, evidence from prior works in children are in agreement with our results converging on early structural abnormalities in three major regions, i.e., interhemispheric bundles (13, 15, 20, 22, 32, 33), main motor pathways (13, 17, 20, 32) and prefrontal and fronto-striatal pathways (18, 35). In adult TS studies, results on DTI changes are also heterogeneous, showing either an increased (36, 37) and a decreased FA (38, 39) of the interhemispheric bundles and WM tracts of sensorimotor regions. The factor age is critical to explain such inconsistencies. The WM undergoes age-dependent changes throughout the lifespan (40), making comparisons between pediatric and adult cohorts hardly affordable. Moreover, given that less than 25% of individuals with TS encounter a significant persistence of tics into adulthood (41), adults may be conceived as a peculiar subpopulation with different or more pronounced neural abnormalities. Comorbidities–which were poorly controlled in former pediatric investigations–represent another relevant aspect contributing to heterogeneous findings. To control for this aspect, we have carefully selected participants based on comorbid conditions and included only treatment-naïve children. To our knowledge, the only previous pediatric study with a similar design is that of Wolff et al. (22), which did not show FA difference in the CC of boys with pure TS compared to controls. Differences in DTI analyses procedure, in particular the selection of a definite TOI, which was further partitioned into five segments (or sub-TOIs), might explain the discrepancy with our results.
The increased FA in the CC of children with TS, that we found, may be indicative of stronger structural connectivity due to increased axonal density, thicker axons, greater myelination, or a combination of these processes. In normal neurodevelopment, there is a linear relationship between interhemispheric connectivity and motor learning/control (42–44). The observed increased connectivity in callosal fibers, which was negatively correlated to tic severity, raises the intriguing possibility of enhanced interhemispheric communication between areas involved in motor control and tic inhibition in children with TS. Moreover, the presence of atypical WM structure in the ATR, which is the major fiber bundle connecting the prefrontal cortex with the thalamus through the anterior limb of the internal capsule, greatly supports the cortico-striatal-thalamic circuitry model of TS, which is considered the leading pathophysiological account of the disorder (45, 46). Lastly, increased FA was also identified in the ILF of TS children. ILF is a large associative bundle connecting the occipital with the temporal lobe, and it is critically involved in visually-guided behaviors and object recognition (47–50). Our finding of increased FA in ILF might explain the enhanced abilities in visuomotor integration and learning in patients with TS, which have been described on both behavioral (51, 52) and neurophysiological levels (53).
Overall, the observation that lower FA was associated with greater tic severity, strongly suggest an inverse relationship between WM organization in TS and disorder expression. Thus, by detecting common structural correlates, in line with our previous findings (11, 12), our data support the conceptualization of TS+OCD as a specific subtype of TS, in which obsessive-compulsive symptoms may be conceived as heterotypic manifestations of a spectrum rather than the expression of a distinct disorder.
WM microstructural alterations in OCD
Children with pure OCD manifested opposite FA changes compared to the TS group, suggesting different WM changes sustain two disorders. To date, several DTI meta-analyses have been conducted in OCD, albeit they have many inconsistencies. In a meta-analysis combining data from pediatric and adult cohorts (54), the most prominent and replicable result was a decreased FA in the genu of CC and left orbito-frontal WM. Conversely, another metanalytic investigation analyzing pediatric and adult studies separately (55) showed no FA alterations in children but decreased FA in the genu and anterior body of CC in adults. In sum, this latter study suggested the existence of different pathophysiological processes in early onset compared to adult OCD. Similarly, no FA changes were detected in a large cohort of children and adolescents in a study from the ENIGMA OCD working group (31). However, a previous study by Fitzgerald et al. (56) specifically addressed the effects of age on FA in children and adolescents with OCD. The authors showed an FA reduction of the CC in children aged 8–11 years and an FA increase in adolescents aged 16–19 years. Thus, WM organization likely undergoes different developmental phases in OCD, which are sustained by differing patterns of pruning and myelination according to age. Therefore, inconsistencies across DTI findings, as well as the null results from metanalyses, might be explained according to neurodevelopmental features, especially when pediatric broad-aged cohorts are considered.
In our study, FA changes in the anterior CC and ATR agrees with previously reported evidence of abnormal function and structure in the cortico-striatal-thalamic circuitry in OCD (57–59). The anterior CC contains WM fibers projecting to the prefrontal regions and connecting the right and left prefrontal, premotor and supplementary motor cortex, all areas repeatedly shown to have aberrant function and volume in OCD [e.g., (60)]. The ATR contains WM pathways connecting the frontal lobe and thalamus. Thus, anterior CC and ATR collectively represent key traits of cortico-striatal-thalamic circuitry, on which pathophysiological models of OCD are built. Compared to controls, OCD children also exhibited decreased FA within the ILF. Such findings indirectly suggest abnormalities in WM areas relevant for visuo-spatial abilities, which have been repeatedly reported to be impaired in OCD patients (61) and have been associated with the persistence of OCD in adult age (62).
Overall, the observed negative correlation of obsessive-compulsive symptoms with FA values in the CST, ATR, CC, and ILF indicate that more severe clinical phenotypes are underpinned by less organized WM tracts in OCD children.
From the original OCD group, a few children were affected by head motion during MRI acquisition leading to question whether FA reductions were due to motion artifacts (63). However, visual inspection procedure was performed to identify artifacts on all subjects, irrespectively of group membership (TS, OCD, or controls). Furthermore, we might expect that, unlike TS, the scans of OCD patients were less likely corrupted by motion effects as a consequence of the clinical symptoms and showed the same probability of motion artifacts as the pediatric healthy population with similar age. As a matter of fact, TS patients showed FA increments in respect to controls, therefore FA reduction in patients with OCD appears to be inherently associated with the disorder, rather than to head motion.
Differences in WM microstructural alterations in TS and OCD
When looking at the specific difference between clinical groups, a dichotomic pattern of WM abnormalities within ATR, CST, and ILF emerged in TS vs. OCD, i.e., increased FA in TS as opposed to decreased FA in OCD. This further supports the hypothesis that pure OCD represents a different entity from OCD in the context of TS, reflecting independent neuroadaptive processes. Interestingly, the TS group showed a negative correlation between FA and YGTSS, pointing to the idea that an early increase in axons, fiber density bundles and/or myelination in TS may be indicative of a compensatory reorganization in response to the disease pathophysiology. In the OCD group, the clinical-neuroradiological correlations suggest opposite considerations, as reduced fiber myelination and organization are associated with the overall disease burden.
A final comment concerns the DTI indexes other than FA such as MD, AD, and RD, collected in the TS and OCD cohorts. It is known that FA reflects myelination and organization of axon fibers (64, 65). Decreased FA would point to less myelinated and less compact WM tracts and vice versa (66, 67). However, given that the interpretation of FA changes is not univocal (68), additional DTI parameters such as MD, RD, and AD might help characterize WM microstructural abnormalities. RD and AD provide measures of myelin and axonal integrity, respectively (69, 70). In the present study, MD, RD, and AD scaled in the same direction within each clinical cohort, i.e., they all resulted decreased in TS and increased in OCD. Overall, these results further support our hypothesis of a differential organization and maturation of WM fibers rather than selective damage to axons or myelin in both TS and OCD.
Some limitations of the current study should be stressed. First, differences in sex distribution between patients and controls may have influenced the results. However, the effects of sex on our DTI findings have been minimized by including such variable as a nuisance covariate in the analysis. Thus, we believe that sex distribution has unlikely influenced our results. Second, the relatively small number of participants would preclude the immediate generalization of our findings to all children with TS and OCD. However, given the specific participants' age range and their careful selection, we believe that our results truly capture the WM microstructural organization in the early stages of TS and OCD, uninfluenced by other comorbidities, long disease duration and pharmacological treatment.
Conclusion
This is the first study to characterize and compare WM microstructure in drug-naïve TS and OCD children. We highlight a shared pattern of WM microstructural changes in pure TS and combined TS+OCD as opposed to pure OCD, pointing to the conceptualization of TS+OCD as a peculiar subtype of TS. Compared to the normative population, the overall TS group showed a unique pattern of increased FA in callosal WM and in tracts linking the frontal, occipital and temporal cortices with each other and with the thalamus. The increased WM connectivity–which inversely correlated to tic severity–may represent an adaptive reorganization to aberrant or overactive sensory-motor processing in TS, possibly allowing partial compensation of tics. Conversely, children with OCD showed widespread reduced WM connectivity of callosal, temporo-occipital, and fronto-thalamic bundles, which were all related to greater disease severity and appear to play a role in the disease pathophysiology since an early stage.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Azienda Ospedaliero-Universitaria Policlinico Umberto I. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author contributions
KB: investigation, writing–original draft, writing–review and editing, MRI acquisition, data analysis, and data curation. GC: investigation, writing–original draft, writing–review and editing, and data curation. ST and CG: writing–review and editing, MRI acquisition, and data analysis. AS and GM: investigation and writing–review and editing. FC: investigation, writing–review and editing, data curation, conceptualization, design of the study, and study supervision. PP: investigation, writing–review and editing, MRI acquisition, data analysis, conceptualization, design of the study, and study supervision. All authors contributed to the article and approved the submitted version.
Acknowledgments
We thank all patients and their families for taking part in the study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer LM declared a past collaboration with one of the authors AS to the handling editor.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2022.960979/full#supplementary-material
Supplementary Figure 1. Mean diffusivity (MD) differences between (A) TS-pure and controls, (B) TS+OCD and controls, (C) OCD and controls (D) TS and OCD at anterior thalamic radiation (ATR), corpus callosum (CC), corticospinal tract (CST), inferior longitudinal fasciculus (ILF). (A): lower MD in TS-pure than in controls, (B): lower MD in TS+OCD than in controls, (C): higher MD in OCD than in controls, (D): lower MD in TS than in OCD. Results were obtained within the mask of ATR, CC, CST, and ILF. Results were presented in the whole brain FA skeleton mask derived from the complete set of participants. MD results were corrected for multiple comparisons at the false discovery rate (FDR) of p < 0.05. Red: Higher MD differences, Blue: Lower MD differences, TS-pure: participants with pure Tourette syndrome (TS), OCD: participants with obsessive compulsive disorder, TS+OCD: TS participants with comorbid condition, TS: participants with TS-pure and TS+OCD.
Supplementary Figure 2. Axial diffusivity (AD) differences between (A) TS-pure and controls, (B) TS+OCD and controls, (C) OCD and controls (D) TS and OCD at anterior thalamic radiation (ATR), corpus callosum (CC), corticospinal tract (CST), inferior longitudinal fasciculus (ILF). (A): lower AD in TS-pure than in controls, (B): lower AD in TS+OCD than in controls, (C): higher AD in OCD than in controls, (D): lower AD in TS than in OCD. Results were obtained within the mask of ATR, CC, CST, and ILF. Results were presented in the whole brain FA skeleton mask derived from the complete set of participants. AD results were corrected for multiple comparisons at the false discovery rate (FDR) of p < 0.05. Red: Higher AD differences, Blue: Lower AD differences, TS-pure: participants with pure Tourette syndrome (TS), OCD: participants with obsessive compulsive disorder, TS+OCD: TS participants with comorbid condition, TS: participants with TS-pure and TS+OCD.
Supplementary Figure 3. Radial diffusivity (RD) differences between (A) TS-pure and controls, (B) TS+OCD and controls, (C) OCD and controls (D) TS and OCD at anterior thalamic radiation (ATR), corpus callosum (CC), corticospinal tract (CST), inferior longitudinal fasciculus (ILF). (A): lower RD in TS-pure than in controls, (B): lower RD in TS+OCD than in controls, (C): higher RD in OCD than in controls, (D): lower RD in TS than in OCD Results were obtained within the mask of ATR, CC, CST, and ILF. Results were presented in the whole brain FA skeleton mask derived from the complete set of participants. RD results were corrected for multiple comparisons at the false discovery rate (FDR) of p < 0.05. Red: Higher RD differences, Blue: Lower RD differences, TS-pure: participants with pure Tourette syndrome (TS), OCD: participants with obsessive compulsive disorder, TS+OCD: TS participants with comorbid condition, TS: participants with TS-pure and TS+OCD.
Supplementary Figure 4. Clinical correlations in TS and OCD with YGTSS and CYBOCS at the anterior thalamic radiation (ATR), corpus callosum (CC), corticospinal tract (CST), and inferior longitudinal fasciculus (ILF). (A): positive correlation in TS between MD and YGTSS, (B): positive correlation in OCD between MD and CYBOCS, (C): positive correlation in TS between AD and YGTSS, (D): positive correlation in OCD between AD and CYBOCS, (E): positive correlation in TS between RD and YGTSS, (F): positive correlation in OCD between RD and CYBOCS. Results were presented in the whole brain FA skeleton mask derived from the complete set of participants. Results were corrected for multiple comparisons at the false discovery rate (FDR) of p < 0.05. TS: participants with pure Tourette and comorbid condition [TS+(TS+OCD)]. OCD: participants with obsessive compulsive disorder, YGTSS: Yale Global Tic Severity Scale, CYBOCS: Children's Yale Brown Obsessive-Compulsive Scale.
References
1. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (DSM-5®). Washington, DC: American Psychiatric Publishing, Inc. (2013).
2. Robertson MM, Eapen V, Singer HS, Martino D, Scharf JM, Paschou P, et al. Gilles de la tourette syndrome. Nat Rev Dis Primers. (2017) 3:16097. doi: 10.1038/nrdp.2016.97
3. Kloft L, Steinel T, Kathmann N. Systematic review of co-occurring OCD and TD: Evidence for a tic-related OCD subtype? Neurosci Biobehav Rev. (2018) 95:280–314. doi: 10.1016/j.neubiorev.2018.09.021
4. Dell'Osso B, Marazziti D, Albert U, Pallanti S, Gambini O, Tundo A, et al. Parsing the phenotype of obsessive-compulsive tic disorder (OCTD): a multidisciplinary consensus. Int J Psychiatry Clin Pract. (2017) 21:156–9. doi: 10.1080/13651501.2017.1291822
5. Taylor S. Early versus late onset obsessive–compulsive disorder: evidence for distinct subtypes. Clin Psychol Rev. (2011) 31:1083–100. doi: 10.1016/j.cpr.2011.06.007
6. Huisman-van Dijk HM, van de Schoot R, Rijkeboer MM, Mathews CA, Cath DC. The relationship between tics, OC, ADHD and autism symptoms: a cross-disorder symptom analysis in Gilles de la Tourette syndrome patients and their family members. Psychiatry Res. (2016) 237:138–46. doi: 10.1016/j.psychres.2016.01.051
7. Mathews CA, Grados MA. Familiality of Tourette syndrome, obsessive-compulsive disorder, and attention-deficit/hyperactivity disorder: heritability analysis in a large sib-pair sample. J Am Acad Child Adolesc Psychiatry. (2011) 50:46–54. doi: 10.1016/j.jaac.2010.10.004
8. Schlemm E, Cheng B, Fischer F, Hilgetag C, Gerloff C, Thomalla G. Altered topology of structural brain networks in patients with Gilles de la Tourette syndrome. Sci Rep. (2017) 7:10606. doi: 10.1038/s41598-017-10920-y
9. Wen H, Liu Y, Rekik I, Wang S, Zhang J, Zhang Y, et al. Disrupted topological organization of structural networks revealed by probabilistic diffusion tractography in Tourette syndrome children. Hum Brain Mapp. (2017) 38:3988–4008. doi: 10.1002/hbm.23643
10. Yun JY, Boedhoe PSW, Vriend C, Jahanshad N, Abe Y, Ameis SH, et al. Brain structural covariance networks in obsessive-compulsive disorder: a graph analysis from the ENIGMA consortium. Brain. (2020) 143:684–700. doi: 10.1093/brain/awaa001
11. Tikoo S, Cardona F, Tommasin S, Giannì C, Conte G, Upadhyay N, et al. Resting-state functional connectivity in drug-naive pediatric patients with Tourette syndrome and obsessive-compulsive disorder. J Psychiatr Res. (2020) 129:129–40. doi: 10.1016/j.jpsychires.2020.06.021
12. Tikoo S, Suppa A, Tommasin S, Giannì C, Conte G, Mirabella G, et al. The cerebellum in drug-naive children with Tourette syndrome and obsessive–compulsive disorder. Cerebellum. (2021). doi: 10.1007/s12311-021-01327-7
13. Govindan RM, Makki MI, Wilson BJ, Behen ME, Chugani HT. Abnormal water diffusivity in corticostriatal projections in children with Tourette syndrome. Hum Brain Mapp. (2010) 31:1665–74. doi: 10.1002/hbm.20970
14. Gruner P, Vo A, Ikuta T, Mahon K, Peters BD, Malhotra AK, et al. White matter abnormalities in pediatric obsessive-compulsive disorder. Neuropsychopharmacology. (2012) 37:2730–9. doi: 10.1038/npp.2012.138
15. Jackson SR, Parkinson A, Jung J, Ryan SE, Morgan PS, Hollis C, et al. Compensatory neural reorganization in Tourette syndrome. Curr Biol. (2011) 21:580–5. doi: 10.1016/j.cub.2011.02.047
16. Jayarajan RN, Venkatasubramanian G, Viswanath B, Janardhan Reddy YC, Srinath S, Vasudev MK, et al. White matter abnormalities in children and adolescents with obsessive-compulsive disorder: a diffusion tensor imaging study. Depress Anxiety. (2012) 29:780–8. doi: 10.1002/da.21890
17. Liu Y, Miao W, Wang J, Gao P, Yin G, Zhang L, et al. Structural abnormalities in early Tourette syndrome children: a combined voxel-based morphometry and tract-based spatial statistics study. PLoS ONE. (2013) 8:e76105. doi: 10.1371/journal.pone.0076105
18. Makki MI, Munian Govindan R, Wilson BJ, Behen ME, Chugani HT. Altered fronto-striato-thalamic connectivity in children with Tourette syndrome assessed with diffusion tensor mri and probabilistic fiber tracking. J Child Neurol. (2009) 24:669–78. doi: 10.1177/0883073808327838
19. Piras F, Piras F, Chiapponi C, Girardi P, Caltagirone C, Spalletta G. Widespread structural brain changes in OCD: a systematic review of voxel-based morphometry studies. Cortex. (2015) 62:89–108. doi: 10.1016/j.cortex.2013.01.016
20. Sigurdsson HP, Pépés SE, Jackson GM, Draper A, Morgan PS, Jackson SR. Alterations in the microstructure of white matter in children and adolescents with Tourette syndrome measured using tract-based spatial statistics and probabilistic tractography. Cortex. (2018) 104:75–89. doi: 10.1016/j.cortex.2018.04.004
21. Silk T, Chen J, Seal M, Vance A. White matter abnormalities in pediatric obsessive-compulsive disorder. Psychiatry Res Neuroimaging. (2013) 213:154–60. doi: 10.1016/j.pscychresns.2013.04.003
22. Wolff N, Luehr I, Sender J, Ehrlich S, Schmidt-Samoa C, Dechent P, et al. A DTI study on the corpus callosum of treatment-naïve boys with ‘pure’ Tourette syndrome. Psychiatry Res Neuroimaging. (2016) 247:1–8. doi: 10.1016/j.pscychresns.2015.12.003
23. Zarei M, Mataix-Cols D, Heyman I, Hough M, Doherty J, Burge L, et al. Changes in gray matter volume and white matter microstructure in adolescents with obsessive-compulsive disorder. Biol Psychiatry. (2011) 70:1083–90. doi: 10.1016/j.biopsych.2011.06.032
24. Leckman JF, Riddle MA, Hardin MT, Ort SI, Swartz KL, Stevenson J, et al. The yale global tic severity scale: initial testing of a clinician-rated scale of tic severity. J Am Acad Child Adolesc Psychiatry. (1989) 28:566–73. doi: 10.1097/00004583-198907000-00015
25. Goodman WK, Price LH, Rasmussen SA, Mazure C, Fleischmann RL, Hill CL, et al. The Yale-Brown Obsessive Compulsive Scale. I development, use, and reliability. Arch Gen Psychiatry. (1989) 46:1006–11. doi: 10.1001/archpsyc.1989.01810110048007
26. Bryden MP. Measuring handedness with questionnaires. Neuropsychologia. (1977) 15:617–24. doi: 10.1016/0028-3932(77)90067-7
27. Kaufman J, Birmaher B, Brent D, Rao U, Flynn C, Moreci P, et al. Schedule for affective disorders and schizophrenia for school-age children-present and lifetime version (K-SADS-PL): initial reliability and validity data. J Am Acad Child Adolesc Psychiatry. (1997) 36:980–8. doi: 10.1097/00004583-199707000-00021
28. Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, Mackay CE, et al. Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage. (2006) 31:1487–505. doi: 10.1016/j.neuroimage.2006.02.024
29. Mori S, Zhang J. Principles of diffusion tensor imaging and its applications to basic neuroscience research. Neuron. (2006) 51:527–39. doi: 10.1016/j.neuron.2006.08.012
30. Piras F, Piras F, Caltagirone C, Spalletta G. Brain circuitries of obsessive compulsive disorder: a systematic review and meta-analysis of diffusion tensor imaging studies. Neurosci Biobehav Rev. (2013) 37:2856–77. doi: 10.1016/j.neubiorev.2013.10.008
31. Piras F, Piras F, Abe Y, Agarwal SM, Anticevic A, Ameis S, et al. White matter microstructure and its relation to clinical features of obsessive–compulsive disorder: findings from the ENIGMA OCD working group. Transl Psychiatry. (2021) 11:173. doi: 10.1038/s41398-021-01276-z
32. Bruce AB, Yuan W, Gilbert DL, Horn PS, Jackson HS, Huddleston DA, et al. Altered frontal-mediated inhibition and white matter connectivity in pediatric chronic tic disorders. Exp Brain Res. (2021) 239:955–65. doi: 10.1007/s00221-020-06017-0
33. Plessen KJ, Grüner R, Lundervold A, Hirsch JG, Xu D, Bansal R, et al. Reduced white matter connectivity in the corpus callosum of children with Tourette syndrome. J Child Psychol Psychiatry. (2006) 47:1013–22. doi: 10.1111/j.1469-7610.2006.01639.x
34. Jeppesen SS, Debes NM, Simonsen HJ, Rostrup E, Larsson HB, Skov L. Study of medication-free children with Tourette syndrome do not show imaging abnormalities. Mov Disord. (2014) 29:1212–6. doi: 10.1002/mds.25858
35. Debes N, Jeppesen S, Raghava JM, Groth C, Rostrup E, Skov L. Longitudinal magnetic resonance imaging (MRI) analysis of the developmental changes of tourette syndrome reveal reduced diffusion in the cortico-striato-thalamo-cortical pathways. J Child Neurol. (2015) 30:1315–26. doi: 10.1177/0883073814560629
36. Neuner I, Kupriyanova Y, Stöcker T, Huang R, Posnansky O, Schneider F, et al. White-matter abnormalities in Tourette syndrome extend beyond motor pathways. Neuroimage. (2010) 51:1184–93. doi: 10.1016/j.neuroimage.2010.02.049
37. Thomalla G, Siebner HR, Jonas M, Baumer T, Biermann-Ruben K, Hummel F, et al. Structural changes in the somatosensory system correlate with tic severity in Gilles de la Tourette syndrome. Brain. (2009) 132:765–77. doi: 10.1093/brain/awn339
38. Cavanna AE, Stecco A, Rickards H, Servo S, Terazzi E, Peterson B, et al. Corpus callosum abnormalities in Tourette syndrome: an MRI-DTI study of monozygotic twins. J Neurol Neurosurg Psychiatry. (2010) 81:533–5. doi: 10.1136/jnnp.2009.173666
39. Draganski B, Martino D, Cavanna AE, Hutton C, Orth M, Robertson MM, et al. Multispectral brain morphometry in Tourette syndrome persisting into adulthood. Brain. (2010) 133:3661–75. doi: 10.1093/brain/awq300
40. Lebel C, Deoni S. The development of brain white matter microstructure. Neuroimage. (2018) 182:207–18. doi: 10.1016/j.neuroimage.2017.12.097
41. Groth C, Mol Debes N, Rask CU, Lange T, Skov L. Course of tourette syndrome and comorbidities in a large prospective clinical study. J Am Acad Child Adolesc Psychiatry. (2017) 56:304–12. doi: 10.1016/j.jaac.2017.01.010
42. Ciechanski P, Zewdie E, Kirton A. Developmental profile of motor cortex transcallosal inhibition in children and adolescents. J Neurophysiol. (2017) 118:140–8. doi: 10.1152/jn.00076.2017
43. Sisti HM, Geurts M, Gooijers J, Heitger MH, Caeyenberghs K, Beets IAM, et al. Microstructural organization of corpus callosum projections to prefrontal cortex predicts bimanual motor learning. Learn Mem. (2012) 19:351–7. doi: 10.1101/lm.026534.112
44. Takeuchi N, Oouchida Y, Izumi SI. Motor control and neural plasticity through interhemispheric interactions. Neural Plast. (2012) 2012:823285. doi: 10.1155/2012/823285
45. Albin RL, Mink JW. Recent advances in Tourette syndrome research. Trends Neurosci. (2006) 29:175–82. doi: 10.1016/j.tins.2006.01.001
46. Mink JW. Neurobiology of basal ganglia and Tourette syndrome: basal ganglia circuits and thalamocortical outputs. Adv Neurol. (2006) 99:89–98.
47. Ashtari M. Anatomy and functional role of the inferior longitudinal fasciculus: a search that has just begun. Dev Med Child Neurol. (2012) 54:6–7. doi: 10.1111/j.1469-8749.2011.04122.x
48. Ortibus E, Verhoeven J, Sunaert S, Casteels I, De Cock P, Lagae L. Integrity of the inferior longitudinal fasciculus and impaired object recognition in children: a diffusion tensor imaging study. Dev Med Child Neurol. (2012) 54:38–43. doi: 10.1111/j.1469-8749.2011.04147.x
49. Latini F, Mårtensson J, Larsson EM, Fredrikson M, Åhs F, Hjortberg M, et al. Segmentation of the inferior longitudinal fasciculus in the human brain: a white matter dissection and diffusion tensor tractography study. Brain Res. (2017) 1675:102–15. doi: 10.1016/j.brainres.2017.09.005
50. Zemmoura I, Burkhardt E, Herbet G. The inferior longitudinal fasciculus: anatomy, function and surgical considerations. J Neurosurg Sci. (2021) 65:590–604. doi: 10.23736/S0390-5616.21.05391-1
51. Kim S, Jackson SR, Groom M, Jackson GM. Visuomotor learning and unlearning in children and adolescents with tourette syndrome. Cortex. (2018) 109:50–9. doi: 10.1016/j.cortex.2018.08.007
52. Takács Á, Kóbor A, Chezan J, Élteto N, Tárnok Z, Nemeth D, et al. Is procedural memory enhanced in Tourette syndrome? evidence from a sequence learning task. Cortex. (2018) 100:84–94. doi: 10.1016/j.cortex.2017.08.037
53. Kleimaker M, Takacs A, Conte G, Onken R, Verrel J, Bäumer T, et al. Increased perception-action binding in Tourette syndrome. Brain. (2020) 143:1934–45. doi: 10.1093/brain/awaa111
54. Hu X, Zhang L, Bu X, Li H, Gao Y, Lu L, et al. White matter disruption in obsessive-compulsive disorder revealed by meta-analysis of tract-based spatial statistics. Depress Anxiety. (2020) 37:620–31. doi: 10.1002/da.23008
55. Li Q, Zhao Y, Huang Z, Guo Y, Long J, Luo L, et al. Microstructural white matter abnormalities in pediatric and adult obsessive-compulsive disorder: a systematic review and meta-analysis. Brain Behav. (2020) 11:e01975. doi: 10.1002/brb3.1975
56. Fitzgerald KD, Liu Y, Reamer EN, Taylor SF, Welsh RC. Atypical frontal-striatal-thalamic circuit white matter development in pediatric obsessive-compulsive disorder. J Am Acad Child Adolesc Psychiatry. (2014) 53:1225–33, 1233.e1–9. doi: 10.1016/j.jaac.2014.08.010
57. Lochner C, Fouché JP, du Plessis S, Spottiswoode B, Seedat S, Fineberg N, et al. Evidence for fractional anisotropy and mean diffusivity white matter abnormalities in the internal capsule and cingulum in patients with obsessive–compulsive disorder. J Psychiatry Neurosci. (2012) 37:193–9. doi: 10.1503/jpn.110059
58. Menzies L, Chamberlain SR, Laird AR, Thelen SM, Sahakian BJ, Bullmore ET. Integrating evidence from neuroimaging and neuropsychological studies of obsessive-compulsive disorder: the orbitofronto-striatal model revisited. Neurosci Biobehav Rev. (2008) 32:525–49. doi: 10.1016/j.neubiorev.2007.09.005
59. Posner J, Marsh R, Maia TV, Peterson BS, Gruber A, Simpson HB. Reduced functional connectivity within the limbic cortico-striato-thalamo-cortical loop in unmedicated adults with obsessive-compulsive disorder. Hum Brain Mapp. (2014) 35:2852–60. doi: 10.1002/hbm.22371
60. Robbins TW, Vaghi MM, Banca P. Obsessive-compulsive disorder: puzzles and prospects. Neuron. (2019) 102:27–47. doi: 10.1016/j.neuron.2019.01.046
61. Shin NY, Lee TY, Kim E, Kwon JS. Cognitive functioning in obsessive-compulsive disorder: a meta-analysis. Psychol Med. (2014) 44:1121–30. doi: 10.1017/S0033291713001803
62. Bloch MH, Sukhodolsky DG, Dombrowski PA, Panza KE, Craiglow BG, Landeros-Weisenberger A, et al. Poor fine-motor and visuospatial skills predict persistence of pediatric-onset obsessive-compulsive disorder into adulthood. J Child Psychol Psychiatry. (2011) 52:974–83. doi: 10.1111/j.1469-7610.2010.02366.x
63. Alhamud A, Tisdall MD, Hess AT, Hasan KM, Meintjes EM, van der Kouwe AJW. Volumetric navigators for real-time motion correction in diffusion tensor imaging. Magn Reson Med. (2012) 68:1097–108. doi: 10.1002/mrm.23314
64. Johansen-Berg H, Rushworth MFS. Using diffusion imaging to study human connectional anatomy. Annu Rev Neurosci. (2009) 32:75–94. doi: 10.1146/annurev.neuro.051508.135735
65. Toga AW, Thompson PM, Sowell ER. Mapping brain maturation. Trends Neurosci. (2006) 29:148–59. doi: 10.1016/j.tins.2006.01.007
66. Beaulieu C. The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed. (2002) 15:435–55. doi: 10.1002/nbm.782
67. Mori S, Wakana S, Zijl PCM van, Nagae-Poetscher LM. MRI Atlas of Human White Matter. Elsevier (2005).p. 249.
68. Alexander AL, Lee JE, Lazar M, Field AS. Diffusion tensor imaging of the brain. Neurotherapeutics. (2007) 4:316–29. doi: 10.1016/j.nurt.2007.05.011
69. Song SK, Sun SW, Ramsbottom MJ, Chang C, Russell J, Cross AH. Dysmyelination revealed through MRI as increased radial (but unchanged axial) diffusion of water. Neuroimage. (2002) 17:1429–36. doi: 10.1006/nimg.2002.1267
Keywords: Tourette syndrome, obsessive-compulsive disorder (OCD), white matter (WM) microstructural organization, drug-naïve, diffusion tensor imaging (DTI), fractional anisotropy (FA)
Citation: Bharti K, Conte G, Tommasin S, Giannì C, Suppa A, Mirabella G, Cardona F and Pantano P (2022) White matter alterations in drug-naïve children with Tourette syndrome and obsessive-compulsive disorder. Front. Neurol. 13:960979. doi: 10.3389/fneur.2022.960979
Received: 03 June 2022; Accepted: 29 August 2022;
Published: 03 October 2022.
Edited by:
Huifang Shang, Sichuan University, ChinaReviewed by:
Joelle Sarlls, National Institutes of Health (NIH), United StatesLuca Marsili, University of Cincinnati, United States
Copyright © 2022 Bharti, Conte, Tommasin, Giannì, Suppa, Mirabella, Cardona and Pantano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Giulia Conte, giulia.conte@uniroma1.it
†These authors have contributed equally to this work and share first authorship