0796

Aberrant Functional Brain Network Topology for Classification between Major Depressive Disorder and Healthy Controls
Yael Jacob1, Laurel S Morris1, Kuang-Han Huang1, Molly Schneider1, Gaurav Verma1, James W Murrough1, and Priti Balchandani1
1Icahn School of Medicine at Mount Sinai, New York, NY, United States

Synopsis

Currently, diagnosis for major depressive disorder (MDD) is largely reliant on self-reported symptoms. The ability to identify MDD without self-report is greatly needed. Implementing a graph-theoretical analysis on resting state fMRI (rsfMRI), we tested whether whole-brain network topology can be used as predictors of MDD using a machine learning algorithm. We found that MDD patients exhibit aberrant network centrality measures within the right hippocampus, supramarginal and parsopercularis. Using these as predictors in a machine learning algorithm we were able to classify MDD and controls with total accuracy of 81%, demonstrating the applicability of rsfMRI for diagnostics of MDD.

Introduction

Major depressive disorder (MDD) is one of the world’s largest health problems1. Currently, the diagnosis entirely depends on clinical symptoms, whereas the underlying brain pathology remains largely unclear. Recently, the mathematical field of graph theory emerged as a tool for characterizing brain network features that can distinguish between healthy and pathological states2-4. The functional brain network graph is composed of nodes, representing regions, and edges, representing connections. Based on graph theory analyses, various studies have reported aberrant topological organization among MDD patients5. Yet, it is currently unclear how these whole-brain network features can be used with data-driven machine learning tools for the classification of MDD. Implementing a graph theoretical analysis based on ultra-high field 7-Tesla functional MRI data, we tested whether whole brain network connectivity hierarchies during resting state can distinguish between MDD patients and controls, and if these can be used as predictors in a machine learning algorithm.

Methods

21 MDD patients and 21 controls underwent resting state fMRI. Data were acquired on a Siemens Magnetom 7T MRI scanner (Erlangen, Germany). Functional images were processed using the multi-echo independent component analysis implemented in the AFNI meica.py toolbox6. Each subjects anatomical T1 weighted brain image was segmented into the Desikan-Killiany Atlas7 applied in FreeSurfer. The segmentation resulted in 84 regions of interest each representing a node of the network (Figure 1A). To define the network edges, we calculated the partial correlation coefficients between the regional mean time-series of all pairwise regions (excluding the effects of all other regions) (Figure 1C). To enable comparison across participants, we used a sparsity threshold S, which retains S% of the top connections for each participant (Figure 1D). Using the Brain Connectivity Toolbox8, we examined common local nodal centrality features; 1) node strength, which is the sum of weights of connections to the node (Figure 1E); 2) betweenness centrality, which is a measure of the number of shortest paths that traverse a given node; and 3) local efficiency, which is the efficiency computed on the neighbors of a given node. We examined the network features across a range of thresholds (10% < S < 30% in steps of 1%)9. We then calculated the area under the curve for each network feature, which provides a summarized measure independent of single threshold selection9 (Figure 1F). Then we conducted a between-group t-test for each region's local feature. All results were corrected for multiple comparisons (number of nodes), using false discovery rate (FDR)10. Next, we conducted Pearson correlations to assess the association between the significant results and disease duration and severity among the MDD patients. Partial correlation was used to control for age, gender and region volume as covariates. Quadratic discriminant analysis (QDA)11 was used for between group classification using Matlab Classification Learner application (The MathWorks Inc., Natick, MA). All the statistically significant network features were used as predictors for the machine learning classification. We then assessed the generalization ability through fivefold cross-validation. For statistical significance 1,000 repetitions were performed and then for permutation tests another 1000 repetitions of shuffled datasets. The overall accuracy, sensitivity, specificity, and positive predictive value (PPV) were computed for each classifier.

Results

The results revealed that compared to controls, MDD patients exhibited both increased and decreased centrality measures. Specifically, there was increased betweenness centrality within the right hippocampus (t=4.45, p<0.00008) and right supramarginal (t=2.99, p<0.001), which also presented increased node strength (t=3.39, p<0.001). There was decreased node strength of the right parsopercularis (t=3.25, p<0.001). All were qFDR<0.05 corrected (Figure 2). The right hippocampus also showed a significant association with the duration of the current episode showing, whereby the higher the betweenness centrality, the longer is the depressive episode (r=0.52, p<0.04) (Figure 3). The machine learning QDA algorithm using the significant network feature values (right parsopercularis strength, right supramarginal strength and betweenness, and right hippocampus betweenness, total of 4 features) revealed significant classification power (Figure 4). The overall accuracy was 81% (0.81 ± 0.037, p<0.001), with sensitivity of 86% (0.86 ± 0.047, p<0.001) and specificity of 76% (0.76 ± 0.053, p<0.02). The PPV for the MDD were 84% (0.84 ± 0.045, p<0.001) and for the controls, 78% (0.78 ± 0.039, p<0.01).

