Share

October 2018; 4 (5) ArticleOpen Access

Protein network analysis reveals selectively vulnerable regions and biological processes in FTD

Luke W. Bonham, Natasha Z.R. Steele, Celeste M. Karch, Claudia Manzoni, Ethan G. Geier, Natalie Wen, Aaron Ofori-Kuragu, Parastoo Momeni, John Hardy, Zachary A. Miller, Christopher P. Hess, Patrick Lewis, Bruce L. Miller, William W. Seeley, Sergio E. Baranzini, Rahul S. Desikan, Raffaele Ferrari, Jennifer S. Yokoyama, on behalf of the International FTD-Genomics Consortium (IFGC)
First published October 1, 2018, DOI: https://doi.org/10.1212/NXG.0000000000000266
Full PDF
Citation
Permissions

Make Comment

See Comments

Downloads
458

Share

Abstract

Objective The neuroanatomical profile of behavioral variant frontotemporal dementia (bvFTD) suggests a common biological etiology of disease despite disparate pathologic causes; we investigated the genetic underpinnings of this selective regional vulnerability to identify new risk factors for bvFTD.

Methods We used recently developed analytical techniques designed to address the limitations of genome-wide association studies to generate a protein interaction network of 63 bvFTD risk genes. We characterized this network using gene expression data from healthy and diseased human brain tissue, evaluating regional network expression patterns across the lifespan as well as the cell types and biological processes most affected in bvFTD.

Results We found that bvFTD network genes show enriched expression across the human lifespan in vulnerable neuronal populations, are implicated in cell signaling, cell cycle, immune function, and development, and are differentially expressed in pathologically confirmed frontotemporal lobar degeneration cases. Five of the genes highlighted by our differential expression analyses, BAIAP2, ERBB3, POU2F2, SMARCA2, and CDC37, appear to be novel bvFTD risk loci.

Conclusions Our findings suggest that the cumulative burden of common genetic variation in an interacting protein network expressed in specific brain regions across the lifespan may influence susceptibility to bvFTD.

Glossary

AD=
Alzheimer disease;
BP=
biological process;
CSEA=
cell-type SEA;
EGFR=
epidermal growth factor;
FDR=
false-discovery rate;
bvFTD=
behavioral variant frontotemporal dementia;
FTLD=
frontotemporal lobar degeneration;
GO=
gene ontology;
GWAS=
genome-wide association study;
IFGC=
International FTD-Genomics Consortium;
MAPK=
mitogen-activated protein kinase;
PD=
Parkinson disease;
PDGF=
platelet-derived growth factor;
PINBPA=
protein interaction network-based pathway analysis;
SEA=
specific enrichment analysis;
VEN=
Von Economo neuron;
W-PPI-NA=
weighted protein-protein interaction network analysis

Although defined by a clear clinical syndrome, behavioral variant frontotemporal dementia (bvFTD) is the most pathologically diverse form of frontotemporal lobar degeneration (FTLD). Despite variable pathology, we and others have shown that a shared, selectively vulnerable brain network, composed of disparate brain regions functionally connected throughout the lifespan, degenerates in bvFTD.1 The fact that different protein pathologies converge on a relatively consistent set of neuroanatomical regions suggests intrinsic molecular properties may predispose this network to disease. Nevertheless, despite extensive advances in our ability to clinically and pathologically diagnose bvFTD, little is known about the protein networks (“nexopathies”) that drive the biological processes (BPs) underlying this selective vulnerability.2,3

New techniques enable evaluation of experimentally determined protein networks in disease by aggregating single-nucleotide polymorphism (SNP)-level risk metrics across an entire gene. This technique increases statistical power, overcomes locus heterogeneity likely to occur in clinical populations, and allows for better detection of multiple variants contributing to polygenic disease risk.4 Combining this information with existing protein interaction data has revealed new genetic risk loci and helped unravel the pathophysiology of complex diseases like multiple sclerosis.5

In this study, we first generated networks of genes underlying bvFTD risk to identify common pathophysiologic processes underlying disease biology. Second, based on prior studies showing that the brain network affected in bvFTD is present and functioning throughout life,6 we evaluated spatial and temporal risk gene expression patterns in human brain across the lifespan. Finally, we sought to validate risk network genes using brain tissue from patients pathologically diagnosed with FTLD.

Methods

Standard protocol approvals, registrations, and patient consents

Written informed consent was obtained from all participants or their guardians at the site of sample collection, and relevant institutional review boards approved all aspects of the study.

Genome-wide association study study participants

This study used SNP level metadata from phase 1 FTD-GWAS data from the International FTD-Genomics Consortium (IFGC), which was composed of 1,377 bvFTD cases and 4,308 controls with genotypes or imputed data at 6,026,384 SNPs. Raw genetic data for each participant was not available for analysis at the inception of this study. Cases and controls within the cohort were diagnosed according to published criteria for bvFTD7 by a trained neurologist or pathologic examination.8

Gene network generation

We calculated gene-level significance values using Versatile Gene-based Association Study (VEGAS)9 and summary data from a recent bvFTD genome-wide association study (GWAS).8 VEGAS assigns SNPs to their respective autosomal genes using their position on the University of California, Santa Cruz Genome Browser. To capture the effects of regulatory regions and SNPs in linkage disequilibrium,10 gene boundaries were defined as 50 kb beyond the 5′ and 3′ untranslated regions of each gene. Empirical p-values for each gene were calculated using Monte-Carlo simulations.5,9 Association blocks were defined as groups of sequential genes with p < 0.05.

We used protein interaction network-based pathway analysis (PINBPA) to compute first-order interactions by filtering a reference network containing 8,960 proteins and 27,724 interactions so that only the genes (and their protein products) with VEGAS p < 0.05 were retained as significant. We evaluated network strength using simulations, which assigned p-values to genes at random from the parent network to create a simulated distribution of similarly sized networks. The empirical bvFTD network's characteristics were compared against the simulated networks.

