Skip to main content

Umbrella menu

  • SfN.org
  • eNeuro
  • The Journal of Neuroscience
  • Neuronline
  • BrainFacts.org

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Latest Articles
    • Issue Archive
    • Editorials
    • Research Highlights
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • EDITORIAL BOARD
  • BLOG
  • ABOUT
    • Overview
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SfN.org
  • eNeuro
  • The Journal of Neuroscience
  • Neuronline
  • BrainFacts.org

User menu

  • My alerts

Search

  • Advanced search
eNeuro
  • My alerts

eNeuro

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Latest Articles
    • Issue Archive
    • Editorials
    • Research Highlights
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • EDITORIAL BOARD
  • BLOG
  • ABOUT
    • Overview
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
PreviousNext
Research ArticleNew Research, Disorders of the Nervous System

Mega-Analysis of Gene Expression in Mouse Models of Alzheimer’s Disease

Beryl Zhuang, B. Ogan Mancarci, Lilah Toker and Paul Pavlidis
eNeuro 25 November 2019, 6 (6) ENEURO.0226-19.2019; DOI: https://doi.org/10.1523/ENEURO.0226-19.2019
Beryl Zhuang
1Graduate Program in Bioinformatics, University of British Columbia, Vancouver V6T1Z4, Canada
2Department of Psychiatry, University of British Columbia, Vancouver V6T1Z4, Canada
3Michael Smith Laboratories, University of British Columbia, Vancouver V6T1Z4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
B. Ogan Mancarci
1Graduate Program in Bioinformatics, University of British Columbia, Vancouver V6T1Z4, Canada
2Department of Psychiatry, University of British Columbia, Vancouver V6T1Z4, Canada
3Michael Smith Laboratories, University of British Columbia, Vancouver V6T1Z4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for B. Ogan Mancarci
Lilah Toker
2Department of Psychiatry, University of British Columbia, Vancouver V6T1Z4, Canada
3Michael Smith Laboratories, University of British Columbia, Vancouver V6T1Z4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul Pavlidis
2Department of Psychiatry, University of British Columbia, Vancouver V6T1Z4, Canada
3Michael Smith Laboratories, University of British Columbia, Vancouver V6T1Z4, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Paul Pavlidis
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

While multiple studies have been conducted of gene expression in mouse models of Alzheimer’s disease (AD), their findings have not reached a clear consensus and have not accounted for the potentially confounding effects of changes in cellular composition. To help address this gap, we conducted a re-analysis based meta-analysis (mega-analysis) of ten independent studies of hippocampal gene expression in mouse models of AD. We used estimates of cellular composition as covariates in statistical models aimed to identify genes differentially expressed (DE) at either early or late stages of progression. Our analysis revealed changes in gene expression at early phases shared across studies, including dysregulation of genes involved in cholesterol biosynthesis and the complement system. Expression changes at later stages were dominated by cellular compositional effects. Thus, despite the considerable heterogeneity of the mouse models, we identified common patterns that may contribute to our understanding of AD etiology. Our work also highlights the importance of controlling for cellular composition effects in genomics studies of neurodegeneration.

  • Alzheimer’s
  • bioinformatics
  • genomics
  • metaanalysis
  • transcriptome

Significance Statement

A molecular understanding of Alzheimer’s disease (AD) is important to the development of treatments. Because of the difficulty of studying brain tissue in humans, especially at very early stages of progression, many groups have performed transcriptomic studies of rodent models. Our study identified changes in gene expression that are strikingly consistent across multiple such studies, providing substantial insight into molecular changes present at early stages of these disease models and which may be of importance in the human condition. Our study also demonstrates the importance of accounting for the complex effects of neurodegeneration on brain tissue in data analysis.

Introduction

There is intense interest in understanding the molecular mechanisms that contribute to Alzheimer’s disease (AD), which involve complex interplays of genetic and environmental factors. However, the early changes in the brain before the onset of cognitive impairment in AD are still poorly understood (Masters et al., 2015). In AD patients, pathologic changes in the brain often precede the occurrence of clinical symptoms by years (Pennanen et al., 2004; Webster et al., 2014; Wirz et al., 2014). Postmortem brain samples of AD patients represent a late phase of the disorder, and it is challenging to distinguish causes from effects of neurodegeneration. Studying early molecular changes in mouse models improves our understanding of human disease pathogenesis. Various transgenic mouse models have been developed based on known mutations in familial AD cases and other AD related genes (Webster et al., 2014). Despite the differences in mechanisms, many of these animal models show similar age-dependent progression in decline of cognitive functions that are similar to AD symptoms observed in human (Webster et al., 2014). Identifying genes that are differentially expressed (DE) during early disease phase in the mouse models versus wild type could help understanding early responses to disease initiation, which may be neuroprotective or contribute to disease progression.

Recent expression profiling studies of AD mouse models have reported up-regulation of genes in classical complement cascade and changes of synapse-related genes in early AD disease progression (D’Onofrio et al., 2011; Hong et al., 2013, 2016). Studies of AD mouse models in later phases of progression (Matarin et al., 2015; Saura et al., 2015) often report expression changes of genes related to inflammation and the immune system. However, their findings are not completely concordant. Differences in mouse models can be one of the contributing factors, but differences in experiment design, sample size, microarray platforms, and data processing methods can also influence the results. Since these mouse models develop similar phenotypes, they could share molecular commonalities. A combined analysis of gene expression profiles of mouse models (which we refer to as a mega-analysis, to distinguish it from a meta-analysis, which traditionally refers to pooling of reported findings from studies rather than re-analyzing the raw data) has the potential to identify cross-mouse-model and disease-phase-specific transcriptional changes.

One of the challenges of bulk tissue transcriptome profiling in AD studies is that samples represent a weighted average of cellular composition, and cell counts of different cell types are typically not directly assessed. As a neurodegenerative disorder (ND), AD is characterized by neuronal loss and neuroinflammation (Serrano-Pozo et al., 2011). In AD patients and in some mouse models, loss of pyramidal neurons and other types of neurons in the hippocampus contributes to neuronal loss (Saxena and Caroni, 2011; Serrano-Pozo et al., 2011). Cell-type proportion changes have been identified as one of the driving forces for expression changes in bulk tissues of an AD mouse model (Srinivasan et al., 2016). Therefore, both cellular composition differences and cell type-specific expression differences can contribute to the overall expression changes in bulk tissue samples between disease models and controls, especially at late stages of progression. In contrast to late stages, early effects are expected to be relatively subtle but perhaps more informative about processes leading to pathology.

