Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 30 June 2020
Sec. Alzheimer's Disease and Related Dementias
This article is part of the Research Topic Translational Advances in Alzheimer's, Parkinson's, and other Neurodegenerative Dementias View all 70 articles

Default Mode Network Analysis of APOE Genotype in Cognitively Unimpaired Subjects Based on Persistent Homology

\r\nLiqun Kuang*Liqun Kuang1*Jiaying JiaJiaying Jia1Deyu ZhaoDeyu Zhao1Fengguang XiongFengguang Xiong1Xie HanXie Han1Yalin Wang*Yalin Wang2*for the Alzheimer&#x;s Disease Neuroimaging Initiativefor the Alzheimer’s Disease Neuroimaging Initiative
  • 1School of Data Science and Technology, North University of China, Taiyuan, China
  • 2School of Computing, Informatics, and Decision Systems Engineering, Arizona State University, Tempe, AZ, United States

Current researches on default mode network (DMN) in normal elderly have mainly focused on finding some dysfunctional areas with decreased or increased connectivity. The global network dynamics of apolipoprotein E (APOE) e4 allele group is rarely studied. In our previous brain network study, we have demonstrated the advantage of persistent homology. It can distinguish robust and noisy topological features over multiscale nested networks, and the derived properties are more stable. In this study, for the first time we applied persistent homology to analyze APOE-related effects on whole-brain functional network. In our experiments, the risk allele group exhibited lower network radius and modularity in whole brain DMN based on graph theory, suggesting the abnormal organization structure. Moreover, two suggested measures from persistent homology detected significant differences between groups within the left hemisphere and in the whole brain in two datasets. They were more statistically sensitive to APOE genotypic differences than standard graph-based measures. In summary, we provide evidence that the e4 genotype leads to distinct DMN functional alterations in the early phases of Alzheimer’s disease using persistent homology approach. Our study offers a novel insight to explore potential biomarkers in healthy elderly populations carrying APOE e4 allele.

Introduction

Alzheimer,s disease (AD) (Lane et al., 2018) is the most common form of dementia among the elderly and the sixth leading cause of death in the United States. There are more than 50 million patients worldwide in 2018, and it is expected to reach a staggering 152 million by 2050 (Patterson, 2018). It is crucial to develop the AD-related biomarkers early in the aging process before the onset of overt cognitive impairment and irreversible brain damage (Korthauer et al., 2018). One hypothesis (Reiman et al., 2009; Lambert et al., 2013; Yu et al., 2019) for the pathogenesis (Karch and Goate, 2015) of AD indicate the apolipoprotein E (APOE) e4 allele (Lane-Donovan and Herz, 2017) involves the accumulation of Amyloid-β (Caselli et al., 2010), leading to increasing neuronal atrophy and synapse loss. To date, APOE is a major genetic risk factor for developing AD (Thompson et al., 2013; Zhu et al., 2019). Functional neuroimaging genetics provides an effective strategy for characterizing the intermediate phenotype of AD and identifying genes that contribute to functional alterations in brain networks (Chiesa et al., 2017, Chiesa et al., 2019). In particular, recent research has demonstrated that default mode network (DMN) (Raichle, 2015) is associated with progressive brain dysfunction and is susceptible to APOE genotype (Song et al., 2015; Ma et al., 2016; Yuan et al., 2016; Palmqvist et al., 2017; Staffaroni et al., 2018; Chiesa et al., 2019).

Graph theory has increasingly been used as a theoretical framework for studying brain network characteristics. At the level of regional connection, functional connectivity between brain nodes as an important biomarker can identify early brain function alteration related to AD pathophysiology (Bokde et al., 2009). It is dedicated to investigating the distinct connectivity within the DMN that could represent the progressive biomarker. However, the results of resting state functional MRI (rs-fMRI) (Teipel et al., 2015) and APOE studies have reported mixed results (Cai et al., 2017; Luo et al., 2017; Caldwell et al., 2019; Zhu et al., 2019). Some reported decreased functional connectivity (Yan et al., 2015) in APOE e4 allele carriers (APOE4+) compared with non-carries (APOE4-), the others found some increased functional connectivity (Song et al., 2015; Zhu et al., 2018), while others didn’t found any differences (Chiesa et al., 2019). At the global whole-brain level, some neurobiologically meaningful graph-theoretic properties have become important indicators for measuring brain functional networks, through which we can understand the altered network architecture in those carrying risk genotype, including a loss of small-world network (Korthauer et al., 2018), a redistribution of hubs (Wink et al., 2018), and a disrupted modular organization (Li et al., 2019). However, there are currently few network measures based on graph theory have been studied in cognitively unimpaired elderly (Seo et al., 2013; Luo et al., 2017; Pietzuch et al., 2019) and some of their results were reported as inconsistent (Seo et al., 2013; Qiu et al., 2016; Luo et al., 2017). For instance, Wink et al. (2018) found decreased centrality of DMN in APOE4+ comparing to non-carriers, while Wang et al. (2017) didn’t find such genotype difference of centrality in normal elderly. Overall, graph-theoretic methods cannot consistently demonstrate functional DMN difference between APOE4+ and APOE4- in normal elderly, and the reason has been debated in the literature (Chiesa et al., 2017).