To increase resolution on the candidate proteins isolated by the PINBPA method, we generated a second independent network with recently developed weighted protein-protein interaction network analysis (W-PPI-NA)11 using the 63 bvFTD genes prioritized by PINBPA as seeds. Of note, the genes prioritized by the PINBPA analysis (CALM2 and CALM3) were combined into one unique entry, CALM1, through this step. We removed TrEMBL, non-protein interactors (e.g., chemicals), obsolete Entrez, and Entrez matching to multiple Swiss-Prot identifiers.

Profiling of gene expression during development and across the human lifespan

To assess whether genes from the PINBPA networks showed enriched expression in human brain tissues of all ages and life stages, we used the R package ABAEnrichment and specific enrichment analysis (SEA). ABAEnrichment performs ontology gene set enrichment analyses for 26 brain regions ranging in age from 8 postconception weeks to 40 years old from the BrainSpan data set distributed by the Allen Brain Institute. Forty-five percent of participants were Caucasian, 33% were African American, 10% were Hispanic, 5% were mixed Caucasian/African American, and 5% had no ethnicity information available. All participants were screened for neurologic conditions and for large-scale chromosomal abnormalities. Samples were grouped by age into 5 groups: prenatal (all prenatal stages; n = 20), infant (0–2 years; n = 6), child (3–11 years; n = 6), adolescent (12–19 years; n = 4), and adult (older than 19 years; n = 6). Additional demographic data are provided in table e-1 (links.lww.com/NXG/A72). All genes evaluated in enrichment analyses were expressed above the 10th percentile when ranked by expression. SEA allows for expression testing across 6 brain regions and 10 different developmental periods. Fisher's exact test was used to assess statistical enrichment, and false-discovery rate (FDR) correction was applied.

Cell-type-specific gene expression profiling

We used 2 RNA expression data sets examining several cell types commonly found in the CNS, one to conduct descriptive representations of the cells that most highly express our genes of interest and the other to specifically test for enrichment.

The first data set used gene expression from human temporal lobe samples.12 For each gene implicated in the bvFTD network, we queried the database and recorded the cell type (fetal astrocytes, mature astrocytes, neurons, oligodendrocytes, microglia/macrophages, and endothelial cells) that most highly expressed the gene of interest. The data set represents expression across the human lifespan (18 postconception weeks [pcw] for fetal samples and 8–63 years old for non-fetal samples).12

The second data set used RNAseq data from the BrainSpan atlas to characterize cell-specific expression patterns. Enrichment for specific cell types was tested using cell-type SEA (CSEA),13 which employs expression cutoffs and tests for cell-type-specific enrichment (CSEA employs methodology similar to SEA). Cell types are determined using well-validated markers,14 and enrichment is tested using Fisher's exact test with FDR correction.

Biological pathways and processes analysis

We classified the biological function of bvFTD network genes using PANTHER (pantherdb.org) and Reactome (reactome.org) to reduce the likelihood of spurious findings from any one source. Statistical tests for overrepresentation with Bonferroni correction were performed using the PANTHER pathways annotation and Gene Ontology (GO) BPs databases. We applied functional annotation analysis to the network built using the 62 PINBPA-prioritized genes as seeds by performing GO BP enrichment analyses through g:Profiler (g:GOSt, biit.cs.ut.ee/gprofiler/15) (62 genes were used because CALM2 and CALM3 were combined into one entry as previously described). A Fisher's one-tailed test was used to assess enrichment; the set counts and sizes threshold was applied as a multiple test correction; statistical domain size was only annotated genes; and no hierarchical filtering was included. We then grouped enriched GO-BP terms into custom-made “semantic classes.” We removed general, thus negligible, semantic classes, such as general, metabolism, enzymes, protein modification, and physiology. Semantic classes were further grouped by similarity in more general classes called functional blocks.

Gene expression in pathologically confirmed FTLD cases

We analyzed gene expression data from 2 cohorts of pathologically confirmed FTLD along with Alzheimer disease (AD) and Parkinson disease (PD) cases compared with normal older adult controls. The first FTLD cohort contained 10 individuals with pathologically confirmed sporadic FTLD, 7 individuals with pathologically confirmed FTLD due to GRN mutations, and 11 pathologically confirmed controls (GSE13162). The second contained 6 FTLD cases with primarily tau pathology (Constantinides group A), 6 FTLD cases with primarily temporal atrophy and minimal tau pathology (Constantinides group C1), and 5 normal controls (E-MEXP-2280). The AD cohort included temporal cortex samples from 84 AD and 80 control brains (syn5550404).16 The PD cohort included substantia nigra samples from 22 PD and 23 control brains (GSE7621).17 In linear regression analyses, we analyzed FTLD cases together in a combined analysis and then performed separate subgroup analyses. Following this, we analyzed the AD and PD groups separately. In all analyses, we covaried for age, sex, and postmortem interval.

Results

bvFTD GWAS data reveals a large protein interaction network driven by risk-associated alleles

Using gene-based p-values for 17,463 genes derived from the bvFTD GWAS cohort, we identified 475 association blocks containing a total of 1,104 genes. We next generated PINBPA networks for the genes implicated in bvFTD. The bvFTD network contained 63 nodes and 72 edges. Permutation testing revealed that the bvFTD gene network was within the top 10th percentile for nodes and edges as compared with a sample of 1,000 randomly generated networks (figure 1 and table e-2, links.lww.com/NXG/A73).

Figure 1
Figure 1 bvFTD PINBPA network

The bvFTD PINBPA network is shown. Nodes are color coded according to their gene-based p-value, with warmer colors indicating lower p-values and cooler colors representing p-values closer to 0.05. Node size represents closeness centrality, a measure of a node's nearness to other nodes within a network. Edge thickness represents edge betweenness, a measure of the number of paths that go through an edge in a network.