A recent study examined expression changes of neuronal-specific marker genes in whole brain tissues of human and mouse models of AD and revealed that expression profiles reflect cell population changes (Hokama et al., 2014). To address the cellular composition problem, methods for estimating cell-type proportions from bulk tissue transcriptome profiles have been developed, sometimes referred to as cell-type deconvolution. A commonly used method to deconvolute effects of cell-type proportion changes of bulk tissue is to estimate cell-type proportions by cell type-specific marker genes (Gaujoux and Seoighe, 2013; Shen-Orr and Gaujoux, 2013; Chikina et al., 2015; Mancarci et al., 2017). While not without caveats, accounting for cell type-specific expression changes improves the interpretability of bulk tissue datasets (Mancarci et al., 2017; Toker et al., 2018).

In this study, we re-analyzed publicly available expression data from the hippocampus of multiple AD mouse models, accounting for between-study differences and cell type-specific effects, with the goal of identifying consistent cross-mouse-model transcriptional changes, especially at early phases.

Materials and Methods

Analyses were implemented in R (Team, 2016) except where noted. The data and scripts used are available at https://github.com/PavlidisLab/AlzMouseModelMeta.

Data pre-processing and quality control

We retrieved gene expression profiling studies of mouse models of AD with the keyword “Alzheimer” from Gene Expression Omnibus (GEO; RRID:SCR_005012) and ArrayExpress (Barrett et al., 2013; Kolesnikov et al., 2015; RRID:SCR_002964). We further filtered the datasets and selected studies that have at least two biological replicate per condition and contained hippocampal samples. Initially 11 independent studies met the selection criteria. All were conducted on microarray platforms. We downloaded expression data and experimental design meta data from GEO. Probeset and gene annotations of the corresponding Affymetrix, Illumina and Agilent platforms were obtained from Gemma (Zoubarev et al., 2012; RRID:SCR_008007). On further examination, dataset GSE36981 was removed because genotype was confounded with sample batches. This left a final group of 10 datasets. Four major types of AD mouse models are among the 10: amyloid transgenic models, tau transgenic models, knock-out (KO) models and anti-NGF AD11 (categorized as “other”). Individual studies varied in their use of male or female animals, with both sexes included in the aggregated data. For overviews of the datasets, see Tables 1, 2.

View this table:
  • View inline
  • View popup
Table 1.

Summary of selected gene expression profiling studies for AD

View this table:
  • View inline
  • View popup
Table 2.

Details of analyzed AD mouse model studies

The quality of the raw expression data of Affymetrix arrays was evaluated as described in (Rogic et al., 2016). All hippocampal samples using Affymetrix arrays passed the quality control procedures. Such raw data quality control procedures were not available for Agilent and Illumina arrays. To standardize data processing, Affymetrix and Agilent arrays were robust multi-array average (RMA) background corrected by affy (Gautier et al., 2004) and limma (Ritchie et al., 2015) R packages (RRID:SCR_012835, RRID:SCR_010943 respectively) followed by quantile normalization and log2 transformation. Illumina arrays were quantile normalized and log2 transformed. Samples with brain tissues other than hippocampus were discarded after normalization. Samples that were outliers (two standard deviations away from the mean sample-to-sample Pearson correlation within a dataset) were removed and the remaining samples were batch-corrected for each dataset by ComBat (Johnson et al., 2007; RRID:SCR_010974) if batch information was available.

The time points when mouse models first develop phenotypes that are similar to the earliest clinical symptoms for diagnosis in AD were used to define early and late phases of AD. The mouse phenotypes of AD mouse models analyzed were based on the behavior data from original publications, or the publications cited in the original paper. Mild cognitive decline is an early diagnostic symptom in AD patients (Webster et al., 2014). Cognitive impairment is often assessed by water maze test in AD mouse models. Therefore, AD mouse samples that did not display impairment in memory and learning measured by water maze tests were categorized as early phase AD samples, while the rest as late phase AD samples. The final dataset constituted data for early time points with controls (116 samples) and late time points with controls (96 samples).

To allow cross-platform comparison, within each dataset, we removed non-specific probes (i.e., probes that mapped to multiple genes), probes that did not map to any genes, and probes that contained missing expression values in one or more samples. When more than one probe mapped to a gene, we retained only the probe with the highest median expression value to represent the mapped gene. Not all the genes are available on all platforms used by the studies; we selected genes that were present in more than at least 2/3 of the platforms as a compromise between maximizing the number of genes in the analysis and the requirement to have multiple measurements to perform a mega-analysis. For each disorder, two integrated datasets were created by combining samples across studies from each disease phase. Within each integrated dataset, gene expression values were quantile normalized to harmonize scales across studies. We then filtered each dataset to remove non-expressed genes. To set the threshold for filtering, we were guided by the expression level of sex-specific genes (Toker et al., 2016). The signal for sex-specific genes in the non-expressing sex (e.g., Y-linked genes in females) can be taken as a rough indicator of background levels. The median expression value of non-expressed sex-specific genes from all samples was 5.2, and thus, we filtered genes with expression value lower than 6 as a more stringent threshold. For the number of genes in each disease phase after gene filtering, see Table 1.

Most of AD mouse models analyzed in this project were transgenic mouse models with transgenes under the control of murine Thy1 tissue-specific regulatory elements. The microarray probesets mapped to these transgenes and the endogenous copy, which artificially increased the measured expression of Thy1. Therefore, Thy1 was removed from the mega-analysis in AD mouse models.

Estimation of cell-type proportion changes

Cell-type proportions of three glial cell types (microglia, astrocytes, oligodendrocytes) and three neuronal cell types (pyramidal cells, dentate granule cells, GABAergic cells) were estimated by marker gene profiles (MGPs) using pre-selected markers specific to the murine hippocampus (Mancarci et al., 2017). Expression of the marker genes were first corrected for between-study variation for each disease phase and then used as input for MGP estimation. For presentation MGPs were normalized to a range between 0 and 1, where the sample with the highest profile was assigned to 1 and the lowest was assigned to 0. Wilcoxon rank-sum test and computation of false discovery rates (FDRs) using the Benjamini–Hochberg procedure (Benjamini and Hochberg, 1995) were applied to test whether the profiles were significantly different between disease mouse models and controls for each cell type.

Fitting linear mixed-effects models (LMMs) and jackknife procedure to rank genes

LMMs allow modeling multiple sources of variation, such as mouse model-specific effects, laboratory effects, and difference between disease models and controls (Demirkale et al., 2010). For each disease phase group, two LMMs were fitted for each gene using the “lmer” function in R package lme4 version 1.1-12 (RRID:SCR_015654), via maximum likelihood estimation [lmer(REML = F); Bates et al., 2015]. The first LMM corrected for between-study variations without correction for the MGPs; the second LMM corrected for both between-study variations and MGPs. The p value for the significance of the fixed effect of disease state (i.e., disease and normal states) was obtained by the “anova” function in R package stats version 3.3.1 (Chambers and Hastie, 1991). Ranking of upregulated or downregulated genes was based on p values in ascending order and the direction of expression changes between mouse models and controls. FDR estimated by the Benjamini–Hochberg procedure was computed for each gene. Significantly DE genes had FDR < 0.05. A jackknife procedure was applied to yield more robust gene rankings.