Recently, persistent homology (Edelsbrunner and Harer, 2010) from algebraic topology has been adopted for the analysis of brain network. It uses graph filtration to construct a multiscale brain network with all possible thresholds wherever the persistent topological features over the network dynamics are identified (Giusti et al., 2016). This method can distinguish robust and noisy topological features over a wide range of filtration values in measuring global brain network organization. The typical approach of persistent homology is Betti number plot (BNP) (Edelsbrunner and Harer, 2010; Lee et al., 2012), which has successfully applied to the brain network research on epilepsy (Choi et al., 2014), autism spectrum disorder and attention-deficit hyperactivity disorder (Lee et al., 2012, 2017), etc. In our previous works (Kuang et al., 2019a, b), we have developed some network properties based on persistent homology and have successfully applied them to measure the metabolic and functional networks of AD and MCI patients. Although the persistent homology works well in cognitively impaired elderly, it has never been applied to study the genetic influence on brain network yet, especially in unimpaired individuals.

In this paper, we study the effect of APOE genotype on functional DMN in cognitively unimpaired subjects. We hypothesized the topological properties of persistent homology may reveal the APOE-related alteration in DMN even before clinical symptoms appear better than graph-theoretic approaches. Using the cross-sectional rs-fMRI imaging data of 27 APOE4+ and 31 APOE4- normal elderly, we test this hypothesis by computing two persistent homology-based properties and measuring the differences between APOE4+ and APOE4- groups. We further run the statistical inference to validate their powers and compare them with some graph-theoretic methods.

Materials and Methods

We summarize the pipeline of our framework in Figure 1. The rs-fMRI (Teipel et al., 2015) data of each subject are preprocessed and the blood oxygen level dependent (BOLD) signals within each region-of-interest (ROI) are obtained. Then we construct one weighted DMN per subject and quantify its global topological structure using graph theory and persistent homology. The details are described in following subsections.

FIGURE 1
www.frontiersin.org

Figure 1. The flow of measuring DMN topological structure based on graph theory and persistent homology in cognitively unimpaired subjects using rs-fMRI data from ADNI.

Participants

Data used in the preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database1 (Jack et al., 2008; Jagust et al., 2010). The ADNI was launched in 2003 as a public-private partnership, led by Principal Investigator Michael W. Weiner, MD. The primary goal of ADNI has been to test whether serial magnetic resonance imaging (MRI), positron emission tomography (PET), other biological markers, and clinical and neuropsychological assessment can be combined to measure the progression of mild cognitive impairment (MCI) and early Alzheimer’s disease (AD).

There were only 38 Normal Controls (NC) between the ages of 60 and 90 from the ADNI-2 who had available rs-fMRI and APOE data. Due to the small sample size of NC, we also introduced Subjective Memory Complaints (SMC), producing a dataset of cognitively unimpaired subjects in this study. The only difference from NC is that SMC reported memory problems by themselves. Individuals with two copies of the apolipoprotein e2 allele (APOE e2/2) were excluded due to its possible protective effects (Suri et al., 2013). Finally, individuals carrying at least one APOE ε4 allele (genotype e3/e4 and e4/e4) were classified as APOE4+, while individuals with genotype e3/e3 were classified as APOE4-.

Data Acquisition and Preprocessing

The experimental dataset were acquired at multiple ADNI sites using 3.0 T Philips MRI scanners. All rs-fMRI data were obtained using an echo-planar imaging (EPI) sequence and the parameters included repetition time (TR) = 3000 ms, echo time (TE) = 30 ms, flip angle = 80°, number of slices = 48, slice thickness = 3.3 mm, voxel size = 3 mm × 3 mm × 3 mm, voxel matrix = 64 × 64, and time points = 140.

All functional images were pre-processed using SPM8 toolbox2, DPARSF3 (Yan and Zang, 2010), and REST3 (Song et al., 2011) according to well-accepted pipelines, the same as our prior work (Kuang et al., 2019a). Briefly, the first ten time points were removed before temporal correction and spatial normalization. Then image smoothing, linear trend adjustment and band-pass filter were performed sequentially.

Construction of DMN

First, the whole brain is divided into 90 functional ROI using standard automated anatomical labeling atlas (AAL90) (Tzourio-Mazoyer et al., 2002). Then 26 areas (Vriend et al., 2018) in AAL90 are identified as the ROI of DMN, as shown in Table 1. There are 13 ROI per hemisphere and each ROI is considered as a network node of DMN.

TABLE 1
www.frontiersin.org

Table 1. The division of twenty six ROI nodes in DMN based on AAL90 atlas.

The average timing BOLD signal serial Ti = (Ti1, Ti2,…Tik) within the i-th ROI node is used as its measurement (Step 2 in Figure 1). We define the functional connectivity (i.e., edge weight) between any pair of ROI as 1-Pearson coefficient of their BOLD signal serials, i.e.

Wij=1cov(Ti,Tj)σTiσTj=1p=1k(TipT¯i)(TjpT¯j)p=1k(TipT¯i)2p=1k(TjpT¯j)2(1)