We then sought to analyze the topological details of the 63 PINBPA-prioritized genes by generating protein interaction networks with our W-PPI-NA to highlight key functional players within the identified gene set. After filtering and scoring the protein network, the result was composed of 1,913 nodes and 3,212 edges, where all but one node, PBX2, were interconnected (figure 2, A). We identified the inter-interactome hubs (IIHs; n = 10)—the proteins characterized by the highest connection degree distribution (figure 2, B and C)—and we used them to define the core of the network with the highest interconnectivity (figure 2, D). By comparing the core of the network with randomly sampled parts of the network, we verified that the IIHs-driven network is the most densely connected. The core of the network was made of 52 nodes (10 IIHs and their interactors) and 121 edges; these were strongly interconnected (average number of neighbors = 4.5), representing the proteins that keep the cohesion of >16% of the initial seeds (n = 63; figure 2, C).

Figure 2
Figure 2 Interactome analyses in bvFTD

(A) Protein network built through the W-PPI-NA pipeline around the 63 genes prioritized by the PINBPA protocol. The seeds are depicted in pink and their interactors in blue. (B) Inter-interactome degree distribution curve reporting the number of nodes (x-axis) able to bridge an N number of seeds (y-axis). The IIHs are the nodes marked by the rectangle. (C) List of IIHs reporting the number of seeds that they bridge and associated %. *UBC has been reported yet ignored as it is likely a false positive (as it may indicate unspecific binding of ubiquitin to proteins tagged for degradation). (D) Core of the network around the IIHs, which are depicted in pink.

bvFTD network genes show enriched expression in bvFTD-relevant brain regions across the lifespan

We next tested whether genes within the bvFTD network showed enriched or specific expression in distinct brain regions of healthy individuals using ABAEnrichment. Hypergeometric tests revealed enriched expression for our bvFTD network genes in neuroanatomical regions affected in bvFTD, including the dorsolateral prefrontal cortex, frontal cortex, and temporal cortex in most age groups (family-wise error corrected p < 0.05 for each region, figure 3, A–D and table e-3, links.lww.com/NXG/A74).

Figure 3
Figure 3 Regional gene expression enrichment across the lifespan

For each of the 5 age groups used in the ABAEnrichment analyses, we counted the number of times each of the 26 available brain regions showed enrichment of the bvFTD network genes. The superior temporal cortex was the most common region of enriched expression with 5 observations (i.e., in every age grouping). Generally, the FTD cohort showed enriched expression in both frontal (DFC, VFC, M1C, PFC, OFC) and temporal (A1C, ITC, STC, MFC) regions. (A and B) The number of times each region was associated is illustrated on a model brain. Brain regions were mapped onto the illustration using information provided in Bahl et al. 2017. When regional overlap was detected (e.g., IPC is contained within PCx or DFC; VFC and OFC are contained within PFC), the more specific region(s) was chosen for presentation. (C) A graphical depiction of the data shown in (A and B). The number of times each brain region was associated from the 5 age groupings is shown. (D) Detailed results are shown in a table format with age groupings as rows and brain regions as columns. Regions shaded in black were statistically associated while unshaded regions were not. A1C = primary auditory cortex; AMY = amygdaloid complex; CN = cerebral nuclei; CB = cerebellum; CBC = cerebellar cortex; DFC = dorsolateral prefrontal cortex; HIP = hippocampus; IPC = posteroventral (inferior) parietal cortex; ITC = inferolateral temporal cortex; M1C = primary motor cortex; MFC = anterior (rostral) cingulate (medial prefrontal) cortex; OFC = orbital frontal cortex; PCx = parietal neocortex; PFC = prefrontal cortex; S1C = primary somatosensory cortex; STC = posterior (caudal) superior temporal cortex; STR = striatum; TCx = temporal neocortex; THA = thalamus; V1C = primary visual cortex; VFC = ventrolateral prefrontal cortex.

Our ABAEnrichment analyses suggested that enrichment was most common during early life (table e-3, links.lww.com/NXG/A74). However, they did not provide information on the role of gene expression within the early life stages when gene expression patterns can vary temporally and influence neurodevelopmental processes.18 Given this, we performed exploratory analyses using SEA and found that there were several nominally associated regions of gene expression enrichment concentrated within the fetal period. At a specificity index probability level of 0.05, the amygdala (late fetal; p = 0.004), thalamus (early fetal; p = 0.006), hippocampus (late fetal; p = 0.007), and cortex (mid-late fetal; p = 0.014) showed nominal enrichment levels but did not withstand FDR correction (table e-4, links.lww.com/NXG/A75).

bvFTD network genes are most highly expressed in fetal astrocytes, mature astrocytes, microglia, and cortical layer 5b neurons

Analysis of whole transcriptome data from specific cell types isolated from human brain tissue12,19 revealed that the bvFTD network genes clustered into 3 main brain cell groupings. The most commonly observed cell type was fetal astrocytes (37%), followed closely by neurons (35%), and microglia (20%). bvFTD network genes were present at a lower frequency in mature astrocytes (5%), oligodendrocytes (5%), and endothelial cells (5%) (figure e-1, links.lww.com/NXG/A79). Network genes represented only 3 of 100 and 11 of 500 of the top-expressed genes in each cell type, suggesting that the cell type enrichment was specific to bvFTD and not a function of detecting genes that were most highly expressed in these cell types. IL2 and SLA2 were not expressed at a level that allowed for reliable cell type determination.

Given our findings in the cell-type-specific expression data set, we next used CSEA to test for enrichment of our bvFTD network genes in specific cell types and populations. The bvFTD network genes demonstrated enriched expression in layer 5b neuronal cells (p = 0.038; figure e-1, links.lww.com/NXG/A79 and table e-5, links.lww.com/NXG/A76); however, this association did not withstand FDR correction. Von Economo neurons (VENs) and fork cells, which are among the earliest cells to degenerate in bvFTD, are located within layer 5b.20

Ontological analyses reveal overrepresentation of immune signaling in bvFTD risk genes

