# Brain Connectivity and Information-Flow Breakdown Revealed by a Minimum Spanning Tree-Based Analysis of MRI Data in Behavioral Variant Frontotemporal Dementia

^{1}Medical and Genomic Statistics Unit, Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy^{2}Neurology Unit, Department of Clinical and Experimental Sciences, Centre for Neurodegenerative Disorders, University of Brescia, Brescia, Italy^{3}Alzheimer’s Research Unit, IRCCS Fatebenefratelli, Brescia, Italy^{4}Neuroradiology Unit, University of Brescia, Brescia, Italy

Brain functional disruption and cognitive shortfalls as consequences of neurodegeneration are among the most investigated aspects in current clinical research. Traditionally, specific anatomical and behavioral traits have been associated with neurodegeneration, thus directly translatable in clinical terms. However, these qualitative traits, do not account for the extensive information flow breakdown within the functional brain network that deeply affect cognitive skills. Behavioural variant Frontotemporal Dementia (bvFTD) is a neurodegenerative disorder characterized by behavioral and executive functions disturbances. Deviations from the physiological cognitive functioning can be accurately inferred and modeled from functional connectivity alterations. Although the need for unbiased metrics is still an open issue in imaging studies, the graph-theory approach applied to neuroimaging techniques is becoming popular in the study of brain dysfunction. In this work, we assessed the global connectivity and topological alterations among brain regions in bvFTD patients using a minimum spanning tree (MST) based analysis of resting state functional MRI (rs-fMRI) data. Whilst several graph theoretical methods require arbitrary criteria (including the choice of network construction thresholds and weight normalization methods), MST is an unambiguous modeling solution, ensuring accuracy, robustness, and reproducibility. MST networks of 116 regions of interest (ROIs) were built on wavelet correlation matrices, extracted from 41 bvFTD patients and 39 healthy controls (HC). We observed a global fragmentation of the functional network backbone with severe disruption of information-flow highways. Frontotemporal areas were less compact, more isolated, and concentrated in less integrated structures, respect to healthy subjects. Our results reflected such complex breakdown of the frontal and temporal areas at both intra-regional and long-range connections. Our findings highlighted that MST, in conjunction with rs-fMRI data, was an effective method for quantifying and detecting functional brain network impairments, leading to characteristic bvFTD cognitive, social, and executive functions disorders.

## Introduction

The bvFTD is clinically defined by personality changes and behavioral disturbances, impairment of executive functions and emotional blunting (Gorno-Tempini et al., 2011; Rascovsky et al., 2011). Abnormal intracellular accumulation of either tau or TDP-43 protein is found in most cases (Mann and Snowden, 2017).

Recently, the increasing interest in unraveling functional and structural features of the brain, has benefitted from complex network analyses, such as graph theory, a multidisciplinary approach that allows to analyse complex systems in a straightforward computable way and to describe cerebral areas as nodes, and their connections as edges (Rubinov and Sporns, 2010). Both structural (anatomical) and functional (statistical relationship between two nodes) connectivity can be assessed (Zhang et al., 2016).

Applying graph theoretical methods to neuroimaging techniques is becoming popular in the study of brain dysfunction (Bullmore and Sporns, 2009; Drakesmith et al., 2015; Chiang et al., 2016). Recently, a few studies have considered graph theory analysis applied to rs-fMRI data in patients with bvFTD and have better described FTD-related brain changes (Agosta et al., 2013; Filippi et al., 2017). However, although conventional graph theoretical analyses are helpful in dissecting disease mechanisms (Bullmore and Sporns, 2012), the methodology is significantly hampered by a number of arbitrary choices. Descriptive metrics and their normalization, network type (weighted or unweighted networks), threshold value (fixed cut-off, fixed average degree, fixed edge density, or variable threshold over a range of values) are some of the critical points making network results difficult to reproduce (van Wijk et al., 2010; Telesford et al., 2011; Stam et al., 2014; Drakesmith et al., 2015; Yu et al., 2016). In addition, several network metrics and node centrality indices may assume different importance at either local or global scale (Telesford et al., 2011; Antonenko et al., 2018), whether the graph model accounts for time variant (i.e., dynamic) or invariant (i.e., static) connectivity (Rashid et al., 2014; Park et al., 2018), and the parcellation type, according to Independent Component Analysis (ICA; McKeown et al., 2003; Griffanti et al., 2014) and specific atlases (see Materials and Methods, for data pre-processing in this work).

Minimum Spanning Tree, a unique acyclic subgraph that connects *N* nodes with (*N*-1) edges, and maximizing synchronization between brain areas (i.e., minimizes edge connections), is a promising unambiguous solution to describe complex brain networks (Stam et al., 2014). The use of MST avoids methodological confounding thanks to an efficient integration of topological properties and functional connectivity information (Tewarie et al., 2014; van Diessen et al., 2016; van Lutterveld et al., 2017), ensuring network robustness and reproducibility respect to classical graph analytical approaches (Otte et al., 2015; Tewarie et al., 2015). MST is a tree which has the minimum total edge weight of all possible spanning trees of the original graph. If the brain network can be interpreted as a kind of transport network, an MST might represent the critical backbone of information flow in weighted networks (i.e., contains with high probability all the shortest paths in the network). Given MSTs efficiency and high sensitivity to small fluctuations of connection weights (Van Mieghem and Magdalena, 2005), they can intrinsically provide an accurate representation of subtle and critical topological perturbations, at local scale. On the other hand, MST sparseness could raise issues at global scale. However, it has been demonstrated that if the MST weight distribution is consistent with a power law with sufficiently small exponent value, the global information flow of the underlying network follows entirely MST paths (Van Mieghem and Magdalena, 2005; Van Mieghem and van Langen, 2005; Meier et al., 2015). Moreover, TOMs can be applied to the original adjacency matrix to modulate neighborhood characteristics in MST nodes (Ravasz et al., 2002). In addition, we used wavelet decomposition and correlation to obtain noise-free and robust functional relationships between brain areas (Achard et al., 2006; Zhang et al., 2016).

Collectively, these aspects enable data-driven network comparison of healthy and diseased groups, without normalization or standardization steps, as recently illustrated in EEG and MEG data (Stam et al., 2014; van Dellen et al., 2014; Numan et al., 2017). Furthermore, network robustness and reproducibility should ensure univocal results, minimizing room for ambiguous interpretations. Besides graph theoretical aspects, data acquisition and pre-processing issues may affect results, including brain parcellation and data acquisition technologies. These aspects should follow the principles of common usage, availability, cost effectiveness, and non-invasivity, that can secondarily affect methodological choices and issues (Hohenfeld et al., 2018). Despite objective difficulties in generating a consensus functional brain map, especially for rare disorders, a set of reference resting state functional networks have been replicated in many different studies (Hohenfeld et al., 2018). According to well-established rs-fMRI literature, three reference networks are mainly involved in bvFTD functional breakdown: the default-mode network (DMN), the salience network (SN), and the executive network (EN), accounting for the cognitive, emotional, and social impairments characterizing this pathology (Raichle, 2015; Trojsi et al., 2015; Sedeño et al., 2016; Hohenfeld et al., 2018).

In this study, we leverage MST-based analysis of rs-fMRI data to investigate large-scale functional network alterations, inspecting global and local network properties in bvFTD patients, compared to a group of HC, and providing a novel MST procedure, combining individual tree-based global evidences and two-group topological aspects.

## Materials and Methods

The analysis workflow of the rs-fMRI data is illustrated in Figure 1 and described in detail below.

**Figure 1.** General study workflow. The most important steps of connectome extraction. MRI signals acquisition and brain parcellation represent the first phases of resting state functional magnetic resonance imaging data pre-processing. Wavelet correlation matrix and minimum spanning tree (MST) network calculation are included in the analysis phase. Wavelet transformation was applied to the average voxel time series mapped to AAL 116 brain regions of interest (ROI), allowing to obtain the statistical relationships between nodes (i.e., wavelet correlation matrices). Each matrix is the input of the MST algorithm, from which parameters and topological features (including edges partition and nodes cluster) have been calculated.

### Subjects

Forty-one patients with a probable bvFTD diagnosis, according to current criteria (Gorno-Tempini et al., 2011), were recruited at the Center for Neurodegenerative Disorders, University of Brescia, Italy All patients underwent an extensive neuropsychological assessment, as previously published (Gazzina et al., 2016), genetic screening for the most frequent monogenic causes of FTD (i.e., *Granulin, C9orf72*, and *Microtubule Associated Protein* Tau; Cosseddu et al., 2018) and brain MRI structural imaging study. In the present study, none of the bvFTD cases carried pathogenic mutations of monogenic bvFTD.

Thirty-nine HC, recruited from voluntary individuals, were used as control group. HC underwent a brief standardized neuropsychological assessment (Mini-Mental State Examination; MMSE >=27). Table 1 shows the demographic information of the participants. The study, in conformity with the Helsinki Declaration, was approved by the Brescia Hospital Ethics Committee. Informed consent was obtained from all participants.

### MRI Acquisition