Discussion

Our results showed that the MDD patients exhibit aberrant network topology as compared to controls during resting state. We also showed that the longer the duration of the depressive episode, the higher the patients’ right hippocampus betweenness centrality. Increased betweenness centrality indicates the enhanced rapid information flow through the hippocampus, which might imply of a maladaptive functional behavior of the network. These results replicate previous studies5,12. Furthermore, using these results as predictors in a machine learning algorithm we were able to classify MDD versus controls with total accuracy of 81%. This result demonstrates the applicability of rsfMRI for the diagnostics of MDD. Further research is needed to replicate this finding and improve the generalization and predictive capacity for robust clinical translation.

Acknowledgements

Funding was provided by NIH R01 MH109544. Additional support was provided by the Icahn School of Medicine Capital Campaign, Translational and Molecular Imaging Institute and Department of Radiology, Icahn School of Medicine at Mount Sinai and Siemens Healthcare.

References

1. Collins, Pamela Y., et al. "Grand challenges in global mental health." Nature 475.7354 (2011): 27.

2. Bassett, Danielle S., et al. "Hierarchical organization of human cortical networks in health and schizophrenia." Journal of Neuroscience 28.37 (2008): 9239-9248.

3. Bullmore, Ed, and Olaf Sporns. "Complex brain networks: graph theoretical analysis of structural and functional systems." Nature reviews neuroscience 10.3 (2009): 186.

4. Sporns, Olaf. Networks of the Brain. MIT press, 2010.

5. Gong, Qiyong, and Yong He. "Depression, neuroimaging and connectomics: a selective overview." Biological psychiatry 77.3 (2015): 223-235.

6. Kundu, Prantik, et al. "Differentiating BOLD and non-BOLD signals in fMRI time series using multi-echo EPI." Neuroimage 60.3 (2012): 1759-1770.

7. Desikan, Rahul S., et al. "An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest." Neuroimage 31.3 (2006): 968-980.

8. Rubinov, Mikail, and Olaf Sporns. "Complex network measures of brain connectivity: uses and interpretations." Neuroimage 52.3 (2010): 1059-1069.

9. Korgaonkar, Mayuresh S., et al. "Abnormal structural networks characterize major depressive disorder: a connectome analysis." Biological psychiatry 76.7 (2014): 567-574.

10. Benjamini, Yoav, and Yosef Hochberg. "Controlling the false discovery rate: a practical and powerful approach to multiple testing." Journal of the Royal statistical society: series B (Methodological) 57.1 (1995): 289-300.

11. Zhang, Michael Q. "Identification of protein coding regions in the human genome by quadratic discriminant analysis." Proceedings of the National Academy of Sciences 94.2 (1997): 565-568.

12. Zhang, Junran, et al. "Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder." Biological psychiatry 70.4 (2011): 334-342.

Figures

Functional network graph theory analysis. (A) Each subjects anatomical brain image was segmented into 84 regions. (B) We then extracted the averaged fMRI time-course per region. (C) To define network edges, we calculated the partial correlation coefficients between all pairwise regions. (D) A threshold of the top connections for each network was applied across a range from 0.1-0.3. (E) We then examined local network properties of node strength, betweenness centrality, and local efficiency.

Local network properties differentiating between MDD and HC. Applying graph theory analysis we found several regions’ centrality measures that are significantly different between the groups. (A) Increased betweenness centrality property of the right hippocampus and the right supramarginal among the MDD patients. Suggesting that these regions become more critical for the information flow among MDD. (B) Increased strength of the right supramarginal indicating its increased connectivity and reduced connectivity of the parsopercularis among MDD.

Graph theory measures association with duration of current episode. Additional exploration revealed that increased right hippocampus betweenness centrality measure was correlated to greater duration of current episode. *p<0.05

Classification between MDD and HC using machine learning algorithm. Quadratic discriminant analysis algorithm was performed on the graph theory results. Permutation tests were conducted 1,000 times for the overall accuracy (A), sensitivity (B), specificity (C) and positive predictive value (PPV) of HC (D) and MDD (E). The dashed red lines represent the accuracy of the original dataset. The bar graph demonstrates the averaged classification features calculated on 1,000 repetitions (F). This analysis revealed significant classification power with an overall accuracy of 81%.

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