Functional (ontological) characterization of the 63 bvFTD network genes indicated overrepresentation of immune-mediated pathways, cell signaling, cell cycle, and development (table 1). Reactome analysis indicated the greatest degree of overrepresentation in pathways known to be involved with microglia-initiated inflammatory responses, such as epidermal growth factor (EGFR) signaling, mitogen-activated protein kinases (MAPKs), and pro-inflammatory cytokines (IL-3, IL-5, granulocyte-macrophage colony-stimulating factor [GM-CSF]) (table e-6, links.lww.com/NXG/A77). PANTHER analysis found similar overrepresentation for interleukin signaling pathways, EGFR signaling, and platelet-derived growth factor (PDGF) signaling.

View this table:
Table 1

Pathway analysis in bvFTD network genes

We evaluated additional GO BPs via g:Profiler for the second protein network built by W-PPI-NA and focused on the interactomes of the 63 genes prioritized through the PINBPA protocol. This first functional enrichment was followed by a second iteration performed only on the core of the network—i.e., defined by the IIHs and their interactors—thus containing the most densely connected proteins. Functional annotation of the latter indicated a list of semantic classes that were a subset of the former: of interest, the subset terms (percentage of retention >15%) pointed to the following functional blocks: (1) “development”; (2) “motility”; (3) “protein modification”; and (4) “cell signaling” (figure e-2, links.lww.com/NXG/A80). These appeared as the common functions characterizing the part of the protein network with strongest cohesion among the initial seeds. In particular, the functional block “development” contained semantic classes referring to cell differentiation and proliferation as well as glial and neuronal cell development, while the functional block “cell signaling” contained semantic classes referring to the previously mentioned MAPK, PDGF, and immune-specific signaling in addition to terms related to Fc (typically relevant for immune cells), ERBB (relevant for both developmental and immune processes), growth factors, and hormone-based signaling events. These findings confirmed the PANTHER and Reactome results, using a completely independent approach.

bvFTD network genes are dysregulated in postmortem brain tissue

To examine whether the risk-associated genes identified by the network analysis are differentially expressed in human disease, we assessed the expression patterns of our bvFTD genes in postmortem samples from 2 pathologically confirmed cohorts of FTLD cases along with AD16 and PD17 cases compared with pathology-free controls. In the first FTLD cohort, GSE13162, 61 bvFTD genes had expression data available. In the second FTLD cohort, E-MEXP-2280, there were 58 bvFTD genes with expression data available. Because the 2 data sets represent a mix of sporadic and autosomal dominant (GRN) postmortem cases (GSD13162) or sporadic cases with tau and non-tau pathology (E-MEXP-2280), we analyzed data in aggregate and by gene status or pathologic subtype (table e-7, links.lww.com/NXG/A78). In GSE13162, 19 bvFTD genes were differentially expressed between cases vs controls (PFDR < 0.05). In E-MEXP-2280, 26 bvFTD genes were differentially expressed (PFDR < 0.05). Of these top genes, 9 (BAIAP2, CALD1, CDC37, ERBB3, GSK3A, MAP3K5, POU2F2, SMARCA2, and TGFB1|1) showed dysregulation in both cohorts (table 2 and table e-7). Several of these genes appear to be novel risk factors for bvFTD (table 2). In AD cases, 7 of 9 bvFTD network genes were differentially expressed (table 2). In contrast, 0 of 9 network genes were differentially expressed in PD cases (table 2).

View this table:
Table 2

Differentially expressed genes enriched in both FTLD pathologic cohorts

Discussion

We identified a network of interacting bvFTD risk genes that demonstrate enriched expression in bvFTD-affected brain regions and cell types across the human lifespan. Differential expression analyses in pathologically confirmed cases of FTLD and controls showed that many of our network genes were dysregulated in pathologically confirmed FTLD, and 5 potentially novel bvFTD risk genes showed altered expression in both pathologic cohorts. Taken together, our results provide evidence that the pathologic heterogeneity seen in bvFTD may occur as the product of multiple dysregulated cellular pathways that converge on specific brain networks across the lifespan. Our findings also contribute to the molecular understanding of selective vulnerability in bvFTD and suggest several potential mechanisms by which genetic variation and gene expression abnormalities increase susceptibility to neurodegenerative diseases like bvFTD. Our analyses suggest that this bvFTD risk gene network may overlap with AD, a related disorder sharing genetic risk factors23 and tau pathology24 that can also present with a bvFTD-like syndrome.25 Additional analyses will be useful to elucidate the specificity of these findings to bvFTD compared with AD and determine any shared mechanisms conferring risk across other neurodegenerative diseases.

By integrating GWAS with protein interaction data, gene expression data from across the lifespan, and postmortem gene expression data from confirmed FTLD cases, we highlight the importance of bvFTD risk genes during development, aging, and disease. Our network genes showed enrichment in regions of the brain commonly affected in bvFTD such as the dorsolateral and ventrolateral prefrontal cortices, and inferior temporal cortex (figure 3). Because of the scarcity of pathologic samples available for enrichment analyses, some regions of the brain that have known roles in bvFTD, such as the cingulate and insula, were not available for analysis. Of interest, our genes also showed enrichment in brain regions typically spared in bvFTD patients such as superior temporal cortex, occipital cortex, and sensory cortex. Although the reason for this enrichment remains to be determined, it is possible that spared regions contain similar cell types as vulnerable ones but remain unaffected because they are functionally disconnected from pathologically affected regions in bvFTD and remain unexposed to pathologic aggregates that may disseminate through functionally connected neural networks.26,27 Alternatively, and perhaps counterintuitively, consistently enriched expression (and thereby greater expression across the lifespan) of bvFTD network genes in regions like cerebellum and superior temporal cortex could potentially provide protection against neurodegenerative processes. In this framework, reduced expression of the same genes in other, less-enriched, brain regions may result in vulnerability to disease pathology because less-enriched brain regions have less “reserve” expression to rely on during aging, potentially reducing the threshold within which variable expression is tolerated by relevant biological systems.