All imaging was obtained using a 1.5T Siemens Avanto MRI scanner (Siemens, Erlangen, Germany), equipped with a circularly polarized transmit-receive coil. In a single session, the following scans were collected from each subject:

(i) dual-echo TSE [repetition time = 2500 ms, echo time (TE) = 50 ms], to exclude presence of macroscopic brain abnormalities, according to exclusion criteria;

(ii) 3D MPRAGE T1-weighted scan (TR = 2050 ms, TE = 2.56 ms, matrix = 1 × 1 × 1, in-plane field of view [FOV] = 256 × 256 mm^{2}, slice thickness = 1 mm, flip angle = 15°);

(iii) T2^{∗}-weighted EPI sensitized to BOLD contrast [TR = 2500 ms, TE = 50 ms, 29 axial slices parallel to anterior commissure-posterior commissure line (AC-PC) line, matrix = 64 × 64, field of view = 224 mm, slice thickness = 3.5 mm], gap between slices 1.75 mm for rs-fMRI.

Echo planar images were collected during rest for an 8-min period, resulting in a total of 195 volumes.

### Neuroimaging Pre-processing

Functional data were pre-processed using FSL 5.0.8 neuroimaging software, as reported in Jenkinson et al. (2012): (i) the first two volumes were removed to allow signal stabilization; each volume was motion-corrected to a reference volume using MCFLIRT; (ii) non-brain structures were removed using Brain Extraction Tool (BET); (iii) the effect of TR during slice acquisition was reduced using slice-timing correction and the data were spatially smoothed applying the Gaussian kernel with a full width and half maximum (FWHM) of about 7 mm; (iv) grand-mean intensity of the entire data was adjusted by a single multiplicative factor; (v) high-pass temporal filtering Gaussian-weighted least-squares straight line fitting (100 s) was applied; (vi) functional data were co-registered into equivalent native-space T1 weighted image using Boundary-Based Registration (BBR); and (vii) each T1 weighted image was co-registered into standard space template MNI152 using linear (affine with 12 degree of freedom) brain image registration (FLIRT).

To control for motion effect, we included in the further steps only subjects with head motion in a range below or equal to 1 mm (translation) and one degree (rotation). Motion parameters were derived from six degree of freedom registration using MCFLIRT (Jenkinson et al., 2012).

After pre-processing, we applied an automatic approach called “FMRIB’s ICA-based X-noiseifier” (FIX) to detect non-signal components in resting state images which combines the classifiers approach and the Independent Component Analysis (ICA) in a MATLAB environment (Griffanti et al., 2014, 2015). We applied FIX procedure in three steps. First, for each subject we estimate the amount of Gaussian noise of the true dimensionality of the data, i.e., the number of activation and non-Gaussian noise sources using a probabilistic ICA approach implemented in MELODIC (Beckmann and Smith, 2004). Then, we made a random selection of subjects (10 healthy subjects and 10 bvFTD subjects, covering approximately 25% of the whole sample set) to create a subsample to train the FIX’s multi-level classifier. For these subjects, we manually selected the components (white matter, susceptibility artifact, head motion, cardiac pulsation) looking into the thresholded spatial map estimated from single-ICA and the power spectrum of the time series for each component (Griffanti et al., 2014). Lastly, to test if FIX successfully detected the noise components and regressed out the variance (including the six motion parameters derived from MCFLIRT), we looked into a sample of subjects to confirm classification of bad components. The “cleaned” rs-fMRI images (after noise and motion variance regression) of each subject were subsequently used for brain parcellation.

### Brain Parcellation

The ALL atlas (AAL; Tzourio-Mazoyer et al., 2002) was used to parcel brain into 116 (90 cortical and subcortical, and 26 cerebellar) regions of interest (ROIs, Supplementary Table 1). Mean time series were extracted from each ROI by averaging the signal from all voxels within each region, using Marsbar software (^{1}Brett et al., 2002). A subsequent descriptive aggregation was applied to obtain 8 (right and left) macro-ROI (referred to as lobes or “macro-regions”): Frontal, Insular, Limbic, Occipital, Parietal, Subcortical Gray Matter (SCGM; including Thalamus), Temporal, and Cerebellum (including Vermis).

### Wavelet Correlation Analysis

For each subject, the final dataset (Saba et al., 2018) was composed by 193 mean time-series extracted from each brain region. Different correlation estimates define statistical relationship between brain region pairs (Zalesky et al., 2012): in the present study, we used wavelet correlation.

Each temporal series was decomposed using wavelet analysis and characterized by weighted coefficients, proportional to the total amount of energy emitted from the system, relative to a specific scale and brain location. Considering total energy as a frequency-time, wavelet decomposition enables data processing at different hierarchical scale resolutions. Indeed, low-frequency components correspond to coefficients of approximative scale, while high-frequency components correspond to finer scale coefficients (Bullmore et al., 2004). The high flexibility for non-stationary characteristics of the data in each decomposition scale favors wavelet multi-modularity application to fMRI data. Therefore, wavelet correlation results in a higher robustness and noise reduction, with a more homogeneous representation of the original time series and their transformations (Bullmore et al., 2004; Zhang et al., 2016).

Consistent with (Zhang et al., 2016) guidelines, the maximal overlap discrete wavelet transform (MODWT; Achard et al., 2006) with a Daubechies wavelet filter (length equal to 8), was used to band-pass filter on mean time series, and extract wavelet coefficients for the wavelet scales. Given our (TR = 2500 ms, Nyquist frequency = 0.2 Hz), four frequency bands (i.e., scales) were used; scale one: 0.2–0.1 Hz; scale two: 0.05–0.1 Hz; scale three: 0.025–0.05 Hz; and scale four: 0.013–0.025 Hz. Then, we defined a correlation matrix whose *ij*th elements were set by the estimated wavelet correlations between brain regions *i*$\left({\text{\lambda}}_{\text{i}}^{\left(s\right)}\right)$ and region *j*$\left({\text{\lambda}}_{\text{j}}^{\left(s\right)}\right)$:

we focused on wavelet decomposition scale two (*s* = 2), as only this scale reached the significance in three distinct bvFTD/HC connectivity-based Wilcoxon rank sum tests (Supplementary Table 2). Specifically, we quantified bivariate connectivity of the wavelet correlation matrix for each subject (Lynall et al., 2010), with three global measures: (i) *strength*, defined as the average of columns mean; (ii) *diversity*, defined as the average of the columns variance, and (iii) *zero correlation*, defined as the number of correlations with *P* > 0.05, testing the null hypothesis, H_{0} : ρij = 0.

### Minimum Spanning Tree (MST) Structures and Graph Metrics

The MST method overcomes issues concerning arbitrary threshold selection in weighted connected graphs, by joining edges with minimum weight. In other words, it assembles connections minimizing the sum of edge-weights, excluding edges that form a cycle. MST is an extremely efficient binary representation of a full graph *G* (*N, E*) characterized by *N* nodes and *E* = (*N*-1) edges. It yields perfectly comparable networks among different samples, without dependences from vertices and edges number variability, and a best possible synthesis of the original graph information, achieved through the most important subgraph (Stam et al., 2014). Among existing methods for MST search, we applied the Prim’s algorithm (Cormen et al., 2001), to obtain a MST from each subjects’ wavelet correlation matrix.

The MST method handles very different network configurations with their extremes represented by linear or star shapes. In the former, each node has a maximum of two edges (i.e., path-like tree) and two extreme leaves (i.e., nodes with only one link). Whereas, in a star, all nodes are leaves, except the central node to which the other ones are connected (Stam et al., 2014).

Every network can be described through a set of graph metrics (e.g., topological indices) characterized according to its configuration (either linear or star) and its specific graph metrics values. We used global MST measures providing information on graph centrality (maximum degree *k _{max}*, maximum betweenness

*B*

_{max}), distance (diameter

*d*, eccentricity

*Ecc*), association (assortativity

*Ass*), and topological aspects (degree divergence

*K*, leaf fraction

*Lf)*. Their definitions are given in Table 2.

Short distances and overload prevention aspects suggest a good tree configuration or network integration. If a tree has a star-like (i.e., highly connected) topology, it will be characterized by a greater information exchange capacity (i.e., spread of information across the tree), yet a greater probability of central node overload. The opposite behavior is true if the graph has a line-like topology. In particular, an increase of maximum betweenness *B*_{max} and leaf fraction *Lf*, with a decrease of diameter d and eccentricity *Ecc*, tend to have a star-type configuration and a better network integration (van Lutterveld et al., 2017). In addition, maximum betweenness *B*_{max} and degree divergence *K* correlate positively with the presence of some high-degree tree nodes (hub communication) in networks with a scale-free degree distribution (a scale-free network has a large number of nodes with a lower degree and few highly connected hubs; Albert and Barabási, 2002; Mears and Pollard, 2016).

Finally, positive assortativity (*Ass* > 0) indicates that nodes are likely to be connected to other nodes with the same degree, and therefore that the high degree nodes (hubs) tend to be connected to each other (Bullmore and Sporns, 2009). Negative assortativity (*Ass <* 0) is typical of biological networks with hierarchical structure where hubs are connected to nodes with lower degree nodes (Newman, 2003).