Functional enrichment analysis

Functional enrichment analysis was performed using the threshold-free precision-recall algorithm in ErmineJ version 3.0.2 (Gillis et al., 2010; RRID:SCR_006450) to determine enrichment of Gene Ontology (GO) terms for the ranked list of genes from the jackknife procedure. We used the multifunctionality-adjusted enrichment rankings provided by ErmineJ to reduce the distorting effect of highly annotated genes (Ballouz et al., 2016).

Results

We reanalyzed a total of 208 gene expression profiles from 10 AD mouse model datasets (Tables 1, 2). With the goal of identifying shared transcriptional alterations among different mouse models in the early and late phases, two separate mega-analyses were conducted: analysis of samples from “early” time points (before most pathologic effects), and a later time point. An important feature of our analysis is the use of MGPs. Used directly, MGPs can be used as an estimate of relative cellular proportions (Mancarci et al., 2017). We also used MGPs as covariates in our linear models to help isolate changes due to regulation as opposed to changes in cellular makeup of the tissue.

Estimation of cell-type proportion changes

The results of MGP estimation were largely consistent with the expected cell-type changes in the early and late phase mouse samples. Relative to controls, mouse model samples in the early phase were predicted to show minimal or small cellular composition changes while late phase samples would show more substantial changes. In agreement, with this hypothesis, during the early phase, MGP analysis indicated no significant changes in neurons and glial cells (Fig. 1A,B). During the late phase, AD mouse models had estimated reduced dentate granule cells and pyramidal cells (Fig. 1A), and increased astrocytes and microglia (Fig. 1B).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

MGPs in AD mouse models. Each row of plots represents MGPs for one cell type as labeled at right. Each point represents one sample. Boxplots representing the interquartile range and the sample means are also shown. The vertical axis is normalized to the range 0 and 1. *** FDR < 0.01. A, Neurons B, Glial cells.

Mega-analysis of gene expression in AD mouse models

Using a mixed effect modeling approach with MGPs as covariates (see Materials and Methods), we identified expression changes associated with early and late phases of the AD mouse models. As expected, early changes were subtle. The mega-analysis of 116 gene expression profiles of mouse models and controls in the early phase revealed small but consistent expression changes (Fig. 2). Only four upregulated and three downregulated genes of the top genes were significant at an FDR of 0.05. A few of the top genes were functionally related based on the enrichment analysis. Three of the top 40 upregulated genes (Sqle, Msmo1, and Nsdhl ranked 2, 3, and 7; FDR < 0.05 except Nsdh1, which met an FDR of 0.08) are involved in cholesterol biosynthesis. C1qa, C1qb, and C1qc, which are involved in the classical complement cascade, were among the top 40 upregulated gene in the early phase (ranks 28, 12, and 29, meeting an FDR of 0.12).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Top 20 upregulated and downregulated genes for late phase. Expression values are corrected for studies and MGPs. Each column is a brain sample, each row is a gene. For display purposes, each row is z score transformed. Gray cells represent missing values. A gap separates the upregulated genes (top) and the downregulated genes (bottom).

In contrast to the subtle changes in the early phase, the expression change signals in the late phase were stronger (Fig. 3). The analysis also detected down-regulation of a known AD risk gene, Trem2, which had also been reported in the original publications of two studies analyzed. Several GO terms were found significantly enriched in the top downregulated genes, including “neuronal cell body” and “regulation of neurogenesis.”

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Top 20 upregulated and downregulated genes for early phase. As in Figure 2, but for the early phase.

We note that using MGPs as covariates have important effects on the results, which we documented by comparing the results of fitting models in which MGPs are not included. First, after MGP correction, as expected the marker genes included in the MGPs themselves are down-ranked, especially in the late phase. The number of significant DE genes (FDR < 0.05) in the late phase was greatly reduced after MGP correction:13.64% of the total genes were DE genes before correction and dropped to 0.22% after correction. We anticipated that the additional parameters included in the models would reduce power, so the fact that fewer genes met the FDR threshold was not in itself surprising. However, the gene rankings changed considerably, especially in the late phase. The Spearman correlation coefficients of before and after correction in the late phase is 0.38, compared to 0.80 in the early phase. Thus, MGP correction affected both the number of DE genes and the gene rankings, which indicated that cellular composition changes (as measured by MGPs) contributed to bulk tissue gene expression changes, especially in the late disease phase.

The gene rankings of early and late AD were quite different, indicating phase-specific effects. The Spearman correlation coefficient for gene ranking of all the genes between early and late phase was 0.24 before MGPs correction and 0.25 after. There were very few genes in common among the top 100 dysregulated genes across the phases. Specifically, only the top genes involved in cholesterol biosynthesis (Sqle, Msmo1, Nsdhl) were dysregulated in concordance in both phases (FDR < 0.05) for both before and after MGPs correction.

Discussion

We identified consistent transcriptomic alterations across mouse models of AD. By categorizing samples into early and late disease phases, gene expression changes specific to the disease phase were revealed. This revealed subtle but biologically interpretable changes shared across mouse models in the early phase, which may reflect reveal early disease mechanisms. Changes in the late phase were stronger, and as expected were more associated with changes in cellular populations. The top-ranked genes in the early phase were not always affected in the late phase, and vice versa, indicating phase-specific expression changes.

Consistent transcriptomic alterations were identified across different mouse models

Candidate genes identified in the original published reports were often inconsistent, though some genes are reported in more than one study. This may be due to differences in analysis methods and thresholds used, but doing a mega-analysis has the advantage of facilitating detection of weaker signals. Thus our mega-analysis was able to capture consistent signals that were not reported in the original source studies. For example, in the late AD phase, two out of seven studies reported the AD risk gene, Trem2, in their hit lists. However, we observed consistent up-regulation of Trem2 in other studies (Fig. 4A). Similarly, we identified genes that had not been reported in any of the original studies, included a top upregulated gene, Msmo1, in the late AD phase (Fig. 4).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Examples of gene expression patterns identified. Each point represents one sample. Control samples (WT) are marked red and AD mouse model samples are blue. “After MGP” reflects analyses in which MGPs were included as covariates in the model; this correction rescales the expression levels to arbitrary units. Studies (datasets) that are marked with an asterisk have reported the gene as top hit in the original publications. A, Expression of Trem2 in the late AD phase before (top row) and after (bottom row) MGP correction. B, Expressions of Msmo1 in the late AD phase before (top) and after MGP correction (bottom).