An exploratory analysis revealed that our bvFTD network genes are most highly expressed in fetal astrocytes, neurons, and microglia/macrophages. The importance of fetal astrocytes is particularly interesting in the context of our pathway findings as they are implicated in neurodevelopmental processes, complement signaling, and forming healthy neuronal synapses.28 Many of the identified bvFTD risk genes were most highly expressed in neuronal populations, which participate in many of the cell signaling pathways highlighted in our GO analyses.29 Neurons may also be the target of abnormal expression of cell-cycle regulation proteins, which have been shown in studies of neurodegenerative diseases to result in apoptosis when activated.30,31 The third major cell type expressing our genes of interest—microglia/macrophages—play an increasingly appreciated role in the pathogenesis of bvFTD and other neurodegenerative disorders.32 Remarkably, our CSEA analyses highlighted layer 5b cortical neurons as sites of bvFTD network gene expression enrichment. Pathologic evidence in human brain tissue suggests this cortical layer as a likely site of selectively vulnerable neuronal populations in bvFTD (VENs and fork cells).33 While this finding was not significant after correction for multiple testing, it is striking for its convergence on the site of bvFTD neurodegeneration in the cortex and provides a strong impetus for further study of bvFTD risk gene expression patterns in these specific cell populations.

Limitations of this study include lack of a suitable case-control data set of sufficient size to confirm whether our network and risk loci are associated with disease risk and lack of access to raw IFGC genotype data to perform subgroup analyses. Phase III data of the IFGC (on over 2,000 new cases) will be released in the near future and that cohort could represent a suitable cohort for validation of these findings and further characterization efforts. Our data-sharing agreement with the IFGC provided SNP-level GWAS results only and we did not have raw genetic data, which precludes potentially more powerful study designs, such as subsetting the cohort into training and testing data sets. Beyond this, our network generation methodology relies on already existing protein interaction data and thus could bias our ontological findings toward more thoroughly characterized biological pathways and processes. An additional limitation of our gene expression analyses stems from Constantinidis diagnoses, which are not easily converted into more modern pathologic diagnoses.

In our analyses, we identified and bioinformatically characterized a gene network linked at the protein level that shows enriched expression across the human lifespan in brain areas most affected by bvFTD and shows dysregulation in pathologically diagnosed FTLD cases. These findings suggest that multiple distinct biological pathways are altered in sporadic bvFTD, including signal transduction, cell cycle regulation, immune/inflammation, and neurodevelopmental processes. Together, these pathway enrichments raise an important and unanswered question regarding the temporal relationship between genetic variation and biological risk for disease. It is plausible that a subset of risk variants promotes disease risk specifically during development, while others promote risk through nuanced changes in biological function in aging. It is also possible that genetic variation may be helpful during development but detrimental during aging if activated inappropriately (e.g., synaptic pruning by the innate immune system). Future longitudinal studies of at-risk individuals will be required to fully address this question. Although no single genetic variant used to generate our network directly causes bvFTD, the cumulative burden of multiple sub-GWAS-associated variants within these pathways may alter the cellular landscape on which development and aging occurs, thereby altering the amount of metabolic stress, inflammation, and apoptotic tendencies of particular cell types in defined brain networks. Our results support the utility of polygenic scores that incorporate the effects of multiple genetic loci for clinical risk assessment and prevention study enrichment, and genetic overlap with AD suggests shared underlying biology across multiple neurodegenerative disorders. Understanding the temporal tropism of these risk factors on disease vulnerability will also be critical for translating these results into therapeutic targets for clinical populations, especially for processes such as metabolic stress and immune dysfunction that are more amenable to intervention. If confirmed, our findings represent a valuable step toward a genetic understanding of selective vulnerability in neurodegenerative disease.

Authors contributions

L.W. Bonham: design and conceptualization of the study, analysis and interpretation of the data, drafting and revising the manuscript for intellectual content. N.Z.R. Steele: analysis and interpretation of the data, drafting and revising the manuscript for intellectual content. C.M. Karch: analysis and interpretation of the data, revising the manuscript for intellectual content. C. Manzoni: interpretation of the data, revising the manuscript for intellectual content. E.G. Geier: interpretation of the data, revising the manuscript for intellectual content. N. Wen, interpretation of the data, revising the manuscript for intellectual content. A. Ofori-Kuragu: analysis and interpretation of the data, revising the manuscript for intellectual content. P. Momeni: design and conceptualization of the study, revising the manuscript for intellectual content. J. Hardy: design and conceptualization of the study, revising the manuscript for intellectual content. Z.A. Miller: design and conceptualization of the study, revising the manuscript for intellectual content. C.P. Hess: design and conceptualization of the study, revising the manuscript for intellectual content. P. Lewis: design and conceptualization of the study, revising the manuscript for intellectual content. B.L. Miller: design and conceptualization of the study, revising the manuscript for intellectual content. W.W. Seeley: design and conceptualization of the study, interpretation of the data, revising the manuscript for intellectual content. S.E. Baranzini: design and conceptualization of the study, interpretation of the data, revising the manuscript for intellectual content. R.S. Desikan: design and conceptualization of the study, interpretation of the data, revising the manuscript for intellectual content. R. Ferrari: design and conceptualization of the study, analysis and interpretation of the data, revising the manuscript for intellectual content. J.S. Yokoyama: design and conceptualization of the study, analysis and interpretation of the data, drafting and revising the manuscript for intellectual content.

Study funding

