Data-driven functional sub-division of the sensory-motor network using hierarchical clustering for resting-state fMRI data.
Yanlu Wang1 and Tie-Qiang Li1,2

1Clinical Sciences, Intervention and Technology, Karolinska Institute, Stockholm, Sweden, 2Medical Physics, Karolinska University Hospital, Stockholm, Sweden

Synopsis

A data-driven analysis method based on hierarchical clustering was used to analyze the sensory-motor resting-state network from resting-state fMRI data. It was used to analyze the network’s functional sub-division, and intra-network functional organization, in different levels of detail. Sub-network for the sensory-motor network as obtained by hierarchical clustering is anatomically and functionally sensible. Further sub-division of the paracentral lobule network hub successfully revealed its functional sub-division in great detail. The intra-network organization of intrinsic functional connectivity derived from spontaneous activity of the brain at rest reflects consistently, the functional and neural anatomic connectivity topography of the resting-state network.

Introduction

We have previously shown that hierarchical clustering can also be used to extract functional connectivity networks, and that it can be used to investigate the intra-network organization of resting-state networks through stratifying data into a hierarchical structure naturally. In this study, we devise a data-driven scheme to analyze resting-state networks using hierarchical clustering and applied this scheme to analyze the sensory-motor network.

Aim

The purpose of this study is to apply the data-driven hierarchical clustering analysis scheme to the sensory-motor network to analyze its functional sub-division, and intra-network functional organization.

Materials and Methods

The sensory-motor was previously extracted using full-brain hierarchical clustering at voxel level from 86 health subjects1. The dendrogram for the network, as obtained though hierarchical clustering, was cut at increasing cut levels (k). At each cut, the new clusters were evaluated for significance (p<0.01). The significance criterion is based on both intra-cluster correlation and cluster size, obtained through group-level statistics from intra-cluster correlation values and Monte-Carlo simulations of cluster formation from Gaussian noise (fig. 1). Clusters that are determined to be not meaningful in terms of functional subdivision were removed from the dendrogram. This process was repeated up to cut level k=80 whereupon enough sub-networks (19) were obtained to warrant figure display. A cluster stands out to be the largest among the resulting sub-networks, and appears to be the central hub of the network2. It was chosen for further investigation to illustrate the frameworks capability to analyze resting-state networks greater detail. The same procedure was carried out as previously. This further sub-division corresponds to cutting the dendrogram at nodes corresponding to a plateau in the inconsistency coefficient, indicating that this is a natural termination point for further sub-divisions.

Results

The intra-network hierarchy of the extracted sub-networks and their neuroanatomical locations are depicted in figure 2. The extracted sub-networks include 4 groups of bilateral sub-networks: Insula-auditory sub-networks (fig. 2, clusters 1-3); paracentral lobule and cingulate-motor sub-networks (fig. 2, clusters 4 and 5); facial expression control (fig. 2, clusters 11-13); and hand movement control (fig. 2, clusters 6, 7, 9, and 10). There exist also two groups of more loosely connected unilateral sub-networks: The right parietal sub-network group (fig. 2, clusters 15 and 16) and the left insular-STG (superior temporal gyrus) sub-network group (fig. 2, clusters 17-19). This intra-network organization of intrinsic functional connectivity derived from spontaneous activity of the brain at rest reflects consistently the functional and neural anatomic connectivity topography network, which consists of S1, M1, and some pre/post- central gyrus areas divided into dorsal and ventral subgroups in addition to the parietal operculum and the auditory cortex3.

The functional sub-division of the paracentral lobule effectively distinguishes the supplementary motor area (SMA) (fig. 3, cluster 4), areas in the somatosensory cortex responsible for the legs and feet4 (fig. 3, clusters 1-3), and higher order motor-learning areas5 bordering the cingulate cortex (fig.3, cluster 5).

Discussion

The implemented hierarchical clustering scheme is a conscientious method where certain internal characteristics of the data, such as cluster size and intra-cluster correlation. Another characteristic of the data, the inconsistency coefficient, created a natural end-point to the sub-division. Further sub-divisions (up to k=200) yielded no more significant clusters given our cluster evaluation criteria. Hierarchical clustering nevertheless produces a large amount of information that can be difficult to survey even under the current scheme. This visualization problem is overcome by presenting only manageable sized sub-sets at a time. To demonstrate this, we further sub-divided the cluster corresponding to the paracentral lobule as an illustrative example. While our choice is in some ways natural due to its comparatively larger size, and central role in functional modulation within the network, it should be noted that this can be done equally well for any other sub-networks

Conclusion

We have shown that hierarchical clustering is capable of analyzing the functional organization of resting-state networks in different levels of detail. The functional hierarchy constructed through voxel-wise hierarchical clustering corroborate consistently with the known modular organizations from previous clinical and experimental studies.

Acknowledgements

No acknowledgement found.

References

1. Wang Y, Li TQ. Analysis of whole-brain resting-state FMRI data using hierarchical clustering approach. PloS one 2013; 8(10): e76315.

2. Bianchi MT, Danielle S, Matt T, et al. Network Approaches to Diseases of the Brain; 2012.

3. Power JD, Cohen AL, Nelson SM, et al. Functional network organization of the human brain. Neuron 2011; 72(4): 665-78.

4. Fisicaro RA, Jiao RX, Stathopoulos C, Petrovich Brennan NM, Peck KK, Holodny AI. Challenges in Identifying the Foot Motor Region in Patients with Brain Tumor on Routine MRI: Advantages of fMRI. AJNR American journal of neuroradiology 2015; 36(8): 1488-93.

5. Rubia K, Overmeyer S, Taylor E, et al. Hypofrontality in attention deficit hyperactivity disorder during higher-order motor control: a study with functional MRI. The American journal of psychiatry 1999; 156(6): 891-6.

Figures

Figure 1: Plot of the cluster evaluation criterion for significance (p<0.01) based a clusters size and intra-cluster correlation values. Clusters with these parameters that fall below this curve is filtered out. The curve indicates that clusters of a smaller size requires higher intra-cluster correlation values to be considered significant.

Figure 2: The intra-network hierarchical organization of the sub-networks extracted from the sensory-motor network. The bottom panel shows the color-coded sub-networks imposed on the smoothed white matter surface. The top panels shows the functional hierarchy of the color-coded sub-networks extracted (left), and the inconsistency coefficient plot the network of nodes of descending height (right). The vertical line indicates where sub-division was stopped to produce the figure.

Figure 3: The intra-network hierarchical organization of the cluster corresponding to the paracentral lobule. EThe top panel shows the color-coded sub-networks imposed on the smoothed white matter surface. The top panels shows the functional hierarchy for the color-coded sub-networks (left), and the inconsistency coefficient plot the network of nodes of descending height (right). The inconsistency coefficient plot plateau is a natural place to terminate the sub-division.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1666