Applying cell population proportion correction revealed transcriptional changes of cell type-specific regulatory events

We used MGPs to adjust the expression profiles for the effects of cellular proportion changes (Mancarci et al., 2017). We note that MGP analysis should interpreted as reflecting cellular proportions with caution, as it is an indirect measure. A more conservative interpretation of an MGP change is a cell type-specific change in expression. However, given the known neurodegeneration in these mouse models, changes in cellular proportions are the likely cause in MGP shifts in our study (Serrano-Pozo et al., 2011). Indeed, the cell-type proportion changes as estimated by MGPs in AD mouse models were consistent with previous mouse and human studies (Schmitz et al., 2004; Serrano-Pozo et al., 2011; Hokama et al., 2014), and indicated neuronal loss and gliosis in the hippocampus in the late phase (Fig. 1A,B). The lack of significant proportion changes of neurons and glial cells before the occurrence of cognitive impairment in AD mouse models (Fig. 1A,B) also agreed with a report that cognitive impairment, neuronal loss and gliosis occur concurrently in the hippocampus of AD mouse model J20 (Wright et al., 2013). Therefore, we interpret the MGP shifts as cellular proportion changes while keeping the indirect nature of the measure in mind.

The changes of gene rankings, and decreased DE signal after MGPs correction, implies that the gene expression changes, especially in the late disease phase, were substantially driven by changes in cell-type proportions. Adjusting MGPs can reveal cell type-specific transcriptional changes. Some of the markers remained top-ranked of the same direction of regulation even after correction, such as top-up-ranked microglia markers Cd68 and Tyrobp, and astrocyte marker Slc14a1 in the late AD phase remained top-up-ranked after correction. Because the dysregulation of these marker genes cannot be fully explained by the MGP changes, they may indicate changes at transcriptional regulation level within microglia that contribute to disease pathophysiology. However, because the changes in composition are confounded with the experimental condition, it is also possible that some residual effects of cellular proportion remain even after including MGPs as covariates. Ultimately resolution of cell type-specific changes in gene expression in these models, especially at late stages, will require cell type-specific transcriptomic studies.

Shared gene expression changes revealed upregulation of genes in cholesterol biosynthesis and classical complement cascade in AD mouse models

Our mega-analysis of the early phase revealed up-regulation of genes in cholesterol biosynthesis and classical complement cascade. The changes were subtle in the early phase compared to the late phase. However, the results had some biological coherence as suggested by the GO enrichment analysis. Thus, the enrichment of genes involved in cholesterol biosynthesis in both early and late phases, even after MGP correction, suggested chronic dysregulation in cholesterol biosynthesis. Up-regulation of genes in cholesterol biosynthesis has been also observed in an AD amyloid transgenic mouse model, APP23 (Tseveleki et al., 2010). Many lines of evidence have linked cholesterol to AD, and to Aβ production in particular (Puglielli et al., 2003; Wood et al., 2014), though the interpretation and implications are still unclear. In mouse and cell culture studies, decreased brain cholesterol levels can reduce Aβ abundance (Di Paolo and Kim, 2011; Wood et al., 2014). Cleavage of APP by β-secretase and γ-secretase (the amyloidogenic pathway) mainly occurs in lipid rafts of the plasma membrane, whereas α-secretase of the non-amyloidogenic pathway tend to localize at the non-lipid-raft sites (Kim et al., 2016). Lipid rafts have high concentration of cholesterol (Kim et al., 2016). Increased level of cholesterol enhances localization of APP, β-secretase and γ-secretase to the lipid rafts, and subsequently promotes Aβ production (Marquer et al., 2011; Kim et al., 2016). However, the mechanisms that initiate the increased gene expression of cholesterol biosynthesis genes in mouse models are not known. The majority of samples analyzed here were amyloid transgenic models, which express human transgenes (APP, PSEN1, PSEN2) with known AD associated mutations that promote APP processing through the amyloidogenic pathway. Expression of these transgenes could play a role in promoting the expression of cholesterol biosynthesis genes, which lead to increased cholesterol level that can accelerate Aβ production. The most established AD risk factor is apolipoprotein E (ApoE, encoded by APOE), which mediates cholesterol metabolism in the brain and is found in Aβ plaques and neurofibrillary tangles (Yang and Song, 2013; Ries and Sastre, 2016). It has been shown that the high-risk isoform APOE4 has lower efficiency in transporting cholesterol from astrocytes to neurons compared to the neutral isoform APOE3, and lead to synaptic dysfunctions (Bu, 2009). Our results suggest cholesterol may play a role in AD-like process initiation and progression in mouse models.

Complement pathway genes, C1qa, C1qb, and C1qc, were among the top-ranked upregulated genes in the early phase, suggested up-regulation of genes in the C1q pathway in early disease progression. The up-regulation of C1q pathway, which initiates classical complement cascade, has been linked to early synaptic loss before Aβ accumulation, and is a response to injury in mouse model, and could be neuroprotective against misfolded proteins (Benoit et al., 2013; Hong et al., 2016). C1q pathway genes were still differentially upregulated in the late phase, however, they were no longer among the top-ranked genes after MGP correction. This suggests that cellular composition changes are more responsible for RNA levels changes of C1q genes in the late phase.

In the late phase, the GO term regulation of neurogenesis was enriched in the downregulated genes. The result was consistent with previous studies (Noh et al., 2014; Matarin et al., 2015), reflecting neuronal dysfunction in the AD mouse brain, though as noted above this could also potentially reflect imperfect correction for proportion effects. Interestingly, GO terms related to inflammation and the immune system were not enriched in upregulated genes. These results are contrary to other studies, which report up-regulations of genes related to inflammation and immune system (Matarin et al., 2015; Saura et al., 2015). This seeming discrepancy can be explained as a combination of effects of cell-type proportion considerations and gene function annotations. Without MGPs correction, our results did indicate up-regulation of genes in the immune response pathways. However, closer inspection revealed that some of the microglia markers are annotated with GO terms “immune response” and “inflammatory response.” This further complicates separating cell-type proportion changes regulatory effects. On the other hand, the dysregulation of a few astrocytes and microglia markers could not be fully explained by the cell-type population changes and remained top-ranked dysregulated genes after MGPs correction. For example, Cd68, a marker for microglial activation and is correlated with Aβ42 load (Zotova et al., 2011), was upregulated in the late phase. Another microglia marker, TYRO protein tyrosine kinase-binding protein gene (Tyrobp) and its receptor Trem2 were among the top upregulated genes. A previous study reported that TYROBP binds to TREM2 and promotes microglial activation (Kobayashi et al., 2016). These results may provide some evidence of microglial activation at the cell type-specific level in AD mouse models.