Primary support for data analyses was provided by the Radiologic Society of North America (RSNA) RMS1741 (L.W.B.), Larry L. Hillblom Foundation 2012-A-015-FEL and 2016-A-005-SUP (J.S.Y.), AFTD Susan Marcus Memorial Fund Clinical Research Grant (J.S.Y.), NIA K01 AG049152 (J.S.Y.), Bluefield Project to Cure FTD (J.S.Y.), Tau Consortium (J.S.Y.), NIA K01 AG046374 (C.M.K.), U24DA041123 (R.S.D.), National Alzheimer's Coordinating Center (NACC) Junior Investigator Award (R.S.D.), RSNA Resident/Fellow Grant (R.S.D.), Foundation of ASNR Alzheimer's Imaging Grant (R.S.D.), and Alzheimer's Society Grant 284 (R.F.). Additional support was provided by an MRC Programme grant (MR/N026004/1; J.H. and P.A.L.), NIA P01 AG1972403 (B.L.M.), NIA P50 AG023501 (B.L.M.), and a MRC New Investigator Research Grant (MR/L010933/1; P.A.L.). P.A. Lewis is a Parkinson's UK research fellow (grant F1002). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Disclosure

L.W. Bonham has received research support from the Radiologic Society of North America. N.Z.R. Steele, C.M. Karch, N. Wen, A. Ofori-Kuragu, P. Momeni report no disclosures. C. Manzoni has received a bursary from NIH; serves on the editorial board of Frontiers; and has receives research support/salary from the Medical Research Council. E.G. Geier holds a patent (pending) for Targeted platinum anticancer agents. J. Hardy has served on scientific advisory boards for Eisai and Lilly; has been a consultant for Eisai and Ceracuity; and has received research support from the Medical Research Council. Z.A. Miller holds patents for Dynamic and adjustable portable X-ray filter grid and Dynamic and adjustable filter grids; and has received research support from NIH. C.P. Hess serves on the editorial boards of the American Journal of Neuroradiology and PLoS One; has been a consultant for Imaging Endpoints; has received research support from General Electric and NIH; and has provided expert witness for various legal trials. P. Lewis has received travel funding from the Biochemical Society; serves on the editorial boards of Biology Direct, PLoS One, and Frontiers in Neurodegenerative Disease; and has received research support from the Medical Research Council, the Biotechnology and Biological Sciences Research Council, Diamond Light Source, and Parkinson's UK. B.L. Miller has served on the scientific advisory boards of the Tau Consortium, the John Douglas French Foundation, the Larry L. Hillblom Foundation, the Consortium for Frontotemporal Dementia Research, the Global Brain Health Institute (GBHI), University of Washington ADRC, Stanford University ADRC, the Arizona Alzheimer's Disease Center (ADC), the Massachusetts Alzheimer Disease Research Center, and the International Society of FTD; has served on the editorial boards of Cambridge University Press, Guilford Publications Inc, Oxford University Press, Neurocase, Elsevier, Up-To-Date, and Frontiers; receives publishing royalties from Cambridge University Press, Elsevier, Guilford Publications, and Oxford University Press; and has received research support from Quest Diagnostics Incorporated, and NIA. W.W. Seeley serves on the editorial boards of Annals of Neurology, Acta Neuropathologica, and Neuroimage: Clinical; has been a consultant for Bristol-Myers Squibb, Merck, and Biogen Idec; has received research support from NIH; and has provided expert witness for a defendant. S.E. Baranzini has served on scientific advisory boards for Novartis, EMD Serono, and Sanofi-Aventis; has received gifts, travel funding, and speaker honoraria from Novartis; serves on the editorial board of the Multiple Sclerosis Journal; holds a patent (pending) for a gene expression signature that could identify patients at high risk of developing multiple sclerosis; has been a consultant for Novartis, EMD Serono, and TEVA Neuroscience; and has received research support from NIH, the Department of Defense, and NMSS. R.S. Desikan has received research support from NIH, the National Alzheimer's Coordinating Center (NACC) Junior Investigator Award, RSNA, and the Foundation of ASNR. R. Ferrari has received travel funding/registration fee reimbursement from Sindem4Juniors and the Alzheimer's Association International Conference; and has received research support from the Alzheimer's Society. J.S. Yokoyama has served on the editorial boards of the Journal of Alzheimer's Disease and Frontiers in Neurology; and has received research support from Quest Diagnostics, NIH, the Association for Frontotemporal Degeneration, the Larry L. Hillblom Foundation, the Rainwater Charitable Foundation, the John Douglas French Alzheimer's Foundation, and the Bluefield Project to Cure FTD. Full disclosure form information provided by the authors is available with the full text of this article at Neurology.org/NG.

Acknowledgment

The authors thank the International FTD-Genomics Consortium (IFGC) for providing phase I summary statistics data for these analyses. IFGC acknowledgements and full collaborator list are provided in appendix e-1 (links.lww.com/NXG/A71).

Footnotes

  • Funding information and disclosures are provided at the end of the article. Full disclosure form information provided by the authors is available with the full text of this article at Neurology.org/NG.

  • The Article Processing Charge was funded by the authors.

  • Received January 29, 2018.
  • Accepted in final form July 31, 2018.

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND), which permits downloading and sharing the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.