where Tip represents the average BOLD signal within the i-th ROI at p-th time point and K = 130 is total number of time points of the rs-fMRI data. Thus, the functional connection matrix (N × N) per subject is obtained (Step 3 in Figure 1) and each subject’s DMN is constructed. Here N = 26 if the DMN of entire brain is studied, otherwise N = 13 if only one hemispheric DMN is studied.

Measuring DMN Using Graph Theory

In the past decade, the neurobiologically meaningful network properties based on graph theory have become important indicators in measuring brain functional networks. We validate some widely used graph measures in this study (Right part of step 5 in Figure 1), including characteristic path length (CPL) (Li et al., 2019), global efficiency (GC) (Shu et al., 2015), network radius (NR) (Fujita et al., 2017), modularity (Mod) (Li et al., 2019), and eigenvector centrality (EC) (Luo et al., 2017). Briefly, the average shortest path length between all pairs of nodes in the network is CPL, while the average inverse shortest path length is called GC. Then NR is the minimum eccentricity of all nodes in the network and nodal eccentricity is the greatest distance between this node and any other nodes. Further, Mod measures the extent to which the network can be subdivided into clearly delineated and non-overlapping groups, and EC computes the sum of centralities of the node’s direct neighbors. All these network measures were calculated by Brain Connectivity Toolbox (BCT)4 in Matlab R2017a.

Measuring DMN Using Persistent Homology

Persistent homology (Edelsbrunner and Harer, 2010) is a mathematical concept derived from algebraic topology and is used to characterize topological features in complex data. There is an important tool, graph filtration (Giusti et al., 2016), in persistent homology that constructs a family of nested networks along an axis at their threshold values by thresholding original weighted network at every possible entry (Step 4 in Figure 1). Thus, it can distinguish robust and noisy topological characteristics in a wide range and enables reasonable inferences regarding the underlying organization. The classic network property based on persistent homology is BNP which detects the dynamic of the zeroth Betti number (i.e., the number of connected components) over all filtration values (Lee et al., 2012). It has been successfully applied to the some studies (Lee et al., 2012, 2017; Choi et al., 2014) of brain network in neurodegenerative diseases.

In our previous work (Kuang et al., 2019a), we proposed an integrated persistent feature (IPF) based on BNP, which introduced a connected component aggregation cost into the zeroth Betti number and thus achieves a holistic description of network dynamics. The IPF at filtration λi is defined as (Kuang et al., 2019a).