In conclusion, our mega-analysis of gene expression in mouse models revealed consistent and disease phase-specific transcriptional changes and cell type-specific regulatory events, despite the considerable heterogeneity of the mouse models of AD. The identification of shared gene expression changes in the early phase increases our understanding of disease initiation and progression. Prioritized top-ranked genes in the early phase can be candidate genes to study mechanisms in disease initiation and worthy of further investigation.

Acknowledgments

Acknowledgements: We thank Weihong Song, Gabriela Cohen Freue, and the Pavlidis lab for discussion. We also thank all the authors who contributed their data to the public repositories, which enabled this study.

Footnotes

  • The authors declare no competing financial interests.

  • This work was supported by the Canadian Institutes of Health Research Joint Program on Neurodegenerative Disease Research, the National Institutes of Health Grant MH111099, the Natural Sciences and Engineering Research Council of Canada, and the University of British Columbia Graduate Program in Bioinformatics.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Ballouz S, Pavlidis P, Gillis J (2016) Using predictive specificity to determine when gene set analysis is biologically meaningful. Nucleic Acids Res 45:e20.
    OpenUrl
  2. ↵
    Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, Robertson CL, Serova N, Davis S, Soboleva A (2013) NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res 41:D991–D995. doi:10.1093/nar/gks1193 pmid:23193258
    OpenUrlCrossRefPubMed
  3. ↵
    Bates D, Mächler M, Bolker B, Walker S (2015) Fitting linear mixed-effects models using lme4. J Stat Softw 67:1–48. doi:10.18637/jss.v067.i01
    OpenUrlCrossRef
  4. ↵
    Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol 57:289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
    OpenUrlCrossRef
  5. ↵
    Benoit ME, Hernandez MX, Dinh ML, Benavente F, Vasquez O, Tenner AJ (2013) C1q-induced LRP1B and GPR6 proteins expressed early in Alzheimer disease mouse models, are essential for the C1q-mediated protection against amyloid-β neurotoxicity. J Biol Chem 288:654–665. doi:10.1074/jbc.M112.400168 pmid:23150673
    OpenUrlAbstract/FREE Full Text
  6. ↵
    Bu G (2009) Apolipoprotein E and its receptors in Alzheimer’s disease: pathways, pathogenesis and therapy. Nat Rev Neurosci 10:333–344. doi:10.1038/nrn2620 pmid:19339974
    OpenUrlCrossRefPubMed
  7. ↵
    Chambers JM, Hastie TJ (1991) Statistical models in S. Boca Raton, FL: Chapman and Hall/CRC.
  8. ↵
    Chikina M, Zaslavsky E, Sealfon SC (2015) CellCODE: a robust latent variable approach to differential expression analysis for heterogeneous cell populations. Bioinformatics 31:1584–1591. doi:10.1093/bioinformatics/btv015
    OpenUrlCrossRefPubMed
  9. Cummings DM, Liu W, Portelius E, Bayram S, Yasvoina M, Ho S-H, Smits H, Ali SS, Steinberg R, Pegasiou C-M, James OT, Matarin M, Richardson JC, Zetterberg H, Blennow K, Hardy JA, Salih DA, Edwards FA (2015) First effects of rising amyloid-β in transgenic mouse brain: synaptic transmission and gene expression. Brain J Neurol 138:1992–2004. doi:10.1093/brain/awv127 pmid:25981962
    OpenUrlCrossRefPubMed
  10. ↵
    Demirkale CY, Nettleton D, Maiti T (2010) Linear mixed model selection for false discovery rate control in microarray data analysis. Biometrics 66:621–629. doi:10.1111/j.1541-0420.2009.01286.x pmid:19522873
    OpenUrlCrossRefPubMed
  11. ↵
    Di Paolo G, Kim T-W (2011) Linking lipids to Alzheimer’s disease: cholesterol and beyond. Nat Rev Neurosci 12:284–296. doi:10.1038/nrn3012 pmid:21448224
    OpenUrlCrossRefPubMed
  12. ↵
    D’Onofrio M, Arisi I, Brandi R, Di Mambro A, Felsani A, Capsoni S, Cattaneo A (2011) Early inflammation and immune response mRNAs in the brain of AD11 anti-NGF mice. Neurobiol Aging 32:1007–1022. doi:10.1016/j.neurobiolaging.2009.05.023 pmid:19604602
    OpenUrlCrossRefPubMed
  13. ↵
    Gaujoux R, Seoighe C (2013) CellMix: a comprehensive toolbox for gene expression deconvolution. Bioinformatics 29:2211–2212. doi:10.1093/bioinformatics/btt351 pmid:23825367
    OpenUrlCrossRefPubMed
  14. ↵
    Gautier L, Cope L, Bolstad BM, Irizarry RA (2004) affy–analysis of Affymetrix GeneChip data at the probe level. Bioinforma Oxf Engl 20:307–315. doi:10.1093/bioinformatics/btg405 pmid:14960456
    OpenUrlCrossRefPubMed
  15. ↵
    Gillis J, Mistry M, Pavlidis P (2010) Gene function analysis in complex data sets using ErmineJ. Nat Protoc 5:1148–1159. doi:10.1038/nprot.2010.78 pmid:20539290
    OpenUrlCrossRefPubMed
  16. ↵
    Hokama M, Oka S, Leon J, Ninomiya T, Honda H, Sasaki K, Iwaki T, Ohara T, Sasaki T, LaFerla FM, Kiyohara Y, Nakabeppu Y (2014) Altered expression of diabetes-related genes in Alzheimer’s disease brains: the Hisayama study. Cereb Cortex 24:2476–2488. doi:10.1093/cercor/bht101 pmid:23595620
    OpenUrlCrossRefPubMed
  17. ↵
    Hong H, Hong Q, Liu J, Tong W, Shi L (2013) Estimating relative noise to signal in DNA microarray data. Int J Bioinforma Res Appl 9:433–448. doi:10.1504/IJBRA.2013.056085 pmid:24001721
    OpenUrlCrossRefPubMed
  18. ↵
    Hong S, Beja-Glasser VF, Nfonoyim BM, Frouin A, Li S, Ramakrishnan S, Merry KM, Shi Q, Rosenthal A, Barres BA, Lemere CA, Selkoe DJ, Stevens B (2016) Complement and microglia mediate early synapse loss in Alzheimer mouse models. Science 352:712–716. doi:10.1126/science.aad8373 pmid:27033548
    OpenUrlAbstract/FREE Full Text
  19. ↵
    Johnson WE, Li C, Rabinovic A (2007) Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8:118–127. doi:10.1093/biostatistics/kxj037 pmid:16632515
    OpenUrlCrossRefPubMed
  20. ↵
    Kim Y, Kim C, Jang HY, Mook-Jung I (2016) Inhibition of cholesterol biosynthesis reduces γ-secretase activity and amyloid-β generation. J Alzheimers Dis 51:1057–1068. doi:10.3233/JAD-150982 pmid:26923021
    OpenUrlCrossRefPubMed
  21. Kleiman RJ, Lanz TA, Finley JE, Bove SE, Majchrzak MJ, Becker SL, Carvajal-Gonzales S, Kuhn AM, Wood KM, Mariga A, Nelson FR, Verhoest PR, Seymour PA, Stephenson DT (2010) Dendritic spine density deficits in the hippocampal CA1 region of young Tg2576 mice are ameliorated with the PDE9A inhibitor PF-04447943. Alzheimers Dement 6:S563–S564. doi:10.1016/j.jalz.2010.05.1922
    OpenUrlCrossRef
  22. ↵
    Kobayashi M, Konishi H, Sayo A, Takai T, Kiyama H (2016) TREM2/DAP12 signal elicits proinflammatory response in microglia and exacerbates neuropathic pain. J Neurosci 36:11138–11150. doi:10.1523/JNEUROSCI.1238-16.2016
    OpenUrlAbstract/FREE Full Text
  23. ↵
    Kolesnikov N, Hastings E, Keays M, Melnichuk O, Tang YA, Williams E, Dylag M, Kurbatova N, Brandizi M, Burdett T, Megy K, Pilicheva E, Rustici G, Tikhonov A, Parkinson H, Petryszak R, Sarkans U, Brazma A (2015) ArrayExpress update–simplifying data submissions. Nucleic Acids Res 43:D1113–D1116. doi:10.1093/nar/gku1057 pmid:25361974
    OpenUrlCrossRefPubMed
  24. Li H, Wang B, Wang Z, Guo Q, Tabuchi K, Hammer RE, Südhof TC, Zheng H (2010) Soluble amyloid precursor protein (APP) regulates transthyretin and Klotho gene expression without rescuing the essential function of APP. Proc Natl Acad Sci USA 107:17362–17367. doi:10.1073/pnas.1012568107 pmid:20855613
    OpenUrlAbstract/FREE Full Text
  25. ↵
    Mancarci BO, Toker L, Tripathy SJ, Li B, Rocco B, Sibille E, Pavlidis P (2017) Cross-laboratory analysis of brain cell type transcriptomes with applications to interpretation of bulk tissue data. eNeuro 4:ENEURO.0212-17.2017. doi:10.1523/ENEURO.0212-17.2017 pmid:29204516
    OpenUrlAbstract/FREE Full Text
  26. ↵
    Marquer C, Devauges V, Cossec JC, Liot G, Lécart S, Saudou F, Duyckaerts C, Lévêque-Fort S, Potier MC (2011) Local cholesterol increase triggers amyloid precursor protein-Bace1 clustering in lipid rafts and rapid endocytosis. FASEB J 25:1295–1305. doi:10.1096/fj.10-168633 pmid:21257714
    OpenUrlCrossRefPubMed
  27. ↵
    Masters CL, Bateman R, Blennow K, Rowe CC, Sperling RA, Cummings JL (2015) Alzheimer’s disease. Nat Rev Dis Primer 1:15056.
    OpenUrl
  28. ↵
    Matarin M, Salih DA, Yasvoina M, Cummings DM, Guelfi S, Liu W, Nahaboo Solim MA, Moens TG, Paublete RM, Ali SS, Perona M, Desai R, Smith KJ, Latcham J, Fulleylove M, Richardson JC, Hardy J, Edwards FA (2015) A genome-wide gene-expression analysis and database in transgenic mice during development of amyloid or tau pathology. Cell Rep 10:633–644. doi:10.1016/j.celrep.2014.12.041 pmid:25620700
    OpenUrlCrossRefPubMed
  29. Nagahara AH, Merrill DA, Coppola G, Tsukada S, Schroeder BE, Shaked GM, Wang L, Blesch A, Kim A, Conner JM, Rockenstein E, Chao MV, Koo EH, Geschwind D, Masliah E, Chiba AA, Tuszynski MH (2009) Neuroprotective effects of brain-derived neurotrophic factor in rodent and primate models of Alzheimer’s disease. Nat Med 15:331–337. doi:10.1038/nm.1912 pmid:19198615
    OpenUrlCrossRefPubMed
  30. ↵
    Noh H, Park C, Park S, Lee YS, Cho SY, Seo H (2014) Prediction of miRNA-mRNA associations in Alzheimer’s disease mice using network topology. BMC Genomics 15:644. doi:10.1186/1471-2164-15-644 pmid:25086961
    OpenUrlCrossRefPubMed
  31. Paesler K, Xie K, Hettich MM, Siwek ME, Ryan DP, Schröder S, Papazoglou A, Broich K, Müller R, Trog A, Garthe A, Kempermann G, Weiergräber M, Ehninger D (2015) Limited effects of an eIF2αS51A allele on neurological impairments in the 5xFAD mouse model of Alzheimer’s disease. Neural Plast 2015:825157. doi:10.1155/2015/825157 pmid:25883808
    OpenUrlCrossRefPubMed
  32. ↵
    Pennanen C, Kivipelto M, Tuomainen S, Hartikainen P, Hänninen T, Laakso MP, Hallikainen M, Vanhanen M, Nissinen A, Helkala E-L, Vainio P, Vanninen R, Partanen K, Soininen H (2004) Hippocampus and entorhinal cortex in mild cognitive impairment and early AD. Neurobiol Aging 25:303–310. doi:10.1016/S0197-4580(03)00084-8 pmid:15123335
    OpenUrlCrossRefPubMed
  33. Pereson S, Wils H, Kleinberger G, McGowan E, Vandewoestyne M, Van Broeck B, Joris G, Cuijt I, Deforce D, Hutton M, Van Broeckhoven C, Kumar-Singh S (2009) Progranulin expression correlates with dense-core amyloid plaque burden in Alzheimer disease mouse models. J Pathol 219:173–181. doi:10.1002/path.2580 pmid:19557827
    OpenUrlCrossRefPubMed
  34. Polito VA, Li H, Martini-Stoica H, Wang B, Yang L, Xu Y, Swartzlander DB, Palmieri M, di Ronza A, Lee VY, Sardiello M, Ballabio A, Zheng H (2014) Selective clearance of aberrant tau proteins and rescue of neurotoxicity by transcription factor EB. EMBO Mol Med 6:1142–1160. doi:10.15252/emmm.201303671 pmid:25069841
    OpenUrlAbstract/FREE Full Text
  35. ↵
    Puglielli L, Tanzi RE, Kovacs DM (2003) Alzheimer’s disease: the cholesterol connection. Nat Neurosci 6:345–351. doi:10.1038/nn0403-345 pmid:12658281
    OpenUrlCrossRefPubMed
  36. ↵
    Ries M, Sastre M (2016) Mechanisms of Aβ clearance and degradation by glial cells. Front Aging Neurosci 8:160.
    OpenUrlCrossRefPubMed
  37. ↵
    Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK (2015) limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 43:e47.
    OpenUrlCrossRefPubMed
  38. ↵
    Rogic S, Wong A, Pavlidis P (2016) Meta-Analysis of Gene Expression Patterns in Animal Models of Prenatal Alcohol Exposure Suggests Role for Protein Synthesis Inhibition and Chromatin Remodeling. Alcohol Clin Exp Res 40:717–727.
    OpenUrlCrossRef
  39. ↵
    Saura CA, Parra-Damas A, Enriquez-Barreto L (2015) Gene expression parallels synaptic excitability and plasticity changes in Alzheimer’s disease. Front Cell Neurosci 9:318.
    OpenUrl
  40. ↵
    Saxena S, Caroni P (2011) Selective neuronal vulnerability in neurodegenerative diseases: from stressor thresholds to degeneration. Neuron 71:35–48. doi:10.1016/j.neuron.2011.06.031 pmid:21745636
    OpenUrlCrossRefPubMed
  41. ↵
    Schmitz C, Rutten BPF, Pielen A, Schäfer S, Wirths O, Tremp G, Czech C, Blanchard V, Multhaup G, Rezaie P, Korr H, Steinbusch HWM, Pradier L, Bayer TA (2004) Hippocampal neuron loss exceeds amyloid plaque load in a transgenic mouse model of Alzheimer’s disease. Am J Pathol 164:1495–1502. doi:10.1016/S0002-9440(10)63235-X pmid:15039236
    OpenUrlCrossRefPubMed
  42. ↵
    Serrano-Pozo A, Frosch MP, Masliah E, Hyman BT (2011) Neuropathological alterations in Alzheimer disease. Cold Spring Harb Perspect Med 1:a006189. doi:10.1101/cshperspect.a006189 pmid:22229116
    OpenUrlAbstract/FREE Full Text
  43. ↵
    Shen-Orr SS, Gaujoux R (2013) Computational deconvolution: extracting cell type-specific information from heterogeneous samples. Curr Opin Immunol 25:571–578. doi:10.1016/j.coi.2013.09.015 pmid:24148234
    OpenUrlCrossRefPubMed
  44. ↵
    Srinivasan K, Friedman BA, Larson JL, Lauffer BE, Goldstein LD, Appling LL, Borneo J, Poon C, Ho T, Cai F, Steiner P, van der Brug MP, Modrusan Z, Kaminker JS, Hansen DV (2016) Untangling the brain’s neuroinflammatory and neurodegenerative transcriptional responses. Nat Commun 7:11295. doi:10.1038/ncomms11295 pmid:27097852
    OpenUrlCrossRefPubMed
  45. Stein TD, Anders NJ, DeCarli C, Chan SL, Mattson MP, Johnson JA (2004) Neutralization of transthyretin reverses the neuroprotective effects of secreted amyloid precursor protein (APP) in APPSW mice resulting in tau phosphorylation and loss of hippocampal neurons: support for the amyloid hypothesis. J Neurosci 24:7707–7717. doi:10.1523/JNEUROSCI.2211-04.2004
    OpenUrlAbstract/FREE Full Text
  46. ↵
    Team RC (2016) R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.
  47. ↵
    Toker L, Feng M, Pavlidis P (2016) Whose sample is it anyway? Widespread misannotation of samples in transcriptomics studies. F1000Research 5:2103. doi:10.12688/f1000research.9471.1
    OpenUrlCrossRef
  48. ↵
    Toker L, Mancarci BO, Tripathy S, Pavlidis P (2018) Transcriptomic evidence for alterations in astrocytes and parvalbumin interneurons in subjects with bipolar disorder and schizophrenia. Biol Psychiatry 84:787–796. doi:10.1016/j.biopsych.2018.07.010 pmid:30177255
    OpenUrlCrossRefPubMed
  49. ↵
    Tseveleki V, Rubio R, Vamvakas S-S, White J, Taoufik E, Petit E, Quackenbush J, Probert L (2010) Comparative gene expression analysis in mouse models for multiple sclerosis, Alzheimer’s disease and stroke for identifying commonly regulated and disease-specific gene changes. Genomics 96:82–91. doi:10.1016/j.ygeno.2010.04.004 pmid:20435134
    OpenUrlCrossRefPubMed
  50. ↵
    Webster SJ, Bachstetter AD, Nelson PT, Schmitt FA, Van Eldik LJ (2014) Using mice to model Alzheimer’s dementia: an overview of the clinical disease and the preclinical behavioral changes in 10 mouse models. Front Genet 5:88.
    OpenUrlCrossRefPubMed
  51. ↵
    Wirz KTS, Keitel S, Swaab DF, Verhaagen J, Bossers K (2014) Early molecular changes in Alzheimer disease: can we catch the disease in its presymptomatic phase? J Alzheimers Dis 38:719–740. doi:10.3233/JAD-130920 pmid:24072070
    OpenUrlCrossRefPubMed
  52. ↵
    Wood WG, Li L, Müller WE, Eckert GP (2014) Cholesterol as a causative factor in Alzheimer’s disease: a debatable hypothesis. J Neurochem 129:559–572. doi:10.1111/jnc.12637 pmid:24329875
    OpenUrlCrossRefPubMed
  53. ↵
    Wright AL, Zinn R, Hohensinn B, Konen LM, Beynon SB, Tan RP, Clark IA, Abdipranoto A, Vissel B (2013) Neuroinflammation and neuronal loss precede Aβ plaque deposition in the hAPP-J20 mouse model of Alzheimer’s disease. PLoS One 8:e59586. doi:10.1371/journal.pone.0059586 pmid:23560052
    OpenUrlCrossRefPubMed
  54. ↵
    Yang Y, Song W (2013) Molecular links between Alzheimer’s disease and diabetes mellitus. Neuroscience 250:140–150. doi:10.1016/j.neuroscience.2013.07.009 pmid:23867771
    OpenUrlCrossRefPubMed
  55. ↵
    Zotova E, Holmes C, Johnston D, Neal JW, Nicoll J. a R, Boche D (2011) Microglial alterations in human Alzheimer’s disease following Aβ42 immunization. Neuropathol Appl Neurobiol 37:513–524. doi:10.1111/j.1365-2990.2010.01156.x pmid:21166690
    OpenUrlCrossRefPubMed
  56. ↵
    Zoubarev A, Hamer KM, Keshav KD, McCarthy EL, Santos JRC, Van Rossum T, McDonald C, Hall A, Wan X, Lim R, Gillis J, Pavlidis P (2012) Gemma: a resource for the reuse, sharing and meta-analysis of expression profiling data. Bioinforma Oxf Engl 28:2272–2273. doi:10.1093/bioinformatics/bts430 pmid:22782548
    OpenUrlCrossRefPubMed