References

  1. 1.↵
    1. Seeley WW,
    2. Crawford RK,
    3. Zhou J,
    4. Miller BL,
    5. Greicius MD
    . Neurodegenerative diseases target large-scale human brain networks. Neuron 2009;62:42–52.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Warren JD,
    2. Rohrer JD,
    3. Schott JM,
    4. Fox NC,
    5. Hardy J,
    6. Rossor MN
    . Molecular nexopathies: a new paradigm of neurodegenerative disease. Trends Neurosci 2013:561–569.
  3. 3.↵
    1. Zhou J,
    2. Gennatas ED,
    3. Kramer JH,
    4. Miller BL,
    5. Seeley WW
    . Predicting regional neurodegeneration from the healthy brain functional connectome. Neuron 2012;73:1216–1227.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Boyle EA,
    2. Li YI,
    3. Pritchard JK
    . An expanded view of complex traits: from polygenic to omnigenic. Cell 2017;169:1177–1186.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Baranzini S,
    2. Khankhanian P,
    3. Patsopoulos N, et al
    . Network-based multiple sclerosis pathway analysis with GWAS data from 15,000 cases and 30,000 controls. Am J Hum Genet 2013;92:854–865.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Wang C,
    2. Ji F,
    3. Hong Z, et al.
    Disrupted salience network functional connectivity and white-matter microstructure in persons at risk for psychosis: findings from the LYRIKS study. Psychol Med 2016;46:2771–2783.
    OpenUrl
  7. 7.↵
    1. Rascovsky K,
    2. Hodges JR,
    3. Kipps CM, et al
    . Diagnostic criteria for the behavioral variant of frontotemporal dementia (bvFTD): current limitations and future directions. Alzheimer Dis Assoc Disord 2007;21:S14–S18.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Ferrari R,
    2. Hernandez DG,
    3. Nalls MA, et al
    . Frontotemporal dementia and its subtypes: a genome-wide association study. Lancet Neurol 2014;13:686–699.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Liu JZ,
    2. McRae AF,
    3. Nyholt DR, et al
    . A versatile gene-based test for genome-wide association studies. Am J Hum Genet 2010;87:139–145.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Frazer KA,
    2. Ballinger DG,
    3. Cox DR
    ; The International HapMap Corsortium. A second generation human haplotype map of over 3.1 million SNPs. Nature 2007;449:851–861.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Ferrari R,
    2. Lovering RC,
    3. Hardy J,
    4. Lewis PA,
    5. Manzoni C
    . Weighted protein interaction network analysis of frontotemporal dementia. J Proteome Res 2017;16:999–1013.
    OpenUrlCrossRef
  12. 12.↵
    1. Zhang Y,
    2. Sloan SA,
    3. Clarke LE, et al
    . Purification and characterization of progenitor and mature human astrocytes reveals transcriptional and functional differences with mouse. Neuron 2016;89:37–53.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Sheng J,
    2. Su L,
    3. Xu Z,
    4. Chen G
    . Progranulin polymorphism rs5848 is associated with increased risk of Alzheimer's disease. Gene 2014;542:141–145.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Schmidt EF,
    2. Warner-Schmidt JL,
    3. Otopalik BG,
    4. Pickett SB,
    5. Greengard P,
    6. Heintz N
    . Identification of the cortical neurons that mediate antidepressant responses. Cell 2012;149:1152–1163.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Reimand J,
    2. Arak T,
    3. Adler P, et al
    . g:Profiler: a web server for functional interpretation of gene lists (2016 update). Nucleic Acids Res 2016;44:W83–W89.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Allen M,
    2. Carrasquillo MM,
    3. Funk C, et al
    . Human whole genome genotype and transcriptome data for Alzheimer's and other neurodegenerative diseases. Sci data 2016;3:160089.
    OpenUrl
  17. 17.↵
    1. Lesnick TG,
    2. Papapetropoulos S,
    3. Mash DC, et al
    . A genomic pathway approach to a complex disease: axon guidance and Parkinson disease. PLoS Genet 2007;3:0984–0995.
    OpenUrl
  18. 18.↵
    1. Kang HJ,
    2. Kawasawa YI,
    3. Cheng F, et al
    . Spatio-temporal transcriptome of the human brain. Nature 2011;478:483–489.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Bennett ML,
    2. Bennett FC,
    3. Liddelow SA, et al
    . New tools for studying microglia in the mouse and human CNS. Proc Natl Acad Sci 2016;113:E1738–E1746.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Dijkstra AA,
    2. Lin LC,
    3. Nana AL,
    4. Gaus SE,
    5. Seeley WW
    . Von economo neurons and fork cells: a neurochemical signature linked to monoaminergic function. Cereb Cortex 2016;45:846–860.
    OpenUrl
  21. 21.
    1. Chen-Plotkin AS,
    2. Geser F,
    3. Plotkin JB, et al
    . Variations in the progranulin gene affect global gene expression in frontotemporal lobar degeneration. Hum Mol Genet 2008;17:1349–1362.
    OpenUrlCrossRefPubMed
  22. 22.
    1. Bronner IF,
    2. Bochdanovits Z,
    3. Rizzu P, et al
    . Comprehensive mRNA expression profiling distinguishes tauopathies and identifies shared molecular pathways. PLOS ONE 2009;4:e6826.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Ferrari R,
    2. Wang Y,
    3. Vandrovcova J, et al
    . Genetic architecture of sporadic frontotemporal dementia and overlap with Alzheimer’s and Parkinson’s diseases. J Neurol Neurosurg Psychiatry 2017;88:152–164.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Iqbal K,
    2. Liu F,
    3. Gong CX
    . Tau and neurodegenerative disease: the story so far. Nat Rev Neurol 2015;12:15–27.
    OpenUrl
  25. 25.↵
    1. Sawyer RP,
    2. Rodriguez-Porcel F,
    3. Hagen M,
    4. Shatz R,
    5. Espay AJ
    . Diagnosing the frontal variant of Alzheimer’s disease: a clinician’s yellow brick road. J Clin Mov Disord 2017;4:2.
    OpenUrl
  26. 26.↵
    1. Brettschneider J,
    2. Del Tredici K,
    3. Lee VM,
    4. Trojanowski JQ
    . Spreading of pathology in neurodegenerative diseases: a focus on human studies. Nat Rev Neurosci 2015;16:109–120.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Cohen TJ,
    2. Lee VM,
    3. Trojanowski JQ
    . TDP-43 functions and pathogenic mechanisms implicated in TDP-43 proteinopathies. Trends Mol Med 2011;17:659–667.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Yang Y,
    2. Higashimori H,
    3. Morel L
    . Developmental maturation of astrocytes and pathogenesis of neurodevelopmental disorders. J Neurodev Disord 2013;5:22.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Ferrari R,
    2. Wang Y,
    3. Vandrovcova J, et al
    . Genetic architecture of sporadic frontotemporal dementia and overlap with Alzheimer’s and Parkinson’s diseases. J Neurol Neurosurg Psychiatry 2017;88:152–164.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Becker EB,
    2. Bonni A
    . Cell cycle regulation of neuronal apoptosis in development and disease. ProgNeurobiol 2004;72:1–25.
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. Copani A,
    2. Caraci F,
    3. Hoozemans JJM,
    4. Calafiore M,
    5. Angela Sortino M,
    6. Nicoletti F
    . The nature of the cell cycle in neurons: focus on a “non-canonical” pathway of DNA replication causally related to death. Biochim Biophys Acta 2007;1772:409–412.
    OpenUrlPubMed
  32. 32.↵
    1. Minami SS,
    2. Min SW,
    3. Krabbe G, et al
    . Progranulin protects against amyloid β deposition and toxicity in Alzheimer’s disease mouse models. Nat Med 2014;20:1157–1164.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Santillo AF,
    2. Englund E
    . Greater loss of von Economo neurons than loss of layer II and III neurons in behavioral variant frontotemporal dementia. Am J Neurodegener Dis 2014;3:64–71.
    OpenUrl
  34. 34.
    1. Chow F,
    2. Gong Y,
    3. Lippa CF
    . The potential role of insulin on the shank-postsynaptic platform in neurodegenerative diseases involving cognition. Am J Alzheimer’s Dis Other Demen 2014;29:303–310.
    OpenUrlCrossRefPubMed
  35. 35.
    1. Kang J,
    2. Park H,
    3. Kim E
    . IRSp53/BAIAP2 in dendritic spine development, NMDA receptor regulation, and psychiatric disorders. Neuropharmacology 2016;100:27–39.
    OpenUrl
  36. 36.
    1. Narayanan RK,
    2. Mangelsdorf M,
    3. Panwar A,
    4. Butler TJ,
    5. Noakes PG,
    6. Wallace RH
    . Identification of RNA bound to the TDP-43 ribonucleoprotein complex in the adult mouse brain. Amyotroph Lateral Scler Front Degener 2013;14:252–260.
    OpenUrl
  37. 37.
    1. Lue LF,
    2. Schmitz CT,
    3. Snyder NL, et al
    . Converging mediators from immune and trophic pathways to identify Parkinson disease dementia. Neurol Neuroimmunol Neuroinflamm 2016;3:e193.
    OpenUrl
  38. 38.
    1. Hooper C,
    2. Killick R,
    3. Lovestone S
    . The GSK3 hypothesis of Alzheimer’s disease. J Neurochem 2008;104:1433–1439.
    OpenUrlCrossRefPubMed
  39. 39.
    1. Trojanowski JQ,
    2. Duff K,
    3. Fillit H, et al
    . New directions for frontotemporal dementia drug discovery. Alzheimers Dement 2008;4:89–93.
    OpenUrlPubMed
  40. 40.
    1. Sekine Y,
    2. Takeda K,
    3. Ichijo H
    . The ASK1-MAP kinase signaling in ER stress and neurodegenerative diseases. Curr Mol Med 2006;6:87–97.
    OpenUrlCrossRefPubMed
  41. 41.
    1. Sousa SB,
    2. Hennekam RC,
    3. Abdul-Rahman O, et al
    . Phenotype and genotype in Nicolaides-Baraitser syndrome. Am J Med Genet C Semin Med Genet 2014;166:302–314.
    OpenUrl
  42. 42.
    1. Zhang YJ,
    2. Gendron TF,
    3. Xu YF,
    4. Ko LW,
    5. Yen SH,
    6. Petrucelli L
    . Phosphorylation regulates proteasomal-mediated degradation and solubility of TAR DNA binding protein-43 C-terminal fragments. Mol Neurodegener 2010;5:1–13.
    OpenUrlCrossRefPubMed
  43. 43.
    1. Moradifard S,
    2. Hoseinbeyki M,
    3. Ganji SM,
    4. Minuchehr Z
    . Analysis of microRNA and gene expression profiles in Alzheimer’s disease: a meta-analysis approach. Sci Rep 2018;8:4767.
    OpenUrl
  44. 44.
    1. Dobson-Stone C,
    2. Luty AA,
    3. Thompson EM, et al
    . Frontotemporal dementia-amyotrophic lateral sclerosis syndrome locus on chromosome 16p12.1-q12.2: genetic, clinical and neuropathological analysis. Acta Neuropathol 2013;125:523–533.
    OpenUrl