### Topological Overlap Measures (TOM) and Shortest Path Tree (SPT)

In the present study, we implemented pre-processing procedures to avoid some of the intrinsic limitations of the MST approach: the absence of triangular connections (i.e., absence of clustering metric) e graph sparseness with a limited number of edges of the resulted network.

First, we applied *topological overlap measures (TOM*; Ravasz et al., 2002; Zhang and Horvath, 2005) to the network adjacency matrix to compensate the absence of triangular connections in the MST. This allowed to analyse MSTs in terms of aggregation or clustering without biases. Specifically:

where l_{ij} = ∑ _{u}a_{iu}a_{ju} is the overlap estimate between two nodes neighborhoods, a_{ij} is the ijth-element of adjacency matrix, k_{i} and k_{j} are the degree measures of ith-node and jth-node, expressed as ki = ∑ _{u≠i}a_{iu}; k_{j} = ∑ _{u≠j}a_{ju}. A higher overlap is associated with a greater relationship between two nodes, and a greater relationship with their common nodes (Mumford et al., 2010). Therefore, TOM modulates neighborhood characteristics of nodes, by quantifying the topological overlap between two nodes against all other nodes in the network. High TOM values identify nodes that constitute a neighborhood (Ravasz et al., 2002).

Secondly, MST is characterized by (*N*-1) edges, resulting a sparse graph with a limited number of edges. Therefore, we evaluate the network performance of MST through the SPT problem (Van Mieghem and Magdalena, 2005; Meier et al., 2015), quantifying if MST is a good representation of the whole graph (i.e., whether the sparse tree may be considered as a critical backbone of original network; Tewarie et al., 2014).

The shortest path is the path with minimum sum of weights from a source to a destination node, and a SPT is the union of the shortest paths from a source node to all other nodes in the graph. An SPT is mainly sensitive to the small, non-negative link weights, around zero. Starting with a full graph *G*, the probability distribution for the link weights of G around zero can be described by a power distribution: F(x) = Pr (X ≤ x) ~ x^{α}, where x ∈ [0,1] represents weights, and the exponent α > 0, defines the *extreme value index* of the probability distribution (Van Mieghem and Magdalena, 2005).

Three specials α-trees correspond to precise α ranges. The α →∞ regime matches a unique weight for all links (i.e., *w* = 1). In the α = 1 regime, link weights result to be uniformly distributed. Finally, in the α → 0 regime, link weights show strong fluctuations. For α → 0, defined as the *stronger disorder condition*, the SPT of the full graph coincide with a MST (Van Mieghem and Magdalena, 2005; Van Mieghem and van Langen, 2005), and the information flow within the network follows only MST links (Meier et al., 2015).

Therefore, we tested the α → 0 regime (i.e., MST) of the weighted full (i.e., original) graph *G* for each subject in three steps:

(1) Define a full graph *G* (*N, E*) by fixing to zero the not statistically significant (*P* > 0.05) wavelet correlations, r_{ij}. The null hypothesis: H_{0} : ρij = 0 was evaluated with the Mutual Information test, $MI=-N\text{\hspace{0.17em}}\text{ln}\left(1-{r}_{\text{ij}}^{2}\right)/2$ (Edwards et al., 2010).