Synthesis

Reviewing Editor: Menno Witter, Norwegian University of Science and Technology

Decisions are customarily a result of the Reviewing Editor and the peer reviewers coming together and discussing their recommendations until a consensus is reached. When revisions are invited, a fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision will be listed below. The following reviewer(s) agreed to reveal their identity: Peng Lei, Lene Olsen.

The manuscript entitled “Meta-analysis of gene expression in mouse models of Alzheimer's disease” has provided an interesting attempt to use bioinformatics for the understanding of pathogenesis of Alzheimer's disease. Such a meta-analysis of transcription data from the hippocampus of multiple animal models of Alzheimer's disease has not been done before. The authors report that consistent transcriptomic alterations across mouse early and late models of AD, identified genes that had not been reported yet and pointed the importance of MGPs.

The use of multiple types of animal models allows identification of common pathways involved resulting in similar pathology. The study differentiates sampling before and after onset of cognitive symptoms in the animals, which allows for identification of stage-specific molecular processes. The authors attempt to correct for the cellular composition effect, which is a well-known problem with bulk transcriptomics data on the disease due to loss of neurons and gliosis during disease progression. Indeed, they found that cell type composition dominates the transcription data at the later stages.

The main advance this study provides for the field is that a such systematic meta-study, if done correctly, might increase the power to detect differences between transgenic and wildtype models, and improve the validity of the findings of the studies included. In addition, the consistently differentially expressed genes they identified can be used as a resource for others in future studies on Alzheimer's disease.