View Abstract

Letters: Rapid online correspondence

No comments have been published for this article.
Comment

NOTE: All contributors' disclosures must be entered and current in our database before comments can be posted. Enter and update disclosures at http://submit.ng.neurology.org. Exception: replies to comments concerning an article you originally authored do not require updated disclosures.

  • Stay timely. Submit only on articles published within the last 8 weeks.
  • Do not be redundant. Read any comments already posted on the article prior to submission.
  • 200 words maximum.
  • 5 references maximum. Reference 1 must be the article on which you are commenting.
  • 5 authors maximum. Exception: replies can include all original authors of the article.
  • Submitted comments are subject to editing and editor review prior to posting.

More guidelines and information on Letters

Compose Comment

More information about text formats

Plain text

  • No HTML tags allowed.
  • Web page addresses and e-mail addresses turn into links automatically.
  • Lines and paragraphs break automatically.
Author Information
NOTE: The first author must also be the corresponding author of the comment.
First or given name, e.g. 'Peter'.
Your last, or family, name, e.g. 'MacMoody'.
Your email address, e.g. higgs-boson@gmail.com
Your role and/or occupation, e.g. 'Orthopedic Surgeon'.
Your organization or institution (if applicable), e.g. 'Royal Free Hospital'.
Publishing Agreement
NOTE: All authors, besides the first/corresponding author, must complete a separate Letters Submission Form and provide via email to the editorial office before comments can be posted.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.

Vertical Tabs

You May Also be Interested in

Back to top

Related Articles

  • No related articles found.

Topics Discussed

Alert Me

Advertisement