1971

Altered transient connectivity states at rest in Tourette’s syndrome
Shukti Ramkiran1,2,3, Ravichandran Rajkumar1,2,3, N. Jon Shah1,3,4,5, and Irene Neuner1,2,3
1Institute of Neuroscience and Medicine - 4 (INM-4), Forschungszentrum Juelich, Juelich, Germany, 2Department of Psychiatry, Psychotherapy and Psychosomatics, RWTH Aachen University, Aachen, Germany, 3JARA – BRAIN – Translational Medicine, Juelich, Germany, 4Department of Neurology, RWTH Aachen University, Aachen, Germany, 5Institute of Neuroscience and Medicine - 11 (INM-11), JARA, Forschungszentrum Juelich, Juelich, Germany

Synopsis

Analysis of dynamic functional connectivity enables identification of spatio-temporal variations in functional connectivity. This is useful in identifying transient states of abnormalities in various brain disorders. Tourette syndrome (TS) is a neurodevelopmental disorder characterized by multiple motor and at least one vocal or phonic tic. We identified eight distinct states in TS patients and healthy controls (HC) using spectral clustering. Our results show preferential state dominance and increased inter-state variability in TS patients as compared to healthy controls.

Introduction

Dynamic functional connectivity (DFC) refers to the analysis of functional connectivity changes over short periods of time using a sliding-window approach. It offers the advantage over classical static functional connectivity of capturing temporal variations in spatial connectivity. This is particularly interesting as the mind does not remain constant and switches between mental-states at rest1. Several studies have shown this tool to be useful in identifying transient states of mindfulness and mind-wandering1, aberrant transient states in schizophrenia2 and traumatic brain injury3. Tourette syndrome (TS) is a neurodevelopmental disorder characterized by multiple motor and at least one vocal or phonic tic4, often accompanied by other comorbidities5. Several studies have shown abnormal patterns of static connectivity in TS5–7 and this study aims to investigate abnormalities in patterns of transient connectivity.

Methods

Resting-state fMRI data from 28 TS patients and 28 age-matched HC, recorded in a 1.5T MR system (Sonata, Siemens, Germany) were used for this study. Details of the participants, the acquisition protocol, the structural and functional preprocessing and denoising can be found in5. For performing dynamic functional connectivity analyses, the entire scan session was decomposed into 24 sliding windows of length 100s and overlap of 75s. Each window contained 35 fMRI volumes. The sliding window decomposition was followed by 1st level (subject-level) analysis of functional connectivity in each window. This was done by specifying a within-subject ROI-ROI bivariate correlation model (132 whole-brain ROIs: Harvard-Oxford maximum likelihood cortical and subcortical atlas, AAL cerebellar parcellations), using CONN18b8 . The connectivity matrices were imported into MATLAB and checked for data inconsistencies. Subjects with missing connections in any given time window were excluded from further analyses. This was done as the un-thresholded connectivity analyses should yield fully connected matrices. The reason for these missing connections could be the noisy data (scrubbing or improper segmentation in a particular ROI leading to 0 voxels in turn leading to missing connections), but this needs to be investigated further. Visual inspection showed that the structural MRI data of some of these subjects were corrupted by motion artefacts. After exclusion, the data consisted of 23 TS patients and 22 HC. Average connectivity matrices were computed for each time window separately for TS patients and HC. Spectral clustering was applied to the time windows to identify distinct states from the 24 windows. Using the Calinski Harabasz criteria, the optimal number of clusters was identified as 8 for both TS patients and HC. Thus, 8 distinct functional connectivity states were obtained separately in TS patients and HC. The dwell time or the time spent in each state (i.e the size of the cluster or the number of windows belonging to each state) was calculated.

Results

The 8 distinct states of functional connectivity obtained, along with their dwell times, are given in figures 1 and 2 for TS patients and HC respectively. While the distribution of states is more uniform in HC (most of the states have DT = 3, mean DT = 3, SD = 0.53), the states in TS patients range from most dominant (state 4, DT = 5) to least dominant (state 8, DT = 1) (mean DT = 3, SD = 1.41). Similarly, the spatial patterns of connectivity appear to be more uniform in HC and more varying in TS patients. Through visual inspection, state 7 of the TS patients resembles the spatial connectivity pattern of HC, while state 1 reflects more positive connectivity and states 2 and 5 reflect more negative connectivity.

Discussion

As we can see from figures 1 and 2, connectivity patterns in the healthy population remain more stable and uniform over time compared to TS patients. TS patients show a preferential state dominance and the most dominant state being state 4, which, based on visual inspection, looks spatially different from the states in the healthy population. On the other hand, the second most dominant state (state 7) closely resembles the states seen in the healthy population. This inter-state variability and transient switching to states closely resembling those of the healthy controls could mean that abnormalities in communication occur in some states, while healthy functioning resumes in other states. Several studies have shown network-level abnormalities, defects in network maturation and alterations in information exchange mechanisms in static connectivity5,7. The transient state changes observed in this preliminary study pave the way to exploring network level and individual subject level changes in transient connectivity.

Conclusion

DFC allows the capturing of spatio-temporal variations in functional connectivity. Application of DFC on TS patients revealed preferential state dominance and increased variability in transient states as compared to healthy volunteers.

Acknowledgements

1. The German Tourette's Association for their support and travel funds.

2. All participating Tourette's patients and their families for their cooperation

3. Ms. Claire Rick for editorial input.

References

1. Mooneyham, B. W. et al. States of mind: Characterizing the neural bases of focus and mind-wandering through dynamic functional connectivity. J. Cogn. Neurosci. 29, 495–506 (2017).

2. Damaraju, E. et al. Dynamic functional connectivity analysis reveals transient states of dysconnectivity in schizophrenia. NeuroImage Clin. 5, 298–308 (2014).

3. Vergara, V. M., Mayer, A. R., Kiehl, K. A. & Calhoun, V. D. Dynamic functional network connectivity discriminates mild traumatic brain injury through machine learning. NeuroImage Clin. 19, 30–37 (2018).

4. Robertson, M. M. et al. Gilles de la Tourette syndrome. Nat. Rev. Dis. Prim. 3, 16097 (2017).

5. Ramkiran, S., Heidemeyer, L., Gaebler, A., Shah, N. J. & Neuner, I. Alterations in basal ganglia-cerebello-thalamo-cortical connectivity and whole brain functional network topology in Tourette’s syndrome. NeuroImage Clin. 24, (2019).

6. Wen, H. et al. Combining Disrupted and Discriminative Topological Properties of Functional Connectivity Networks as Neuroimaging Biomarkers for Accurate Diagnosis of Early Tourette Syndrome Children. Molecular Neurobiology 1–19 (2017). doi:10.1007/s12035-017-0519-1

7. Worbe, Y. et al. Functional immaturity of cortico-basal ganglia networks in Gilles de la Tourette syndrome. Brain 135, 1937–1946 (2012).

8. Whitfield-Gabrieli, S. & Nieto-Castanon, A. Conn: A Functional Connectivity Toolbox for Correlated and Anticorrelated Brain Networks. Brain Connect. 2, 125–141 (2012).

Figures

Figure 1: DFC states in TS patients – states 1 to 8 in the order left to right and top to bottom. The rows and columns indicate the different ROIs and the colour bar indicates the connectivity value. The graph on the bottom right shows the dwell time of each of the 8 states.

Figure 2: DFC states in healthy controls – states 1 to 8 in the order left to right and top to bottom. The rows and columns indicate the different ROIs and the colour bar indicates the connectivity value. The graph in the bottom right shows the dwell time of each of the 8 states.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
1971