We thus recognize the importance of meta-analysis on varies AD models but there are considerable doubts on how these models have been grouped and some other methodological issues raised by the reviewers.

1. The process of meta-analysis needs to be standardized. Although the authors gave the descriptions of data pre-processing and quality control, using flow charts to show the inclusion and exclusion process of original research for meta-analysis is necessary. They also need to state the methods for meta-analysis. Are they using combining P values and combining effect sizes? If properly described, the applied statistical data might be appropriate for the study and ensure robust findings across the data sets. However, due to journal policy, the code was not provided, so the reviewers could not determine that the analyses were conducted in an appropriate manner. Therefore, a short description in the methods section might be appropriate.

2. The current manuscript has presented the differences between early and late AD mouse model, which is interesting. However, it is important to provide a clear definition of early and late stage of AD.

3. The authors grouped varies models based on the presentation of pathology. However, what is the rational to combine amyloid pathology and tau pathology for any analysis? To perform such meta-analysis, the authors made the assumption that the pathologies were caused by the same process, which is not supported by the experimental evidences in the literature.

4. Many of ‘AD’ models used in this study is questionable. For example, AD11 anti-NGF mice is rarely considered as AD models. Introducing such model into the analysis could easily introduce more inflammatory factors that may not be of importance. Similarly, APP KO can hardly be considered as early stage of AD, nor does APLP2 KO.