(2) Apply topological overlap measures (i.e., TOM) to the${\text{r}}_{\text{ij}}^{\left(\text{0}\right)}$ = r_{ij} if *P* < 0.05 and ${\text{r}}_{\text{ij}}^{\left(\text{0}\right)}$ = 0 if *P* ≥ 0.05. MST was created by edges weights defined as w_{ij} = min(1/*TOM*_{ij}; 100)/100. This transformation ensures that weights are enclosed in the range (0, 1), and the most important edges (with small weights, i.e., large TOM values) represent the strongest neighborhood connections.

(3) Ranking weights w_{ij} in descending order and estimating the α exponent from the power functionF (x) = c ⋅ x^{α}, by plotting *Y* = log10 [*F*(*x*)] versus *X* = log10 [(x)]. A straight line with an R-squared index approaching 1 is indicative of good fit to the power function, and slope indicates SPT coincidence with MST.

### Edges Partition and Nodes Clustering in MST

Two-group comparisons of topological MST properties were applied to the individual wavelet correlation matrices, between bvFTD and HC, averaged over subjects. Topological structure of nodes and links in an MST have a key role to capture paths with higher importance in information flow (Meier et al., 2015). It is possible to identify an MST subnetwork showing a higher average node (or link) betweenness centrality compared to the rest of MST (Wu et al., 2006). In other terms, this subset of nodes (or links) is used more often than others, and their paths can be considered as a set of *superhighways* (SHW) in MST, i.e., the most important “roads.” Other links in MST constitute “secondary roads.” Identifying SHW enabled to subdivide MST edges (links) in two distinct components with significantly different transport properties.

Based on SHW’s definition, we applied the method suggested by Wu et al. (2006), on both scale-free and Erdos–Renyi (i.e., random) networks, for edges partition in MST. Briefly, considering the fully connected network with the previous TOM-based link weights, we extracted one MST for each case-control group. Through an iterative process, we removed links in descending order of their weights and calculated the degree divergence value (K), which decreases in each cycle with link removals. As demonstrated by Braunstein et al. (2007), the process ends when *K* < 2, and the largest remaining component is the SHW set. Lastly, we measured node betweenness differences (i.e., bvFTD – HC) to assess the amplitude of information flow connectivity gain or loss, and we defined a threshold *b* as the non-zero median of the absolute betweenness differences. Areas showing differences above *b* (or below -*b*), will be taken as markers of functional connectivity gain (or loss).

Clustering nodes in a tree is more complex respect to other graphs. Several types of clustering algorithms have been developed which revealed to be limited and unsuitable to MST characteristics (Yu et al., 2015). Here, we used a novel approach for MST clustering suggested by Yu et al. (2015) based on the geodesic distance matrix D, where the ijth-elements of D represent the number of links of a shortest path between two nodes. As proposed by the authors, we computed vector similarities as the Spearman’s distance, d_{S} = 1 - r_{S}, where r_{S} is the Spearman’s rank correlation between all row pairs of D. Next, we applied the iterative hierarchical clustering algorithm, using d_{S} as input and *average-linkage* method (Gordon et al., 2016) to define the distance between two clusters. This method merges node pairs into corresponding clusters by decreasing similarity until all nodes are merged into one cluster. The different stages of the algorithm were represented in the form of a dendrogram. We partitioned the case-control MSTs in an equal number of clusters using the same cut-off value (0.2) on each case-control dendrogram.

### Statistical Analysis

Global network connectivity (strength, diversity, and zero correlation), and MST global parameters (maximum degree, maximum betweenness, diameter, eccentricity, assortativity, degree divergence, leaf fraction, extreme value index, R-squared), provided a dataset of 80 rows (subjects) and 12 columns (parameters) for further statistical analysis. Since data were generally not Gaussian, non-parametric, Wilcoxon signed rank test was used to compare FTD and HC groups *P*-values were adjusted using Benjamini–Hochberg correction, and fixing the significance threshold at *P* < 0.05 (two sided). Sex and age showed significant differences between groups, therefore data were corrected for age, sex and age^{∗}sex interaction.

Global measures from the two case-control average correlation matrices were also compared by permutation tests as follows: (i) by computing the observed absolute difference between the global measure in FTD and HC groups, (ii) by permuting group assignments of the individuals’ values of the global measures for FTD and HC groups (*B* = 10000 iterations), and (iii) by repeating step (i) to obtain *B* = 10000 sampled permutations of the absolute differences between FTD and HC groups. Then, *P*-values were obtained using the sample permutation distribution with the same significant threshold of the individual signed rank tests.

### Methodological Comparison

Minimum spanning tree properties were compared with other conventional graph theory approaches, to measure the extent of network metrics reproducibility and wavelet scale specificity. The first method we applied uses the efficiency cost optimization (ECO; De Vico Fallani et al., 2017) criterion, imposing a fixed edge density threshold, based on the trade-off between network efficiency and wiring cost. We further applied two methods defining per subject optimal correlation thresholds, based on the extended Bayesian information criterion (EBIC; Chen and Chen, 2008), and spectral analysis (Perkins and Langston, 2009). Finally, we tested the performances of a scale free model-based method that chooses the correlation threshold optimizing power law fitting for each subject (Mumford et al., 2010). Differently from MST topology, these networks include triangular connections and cycles (i.e., they are not acyclic graphs), where classical clustering-based indices can be calculated. MST-specific indices are leaf fraction and alpha, while non-MST indices include clustering coefficient, average path length, and efficiency. Common metrics include maximum degree, maximum betweenness, degree divergence, diameter, eccentricity, and assortativity. Wilcoxon rank sum tests were calculated for every method at each wavelet scale.

### Software

Network analyses, graph visualization, and statistical analyses were performed in R (R Core team, 2018), using packages *igraph* (Csardi and Nepusz, 2006), *WGCNA* (Langfelder and Horvath, 2008), *brainwaver* (Achard, 2015), *brainGraph* (Watson, 2018), and custom R functions.

## Results

### Global Functional Connectivity of rs-fMRI Data

For each subject, we examined the complexity of rs-fMRI data (Saba et al., 2018) using bivariate measures (i.e., strength, diversity, and zero correlation) computed as summary regional values of the wavelet correlations (Supplementary Table 2). Scale 2 was the only reaching significance at these three tests (Figure 2 and Supplementary Table 2). Diversity and zero correlation showed a significant increase in bvFTD respect to HC, while strength showed significant decrease (*P* < 0.05 for every test). The significant increment of wavelet correlation diversity and its percentage of zeros (median HC: 25% vs. median bvFTD: 35%) indicated a decreased heterogeneity and increased null functional connectivity between brain regions in bvFTD compared to HCs. In addition, strength decrease denounces a generalized connectivity weakening in bvFTD respect to HCs.

**Figure 2.** Box-plots and *P*-values of Wilcoxon rank sum test for global connectivity measures and minimum spanning tree global parameters. The null hypothesis of the test is that the distributions of the HC and bvFTD do not differ (i.e., true location shift equal to 0). Data were corrected by age, sex, and age^{∗}sex interaction, and *P*-values were adjusted by Benjamini–Hochberg correction. Legend: rzero, zero correlation; kmax, maximum degree; Bmax, maximum betweenness; d, diameter; Ecc, eccentricity; Ass, assortativity; K, degree divergence; Lf, leaf fraction; alpha, extreme value index; R2 alpha, R-square of the extreme value index.

### Stronger Disorder Limit

By analyzing TOM-based weighted fully connected graphs for each subject (see Materials and Methods), we obtained good fitting for the power function (Figure 2 and Supplementary Table 2). High *R*^{2} (R-squared) indices were observed in both groups (0.93–0.97 for HC and 0.90–0.99 for bvFTD). When the extreme value index was estimated ($\widehat{\text{\alpha}}$), excluding outliers, low values (0.14–0.57 for HC and 0.10–0.42 for bvFTD) were found. Values of $\widehat{\text{\alpha}}$ less than one indicated a strong disorder limit tendency.

### MST Global Graph Metrics

Minimum spanning trees indicated brain connections alterations in bvFTD patients when compared to controls (Figure 2 and Supplementary Table 2). Data suggested a reduction of the degree centrality and leaf fraction, and an increment of distance metrics in the bvFTD group. Their trees were composed by nodes with a lower maximum degree and number of leafs. Moreover, the trees were characterized by a higher inter-distance, translating in a higher diameter and eccentricity in bvFTD compared to HC trees. These measures indicated brain impairments in bvFTD, highlighting less node-connections and loss of efficiency in exchange information capacity, that support a linear-shaped configuration network. Conversely, HC tree metrics showed a better network integration, characterized by parameter values that tend to a star-type configuration (i.e., an increase in the number of hubs and leaf points, and a decrease of diameter and eccentricity). Although assortativity index was not significant, it confirmed the topological hierarchy and biological nature of all MSTs, with assortativity values less than 0, for both groups. Lastly, permutation testing of two-group differences on the average correlation matrices yielded similar results (data not shown).

### MST Topological Two-Group Comparison

Cluster, spatial, and anatomical data were integrated through a network representation, showing the topological properties of the two-group average for HC and bvFTD graphs (see Materials and Methods section). The results of edge partitioning and node clustering are shown in Figures 3A,B, 4A,B. Every node in Figures 3A,B, colored by cluster membership, correspond to a single brain area belonging to a specific macro-region or lobe (i.e., node name), and traversed by two kinds of functional connections: (i) superhighways (bold gray), and (ii) secondary functional routes (thin gray). Finally, Figures 5A,B displays lobe partition (node color), node degree centrality (node size), and superhighways information flow connectivity (edge thickness) for bvFTD and HC groups, respectively.

**Figure 3.** **(A,B)** Clusters and superhighways differences in brain macro-regions. MSTs for both bvFTD and HCs are shown on the left **(A)** and right **(B)** panel, respectively. Clusters are defined by node colors, while thick edges show superhighway paths. Hubs (i.e., nodes having degree centrality >5) are marked with a black dot. Nodes are labeled according to macro-regions (i.e., lobes) membership (see Supplementary Table 1 for label encodings).

**Figure 4.** **(A,B)** Heatmaps and dendrograms of brain areas clustering elaboration. Hierarchical subdivision of brain areas and clustering for both bvFTD and HCs are shown on the left **(A)** and right **(B)** panel, respectively. The same number of clusters for both groups is obtained through the application of a cutting-height on dendrograms. We applied a cut-off equal to 0.2 to obtain four clusters (yellow-to-red squares within clustering areas). The four clusters (showed as colored boxes), reveal a clearer subdivision in the bvFTD group **(A)** respect to HCs **(B)**.

**Figure 5.** **(A,B)** Node degree centrality and superhighway routes in brain macro-regions. MSTs for both bvFTD and HCs are shown on the left **(A)** and right **(B)** panel, respectively. Macro-regions (i.e., lobes) are defined by node colors. Thick edges show superhighway paths and node size is proportional to degree centrality. Nodes are labeled by AAL 116 region identifier. Identifiers and corresponding areas and macro-regions can be found in Supplementary Table 1. Macro-region color code: Red, Frontal; Green, Temporal; Yellow, Parietal; Orange, Limbic; Cyan, Occipital; Brown, Subcortical Gray Matter; Light pink, Insula; Blue, Cerebellum.

In Figures 3A,B, 5A,B emerged the HC multiple-star structure, whose central nodes, namely Lingual-L, Occipital-Inf-R, Precuneus-L, and Cerebellum-6-L regions (nodes #47, #54, #67, and #99 in Figure 5B), serve as starting points for superhighways and bridge components in the tree backbone. This architecture was impaired in the MST of bvFTD patients, where few conserved stars, namely Precuneus-L, and Cerebellum-6-L (nodes #67 and #99 in Figure 5A), maintained a reduced functional connectivity, leading to a general isolation of Frontal and Temporal areas from central nodes. Network complexity reduction from HC to bvFTD is evident from bvFTD hierarchical clustering in Figure 3A, where Frontal and Temporal areas (green and red clusters, respectively) are almost completely separated from Parietal (yellow cluster) and Occipital-Cerebellar regions (cyan cluster). On the other hand, HC heatmap in Figure 3B, shows two highly connected network communities (green-brown and red-cyan clusters). Specifically, two different groups of Frontal-Parietal regions (brown and red clusters) are highly connected to Frontal (green cluster) and Temporal-Occipital-Cerebellar regions (cyan cluster), respectively. Notably, network star nodes (i.e., Lingual-L, Occipital-Inf-R, Precuneus-L, and Cerebellum-6-L) traverse and integrate these two highly connected communities in HC (Figure 5B). Conversely, bvFTD clusters are much more homogeneous (i.e., areas from the same lobe tend to cluster together), involving fewer and isolated stars. This is evident in Figure 5, where the bvFTD network panel (A) shows a clear lobe segregation (i.e., node color), especially for Frontal (red), Parietal (yellow), Occipital (light blue), and Cerebellum (blue) lobes, while the HC network panel (B) shows a much higher level of integration.

The survival ratio (i.e., the intersection between graphs calculated as the fraction of links found common in two MSTs) was equal to 42%, defining dissimilar topological structures of MSTs in the two groups. These differences are highlighted in Figure 6, where edges present in the bvFTD graph but not in the HC one panel (A), and vice versa panel (B), are shown. More specifically, highly connected star-clusters seen in the HC group were absent in the bvFTD graph. Conversely, leafs and short linear structures seen in the bvFTD group were absent in the HC graph.

**Figure 6.** **(A,B)** Graph-edges difference between the two groups. MSTs differences for both bvFTD and HCs are shown. Left **(A)** panel shows bvFTD-HC residual graph, while right **(B)** panel reports the difference HC-bvFTD. Node colors follow cluster membership, as in Figure 3. Nodes are labeled according to macro-regions membership (see Supplementary Table 1 for label encodings).

Nodes and edges in the axial orientation (x-y)-coordinates of the anatomical automatic labeling (AAL 116) brain atlas are shown in Figures 7A,B, providing a complete frontotemporal brain state representation, through the visualization of type, number, and origin of connections, confirming the abnormalities found in bvFTD, compared to HC. The most evident feature is a massive grouping of Frontal areas (red nodes in Figure 7A), disconnected from both Temporal lobes (green nodes in Figure 7A) and other areas (yellow nodes in Figure 7A), including the conserved star nodes Precuneus-L and Cerebellum-6-L (nodes #67 and #99 in Figure 7A). Strikingly, the massive frontal aggregation in bvFTD belongs entirely to a single cluster (red nodes in Figures 5A, 7A), at the center of which is present a new bvFTD-specific star node (Frontal Sup-Medial-L, node #23 in Figures 5A, 7A), indicative of a new isolated functional macro-region in bvFTD. Notably, comparing superhighways distribution between groups (Figures 7A,B), HCs show a deeply intertwined connectivity that integrates star nodes with Frontal and Temporal areas, linking them each other. Conversely, bvFTD superhighway connectivity collapses around conserved hubs and within the Frontal macro-region, causing their isolation.

**Figure 7.** **(A,B)** Spatial location of areas and superhighway routes according to x-y AAL 116 coordinates. MSTs for both bvFTD and HCs are shown on the left **(A)** and right **(B)** panel, respectively. Thick edges show superhighway paths and node size is proportional to degree centrality. Red nodes correspond to Frontal areas, while Temporal areas are colored in green. Hubs (i.e., nodes having degree centrality >5) are labeled by AAL 116 region identifier. Left MST **(A)** reports conserved hubs (areas #67 and #99), lost hubs (areas #47 and #54), and the gained frontal hub (area #23) in bvFTD, respect to HCs.

### Methodological Comparison

Minimum spanning tree network metrics were compared to other classical graph theory approaches (see Materials and Methods), to assess network properties reproducibility and wavelet scale specificity. Network metrics and Wilcoxon rank sum test results are reported in Supplementary Table 2. Beside the MST, only the efficiency/cost trade-off based method (i.e., ECO) showed significant bvFTD/HC specifically at scale 2, while all the other scales were non-informative (i.e., non-significant bvFTD/HC differences). The EBIC-based method reported few significant indices not directly related to information exchange efficiency (i.e., maximum degree, degree divergence, and assortativity), for both scales 1 and 2, indicating poor discriminant power. Lastly, spectral analysis and scale free model-based methods did not show significant case/control differences at any wavelet scale. Notably, both MST and ECO are based on the optimization of the information flow exchange across the functional network, maximizing graph integration and minimizing wiring costs. Common metrics between MST and ECO showing significant results at scale 2 were the same in both methods. However, differently from MST, ECO yielded several disconnected nodes (up to 74) per subject.

## Discussion

In the present study, we evaluated the power of MST representation within the framework of bvFTD, combining rs-fMRI data and graph theory analysis. Different neuroimaging methods have been proposed to highlight brain damage in bvFTD (Whitwell et al., 2011; Rohrer et al., 2015). However, since the description of the BOLD signal (Ogawa et al., 1990), functional neuroimaging allowed to go beyond the mere anatomical description of brain connectivity, identifying functionally connected (i.e., synchronized time-dependent fluctuations of the BOLD signal) networks of cortical and subcortical regions (Seeley et al., 2008). In bvFTD, multiple independent studies identified the frontal brain regions as the affected core networks (Seeley et al., 2009; Zhou et al., 2010; Whitwell et al., 2011; Borroni et al., 2012; Farb et al., 2013; Lee et al., 2014).

Graph theory has already been applied in a few MRI studies on bvFTD, demonstrating disruption of the global topologic organization, increased path length and assortativity, with loss of cortical hubs and network centrality, extending the list of sensible target regions to salience and executive functions (Agosta et al., 2013; Trojsi et al., 2015; Sedeño et al., 2016; Hohenfeld et al., 2018). Despite the wide use of graph theory in the analysis of neurodegenerative disorders, many methodological aspects are left to the experimenter choice, leading to possible confounding results.

When comparing graphs with the same number of nodes (*N*) and edges (*N*-1), the MST is an unambiguous method for brain network analysis, allowing to avoid methodological biases (Tewarie et al., 2015). Although it does not completely replace traditional graph theory approaches, MST remains the simplest and most effective representation of a full graph (Stam et al., 2014), where minimum weight (i.e., maximum connectivity) links constitute significant information flow paths. Moreover, it has been recently demonstrated how MST provides robust network estimates, with results in accordance to classical network analytical data (Otte et al., 2015; Tewarie et al., 2015). As indicated by our methodological comparison (see Results), algorithms based on information flow optimization and wiring cost minimization (including MST and ECO methods) achieved the best bvFTD/HC separation performances (Supplementary Table 2). On the other hand, fixed-correlation threshold and scale free model-based graphs showed poor or non-significant results, indicating that enforcing an arbitrary correlation threshold or network nature (i.e., scale free power law distribution) strongly limit network descriptive power. Conversely, MST has the advantage of requiring no correlation thresholds, network density or *a priori* distribution, ensuring full reproducibility and robustness in different conditions. Furthermore, since our results indicated a strong disorder limit tendency in both groups, the MST was able to preserve the connectivity features of the underlying functional networks (Tewarie et al., 2014). Through our TOM-adjusted edge weights, MSTs and relative parameters conserved their neighborhood node characteristics, highlighting nodes aggregation (i.e., star-type configurations), and providing a valid method to identify a set of shortest paths in MSTs that may be considered as the “critical backbone” of original graph (Tewarie et al., 2014).

Our MSTs suggest that the differences between groups may be attributed to functional alterations of such major organization, since they represent the information-flow highways of the fully connected network. Shape-linear configuration tendency in bvFTD graphs highlights different impairments: high distance between nodes, low centrality parameter values, and a low exchange information capacity (i.e., low network integration). Connection efficiency loss is particularly evident in Figure 7, where the superhighway system in HCs, linking hubs to Frontal and Temporal brain areas, is replaced by a local (i.e., isolated) network surrounding conserved hubs. Functional isolation is a generalized process in bvFTD, where brain areas tend to interact within lobes (i.e., colors in Figure 5), showing a homogeneous brain area distribution, longer distances between hubs, and longer within-lobe superhighways. Therefore, bvFTD functional breakdown is not merely described by connectivity loss, but though disease-specific reorganizations and regularization of the information flow. This contrasts with the marked integration of a healthy functional network, where superhighways serve as shortcuts to connect areas from different brain macro-regions. Network regularization has already been observed as a distinctive FTD trait, respect to other neurodegenerative disorders, including Alzheimer’s disease (de Haan et al., 2009; Zhou et al., 2010). We further investigated this aspect by gathering evidences from both global and local network metrics.

Although global functional parameters (i.e., strength, diversity, and zero correlation) showed a significantly weaker and reduced connectivity in bvFTD, edge-level and node-level features (i.e., superhighways, and node degree and betweenness centrality), highlighted a more complex scenario, explaining some of the key dysfunctions observed in large scale resting-state functional networks, including the DMN, SN, and EN networks (Zhou et al., 2010; Trojsi et al., 2015; Sedeño et al., 2016; Hohenfeld et al., 2018). The first evidence from our data is the formation of a new FTD-specific hub (area #23: Frontal_Sup_Medial_L, Figure 5A), absent in HCs. This hub is the starting point of a huge superhighway, fully extending within frontal lobe, clearly derived from an elongation of the original route in HCs (Figure 5B). Specifically, our data showed the involvement of regions: Frontal_Sup_Medial_L/R, Frontal_Sup_Orb_L/R, Frontal_Mid_Orb_L/R, Frontal_Med_Orb_L/R, Rectus_L/R, and Olfactory_L/R (areas #23, #24, #5, #6, #9, #10, #25, #26, #27, #28, #21, and #22 in Supplementary Table 1 and Figure 5A). Notably, all the elements of this superhighway are part of the DMN (Zhou et al., 2010; Trojsi et al., 2015; Hohenfeld et al., 2018), strongly supporting the evidence of a compensation mechanism and frontal DMN decreased connectivity with areas from other lobes. However, this process is not exclusive of the Frontal lobe. Within-lobe superhighway formation, and consequent isolation, involves also the two conserved nodes in bvFTD: Precuneus_L (area #67, parietal lobe in yellow in Figure 5A) and Cerebellum_6_L (area #99, blue in Figure 5A), where the former plays a key role in the DMN network. Betweenness difference (*b*) is a good local indicator of these impairments at node level (threshold set at *b* = ± 435, (see Supplementary Table 3). Notably, while areas #23 and #67 show a strongly increased node betweenness (*b _{23}* = 1545, and

*b*= 879) from bvFTD to HCs, area #99 loses a great portion of its centrality (

_{67}*b*= -1534), supporting the evidence of a functional deterioration of the cerebellar lobule VI, that has been associated with cerebellar atrophy in both bvFTD and Alzheimer’s disease (Guo et al., 2016; Schmahmann, 2016). Nevertheless, area #99 is the center of an enlarged intra-cerebellar superhighway system, showing that superhighway elongation is due to a generalized network centrality reorganization, rather than a region-specific impairment. It has been recently suggested that long-distance connections have an important role in integrating distinct brain areas, leading to a greater functional diversification, robustness, and specialization (Betzel and Bassett, 2018). However, long-range connections number and length is strictly controlled by their metabolic demand. In contrast, bvFTD long within-lobe superhighways seem not to contribute to the overall functional integration, but rather being the result of a compensatory mechanism, in response to a generalized functional deterioration. To verify this hypothesis, we focused our attention to those areas showing a strong betweenness centrality loss (

_{99}*b*< -435) in bvFTD respect to HCs. The clearest example in our data is given by three connected areas (Figure 5B): Angular_R (area #66,

*b*= -1552), Frontal_Sup_R (area #4,

_{66}*b*= -1359), and Cingulum_Ant_R (area #32,

_{4}*b*= -1299), being involved in DMN, EN, and EN/SN, respectively. As shown in Figure 5B, these three areas connect the Precuneus_L (parietal area #67) and its superhighway system to the frontal superhighway starting from area #23 (Frontal_Sup_Medial_L), that in bvFTD is markedly enlarged (Figure 5A). In the bvFTD network, area #23 is a new hub, and the #23–67 connection is now a much longer and linear path, suggesting that superhighway elongation could be a compensatory reaction to a less efficient network integration. Notably, Cingulum_Ant_R (area #32) is part of the limbic system, involved in both DMN and salience/executive functions, suggesting that limbic system failure could be an underlying cause of the global network rearrangement observed in bvFTD subjects.

_{32}A further support to this hypothesis is represented by the disruption of two HC network hubs (Figure 5B): Occipital_Inf_R (area #54, *b _{54}* = -4054) and Lingual_L (area #47,

*b*= -2913), replaced by a long linear-shaped sequence of occipital areas (Figure 5A), denouncing a massive loss of degree and betweenness centrality also in this part of the network. Beside hubs, two adjacent connector areas experienced huge loss of betweenness: Fusiform_R (area #56,

_{47}*b*= -531) and ParaHippocampal_R (area #40,

_{56}*b*= -540). Strikingly, while the former is a direct bridge to the hub group, the latter is the first node of a limbic superhighway (Figure 5B), whose nodes occupy distinct peripheral positions in bvFTD (Figure 5A), including: Hippocampus_L/R (area #37–38), ParaHippocampal_L/R (areas #39–40), and Amygdala_L/R (areas #41–42).

_{40}Collectively, these evidences show an underlying involvement of the limbic system in the observed bvFTD functional deterioration, associated to the well-studied impairments affecting emotion recognition, social inference, and executive functions typical of this neurodegenerative disorder (Zhou et al., 2010; Trojsi et al., 2015; Sedeño et al., 2016; Hohenfeld et al., 2018).

## Conclusion

The present work had the primary goal of highlighting alterations in brain connectivity of bvFTD subjects, providing at the same time a detailed description of the observed functional impairments, and insights about their possible causes. According to the most recent fMRI literature (Guo et al., 2017; Cui et al., 2018), we applied an MST model to wavelet correlation matrices from bvFTD and HC subjects, exploiting three strongpoints of MST-based methods: (i) assumption-free network construction and reproducibility, (ii) independence from node and edge number during network comparison, and (iii) simplicity of representation (i.e., three branching is directly interpretable in terms of most efficient shortest paths). On the other hand, MSTs have one main limitation: the resulting network is a sparse representation, implicitly excluding triangular connections, thus causing non-applicability of some common clustering metrics (e.g., transitivity and coreness) and evaluation of the network small-worldness. However, we coped with this issue by using TOM (Zhang and Horvath, 2005) and the extreme value index evaluation (Van Mieghem and Magdalena, 2005).

The combination of this theoretical model with rs-fMRI data allowed us not only to generate a clear picture of the functional divergence of bvFTD from HCs, but also to shed light on the possible causes of topological and functional rearrangements, and compensatory mechanisms, underlying cognitive, social, and executive impairments characterizing bvFTD phenotype.

Further developments to the present work, that now constitute main limitations, are represented by: (i) the lack of clinical and/or metabolic parameters that could confirm or reveal new causal hypotheses, and (ii) a combination with anatomical variables (e.g., gray matter mass), to achieve a better model resolution, and associate degenerative processes to functional deterioration.

Collectively, the application of MST-based analysis to rs-fMRI data looks a promising way to clarify the role of degenerative processes involved in FTD functional breakdown, improving the discovery of new fMRI biomarkers.

## Data Availability

Pre-processed and parceled data (ItalianFTD-fMRI_May2018) for this study can be found in the SOURCEFORGE repository at: https://sourceforge.net/projects/bionet-finder/.

## Author Contributions

VS performed the literature research, drafted the manuscript, performed the data analysis and R coding, and wrote, reviewed, and edited the manuscript. EP conceived and designed the study, acquired the data, analyzed the image, and wrote, reviewed, and edited the manuscript. VC performed the literature research and analyzed the image. FP revised the data analysis, results, and discussion, and wrote, reviewed, and edited the manuscript. SG, OZ, RG, and AP acquired the data. BB conceived and designed the study, recruited the patients, drafted the manuscript, and supervised the study. MG conceived and designed the study, administered the project, performed the data analysis and R coding, and wrote, reviewed, and edited the manuscript. All authors have read and approved the final manuscript.

## Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors (Materials and Methods section, Subjects).

## Conflict of Interest Statement

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.

## Acknowledgments

The authors would like to thank patients, their families and all volunteers for allowing and supporting this research. The authors would also like to acknowledge the support of Raffaele Ferrari (Department of Molecular Neuroscience, University College London, London, United Kingdom) for scientific discussion and revision of the article.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2019.00211/full#supplementary-material

## Abbreviations

AAL, automated anatomical labeling; bvFTD, behavioral variant Frontotemporal Dementia; HC, healthy controls; MST, minimum spanning tree; ROI, regions of interest; rs-fMRI, resting state functional MRI; SPT, shortest path tree; TOM, topological overlap measures.

## Footnotes

## References

Achard, S. (2015). *Basic Wavelet Analysis of Multivariate Time Series With a Visualisation and Parametrisation Using Graph Theory. CRAN Repository.* Available at: https://cran.r-project.org/web/packages/brainwaver/brainwaver.pdf

Achard, S., Salvador, R., Whitcher, B., Suckling, J., and Bullmore, E. (2006). A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. *J. Neurosci.* 26, 63–72. doi: 10.1523/JNEUROSCI.3874-05.2006

Agosta, F., Sala, S., Valsasina, P., Meani, A., Canu, E., Magnani, G., et al. (2013). Brain network connectivity assessed using graph theory in frontotemporal dementia. *Neurology* 81, 134–143. doi: 10.1212/WNL.0b013e31829a33f8

Albert, R., and Barabási, A. L. (2002). Statistical mechanics of complex networks. *Rev. Mod. Phys.* 74, 47–97. doi: 10.1103/RevModPhys.74.47

Antonenko, D., Külzow, N., Sousa, A., Prehn, K., Grittner, U., and Flöel, A. (2018). Neuronal and behavioral effects of multi-day brain stimulation and memory training. *Neurobiol. Aging* 61, 245–254. doi: 10.1016/j.neurobiolaging.2017.09.017

Beckmann, C. F., and Smith, S. M. (2004). Probabilistic independent component analysis for functional magnetic resonance imaging. *IEEE Trans. Med. Imaging* 23, 137–152. doi: 10.1109/TMI.2003.822821

Betzel, R. F., and Bassett, D. S. (2018). Specificity and robustness of long-distance connections in weighted, interareal connectomes. *PNAS* 115, E4880–E4889. doi: 10.1073/pnas.1720186115

Borroni, B., Alberici, A., Cercignani, M., Premi, E., Serra, L., Cerini, C., et al. (2012). Granulin mutation drives brain damage and reorganization from preclinical to symptomatic FTLD. *Neurobiol. Aging* 33, 2506–2520. doi: 10.1016/j.neurobiolaging.2011.10.031

Braunstein, L. A., Wu, Z., Chen, Y., Buldyrev, S. V., Kalisky, T., Sreenivasan, S., et al. (2007). Optimal path and minimal spanning trees in random weighted networks. *Int. J. Bifurcat. Chaos* 17, 2215–2255. doi: 10.1142/S0218127407018361

Brett, M., Anton, J.-L., Valabregue, R., and Poline, J.-B. (2002). Region of interest analysis using an SPM toolbox. *Neuroimage* 16, 1140–1141. doi: 10.1016/S1053-8119(02)90013-3

Bullmore, E., Fadili, J., Maxim, V., Sendur, L., Whitcher, B., Suckling, J., et al. (2004). Wavelets and functional magnetic resonance imaging of the human brain. *Neuroimage* 23(Suppl. 1), S234–S249. doi: 10.1016/j.neuroimage.2004.07.012

Bullmore, E., and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. *Nat. Rev. Neurosci.* 10, 186–198. doi: 10.1038/nrn2575

Bullmore, E., and Sporns, O. (2012). The economy of brain network organization. *Nat. Rev. Neurosci.* 13, 336–349. doi: 10.1038/nrn3214

Chen, J., and Chen, Z. (2008). Extended Bayesian information criteria for model selection with large model spaces. *Biometrika* 95, 759–771. doi: 10.1093/biomet/asn034

Chiang, S., Cassese, A., Guindani, M., Vannucci, M., Yeh, H. J., Haneef, Z., et al. (2016). Time-dependence of graph theory metrics in functional connectivity analysis. *Neuroimage* 125, 601–615. doi: 10.1016/j.neuroimage.2015.10.070

Cormen, T. H., Leiserson, C. E., Rivest, R. L., and Stein, C. (2001). *Introduction to Algorithms*, 2nd Edn. Cambridge, MA: MIT Press.

Cosseddu, M., Benussi, A., Gazzina, S., Turrone, R., Archetti, S., Bonomi, E., et al. (2018). Mendelian forms of disease and age at onset affect survival in frontotemporal dementia. *Amyotroph. Lateral Scler. Frontotemporal Degener.* 19, 87–92. doi: 10.1080/21678421.2017.1384020

Csardi, G., and Nepusz, T. (2006). The igraph software package for complex network research. *InterJournal Complex Syst.* 1695, 1–9.

Cui, X., Xiang, J., Guo, H., Yin, G., Zhang, H., Lan, F., et al. (2018). Classification of Alzheimer’s disease, mild cognitive impairment, and normal controls with subnetwork selection and graph kernel principal component analysis based on minimum spanning tree brain functional network. *Front. Comput. Neurosci.* 12:31. doi: 10.3389/fncom.2018.00031

de Haan, W., Pijnenburg, Y. A., Strijers, R. L., van der Made, Y., van der Flier, W. M., Scheltens, P., et al. (2009). Functional neural network analysis in frontotemporal dementia and Alzheimer’s disease using EEG and graph theory. *BMC Neurosci.* 10:101. doi: 10.1186/1471-2202-10-101

De Vico Fallani, F., Latora, V., and Chavez, M. (2017). A topological criterion for filtering information in complex brain networks. *PLoS Comput. Biol.* 13:e1005305. doi: 10.1371/journal.pcbi.1005305

Drakesmith, M., Caeyenberghs, K., Dutt, A., Lewis, G., David, A. S., and Jones, D. K. (2015). Overcoming the effects of false positives and threshold bias in graph theoretical analyses of neuroimaging data. *Neuroimage* 118, 313–333. doi: 10.1016/j.neuroimage.2015.05.011

Edwards, D., de Abreu, G. C. G., and Labouriau, R. (2010). Selecting high-dimensional mixed graphical models using minimal AIC or BIC forests. *BMC Bioinformatics* 11:18. doi: 10.1186/1471-2105-11-18

Farb, N. A. S., Grady, C. L., Strother, S., Tang-Wai, D. F., Masellis, M., Black, S., et al. (2013). Abnormal network connectivity in frontotemporal dementia: evidence for prefrontal isolation. *Cortex* 49, 1856–1873. doi: 10.1016/j.cortex.2012.09.008

Filippi, M., Basaia, S., Canu, E., Imperiale, F., Meani, A., Caso, F., et al. (2017). Brain network connectivity differs in early-onset neurodegenerative dementia. *Neurology* 89, 1764–1772. doi: 10.1212/WNL.0000000000004577

Gazzina, S., Premi, E., Zanella, I., Biasiotto, G., Archetti, S., Cosseddu, M., et al. (2016). Iron in frontotemporal lobar degeneration: a new subcortical pathological pathway? *Neurodegener Dis.* 16, 172–178. doi: 10.1159/000440843

Gordon, E., Rohrer, J. D., and Fox, N. C. (2016). Advances in neuroimaging in frontotemporal dementia. *J. Neurochem.* 138(Suppl. 1), 193–210. doi: 10.1111/jnc.13656

Gorno-Tempini, M. L., Hillis, A. E., Weintraub, S., Kertesz, A., Mendez, M., Cappa, S. F., et al. (2011). Classification of primary progressive aphasia and its variants. *Neurology* 76, 1006–1014. doi: 10.1212/WNL.0b013e31821103e6

Griffanti, L., Dipasquale, O., Laganà, M. M., Nemni, R., Clerici, M., Smith, S. M., et al. (2015). Effective artifact removal in resting state fMRI data improves detection of DMN functional connectivity alteration in Alzheimer’s disease. *Front. Hum. Neurosci.* 9:449. doi: 10.3389/fnhum.2015.00449

Griffanti, L., Salimi-Khorshidi, G., Beckmann, C. F., Auerbach, E. J., Douaud, G., Sexton, C. E., et al. (2014). ICA-based artefact removal and accelerated fMRI acquisition for improved resting state network imaging. *Neuroimage* 95, 232–247. doi: 10.1016/j.neuroimage.2014.03.034

Guo, C. C., Tan, R., Hodges, J. R., Hu, X., Sami, S., and Hornberger, M. (2016). Network-selective vulnerability of the human cerebellum to Alzheimer’s disease and frontotemporal dementia. *Brain* 139, 1527–1538. doi: 10.1093/brain/aww003

Guo, H., Liu, L., Chen, J., Xu, Y., and Jie, X. (2017). Alzheimer classification using a minimum spanning tree of high-order functional network on fMRI dataset. *Front. Neurosci.* 11:639. doi: 10.3389/fnins.2017.00639

Hohenfeld, C., Werner, C. J., and Reetz, K. (2018). Resting-state connectivity in neurodegenerative disorders: is there potential for an imaging biomarker? *Neuroimage Clin.* 18, 849–870. doi: 10.1016/j.nicl.2018.03.013

Jenkinson, M., Beckmann, C. F., Behrens, T. E. J., Woolrich, M. W., and Smith, S. M. (2012). FSL. *Neuroimage* 62, 782–790. doi: 10.1016/j.neuroimage.2011.09.015

Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. *BMC Bioinformatics* 9:559. doi: 10.1186/1471-2105-9-559

Lee, S. E., Khazenzon, A. M., Trujillo, A. J., Guo, C. C., Yokoyama, J. S., Sha, S. J., et al. (2014). Altered network connectivity in frontotemporal dementia with C9orf72 hexanucleotide repeat expansion. *Brain* 137, 3047–3060. doi: 10.1093/brain/awu248

Lynall, M.-E., Bassett, D. S., Kerwin, R., McKenna, P. J., Kitzbichler, M., Muller, U., et al. (2010). Functional connectivity and brain networks in schizophrenia. *J. Neurosci.* 30, 9477–9487. doi: 10.1523/JNEUROSCI.0333-10.2010

Mann, D. M. A., and Snowden, J. S. (2017). Frontotemporal lobar degeneration: pathogenesis, pathology and pathways to phenotype. *Brain Pathol.* 27, 723–736. doi: 10.1111/bpa.12486

McKeown, M. J., Hansen, L. K., and Sejnowsk, T. J. (2003). Independent component analysis of functional MRI: what is signal and what is noise? *Curr. Opin. Neurobiol.* 13, 620–629. doi: 10.1016/j.conb.2003.09.012

Mears, D., and Pollard, H. B. (2016). Network science and the human brain: using graph theory to understand the brain and one of its hubs, the amygdala, in health and disease. *J. Neurosci. Res.* 94, 590–605. doi: 10.1002/jnr.23705

Meier, J., Tewarie, P., and Van Mieghem, P. (2015). The union of shortest path trees of functional brain networks. *Brain Connect.* 5, 575–581. doi: 10.1089/brain.2014.0330

Mumford, J. A., Horvath, S., Oldham, M. C., Langfelder, P., Geschwind, D. H., and Poldrack, R. A. (2010). Detecting network modules in fMRI time series: a weighted network analysis approach. *Neuroimage* 52, 1465–1476. doi: 10.1016/j.neuroimage.2010.05.047

Newman, M. E. J. (2003). Mixing patterns in networks. *Phys. Rev. E* 67(2 Pt 2):026126. doi: 10.1103/PhysRevE.67.026126

Numan, T., Slooter, A. J. C., van der Kooi, A. W., Hoekman, A. M. L., Suyker, W. J. L., Stam, C. J., et al. (2017). Functional connectivity and network analysis during hypoactive delirium and recovery from anesthesia. *Clin. Neurophysiol.* 128, 914–924. doi: 10.1016/j.clinph.2017.02.022

Ogawa, S., Lee, T. M., Kay, A. R., and Tank, D. W. (1990). Brain magnetic resonance imaging with contrast dependent on blood oxygenation. *Proc. Natl. Acad. Sci. U.S.A.* 87, 9868–9872. doi: 10.1073/pnas.87.24.9868

Otte, W. M., van Diessen, E., Paul, S., Ramaswamy, R., Subramanyam Rallabandi, V. P., Stam, C. J., et al. (2015). Aging alterations in whole-brain networks during adulthood mapped with the minimum spanning tree indices: the interplay of density, connectivity cost and life-time trajectory. *Neuroimage* 109, 171–189. doi: 10.1016/j.neuroimage.2015.01.011

Park, B.-Y., Moon, T., and Park, H. (2018). Dynamic functional connectivity analysis reveals improved association between brain networks and eating behaviors compared to static analysis. *Behav. Brain Res.* 337, 114–121. doi: 10.1016/j.bbr.2017.10.001

Perkins, A. D., and Langston, M. A. (2009). Threshold selection in gene co-expression networks using spectral graph theory techniques. *BMC Bioinformatics* 10:S4. doi: 10.1186/1471-2105-10-S11-S4

R Core team (2018). *R: A Language and Environment for Statistical Computing.* Vienna: The R Foundation for Statistical Computing.

Raichle, M. E. (2015). The brain’s default mode network. *Annu. Rev. Neurosci.* 38, 433–447. doi: 10.1146/annurev-neuro-071013-014030

Rascovsky, K., Hodges, J. R., Knopman, D., Mendez, M. F., Kramer, J. H., Neuhaus, J., et al. (2011). Sensitivity of revised diagnostic criteria for the behavioural variant of frontotemporal dementia. *Brain* 134, 2456–2477. doi: 10.1093/brain/awr179

Rashid, B., Damaraju, E., Pearlson, G. D., and Calhoun, V. D. (2014). Dynamic connectivity states estimated from resting fMRI Identify differences among Schizophrenia, bipolar disorder, and healthy control subjects. *Front. Hum. Neurosci.* 8:897. doi: 10.3389/fnhum.2014.00897

Ravasz, E., Somera, A. L., Mongru, D. A., Oltvai, Z. N., and Barabási, A. L. (2002). Hierarchical organization of modularity in metabolic networks. *Science* 297, 1551–1555. doi: 10.1126/science.1073374

Rohrer, J. D., Nicholas, J. M., Cash, D. M., van Swieten, J., Dopper, E., Jiskoot, L., et al. (2015). Presymptomatic cognitive and neuroanatomical changes in genetic frontotemporal dementia in the genetic frontotemporal dementia initiative (GENFI) study: a cross-sectional analysis. *Lancet Neurol.* 14, 253–262. doi: 10.1016/S1474-4422(14)70324-2

Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. *Neuroimage* 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003

Saba, V., Premi, E., Cristillo, V., Gazzina, S., Palluzzi, F., Zanetti, O., et al. (2018). *ItalianFTD-fMRI_May2018. SourceForge.* Available at: https://sourceforge.net/projects/bionet-finder/

Schmahmann, J. D. (2016). Cerebellum in Alzheimer’s disease and frontotemporal dementia: not a silent bystander. *Brain* 139, 1314–1318. doi: 10.1093/brain/aww064

Sedeño, L., Couto, B., García-Cordero, I., Melloni, M., Baez, S., Morales Sepúlveda, J. P., et al. (2016). Brain network organization and social executive performance in frontotemporal dementia. *J. Int. Neuropsychol. Soc.* 22, 250–262. doi: 10.1017/S1355617715000703

Seeley, W. W., Crawford, R., Rascovsky, K., Kramer, J. H., Weiner, M., Miller, B. L., et al. (2008). Frontal paralimbic network atrophy in very mild behavioral variant frontotemporal dementia. *Arch. Neurol.* 65, 249–255. doi: 10.1001/archneurol.2007.38

Seeley, W. W., Crawford, R. K., Zhou, J., Miller, B. L., and Greicius, M. D. (2009). Neurodegenerative diseases target large-scale human brain networks. *Neuron* 62, 42–52. doi: 10.1016/j.neuron.2009.03.024

Stam, C. J., Tewarie, P., Van Dellen, E., van Straaten, E. C. W., Hillebrand, A., and Van Mieghem, P. (2014). The trees and the forest: characterization of complex brain networks with minimum spanning trees. *Int. J. Psychophysiol.* 92, 129–138. doi: 10.1016/j.ijpsycho.2014.04.001

Telesford, Q. K., Joyce, K. E., Hayasaka, S., Burdette, J. H., and Laurienti, P. J. (2011). The ubiquity of small-world networks. *Brain Connect.* 1, 367–375. doi: 10.1089/brain.2011.0038

Tewarie, P., Hillebrand, A., Schoonheim, M. M., van Dijk, B. W., Geurts, J. J. G., Barkhof, F., et al. (2014). Functional brain network analysis using minimum spanning trees in multiple sclerosis: an MEG source-space study. *Neuroimage* 88, 308–318. doi: 10.1016/j.neuroimage.2013.10.022

Tewarie, P., van Dellen, E., Hillebrand, A., and Stam, C. J. (2015). The minimum spanning tree: an unbiased method for brain network analysis. *Neuroimage* 104, 177–188. doi: 10.1016/j.neuroimage.2014.10.015

Trojsi, F., Esposito, F., de Stefano, M., Buonanno, D., Conforti, F. L., Corbo, D., et al. (2015). Functional overlap and divergence between ALS and bvFTD. *Neurobiol. Aging* 36, 413–423. doi: 10.1016/j.neurobiolaging.2014.06.025

Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. *Neuroimage* 15, 273–289. doi: 10.1006/nimg.2001.0978

van Dellen, E., Douw, L., Hillebrand, A., de Witt Hamer, P. C., Baayen, J. C., Heimans, J. J., et al. (2014). Epilepsy surgery outcome and functional network alterations in longitudinal MEG: a minimum spanning tree analysis. *Neuroimage* 86, 354–363. doi: 10.1016/j.neuroimage.2013.10.010

van Diessen, E., Otte, W. M., Stam, C. J., Braun, K. P. J., and Jansen, F. E. (2016). Electroencephalography based functional networks in newly diagnosed childhood epilepsies. *Clin. Neurophysiol.* 127, 2325–2332. doi: 10.1016/j.clinph.2016.03.015

van Lutterveld, R., van Dellen, E., Pal, P., Yang, H., Stam, C. J., and Brewer, J. (2017). Meditation is associated with increased brain network integration. *Neuroimage* 158, 18–25. doi: 10.1016/j.neuroimage.2017.06.071

Van Mieghem, P., and Magdalena, S. M. (2005). Phase transition in the link weight structure of networks. *Phys. Rev. E Stat. Nonlin. Soft Matter Phys.* 72:056138. doi: 10.1103/PhysRevE.72.056138

Van Mieghem, P., and van Langen, S. (2005). Influence of the link weight structure on the shortest path. *Phys. Rev. E Stat. Nonlin. Soft Matter Phys.* 71:056113. doi: 10.1103/PhysRevE.71.056113

van Wijk, B. C. M., Stam, C. J., and Daffertshofer, A. (2010). Comparing brain networks of different size and connectivity density using graph theory. *PLoS One* 5:e13701. doi: 10.1371/journal.pone.0013701

Watson, C. G. (2018). *Graph Theory Analysis of Brain MRI Data. CRAN Repository.* Available at: https://github.com/cwatson/brainGraph

Whitwell, J. L., Josephs, K. A., Avula, R., Tosakulwong, N., Weigand, S. D., Senjem, M. L., et al. (2011). Altered functional connectivity in asymptomatic MAPT subjects: a comparison to bvFTD. *Neurology* 77, 866–874. doi: 10.1212/WNL.0b013e31822c61f2

Wu, Z., Braunstein, L., Havlin, S., and Stanley, H. (2006). Transport in weighted networks: partition into superhighways and roads. *Phys. Rev. Lett.* 96:148702. doi: 10.1103/PhysRevLett.96.148702

Yu, M., Gouw, A. A., Hillebrand, A., Tijms, B. M., Stam, C. J., van Straaten, E. C. W., et al. (2016). Different functional connectivity and network topology in behavioral variant of frontotemporal dementia and Alzheimer’s disease: an EEG study. *Neurobiol. Aging* 42, 150–162. doi: 10.1016/j.neurobiolaging.2016.03.018

Yu, M., Hillebrand, A., Tewarie, P., Meier, J., van Dijk, B., Van Mieghem, P., et al. (2015). Hierarchical clustering in minimum spanning trees. *Chaos* 25:023107. doi: 10.1063/1.4908014

Zalesky, A., Fornito, A., and Bullmore, E. (2012). On the use of correlation as a measure of network connectivity. *Neuroimage* 60, 2096–2106. doi: 10.1016/j.neuroimage.2012.02.001

Zhang, B., and Horvath, S. (2005). A general framework for weighted gene co-expression network analysis. *Stat. Appl. Genet. Mol. Biol.* 4:Article17. doi: 10.2202/1544-6115.1128

Zhang, Z., Telesford, Q. K., Giusti, C., Lim, K. O., and Bassett, D. S. (2016). Choosing wavelet methods, filters, and lengths for functional brain network construction. *PLoS One* 11:e0157243. doi: 10.1371/journal.pone.0157243

Keywords: functional connectivity, functional magnetic resonance imaging, resting state, minimum spanning tree, graph theory, behavioral variant frontotemporal dementia, neurodegeneration

Citation: Saba V, Premi E, Cristillo V, Gazzina S, Palluzzi F, Zanetti O, Gasparotti R, Padovani A, Borroni B and Grassi M (2019) Brain Connectivity and Information-Flow Breakdown Revealed by a Minimum Spanning Tree-Based Analysis of MRI Data in Behavioral Variant Frontotemporal Dementia. *Front. Neurosci.* 13:211. doi: 10.3389/fnins.2019.00211

Received: 17 September 2018; Accepted: 25 February 2019;

Published: 14 March 2019.

Edited by:

Zhi Yang, Shanghai Mental Health Center, ChinaReviewed by:

Hyunjin Park, Sungkyunkwan University, South KoreaSunghyon Kyeong, Yonsei University College of Medicine, South Korea

Copyright © 2019 Saba, Premi, Cristillo, Gazzina, Palluzzi, Zanetti, Gasparotti, Padovani, Borroni and Grassi. 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: Valentina Saba, valentina.saba01@universitadipavia.it