Abstract
A network neuroscience perspective can significantly advance the understanding of neurodegenerative disorders, particularly frontotemporal dementia (FTD). This study employed an innovative multiplex connectomics approach, integrating cortical thickness (CTH) and fluorodeoxyglucose-positron emission tomography (FDG-PET) in a dual-layer model to investigate network alterations in FTD subtypes across two geographically distinct sites. The cohort included groups of behavioral variant FTD (bvFTD), primary progressive aphasia (PPA), mild cognitive impairment (MCI), and cognitively normal (CN) individuals who were analyzed from two separate sites. Site 1 included 28 bvFTD, 20 PPA, and 27 MCI participants, whereas Site 2 included 26 bvFTD, 43 PPA, and 43 CN individuals, respectively. Utilizing CTH and FDG-PET data after standard preprocessing, a multiplex network pipeline in BRAPH2 toolbox was used to derive multiplex participation coefficient (MPC) between the groups. The analysis revealed an increase in global MPC as an indicator of disease in PPA at both sites. Additionally, nodal MPC alterations in the anterior cingulate, frontal, and temporal lobes in PPA were compared with bvFTD. Comparisons with the CN showed that nodal MPC alterations were more extensive in PPA when compared with bvFTD. These findings underscore the potential utility of multiplex connectomes for identifying network disruptions in neurodegenerative disorders, offering promising implications for future research and clinical applications.
Author Summary
This study uses multiplex connectomics, a method combining cortical thickness and glucose metabolism data, to investigate frontotemporal dementia (FTD) across two sites. Behavioral variant FTD (bvFTD) and primary progressive aphasia showed global and nodal brain network alterations. PPA exhibited more widespread disruptions in the anterior cingulate, frontal, and temporal lobes than bvFTD, suggesting greater network changes. bvFTD showed distinct alterations in the left caudal anterior cingulate and temporal pole compared with mild cognitive impairment and more extensive changes compared with cognitively normal individuals. These findings highlight the value of multiplex connectomics in identifying complex brain network changes in FTD and emphasize the need for further research to uncover the mechanisms behind these alterations.
INTRODUCTION
The rapid growth of India’s aging population, coupled with a rise in modifiable risk factors—air pollution, hypertension, traumatic brain injury, and obesity—presents a significant challenge for the nation’s public health and social care systems (Livingston et al., 2020). Neurodegenerative disorder is a major cause of cognitive disability and dependency among the elderly and is expected to affect millions of Indian families over the next 2 decades (Ravindranath & Sundarakumar, 2021). Frontotemporal dementia (FTD) is an early-onset neurodegenerative disorder that primarily affects the frontal and temporal lobes of the brain leading to a range of behavioral and language impairments (Rascovsky et al., 2011; Warren et al., 2013). FTD group has a behavioral variant FTD (bvFTD) presenting with alterations in personality, cognition, and behavior, and a primary progressive aphasia (PPA) variant predominantly presenting with disorders in language and speech (Gorno-Tempini et al., 2011; Rascovsky et al., 2011). Survival differs across these clinical groups, with patients diagnosed with the semantic variant of PPA having an average survival of up to 12.2 years, while those with bvFTD have a survival of around 8.2 years (Kansal et al., 2016). Pathologically, FTD has the intracellular deposition of Tau, TAR DNA binding protein-43, or fused in sarcoma (RNA binding protein; Mackenzie et al., 2010) in the majority of patients. It is important to note the lack of correlation between pathological subtypes with clinical syndromes, wherein one clinical subtype could have several pathological subtypes and one pathological subtype could have varying clinical presentations.
Multilayer network approaches integrate structural and metabolic connectivity, providing a more comprehensive understanding of the brain’s network alterations than single-layer methods. Studies have shown that multilayer approaches can reveal hidden network properties related to disease progression and cognitive decline such as Alzheimer’s disease (AD) (Canal-Garcia et al., 2022). In the context of bvFTD and PPA, the added value lies in the potential to capture complex alterations that may not be detectable in single-layer metrics. These multilayer analyses can help in identifying specific biomarkers for differentiating these syndromes, particularly through the examination of region-specific nodal changes.
Neuroimaging is a key component in the diagnostic workup of patients with FTD and reveals characteristic anterior cingulate and fronto-insular hypometabolism or atrophy in bvFTD and asymmetric temporal prominent atrophy in the language variant of PPA supporting the clinical syndromes of FTD (Cerami et al., 2017; Gorno-Tempini et al., 2004; Rogalski et al., 2011). These regions are consistent with disease epicenters representing neurons with heightened vulnerability to disease especially in the early stages of the disease onset (Seeley et al., 2009). Existing research using whole-brain, data-driven, graph theory methods (Bullmore & Sporns, 2009) utilizing measures of segregation (clustering coefficient, local efficiency, modularity), integration (path-length, global efficiency, participation coefficient), and centrality (betweenness centrality) when compared with healthy controls have revealed significant anterior predominant alterations in frontal, temporal, and limbic regions in bvFTD (Nigro et al., 2021) and bilateral left more than right frontotemporal regions in PPA (Nigro et al., 2022; Tao et al., 2020) reiterating the epicenter concept of the disease. Several studies also have investigated network differences when comparing disease-matched controls. The comparison of AD with bvFTD has shown a higher participation coefficient in the left frontomedial cortices in AD compared with bvFTD patients (Ng et al., 2021). Resting-state functional MRI (rs-fMRI) studies have shown significant network alterations in the parietal lobe in early-onset AD compared with bvFTD (Filippi et al., 2017). Unilayer networks in neuroimaging refer to the analysis within a single layer of data, either from fluorodeoxyglucose-positron emission tomography (FDG-PET), rs-fMRI, or structural MRI, measuring regional variations of glucose metabolism, BOLD changes, or cortical thickness (CTH). A multilayer network, however, consists of several distinct unilayer networks, which use identical regions of interest (ROI) as nodes and with the concept of adjacency matrix (De Domenico, 2017), which measures edges between layers to identify regions with similar network measures across modalities (Battiston et al., 2014). Several network measures, like multiplex participation coefficient (MPC), multiplex clustering coefficient, and multiplex overlapping degree, can be derived from such networks (Canal-Garcia et al., 2022). Nodal MPC is an integration measure that can be used to detect highly connected hubs in both layers, and increasing nodal MPC indicates densely connected hubs that are similar in both layers. A higher MPC in multiplex terminology does not reflect increased hubness of the nodes; it reflects the differential contribution of the layers as the edges define the participation coefficient strictly between layers. Multiplex clustering coefficient is a segregation measure of the number of triangles formed by the nodes across layers, indicating their cliqueness. Multiplex overlapping degrees is another measure of segregation, which is the sum of the degree of the node in both layers. Few studies have used a multiplex model using rs-fMRI, electroencephalogram, or magnetoencephalography in dementia (Klepl et al., 2023; Wang et al., 2021). A recent study using CTH and amyloid PET using multiplex graph measures in AD with amyloid beta (Aβ+) patients found alterations in MPC in several regions including the medial temporal lobe, cuneus, and anterior cingulate regions (Canal-Garcia et al., 2022). Another study using multiple frequency bands in rs-fMRI study as different layers (Wang et al., 2021) revealed reduced MPC in AD compared with controls.
Given the varying clinical profile, imaging, and network features in patients with bvFTD and PPA, in the current study, we used the advanced imaging modalities in the clinical care of patients—CTH and FDG-PET—within a multiplex connectome to explore differences in the multilayer graph measures between bvFTD and PPA. Our study used mild cognitive impairment (MCI) to compare with dementia patients. We also used publicly available Neuroimaging in Frontotemporal Dementia (NIFD) data in patients with bvFTD, PPA, and cognitively normal (CN) to evaluate the consistency of the pipeline. We hypothesized that bvFTD will show significant differences in multiplex network measures compared with PPA based on the existing clinical and imaging evidence.
METHODS
Demographic Information of Participants
In this ambispective comparative cohort study of consecutive patients referred for positron emission tomography-magnetic resonance imaging between 2016 and 2022, we examined a total of 75 patients (from Site 1: 28 bvFTD, 20 PPA, and 27 MCI). The severity of dementia was assessed using the Clinical Dementia Rating (CDR) scale. Cognitive evaluations were conducted using Addenbrooke’s Cognitive Examination-III (ACE-III; Mekala et al., 2020). The diagnosis of FTD was made in patients using the standard criteria (Gorno-Tempini et al., 2011; Rascovsky et al., 2011), while the diagnosis of MCI was made based on the modified Petersen criteria (Petersen, 2004). In the retrospective group from Site 1 (31 subjects), CDR scales and ACE-III scores were extracted from the patient case files, and imaging data were obtained from the imaging data bank. For all the subjects in the prospective group (44 subjects), informed consent, and clinical and behavioral scores were obtained in person. Additionally, for retrospective cases, we received the approval from Institute Ethics Committee for waiver of consent. The study protocol was approved by this committee (Ethics NO. NIMH/DO/[BS&NS DIV.] 2019–20). All methods were carried out in compliance with the approved Indian Council of Medical Research guidelines. For the purpose of comparison with global data, we used the 112 patients from NIFD to replicate the pipeline (site 2, 26 bvFTD, 43 PPA, and 43 CN; https://adni.loni.usc.edu). Since ACE-III was unavailable for the Site 2 data, Mini-Mental State Examination (MMSE) scores were used for cognitive evaluation.
Comparing the groups, the MCI from Site 1 was comparable to bvFTD and PPA (Table 1) in age, gender, and disease duration. bvFTD and PPA revealed significant differences in ACE-III score compared with MCI and bvFTD had significantly lower years of education compared with MCI (bvFTD = 12.86 ± 3.99, MCI = 15.56 ± 3.10; Table 1). In Site 2, there were no significant differences in gender and years of education between these groups; however, the bvFTD was significantly younger than PPA and CN (bvFTD = 61.85 ± 6.56, PPA = 66.81 ± 8.04, and CN = 62.16 ± 7.66 years). MMSE scores exhibited no significant difference between bvFTD and PPA (Table 2). The effects of age, gender, and education were regressed from the analysis at the time of network construction, to ensure that the results were unaffected by these confounds.
Demographic, clinical characteristics MCI, bvFTD, and PPA groups from Site 1
Characteristic . | MCI . | FTD . | . | F . | PR(> F) . | Post hoc (Bonferroni) . | ||||
---|---|---|---|---|---|---|---|---|---|---|
bvFTD . | PPA . | |||||||||
N | 27 | 28 | 20 | -- | -- | -- | ||||
Gender (M:F) | 20:07 | 18:10 | 13:07 | 7.61 | 0.02 | NS | ||||
Mean | SD | Mean | SD | Mean | SD | |||||
Age (years) | 63.52 | 9.02 | 63.07 | 9.08 | 59.20 | 6.74 | 1.65 | 0.20 | NS | |
Education (years) | 15.56 | 3.05 | 12.86 | 3.99 | 14.65 | 3.99 | 3.67 | 0.03 | MCI > bvFTD | |
Duration of illness (years) | 2.24 | 1.89 | 3.02 | 1.40 | 3.10 | 1.17 | 2.33 | 0.10 | NS | |
ACE III | 81.13 | 11.41 | 58.04 | 13.92 | 39.55 | 21.84 | 39.68 | <0.001 | MCI > bvFTD, MCI > PPA & bvFTD > PPA | |
CDR | Very mild | 27 | 05 | 05 | 23.38 | <0.001 | PPA > MCI & bvFTD > MCI | |||
Mild | -- | 15 | 09 | |||||||
Moderate | -- | 06 | 06 | |||||||
Severe | -- | 02 | -- |
Characteristic . | MCI . | FTD . | . | F . | PR(> F) . | Post hoc (Bonferroni) . | ||||
---|---|---|---|---|---|---|---|---|---|---|
bvFTD . | PPA . | |||||||||
N | 27 | 28 | 20 | -- | -- | -- | ||||
Gender (M:F) | 20:07 | 18:10 | 13:07 | 7.61 | 0.02 | NS | ||||
Mean | SD | Mean | SD | Mean | SD | |||||
Age (years) | 63.52 | 9.02 | 63.07 | 9.08 | 59.20 | 6.74 | 1.65 | 0.20 | NS | |
Education (years) | 15.56 | 3.05 | 12.86 | 3.99 | 14.65 | 3.99 | 3.67 | 0.03 | MCI > bvFTD | |
Duration of illness (years) | 2.24 | 1.89 | 3.02 | 1.40 | 3.10 | 1.17 | 2.33 | 0.10 | NS | |
ACE III | 81.13 | 11.41 | 58.04 | 13.92 | 39.55 | 21.84 | 39.68 | <0.001 | MCI > bvFTD, MCI > PPA & bvFTD > PPA | |
CDR | Very mild | 27 | 05 | 05 | 23.38 | <0.001 | PPA > MCI & bvFTD > MCI | |||
Mild | -- | 15 | 09 | |||||||
Moderate | -- | 06 | 06 | |||||||
Severe | -- | 02 | -- |
Demographic, clinical characteristics of CN, bvFTD, and PPA groups from Site 2
Characteristic . | CN . | . | FTD . | F . | PR(> F) . | Post hoc (Bonferroni) . | ||||
---|---|---|---|---|---|---|---|---|---|---|
bvFTD . | PPA . | |||||||||
N | 43 | 26 | 43 | -- | -- | -- | ||||
Gender (M:F) | 17:26 | 17:09 | 24:19 | 4.79 | 0.09 | NS | ||||
Mean | SD | Mean | SD | Mean | SD | |||||
Age (years) | 62.16 | 7.66 | 61.85 | 6.56 | 66.81 | 8.04 | 5.13 | 0.01 | PPA > bvFTD & PPA > CN | |
Education (years) | 20.93 | 17.36 | 15.65 | 3.17 | 16.02 | 2.80 | 2.72 | 0.07 | NS | |
MMSE | 28.51 | 5.24 | 21.23 | 7.39 | 17.65 | 9.22 | 22.69 | <0.001 | CN > bvFTD & CN > PPA | |
CDR | None | 43 | -- | -- | 92.91 | <0.001 | bvFTD > CN, PPA > CN & bvFTD > PPA | |||
Very mild | -- | 01 | 24 | |||||||
Mild | -- | 09 | 10 | |||||||
Moderate | -- | 16 | 09 | |||||||
Severe | -- | -- | -- |
Characteristic . | CN . | . | FTD . | F . | PR(> F) . | Post hoc (Bonferroni) . | ||||
---|---|---|---|---|---|---|---|---|---|---|
bvFTD . | PPA . | |||||||||
N | 43 | 26 | 43 | -- | -- | -- | ||||
Gender (M:F) | 17:26 | 17:09 | 24:19 | 4.79 | 0.09 | NS | ||||
Mean | SD | Mean | SD | Mean | SD | |||||
Age (years) | 62.16 | 7.66 | 61.85 | 6.56 | 66.81 | 8.04 | 5.13 | 0.01 | PPA > bvFTD & PPA > CN | |
Education (years) | 20.93 | 17.36 | 15.65 | 3.17 | 16.02 | 2.80 | 2.72 | 0.07 | NS | |
MMSE | 28.51 | 5.24 | 21.23 | 7.39 | 17.65 | 9.22 | 22.69 | <0.001 | CN > bvFTD & CN > PPA | |
CDR | None | 43 | -- | -- | 92.91 | <0.001 | bvFTD > CN, PPA > CN & bvFTD > PPA | |||
Very mild | -- | 01 | 24 | |||||||
Mild | -- | 09 | 10 | |||||||
Moderate | -- | 16 | 09 | |||||||
Severe | -- | -- | -- |
Acquisition Parameters
Magnetization prepared rapid acquisition gradient-recalled echo (MPRAGE).
At Site 1, subjects underwent structural MRI on a 3 Tesla (Biograph mMR, Siemens, Erlangen, Germany) scanner using a 16-channel transmit receive head coil. Anatomic images were obtained using a three-dimensional T1-weighted MPRAGE sequence; additional conventional sequences, fluid-attenuated inversion recovery, T2, and susceptibility weighted imaging were acquired; and structural and unrelated metabolic abnormalities were excluded in all patients by domain experts (e.g., Dr. R.D.B., Dr. S.M., Dr. C.N.). At Site 2, subjects were imaged on a 3 Tesla scanner (Discovery MR 750, GE Medical Systems), where 3D T1-weighted MPRAGE sequences were also obtained. Detailed acquisition parameters are provided in Supporting Information Table S1.
FDG-PET.
Subjects were administered FDG 185 ± 10%, Megabecquerel (5.0 mCi) intravenously 30 min before acquiring PET scan in list mode for 15 min. Standard corrections for attenuation, scatter, random coincidences, decay, and a 5- mm Gaussian postreconstruction filter were applied. The acquisition parameters in Site 2 were identical, and a comparative parameter table is provided in Supporting Information Table S1. There were no differences in the acquisition parameters between the retrospective and prospective groups of patients. Since Site 2 had 30 min of PET data, only the first 15 min of data was used for analysis to make it comparable with Site 1 data in the current study.
Preprocessing
CTH.
The process of reconstructing the cortex and determining CTH was carried out using the FreeSurfer image analysis tool version 6.05 (Fischl, 2012; https://surfer.nmr.mgh.harvard.edu). The pipeline within the FreeSurfer suite involves an automated process for removing the skull, transforming images into Talairach space, segmenting subcortical gray and white matter, and delineating the boundaries between gray and white matter. Subsequently, the cortical surface underwent reconstruction through precise intersubject alignment methods, resulting in a detailed measurement of CTH and categorization into various predefined cortical brain regions as per different atlases (Dale et al., 1999). CTH was determined as the minimum distance between the junction of gray and white matter and the pial surface. The results of these preliminary processing stages were visually examined to confirm that the analyses were conducted accurately. To ensure consistency with PET segmentation, the cortical surface was divided into unique areas based on gyral and sulcal structures using the Desikan-Killiany (DK) atlas, which identifies 34 regions in each hemisphere (Desikan et al., 2006).
FDG-PET.
PETSurfer analysis tool v6.05 was used for FDG-PET image analysis (Greve et al., 2014). PETSurfer performs PET processing by utilizing the result from FreeSurfer, involving registration into the Montreal Neurological Institute space. Further, it provides partial volume correction using geometric transfer matrix methods for the partial volume effect, which can impact the accuracy of PET measurements (Greve et al., 2016). After segmentation and registration, the pipeline applies PET Surfer and calculates the standardized uptake value ratio (SUVR) values with pons as a reference region, volume size in the region, and voxel variance. The DK atlas was used to parcellate the cortical surface into 34 regions per hemisphere (Desikan et al., 2006; Khokhar et al., 2023). Additionally, we also utilized the Destrieux atlas, which comprises 74 regions per hemisphere, to assess the stability of our findings and to explore whether consistent results could be obtained using two different atlases. The results based on the Destrieux atlas are included in the supplementary material for reference. Further, we quantified the metabolic index for FDG-PET using the approach described by Thibaut et al. (2021) to generate a mean SUV image representing FDG uptake (Supporting Information Figure S3).
Multilayer Graph Theory
The network construction procedure involves the following steps:
Schematic representation of the image processing and multilayer graph theory analysis pipeline. CTH and FDG-PET values are extracted for each ROI. Pearson’s correlation coefficients are calculated between all ROI pairs within each group, generating group-level connectivity matrices for each modality (CTH and FDG-PET). These intralayer matrices are then integrated into a supra-adjacency matrix to facilitate multiplex graph analysis. (N represents number of nodes in matrix.)
Schematic representation of the image processing and multilayer graph theory analysis pipeline. CTH and FDG-PET values are extracted for each ROI. Pearson’s correlation coefficients are calculated between all ROI pairs within each group, generating group-level connectivity matrices for each modality (CTH and FDG-PET). These intralayer matrices are then integrated into a supra-adjacency matrix to facilitate multiplex graph analysis. (N represents number of nodes in matrix.)
Overlapping degree (global).
MPC (nodal).
MPC (global).
Global MPC represents the mean value of the nodal participation coefficients across the various layers of a network.
Statistical Analysis
To assess differences between groups in demographic variables, statistical analysis was performed with “scipy” library in Python. Differences between groups in age, education, and ACE-III/MMSE were tested using an analysis of variance for continuous variables, while gender differences between groups were tested using a chi-square test. Post hoc analyses using Bonferroni-corrected pairwise t tests, following a Kruskal–Wallis test for the ordinal variable CDR, were conducted into pairwise group differences with a significance level of 0.05.
For multiplex analysis, we carried out a nonparametric permutation test with 5,000 permutations. Group differences were considered significant for a two-tailed p value based on 95% confidence intervals. FDR corrections (Benjamini & Hochberg, 1995) were applied on multiplex graph measures to control for multiple comparisons (p < 0.05). The multiplex analysis results were visualized with “fsbrain” library, a visualization for brain statistics R-package (Schäfer & Ecker, 2020) using the DK atlas.
RESULTS
In this multisite analysis of multiplex network from sMRI and FDG-PET, PPA revealed significantly increased global MPC when compared with MCI and CN in Site 1 and Site 2, respectively. Nodal MPC revealed more widespread alterations when compared with controls with significant alterations in anterior cingulate and temporal lobes in PPA when compared with bvFTD in both Site 1 and Site 2.
Multiplex Global Measures
Overlapping degree average (global).
There was a significantly reduced global overlapping degree in patients with PPA when compared with CN in Site 2 (d = 0.7, diff = 0.13). However, when comparing bvFTD to MCI, bvFTD to CN, PPA to MCI, and bvFTD to PPA (Sites 1 & 2), no significant differences were observed in the overlapping degree.
MPC (global).
There was an increase in the global MPC in patients with bvFTD when compared with PPA in Site 2 (d = 0.4, 0.45, 0.5, 0.55, diff = −0.20, −0.24, −0.25, −0.27). PPA, when compared with MCI in Site 1 (d = 0.6, 0.65, 0.7, diff = 0.44, 0.47, 0.40), showed increase global MPC, also seen in bvFTD and PPA when compared with CN (d = 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7; diff = 0.22, 0.26, 0.32, 0.38, 0.45, 0.65, 0.7 and d = 0.7; diff = 0.23, respectively). When comparing bvFTD to MCI, and bvFTD to PPA in Site 1, no significant differences were observed in the global MPC.
MPC (Nodal).
Nodal MPC identified several alterations in the frontal and temporal lobes in both sites.
bvFTD versus PPA.
In Site 1, bvFTD showed increased MPC in the left entorhinal and parahippocampal regions compared with PPA, while the right rostral and caudal anterior cingulate regions had lower MPC in bvFTD compared with PPA (Figure 2A). In Site 2, the left parahippocampal region, similar to Site 1, and the right-rostral and caudal anterior cingulate regions exhibited increased MPC in bvFTD compared with PPA (Figure 2B). However, the left-inferior temporal, fusiform, and frontal pole regions showed lower MPC in bvFTD compared with PPA. As shown in Figures 2C and 2D, these connectivity differences indicate relative variations in the contribution of different layers between the groups. Figures 2C and 2D, where the increased nodal MPC in bvFTD in the left parahippocampal gyrus was mediated by the increased contribution of the FDG in Site 1 and CTH in Site 2, compared to PPA. When it comes to the right rostral and caudal anterior cingulate, it can be seen that the decreased MPC in Site 1 is mediated by the prominent contribution of the FDG layer and the increased MPC in Site 2 is due to the increased contribution of the CTH layer.
Nodal MPC of significant differences between the bvFTD and PPA from Site 1 (A), and bvFTD versus PPA from Site 2 (B), two layers of connectivity imbalance can be seen for Site 1 and Site 2 (C and D, respectively) between bvFTD and PPA (FDR corrected at p < 0.05).
Nodal MPC of significant differences between the bvFTD and PPA from Site 1 (A), and bvFTD versus PPA from Site 2 (B), two layers of connectivity imbalance can be seen for Site 1 and Site 2 (C and D, respectively) between bvFTD and PPA (FDR corrected at p < 0.05).
PPA versus MCI/CN.
PPA revealed wider areas of altered MPC than bvFTD when compared with MCI/CN. Right-rostral anterior cingulate, caudal anterior cingulate, left lateral orbitofrontal, pars triangularis, rostral middle frontal, and banks of the superior temporal sulcus showed increased MPC in PPA when compared with MCI. The right entorhinal and left temporal pole showed decreased MPC in PPA when compared with MCI. In Site 2 data, PPA revealed increased nodal MPC in right-rostral anterior cingulate, parahippocampal, isthmus cingulate, inferior temporal, middle temporal, fusiform, frontal pole, superior frontal, left-inferior temporal, isthmus cingulate, middle temporal, pars opercularis, rostral middle frontal, caudal middle frontal, superior temporal, fusiform, left banks of the superior temporal sulcus, and frontal pole, superior frontal. Conversely, there were regions, the right—cuneus, pericalcarine, transverse temporal; left—pericalcarine and transverse temporal, that revealed decreased MPC in PPA (Figures 3A and 3B).
Nodal MPC shows the significant differences between the group PPA and MCI (A), PPA and CN (B), bvFTD and MCI (C), and bvFTD and CN (D) (FDR corrected at p < 0.05).
Nodal MPC shows the significant differences between the group PPA and MCI (A), PPA and CN (B), bvFTD and MCI (C), and bvFTD and CN (D) (FDR corrected at p < 0.05).
bvFTD versus MCI/CN.
The left caudal anterior cingulate and temporal pole showed significantly lower MPC in bvFTD in comparison with MCI in Site 1. Increased MPC was more diffuse and bilateral involving the frontal lobes and cingulate gyrus when compared with healthy controls (site 2) than MCI (site 1). Specifically in Site 2, the right-posterior cingulate, rostral anterior cingulate, caudal anterior cingulate, pars opercularis, caudal middle frontal, paracentral, superior frontal, left-lateral orbitofrontal, pars opercularis, pars triangularis, rostral middle frontal, caudal middle frontal, and superior frontal showed significantly higher MPC in bvFTD compared with CN (Figures 3C and 3D).
DISCUSSION
The complex nature of neurodegenerative disorders necessitates a connectomics approach to understand their pathobiology in vivo, and the current study with multiplex global and nodal connectomics across two geographically distinct sites has revealed an increasing global network MPC as a marker of the disease in patients with PPA in both the sites. Nodal MPC using data-driven methods consistently identified alterations of several frontal and temporal lobe regions in both sites, supporting the clinical and pathological concept of disease epicenters. Both PPA and bvFTD show increased global and nodal MPC in similar frontal regions when compared with CN individuals, suggesting a shared network-level disruption in frontal circuits across these disorders. However, PPA showed distinct increases in nodal MPC within temporal regions, distinguishing it from MCI and highlighting the diagnostic potential of temporal network alterations, likely reflecting compensatory network changes associated with disease progression. Additionally, differences in nodal MPC in parahippocampal regions emerged as a critical feature for differentiating bvFTD from PPA.
Only a few studies on neurodegenerative disorders, such as AD, have investigated multilayer connectomes (Canal-Garcia et al., 2022; Wang et al., 2021) and have found decreased nodal participation coefficient in several regions in the medial temporal, anterior cingulate, and cuneus in patients with AD Aβ+ compared with CN Aβ− participants due to differential contributions from amyloid PET and CTH. However, it needs to be noted that the increase in global MPC found in this study does not have the same biological meaning as the increased participation coefficient found in unilayer studies (Ng et al., 2021; Tao et al., 2021). Despite not having the biological correlate, multiplex connectomics is useful as it identifies areas that are involved across both layers providing summative evidence from multiple modalities.
bvFTD revealed significant alterations in the left caudal anterior cingulate and temporal pole when compared with MCI and revealed more widespread alterations involving various parts of the right cingulate gyrus, right paracentral, bilateral superior and middle frontal gyrus, pars opercularis and left pars triangularis and lateral orbitofrontal region when compared with CN. It is important to note that the data-driven method using multilayer connectomics has identified these areas that have been already identified in the literature using unilayer connectomics (Agosta et al., 2013; Mandelli et al., 2018; Reyes et al., 2018; Sedeño et al., 2017), pathology, and clinical methods.
It is interesting to note that PPA revealed more extensive alterations compared with bvFTD in both sites involving the right anterior cingulate and left-lateral orbitofrontal, pars triangularis, rostral middle frontal, and banks of the superior temporal sulcus when compared with the MCI group and involving several regions of bilateral cingulate isthmus, middle and inferior temporal, fusiform, superior frontal and frontal poles, left middle and inferior frontal, pars opercularis, transverse temporal, and pericalcarine regions. Asymmetric left predominant involvement was noted in the left middle and inferior frontal, pars opercularis, frontal pole, and superior temporal sulcus.
Contrasting the bvFTD and PPA with CN revealed more areas with alterations than with MCI. It is logical to consider these differences to be due to varying disease stages spotlighting the longitudinal evolution of the disease with more areas involved as the disease progresses, which is amplified when compared with CN and minimized in comparison with the MCI. Notably, the asymmetric left prominent involvement of the frontal and temporal lobes in PPA in both sites is also interesting. Asymmetric involvement of the left frontal and temporal lobes is one of the imaging markers for patients with PPA with left anterioinferior temporal lobe in semantic variant PPA and left side inferior frontal and insula in progressive nonfluent aphasia (Cerami et al., 2017; Gorno-Tempini et al., 2004; Rogalski et al., 2011). Though the regions involved closely map the existing literature in FTD, it is intriguing to note widespread alterations in PPA compared with bvFTD in both sites. This is counterintuitive since clinically PPA is a milder form of disease with a better prognosis than bvFTD. FDG-PET and conventional MRI also reveal more focal disease in PPA as compared with bvFTD. One explanation for this result could be that the network changes denote both disease-induced changes and compensatory network modulations due to the disease (Liu et al., 2022). Our analysis of bvFTD versus PPA revealed significant differences in the right cingulate gyrus and the left parahippocampal gyrus in both sites. Additional involvement of the left entorhinal, fusiform gyrus, and frontal pole was identified in Site 2. Additional involvement of the left temporal and frontal areas in Site 2 can be partially explained by the differences in the severity of the diseases. CDR scores from Site 2 revealed a higher number of PPA patients to belong to the minor group whereas the majority of bvFTD belonged to the severe group. The CDR scores across groups were more balanced in Site 1. Though it is logical to presume that the differences are reflective of the contributions by PPA, a detailed unilayer network analysis will be required to further understand the biological meaning of this finding as the multiplex methodology depict contributions of the layers and not the disease. Since the derived regions indicate group results, it is difficult to correlate these with clinical and behavioral scores.
Multiplex connectomics provides a platform for the integration of multimodal techniques. Several studies using unilayer CTH and FDG-PET networks have revealed highly heterogeneous results with the involvement of several areas in the literature (Malpetti et al., 2019; Vijverberg et al., 2017). Some studies have also observed heterogeneity between FDG-PET and CTH patterns in FTD patients (Bejanin et al., 2020; Cerami et al., 2016). Some studies reveal increased connectivity, and some others reveal decreased connectivity when compared with CN, justified by differences in the temporal evolution of the disease, and varying image processing pipelines and statistical thresholds used while deriving the results. Studies from geographically different locations also add to the variability necessitating large data samples to comprehend disease pathobiology. It is important to note that the current multiplex connectomics study across two geographically different sites revealed increased global MPC in FTD with altered nodal MPC in several regions involving the frontal and temporal lobes. Though it is difficult to derive the network pathobiology in FTD from this data, since the results provided consistent results across two distinct data sets, multiplex connectomics by its very characteristic of defining edges between layers appears to be an appropriate tool in addressing the heterogeneity in data sets.
Despite the substantial contributions of our study in comprehending network changes in dementia through the innovative application of a multiplex network approach, certain limitations need acknowledgment. First, in the current study, the DK atlas was utilized for the parcellation of CTH and FDG uptake regions. A significant limitation stems from this atlas’s inability to delineate the thickness of subcortical brain regions and the cerebellum. As a result, this confines our understanding of the whole brain’s structural and metabolic interactions to a certain extent. Further, we did an analysis using Destrieux atlas and the findings showed concordance with DK atlas, and we found increasing global MPC in patients with PPA, when compared with MCI and CN in Site 1 and Site 2, respectively, and in bvFTD when compared with CN in and Site 2 with no differences between bvFTD and MCI in site 1. Nodal MPC also revealed significant alterations in bvFTD and PPA when compared with CN and MCI in the frontotemporal region aligning with the DK atlas; however, in bvFTD, when compared with PPA in both Site 1 and Site 2, more regions in the neocortical frontal, occipital, temporal and cingulate, which unlike in DK atlas, were limited to cingulate and temporal regions. These differences could be attributed to increased granularity of the Destrieux Atlas with 148 ROI, underscoring the importance of parcellation in the nodal participation coefficient (supporting information results).
Second, extracting graph measures from structural and PET data on an individual basis proved challenging, preventing correlation with clinical and cognitive scores. Recent developments in FDG-PET time series (functional PET) acquisition, combined with MR-PET hybrid scanners, could enhance this approach at the subject level. Third, the absence of healthy control subjects at Site 1 posed a significant limitation, potentially introducing biases. Fourth, the analytical framework primarily hinged on cross-sectional imaging data for executing multiplex analyses, which inherently restricted the investigation of temporal variations of FDG-PET and gray matter atrophy across groups. Fifth, our findings are limited by a smaller sample size, which may affect their generalizability. Additionally, we have included both the semantic variant and the nonfluent variant of PPA under the PPA category due to smaller numbers. The biological meaning of the network alterations could not be derived because of the interlayer measures employed, thereby limiting the interpretability of the clinical context. Subsequent research should incorporate other imaging techniques, like fMRI, functional FDG-PET and Tau-PET or Translocator Protein-PET, into multiplex analyses to gain insights into the relationship between brain atrophy, pathology, and other pathological processes in FTD. Moreover, longitudinal study designs could be instrumental in exploring the temporal trajectory of tau pathologies vis-a-vis gray matter atrophy. If the pipeline in the future can, in addition, provide within-network measures, it will be more useful in deriving clinical correlates.
Conclusion
The current study using a multiplex connectomics approach revealed widespread changes in the MPC, in FTD subtypes consistently across geographically distinct sites using the most popular imaging tools—CTH and FDG-PET. These findings offer new insights into the interplay between metabolism and brain atrophy over the course of FTD and subtypes and move beyond traditional graph theory analyses based on unilayer brain networks.
Acknowledgments
We would like to acknowledge the participants who took part in this study and the staff at the PET-MRI facility for their assistance in data acquisition. We would like to express our gratitude to the National Institute of Mental Health and Neuro Sciences (NIMHANS) for providing the resources.
Frontotemporal Lobar Degeneration Neuroimaging Initiative (FTLDNI) was funded through the National Institute of Aging and started in 2010. The primary goals of FTLDNI were to identify neuroimaging modalities and methods of analysis for tracking frontotemporal lobar degeneration and to assess the value of imaging versus other biomarkers in diagnostic roles. The Principal Investigator of NIFD was Dr. Howard Rosen, MD at the University of California, San Francisco. The data are the result of collaborative efforts at three sites in North America. For up-to-date information on participation and protocol, please visit https://memory.ucsf.edu/research/studies/nifd. Grant Number: R01 AG032306.
Supporting Information
Supporting information for this article is available at https://doi.org/10.1162/netn_a_00448.
Author Contributions
Sunil Kumar Khokhar: Conceptualization; Data curation; Formal analysis; Methodology; Resources; Software; Visualization; Writing – original draft; Writing – review & editing. Manoj Kumar: Investigation; Supervision; Validation; Writing – review & editing. Faheem Arshad: Data curation; Investigation; Resources; Validation; Writing – review & editing. Sheetal Goyal: Data curation; Investigation; Validation. Megha Tiwari: Data curation; Investigation; Methodology. Nithin Thanissery: Data curation; Methodology; Writing – review & editing. Subasree Ramakrishnan: Data curation; Supervision. Chandana Nagaraj: Supervision; Validation. Rajan Kashyap: Formal analysis; Methodology; Validation. Sandhya Mangalore: Methodology; Resources; Writing – review & editing. Tapan K. Gandhi: Conceptualization; Methodology; Supervision; Validation; Writing – review & editing. Suvarna Alladi: Conceptualization; Data curation; Supervision; Validation; Writing – review & editing. Rose Dawn Bharath: Conceptualization; Investigation; Methodology; Resources; Supervision; Validation; Writing – review & editing.
Funding Statement
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. The authors independently conducted the study, and all necessary resources and materials were provided by National Institute of Mental Health and Neuro Sciences.
Author Disclosure Statement
The authors of this manuscript declare that they have no conflicts of interest to disclose concerning the contents of this research paper.
Data Availability
• Data from Site 1 analyzed during the current study are not publicly available due to conditions of our ethics approval, particularly regarding patient confidentiality, and do not permit public archiving of the original data but are available from the corresponding author on reasonable request.
• Data from Site 2 analyzed during the current study are available in the Frontotemporal Lobar Degeneration Neuroimaging Initiative database repository, available at https://adni.loni.usc.edu.
Author List of FTLDNI/NIFD Consortium
Howard Rosen; University of California, San Francisco (PI). Bradford C. Dickerson; Harvard Medical School and Massachusetts General Hospital. Kimoko Domoto-Reilly; University of Washington School of Medicine. David Knopman; Mayo Clinic, Rochester. Bradley F. Boeve; Mayo Clinic Rochester. Adam L. Boxer; University of California, San Francisco. John Kornak; University of California, San Francisco. Bruce L. Miller; University of California, San Francisco. William W. Seeley; University of California, San Francisco. Maria-Luisa Gorno-Tempini; University of California, San Francisco. Scott McGinnis; University of California, San Francisco. Maria Luisa Mandelli; University of California, San Francisco.
TECHNICAL TERMS
- Multiplex connectomics:
An approach that integrates multiple types of data (such as CTH and FDG-PET) into a unified network model to analyze brain connectivity patterns across different modalities.
- Multiplex network:
A multiplex network is a specific type of multilayer network where interlayer connections exist solely between the corresponding nodes across different layers.
- Multiplex participation coefficient (MPC):
A measure in multiplex network analysis that quantifies the extent to which a node (ROI) participates in different layers (modalities) of the network.
- Overlapping degree:
The overlapping degree of a graph is the sum of the degrees of a node in all layers.
- Supra-adjacency matrix:
A matrix representing the connections in a multiplex network, containing intralayer adjacency matrices on the diagonal and interlayer adjacency matrices on the off-diagonal.
REFERENCES
Competing Interests
Competing Interests: The authors have declared that no competing interests exist.
Author notes
Handling Editor: Alex Fornito