5. Many of the studies included in the analysis only contain 2 animals. Is it statistically correct to include these data?

6. Is the categorization of late and early time points performed according to the phenotype of the particular animal model, or was a general time point set for all?

7. The paper only indicates QC for Affymetrix data. What about the data from other platforms?

8. What did the authors do with multiple probes covering the same gene?

9. It would be clearer for the reader if the authors specified at least once that they only used microarray and not sequencing data.

10. The authors state that only four genes were significantly upregulated at early time points after adjusting for multiple testing. However, they list three genes involved in cholesterol biosynthesis, and three in the complement cascade as being among the top upregulated results. Did they make some sort of cut off apart from the FDR < 0.05?

11. Some of the data sets contain testing of conditions other than just transgene vs. control mice. While the number of samples included in the table 2 indicates removal of such samples, it is not stated in the method section that such samples were removed.

12. A few mistakes in the text that warrant correction:

• There are some mistakes in the references, with the year being listed twice (e.g. on line 45)

• Mistake line 159: The word “Figur” should be deleted.

• In Table 1, the letter ‘c’ should be deleted from the header in the last column.

• Lines 196-199 states that Trem2 is downregulated, but that Msmo1 is upregulated. Figure 4 shows that these two genes are both upregulated in the Alzheimer's model mice compared to controls.

Back to top

In this issue

eneuro: 6 (6)
eNeuro
Vol. 6, Issue 6
November/December 2019
  • Table of Contents
  • Index by author
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Mega-Analysis of Gene Expression in Mouse Models of Alzheimer’s Disease
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Citation Tools
Mega-Analysis of Gene Expression in Mouse Models of Alzheimer’s Disease
Beryl Zhuang, B. Ogan Mancarci, Lilah Toker, Paul Pavlidis
eNeuro 25 November 2019, 6 (6) ENEURO.0226-19.2019; DOI: 10.1523/ENEURO.0226-19.2019

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Mega-Analysis of Gene Expression in Mouse Models of Alzheimer’s Disease
Beryl Zhuang, B. Ogan Mancarci, Lilah Toker, Paul Pavlidis
eNeuro 25 November 2019, 6 (6) ENEURO.0226-19.2019; DOI: 10.1523/ENEURO.0226-19.2019
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • Alzheimer’s
  • bioinformatics
  • genomics
  • metaanalysis
  • transcriptome

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

New Research

  • Postnatal fluoxetine treatment alters perineuronal net formation and maintenance in the hippocampus
  • The Contribution of Environmental Enrichment to Phenotypic Variation in Mice and Rats
  • Unique effects of social defeat stress in adolescent male mice on the Netrin-1/DCC pathway, prefrontal cortex dopamine and cognition (Social stress in adolescent vs. adult male mice)
Show more New Research

Disorders of the Nervous System

  • Postnatal fluoxetine treatment alters perineuronal net formation and maintenance in the hippocampus
  • The Contribution of Environmental Enrichment to Phenotypic Variation in Mice and Rats
  • Unique effects of social defeat stress in adolescent male mice on the Netrin-1/DCC pathway, prefrontal cortex dopamine and cognition (Social stress in adolescent vs. adult male mice)
Show more Disorders of the Nervous System

Subjects

  • Disorders of the Nervous System
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2021 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.