- Research article
- Open Access
Multiscale network modeling of oligodendrocytes reveals molecular components of myelin dysregulation in Alzheimer’s disease
Molecular Neurodegeneration volume 12, Article number: 82 (2017)
Oligodendrocytes (OLs) and myelin are critical for normal brain function and have been implicated in neurodegeneration. Several lines of evidence including neuroimaging and neuropathological data suggest that Alzheimer’s disease (AD) may be associated with dysmyelination and a breakdown of OL-axon communication.
In order to understand this phenomenon on a molecular level, we systematically interrogated OL-enriched gene networks constructed from large-scale genomic, transcriptomic and proteomic data obtained from human AD postmortem brain samples. We then validated these networks using gene expression datasets generated from mice with ablation of major gene expression nodes identified in our AD-dysregulated networks.
The robust OL gene coexpression networks that we identified were highly enriched for genes associated with AD risk variants, such as BIN1 and demonstrated strong dysregulation in AD. We further corroborated the structure of the corresponding gene causal networks using datasets generated from the brain of mice with ablation of key network drivers, such as UGT8, CNP and PLP1, which were identified from human AD brain data. Further, we found that mice with genetic ablations of Cnp mimicked aspects of myelin and mitochondrial gene expression dysregulation seen in brain samples from patients with AD, including decreased protein expression of BIN1 and GOT2.
This study provides a molecular blueprint of the dysregulation of gene expression networks of OL in AD and identifies key OL- and myelination-related genes and networks that are highly associated with AD.
Alzheimer’s disease (AD) causes a progressive dementia that affects approximately 1/6th of people in the US age 75 and above . Although the number one risk factor for AD is advanced age, the reason for this remains unknown [2, 3]. A wealth of evidence has emerged over the past decade to support the role of non-neuronal cells, especially astrocytes and microglia, in amyloid beta (Aβ) processing and AD pathogenesis [4,5,6]. While less well studied, several lines of evidence have suggested that dysregulation of oligodendrocytes (OLs) and associated dysmyelination might be important in AD pathology. For example, human neuroimaging studies have shown that white matter changes occur early in AD and are predictive of disease status [7,8,9]. In particular, MRI studies have detected white matter volume atrophy in multiple brain regions prior to changes in gray matter in AD progression [10,11,12]. Further, post-mortem human pathological studies have demonstrated that the pattern of neurofibrillary tangle deposition in AD is highly correlated with the developmental pattern of myelination, with late-myelinated axonal tracts substantially more vulnerable to degeneration in AD [13, 14]. Recent reports have highlighted the importance of OLs and myelin metabolic function for axonal health and transport capacity [15,16,17,18,19,20], thereby suggesting the possibility that OL-driven axonal damage could precede a secondary demyelination in the pathogenesis of AD. In particular, there are several mouse models of ablation of the OL-associated myelin genes Ugt8 , Cnp [22, 23], and Plp1 [24, 25], in which axonal degeneration occurs in the presence of minimal ultrastructural myelin alterations and therefore are well suited to study altered OL gene expression, presumably leading to myeling dysfunction preceding the onset of neurodegeneration. To investigate the hypothesis that OL dysregulation in AD may be part of the underlying mechanism leading to neurodegeneration, we sought to employ a detailed molecular and systems-level analysis to provide a molecular substrate for the potential role of OLs in mediating the initial axonal damage.
In this study, we systematically examined and validated OL-enriched gene networks to uncover key genes and molecular signaling circuits of OLs in AD. We built upon AD-associated and OL-enriched networks constructed in a previous study of genetic, gene expression, and pathophysiologic data in late-onset AD . We constructed a union of the three OL-enriched modules from a multi-tissue AD co-expression network and found that it was strongly enriched for AD risk factor genes. Our OL-enriched consensus module includes genes encoding proteins associated with Aβ-production PSEN1 and BACE1, as well as the AD risk factor genes BIN1, PSEN1, PICALM, and UNC5C [27,28,29,30]. We next built co-expression networks from a large-scale proteomics data set, identifying a strong loss of coordination among proteins in the most OL-enriched network, an interaction of this dysregulation and dementia status, as well as a down-regulation of key OL network genes, including BIN1. We then used the OL modules to construct regulatory networks and found that the topological structures of our OL-enriched networks were validated through in vitro and in vivo perturbations of the predicted key regulatory genes in the networks. Further, transcriptomic analysis of brain tissue isolated from mice with a genetic ablation of three top key driver genes (Ugt8, Cnp, and Plp1) recapitulated key aspects of the dysregulation in gene pathways related to myelination that are seen in human AD brains. We chose to profile the mice with these genetic ablations at an early stage of development (postnatal day 20), in order to detect alterations in gene pathways occurring during the process of myelination and prior to the onset of widespread axon degeneration in Cnp-KO  or Plp1-KO  mice. This approach allowed us to define alterations in OL gene expression occurring during the “pre-symptomatic” or prodromal phase of the disease, as opposed to reactive changes that might be consequent to axonal degeneration. We found that differentially expressed gene (DEG) signature in the Cnp knockout (KO) mouse mimicked gene expression changes detected in AD brains at the early stages of the pathology, both at the gene pathway and individual protein level, thereby suggesting that dysregulation of OLs in general and CNP in particular may play a key role in driving AD-associated gene expression changes.
Oligodendrocyte network construction from AD data and key driver analysis
We re-analyzed data from the Harvard Brain Tissue Resource Center (HBTRC) consisting of 376 late-onset Alzheimer’s disease brain samples as well as 173 non-demented brain samples, harvested from the cerebellum (CBM), dorsolateral prefrontal cortex (Brodmann area 9; henceforth, PFC), and visual cortex (Brodmann area 7; VC), which has been previously described  (GEO: GSE44772). Based upon the HBTRC AD data, we previously built up a multi-tissue weighted gene co-expression network using WGCNA [26, 31]. In order to find the AD OL-specific modules, we tested the enrichment of each module for the gene signatures specifically expressed in OLs. To do this, we employed datasets on cell enrichment and identified the genes most highly expressed in myelinating OLs compared to the non-OL cell types . Three modules, one with probes primarily derived from each of the three investigated brain regions were found to be the most enriched for the OL gene signature. We therefore refer to these as the “region-specific oligodendrocyte co-expression network” modules. We further inferred an interaction network for each OL-enriched and region-specific gene module by integrating gene expression and DNA genotypic data, as previously described . Specifically, relationships between genes in each co-expression module were inferred based on conditional independence tests of gene expression in a Bayesian framework, and the presence of more expression single nucleotide polymorphisms (eSNPs) associated with a gene were used as priors to break Markov equivalence . We refer to each of these networks as the “region-specific OL-enriched Bayesian interaction networks”. For robustness, we combined the genes in the three region-specific co-expression networks into a core OL gene set (COLGS), merged the three brain region-specific Bayesian interaction networks by a set union of directed links into a core OL-enriched Bayesian interaction network (COLBN). The gene symbols from this previous study were updated to current Hugo Gene Nomenclature Committee (HGNC) standards using the R package HGCNhelper (v. 0.3.1), and the two gene symbols in the OL-enriched networks that were still mapped to multiple symbols were mapped to the more common symbol (MARCH1 and LPAR1) for legibility. In order to identify key regulatory genes in this interaction network, we then performed key driver analysis on COLBN [34, 35], which uses network connectivity to infer regulatory importance scores for genes.
Enrichment of brain co-expression modules in AD GWAS risk factor genes
In order to find the enrichment of COLGS and other co-expression modules in Alzheimer’s risk factor genes, we converted the International Genomics of Alzheimer’s Project (IGAP) GWAS SNP-level data set  to gene-level p-value calls using VEGAS2 , and used the genes with significant association at a nominal p < 0.05 for further analysis. We then measured the enrichment of this 543 AD GWAS risk factor gene set in the 62 overall qualifying multiscale AD modules with at least 50 genes using Fisher’s Exact Test (FET).
Proteomics data analysis from PFC human postmortem brain samples
Grey matter brain samples in 50 mg aliquots were harvested from the prefrontal cortex (PFC; Brodmann Area 10) from the autopsied brains of persons with a wide range of cognitive status at the time of death, ranging from no cognitive impairment to dementia, as well as a wide range of Braak scores, from 0 to 6. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to measure the abundance of peptides in each brain sample, from which a protein-level quantitation was estimated using MaxQuant (v188.8.131.52). We used WGCNA to define modules of proteins in the proteomics data, with a soft-thresholding power coefficient of 3. We next annotated these modules based on their relative FET enrichment in the human homologues of genes specifically expressed in each of the five major mouse brain cell types (see “Estimating brain cell type enrichment”).
Within this OL-enriched module, we performed modular differential connectivity analysis in order to assess the overall difference in between samples classified as non-AD (Braak 0–2) and AD (Braak 5–6), calculated using the mean difference in z-scores option in DGCA  (v. 1.0.1), with 10,000 permutation samples to assess significance. We used a Student’s t-test to measure differences in average expression between conditions, using the q value R package  to estimate false discovery rates in the context of the multiple hypothesis tests.
Generating in vivo mouse genetic perturbation signatures
Use of animals in this research was strictly compliant with the guidelines set forth by the US Public Health Service in their policy on Humane Care and Use of Laboratory Animals, and in the Guide for the Care and Use of Laboratory Animals. All animal procedures received prior approval from the Institutional Animal Care and Use Committee at Icahn School of Medicine at Mount Sinai.
For each of the control and knockout model mice of the three key drivers (Cnp, Ugt8, and Plp1), mice of either sex were sacrificed at postnatal day 20 and tissue was flash-frozen in liquid nitrogen vapors. The frontal cortex and cerebellum were dissected on ice and immediately processed for RNA isolation. Integrity of the RNA was confirmed by measuring RNA Integrity Number (RIN) and only samples with RIN > 8.5 were used for RNA sequencing.
RNA isolation and reverse transcription
RNA was isolated using Trizol reagent (Invitrogen, CA) and cleaned using RNeasy Mini kit (Qiagen, CA). Ribosomal RNA was removed from the samples using Ribo-Zero rRNA Removal Kit (Illumina, CA). Approximately 500 nanograms of total RNA was used in cDNA library construction with the TruSeq RNA Sample Prep Kit (Illumina, CA), followed by RNA sequencing using an Illumina HiSeq2000.
Read mapping and quantification of RNAseq data
RNA-sequencing reads were mapped to the mouse genome (mm10, UCSC assembly) using Bowtie (version 184.108.40.206), TopHat (version 2.0.11), and SamTools (version 0.1.19.0) using a read length of 100. For RNAseq knockout experiments, reads were converted to counts at the gene level using HTSeq  on the BAM files from TopHat2 using the UCSC known genes data set. For each key driver knockout and brain region, genes that mapped less than 100 counts in 80% or more of the samples were filtered out from downstream analysis, because these genes are likely to have especially high variance in expression calls.
Compartmental approach to compare of mouse key driver perturbation signatures with human postmortem AD gene expression signatures
We found the enrichment of both mouse key driver knockout DEG signatures, aggregated across significant DEGs (p < 0.05, FDR < 0.3) from both the frontal cortex and cerebellum, as well as two human AD gene expression signatures in the gene ontology pathways. The prefrontal cortex (PFC) human AD DEG signature was estimated from the HBTRC cohort. We used the list of genes identified as previously identified as having significantly different RNA levels in brain samples from persons with high levels of Alzheimer’s neuropathology (Braak = 5–6) compared brain samples from persons with low levels of Alzheimer’s neuropathology (Braak = 0–2) via a Student’s t-test . The second human AD DEG signature is from a previous study that identified genes with a significant trend in RNA expression changes across the severity spectrum of AD (control, incipient, moderate, and severe) in the hippocampus (HIPP) . In this data set, we selected the genes as having an increasing or decreasing trend in expression across AD severity stages and an ANOVA p-value <0.05 as the HIPP AD DEG signature. We used the moduleGO function in DGCA  (version 1.0.1) to perform gene ontology (GO) enrichment analysis on these DEG sets, which leverages the GOstats (version 2.34)  and org.Hs.eg.db GO annotation (version 3.1.2) R packages. We filtered for those GO terms with less than 800 and greater than 100 gene symbols. We adjusted the enrichment p-values for all GO terms in each DEG set using the Benjamani-Hochberg method. In order to determine whether there was a similar pattern of dysregulation in the mouse key driver knockout models as in human AD, we next found the degree to which the compartment overlaps intersected, using the R package SuperExactTest (version 0.99.2) .
In vivo validation of gene expression
The Cnp-cre knock-in mouse line has been described previously . For the sake of clarity, in this manuscript, Cnp +/+ mice are called Cnp-WT and Cnp cre/cre mice are called Cnp-KO.
Mice of either sex were sacrificed at postnatal day 60 and brains were cut into 1 mm coronal sections using a refrigerated brain matrix. Corpus callosum (CC) was dissected from coronal sections on ice using a light microscope. The tissue was then immediately processed for protein extraction. Protein lysates (50 μg) were separated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE) and transferred onto a PolyVinylidene DiFluoride (PVDF) (Millipore, Billerica, MA, USA) membrane using a buffer containing 25 mM Tris base, pH 8.3, 192 mM glycine, 20% (vol/vol) methanol for 1 h at 100 V at 4 °C. Membranes were blocked for 1 h in 10% Milk/0.1% Tween/TBS, then incubated overnight at 4 °C with the primary antibody diluted 1:1000 in 5% BSA/0.02% sodium azide/0.1% Tween/TBS. Antibodies used are: mouse anti-CNPase (Sternberger Inc., SMI-91, 1:5,00), mouse anti-alpha-TUBULIN (Calbiochem, CP06, 1:5000), rabbit anti-BIN1 (Abcam, ab185950, 1:2500), rabbit anti-GOT2 (Abcam, ab171739, 1:2500). After rinsing with 0.1% Tween/TBS, membranes were incubated 2 h at room temperature with the secondary light-chain specific antibody (Jackson Immunoresearch, 1:10,000) in 10% Milk/0.1% Tween/TBS. After rinsing, membranes were incubated with ECL (Amersham) for 3 min and then revealed. Quantification was carried out on three biological and technical replicates per genotype, using ImageJ. Protein expressions were normalized to alpha-TUBULIN expression, then compared to their respective expression levels in Cnp-WT samples.
Mice were perfused with 4% paraformaldehyde and post-fixed overnight in the same solution at 4 °C. Tissue samples were then transferred to 70% ethanol, sequentially dehydrated and embedded in paraffin. Four-micrometer sections were cut, deparaffinized and rehydrated. Antigen retrieval was performed by incubating slides in sub-boiling (94 °C) citrate buffer (pH 6.0) for 15 min. Slides were incubated in blocking buffer (20% Normal Goat Serum / 1% BSA / PBS 1X) for 1 h at room temperature and then incubated overnight at 4 °C with the primary antibodies in 1% BSA / PBS 1X. Antibodies used are: mouse anti-OLIG2 (Millipore, MABN50, 1:200), mouse anti-NeuN (Millipore, MAB377, 1:200), mouse anti-GFAP (Sternberger Inc., SMI-22, 1:200), rabbit anti-BIN1 (Abcam, ab185950, 1:200), rabbit anti-GOT2 (Abcam, ab171739, 1:200). After rinsing with Tris-buffer / 2% milk, sections were processed with he appropriate Alexa Fluor conjugated secondary antibodies (1:1000 in 1% BSA / PBS 1X, Invitrogen), washed with Tris-buffer, and mounted using Fluoromount-G with DAPI. All images were acquired using a Zeiss Observer A1 fluorescent microscope or Zeiss LSM780 upright Confocal. Quantification was carried out on two-three sections per mouse and three-four mice for each genotype evaluated using ImageJ.
A robust myelination- and oligodendrocyte-enriched gene module is strongly associated with genetic data in late-onset AD
The primary goal of this study was to interrogate OL-enriched multiscale gene networks constructed from human AD postmortem brain samples (Fig. 1). We first identified co-expressed human gene modules, that were built from data obtained from three regions [e.g. prefrontal cortex (PFC), visual cortex (VC), and cerebellum (CBM)] from postmortem late-onset AD human brains samples from a large cohort of patients . We further annotated these gene modules with additional gene sets to investigate their associations with OLs and AD. We then constructed Bayesian gene regulatory networks for the OL-enriched co-expression modules, identified key regulatory genes in the regulatory networks, and systematically validated the topological structures of the regulatory networks based on a series of gene perturbation experiments in vitro and in vivo. Specifically, in order to identify OL co-expression networks in AD, we tested each of the 62 co-expression modules with at least 50 genes for the enrichment of genes expressed specifically in each of five major brain cell types, i.e., astrocytes, endothelial cells, microglia, neurons, and myelinating OLs  (Fig. 2). We identified three co-expression modules with the strongest enrichment of OL genes (Fig. 2), which we subsequently combined, leading to a set of 1631 unique gene symbols, which we henceforth refer to as the core oligodendrocyte gene set (COLGS). Notably, because the three OL-enriched co-expression modules were each primarily identified in one of the three brain regions from the multi-tissue experiment, by combining them we sought to obtain a more sensitive measure of genes whose expression is associated with OLs in the context of AD across brain regions. COLGS is highly enriched for genes encoding proteins identified in the myelin proteome  (Fold Enrichment (FE) = 1.92, Fisher’s Exact Test (FET) p = 2.4e-15), and is also enriched for genes specifically expressed in each of oligodendrocyte precursor cells (FE = 2.8, p = 4.2e-17), newly formed oligodendrocytes (FE = 6.2, p = 1.0e-84), and myelinating oligodendrocytes  (FE = 7.4, p = 2.5e-87), indicating that genes in COLGS capture a wide spectrum of OL functions.
To systematically evaluate the AD genetics of COLGS, we identified 543 genes with nominally significant (p < 0.05) gene-level associations with AD based on a meta-analysis of AD genome-wide association study (GWAS) data by the International Genomics of Alzheimer’s Project (IGAP) . We identified a significant enrichment of the AD risk genes in COLGS (FE = 1.71, p-value = 0.0004), due to the presence of BIN1, PICALM, NME8, SNX1, and other genes. Of all the individual co-expression modules, one OL-enriched module had the second strongest enrichment for AD GWAS risk genes, second only to an immune- and microglia-enriched module (Fig. 2, Table 1). We further highlighted the specific genes in COLGS that were also identified as AD risk genes (Fig. 3, Table 2). These results suggest that the COLGS contains a relatively large proportion of AD risk genes.
Next, we sought to examine the robustness of COLGS to a variety of sources of variance. We first used data from an independent study that also identified coexpression modules in human AD postmortem brain RNA expression samples . Despite originating from a different brain region (the hippocampus), we found that 83% of the genes in the module with the strongest OL-enrichment in this data set overlapped with the members of COLGS (FE = 18.3, p = 9.1e-79). In order to measure the robustness of the co-expression network with respect to an alternative gene expression modality, we utilized a large AD proteomic data set from the prefrontal cortex (PFC; Brodmann area 10), which we corrected for batch, age of death, and sex. We used WGCNA  to identify gene modules in this proteomic data set, and found that the module with the strongest OL-enrichment (FE = 1.7, p = 0.0053), has 50% overlap with COLGS (FE = 10.7, p = 6.6e-58), including the AD risk factor BIN1 (Additional file 1: Figure S1, Additional file 2). These results suggest a robust co-regulation of the COLGS genes at both the transcript and protein levels in AD.
Dysregulation of an oligodendrocyte-associated module in AD at the protein level
Our previous analysis showed that the three OL-enriched modules were each among the modules with the strongest loss of connectivity in AD samples on the RNA expression level . We sought to extend these results by interrogating gene expression dysregulation in the OL-enriched proteomics module, which includes 150 proteins and is the most OL-enriched module. We found that this module has a significant decrease in correlation in AD samples (mean difference in z-transformed correlation = −0.466, empirical p-value = 0.0489; Fig. 4a), thus validating the loss of coordination among OL network genes in AD at the protein level. To further explore gene expression changes in this module in AD, we measured differences in protein expression of the module members, identifying 17 proteins up-regulated in AD and 7 proteins down-regulated in AD at FDR < 0.3 and p-value <0.05 (Additional file 3). Notably, we found that a BIN1 protein isoform was down-regulated in AD (t = −2.4, p = 0.019, FDR = 0.19), as well as an MBP protein isoform (t = −2.3, p = 0.023, FDR = 0.19). Therefore, we found that there are both variable changes in expression levels for individual proteins as well as a loss of overall coordination of OL network protein expression in AD.
Construction and validation of an oligodendrocyte gene regulatory network in AD
We next sought to predict the gene-gene regulatory relationships among the COLGS genes in the RNA expression data, which profiles a wider set of genes than the proteomics data. Specifically, we constructed a Bayesian gene regulatory network for each of the three OL-enriched co-expression modules we identified by integrating RNA expression and genotype data from autopsied brain samples of persons with AD. As detailed in our previous studies [33, 34], Bayesian networks are a type of probabilistic causal network, providing a natural framework for integrating highly dissimilar genetic and gene expression data to predict regulatory relationships. We then combined the three OL Bayesian networks by a set union of directed links, leading to a more robust core oligodendrocyte-enriched Bayesian network (COLBN), and performed key driver analysis on COLBN to identify master regulatory OL genes in the context of AD (Fig. 1b; Table 3).
To interrogate the topology of the COLBN as well as how the dysregulation of COLBN key drivers could relate to AD, we identified an in vitro experiment that perturbed a key driver gene in the COLBN, MYRF, and examined how the predicted network structures correspond to its identified experimental targets. Specifically, we used data from a previous study that performed transcriptional profiling of cultured mouse OLs with a deletion of the myelination transcription factor Myrf (also known as C11orf9) . The set of genes differentially expressed in the cells with a Myrf deletion compared with the control was significantly enriched in the 5-layer downstream neighborhood of MYRF in the COLBN (FE = 2.5, p = 7.5e-33; Fig. 5a, Fig. 6a-b), thus validating the topology of the subnetwork regulated by MYRF. The validated downstream targets of MYRF include PLD1, which encodes a protein that regulates the shuttling of APP and is associated with APP in AD brains [46, 47], as well as APLP1, which encodes a protein that has been shown to accumulate within neuritic plaques .
Next, we validated the network structure of COLBN in vivo using knockout mouse models of three of the top 40 key drivers in COLBN, Ugt8, Cnp, and Plp1, each of which have been found to cause axon pathology without substantial myelin structural alterations [21, 23, 49]. We performed RNAseq (GEO: GSE80437) on tissue samples of postnatal day 20 mice from the frontal cortex (FC) and cerebellum (CBM) in these three mouse models to determine their KO DEG gene signatures at FDR < 0.3 (Additional files 4 and 5, Additional file 1: Figure S2). We profiled RNA from the FC and CBM because these two regions were also profiled in the human AD postmortem brain tissue study, from which the COLBN was generated. We selected this time point because we wanted to focus on gene changes that occurred prior to the development of axonal pathology and therefore more directly would allow us to determine gene changes characteristic of the prodromic stage of the disease.We identified no DEGs for ΔCnp in the CBM and ΔPlp1 in the FC, while we found that all of the other DEG signatures were significantly enriched in the downstream neighborhoods of their corresponding driver genes (Fig. 6c-d). UGT8 encodes an enzyme that transfers galactose to ceramide to generate galactosylceramide, which makes up approximately one-fourth of myelin lipid dry mass . UGT8’s 5-layer neighborhood is most enriched for the Ugt8 KO signature in the CBM (FE = 2.1, p = 2.4e-9). One validated downstream target of UGT8 is LIPA, which encodes a lysosomal cholesterol-metabolizing enzyme associated with genetic polymorphisms that affect plasma 24S–hydroxycholesterol/cholesterol levels in AD patients  (Fig. 5d). The key driver CNP encodes a protein that plays a role in OL process outgrowth and promoting axon survival [23, 52]. CNP’s 5-layer neighborhood is significantly enriched for the Cnp KO signature in the FC (FE = 1.7, p = 6.4e-6). One of the validated downstream targets of CNP is SEC14L5, which has been previously found to be down-regulated in CA1 in AD  (Fig. 5c). Finally, the key driver PLP1 encodes a protein that is the most abundant protein in CNS myelin sheaths , regulates OPC process outgrowth , and promotes axonal integrity . PLP1’s 5-layer neighborhood is significantly enriched for the Plp1 KO signature in the CBM (FE = 2.0, p = 1.2e-4). A validated downstream target of Plp1 is Fgf1 (Fig. 5b), which encodes a protein that has been found to stimulate neuronal growth and promote remyelination  and has been found to have increased levels in the CSF in AD .
Overall, as with the in vitro perturbation data, the in vivo perturbation experiments strongly support the network structure of COLBN and pinpoint specific downstream targets to explain how dysregulation of the key drivers may help mediate AD pathology.
Perturbations of key OL network drivers in mice mimic aspects of dysregulation in human AD brain samples
Since the key drivers of COLBN are predicted to orchestrate the expression of a network with a strong genetic association with AD, we hypothesized that the KO signatures of COLBN key drivers in vivo would mimic their gene expression dysregulation in human AD brains. In order to test this hypothesis, we performed gene ontology analysis on DEG signatures from postmortem brain tissue from both key driver KO mice and patients with AD. The DEG signatures in AD were derived based on samples from the PFC in the Harvard Brain Tissue Resource Center cohort  and samples from the hippocampus (HIPP) in the University of Kentucky Brain Bank cohort , because these regions are among those with the strongest AD pathology. Overall, the enrichment of DEGs from the human AD cases in these gene compartments show several similar dysregulation patterns as that seen in the key driver knockouts (Fig. 7a). For example, we identified a strong enrichment for genes in the gene ontology (GO) category “mitochondrial protein complex” in the down-regulated ΔCnp signature (FE = 10.6, p = 8.2e-22) and the down-regulated AD HIPP signature (FE = 6.7, p = 7.5e-11). We next found that the down-regulated ΔCnp and AD HIPP signatures intersected along with the mitochondrial gene set substantially more than expected due to chance (FE = 34, p = 2.1e-9; Fig. 7b), suggesting that the actual genes dysregulated in the mitochondria are similar between Cnp-KO mice and in the hippocampus of AD patients. For example, COX6A1, a mitochondrial-associated gene in which mutations are causative of peripheral neuropathy , was found to be down-regulated in both Cnp-KO mice and in the hippocampus of AD patients. We also identified a strong enrichment for genes in the GO category “ribosome” in the down-regulated ΔCnp signature (FE = 4.8, p = 3e-11), the down-regulated ΔPlp1 signature (FE = 8.1, p = 4e-18), and the down-regulated AD PFC signature (FE = 3.0, p = 0.002). We identified a significant overlap of all three of these down-regulated signatures and the ribosomal gene set (FE = 76, p = 0.013; Fig. 7c), which has limited power due to the small size of the overall ribosomal gene set. Since the Cnp and Plp1 knockout RNAseq profiling occurred prior to the typical age of onset of axon degeneration in these mice, the mitochondrial and ribosomal gene set dysregulation seen is suggestive of “presymptomatic” pathways that precede subsequent axonal and neuronal degeneration.
Finally, we also detected that there was a significant enrichment for genes in myelin sheath GO category in the up-regulated ΔPlp1 signature (FE = 4.3, p = 3e-8), the down-regulated ΔCnp signature (FE = 2.5, p = 0.002), and the down-regulated AD HIPP signature (FE = 5.2, p = 2e-9). However, the myelin sheath GO signature only contains 161 genes, which makes overlap analysis difficult. In order to use a more well-powered gene set, we utilized a larger set of 1778 genes that have been previously reported to be present in the myelin proteome . In this gene set, we found a strong enrichment for the intersection of genes in the myelin proteome with down-regulated genes in Cnp-KO mice and human AD hippocampus (FE = 6.0, p = 2.1e-8; Fig. 7d). This shared set of down-regulated myelin genes includes CDK5, which encodes an enzyme whose activity is associated with tau pathology  and is essential for myelination . Both Cnp-KO mice and patients with AD undergo axon damage in the absence dramatic ultrastructural changes in myelin structure , and our data shows that both conditions share a down-regulation of myelin-associated genes that may be associated with dysmyelination  and subsequent axon damage.
In order to validate the decreased expression of key myelin and mitochondrial genes in Cnp-KO mice, we measured the protein expression of BIN1 and GOT2 in Cnp-KO mice compared to WT mice. BIN1 is primarily expressed in mature OLs and white matter in AD patients . Interestingly, the longer, central nervous system (CNS)-specific isoform of Bin1 has decreased expression levels in AD, although overall Bin1 transcript levels show increased expression . This may reflect a proliferation of monocyte lineage cells in AD and a concomitant decrease in expression of Bin1 in other cell types, including OLs. BIN1 protein expression has been found to decrease in the myelin proteome of Cnp-KO mice compared to WT mice, but not in Plp-KO, Mag-KO, or Sept8-KO mice . Consistent with this previous data, we found a significant down-regulation of BIN1 expression in the corpus callosum in Cnp-KO mice compared to WT mice (p = 0.0041 and p = 0.0007, Fig. 8a-b). We selected to profile the corpus callosum because it is an OL- and myelin-rich tissue, and therefore we reasoned that protein expression in this region more accurately would reflect the effect of Cnp-KO on OLs protein expression. No changes in the subcellular distribution of BIN1 were detected in Cnp-KO mice, which was primarily expressed in the cytoplasm of cells that were also stained by antibodies specific for the nuclear pan-OL lineage marker OLIG2 (Fig. 8c-e). Next, we measured the protein expression of GOT2, which is a mitochondrial protein that has been shown to be down-regulated in many AD gene expression studies , and has down-regulated transcript levels in Cnp-KO mice in our RNAseq expression profiles. Consistent with this data, we found a down-regulation of GOT2 protein levels in Cnp-KO mice (p = 0.0332; Fig. 8a-b). Although GOT2 has been identified in the myelin proteome, we found that it was primarily co-expressed with neurons (Fig. 8f–h), suggesting that the loss of Cnp expression may also affect mitochondrial gene expression within neurons. Taken together, these results validate that Cnp-KO mice have gene dysregulation in key myelin and mitochondrial proteins that is similar to that seen in AD brain samples.
In this study, we employed an unbiased systems biology approach to characterize OL-enriched and AD-associated co-expression and regulatory molecular networks. We derived a core OL-enriched gene set (COLGS) that was highly enriched in AD GWAS genes. We found that this set of genes strongly overlapped with the corresponding genes from a protein co-expression module, which we found to be highly dysregulated in AD. We next constructed a core OL Bayesian regulatory network (COLBN) to dissect the causal relationships among the genes in COLGS. By employing a series of in vitro and in vivo perturbations of key driver genes in COLBN, we validated the predicted network structures in COLBN. We further showed that the knockouts of key drivers of COLBN in mice mimic the dysregulation of OL-associated compartments of genes in human postmortem AD brains. In particular, our data revealed a surprisingly strong, convergent gene expression effect of the knockouts Cnp and Plp1 on organelle-associated gene expression pathways, specifically in genes annotated for mitochondrial and ribosome functions. These organelles have also been reported to be dysregulated in axons in AD brains [17, 67,68,69], suggesting that altered OL-axon communication may lead to dysregulated expression of ribosomal and mitochondrial genes and possibly play a role in contributing to AD axonal pathology. It is difficult to study prodromal changes in late-onset AD brains, because we are not able to determine a priori the individuals who would progress to AD. However, in mice with well-defined ablation of myelin genes also detected as key drivers of human gene network in AD, it is possible to define alterations that occur in OL and that precede frank neurodegeneration. For this reason the detection of similar gene changes in brain tissue samples of Cnp-KO and Plp1-KO mice prior to the development of any axonal pathology [22, 24], allowed us to infer that the dysregulated expression of similar genes in murine samples and in human AD brains is suggestive of similar events occurring during the early part of the pathological cascade. Overall, this study improves our understanding of the molecular underpinnings of myelination and OLs in AD by identifying biologically relevant pathways, dissecting the causal relationships among the OL- and myelin-related genes, and implicating key driver genes in AD pathogenesis.
At the individual gene level, many of the genes in COLGS have been described as genetic risk factors associated with late-onset AD, including BIN1 , PICALM , NME8 , UNC5C , and PSEN . Notably, BIN1 is the nearest protein-coding gene to the SNP with the second-strongest GWAS signal for AD, following APOE . Histologically, BIN1 is primarily found in the brain at the nodes of Ranvier, consistent with its high RNA expression in OLs and its presence in the myelin proteome [43, 70, 71]. In COLBN, BIN1 is downstream of ABCA2, a cholesterol transporter that has been associated with the risk of AD in many study populations [28, 72,73,74,75]. Much of the literature about the role of BIN1 in AD has focused on its roles in neurons and microglia , but our data, in addition to another recent study , suggest that its role in OLs should be explored further. In addition to genetic risk factors, our AD OL network also contains many genes encoding proteins that have been associated with AD pathophysiology (e.g., via Aβ production) including PSEN1 , BACE1 , PLD1 [46, 78], and APLP1 . Consistent with the important role of BACE1 in OLs suggested by our network, BACE1 has been shown to play a key role in myelination [80, 81]. The role of BACE1 in OLs is of high relevance to AD, as mutations in the BACE1-cleaving region of APP have been associated with a decreased risk of AD , and β-secretase inhibitors intended to treat AD may have side-effects of myelin defects . A focus on the interaction targets of both PSEN1 and BACE1 within OLs using regulatory networks may be a fruitful avenue to identify treatment modalities that decrease deleterious Aβ production without causing off-target effects.
At the gene set level, our enrichment analysis of the AD co-expression modules shows that the three modules significantly enriched for AD risk genes are associated with three different cell types, i.e., microglia, OLs, and neurons (Table 1). Note that gene modules were identified based on the correlation between gene expression profiles in postmortem human AD brains, and they don’t necessarily correspond to any particular known cell types or biological processes due to interactions among cell types and biological processes. Therefore, we denote the co-expression modules by randomly selected color names in addition to their most enriched gene ontology term, to emphasize their multifaceted and highly context-dependent functions. The top ranked module was enriched for immune (microglia/macrophage) genes, consistent with recent reports that immune cells and in particular innate immunity plays a critical role in promoting AD [76, 84,85,86]. However, it is imprudent to focus on a single cell type and ignore the interactions with other cells. For example, TREM2, an established AD risk factor that is primarily expressed in immune cells [87, 88], is also the causative gene of Nasu-Hakola disease, an early-onset subcortical dementia that presents with white matter demyelination . Mice lacking Trem2 have been shown to have delayed myelin debris clearance, which may lead to increased microglia activation and thus demyelination and neuronal death . The dysregulation of myelin proteome genes that we observed in AD may contribute to pathologic inflammation, by increasing the available lipid pool for scavenging by microglia, which can activate microglia into a pro-inflammatory state [91, 92]. Further investigation of cell type interactions in AD via network biology is a promising approach in addressing the underlying causes of AD.
Existing mouse models of AD tend to focus on Aβ and/or neuronal deficits in AD. For example, several mouse models of AD express genes with familial AD-causative mutations under the Thy1 promoter [93,94,95], which is a neuronal marker and will serve to restrict the pathologic changes to neurons. However, the data presented in this study and others suggest that the dysregulation of other cell types, including OLs, may play a role in AD. This opens up a need for mouse models of AD that can recapitulate the OL- and myelin-associated dysregulation in AD. In this study, we found that a mouse knockout model of one of the key drivers in the OL network, Cnp, demonstrates a strikingly similar myelin and mitochondrial dysregulation pattern as is seen in brain samples of patients with AD. Notably, we did expression profiling prior to the onset of axon degeneration in Cnp-KO mice , to minimize the possibility that the myelin gene expression changes observed would be reactive to as opposed to premonitory for axon damage. Taken together, our data suggest that the Cnp-KO mice may be a good model of the OL network gene dysregulation and dysmyelination that occurs in the brains of patients with AD. Therefore, therapeutic agents that are able to mitigate and/or prevent the dysmyelination and axon degeneration seen in Cnp-KO mice are worthwhile of investigation as potential therapeutic agents for ameliorating cognitive deficits in patients with AD.
In this study, we focused on the molecular networks in AD in a brain cell type, OLs, which have not been widely studied in AD. Our network modeling approach uncovered a network of OL-associated genes that is enriched for AD GWAS genes, pathways through which dysregulation of this OL network may promote AD pathology, and key driver genes that orchestrate these pathways. Further work on our model for the role of OLs in AD may help to address why aging is the major risk factor for AD, since myelin maintenance and plasticity are known to become progressively less robust in normal aging [62, 96]. In particular, preventing or reversing the dysregulation of key OL driver genes such as CNP and downstream targets such as BIN1 deserve further research as interventions to help to alleviate the progression of cognitive deficits in AD.
This study systematically identifies and validates a comprehensive molecular blueprint of OLs in the context of AD. Our gene network analysis of large-scale genetic and genomic data from AD brains reveals that OL/myelin-enriched subnetworks in multiple brain regions are strongly associated with clinical and neuropathologic traits in AD. These OL/myelin-enriched subnetworks harbor not only Aβ production-related genes, but also several genes involved in myelin biology. Our network analysis further identifies key causal regulators of the OL/myelin-enriched networks, including UGT8, CNP, MYRF, and PLP1. These predicted network structures can be validated by gene perturbation signatures of these drivers. Mice with genetic ablations of Cnp mimicked components of organelle and myelin sheath gene expression dysregulation seen in brain samples from patients with AD, including decreased protein expression of BIN1 and GOT2. Overall, our network models of OLs in AD and the comprehensive validation experiments reveal details of the molecular mechanisms of OL dysregulation in AD and thus pave a way for the development of novel AD therapies.
Central nervous system
Core oligodendrocyte-enriched Bayesian network
Core oligodendrocyte gene set
Differentially expressed gene
Expression single nucleotide polymorphisms
False discovery rate
Fisher’s exact test
Gene expression omnibus
Genome-wide association study
Harvard brain tissue resource center
Hugo Gene nomenclature committee
International genomics of Alzheimer’s Project
Mount Sinai Brain Bank
RNA integrity number
Weighted gene coexpression network analysis
Hebert LE, Weuve J, Scherr PA, Evans DA. Alzheimer disease in the United States (2010-2050) estimated using the 2010 census. Neurology. 2013;80:1778–83.
Plassman BL, Langa KM, Fisher GG, Heeringa SG, Weir DR, Ofstedal MB, Burke JR, Hurd MD, Potter GG, Rodgers WL, et al. Prevalence of dementia in the United States: the aging, demographics, and memory study. Neuroepidemiology. 2007;29:125–32.
Jack CR, Wiste HJ, Weigand SD, Knopman DS, Vemuri P, Mielke MM, Lowe V, Senjem ML, Gunter JL, Machulda MM, et al. Age, sex, and APOE ε4 effects on memory, brain structure, and β-Amyloid across the adult life span. JAMA Neurol. 2015;72:511–9.
Veeraraghavalu K, Zhang C, Zhang X, Tanzi RE, Sisodia SS. Age-dependent, non-cell-autonomous deposition of amyloid from synthesis of β-amyloid by cells other than excitatory neurons. J Neurosci. 2014;34:3668–73.
Li C, Zhao R, Gao K, Wei Z, Yin MY, Lau LT, Chui D, Hoi Yu AC. Astrocytes: implications for neuroinflammatory pathogenesis of Alzheimer’s disease. Curr Alzheimer Res. 2011;8:67–80.
Griciuc A, Serrano-Pozo A, Parrado AR, Lesinski AN, Asselin CN, Mullin K, Hooli B, Choi SH, Hyman BT, Tanzi RE. Alzheimer’s disease risk gene CD33 inhibits microglial uptake of amyloid beta. Neuron. 2013;78:631–43.
Salat DH, Greve DN, Pacheco JL, Quinn BT, Helmer KG, Buckner RL, Fischl B. Regional white matter volume differences in nondemented aging and Alzheimer’s disease. NeuroImage. 2009;44:1247–58.
Cash DM, Ridgway GR, Liang Y, Ryan NS, Kinnunen KM, Yeatman T, Malone IB, Benzinger TLS, Jack CR, Thompson PM, et al. The pattern of atrophy in familial Alzheimer disease: volumetric MRI results from the DIAN study. Neurology. 2013;81:1425–33.
Salat DH, Kaye JA, Janowsky JS. Prefrontal gray and white matter volumes in healthy aging and Alzheimer disease. Arch Neurol. 1999;56:338–44.
Hua X, Hibar DP, Lee S, Toga AW, Jack CR, Weiner MW, Thompson PM. Sex and age differences in atrophic rates: an ADNI study with n=1368 MRI scans. Neurobiol Aging. 2010;31:1463–80.
Lee GJ, Lu PH, Hua X, Lee S, Wu S, Nguyen K, Teng E, Leow AD, Jack CR, Toga AW, et al. Depressive symptoms in mild cognitive impairment predict greater atrophy in Alzheimer’s disease-related regions. Biol Psychiatry. 2012;71:814–21.
Lee GJ, Lu PH, Medina LD, Rodriguez-Agudelo Y, Melchor S, Coppola G, Braskie MN, Hua X, Apostolova LG, Leow AD, et al. Regional brain volume differences in symptomatic and presymptomatic carriers of familial Alzheimer’s disease mutations. J Neurol Neurosurg Psychiatry. 2013;84:154–62.
Braak H, Braak E. Development of Alzheimer-related neurofibrillary changes in the neocortex inversely recapitulates cortical myelogenesis. Acta Neuropathol. 1996;92:197–201.
Braak H, Del Tredici K. The preclinical phase of the pathological process underlying sporadic Alzheimer’s disease. Brain. 2015;138((Pt 10)):2814–33.
Kim JY, Shen S, Dietz K, He Y, Howell O, Reynolds R, Casaccia P. HDAC1 nuclear export induced by pathological conditions is essential for the onset of axonal damage. Nat Neurosci. 2010;13:180–9.
Werner HB, Kuhlmann K, Shen S, Uecker M, Schardt A, Dimova K, Orfaniotou F, Dhaunchak A, Brinkmann BG, Möbius W, et al. Proteolipid protein is required for transport of sirtuin 2 into CNS myelin. J Neurosci. 2007;27:7717–30.
Nave K-A. Myelination and the trophic support of long axons. Nat Rev Neurosci. 2010;11:275–83.
Nave K-A, Werner HB. Myelination of the nervous system: mechanisms and functions. Annu Rev Cell Dev Biol. 2014;30:503–33.
Vossel KA, Zhang K, Brodbeck J, Daub AC, Sharma P, Finkbeiner S, Cui B, Mucke L. Tau reduction prevents Abeta-induced defects in axonal transport. Science (New York, NY). 2010;330:198.
Zhang B, Maiti A, Shively S, Lakhani F, McDonald-Jones G, Bruce J, Lee EB, Xie SX, Joyce S, Li C, et al. Microtubule-binding drugs offset tau sequestration by stabilizing microtubules and reversing fast axonal transport deficits in a tauopathy model. Proc Natl Acad Sci U S A. 2005;102:227–31.
Coetzee T, Fujita N, Dupree J, Shi R, Blight A, Suzuki K, Popko B. Myelination in the absence of galactocerebroside and sulfatide: normal structure with abnormal function and regional instability. Cell. 1996;86:209–19.
Edgar JM, McLaughlin M, Werner HB, McCulloch MC, Barrie JA, Brown A, Faichney AB, Snaidero N, Nave K-A, Griffiths IR. Early ultrastructural defects of axons and axon-glia junctions in mice lacking expression of Cnp1. Glia. 2009;57:1815–24.
Lappe-Siefke C, Goebbels S, Gravel M, Nicksch E, Lee J, Braun PE, Griffiths IR, Nave K-A. Disruption of Cnp1 uncouples oligodendroglial functions in axonal support and myelination. Nat Genet. 2003;33:366–74.
Edgar JM, McLaughlin M, Yool D, Zhang SC, Fowler JH, Montague P, Barrie JA, McCulloch MC, Duncan ID, Garbern J, et al. Oligodendroglial modulation of fast axonal transport in a mouse model of hereditary spastic paraplegia. J Cell Biol. 2004;166(1):121–31.
Griffiths I, Klugmann M, Anderson T, Yool D, Thomson C, Schwab MH, Schneider A, Zimmermann F, McCulloch M, Nadon N, et al. Axonal swellings and degeneration in mice lacking the major proteolipid of myelin. Science. 1998;280(5369):1610–3.
Zhang B, Gaiteri C, Bodea L-G, Wang Z, McElwee J, Podtelezhnikov AA, Zhang C, Xie T, Tran L, Dobrin R, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer’s disease. Cell. 2013;153:707–20.
Cruchaga C, Haller G, Chakraverty S, Mayo K, Vallania FLM, Mitra RD, Faber K, Williamson J, Bird T, Diaz-Arrastia R, et al. Rare variants in APP, PSEN1 and PSEN2 increase risk for AD in late-onset Alzheimer’s disease families. PLoS One. 2012;7:e31039.
Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, DeStafano AL, Bis JC, Beecham GW, Grenier-Boley B, et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45:1452–8.
Ahmed RR, Holler CJ, Webb RL, Li F, Beckett TL, Murphy MP. BACE1 and BACE2 enzymatic activities in Alzheimer’s disease. J Neurochem. 2010;112:1045–53.
Wetzel-Smith MK, Hunkapiller J, Bhangale TR, Srinivasan K, Maloney JA, Atwal JK, Sa SM, Yaylaoglu MB, Foreman O, Ortmann W, et al. A rare mutation in UNC5C predisposes to late-onset Alzheimer’s disease and increases neuronal cell death. Nat Med. 2014;20:1452–7.
Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17.
Zhang Y, Chen K, Sloan SA, Bennett ML, Scholze AR, O’Keeffe S, Phatnani HP, Guarnieri P, Caneda C, Ruderisch N, et al. An RNA-sequencing Transcriptome and splicing database of Glia, neurons, and vascular cells of the cerebral cortex. J Neurosci. 2014;34:11929–47.
Zhu J, Wiener MC, Zhang C, Fridman A, Minch E, Lum PY, Sachs JR, Schadt EE. Increasing the power to detect causal associations by combining genotypic and expression data in segregating populations. PLoS Comput Biol. 2007;3(4):e69.
Zhu J, Zhang B, Smith EN, Drees B, Brem RB, Kruglyak L, Bumgarner RE, Schadt EE. Integrating large-scale functional genomic data to dissect the complexity of yeast regulatory networks. Nat Genet. 2008;40:854–61.
Wang I-M, Zhang B, Yang X, Zhu J, Stepaniants S, Zhang C, Meng Q, Peters M, He Y, Ni C, et al. Systems Analysis of Eleven Rodent Disease Models Reveals an Inflammatome Signature and Key Drivers. Mol Syst Biol 2012, (in press).
Mishra A, Macgregor S. VEGAS2: software for more flexible gene-based testing. Twin Res Hum Genet. 2015;18:86–91.
McKenzie AT, Katsyv I, Song WM, Wang M, Zhang B. DGCA: a comprehensive R package for differential gene correlation analysis. BMC Syst Biol. 2016;10(1):106.
Storey J. Qvalue: Q-value estimation for false discovery rate control. 2015.
Anders S, Reyes A, Huber W. Detecting differential usage of exons from RNA-seq data. Genome Res. 2012;22:2008–17.
Blalock EM, Geddes JW, Chen KC, Porter NM, Markesbery WR, Landfield PW. Incipient Alzheimer's disease: microarray correlation analyses reveal major transcriptional and tumor suppressor responses. Proc Natl Acad Sci U S A. 2004;101:2173–8.
Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics (Oxford, England). 2007;23:257–8.
Wang M, Zhao Y, Zhang B. Efficient test and visualization of multi-set intersections. Sci Rep. 2015;5:16923.
Dagley LF, White CA, Liao Y, Shi W, Smyth GK, Orian JM, Emili A, Purcell AW. Quantitative proteomic profiling reveals novel region-specific markers in the adult mouse brain. Proteomics. 2014;14:241–61.
Miller JA, Oldham MC, Geschwind DH. A systems level analysis of transcriptional changes in Alzheimer’s disease and normal aging. J Neurosci. 2008;28:1410–20.
Emery B, Agalliu D, Cahoy JD, Watkins TA, Dugas JC, Mulinyawe SB, Ibrahim A, Ligon KL, Rowitch DH, Barres BA. Myelin gene regulatory factor is a critical transcriptional regulator required for CNS Myelination. Cell. 2009;138:172–85.
Cai D, Zhong M, Wang R, Netzer WJ, Shields D, Zheng H, Sisodia SS, Foster DA, Gorelick FS, Xu H, et al. Phospholipase D1 corrects impaired betaAPP trafficking and neurite outgrowth in familial Alzheimer's disease-linked presenilin-1 mutant neurons. Proc Natl Acad Sci U S A. 2006;103:1936–40.
Jin J-K, Ahn B-H, Na Y-J, Kim J-I, Kim Y-S, Choi E-K, Ko Y-G, Chung KC, Kozlowski PB, Min DS. Phospholipase D1 is associated with amyloid precursor protein in Alzheimer’s disease. Neurobiol Aging. 2007;28:1015–27.
Bayer TA, Paliga K, Weggen S, Wiestler OD, Beyreuther K, Multhaup G. Amyloid precursor-like protein 1 accumulates in neuritic plaques in Alzheimer’s disease. Acta Neuropathol. 1997;94:519–24.
Klugmann M, Schwab MH, Pühlhofer A, Schneider A, Zimmermann F, Griffiths IR, Nave KA. Assembly of CNS myelin in the absence of proteolipid protein. Neuron. 1997;18:59–70.
Morell P, Quarles RH. Characteristic Composition of Myelin. Philadelphia: Lippincott-Raven; 1999.
von Trotha KT, Heun R, Schmitz S, Lutjohann D, Maier W, Kolsch H. Influence of lysosomal acid lipase polymorphisms on chromosome 10 on the risk of Alzheimer’s disease and cholesterol metabolism. Neurosci Lett. 2006;402(3):262–6.
Lee J. Process outgrowth in oligodendrocytes is mediated by CNP, a novel microtubule assembly myelin protein. J Cell Biol. 2005;170:661–73.
Miller JA, Woltjer RL, Goodenbour JM, Horvath S, Geschwind DH. Genes and pathways underlying regional and cell type changes in Alzheimer’s disease. Genome Med. 2013;5(5):48.
de Monasterio-Schrader P, Jahn O, Tenzer S, Wichert SP, Patzig J, Werner HB. Systematic approaches to central nervous system myelin. Cell Mol Life Sci. 2012;69(17):2879–94.
Harlow DE, Saul KE, Culp CM, Vesely EM, Macklin WB. Expression of proteolipid protein gene in spinal cord stem cells and early oligodendrocyte progenitor cells is dispensable for normal cell migration and myelination. J Neurosci. 2014;34:1333–43.
Garbern JY, Yool DA, Moore GJ, Wilds IB, Faulk MW, Klugmann M, Nave K-A, Sistermans EA, van der Knaap MS, Bird TD, et al. Patients lacking the major CNS myelin protein, proteolipid protein 1, develop length-dependent axonal degeneration in the absence of demyelination and inflammation. Brain. 2002;125:551–61.
Mohan H, Friese A, Albrecht S, Krumbholz M, Elliott CL, Arthur A, Menon R, Farina C, Junker A, Stadelmann C, et al. Transcript profiling of different types of multiple sclerosis lesions yields FGF1 as a promoter of remyelination. Acta Neuropathol Commun. 2014;2:168.
Mashayekhi F, Hadavi M, Vaziri HR, Naji M. Increased acidic fibroblast growth factor concentrations in the serum and cerebrospinal fluid of patients with Alzheimer’s disease. J Clin Neurosci. 2010;17(3):357–9.
Tamiya G, Makino S, Hayashi M, Abe A, Numakura C, Ueki M, Tanaka A, Ito C, Toshimori K, Ogawa N, et al. A mutation of COX6A1 causes a recessive axonal or mixed form of Charcot-Marie-tooth disease. Am J Hum Genet. 2014;95:294–300.
Piedrahita D, Hernández I, López-Tobón A, Fedorov D, Obara B, Manjunath BS, Boudreau RL, Davidson B, Laferla F, Gallego-Gómez JC, et al. Silencing of CDK5 reduces neurofibrillary tangles in transgenic alzheimer’s mice. J Neurosci. 2010;30:13966–76.
Yang Y, Wang H, Zhang J, Luo F, Herrup K, Bibb JA, Lu R, Miller RH. Cyclin dependent kinase 5 is required for the normal development of oligodendrocytes and myelin formation. Dev Biol. 2013;378:94–106.
Bartzokis G. Alzheimer’s disease as homeostatic responses to age-related myelin breakdown. Neurobiol Aging. 2011;32:1341–71.
De Rossi P, Buggia-Prevot V, Clayton BL, Vasquez JB, van Sanford C, Andrew RJ, Lesnick R, Botte A, Deyts C, Salem S, et al. Predominant expression of Alzheimer’s disease-associated BIN1 in mature oligodendrocytes and localization to white matter tracts. Mol Neurodegener. 2016;11(1):59.
Holler CJ, Davis PR, Beckett TL, Platt TL, Webb RL, Head E, Murphy MP. Bridging integrator 1 (BIN1) protein expression increases in the Alzheimer’s disease brain and correlates with neurofibrillary tangle pathology. J Alzheimers Dis. 2014;42(4):1221–7.
Patzig J, Erwig MS, Tenzer S, Kusch K, Dibaj P, Mobius W, Goebbels S, Schaeren-Wiemers N, Nave KA, Werner HB. Septin/anillin filaments scaffold central nervous system myelin to accelerate nerve conduction. elife. 2016;5:e17119.
Puthiyedth N, Riveros C, Berretta R, Moscato P. Identification of differentially expressed genes through integrated study of Alzheimer’s disease affected brain regions. PLoS One. 2016;11(4):e0152342.
Calkins MJ, Manczak M, Mao P, Shirendeb U, Reddy PH. Impaired mitochondrial biogenesis, defective axonal transport of mitochondria, abnormal mitochondrial dynamics and synaptic degeneration in a mouse model of Alzheimer’s disease. Hum Mol Genet. 2011;20:4515–29.
Iijima-Ando K, Sekiya M, Maruko-Otake A, Ohtake Y, Suzuki E, Lu B, Iijima KM. Loss of axonal mitochondria promotes tau-mediated neurodegeneration and Alzheimer’s disease-related tau phosphorylation via PAR-1. PLoS Genet. 2012;8:e1002918.
González C, Couve A. The axonal endoplasmic reticulum and protein trafficking: cellular bootlegging south of the soma. Semin Cell Dev Biol. 2014;27:23–31.
Butler MH, David C, Ochoa GC, Freyberg Z, Daniell L, Grabs D, Cremona O, De Camilli P. Amphiphysin II (SH3P9; BIN1), a member of the amphiphysin/Rvs family, is concentrated in the cortical cytomatrix of axon initial segments and nodes of ranvier in brain and around T tubules in skeletal muscle. J Cell Biol. 1997;137:1355–67.
Zhang Y, Sloan Steven A, Clarke Laura E, Caneda C, Plaza Colton A, Blumenthal Paul D, Vogel H, Steinberg Gary K, Edwards Michael SB, Li G, et al. Purification and characterization of progenitor and mature human Astrocytes reveals transcriptional and functional differences with mouse. Neuron. 2015;89:37–53.
Macé S, Cousin E, Ricard S, Génin E, Spanakis E, Lafargue-Soubigou C, Génin B, Fournel R, Roche S, Haussy G, et al. ABCA2 is a strong genetic risk factor for early-onset Alzheimer’s disease. Neurobiol Dis. 2005;18:119–25.
Wollmer MA, Kapaki E, Hersberger M, Muntwyler J, Brunner F, Tsolaki M, Akatsu H, Kosaka K, Michikawa M, Molyva D, et al. Ethnicity-dependent genetic association of ABCA2 with sporadic Alzheimer’s disease. Am J Med Genet B Neuropsychiatr Genet. 2006;141B:534–6.
Xu X, Wang Y, Wang L, Liao Q, Chang L, Xu L, Huang Y, Ye H, Xu L, Chen C, et al. Meta-analyses of 8 polymorphisms associated with the risk of the Alzheimer’s disease. PLoS One. 2013;8:e73129.
Minster RL, DeKosky ST, Kamboh MI. No association of DAPK1 and ABCA2 SNPs on chromosome 9 with Alzheimer’s disease. Neurobiol Aging. 2009;30:1890–1.
Karch CM, Goate AM. Alzheimer’s disease risk genes and mechanisms of disease pathogenesis. Biol Psychiatry. 2015;77:43–51.
Karran E, Mercken M, De Strooper B. The amyloid cascade hypothesis for Alzheimer’s disease: an appraisal for the development of therapeutics. Nat Rev Drug Discov. 2011;10(9):698–712.
Liu Y, Zhang Y-W, Wang X, Zhang H, You X, Liao F-F, Xu H. Intracellular trafficking of presenilin 1 is regulated by beta-amyloid precursor protein and phospholipase D1. J Biol Chem. 2009;284:12145–52.
Zhou F, van Laar T, Huang H, Zhang L. APP and APLP1 are degraded through autophagy in response to proteasome inhibition in neuronal cells. Protein & cell. 2011;2:377–83.
Treiber H, Hagemeyer N, Ehrenreich H, Simons M. BACE1 in central nervous system myelination revisited. Mol Psychiatry. 2012;17:237–9.
Hu X, Hicks CW, He W, Wong P, Macklin WB, Trapp BD, Yan R. Bace1 modulates myelination in the central and peripheral nervous system. Nat Neurosci. 2006;9:1520–5.
Jonsson T, Atwal JK, Steinberg S, Snaedal J, Jonsson PV, Bjornsson S, Stefansson H, Sulem P, Gudbjartsson D, Maloney J, et al. A mutation in APP protects against Alzheimer’s disease and age-related cognitive decline. Nature. 2012;488:96–9.
Willem M, Garratt AN, Novak B, Citron M, Kaufmann S, Rittger A, DeStrooper B, Saftig P, Birchmeier C, Haass C. Control of peripheral nerve myelination by the beta-secretase BACE1. Science (New York, NY). 2006;314:664–6.
Miller JA, Woltjer RL, Goodenbour JM, Horvath S, Geschwind DH. Genes and pathways underlying regional and cell type changes in Alzheimer’s disease. Genome Med. 2013;5:48.
Jones L, Lambert J-C, Wang L-S, Choi S-H, Harold D, Vedernikov A, Escott-Price V, Stone T, Richards A, Bellenguez C, et al. Convergent genetic and expression data implicate immunity in Alzheimer’s disease. Alzheimers Dement. 2015;11(6):658–71.
Seyfried NT, Dammer EB, Swarup V, Nandakumar D, Duong DM, Yin L, Deng Q, Nguyen T, Hales CM, Wingo T, et al. A multi-network approach identifies protein-specific co-expression in asymptomatic and symptomatic Alzheimer’s disease. Cell Syst. 2017;4(1):60–72 e64.
Jin SC, Benitez BA, Karch CM, Cooper B, Skorupa T, Carrell D, Norton JB, Hsu S, Harari O, Cai Y, et al. Coding variants in TREM2 increase risk for Alzheimer’s disease. Hum Mol Genet. 2014;23:5838–46.
Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, Cruchaga C, Sassi C, Kauwe JSK, Younkin S, et al. TREM2 variants in Alzheimer’s disease. N Engl J Med. 2013;368:117–27.
Paloneva J, Manninen T, Christman G, Hovanes K, Mandelin J, Adolfsson R, Bianchin M, Bird T, Miranda R, Salmaggi A, et al. Mutations in two genes encoding different subunits of a receptor signaling complex result in an identical disease phenotype. Am J Hum Genet. 2002;71:656–62.
Poliani PL, Wang Y, Fontana E, Robinette ML, Yamanishi Y, Gilfillan S, Colonna M. TREM2 sustains microglial expansion during aging and response to demyelination. J Clin Invest. 2015;125:2161–70.
Painter MM, Atagi Y, Liu CC, Rademakers R, Xu H, Fryer JD, Bu G. TREM2 in CNS homeostasis and neurodegenerative disease. Mol Neurodegener. 2015;10:43.
Cherry JD, Olschowka JA, O’Banion MK. Neuroinflammation and M2 microglia: the good, the bad, and the inflamed. J Neuroinflammation. 2014;11:98.
Oakley H, Cole SL, Logan S, Maus E, Shao P, Craft J, Guillozet-Bongaarts A, Ohno M, Disterhoft J, Van Eldik L, et al. Intraneuronal beta-amyloid aggregates, neurodegeneration, and neuron loss in transgenic mice with five familial Alzheimer's disease mutations: potential factors in amyloid plaque formation. J Neurosci. 2006;26(40):10129–40.
Radde R, Bolmont T, Kaeser SA, Coomaraswamy J, Lindau D, Stoltze L, Calhoun ME, Jaggi F, Wolburg H, Gengler S, et al. Abeta42-driven cerebral amyloidosis in transgenic mice reveals early and robust pathology. EMBO Rep. 2006;7(9):940–6.
Sturchler-Pierrat C, Abramowski D, Duke M, Wiederhold KH, Mistl C, Rothacher S, Ledermann B, Burki K, Frey P, Paganetti PA, et al. Two amyloid precursor protein transgenic mouse models with Alzheimer disease-like pathology. Proc Natl Acad Sci U S A. 1997;94(24):13287–92.
Krestinina O, Azarashvili T, Baburina Y, Galvita A, Grachev D, Stricker R, Reiser G. In aging, the vulnerability of rat brain mitochondria is enhanced due to reduced level of 2′,3′-cyclic nucleotide-3′-phosphodiesterase (CNP) and subsequently increased permeability transition in brain mitochondria in old animals. Neurochem Int. 2015;80:41–50.
Krämer-Albers E-M, Bretz N, Tenzer S, Winterstein C, Möbius W, Berger H, Nave K-A, Schild H, Trotter J. Oligodendrocytes secrete exosomes containing major myelin and stress-protective proteins: Trophic support for axons? Proteomics Clin Appl. 2007;1:1446–61.
Subramaniam S, Fahy E, Gupta S, Sud M, Byrnes RW, Cotter D, Dinasarapu AR, Maurya MR. Bioinformatics and systems biology of the lipidome. Chem Rev. 2011;111:6452–90.
We would like to thank members of the Zhang and Casaccia labs for many fruitful discussions.
This work was supported by the NIH grants R01AG046170, U01AG052411, RF1AG054014, RF1AG057440, R01AG057907, U01AI111598–01, R01NS067550, P50AG025688, U01AG046161, and F30AG052261. R01AG046170 is a component of the AMP-AD Target Discovery and Preclinical Validation Project. We thank the International Genomics of Alzheimer’s Project (IGAP) for providing summary results data for these analyses. The investigators within IGAP contributed to the design and implementation of IGAP and/or provided data but did not participate in analysis or writing of this report. IGAP was made possible by the generous participation of the control subjects, the patients, and their families. The i–Select chips was funded by the French National Foundation on Alzheimer’s disease and related disorders. EADI was supported by the LABEX (laboratory of excellence program investment for the future) DISTALZ grant, Inserm, Institut Pasteur de Lille, Université de Lille 2 and the Lille University Hospital. GERAD was supported by the Medical Research Council (Grant n° 503,480), Alzheimer’s Research UK (Grant n° 503,176), the Wellcome Trust (Grant n° 082604/2/07/Z) and German Federal Ministry of Education and Research (BMBF): Competence Network Dementia (CND) grant n° 01GI0102, 01GI0711, 01GI0420. CHARGE was partly supported by the NIH/NIA grant R01 AG033193 and the NIA AG081220 and AGES contract N01–AG–12100, the NHLBI grant R01 HL105756, the Icelandic Heart Association, and the Erasmus Medical Center and Erasmus University. ADGC was supported by the NIH/NIA grants: U01 AG032984, U24 AG021886, U01 AG016976, and the Alzheimer’s Association grant ADGC–10–196,728.
Availability of data and materials
The RNA-sequencing data from the mouse key driver knockout experiments are available at Gene Expression Omnibus (GEO) GSE80437. All other relevant data is described elsewhere and available from the authors upon request.
Use of animals in this research was strictly compliant with the guidelines set forth by the US Public Health Service in their policy on Humane Care and Use of Laboratory Animals, and in the Guide for the Care and Use of Laboratory Animals. All animal procedures received prior approval from the Institutional Animal Care and Use Committee at Icahn School of Medicine at Mount Sinai (IACUC 08–676).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sarah Moyon, Patrizia Casaccia and Bin Zhang are Co-senior authors.
Supplementary Experimental Procedures. Figure S1. The topological overlap matrix plot of the protein co-expression network constructed from the proteomics data from the autopsied brains in the MSBB cohort, along with the dendrogram showing the tree cutting process used to define modules (above). Figure S2. Confirmation that the key driver knockouts abrogate gene expression of the key driver in the RNA-seq experiments. For each of the key driver knockouts whose genome-wide gene expression was profiled in this study using RNA-seq, we plotted the log10 counts overlapping that gene in both the wildtype (WT) and knockout (KO) samples. Notably, one of the matched samples from Cnp was detected as an outlier in both the CBM and FC brain regions (red), due to suspected mislabeling. These samples were removed prior to downstream differential expression analysis. (DOCX 161 kb)
Module membership file for the Mount Sinai Brain Bank (MSBB) proteomic coexpression network. The module label, a randomly chosen color name, is in the 1st column, while the protein name is in the 2nd column. (TSV 34 kb)
Differentially expressed proteins in the MSBB proteomics data set between AD cases and controls in the OL-enriched module (which has been randomly assigned the color name “Yellow”). (TSV 9 kb)
Summary of read mapping from the three knockout mouse RNAseq experiments generated by TopHat. (XLSX 52 kb)
Differentially expressed genes found the mouse knockout RNAseq analyses of Ugt8 in the CBM (Data 1) and FC (Data 2), Cnp in the FC (Data 3), and Plp1 in the CBM (Data 4). For these differential expression analyses, we mapped RNAseq reads using TopHat, converted to count space using HTSeq, used voom to transform the read space data to log2 counts per million, and used limma for differential expression analysis. We also used the Ensembl database to identify the human gene with the highest homology percentage based on protein coding region DNA divergence, and report this homology percentage for each gene. Note that the differential expression signatures of Cnp in the CBM and Plp1 in the FC were not found not have any differentially expressed genes at FDR < 0.3, so they are not included here. (ZIP 358 kb)
About this article
- Alzheimer’s disease
- co-expression network
- Causal network
- RNA sequencing
- Differential expression