IPFλi={m-im(m-1)k=i+1m-1λk0im-20i=m-1(2)

Here, m is total number of network nodes and λ0 = 0 < λ0 < λ1 < λ2 < … < λm–1 is the filtration value which is actually the set of weights of minimum spanning tree of the original weighted network. Previous work has proven that the IPF is a monotonically decreasing convergence function over all possible filtration. In summary, when λ increases from zero, the IPF value of the network will decrease to zero accordingly until all nodes are connected into a single connected component. Therefore, the slope of the IPF plot (SIP) can be used as an important network property to quantify the brain network dynamics. Both network measures BNP and SIP can be considered as information diffusion rate or convergence rate of the network. We provided their implementations at http://gsl.lab.asu.edu/software/IPF and applied them in this study (Left part of step 5 in Figure 1).

Results

Demographic Information

In this experiment, 58 subjects without cognitive impairment were selected from ADNI-2, and were divided into two groups, APOE e4 carriers and non-carriers, according to their APOE genotype. Among them, 20 subjects were identified as SMC and remaining 38 subjects were NC. The only difference from NC is that SMC reported memory problems by themselves. We considered both as cognitively impairment subjects in this study. As shown in Table 2, there were no significant differences in age, education, Mini Mental State Examination (MMSE) score, and Clinical Dementia Rating (CDR) global scores between groups. All subjects had MMSE of 24–30, CDR = 0, and were cognitively unimpaired. Furthermore, all studied images did not have excessive head motion (six-parameter rigid body) defined by a displacement of less than 1 mm or an angular rotation of less than in any direction 1°.

TABLE 2
www.frontiersin.org

Table 2. Demographic characteristics of the high-risk (APOE4+) and low-risk (APOE4-) groups.

Multiscale Brain DMN Dynamics

We constructed one original weighted DMN per subject. The 26 network nodes in whole brain DMN were determined according to Table 2 and are visualized in Figure 2 using Brain Net Viewer software (Xia et al., 2013). Then, the edge weights between them were calculated using Eq. (1). We further constructed multiscale networks based on the original DMN using graph filtration tool (Step 4 in Figure 1). As we only observe the zeroth homology in this study, the filtration value λ is actually the set of weights of minimum spanning tree of DMN. Figure 3 shows the multiscale network dynamics for two mean DMN of two groups over some filtration values. Figure 4 shows the change of zeroth Betti number using single linkage dendrogram (Lee et al., 2012). All the nodes on the left are connected to form the larger component on the right, until finally a fully connected network is constructed. The zeroth Betti number starts at 26, and gradually decreases to 1 while more and more nodes are connected.

FIGURE 2
www.frontiersin.org

Figure 2. The ROI in DMN from (A) sagittal view, (B) axial view, and (C) coronal view. The color bar shows the ROI node index predefined in AAL90 atlas.

FIGURE 3
www.frontiersin.org

Figure 3. Multiscale DMN dynamics for two mean networks of (A) APOE4+ and (B) APOE4- groups at six different filtration values 0.20, 0.25, 0.30, …, 0.45. The color bar shows the ROI node index predefined in AAL90 atlas.

FIGURE 4
www.frontiersin.org

Figure 4. The single linkage dendrograms of (A) APOE4 carriers and (B) non-carriers groups show the change of the zeroth Betti number. The color represents the target distance (total edge weight) from current connected component to the full connected component (the rightmost line).

From Figures 3, 4, we intuitively saw that the connected components in APOE+ aggregated slightly faster than APOE-, especially after λ is larger than 0.35. However, it needs to be further quantitatively measured by network properties based on persistent homology.

Brain DMN Properties

We calculated the corresponding Betti number β0 and IPF of the multiscale DMN at all different filtration values for two group means, and plotted them, as shown in Figure 5. We found that the APOE+ curve in both Betti number plot and IPF plot were steeper than the APOE- curve, suggesting the faster aggregation of APOE+, which is consistent with the above observation of multiscale brain dynamics (see Figures 3, 4). All subjects’ values of BNP and SIP properties based on persistent homology were summarized using box plot as shown in Figures 6A,B, separately, where 1 represents APOE4+ and 2 is APOE4-. The distributions of both BNP and SIP property values between groups are obviously different, indicating both persistent features may be able to discriminate APOE4+ from APOE4-.

FIGURE 5
www.frontiersin.org

Figure 5. The persistent feature dynamics over filtrations for APOE4+ and APOE4- by (A) Betti number plot and (B) IPF plot.

FIGURE 6
www.frontiersin.org

Figure 6. The box plots of property values for APOE4 carriers (1) and non-carriers (2) groups using two persistent homology properties, (A) SIP and (B) BNP, and five graph-theoretic properties, (C) CPL, (D) NR, (E) EC, (F) Mod, and (G) GC.

Traditionally, brain network properties have been measured using graph theory methods. In order to compare with our suggested methods based on persistent homology, we also calculated some classical graph theory properties, including CPL, NR, EC, Mod, and GC. The distributions of all attribute values are shown in Figures 6C–G where 1 and 2 represents APOE4+ and APOE4-, respectively. We observed that the between-group differences of SIP, BNP, NR, and Mod are more apparent than those of CPL, EC, and GC.

Statistical Group Difference

In the statistical analysis of differences between groups of APOE4+ and APOE4-, we performed the permutation test of 10,000 permutations on all network properties using Matlab R2017a and calculated their resulting p-values as shown in Table 3. First, the differences between groups in whole brain DMN with 26 ROI nodes were measured. Two persistent features SIP and BNP obtained significant differences at significance level of 0.05, which were p = 0.021 and p = 0.009, respectively. In the statistic inferences for five compared graph theory-based properties, only NR and Mod obtained significant differences with p = 0.024 and p = 0.037, respectively, while there were no significant differences in other three properties, CPL, EC, and GC. Then, we analyzed the group differences within the single hemisphere and only two measures SIP and BNP achieved significant differences (p = 0.027, both) within the left hemisphere. We did not find any differences within the right hemisphere.

TABLE 3
www.frontiersin.org

Table 3. Statistical p-values of different network properties between APOE4+ and APOE4- groups.

In short, our experimental results show that both persistent properties achieved more significant group differences between APOE4+ and APOE4- than traditional graph-theoretic measures, and BNP obtained the most significant difference (p = 0.009) in the study of whole brain DMN.

Discussion

Present Findings

There are three main findings in this study.

First, we found that the e4 allele carriers exhibited lower NR and Mod (p = 0.024, 0.037, respectively) in the study of whole brain DMN using traditional graph-theoretic methods, suggesting the abnormal organization structure in the risk allele group. To our knowledge, there have been few studies (Luo et al., 2017; Pietzuch et al., 2019) on graph theory that have reported APOE genotypic differences in functional network properties of whole brain DMN in normal elderly, although a lot of studies have found differences in functional network properties between AD/MCI and NC. Some studies (Staffaroni et al., 2018; Chiesa et al., 2019) even found no difference between elderly APOE4+ and APOE4- groups in functional DMN. The reason why findings of DMN on unimpaired individuals do not consistently demonstrate differences between APOE4+ and APOE4- is still debated in the literature (Chiesa et al., 2017). In our study, we found two measures could detect their differences significantly, which would further enhance the APOE research based on graph theory.

Second, we introduced two measures from our previous studies based on persistent homology and found they were more statistically powerful than graph-theoretic measures in discriminating APOE4+ from APOE4- in our experiment, and the BNP obtained the most significant difference (p = 0.009) between groups. The persistent homology approach can distinguish robust and noisy topological features over multiscale nested networks, and the obtained properties are more persistent and stable. So far, many studies of brain network based on persistent homology have demonstrated the superiority of the performance. To our knowledge, it is the first time we introduced persistent homology to study the APOE genotype effect on DMN.

Finally, the functional disruption within the left hemisphere may be more pronounced than the right one. All persistent homology-based features SIP and BNP detected the significant differences of DMN in whole brain and left hemisphere. However, no significant differences were found within the right hemispheric DMN using any network measures. This finding is consistent with existing APOE studies documenting the effect of the e4 allele on left hippocampus rather than right hippocampus of non-demented individuals (Shi et al., 2014; Li et al., 2016; Dong et al., 2019).

Verification on Normal Control Subjects

As there were only 16 NC individuals who had e4 allele and available rs-fMRI data in ADNI-2, we expanded our experimental sample size by including SMC subjects. However, some studies (Caldwell et al., 2019) on APOE have only investigated NC individuals from ADNI. Thus, we excluded all SMC subjects (see Table 2) and further repeated our experiment on NC subjects. As shown in Table 4, there were no significant differences in age, education, sex, and MMSE between APOE+ and APOE- in NC dataset. We calculated the differences of DMN between groups using different measures, as shown in Table 5. Again, two measures based on persistent homology obtained stronger statistical power than graph theory methods, and the differences within the left hemispheric DMN is more significant than the right one. Moreover, compared to experiments performed on cognitively unimpaired dataset (Table 3), the SIP obtained very significant difference (p = 0.006) on NC dataset.

TABLE 4
www.frontiersin.org

Table 4. Demographic characteristics of NC dataset.

TABLE 5
www.frontiersin.org

Table 5. Statistical p-values of different network properties on NC between APOE4+ and APOE4- groups.

Limitation and Future Works

Despite the promising results were obtained by applying two suggested network properties SIP and BNP based on persistent homology to discriminate APOE e4 allele carriers from non-carriers in cognitively unimpaired subjects, there are three important caveats.

First, both persistent homology-based properties BNP and SIP adopted in this study only investigated the dynamics of the zeroth persistent homology. Higher-dimensional persistent homology characterizes higher-dimensional topological features, and more complexed topological structures such as circular holes can be detected. Therefore, the performance of network measurement may be further boosted if higher-dimensional homology is applied, especially in the sparse network that tends to have more holes.

Then, although the DMN have been heavily studied and are reported as a promising kind of network to study, current researches have mainly focused on finding some dysfunctional areas with decreased or increased connectivity (Song et al., 2015; Yan et al., 2015; Zhu et al., 2018; Chiesa et al., 2019). There are relatively few studies on global network dynamics of e4 allele group in normal elderly. In this study, we measured the global brain network in two datasets (cognitively unimpaired dataset and its subset NC) and found some statistically powerful measures. In future, we will validate these measures in other independent datasets.

Finally, current findings are achieved based on cross-sectional study. With longitudinal analysis, we may further study the longitudinal trajectories (Chiesa et al., 2019) of functional brain dynamics and the impact of e4 allele on individuals at risk for Alzheimer’s disease by quantifying the difference of their persistent features. In addition, the more aggregated structure in APOE e4 allele carriers may be due to worse development in childhood. Examining longitudinal MRI since childhood, the possibility of this phenomenon could be further assessed.

Conclusion

This work measured the DMN structure of rs-fMRI on cognitively unimpaired e4 allele carriers based on our prior work of persistent homology, which encodes a great deal of dynamic information over all possible scales. The significant differences between APOE4+ and APOE4- are identified within the left hemispheric DMN and in the whole brain DMN in two datasets, providing evidence that the APOE e4 genotype leads to distinct alterations of functional DMN several years before the occurrence of dementia symptoms. Moreover, our suggested approaches of persistent homology are more sensitive to APOE genotypic differences than standard graph-based network measures. To the best of our knowledge, this is the first study applying persistent homology to analyze APOE-related effect on whole-brain functional network. This study offers a novel insight to explore potential biomarkers in healthy elderly populations carrying APOE e4 allele.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found from the ADNI (http://adni.loni.usc.edu).

Ethics Statement

In this study, all subjects were selected from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (http://adni.loni.usc.edu/). All ADNI subjects gave written informed consent at enrollment for data collection, storage, and use for research.

Author Contributions

LK and YW designed the study and revised the manuscript. FX acquired the data. JJ, DZ, and LK analyzed and interpreted the results of the data. LK and XH drafted the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by the National Key Research and Development Program of China (2018YFB2101504 to LK, FX, and XH), the Shanxi Provincial Key Research and Development Project (201803D121081 and 201903D121147 to LK, FX, and XH), the Natural Science Foundation of Shanxi Province of China (201901D111150 to LK, FX, and XH), the Graduate Science and Technology Project at North University of China (20191638 to JJ), the National Institute on Aging (R21AG043760, R21AG049216, and RF1AG051710 to YW), the National Institute of Biomedical Imaging and Bioengineering (R01EB025032 to YW), the National Heart, Lung, and Blood Institute (R01HL128818 to YW), and the National Science Foundation (DMS-1413417 and IIS-1421165 to YW).

Conflict of Interest

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

Data collection and sharing for this project was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI was funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuroimaging at the University of Southern California.

Footnotes

  1. ^ adni.loni.usc.edu
  2. ^ http://www.fil.ion.ucl.ac.uk/spm/
  3. ^ http://www.restfmri.net
  4. ^ https://sites.google.com/site/bctnet

References

Bokde, A. L. W., Ewers, M., and Hampel, H. (2009). Assessing neuronal networks: understanding Alzheimer’s disease. Prog. Neurobiol. 89, 125–133. doi: 10.1016/j.pneurobio.2009.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, S., Jiang, Y., Wang, Y., Wu, X., Ren, J., Lee, M. S., et al. (2017). Modulation on brain gray matter activity and white matter integrity by APOE epsilon 4 risk gene in cognitively intact elderly: a multimodal neuroimaging study. Behav. Brain Res. 322, 100–109. doi: 10.1016/j.bbr.2017.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Caldwell, J. Z. K., Zhuang, X., Leavitt, M. J., Banks, S. J., Cummings, J., Cordes, D., et al. (2019). Sex moderates amyloid and apolipoprotein epsilon 4 effects on default mode network connectivity at rest. Front. Neurol. 10:900. doi: 10.3389/fneur.2019.00900

PubMed Abstract | CrossRef Full Text | Google Scholar

Caselli, R. J., Walker, D., Sue, L., Sabbagh, M., and Beach, T. (2010). Amyloid load in nondemented brains correlates with APOE e4. Neurosci. Lett. 473, 168–171. doi: 10.1016/j.neulet.2010.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiesa, P. A., Cavedo, E., Lista, S., Thompson, P. M., Hampel, H., and Alzheimer Precision Medicine Initiative, (2017). Revolution of resting-state functional neuroimaging genetics in Alzheimer’s disease. Trends Neurosci. 40, 469–480. doi: 10.1016/j.tins.2017.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiesa, P. A., Cavedo, E., Vergallo, A., Lista, S., Potier, M.-C., Habert, M.-O., et al. (2019). Differential default mode network trajectories in asymptomatic individuals at risk for Alzheimer’s disease. Alzheimers Dementia 15, 940–950. doi: 10.1016/j.jalz.2019.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, H., Kim, Y. K., Kang, H., Lee, H., Im, H.-J., Kim, E. E., et al. (2014). Abnormal metabolic connectivity in the pilocarpine-induced epilepsy rat model: a multiscale network analysis based on persistent homology. Neuroimage 99, 226–236. doi: 10.1016/j.neuroimage.2014.05.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, Q., Zhang, W., Wu, J., Li, B., Schronh, E. H., McMahon, T., et al. (2019). Applying surface-based hippocampal morphometry to study APOE-E4 allele dose effects in cognitively unimpaired subjects. Neuroimage Clin. 22:101744. doi: 10.1016/j.nicl.2019.101744

PubMed Abstract | CrossRef Full Text | Google Scholar

Edelsbrunner, H., and Harer, J. (2010). Computational Topology: an Introduction. Providence, RI: American Mathematical Society.

Google Scholar

Fujita, A., Takahashi, D. Y., Balardin, J. B., Vidal, M. C., and Sato, J. R. (2017). Correlation between graphs with an application to brain network analysis. Comput. Stat. Data Anal. 109, 76–92. doi: 10.1016/j.csda.2016.11.016

CrossRef Full Text | Google Scholar

Giusti, C., Ghrist, R., and Bassett, D. S. (2016). Two’s company, three (or more) is a simplex : algebraic-topological tools for understanding higher-order structure in neural data. J. Comput. Neurosci. 41, 1–14. doi: 10.1007/s10827-016-0608-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Jack, C. R., Bernstein, M. A., Fox, N. C., Thompson, P., Alexander, G., Harvey, D., et al. (2008). The Alzheimer’s disease neuroimaging initiative (ADNI): MRI methods. J. Magnet. Reson. Imaging 27, 685–691.

Google Scholar

Jagust, W. J., Bandy, D., Chen, K., Foster, N. L., Landau, S. M., Mathis, C. A., et al. (2010). The Alzheimer’s disease neuroimaging initiative positron emission tomography core. Alzheimer’s Dementia 6, 221–229.

Google Scholar

Karch, C. M., and Goate, A. M. (2015). Alzheimer’s disease risk genes and mechanisms of disease pathogenesis. Biol. Psychiatry 77, 43–51. doi: 10.1016/j.biopsych.2014.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Korthauer, L. E., Zhan, L., Ajilore, O., Leow, A., and Driscoll, I. (2018). Disrupted topology of the resting state structural connectome in middle-aged APOE epsilon 4 carriers. Neuroimage 178, 295–305. doi: 10.1016/j.neuroimage.2018.05.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuang, L., Han, X., and Chen, K. (2019a). A concise and persistent feature to study brain resting-state network dynamics: findings from the Alzheimer’s disease neuroimaging initiative. Hum. Brain Mapp. 40, 1062–1081. doi: 10.1002/hbm.24383

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuang, L., Zhao, D., Xing, J., Chen, Z., Xiong, F., and Han, X. (2019b). Metabolic brain network analysis of FDG-PET in Alzheimer’s disease using kernel-based persistent features. Molecules 24:2301. doi: 10.3390/molecules24122301

PubMed Abstract | CrossRef Full Text | Google Scholar

Lambert, J.-C., Ibrahim-Verbaas, C. A., Harold, D., Naj, A. C., Sims, R., Bellenguez, C., et al. (2013). Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat. Genet. 45, 1452–U1206. doi: 10.1038/ng.2802

PubMed Abstract | CrossRef Full Text | Google Scholar

Lane, C. A., Hardy, J., and Schott, J. M. (2018). Alzheimer’s disease. Eur. J. Neurol. 25, 59–70. doi: 10.1111/ene.13439

PubMed Abstract | CrossRef Full Text | Google Scholar

Lane-Donovan, C., and Herz, J. (2017). ApoE, ApoE receptors, and the synapse in Alzheimer’s disease. Trends Endocrinol. Metab. 28, 273–284. doi: 10.1016/j.tem.2016.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Kang, H., Chung, M. K., Kim, B.-N., and Lee, D. S. (2012). Persistent brain network homology from the perspective of dendrogram. IEEE Trans. Med. Imaging 31, 2267–2277. doi: 10.1109/tmi.2012.2219590

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Kang, H., Chung, M. K., Lim, S., Kim, B. N., and Lee, D. S. (2017). Integrated multimodal network approach to PET and MRI based on multidimensional persistent homology. Hum. Brain Mapp. 38, 1387–1402. doi: 10.1002/hbm.23461

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Shi, J., Gutman, B. A., Baxter, L. C., Thompson, P. M., Caselli, R. J., et al. (2016). Influence of APOE genotype on hippocampal atrophy over time – an N=1925 surface-based ADNI study. PLoS One 11:e0152901. doi: 10.1371/journal.pone.0152901

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Yao, Z. J., Yu, Y., Fu, Y., Zou, Y., Hu, B., et al. (2019). The influence of Cerebrospinal fluid abnormalities and APOE 4 on PHF-Tau protein: evidence from voxel analysis and graph theory. Front. Aging Neurosci. 11:208. doi: 10.3389/fnagi.2019.00208

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, X., Qiu, T., Jia, Y., Huang, P., Xu, X., Yu, X., et al. (2017). Intrinsic functional connectivity alterations in cognitively intact elderly APOE epsilon 4 carriers measured by eigenvector centrality mapping are related to cognition and CSF biomarkers: a preliminary study. Brain Imaging Behav. 11, 1290–1301. doi: 10.1007/s11682-016-9600-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, C., Zhang, Y., Li, X., Chen, Y., Zhang, J., Liu, Z., et al. (2016). The TT allele of rs405509 synergizes with APOE epsilon 4 in the impairment of cognition and its underlying default mode network in non-demented elderly. Curr. Alzheimer Res. 13, 708–717. doi: 10.2174/1567205013666160129100350

PubMed Abstract | CrossRef Full Text | Google Scholar

Palmqvist, S., Scholl, M., Strandberg, O., Mattsson, N., Stomrud, E., Zetterberg, H., et al. (2017). Earliest accumulation of beta-amyloid occurs within the default-mode network and concurrently affects brain connectivity. Nat. Commun. 8:1214. doi: 10.1038/s41467-017-01150-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Patterson, C. (2018). World Alzheimer Report 2018. London: Alzheimer’s Disease International.

Google Scholar

Pietzuch, M., King, A. E., Ward, D. D., and Vickers, J. C. (2019). The influence of genetic factors and cognitive reserve on structural and functional resting-state brain networks in aging and Alzheimer’s disease. Front. Aging Neurosci. 11:30. doi: 10.3389/fnagi.2019.00030

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, T., Luo, X., Shen, Z., Huang, P., Xu, X., Zhou, J., et al. (2016). Disrupted brain network in progressive mild cognitive impairment measured by eigenvector centrality mapping is linked to cognition and cerebrospinal fluid biomarkers. J. Alzheimer’s Dis. 54, 1483–1493. doi: 10.3233/JAD-160403

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Reiman, E. M., Chen, K., Liu, X., Bandy, D., Yu, M., Lee, W., et al. (2009). Fibrillar amyloid-beta burden in cognitively normal people at 3 levels of genetic risk for Alzheimer’s disease. Proc. Natl. Acad. Sci. U.S.A. 106, 6820–6825. doi: 10.1073/pnas.0900345106

PubMed Abstract | CrossRef Full Text | Google Scholar

Seo, E. H., Lee, D. Y., Lee, J.-M., Park, J.-S., Sohn, B. K., Choe, Y. M., et al. (2013). Influence ofAPOE genotype on whole-brain functional networks in cognitively normal elderly. PLoS One 8:e83205. doi: 10.1371/journal.pone.0083205

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, J., Lepore, N., Gutman, B. A., Thompson, P. M., Baxter, L. C., Caselli, R. J., et al. (2014). Genetic influence of apolipoprotein E4 Genotype on hippocampal morphometry: an N=725 surface-based Alzheimer’s disease neuroimaging initiative study. Hum. Brain Mapp. 35, 3903–3918. doi: 10.1002/hbm.22447

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, N., Li, X., Ma, C., Zhang, J., Chen, K., Liang, Y., et al. (2015). Effects of APOE promoter polymorphism on the topological organization of brain structural connectome in nondemented elderly. Hum. Brain Mapp. 36, 4847–4858. doi: 10.1002/hbm.22954

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, H., Long, H., Zuo, X., Yu, C., Liu, B., Wang, Z., et al. (2015). APOE effects on default mode network in Chinese cognitive normal elderly: relationship with clinical cognitive performance. PLoS One 10:e0133179. doi: 10.1371/journal.pone.0133179

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, X.-W., Dong, Z.-Y., Long, X.-Y., Li, S.-F., Zuo, X.-N., Zhu, C.-Z., et al. (2011). REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS One 6:e25031. doi: 10.1371/journal.pone.0025031

PubMed Abstract | CrossRef Full Text | Google Scholar

Staffaroni, A. M., Brown, J. A., Casaletto, K. B., Elahi, F. M., Deng, J., Neuhaus, J., et al. (2018). The longitudinal trajectory of default mode network connectivity in healthy older adults varies as a function of age and is associated with changes in episodic memory and processing speed. J. Neurosci. 38, 2809–2817. doi: 10.1523/jneurosci.3067-17.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Suri, S., Heise, V., Trachtenberg, A. J., and Mackay, C. E. (2013). The forgotten APOE allele: a review of the evidence and suggested mechanisms for the protective effect of APOE epsilon 2. Neurosci. Biobehav. Rev. 37, 2878–2886. doi: 10.1016/j.neubiorev.2013.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Teipel, S., Drzezga, A., Grothe, M. J., Barthel, H., Chételat, G., Schuff, N., et al. (2015). Multimodal imaging in Alzheimer’s disease: validity and usefulness for early detection. Lancet Neurol. 14, 1037–1053. doi: 10.1016/s1474-4422(15)00093-9

CrossRef Full Text | Google Scholar

Thompson, P. M., Ge, T., Glahn, D. C., Jahanshad, N., and Nichols, T. E. (2013). Genetics of the connectome. Neuroimage 80, 475–488. doi: 10.1016/j.neuroimage.2013.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Vriend, C., van den Heuvel, O. A., Berendse, H. W., van der Werf, Y. D., and Douw, L. (2018). Global and subnetwork changes of the structural connectome in de novo Parkinson’s disease. Neuroscience 386, 295–308. doi: 10.1016/j.neuroscience.2018.06.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Dai, Z., Shu, H., Liao, X., Yue, C., Liu, D., et al. (2017). APOE genotype effects on intrinsic brain network connectivity in patients with amnestic mild cognitive impairment. Sci. Rep. 7:397. doi: 10.1038/s41598-017-00432-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wink, A. M., Tijms, B. M., ten Kate, M., Raspor, E., de Munck, J. C., Altena, E., et al. (2018). Functional brain network centrality is related to APOE genotype in cognitively normal elderly. Brain Behav. 8:e01080. doi: 10.1002/brb3.1080

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, M., Wang, J., and He, Y. (2013). BrainNet viewer: a network visualization tool for human brain connectomics. PLoS One 8:e68910. doi: 10.1371/journal.pone.0068910

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C., and Zang, Y. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, F.-X., Wu, C. W., Chao, Y.-P., Chen, C.-J., and Tseng, Y.-C. (2015). APOE-epsilon 4 allele altered the rest-stimulus interactions in healthy middle-aged adults. PLoS One 10:e0128442. doi: 10.1371/journal.pone.0128442

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, L., Petyuk, V. A., Tasaki, S., Boyle, P. A., Gaiteri, C., Schneider, J. A., et al. (2019). Association of cortical beta-amyloid protein in the absence of insoluble deposits with Alzheimer disease. JAMA Neurol. 76, 818–826. doi: 10.1001/jamaneurol.2019.0834

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, B., Xie, C., Shu, H., Liao, W., Wang, Z., Liu, D., et al. (2016). Differential effects of APOE genotypes on the anterior and posterior subnetworks of default mode network in amnestic mild cognitive impairment. J. Alzheimers Dis. 54, 1409–1423. doi: 10.3233/jad-160353

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, L., Shu, H., Liu, D., Guo, Q., Wang, Z., and Zhang, Z. (2018). Apolipoprotein E epsilon 4 specifically modulates the hippocampus functional connectivity network in patients with amnestic mild cognitive impairment. Front. Aging Neurosci. 10:289. doi: 10.3389/fnagi.2018.00289

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Gong, L., He, C., Wang, Q., Ren, Q., Xie, C., et al. (2019). Default mode network connectivity moderates the relationship between the APOE genotype and cognition and individualizes identification across the Alzheimer’s disease spectrum. J. Alzheimers Dis. 70, 843–860. doi: 10.3233/jad-190254

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: APOE, Alzheimer’s disease, persistent homology, resting state functional magnetic resonance imaging, graph theory, network measure

Citation: Kuang L, Jia J, Zhao D, Xiong F, Han X and Wang Y for the Alzheimer’s Disease Neuroimaging Initiative (2020) Default Mode Network Analysis of APOE Genotype in Cognitively Unimpaired Subjects Based on Persistent Homology. Front. Aging Neurosci. 12:188. doi: 10.3389/fnagi.2020.00188

Received: 26 March 2020; Accepted: 02 June 2020;
Published: 30 June 2020.

Edited by:

Jiehui Jiang, Shanghai University, China

Reviewed by:

Takahito Yoshizaki, Keio University, Japan
Dafin F. Muresanu, Iuliu Haţieganu University of Medicine and Pharmacy, Romania

Copyright © 2020 Kuang, Jia, Zhao, Xiong, Han and Wang for the Alzheimer’s Disease Neuroimaging Initiative. 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: Liqun Kuang, kuang@nuc.edu.cn; Yalin Wang, ylwang@asu.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.