# A Microglial Signature Directing Human Aging and Neurodegeneration-Related Gene Networks

Shradha Mukherjee, Christine Klaus, Mihaela Pricop-Jeckstadt, Jeremy A Miller, Felix L Struebing

► To cite this version:

Shradha Mukherjee, Christine Klaus, Mihaela Pricop-Jeckstadt, Jeremy A Miller, Felix L Struebing. A Microglial Signature Directing Human Aging and Neurodegeneration-Related Gene Networks. Frontiers in Neuroscience, 2019, 13, 10.3389/fnins.2019.00002. [https://doi.org/10.3389/fnins.2019.00002](https://doi.org/10.3389/fnins.2019.00002). [hal-04088107](https://HAL.inist.fr/hal-04088107)

HAL Id: hal-04088107

https://hal.science/hal-04088107

Submitted on 3 May 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Microglial Signature Directing Human Aging and Neurodegeneration-Related Gene Networks

OPEN ACCESS

**Edited by:**
Rupert W. Overall,
German Center for Neurodegenerative
Diseases (DZNE), Germany

**Reviewed by:**
Björn Spittau,
Universitätsmedizin Rostock,
Germany
Khyobeni Mozhui,
University of Tennessee Health
Science Center (UTHSC),
United States
Michael Oldham,
University of California,
San Francisco, United States

*Correspondence:
Shradha Mukherjee
smukher2@gmail.com
Felix L. Struebing
felix.struebing@
med.uni-muenchen.de

† These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Neurogenomics,
a section of the journal
Frontiers in Neuroscience

**Received:** 31 August 2018
**Accepted:** 03 January 2019
**Published:** 24 January 2019

**Citation:** Mukherjee S, Klaus C, Pricop-Jeckstadt M, Miller JA and Struebing FL (2019) A Microglial Signature Directing Human Aging and Neurodegeneration-Related Gene Networks. Front. Neurosci. 13:2. doi: 10.3389/fnins.2019.00002

Shradha Mukherjee¹,²,³*†, Christine Klaus⁴†, Mihaela Pricop-Jeckstadt⁵, Jeremy A. Miller⁶ and Felix L. Struebing⁷,⁸,⁹*†

¹ Health Informatics Advanced Science Masters Program, Arizona State University, Tempe, AZ, United States, ² Department of Neurology, University of California, Los Angeles, Los Angeles, CA, United States, ³ Department of Bioinformatics, University of California, Los Angeles, Los Angeles, CA, United States, ⁴ Neural Regeneration Group, Institute of Reconstructive Neurobiology, University of Bonn, Bonn, Germany, ⁵ Institute for Medical Informatics and Biometry, Faculty of Medicine "Carl Gustav Carus", TU Dresden, Dresden, Germany, ⁶ Allen Institute for Brain Science, Seattle, WA, United States, ⁷ Department of Translational Brain Research, German Center for Neurodegenerative Diseases, Munich, Germany, ⁸ Center for Neuropathology and Prion Research, Ludwig Maximilian University of Munich, Munich, Germany, ⁹ Department of Ophthalmology, Emory University, Atlanta, GA, United States

Aging is regarded as a major risk factor for neurodegenerative diseases. Thus, a better understanding of the similarities between the aging process and neurodegenerative diseases at the cellular and molecular level may reveal better understanding of this detrimental relationship. In the present study, we mined publicly available gene expression datasets from healthy individuals and patients affected by neurodegenerative diseases (Alzheimer's disease, Parkinson's disease, and Huntington's disease) across a broad age spectrum and compared those with mouse aging and mouse cell-type specific gene expression profiles. We performed weighted gene co-expression network analysis (WGCNA) and found a gene network strongly related with both aging and neurodegenerative diseases. This network was significantly enriched with a microglial signature as imputed from cell type-specific sequencing data. Since mouse models are extensively used for the study of human diseases, we further compared these human gene regulatory networks with age-specific mouse brain transcriptomes. We discovered significantly preserved networks with both human aging and human disease and identified 17 shared genes in the top-ranked immune/microglia module, among which we found five human hub genes *TYROBP*, *FCER1G*, *ITGB2*, *MYO1F*, *PTPRC*, and two mouse hub genes *Trem2* and *C1qa*. Taken together, these results support the hypothesis that microglia are key players involved in human aging and neurodegenerative diseases, and suggest that mouse models should be appropriate for studying these microglial changes in human.

**Keywords:** aging, neurodegeneration, microglia, Alzheimer, Parkinson, bioinformatics, gene networks, WGCNA
# INTRODUCTION

Healthy aging is turning into a major societal challenge. Worldwide, over 960 million people are 60 years or older and this number is expected to double by the year 2050 (United Nations, Department of Economic and Social Affairs, Population Division, 2017). Furthermore, the risk of developing neurodegenerative disease increases with age (Niccoli and Partridge, 2012). The United States medical health system is already spending over $800 billion dollars for the therapy of disorders related to aging (Gooch et al., 2017), of which approximately one third is used for dementia-related diseases. While great advances have been made in unraveling the underlying molecular processes of the aging brain (Kirkwood, 2011; Mattson and Arumugam, 2018), so far only symptomatic treatments have been established, while curative treatments are still lacking.

Some neurodegenerative disorders are clearly of a heritable nature, but unfortunately, most of them are more complex and caused in a polygenic or even omnigenic manner (Boyle et al., 2017). Although genome-wide association studies (GWAS) have improved our understanding of common risk genes for neurodegenerative diseases (Bertram and Tanzi, 2009; Ramanan and Saykin, 2013; Pihlstrom et al., 2017; Pottier et al., 2018), a clearer understanding of how ‘healthy aging’ turns into neurological impairments (‘poor aging’) is needed. The expression of certain disease-associated genes is known to change with time and there might be a critical transition state before the onset of a neurodegenerative disease (Fjell et al., 2014), however, it is not clear whether complete networks break down or whether certain cell types drive healthy aging into poor aging (Kumar et al., 2013; Carmona-Gutierrez et al., 2016; Plaza-Zabala et al., 2017).

In the recent years, time-effective and low-cost high-throughput sequencing made it possible to analyze whole transcriptomes of humans and rodents even on a single-cell level. Here, we use a systems biology approach to compare publicly available gene expression datasets from healthy humans and patients affected by complex neurodegenerative disorders [Alzheimer's disease (AD), Huntington's disease (HD), and Parkinson's disease (PD)] across a broad spectrum of ages to address the question how 'healthy aging' turns into 'poor aging.' We obtain relevant transcriptional modules based on their co-expression relationship and further demonstrate that some modules are enriched in specific cell types such as microglia.

Since studying the human brain transcriptome is limited to post-mortem tissue analyses, most knowledge about brain aging comes from mouse or lower vertebrate models (Rosenthal and Brown, 2007). However, the cognitive decline in these models, which is the most accurate clinical hallmark for aging and neurodegeneration, is less severe than what is seen in humans. Thus, we further contrast our findings of human aging to data from healthy aging mouse brains and found similarly preserved networks and enrichments in cell type signature genes. Our study highlights microglia signatures in the center of common biological processes in human aging and neurodegenerative disease, as well as mouse aging.

# MATERIALS AND METHODS

## Data Availability

Our computational pipeline used in this publication is open-sourced and available at: https://github.com/smukher2/GithubFrontiersNeurosciDec2018.

## Criteria for Selection and Data Sources

Gene expression profile (microarray and RNA-seq) datasets available on Gene Expression Omnibus (GEO) were selected. For analysis of aging and disease, datasets from studies with multiple time points across age for control and neurodegenerative disease brains were selected. Human AD, HD, and PD disease studies that met these selection criteria, GSE33000 [microarray, prefrontal cortex, brain (Narayanan et al., 2014)] and GSE43490 [microarray, substantia nigra, brain (Corradini et al., 2014)] were analyzed. For cross species comparison, mouse hippocampus RNA-seq datasets from studies with multiple time points across age were selected: GSE61915 (Stilling et al., 2014), GSE73503 (Aaronson and Rosinski, 2015) and GSE83931 (Bundy et al., 2017). To investigate cell type signatures associated with aging and neurodegenerative diseases, an RNA-seq dataset taken from mouse brain (GSE52564, Zhang et al., 2014) encompassing microglia, astrocyte, neuron, oligodendrocyte precursors, newly formed oligodendrocytes, mature oligodendrocytes and endothelial cells was utilized.

## Microarray Data Preparation and Annotation

Two human aging and neurodegenerative disease microarray datasets were downloaded from GEO using GEOquery (Davis and Meltzer, 2007): GEO33000 (dual-color) and GEO43490 (single-color). Gene expression matrices were obtained from GEO, which consisted of normalized log₁₀ (Cy5/Cy3) test/reference for GSE33000 (raw 'CEL' files not available) and normalized log₂ signal intensity for GSE43490. Gene expression matrices were converted to linear space, quantile normalized with Lumi 2.32.0 (Du et al., 2008), and then log₂+1-scaled. ProbeIDs were annotated using the appropriate GPL file downloaded from GEO. The output metadata and gene expression files were merged in the same order. This resulted in a combined sample size of 637 human aging and human neurodegenerative disease data points.

## RNA-Seq Data Preparation and Annotation

Raw fastq RNA-seq reads were obtained for mouse hippocampus aging and mouse cerebral cortex cell types. GSE83931 and GSE73503 were total RNA samples prepared for single-end (NEBnext Ultra mRNA library preparation kit) and paired-end sequencing (standard Illumina library preparation protocols), respectively. GSE61915 and GSE52564 were polyA samples prepared for single-end sequencing using the TruSeq library preparation kit. Fastq files were processed on Cyverse's
Atmosphere cloud-computing platform led by Arizona State University, Tempe, AZ, United States (Goff et al., 2011; Merchant et al., 2016). Briefly, raw reads were obtained with sratoolkit.2.9.1 (Leinonen et al., 2011) and quality was evaluated with fastqc_v0.11.7 (Brown et al., 2017). We employed the tophat 2.1.1 aligner (Trapnell et al., 2009) to map reads to the mouse reference genome mm10 (GRCm38.92) from Ensembl (Zerbino et al., 2018) with the annotation file from GENCODE (GRCm38.92), using recommended/default tophat parameters for paired-end and single-end reads (parameters: -b2 -very-fast -no-coverage-search -no-novel-juncs), which generated bam files from samples aligned at ~90% or more. The mapped bam files were sorted with samtools-0.1.19 (Li et al., 2009). Sorted bam files were processed with HTSeq 0.10.0 (Anders et al., 2015) to obtain counts of reads aligned to genes (parameters: -r pos -t gene_name). The output metadata and gene expression htseq-count files were merged, normalized to the log2+1 scale using edgeR 3.22.5 (Robinson et al., 2009) and filtered for a cpm library size of >1 count in ≥3 samples, which resulted in a combined mouse aging sample size of 71. Additionally, for GSE52564, the sorted bam files were processed with cuffdiff2 from cufflinks_2.1.1-4 (Trapnell et al., 2012) to obtain a list of differentially expressed genes or cell type signature genes as defined by a fold change >20 in comparison of the given cell type to all other cell types. Our list was comparable to the previously published cell type signature gene list obtained from this same GSE52564 dataset (Zhang et al., 2014).

## Adjusting for Surrogate Variables Including Batch Effects

We used the Surrogate Variable analysis (SVA 3.28.0) R package (Leek and Storey, 2007; Leek et al., 2012) for removal of confounding effects from the combined human microarray expression data and combined mouse RNA-seq expression data (SVA model: ~study or batch + gender + disease + age; null model: ~study or batch + gender). Because tissue and study/batch were highly correlated, we only used one of these confounding factors (study/batch) for SVA modeling to avoid multicollinearity. For the estimation of variable effects by SVA the “be” method was used. SVA itself does not give the adjusted or normalized expression values and is therefore used in combination with the linear regression model (LM) fitting function from the Limma 3.36.5 R package (Ritchie et al., 2015). This resulted in SVA+LM normalized expression for age and neurodegenerative disease effects in human and mouse merged data, which was then used as input for weighted gene co-expression network analysis (WGCNA) separately for mouse and human.

## Data Used for WGCNA

Matching of unique microarray gene symbols between the human SVA+LM normalized datasets (GSE33000 and GSE43490) left 7301 common genes that were run through the WGCNA meta-analysis pipeline. For the mouse data, gene symbols were first mapped to human gene symbols using biomaRt

[version 2.36.1, (Kasprzyk, 2011)], then arranged by expression variability in descending order, and finally the same number of genes as in human (7301) were used as input for WGCNA.

## Meta-Analysis of Gene Expression Networks Using WGCNA

We conducted WGCNA meta-analysis (Miller et al., 2010) and utilized other functions of the WGCNA (Langfelder et al., 2008) R package (version 1.66) to compare co-expressed gene networks between (1) human aging and human neurodegenerative disease and (2) mouse aging. For both, genes were hierarchically clustered and groups of co-expressed genes (modules) were generated using the dynamic tree cutting algorithm with the following major parameters: soft power = 12, TOMType = “signed,” deepSplit = 2, pamStage = TRUE, cutHeight = 0.99 and minClusterSize = 30-3*deepSplit. Each module was assigned a unique color label. Hub genes were identified for each module using gene-module association KME values. Statistical significance was determined by regressing the traits onto the eigengene of the module, also called trait analysis (Plaisier et al., 2009), which revealed modules that were significantly associated with age (for mouse and human datasets) and disease (for human datasets).

## Identification of Cell Type Signatures in Gene Modules

To identify cell type signatures enriched in the aging and neurodegenerative disease preserved modules in mouse and human, we utilized the inbuilt hypergeometric test function “userListEnrichment” from the WGCNA 1.66 R package. We found overlapping genes between WGCNA modules and the set of cell type signature genes by pairwise comparison, and determined the significance of the overlaps with *p*-values adjusted with Bonferroni's method for multiple testing. Cell type signatures of microglial hub genes were confirmed using the RNA-Seq Data Navigator tool, which is freely available online as part of the Allen Cell Types Database (© 2015 Allen Institute for Brain Science)<fnref n="1" />. This was done by inputting hub genes in the “Gene Selection” box and building a “Group Plot” using default parameters. Cluster names were modified from the online tool to highlight broad cell types of origin, but no modifications were made to the resulting gene expression signatures. For mouse, cell types were grouped by “subclass” for ease of visualization.

## Cross-Species Comparison of Networks

We then assessed which modules, if any, were preserved between the human and mouse networks using module preservation Z statistics and *p*-values as described in the WGCNA package documentation (Langfelder et al., 2011). Similarly, we used this hypergeometric test described in the previous section to identify modules in the human and mouse networks with a significant number of overlapping genes. Together, these

<footnote n="1">¹http://celltypes.brain-map.org/rnaseq</footnote>
strategies allowed us to identify mouse and human modules that represent equivalent cell types of biological processes.

## KEGG Pathway and Reactome Pathway Analysis of Gene Expression Networks

The modules that were most preserved were annotated and characterized with KEGG pathways and Reactome pathways with the R package enrichR 1.0 (Jawaid, 2017). The top hits of the package output 'combined scores' and adjusted *p*-values were used to determine significance.

## RESULTS

The present article comprises results from the Systems Genetics of Neurodegeneration Summer School held in 2017 and combines analyses from multiple published high-throughput transcriptome datasets (Table 1) that were re-analyzed according to the experimental diagram shown in Figure 1.

We first combined two large microarray expression datasets for three neurodegenerative disorders: AD and HD with their respective controls (Narayanan et al., 2014), and PD with its

<table><caption>TABLE 1 | Overview of the datasets used.</caption><thead><tr><th>Dataset</th><th>Tissue</th><th>Design</th><th>Platform</th><th>GEO accession</th><th>Reference ID #</th><th>Reference</th></tr></thead><tbody><tr><td>DLPFC (BA9) brain tissues of AD patients, HD patients, and non-demented controls samples</td><td>Postmortem prefrontal cortex from Harvard Brain tissue resource center</td><td>624 individual DLPFC samples were profiled against a common DLPFC pool constructed from the same set of samples</td><td>Agilent 44K array (GPL4372)</td><td>GSE33000</td><td>GSM1423780 to GSM1424403</td><td>Narayanan et al., 2014</td></tr><tr><td>Transcriptomic profiles of controls and Parkinson's disease patients</td><td>Postmortem samples of dorsal motor vagal nucleus, locus coeruleus, and substantia nigra from controls and PD subjects from Brain Bank of the Brazilian Aging Brain Study Group, BEHEEC-FMUSP Note: In our present study, we only used the substantia nigra tissue data</td><td>Transcriptomic profiles of controls and Parkinson's disease patients were compared using SAM test for LC or Wilcoxon Mann-Whitney test for SN and VA (p < 0.005 and p < 0.01, respectively) in order to identify differentially expressed transcripts</td><td>Agilent-014850 Whole Human Genome Microarray 4x44K G4112F (Probe Name version)</td><td>GSE43490</td><td>GSM1294118 to GSM1294130</td><td>Corradini et al., 2014</td></tr><tr><td>Transcriptome database of eight cell types of mouse cerebral cortex</td><td>Isolated and purified neurons, astrocytes, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, endothelial cells, and pericytes from mouse cerebral cortex (different mouse lines)</td><td>RNA isolated from purified cell samples using a highly sensitive algorithm to detect alternative splicing in each gene were used to identify cell type enriched genes</td><td>Illumina HiSeq 2000 (Mus musculus) (GPL13112)</td><td>GSE52564</td><td>GSM1269903 to GSM1269916</td><td>Zhang et al., 2014</td></tr><tr><td>Transcriptome data of young and old mouse hippocampus</td><td>Mouse hippocampus of 3, 24, and 29 months old C57BL/6J mice</td><td>polyA-enriched RNA extracted from mouse hippocampus in three different age groups [3M vs. 24M (n = 5-6, single-end sequencing) and 3M vs. 29M (n = 3, paired-end sequencing)]</td><td>Illumina HiSeq 2000 (Mus musculus) (GPL13112)</td><td>GSE61915</td><td>SRR1593496 to SRR1593512</td><td>Stilling et al., 2014</td></tr><tr><td>Transcriptome data of knock-in mouse models of Huntington's disease</td><td>Mouse hippocampus of 2, 6, and 10 months old knock-in mice with CAG lengths of 20, 80, 92, 111, 140, 175 along with littermate control wild-type animals Note: In our present study, we only used WT data</td><td>mRNA expression profile of male and female (n = 4 each) mice of three different age points and seven different conditions (CAG lengths = normal, 20, 80, 92, 111, 140, 175)</td><td>Illumina HiSeq 2000 (Mus musculus) (GPL13112)</td><td>GSE73503</td><td>SRR2531532 to SRR2531555</td><td>Aaronson and Rosinski, 2015</td></tr><tr><td>Transcriptome data of the developing mouse hippocampus</td><td>Mouse hippocampus of 1, 2, and 4 months old B6 mice</td><td>Hippocampal mRNA from 1, 2, and 4 months old male and female B6 mice were analyzed by RNA sequencing of five biological replicates</td><td>Illumina HiSeq 2500 (Mus musculus) (GPL17021)</td><td>GSE83931</td><td>SRR3734796 to SRR3734825</td><td>Bundy et al., 2017</td></tr></tbody></table>
<figure><img src="image_2.png" /><figcaption>**FIGURE 1** | Experimental design and data pre-processing. Two microarray-based gene expression datasets of human patients including three neurodegenerative diseases and age-matched controls were combined and analyzed. Simultaneously, three RNA-sequencing-based gene expression datasets from mouse hippocampus were combined. The human and the mouse datasets were corrected for batch effects and gene networks were constructed with WGCNA. Species-specific networks were checked for overlap with cell type-specific RNA-sequencing data. AD, Alzheimer's disease; HD, Huntington's disease; PD, Parkinson's disease; SN, substantia nigra.</figcaption></figure>

A

Age distribution: PD, AD, HD and controls

<figure><img src="image_3.png" /></figure>

B

Original data: Human

<figure><img src="image_4.png" /></figure>

C

Batch-corrected data: Human

<figure><img src="image_5.png" /></figure>

D

WGCNA modules: Human

<figure><img src="image_6.png" /><figcaption>**FIGURE 2** | Removal of batch effects and gene network construction in human aging and neurodegeneration. (A) Age and gender distribution for the combined human dataset. (B,C) The different microarray datasets showed a strong batch effect, which was corrected using SVA to remove confounders while retaining age and disease as main discriminators (CON, healthy controls; ND, neurodegeneration). (D) WGCNA modules: Each row corresponds to one gene co-expression network (labeled by color). Numbers in the table indicate the Pearson correlation coefficient *r* and the associated *p*-values in parentheses. Coloring of the table encodes the correlation between each phenotype and each module eigengene (scale on the right).</figcaption></figure>
respective controls (Corradini et al., 2014) (Figure 1). These gene expression data were either taken from human post-mortem prefrontal cortex (AD + HD), or human post-mortem substantia nigra (PD). We obtained one large dataset of *n* = 637 patients and controls with balanced gender distribution and an age range of 18–106 years (Figure 2A). Since these data were originally gathered in different laboratories and originated from different brain areas, there was a strong batch effect (Figure 2B). Furthermore, neurodegenerative disease is known to be a strong confounder of aging, and vice-versa (Hung et al., 2010). Thus, we leveraged Surrogate Variable Analysis (Leek et al., 2012) and linear regression models to remove confounders (batch/tissue and gender, see the section “Materials and Methods,” Figure 1). Dimensionality reduction techniques showed that age and neurodegeneration were still retained as a primary components for further analysis, while there were no visible batch effects (Figure 2C). To define and contrast gene networks that were conserved between human aging and neurodegenerative disease, a co-expression gene network was then constructed that consisted of a normalized and combined AD, HD, and PD samples coded as “neurodegeneration” and respective controls. This approach increased power, but had the disadvantage to mask any potential tissue-specific networks.

Hierarchical clustering of expression data using WGCNA partitioned this large human data set into 12 modules (Figure 2D). These modules, which are designated in different colors, can be thought of as functionally different compartments of the human brain transcriptome, forming groups of highly interconnected transcripts that may shape a pathway that is relevant to aging and/or neurodegenerative disease (Fuller et al., 2007; Langfelder et al., 2008; Struebing et al., 2016, 2018). No module had any relationship to the confounding variables study, gender or tissue, validating our batch adjustment strategy (Figure 2D). In order to gain more insight on the biological function of these modules, we overlapped genes within each of the previously found 12 co-expression modules with a signature gene list that was specific for several cortical cell types. These included neurons, astrocytes, oligodendrocyte precursor cells, newly formed oligodendrocytes, myelinating oligodendrocytes, microglia, endothelial cells and pericytes (Zhang et al., 2014). We then tested each module-cell type association for significance using a hypergeometric test. From all mutual comparisons, only three reached significance after adjusting for multiple comparisons (Supplementary Table S1): The human “black” module was related to microglia (*p* = 8.00e-60), whereas the human “blue” module was enriched in endothelial cells

<figure><img src="image_7.png" /></figure>
module also demonstrated a highly significant hypergeometric enrichment (p = 8e-16; **Supplementary Table S3**) with the mouse "yellow" module, which we previously showed was enriched with microglia genes as well (p = 2.28e-37). Functional annotations of the microglia-related human "black" and mouse "yellow" modules showed similar enrichments in immune system-related processes (**Figure 4**). Besides the Reactome identifiers "immune system" (human p = 3.2e-27, mouse p = 9.33e-11) and "innate immune system" (human p = 6.4e-11, mouse p = 4.7e-06), several infectious disease pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) were identified. These results suggest that the main actors driving the transcriptome changes associated with aging in human and mice are immune system-related genes expressed in microglia.

The human "black" module and the mouse "yellow" module had 17 genes in common (**Figure 5A** and **Supplementary Table S3**). All of these overlapping genes were also microglia signature genes, out of which seven were among the top ten hub genes in human or mouse (**Tables 2, 3** and detailed: **Supplementary Tables S3, S4**), meaning that they displayed some of the most representative gene expression patterns of the module as a whole and are likely good markers of microglia or their associated processes. FCER1G (transmembrane signaling receptor), ITGB2 (=CD18, subunit of complement receptor 3), MYO1F (required in cytoskeleton remodeling and migration) (Kalhammer and Bähler, 2000; Chen and Iijima,

2012; Maravillas-Montero and Santos-Argumedo, 2012), PTPRC (=CD45 or B220, signaling molecule) and TYROBP (=DAP12, activatory adaptor protein) were hub genes in the human "black" module, while C1qa (subunit of complement component C1q) and Trem2 (transmembrane signaling receptor) were hub genes within the mouse "yellow" module (**Figure 5A**). By directly interrogating their expression changes during the aging process, we found significant associations of these hub genes with human aging in the "black" module (**Figure 5B** and **Tables 2, 3**). All of these genes were strongly enriched in microglia, as evidenced by comparison with our cell type-specific gene list (**Supplementary Tables S1, S2**). Cross-referencing them with RNA-seq data from the Allen Cell Types Database (**Figure 5C**) further confirmed the cellular identity of these modules. These results suggest that even though gene networks enriched in microglial signatures that potentially direct human and mouse aging show significant overlap and are preserved, they are controlled by different regulatory hubs in the two species.

In conclusion, our analysis identifies co-expressed networks involved in human aging and neurodegeneration that are well preserved in mouse aging and associated with similar molecular pathways. A large part of this conserved regulatory pathway is mediated by microglia. Ultimately, our analyses suggest that the mouse is a valid model system to study changes in microglia during human aging and neurodegeneration.

<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Microglia-specific hub genes in mouse and human and enrichment with the Allen Brain Cell Types Database. (A) The human "black" module and the mouse "yellow" module shared 17 genes, out of which five were hub genes in human and two were hub genes in mouse. All 17 shared genes were microglia signature genes as confirmed with the published cell type signature genes. (B) These hub genes correlated positively with age in both mouse and human, while mean expression is higher in human neurodegeneration (ND) than in controls (CON). Each gray dot represents one sample and gene, whereas the blue (ND) and the red (CON) dots indicate means and the colored lines represent a linear fit with a gray-shaded 95% confidence interval. (C) The hub genes were enriched in human (left) and mouse (right) microglia. Gene symbols are given on the y-axis, while different cell types as defined in the Allen Cell Types Database are displayed on the x-axis. Dot size and color correspond to the fraction of cells in a cluster expressing a given gene and the median gene expression, respectively. CPM, counts per million; MTG, middle temporal gyrus; VISp, primary visual cortex; ALM, anterior lateral motor cortex; Inh, inhibitory neurons; Exc, excitatory neurons; oligodendrocytes, LM, leptomeningeal; IT, intratelencephalic; PT, pyramidal tract; NP, near-projecting; CT, corticothalamic. Microglia are highlighted with a red box.</figcaption></figure>
<table><caption>**TABLE 2** | Overview of top 10 hub genes in human "black" module.</caption><thead><tr><th>Rank</th><th>Gene</th><th>CorrDisease</th><th>PvalDisease</th><th>CorrAge</th><th>PvalAge</th><th>CorrGene</th><th>PvalGene</th></tr></thead><tbody><tr><td>1</td><td>TBXAS1</td><td>0.399</td><td>1.03E-25</td><td>0.200</td><td>3.56E-07</td><td>0.882</td><td>6.87E-210</td></tr><tr><td>2</td><td><strong>TYROBP</strong></td><td>0.397</td><td>2.03E-25</td><td>0.085</td><td>0.032</td><td>0.872</td><td>5.07E-199</td></tr><tr><td>3</td><td><strong>PTPRC</strong></td><td>0.327</td><td>2.35E-17</td><td>0.244</td><td>4.41E-10</td><td>0.871</td><td>3.43E-198</td></tr><tr><td>4</td><td><strong>ITGB2</strong></td><td>0.395</td><td>3.01E-25</td><td>0.194</td><td>8.27E-07</td><td>0.859</td><td>1.76E-186</td></tr><tr><td>5</td><td><strong>MYO1F</strong></td><td>0.380</td><td>2.76E-23</td><td>0.102</td><td>0.010</td><td>0.846</td><td>9.05E-176</td></tr><tr><td>6</td><td>LST1</td><td>0.286</td><td>1.94E-13</td><td>0.157</td><td>7.10E-05</td><td>0.841</td><td>2.36E-171</td></tr><tr><td>7</td><td>CYBA</td><td>0.435</td><td>9.96E-31</td><td>0.282</td><td>4.43E-13</td><td>0.840</td><td>1.59E-170</td></tr><tr><td>8</td><td>SLC7A7</td><td>0.371</td><td>3.69E-22</td><td>0.133</td><td>0.001</td><td>0.837</td><td>4.71E-168</td></tr><tr><td>9</td><td>RNASET2</td><td>0.403</td><td>2.65E-26</td><td>0.189</td><td>1.63E-06</td><td>0.835</td><td>3.72E-167</td></tr><tr><td>10</td><td><strong>FCER1G</strong></td><td>0.456</td><td>5.16E-34</td><td>0.246</td><td>2.95E-10</td><td>0.829</td><td>2.60E-162</td></tr></tbody></table>

Rank, rank of hub gene (ordered descending by CorrGene); gene, gene symbol; CorrDisease, Pearson correlation of human "black" module eigengene with a numeric vector indicating whether the sample comes from a control (0) or disease (1) donor; PvalDisease, P-value of correlation in CorrDisease; CorrAge, Pearson correlation of human "black" module eigengene with donor age; PvalAve, P-values of correlation in CorrAge; CorrGene, Pearson correlation of human "black" module eigengene with expression level of the gene in the "Gene" column; PvalGene, P-values of correlation in CorrGene. Genes labeled in bold are microglia signature genes that are shared by the human "black" module and the mouse "yellow" module.

<table><caption>**TABLE 3** | Overview of top 10 hub genes in mouse "yellow" module.</caption><thead><tr><th>Rank</th><th>Gene</th><th>CorrAge</th><th>PvalAge</th><th>CorrGene</th><th>PvalGene</th></tr></thead><tbody><tr><td>1</td><td>C4B</td><td>0.652</td><td>7.39E-10</td><td>0.940</td><td>6.58E-34</td></tr><tr><td>2</td><td>C4A</td><td>0.652</td><td>7.39E-10</td><td>0.940</td><td>6.58E-34</td></tr><tr><td>3</td><td>CTSS</td><td>0.595</td><td>4.50E-08</td><td>0.900</td><td>1.41E-26</td></tr><tr><td>4</td><td>LGALS3BP</td><td>0.529</td><td>2.12E-06</td><td>0.874</td><td>2.46E-23</td></tr><tr><td>5</td><td>APOD</td><td>0.561</td><td>3.58E-07</td><td>0.874</td><td>2.54E-23</td></tr><tr><td>6</td><td>C1QC</td><td>0.521</td><td>3.18E-06</td><td>0.872</td><td>3.99E-23</td></tr><tr><td>7</td><td><strong>TREM2</strong></td><td>0.567</td><td>2.59E-07</td><td>0.857</td><td>1.37E-21</td></tr><tr><td>8</td><td>C1QB</td><td>0.554</td><td>5.54E-07</td><td>0.855</td><td>2.35E-21</td></tr><tr><td>9</td><td>MPEG1</td><td>0.568</td><td>2.33E-07</td><td>0.852</td><td>4.84E-21</td></tr><tr><td>10</td><td><strong>C1QA</strong></td><td>0.548</td><td>7.71E-07</td><td>0.849</td><td>8.53E-21</td></tr></tbody></table>

Rank, rank of hub gene (ordered descending by CorrGene); gene, gene symbol; CorrAge, Pearson correlation of mouse "yellow" module eigengene with donor age; PvalAve, P-values of correlation in CorrAge; CorrGene, Pearson correlation of mouse "yellow" module eigengene with expression level of the gene in the "Gene" column; PvalGene, P-values of correlation in CorrGene. Genes labeled in bold are microglia signature genes that are shared by the human "black" module and the mouse "yellow" module.

## DISCUSSION

There is growing evidence suggesting that aging and neurodegenerative diseases show common pathological gene regulation patterns (Chakrabarti and Mohanakumar, 2016). While different cell types are clearly affected differently from one another, the degree to which each cell type is involved in the dysregulation of genetic networks is not entirely understood. Therefore, the overall goal of this study was to identify cell type signature(s) of aging and neurodegenerative diseases by comparing their gene regulatory networks. Through a systematic WGCNA-based analysis, we identified microglia as a major cell type directing gene networks related to aging and three neurodegenerative diseases (AD, HD, and PD) in post-mortem human brains. Moreover, we found that this role of microglia in aging was conserved in mice.

Combining transcriptomic data from different sources and different neurodegenerative diseases enabled us to investigate a large dataset that covered different diseases and a representative age range from 18 to 106 years. Nevertheless, using datasets from different sources often introduces unwanted batch effects. To overcome this hurdle, we leveraged a unique combination of batch-correcting statistical methods (SVA and linear regression models) to segregate the effects of age and neurodegeneration in their respective molecular networks.

We first identified three gene co-expression networks that were significantly correlated to aging and neurodegenerative diseases. Each of these was enriched in genes specific for a distinct cell type, namely neurons (human "turquoise" module), endothelial cells (human "blue" module), and microglia (human "black" module). The association with neurons is consistent with the major feature of neuronal damage, while the relation to endothelial cells could implicate neurovascular dysfunction, which is also seen in aging and neurodegeneration (Kisler et al., 2017). Although a contribution of glial cell types to these conditions was expected, it was surprising to find the strongest association with microglia, while there was only non-significant overlap with astrocytes and no association with oligodendrocytes. Functional annotations of the microglia-associated human and mouse modules showed that it was strongly related to the immune system. Taken together, these results suggest that microglia are the major cell type common to aging and three major neurodegenerative diseases (AD, HD, PD).

Cross-species comparison of gene co-expression networks can provide additional information about network convergence and divergence of networks, which can help to identify key members of the network and species-specific differences (Hawrylycz et al., 2012). As mouse models are widely used in the study of human aging and diseases (Klopstock et al., 2011), we mined mouse aging data to test the cross-species relevance of our finding that microglia signature genes were strongly associated with aging as well as neurodegeneration. In our results, the mouse "yellow" module was significantly preserved between mouse aging and human aging. Furthermore, this module also
showed a significant enrichment in microglial genes, suggesting that aging and neurodegenerative disease-associated gene co-expression networks are conserved between mouse and human and associated with microglia.

Comparing the human and mouse microglia-related modules yielded 17 shared genes, of which seven were hub genes: FCER1G, ITGB2, MYO1F, PTPRC, and TYROBP in human and C1qa and Trem2 in mouse. Recent findings support the role of these genes in aging or neurodegeneration: For example, C1qa as part of the complement system mediates early synapse loss in AD, while TYROBP expression is increased in AD and simultaneously increases the expression of AD-related genes such as Cd33 in an AD mouse model, consistent with its role being a hub gene (Zhang et al., 2013; Herms and Dorostkar, 2016; Haure-Mirande et al., 2018). Furthermore, TREM2 mutations are well known to increase AD risk in humans, and increased Trem2 expression in microglia has been found to correlate with amyloid deposition (Brendel et al., 2017). The intersection of hub genes with the Allen Brain Cell Types Database showed that they were indeed strongly expressed in microglia. Consistent with their role as being major gene network drivers, their expression pattern significantly correlated with aging (mouse and human) and neurodegeneration (human). Even though hub genes were present in both modules, they assumed this role in only one of the species and not both. This could indicate a form of network re-wiring (Halfon, 2017), suggesting that the aging- and neurodegeneration-related networks are evolutionarily dynamic. These results are also consistent with previous findings of convergent and divergent pathways in human and mouse brain transcriptomes using a similar approach (Miller et al., 2010).

Several caveats concerning this study such as residual batch effects or tissue heterogeneity should be kept in mind. As is the case for any study that is correlational in nature, inferences that imply a causal relationship must be critically evaluated. Ultimately, biological validation will be required in both mouse and human organoid models of aging, AD, PD, and HD, to ultimately prove their dysfunctional regulation in aging and neurodegeneration.

## CONCLUSION

In conclusion, we have provided a systems biology-based study that supports two ideas: (i) that there is a common molecular mechanism underlying aging, AD, HD, and PD that might be related to microglia, and (ii) that these microglial gene networks associated with neurodegeneration and aging are conserved in human and mice. Specifically, we have identified seven candidate microglia genes that should be the focus of further investigations, for example in mouse models and human

organoid models. Understanding their transition during aging and in neurodegenerative diseases might be of particular interest for future diagnostics.

## AUTHOR CONTRIBUTIONS

The initial idea for the project was formed by all authors during the 2017 Systems Genetics of Neurodegeneration summer school, Germany. SM performed and guided the majority of data analysis. FS assisted in the analysis and facilitated manuscript submission by preparation of figures. CK created strategies to enable continuation of effective collaboration and prepared main tables for the manuscript. SM, CK, and FS jointly wrote the manuscript. MP-J and JM gave statistical and methodological advice. All authors contributed significantly with regular reviews of the experimental design, interpretation of the results and manuscript preparation.

## FUNDING

Cyverse Atmosphere cloud-computing platform was supported by the National Science Foundation under Award Numbers DBI-0735191 and DBI-1265383. FS was supported by NIH grant T32EY007092-30 at Emory University and is now funded under the European Union's Horizon 2020 Research and Innovation Program (grant agreement MSCA-IF-EF-RI 792832).

## ACKNOWLEDGMENTS

We would like to acknowledge all organizers of the 2017 Systems Genetics of Neurodegeneration summer school, namely Rupert Overall, Robert W. Williams, Gerd Kempermann, and Rudi Balling, furthermore the nuns of the Abtei Frauenchiemsee, headed by Sister Scholastica, for hosting us. We would like to thank Moeen Riaz, Victoria Perreau, and Karl Broman for their input in developing the idea for this project. Cloud computing processing time was obtained free of cost with special permission from Arizona State University, Tempe, AZ, United States. We also thank the Allen Institute Founders, Paul G. Allen, and Jody Allen, for their vision, encouragement and support.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2019.00002/full#supplementary-material

## REFERENCES

Aaronson, J., and Rosinski, J. (2015). Transcriptome profiling in knock-in mouse models of Huntington's disease [Hippocampus_mRNA]. Available at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73503

Anders, S., Pyl, P. T., and Huber, W. (2015). HTSeq-a python framework to work with high-throughput sequencing data. *Bioinformatics* 31, 166-169. doi: 10.1093/bioinformatics/btu638

Bertram, L., and Tanzi, R. E. (2009). Genome-wide association studies in Alzheimer's disease. *Hum. Mol. Genet* 18, R137-R145. doi: 10.1093/hmg/ddp406
Boyle, E. A., Li, Y. I., and Pritchard, J. K. (2017). An expanded view of complex traits: from polygenic to omnigenic. *Cell* 169, 1177-1186. doi: 10.1016/j.cell.2017.05.038

Brendel, M., Kleinberger, G., Probst, F., Jaworska, A., Overhoff, F., Blume, T., et al. (2017). Increase of TREM2 during aging of an Alzheimer's disease mouse model is paralleled by microglial activation and amyloidosis. *Front. Aging Neurosci.* 9:8. doi: 10.3389/fnagi.2017.00008

Brown, J., Pirrung, M., and McCue, L. A. (2017). FQC dashboard: integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. *Bioinformatics* 33, 3137-3139. doi: 10.1093/bioinformatics/btx373

Bundy, J. L., Vied, C., and Nowakowski, R. S. (2017). Sex differences in the molecular signature of the developing mouse hippocampus. *BMC Genomics* 18:237. doi: 10.1186/s12864-017-3608-7

Carmona-Gutierrez, D., Hughes, A. L., Madeo, F., and Ruckenstuhl, C. (2016). The crucial impact of lysosomes in aging and longevity. *Ageing Res. Rev.* 32, 2-12. doi: 10.1016/j.arr.2016.04.009

Chakrabarti, S., and Mohanakumar, K. P. (2016). Aging and neurodegeneration: a tangle of models and mechanisms. *Aging Dis.* 7, 111-113. doi: 10.14336/AD.2016.0312

Chen, C.-L., and Iijima, M. (2012). Myosin I: a new pip(3) effector in chemotaxis and phagocytosis. *Commun. Integr. Biol.* 5, 294-296. doi: 10.4161/cib.19892

Corradini, B. R., Iamashita, P., Tampellini, E., Farfel, J. M., Grinberg, L. T., and Moreira-Filho, C. A. (2014). Complex network-driven view of genomic mechanisms underlying Parkinson's disease: analyses in dorsal motor vagal nucleus, locus coeruleus, and substantia nigra. *Biomed. Res. Int.* 2014, 543673. doi: 10.1155/2014/543673

Davis, S., and Meltzer, P. S. (2007). GEOquery: a bridge between the Gene Expression Omnibus (GEO) and bioConductor. *Bioinformatics* 23, 1846-1847. doi: 10.1093/bioinformatics/btm254

Du, P., Kibbe, W. A., and Lin, S. M. (2008). Lumi: a pipeline for processing illumina microarray. *Bioinformatics* 24, 1547-1548. doi: 10.1093/bioinformatics/btn224

Fjell, A. M., McEvoy, L., Holland, D., Dale, A. M., and Walhovd, K. B. (2014). What is normal in normal aging? Effects of aging, amyloid and Alzheimer's disease on the cerebral cortex and the hippocampus. *Prog. Neurobiol.* 117, 20-40. doi: 10.1016/j.pneurobio.2014.02.004

Fuller, T. F., Ghazalpour, A., Aten, J. E., Drake, T. A., Lusis, A. J., and Horvath, S. (2007). Weighted gene coexpression network analysis strategies applied to mouse weight. *Mamm. Genome* 18, 463-472. doi: 10.1007/s00335-007-9043-3

Goff, S. A., Vaughn, M., McKay, S., Lyons, E., Stapleton, A. E., Gessler, D., et al. (2011). The iPlant collaborative: cyberinfrastructure for plant biology. *Front. Plant Sci.* 2:34. doi: 10.3389/fpls.2011.00034

Gooch, C. L., Pracht, E., and Borenstein, A. R. (2017). The burden of neurological disease in the united states: a summary report and call to action. *Ann. Neurol.* 81, 479-484.

Halfon, M. S. (2017). Perspectives on gene regulatory network evolution. *Trends Genet* 33, 436-447. doi: 10.1016/j.tig.2017.04.005

Haure-Mirande, J.-V., Wang, M., Audrain, M., Fanutza, T., Kim, S. H., Heja, S., et al. (2018). Integrative approach to sporadic Alzheimer's disease: deficiency of TYROBP in cerebral Aβ amyloidosis mouse normalizes clinical phenotype and complement subnetwork molecular pathology without reducing Aβ burden. *Mol. Psychiatry* doi: 10.1038/s41380-018-0255-6 [Epub ahead of print].

Hawrylycz, M. J., Lein, E. S., Guillozet-Bongaarts, A. L., Shen, E. H., Ng, L., Miller, J. A., et al. (2012). An anatomically comprehensive atlas of the adult human brain transcriptome. *Nature* 489, 391-399. doi: 10.1038/nature11405

Herms, J., and Dorostkar, M. M. (2016). Dendritic spine pathology in neurodegenerative diseases. *Annu. Rev. Pathol. Mech. Dis.* 11, 221-250. doi: 10.1146/annurev-pathol-012615-044216

Hung, C. W., Chen, Y. C., Hsieh, W. L., Chiou, S. H., and Kao, C. L. (2010). Ageing and neurodegenerative diseases. *Ageing Res. Rev.* 9(Suppl. 1), S36-S46. doi: 10.1016/j.arr.2010.08.006

Jawaid, W. (2017). *enrichR: Provides an R Interface to "Enrichr" R Package Version 1.0.*

Kalhammer, G., and Bähler, M. (2000). Unconventional myosins. *Essays Biochem.* 35, 33-42. doi: 10.1042/bse0350033

Kasprzyk, A. (2011). *BioMart: driving a paradigm change in biological data management*. *Database* 2011:bar049. doi: 10.1093/database/bar049

Kirkwood, T. B. L. (2011). Systems biology of ageing and longevity. *Philos. Trans. R. Soc. B Biol. Sci.* 366, 64-70. doi: 10.1098/rstb.2010.0275

Kisler, K., Nelson, A. R., Montagne, A., and Zlokovic, B. V. (2017). Cerebral blood flow regulation and neurovascular dysfunction in Alzheimer disease. *Nat. Rev. Neurosci.* 18, 419-434. doi: 10.1038/nrn201748

Klopstock, T., Elstner, M., and Bender, A. (2011). Creatine in mouse models of neurodegeneration and aging. *Amino Acids* 40, 1297-1303. doi: 10.1007/s00726-011-0850-1

Kumar, A., Gibbs, J. R., Beilina, A., Dillman, A., Kumaran, R., Trabzuni, D., et al. (2013). Age-associated changes in gene expression in human brain and isolated neurons. *Neurobiol. Aging* 34, 1199-1209. doi: 10.1016/j.neurobiolaging.2012.10.021

Langfelder, P., Horvath, S., Fisher, R., Zhou, X., Kao, M., Wong, W., et al. (2008). WGCNA: an R package for weighted correlation network analysis. *BMC Bioinformatics* 9:559. doi: 10.1186/1471-2105-9-559

Langfelder, P., Luo, R., Oldham, M. C., and Horvath, S. (2011). Is my network module preserved and reproducible? *PLoS Comput. Biol.* 7:e1001057. doi: 10.1371/journal.pcbi.1001057

Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., and Storey, J. D. (2012). The SVA package for removing batch effects and other unwanted variation in high-throughput experiments. *Bioinformatics* 28, 882-883. doi: 10.1093/bioinformatics/bts034

Leek, J. T., and Storey, J. D. (2007). Capturing heterogeneity in gene expression studies by surrogate variable analysis. *PLoS Genet.* 3, 1724-1735. doi: 10.1371/journal.pgen.0030161

Leinonen, R., Sugawara, H., and Shumway, M. (2011). The sequence read archive. *Nucleic Acids Res.* 39, D19-D21. doi: 10.1093/nar/gkq1019

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment / map format and SAMtools. *Bioinformatics* 25, 2078-2079. doi: 10.1093/bioinformatics/btp352

Maravillas-Montero, J. L., and Santos-Argumedo, L. (2012). The myosin family: unconventional roles of actin-dependent molecular motors in immune cells. *J. Leukoc. Biol.* 91, 35-46. doi: 10.1189/jlb.0711335

Mattson, M. P., and Arumugam, T. V. (2018). Hallmarks of brain aging: adaptive and pathological modification by metabolic states. *Cell Metab.* 27, 1176-1199. doi: 10.1016/j.cmet.2018.05.011

Merchant, N., Lyons, E., Goff, S., Vaughn, M., Ware, D., Micklos, D., et al. (2016). The Iplant collaborative: cyberinfrastructure for enabling data to discovery for the life sciences. *PLoS Biol.* 14:e1002342. doi: 10.1371/journal.pbio.1002342

Miller, J. A., Horvath, S., and Geschwind, D. H. (2010). Divergence of human and mouse brain transcriptome highlights Alzheimer disease pathways. *Proc. Natl. Acad. Sci. U.S.A.* 107, 12698-12703. doi: 10.1073/pnas.0914257107

Narayanan, M., Huynh, J. L., Wang, K., Yang, X., Yoo, S., McElwee, J., et al. (2014). Common dysregulation network in the human prefrontal cortex underlies two neurodegenerative diseases. *Mol. Syst. Biol.* 10, 743-743. doi: 10.15252/msb.20145304

Niccoli, T., and Partridge, L. (2012). Ageing as a risk factor for disease. *Curr. Biol.* 22, R741-R752. doi: 10.1016/j.cub.2012.07.024

Oldham, M. C., Horvath, S., and Geschwind, D. H. (2006). Conservation and evolution of gene coexpression networks in human and chimpanzee brains. *Proc. Natl. Acad. Sci. U.S.A.* 103, 17973-17978. doi: 10.1109/APS.1998.690857

Pihlstrom, L., Wiethoff, S., and Houlden, H. (2017). Genetics of neurodegenerative diseases: an overview. *Handb. Clin. Neurol.* 145, 309-323. doi: 10.1016/B978-0-12-802395-2.00022-5

Plaisier, C. L., Horvath, S., Huertas-Vazquez, A., Cruz-Bautista, I., Herrera, M. F., Tusie-Luna, T., et al. (2009). A systems genetics approach implicates USF1, FADS3, and other causal candidate genes for familial combined hyperlipidemia. *PLoS Genet.* 5:e1000642. doi: 10.1371/journal.pgen.1000642

Plaza-Zabala, A., Sierra-Torre, V., and Sierra, A. (2017). Autophagy and microglia: novel partners in neurodegeneration and aging. *Int. J. Mol. Sci.* 18:e598. doi: 10.3390/ijms18030598

Pottier, C., Zhou, X., Perkerson, R. B., Baker, M., Jenkins, G. D., Serie, D. J., et al. (2018). Potential genetic modifiers of disease risk and age at onset in patients with frontotemporal lobar degeneration and GRN mutations: a genome-wide association study. *Lancet Neurol.* 17, 548-558. doi: 10.1016/S1474-4422(18)30126-1
Ramanan, V. K., and Saykin, A. J. (2013). Pathways to neurodegeneration: mechanistic insights from GWAS in Alzheimer's disease, Parkinson's disease, and related disorders. *Am. J. Neurodegener. Dis.* 2, 145-175.

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* 43:e47. doi: 10.1093/nar/gkv007

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2009). edgeR: A bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* 26, 139-140. doi: 10.1093/bioinformatics/btp616

Rosenthal, N., and Brown, S. (2007). The mouse ascending: perspectives for human-disease models. *Nat. Cell Biol.* 9, 993-999. doi: 10.1038/ncb437

Stilling, R. M., Benito, E., Barth, J., Gertig, M., Capece, V., Burckhardt, S., et al. (2014). De-regulation of gene expression and alternative splicing affects distinct cellular pathways in the aging hippocampus. *Front. Cell. Neurosci.* 8:373. doi: 10.3389/fncel.2014.00373

Struebing, F. L., King, R., Li, Y., Chrenek, M. A., Lyuboslavsky, P. N., Sidhu, C. S., et al. (2018). Transcriptional changes in the mouse retina after ocular blast injury: a role for the immune system. *J. Neurotrauma* 35, 118-129. doi: 10.1089/neu.2017.5104

Struebing, F. L., Lee, R. K., Williams, R. W., and Geisert, E. E. (2016). Genetic networks in mouse retinal ganglion cells. *Front. Genet.* 7:169. doi: 10.3389/fgene.2016.00169

Trapnell, C., Pachter, L., and Salzberg, S. L. (2009). TopHat: discovering splice junctions with RNA-Seq. *Bioinformatics* 25, 1105-1111. doi: 10.1093/bioinformatics/btp120

Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments

with TopHat and Cufflinks. *Nat. Protoc.* 7, 562-578. doi: 10.1038/jid.2014.371

United Nations, Department of Economic and Social Affairs, Population Division (2017). *World Population Ageing 2017* (ST/ESA/SER.A/408). Available at: http://www.un.org/en/development/desa/population/publications/pdf/ageing/WPA2017_Report.pdf

Zerbino, D. R., Achuthan, P., Akanni, W., Amode, M. R., Barrell, D., Bhai, J., et al. (2018). Ensembl 2018. *Nucleic Acids Res.* 46, D754–D761. doi: 10.1093/nar/gkx1098

Zhang, B., Gaiteri, C., Bodea, L. G., Wang, Z., McElwee, J., Podtelezhnikov, A. A., et al. (2013). Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease. *Cell* 153, 707-720. doi: 10.1016/j.cell.2013.03.030

Zhang, Y., Chen, K., Sloan, S. A., Bennett, M. L., Scholze, A. R., O'Keeffe, S., et al. (2014). An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. *J. Neurosci.* 34, 11929-11947. doi: 10.1523/JNEUROSCI.1860-14.2014

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Mukherjee, Klaus, Pricop-Jeckstadt, Miller and Struebing. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.