Skip to main content

Ocular hypertension suppresses homeostatic gene expression in optic nerve head microglia of DBA/2 J mice

Abstract

Glaucoma is the leading cause of irreversible vision loss. Ocular hypertension is a major risk factor for glaucoma and recent work has demonstrated critical early neuroinflammatory insults occur in the optic nerve head following ocular hypertension. Microglia and infiltrating monocytes are likely candidates to drive these neuroinflammatory insults. However, the exact molecular identity / transcriptomic profile of microglia following ocular hypertensive insults is unknown. To elucidate the molecular identity of microglia after long-term exposure to ocular hypertension, we used a mouse model of glaucoma (DBA/2 J). We performed RNA-sequencing of microglia mRNA from the optic nerve head at a time point following ocular hypertensive insults, but preceding detectable neurodegeneration (with microglia identified as being CD45lo/CD11b+/CD11c). Furthermore, RNA-sequencing was performed on optic nerve head microglia from mice treated with radiation therapy, a potent therapy preventing neuroinflammatory insults. Transcriptomic profiling of optic nerve head microglia mRNA identifies metabolic priming with marked changes in mitochondrial gene expression, and changes to phagocytosis, inflammatory, and sensome pathways. The data predict that many functions of microglia that help maintain tissue homeostasis are affected. Comparative analysis of these data with data from previously published whole optic nerve head tissue or monocyte-only samples from DBA/2 J mice demonstrate that many of the neuroinflammatory signatures in these data sets arise from infiltrating monocytes and not reactive microglia. Finally, our data demonstrate that prophylactic radiation therapy of DBA/2 J mice potently abolishes these microglia metabolic transcriptomic changes at the same time points. Together, our data provide a unique resource for the community to help drive further hypothesis generation and testing in glaucoma.

Background

Glaucoma is one of the most common neurodegenerations affecting an estimated 80 million people worldwide [1]. It is a complex and multifactorial disease characterised by the progressive dysfunction and loss of retinal ganglion cells and their axons (that make up the optic nerve). A common theme between animal models and human glaucoma is the activation or reactivity of glial cells in the retina, optic nerve, and optic nerve head [2,3,4,5,6,7]. Activated microglia are known to affect the progression of neurodegenerative diseases due to their influence over homeostatic and immune responses. Ongoing microglial responses may be protective, damage neural tissue, and/or lead to chronic inflammation. Given the primary role of microglia in regulating neuroinflammation (such as in Alzheimer’s disease [8,9,10,11]), identifying any dysregulation of microglia is of paramount importance for the development of neuroprotective treatments [12].

Microglia perform a diverse series of functions to support neural activity, including maintenance of synapses and axons, removal of cellular debris, surveillance for injury and pathogens, and co-ordination of neuroinflammatory responses [12,13,14,15,16,17]. These functions require microglia to continuously sense and respond to their environment. Many environmental cues lead to changes in microglial gene expression that support functional specialization. Genome-wide gene expression profiling has been used to identify important functional differences between microglia in normal physiological, neuroinflammatory, and neurodegenerative conditions [18,19,20]. However, microglia-specific transcriptomic profiles in the normal and glaucomatous optic nerve head are yet to be generated. These previous studies support the need for transcriptomic profiling of microglia in glaucoma relevant tissues.

The DBA/2 J mouse develops neurodegeneration due to elevated intraocular pressure (IOP) with hallmark features of an inherited, chronic human glaucoma. An important, human relevant aspect of disease in DBA/2 J mice is elevated IOP-related injury at the optic nerve head (ONH). Elevated IOP is a major risk factor for human glaucoma. In our colony, elevated IOP develops from 6 months of age, and by 9 months of age there are signs of injury at the ONH [21, 22]. The neurodegeneration it triggers is complete in the majority of eyes by 12 months of age (based on retinal ganglion cell soma and axon loss) [22]. We have previously profiled genome wide gene expression in the whole ONH from 8.5 to 10.5 months of age to investigate disease progression due to elevated IOP before and during neurodegeneration [2, 3]. Notably these data pointed to major neuroinflammatory changes at a pre-degenerative stage of disease, including changes in the complement system, endothelin system, and cell-adhesion pathways. Our group and others have experimentally validated molecules in these pathways to demonstrate their importance and role in glaucoma (from mouse to human [2, 12, 23,24,25,26,27,28,29];). In addition, these pathways have been independently identified in other models of glaucomatous insult and retinal ganglion cell death (ocular hypertension, optic nerve crush, optic nerve transection / axotomy [30,31,32]) supporting the utility of DBA/2 J retina and optic nerve head tissue for modelling glaucoma. Despite the importance of neuroinflammation in DBA/2 J glaucoma, it remains unclear how microglia contribute to disease progression.

To investigate the role of microglia in the ONH in response to ocular hypertension, we performed RNA-sequencing of mRNA from optic nerve head microglia at a pre-degenerative stage of disease (9 months of age prior to detectable optic nerve degeneration). Furthermore, we performed transcriptomic profiling of microglia from radiation-treated mice, a potent and robust anti-inflammatory and neuroprotective therapy for DBA/2 J glaucoma [3, 33, 34]. Our results provide new information about potential dysfunction of these cells. We expect these data to be a novel resource for the glaucoma and neurodegeneration community.

Methods

Mouse strain, breeding and husbandry

Mice were housed and fed in a 14 h light / 10 h dark cycle with food and water available ad libitum; all mice used in the study were female. All breeding and experimental procedures were undertaken in accordance with the Association for Research for Vision and Ophthalmology Statement for the Use of Animals in Ophthalmic and Vision Research. The Institutional Biosafety Committee (IBC) and the Animal Care and Use Committee (ACUC) at The Jackson Laboratory approved this study. The DBA/2 J and DBA/2 J-GpnmbR150X (D2-Gpnmb+) strains were utilized and have been described in detail elsewhere [21]. Mice were used at 9 months of age when the majority of eyes have had ongoing IOP elevation but detectable neurodegeneration has yet to occur [2, 21, 35]. In DBA/2 J mice, mutations in two genes (GpnmbR150X and Tyrp1b) drive an iris disease with features of human iris atrophy and pigment dispersion. In this disease, pigment disperses from the iris and induces damage in the drainage structures of the eye. This inhibits aqueous humour outflow and leads to an increase in intraocular pressure [36]. We used D2-Gpnmb+ mice as a control, a non-glaucomatous substrain of DBA/2 J. For radiation treated DBA/2 J mice, mice were placed on a rotating platform and a sub-lethal dose of γ-radiation (7.5 Gy; D2-RAD) was administered using a 137Cesium source in a single dose at 10 weeks of age. Our previous data has demonstrated that this level of treatment does not cause any adverse conditions and does not require bone marrow reconstitution [3, 37]. The optic nerves of all mice used in this study were confirmed to have no detectable nerve damage or axon loss as assessed by PPD staining (data not shown).

FAC sorting

FAC sorting of cells from the optic nerve head was performed as previously described [37]. Prior to cell collection, all surfaces and volumes were cleaned with 70% ethanol and RNaseZap (ThermoFisher Scientific) solution followed by dH20. Mice were euthanized by cervical dislocation, eyes enucleated, and placed immediately into ice-cold HBSS. For single ONH isolation, eyes were enucleated from the globe with curved scissors and the eye placed directly into ice cold HBSS (Gibco). The eye was punctured at the limbus with a 28 G needle and the cornea, iris, and lens removed to leave a posterior eye cup. A single cut was made through the choroid to the tissue surrounding the ONH, and the ONH (including the central retina and some extraocular tissue and choroid) was removed free from the globe. The ONH, representing the transparent pre-myelin transition zone segment of the optic nerve (~ 500 μm in length), was separated from the central retina and myelinated optic nerve and any additional pigmented tissue/extraocular tissue was removed. Single ONHs were placed directly into 100 μl of a custom HBSS, dispase (5 U/ml) (Stemcell Technologies), DNase I (2000 U/ml) (Worthington Biochemical), and SUPERase (1 U/μl) (ThermoFisher Scientific) solution. Samples were incubated for 20 mins at 37 °C and shaken at 350 RPM in an Eppendorf Thermomixer R followed by gentle trituration using a 200 μl pipette. Samples were blocked in 2% BSA, SUPERase (1 U/μl) in HBSS, and stained with secondary conjugated antibodies against CD11b, CD11c, CD34, CD45.2, GFAP (detailed in Table 1), as well as DAPI. This cocktail allowed other cell types to be accurately removed during FACS. FACS was performed on a FACSAria (BD Biosciences). Singlet gating of SSC-H vs SSC-W and FSC-H vs FSC-W gating plots was used to identify single cells. FSC-H vs DAPI was used to identify viable cells and CD11b+/CD45.2lo (and negative for all other markers; Fig. 1a) microglia were sorted into 100 μl buffer RLT + 1% β-ME, vortexed and frozen at − 80 °C until further processing. We have previously performed RNA-sequencing on mRNA of infiltrating monocytes (CD45hi/CD11b+/CD11c+) from the optic nerve head of 9 month DBA/2 J mice [37]. This marker panel was based on our previous findings that identified the majority of infiltrating immune cells as CD11c+ [3]. In these previous flow cytometry experiments < 3% of all myeloid derived cells (CD45+/CD11b+) in 9 month of age DBA/2 J optic nerve head tissue were resident CD11c+ microglia (CD45lo/CD11b+) [3, 37], and thus, CD11c+ microglia make up a negligible proportion of optic nerve head immune cells at this time point. In the current study, we aimed to enrich for resident microglia, as opposed to microglia that may be derived from infiltrating immune cells, and enriched for CD11c microglia for this purpose.

Table 1 Flow cytometry antibody details
Fig. 1
figure1

RNA-sequencing of optic nerve head microglia from D2 and D2-Gpnmb+mice. Microglia were FAC sorted from freshly isolated optic nerve heads and identified by being CD45lo/CD11b+ (a, see also Methods). Following RNA-sequencing of microglia from D2 and D2-Gpnmb+ optic nerve heads, samples were grouped by unsupervised hierarchical clustering (b; blue = strong correlation, red = weak correlation), creating D2 and D2-Gpnmb+ clusters (* denotes outlier excluded from subsequent analysis). c Genes were binned by log2 fold change (bin width 0.2) and coloured to show DE genes (red; q < 0.05). A simple summary is shown in the inset of c. d Scatter plot of all genes by mean log2 counts per million (CPM) for D2 (y) against D2-Gpnmb+ (x) showing DE genes (pink, q < 0.05; red, q < 0.001; non DE genes in grey) with top 20 DE genes annotated. Pathway analysis of DE genes revealed dysregulation of pathways involved in the microglial sensome and metabolism (e and f. Wikipathway analysis (e) showed significant dysregulation of 2 inflammatory pathways (q < 0.05). The number of DE genes within each pathway is shown . Ingenuity pathway analysis (f) also showed dysregulation of a number of sensome and inflammatory pathways, and metabolism pathways. Top 20 pathways sorted by z-score are shown, with the threshold for significant activation or inhibition marked. We queried metabolism dysregulation (g) demonstrating a trend towards an increased ratio of mt-RNA:nu-RNA in D2 microglia (P = 0.23) and DE gene expression (red, q < 0.05) in mitochondrial transcripts from nu-RNA and mt-RNA. DE of genes involved in OXPHOS and glycolysis/gluconeogenesis (GNG). In fPPRs = pattern recognition receptors

RNA-sequencing and analysis

Microglia samples were defrosted on ice and homogenized by syringe in RLT Buffer (total volume 300 μl). Total RNA was isolated using RNeasy micro kits as according to manufacturer’s protocols (Qiagen) including the optional DNase treatment step, and quality was assessed using an Agilent 2100 Bioanalyzer. The concentration was determined using a Ribogreen Assay from Invitrogen. Amplified dscDNA libraries were created using a Nugen Ovation RNA-seq System V2 and a primer titration was performed to remove primer dimers from the sample to allow sample inputs as low as 50 pg RNA. The SPIA dscDNA was sheared to 300 bp in length using a Diogenode Disruptor. Quality control was performed using an Agilent 2100 Bioanalyzer and a DNA 1000 chip assay. Library size produced was analysed using qPCR using the Library Quantitation kit/Illumina GA /ABI Prism (Kapa Biosystems). Libraries were barcoded, pooled, and sequenced 6 samples per lane on a HiSeq 2000 sequencer (Illumina) giving a depth of 30–35 million reads per sample.

Following RNA-sequencing samples were subjected to quality control analysis by a custom quality control python script. Reads with 70% of their bases having a base quality score ≥ 30 were retained for further analysis. Read alignment was performed using TopHat v 2.0.7 and expression estimation was performed using HTSeq with supplied annotations and default parameters against the DBA/2 J mouse genome (build-mm10). Bamtools v 1.0.2 were used to calculate the mapping statistics. Differential gene expression analysis between groups was performed using edgeR v 3.10.5 [38, 39] and the removal of outlier samples and lowly expressed genes was achieved by removing genes at a pre-defined cut-off level. We only included genes that were expressed at > 1 counts per million (CPM) in ≥4 samples across all samples for D2-Gpnmb+ to DBA/2 J (D2) comparison and ≥ 3 samples across all samples for D2 to D2-RAD comparison (chosen based on the size of the smallest group). We used unsupervised hierarchical clustering (HC) to generate clusters of samples with distinct gene expression profiles in which as many control samples were represented in a single cluster. HC was performed in R (1-cor, Spearman’s rho) based on a matrix of all samples, representing all genes post cut-off. Clusters were required to have ≥3 samples in order to compare using statistical testing; clusters that did not meet these criteria were removed from the analysis as outliers. (Although D2 sample D2_S1 was more distant than others within its cluster, analysing the data without this sample made no changes to the outcomes, data not shown.) Adjustment for multiple testing was performed using false discovery rate. DEseq2 and limma were used for differential expression analysis [40, 41]. Genes were considered to be significantly differentially expressed at a false discovery rate (FDR; Benjamini and Hochberg adjusted p values; q) of q < 0.05. Pathway analyses (see Results) were performed in R, WebGestalt (www.webgestalt.org; provides continuously curated, publicly availably pathways for exploring pathway enrichment) [42], and Ingenuity pathway analysis (IPA, Qiagen; which further explores directionality of pathway enrichment). Mitochondria gene lists were taken from published lists [43, 44]. Graphing was performed in R. Complete raw untrimmed count files can be found in Supplementary Data 1.

For comparisons to other published datasets, the microglia dataset generated here was compared to DBA/2 J whole optic nerve head at 8.5 months of age [2] (publically available from Datgan [45]) and DBA/2 J optic nerve head monocytes at 9 months of age [37] (Supplementary Data 2).

For the Supplementary Data 3 and 4, t refers to the t statistic, and B refers to log-odds that the gene is differentially expressed.

Results

Transcriptomic profiling of optic nerve head microglia

We performed RNA-sequencing on CD45lo/CD11b+/CD11c microglia samples from DBA/2 J (D2), control D2-Gpnmb+, and radiation treated DBA/2 J (D2-RAD) mice isolated from single optic nerve heads (ONHs) and enriched through FAC sorting (see Methods, and Fig. 1a). A total of 15 samples were amplified and sequenced (n = 5 for all groups; cell inputs = 165 ± 37 for D2, 182 ± 145 for D2-Gpnmb+, and 162 ± 173 for D2-RAD). To confirm the isolated cells were microglia, genes known to be highly expressed in microglia, astrocytes, neurons, oligodendrocytes, and infiltrating monocytes were analysed (Table 2). Genes associated with microglia were highly expressed in these samples. Low to no expression of other cell-type specific genes was observed, consistent with the samples primarily containing microglia.

Table 2 Cell type specific gene expression

Differences between D2 samples were expected due to the spontaneous nature of IOP insults, with some samples still resembling controls [35]. Hierarchical clustering (HC) [2, 28] demonstrated that the majority of D2 samples were more similar to each other than to normotensive D2-Gpnmb+ samples. HC generated two major clusters of samples containing: (1) four D2 samples (D2_S1–4), (2) five D2-Gpnmb+ samples (D2Gpnmb_S1–5) and one D2 sample (D2_S5) (Fig. 1b). The D2 sample in the second cluster (D2_S5) was removed for further analyses to create distinct disease (cluster 1; n = 4) and control (cluster 2; n = 5) groups. Between groups, 982 genes were differentially expressed (DE, Fig. 1c and d, Supplementary Data 3). Functional analysis of DE genes revealed significant changes in pathways associated with microglial function including surveillance, phagocytosis, metabolism, and inflammation (Fig. 1e and f), which are further expanded on below.

Genes that regulate microglial surveillance and phagocytosis are downregulated by chronic IOP elevation

A majority of the DE genes were downregulated in microglia from D2 eyes (Fig. 1c). These downregulated genes contributed to the enrichment of the TYROBP causal network and microglial pathogen phagocytosis pathways (Fig. 1e). The TYROBP causal network links the activation of TREM2 receptors with gene expression and thus contributes to the overall state of microglial activation. The downregulation of numerous genes in the network is consistent with a decrease in TYROBP signalling [46]. Ingenuity Pathway Analysis (IPA) identified another 19 pathways predicted to be significantly inhibited or less active in D2 microglia compared to D2-Gpnmb+ microglia (Fig. 1f). These pathways span a wide range of functions including phagocytosis, cell movement and shape, receptor-mediated signalling, and inflammation. Common downregulated DE genes within these pathways included C1 complex encoding genes (C1qa, C1qb, C1qc), integrins (Itgam, Itgb2, Cd37), Ig superfamily (Trem1, Trem2, Il10ra, Il13ra), phagocytic components (Nckap1l) and toll-like receptor signalling (Tlr1, Tlr3, Tlr7, and toll-like receptor pro-inflammatory enhancer Themis2 [47]). These data predict that many functions of microglia that help maintain tissue homeostasis are affected, and potentially inhibited, by chronic ocular hypertension.

Metabolism-related transcripts are upregulated in microglia by chronic IOP elevation

Changes in genes with mitochondrial and metabolic functions were identified in both gene and pathway level analyses. Of the top 20 DE genes (sorted by q) in D2 microglia compared to D2-Gpnmb+, 10 were mitochondrial transcriptome derived (mt-RNA, Fig. 1d). Nuclear encoded mitochondrial transcripts (nu-RNA) also differed from controls (22 up, 29 down, Fig. 1g). Ingenuity pathway analysis predicted that these changes promote oxidative phosphorylation (OXPHOS) activity (Fig. 1f). 17.5% of OXPHOS genes were DE and 19/20 DE genes had higher expression in D2 mice (Fig. 1g). Changes in the ratio of mt-RNA:nu-RNA can indicate changes in intracellular signalling between the mitochondria and nucleus. This ratio showed high variation between D2 samples (0.25 ± 0.17) but the mean ratio was not significantly different from D2-Gpnmb+ samples (0.13 ± 0.02; P = 0.23, Student’s t-test, Fig. 1g).

To further understand the metabolic state of microglia, we considered additional changes in metabolic genes. Hif1a, a master regulator of glycolysis and cell stress responses [48], was upregulated consistent with stress or inflammation, and with previous findings in the inner retina during glaucoma [35, 49]. Glycolysis regulator Pfkfb2 [50] was upregulated, as well as other glycolysis genes (Gapdh, Pgam1, Pgk1, Pgm2l1) (Fig. 1g). Slc16a1 (MCT1) was also up-regulated suggesting increased lactate, pyruvate, or ketone bodies transport [51]. The transporter is bi-directional, and as such could reflect either an attempt to increase microglial energy sources, or to increase metabolic support to retinal ganglion cell axons in the optic nerve head. Taken together, these data suggest that optic nerve head microglia in D2s develop an increased capacity to metabolise energy from various sources.

Metabolic switching between oxidative phosphorylation and glycolysis occurs in disease when microglia transition between pro-inflammatory and anti-inflammatory states. This cellular transition is associated with gene expression changes induced by both ROS and cytokine signalling pathways. In optic nerve head microglia, no enrichment was observed in pro or anti-inflammatory pathways based on the number of DE genes (Fig. 1e). Based on the direction of expression changes, three pro-inflammatory pathways were predicted to be inhibited; neuroinflammation signalling, Th1 signalling, and inflammasome signalling (Fig. 1f). We also assessed changes in a list of 20 genes associated with canonical M1 and M2 inflammation phenotypes [52]. Only two genes were differentially expressed, Ccl5 (decreased) and Tgfb2 (increased). TGF-β is an important anti-inflammatory signal for microglia that regulates their morphology, proliferation, and survival [53]. TGF-β signalling has been implicated in both homeostatic and anti-inflammatory signalling in microglia. Overall microglia gene expression in response to chronic IOP elevation was not pro-inflammatory, although they are predicted to be primed metabolically to facilitate rapid changes in phenotype.

Microglia and infiltrating monocytes have distinct phenotypes in pre-degenerative tissue

We compared gene expression changes in the isolated D2 microglia to previously defined gene expression changes from whole ONH tissue [2, 3] and from infiltrating monocytes [37] both from D2 mice at similar ages (Fig. 2, Supplementary Data 2). Howell et al. previously identified 5 molecularly distinct stages of disease in the D2 whole optic nerve head using hierarchical clustering of RNA microarray data, where stages 1–3 show no morphologically detectable neurodegeneration. A further 3 early stages (between stage 1 and 2) were subsequently identified [2, 3]. Based on corresponding upregulation or downregulation of genes, there was partial overlap between expression changes in isolated D2 microglia and each staged group of ONH samples (Fig. 2a). The most overlap with the D2 ONH was at an ocular hypertensive and pre-degenerative stage, representing 0.7% of DE genes in the ONH and 8% in microglia. A similar comparison between infiltrating monocytes and the ONH showed greater overlap, with 4-fold more genes in common. Very few DE genes showed the same directional DE changes across all 3 groups (Fig. 2b). A comparison irrespective of directional regulation showed a larger number of common DE genes between either microglia or monocytes with whole optic nerve head tissue (Fig. 2c), possibly as cell-specific effects may be masked in whole tissue analysis. Taken together, these data suggest that these two myeloid-cell derived populations have largely distinct responses relative to whole tissue changes in the ONH after chronic ocular hypertension.

Fig. 2
figure2

Monocytes, not microglia, are likely drivers of inflammatory signatures in glaucoma gene expression datasets. Comparison between DE genes in microglia (Tribble; current dataset), monocytes (Williams [37];), and whole optic nerve head in the D2 (Howell [2, 3];). Howell et al. [2] previously identified 5 molecularly distinct stages of disease in the D2 whole optic nerve head, where stages 1–2 show no morphologically detectable neurodegeneration. Further study subdivided stage 1 into 3 stages (1a, 1b, and 1c as shown between stages 1 and 2 here) [3]. a Euler diagrams show total number of DE genes within each dataset (number within each circle), the number of shared DE genes (with matching upregulation or downregulation, shown outside the corresponding intersection, and centrally (in red text) for matching genes for all 3 datasets (also displayed in b). c When direction of change is not considered, there are a greater number of shared DE genes for both microglia and monocytes, but a greater percentage of these are of matching upregulation or downregulation in monocytes. This may reflect strength of contribution of monocytes to whole optic nerve head changes over microglia. Further comparisons of neuroinflammatory genes (based on IPA gene sets) or mitochondria-related transcripts (from MitoCarta [43, 44]) demonstrate a number of genes that are uniquely DE in either microglia or monocytes (d; red = DE in microglia, blue = DE in monocytes; log2 CPM from D2 microglia and D2 infiltrating monocytes). Of the few shared DE genes (d; purple) the majority of neuroinflammatory and nu-genes show the same directional changes but the opposite was true for mt derived genes (e; upregulation = red, downregulation = blue)

To further elucidate the neuroinflammatory and metabolic phenotypes, we compared neuroinflammatory and mitochondria gene expression in microglia and monocytes (Fig. 2d and e). Of 144 neuroinflammatory genes annotated by Ingenuity, few were DE in either microglia (n = 18) or monocytes (n = 16) at this pre-degenerative stage of disease (Fig. 2d). The shared neuroinflammatory DE genes (n = 5) showed the same direction of change in both cell types (Fig. 2e). These five genes (Cx3cr1, Tgfb2, Tgfbr1, Trem2, and P2rx7) offer excellent candidates for genetic manipulation to test function of innate immune pathways in the optic nerve head.

Metabolic changes are a prominent feature in cells and tissue affected by elevated IOP based on RNA-sequencing datasets [35, 37, 54]. Changes to transcripts encoding mitochondrial proteins feature in both the microglia and monocyte RNA-sequencing datasets (Fig. 2d and e). Dysregulation of nu-RNA was greater in monocytes (n = 121 genes) than microglia (n = 48 genes). Of the 10 shared DE nu-genes, 7 showed the same directional change (Fig. 2e) but with a greater magnitude of dysregulation in monocytes. For mt-RNA transcripts 12/15 genes were DE in both monocytes and microglia, but with markedly different expression profiles (0/12 being co-up- or co-down- regulated; Fig. 2e). Thus, our data suggest a pro-metabolic status in microglia that is not matched in infiltrating monocytes. Together our data implicate mitochondrial / metabolic changes in optic nerve head immune cells as an early disease feature.

Pre-treatment by irradiation reduces the effects of ocular hypertension on microglia

Low dose irradiation of mice at a young age prevents glaucomatous neurodegeneration in D2 mice without lowering IOP [3]. Altered microglia have been suggested to contribute to the protective effects of radiation [34]. We compared and analysed D2 against RAD-D2 samples. HC generated 2 clusters representing 1) three D2 samples (D2_S1, 2, and 4), and 2) five RAD-D2 samples (D2-RAD_S1–5) and 1 D2 sample (D2_S3) (Fig. 3a). This single D2 sample from cluster 2 was removed by the dendogram cut (Fig. 3a). RAD-D2 microglia exhibited 2246 DE genes compared to D2 microglia (783 upregulated, 1463 downregulated) (Fig. 3b and c, Supplementary Data 4). Pathway analyses (Fig. 3d-f) showed that radiation treatment affects phagocytosis, metabolism, mitochondria, and inflammation related genes, all pathways altered in microglia (see above). There were 579 genes altered by radiation that overlapped with changes in glaucoma (D2 vs. D2-Gpnmb+comparison). For 578 of these genes, radiation treatment corrected the disease-related change (Fig. 3f). Dcbld2 was the single gene DE downregulated in both datasets. Dcbld2 encodes the endothelial and smooth muscle cell-derived neuropilin-like protein (ESDN) which is upregulated in endothelial cells following vascular injury [55]. Its deletion or downregulation impairs retinal angiogenesis [56] and promotes insulin signalling [57]. Its expression is not limited to endothelial cells; with relatively robust expression in microglia, astrocytes, and neurons (brain RNA-sequencing [58]).

Fig. 3
figure3

RNA-sequencing of optic nerve head microglia from D2 and radiation therapy treated D2 mice. Following RNA-sequencing of microglia form D2 and radiation therapy treated (RAD-D2) optic nerve heads, samples were grouped by unsupervised hierarchical clustering (a; blue = strong correlation, red = weak correlation), creating D2 and RAD-D2 clusters (* denotes outlier excluded from subsequent analysis). b Genes were binned by log2 fold change (bin width 0.2) and coloured to show DE genes (red; q < 0.05). A simple summary is shown in the inset of b. (c) Scatter plot of all genes by mean log2 counts per million (CPM) for D2 (x) against RAD-D2 (y), showing DE genes (pink, q < 0.05; red, q < 0.001; non DE genes in grey) with top 20 DE genes annotated. Ingenuity pathway analysis (d) showed dysregulation of a number of sensome and inflammatory pathways, and metabolism pathways, representing a correction of the D2 phenotype in RAD treated animals. Top 20 pathways sorted by z-score are shown, with the threshold for significant activation or inhibition marked. Relevant to metabolism and signaling (e) mt-RNA:nu-RNA ratios are shown. DE gene expression (red, q < 0.05) in mitochondrial transcripts in nu-RNA and mt-RNA are shown. Genes involved in OXPHOS and glycolysis/gluconeogenesis (GNG) are compared. f Overlap of DE genes from D2 vs. D2-Gpnmb+ and D2 vs. RAD-D2 datasets. In ePPRs = pattern recognition receptors

Discussion

Neuroinflammation at the site of the optic nerve head (ONH) may be a critical pathogenic event in glaucoma. Our group and others have identified the ONH as a likely candidate for the initial site of damage in glaucoma across species [22, 59,60,61]. In this manuscript we have identified further transcriptomic changes at the level of ONH microglia. The current study focuses on CD11c microglia, the most prominent microglia subtype in our previous datasets (representing > 97% of microglia identified in the ONH, [3, 37]). We have previously used CD11c as a marker to distinguish myeloid-derived cell subtypes within the ONH; with the majority of infiltrating monocyte-like cells in the ONH being CD45hi and CD11c+. Emerging evidence is demonstrating a role for CD11c+ microglia in neurodegenerative disease progression, especially during demyelinating events [62], and with T-cell interactions in the brain [63]. Given that the ONH is unmyelinated, and that T-cell changes have not been found in the glaucomatous ONH, it is unsurprising to observe so few CD11c+ microglia. Here, we focused on Cd11c microglia, a sub-type more associated with tissue surveillance and inflammation. These microglia were affected by chronic elevated intraocular pressure (IOP) based on changes at a transcriptional level, consistent with previous studies showing microglial activation in glaucoma.

To explore molecular changes that lead to optic nerve degeneration we have previously performed microarray gene expression profiling of the whole ONH [2, 3]. This data set shows changes to inflammatory pathways, but the attributive cells were unknown. To further understand the molecular events that happen in the ONH at a cell-type level we have performed transcriptomic profiling of monocyte-like cells in the same model of glaucoma [37]. These cells were highly pro-inflammatory and express various complement genes and integrins. Targeting the α subunit of complement receptor 3 (genetic ablation of Itgam encoding CD11b) prevents monocyte-like cell extravasation into the ONH and significantly reduced the risk of developing severe glaucomatous neurodegeneration. CD11b is well expressed on microglia [64] and we used cell-surface expression of CD11b to enrich for microglia. In the data presented here we predict that DBA/2 J microglia are initially anti-inflammatory. Itgam is downregulated in DBA/2 J microglia (in this data set) following periods of ocular hypertension. Nevertheless, the protection that results from removing CD11b (Itgam knockout [37];) could, in part, be due to its effects on microglia, but elucidating exactly which cell-type is at play will take definitive testing using cell-type-specific cre-lines. This dichotomy of microglia and monocytes in the ONH offers an exciting avenue for further research, particularly during the stage of monocyte extravasation into the ONH. For example, microglia processes have been shown to reseal injured capillaries following recruitment via P2ry12 mediated chemotaxis [65]. P2ry12 is downregulated in D2 microglia in the present study, which could suggest defects in initial-protective microglial responses that exacerbate monocyte initiated injury.

To date, one of the most protective therapies in DBA/2 J glaucoma has been radiation therapy [3, 6, 33, 34]. A sub-lethal dose of radiation (γ- or X- ray) early in life changes the neuroinflammatory response to ocular hypertension later in life. One reason that radiation treatment is protective is that it increases the expression of GlyCAM1 in the retina and ONH, thereby reducing the entry of pro-inflammatory monocytes [6]. Another type of neuroinflammatory cell-type affected by radiation therapy is microglia [34]. Microglia respond acutely to radiation therapy; however, the long-term effects on microglia need to be more fully defined. As presented here, the majority of genes (59%) dysregulated by elevated IOP in untreated DBA/2 J microglia had a different level of expression in D2-RAD microglia. In addition, microglia in treated mice had a different predicted activation state of TREM1 signalling and more broadly alterations of neuroinflammation signalling and metabolic pathways related to energy production (Fig. 3d). It is not clear yet whether this is a direct effect of the radiation treatment within microglia or a secondary effect of microglial interactions with other cells. TREM1 signalling may be a critical factor in controlling microglial activation in glaucoma and needs to be tested further. Determining the function of TREM1 signalling in glaucoma is likely to uncover more specific mechanisms of damage by elevated IOP and protection by radiation therapy.

Phagocytosis is suggested to be a critical process in the ONH for healthy tissue maintenance [66], but its role in ONH microglia has not been explored. Microglial phagocytosis is generally regulated by the TREM2-TYROBP signalling pathway [67], and the TYROBP signalling network in microglia is among the earliest affected by chronic elevated IOP. Mutations in TYROBP cause Nasu-Hakola disease [68], in which the neurodegenerative pathology is defined as a primary microglial disorder [69]. This underlies the important function of this pathway and microglia toward directly preventing disease. Disruption of the TYROBP signalling network has also been demonstrated in Alzheimer’s disease [46]. Furthermore, the increased risk of Alzheimer’s disease caused by TREM2 mutations [70] has led to research showing many effects of TREM2 on phagocytosis, transcription, metabolism, and inflammation [71,72,73]. These data suggest that this signalling pathway may be a master regulator of other changes observed here. Manipulating TYROBP and TREM2 is a promising strategy to define functions of microglia in glaucoma and a possible avenue for treatment.

Retinal ganglion cell axons remain unmyelinated at the ONH and therefore are particularly vulnerable to glaucoma related stresses (age and elevated intraocular pressure) in addition to metabolic strain. Recently we have demonstrated that a critical metabolic vulnerability exists in retinal ganglion cells in DBA/2 J glaucoma with marked mitochondrial and metabolic changes [35]. Preventing these metabolic events (either systemically or specific to retinal ganglion cells) robustly protects from glaucomatous neurodegeneration [25, 54] although neuroinflammatory features still remain [25]. The neuro-glia-vascular complexes of the ONH form a contained metabolic unit in which glia provide trophic and metabolic support to retinal ganglion cell axons [74]. In our data, microglia are predicted to become more metabolically active, a metabolic shift that is typically anti-inflammatory and pro-supportive to neurons [75]. In this sense, microglia could be early mediators of neuroprotection in glaucoma. More work is needed to explore microglial function at different stages of disease, along with a consideration of microglial subtypes. One subtype, CD11c+ microglia, represent a small proportion of ONH microglia in our early-stage data, but may become important at later stages or in other regions of affected tissue. It is also possible that microglia later become pro-inflammatory. Thus further experiments are needed to resolve how elevated IOP triggers changes in microglia and if polarization of microglia continues to change during disease progression.

Metabolic treatments that protect retinal ganglion cells do not prevent all neuroinflammatory events in the retina and ONH [25, 35]. The long-term effects of this neuroinflammation on survival or function of the optic nerve are not known. Thus, inflammation remains an important target to consider for new therapies. At this pre-degenerative stage of glaucoma, we observed few changes in inflammatory molecules in microglia, and attribute many of the neuroinflammatory changes to infiltrating monocyte-like cells. Recent evidence using single cell RNA-sequencing has demonstrated that microglia display a heterogeneous repertoire of inflammatory responses in diseased tissue. Our study may lack the cellular resolution to detect microglial heterogeneity associated with inflammation, if present. However, reactive microglia in the ONH have been identified immunohistochemically in glaucoma [5, 7] and more precisely characterizing the timing and cellular heterogeneity of such changes will provide deeper insights into glaucoma pathophysiology.

Metabolic regulation is an important aspect of myeloid-derived cell polarity and function. Resident microglia can exhibit similar M1/M2 phenotypes to peripheral macrophages, the former representing an activated, pro-inflammatory phenotype and the latter a resting, anti-inflammatory phenotype [75]. An M1 state is consistent with aerobic glycolysis, where metabolic resources can be directed towards cell proliferation and activation, and where the ROS generated for phagocytic clearance will not interfere with and uncouple electron transport [76]. An M2 polarisation state is consistent with increased glucose metabolism and mitochondrial biogenesis [77]. Microglia in our dataset did not conform to either phenotype; we observed increased expression in glycolysis genes (including Hif1a), but also an upregulation of OXPHOS genes. In fact, the microglia exhibited increased metabolic upregulation from multiple energy sources. RNA sequencing of microglia in neurodegenerative disease has shown that microglial exhibit more nuanced states than these simple M1 and M2 polar opposites [14, 72, 78,79,80]. There is also evidence to suggest that dysregulation of metabolism and accumulation of mt-DNA mutations within microglia may itself be a trigger for microglial activation in neurodegenerative disease [81]. Dysregulated metabolism in microglia may also impair their ability to mitigate neuroinflammation. Trem2−/− 5XFAD mice (an Alzheimer’s disease model) show impaired glycolysis, reduced ATP production, and increased autophagy in microglia. This disruption to metabolism may contribute to the reduced microglial phagocytosis observed in Trem2−/− 5XFAD mice. Microglia in aged mice also display upregulated oxidative phosphorylation and utilisation of ketone energy sources, which may represent a stress response or loss of transcriptional regulation [82]. The dysregulation of metabolism seen in D2 microglia may be indicative of metabolic stress which may perturb microglial responses to neuroinflammation. Therapeutic strategies that target microglia metabolism may be valid targets for glaucoma treatment.

Conclusions

This study identifies alterations in mitochondrial gene expression as well as changes to phagocytosis, inflammatory, and sensome pathways in microglia prior to optic nerve degeneration in DBA/2 J glaucoma. Out of this complex set of changes, the TREM/TYROBP signalling network emerged as a potential master regulator of microglial responses to elevated IOP. This type of response is observed in other neurodegenerative diseases where homeostatic functions of microglia are suppressed by inhibition of TYROBP signalling. This occurred in the absence of a clear pro-inflammatory response by microglia, suggesting that early pro-inflammatory signals may largely derive from infiltrating monocyte-like cells in this model of glaucoma. These data shed new light on the early myeloid changes in the ONH following periods of ocular hypertension and offer novel targets for treatment and further mechanistic exploration.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its supplementary information files].

Abbreviations

ACUC:

Animal Care and Use Committee

D2:

DBA/2 J mouse

DE:

Differentially expressed

FACS:

Fluorescent activated cell sorting

FDR:

False discovery rate

HC:

Hierarchical clustering

IBC:

Institutional Biosafety Committee

IOP:

Intraocular pressure

ONH:

Optic nerve head

PPD:

Paraphenylenediamine

RAD:

Radiation treated / irradiation therapy

References

  1. 1.

    Tham YC, Li X, Wong TY, Quigley HA, Aung T, Cheng CY. Global prevalence of glaucoma and projections of glaucoma burden through 2040: a systematic review and meta-analysis. Ophthalmology. 2014;121(11):2081–90.

    PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Howell GR, Macalinao DG, Sousa GL, Walden M, Soto I, Kneeland SC, et al. Molecular clustering identifies complement and endothelin induction as early events in a mouse model of glaucoma. J Clin Investig. 2011;121(4):1429–44.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Howell GR, Soto I, Zhu X, Ryan M, Macalinao DG, Sousa GL, et al. Radiation treatment inhibits monocyte entry into the optic nerve head and prevents neuronal damage in a mouse model of glaucoma. J Clin Investig. 2012;122(4):1246–61.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Williams PA, Tribble JR, Pepper KW, Cross SD, Morgan BP, Morgan JE, et al. Inhibition of the classical pathway of the complement cascade prevents early dendritic and synaptic degeneration in glaucoma. Mol Neurodegener. 2016;11:26.

  5. 5.

    Bosco A, Steele MR, Vetter ML. Early microglia activation in a mouse model of chronic Glaucoma. J Comp Neurol. 2011;519(4):599–620.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Williams PA, Braine CE, Foxworth NE, Cochran KE, John SWM. GlyCAM1 negatively regulates monocyte entry into the optic nerve head and contributes to radiation-based protection in glaucoma. J Neuroinflammation. 2017;14(1):93.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. 7.

    Neufeld AH. Microglia in the optic nerve head and the region of parapapillary chorioretinal atrophy in glaucoma. Arch Ophthalmol. 1999;117(8):1050–6.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Tsilis AG, Tsilidis KK, Pelidou SH, Kitsos G. Systematic review of the association between Alzheimer's disease and chronic glaucoma. Clin Ophthalmol. 2014;8:2095–104.

    PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Mancino R, Martucci A, Cesareo M, Giannini C, Corasaniti MT, Bagetta G, et al. Glaucoma and Alzheimer disease: one age-related neurodegenerative disease of the brain. Curr Neuropharmacol. 2018;16(7):971–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Colligris P, Perez de Lara MJ, Colligris B, Pintor J. Ocular Manifestations of Alzheimer's and Other Neurodegenerative Diseases: The Prospect of the Eye as a Tool for the Early Diagnosis of Alzheimer's Disease. J Ophthalmol. 2018;2018:8538573.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. 11.

    Fernández-Albarral JA, Salobrar-García E, Martínez-Páramo R, Ramírez AI, de Hoz R, Ramírez JM, et al. Retinal glial changes in Alzheimer's disease - a review. J Optom. 2018;12(3):198-207.

    PubMed  Article  Google Scholar 

  12. 12.

    Williams PA, Marsh-Armstrong N, Howell GR, Lasker IIA. Neuroinflammation in glaucoma: a new opportunity. Exp Eye Res. 2017;157:20–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Hickman S, Izzy S, Sen P, Morsett L, El Khoury J. Microglia in neurodegeneration. Nat Neurosci. 2018;21(10):1359–69.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Bachiller S, Jiménez-Ferrer I, Paulus A, Yang Y, Swanberg M, Deierborg T, et al. Microglia in neurological diseases: a road map to brain-disease dependent-inflammatory response. Front Cell Neurosci. 2018;12:488.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Song WM, Colonna M. The identity and function of microglia in neurodegeneration. Nat Immunol. 2018;19(10):1048–58.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Deczkowska A, Keren-Shaul H, Weiner A, Colonna M, Schwartz M, Amit I. Disease-associated microglia: a universal immune sensor of Neurodegeneration. Cell. 2018;173(5):1073–81.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Haruwaka K, Ikegami A, Tachibana Y, Ohno N, Konishi H, Hashimoto A, et al. Dual microglia effects on blood brain barrier permeability induced by systemic inflammation. Nat Commun. 2019;10(1):5816.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Mukherjee S, Klaus C, Pricop-Jeckstadt M, Miller JA, Struebing FL. A microglial signature directing human aging and Neurodegeneration-related gene networks. Front Neurosci. 2019;13:2.

    PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Bonham LW, Sirkis DW, Yokoyama JS. The transcriptional landscape of microglial genes in aging and neurodegenerative disease. Front Immunol. 2019;10:1170.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Hansen DV, Hanson JE, Sheng M. Microglia in Alzheimer's disease. J Cell Biol. 2018;217(2):459–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Libby RT, Anderson MG, Pang IH, Robinson ZH, Savinova OV, Cosma IM, et al. Inherited glaucoma in DBA/2J mice: pertinent disease features for studying the neurodegeneration. Vis Neurosci. 2005;22(5):637–48.

    PubMed  Article  Google Scholar 

  22. 22.

    Howell GR, Libby RT, Jakobs TC, Smith RS, Phalan FC, Barter JW, et al. Axons of retinal ganglion cells are insulted in the optic nerve early in DBA/2J glaucoma. J Cell Biol. 2007;179(7):1523–37.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Ren L, Danias J. A role for complement in glaucoma? Adv Exp Med Biol. 2010;703:95–104.

    PubMed  Article  Google Scholar 

  24. 24.

    Stasi K, Nagel D, Yang XY, Wang RF, Ren LZ, Podos SM, et al. Complement component 1Q (C1Q) upregulation in retina of murine, primate, and human glaucomatous eyes. Invest Ophthalmol Vis Sci. 2006;47(3):1024–9.

    PubMed  Article  Google Scholar 

  25. 25.

    Harder JM, Braine CE, Williams PA, Zhu X, MacNicoll KH, Sousa GL, et al. Early immune responses are independent of RGC dysfunction in glaucoma with complement component C3 being protective. Proc Natl Acad Sci U S A. 2017;114(19):E3839–E48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Bosco A, Anderson SR, Breen KT, Romero CO, Steele MR, Chiodo VA, et al. Complement C3-targeted gene therapy restricts onset and progression of Neurodegeneration in chronic mouse Glaucoma. Mol Ther. 2018;26(10):2379–96.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Wang L, Fortune B, Cull G, Dong J, Cioffi GA. Endothelin B receptor in human glaucoma and experimentally induced optic nerve damage. Arch Ophthalmol. 2006;124(5):717–24.

    CAS  PubMed  Article  Google Scholar 

  28. 28.

    Minton AZ, Phatak NR, Stankowska DL, He S, Ma HY, Mueller BH, et al. Endothelin B receptors contribute to retinal ganglion cell loss in a rat model of glaucoma. PLoS One. 2012;7(8):e43199.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Howell GR, MacNicoll KH, Braine CE, Soto I, Macalinao DG, Sousa GL, et al. Combinatorial targeting of early pathways profoundly inhibits neurodegeneration in a mouse model of glaucoma. Neurobiol Dis. 2014;71:44–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Johnson EC, Jia L, Cepurna WO, Doser TA, Morrison JC. Global changes in optic nerve head gene expression after exposure to elevated intraocular pressure in a rat glaucoma model. Invest Ophthalmol Vis Sci. 2007;48(7):3161–77.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Yasuda M, Tanaka Y, Ryu M, Tsuda S, Nakazawa T. RNA sequence reveals mouse retinal transcriptome changes early after axonal injury. PLoS One. 2014;9(3):e93258.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  32. 32.

    Yasuda M, Tanaka Y, Omodaka K, Nishiguchi KM, Nakamura O, Tsuda S, et al. Transcriptome profiling of the rat retina after optic nerve transection. Sci Rep. 2016;6:28736.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Anderson MG, Libby RT, Gould DB, Smith RS, John SWM. High-dose radiation with bone marrow transfer prevents neurodegeneration in an inherited glaucoma. Proc Natl Acad Sci U S A. 2005;102(12):4566–71.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Bosco A, Crish SD, Steele MR, Romero CO, Inman DM, Horner PJ, et al. Early reduction of microglia activation by irradiation in a model of chronic glaucoma. PLoS One. 2012;7(8):e43602.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Williams PA, Harder JM, Foxworth NE, Cochran KE, Philip VM, Porciatti V, et al. Vitamin B-3 modulates mitochondrial vulnerability and prevents glaucoma in aged mice. Science. 2017;355(6326):756–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Howell GR, Libby RT, Marchant JK, Wilson LA, Cosma IM, Smith RS, et al. Absence of glaucoma in DBA/2J mice homozygous for wild-type versions of Gpnmb and Tyrp1. BMC Genet. 2007;8:45.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  37. 37.

    Williams PA, Braine CE, Kizhatil K, Foxworth NE, Tolman NG, Harder JM, et al. Inhibition of monocyte-like cell extravasation protects from neurodegeneration in DBA/2J glaucoma. Mol Neurodegener. 2019;14(1):6.

    PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–97.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  41. 41.

    Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  42. 42.

    Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: a more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Res. 2017;45(W1):W130–W7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Calvo SE, Clauser KR, Mootha VK. MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins. Nucleic Acids Res. 2016;44(D1):D1251–7.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Pagliarini DJ, Calvo SE, Chang B, Sheth SA, Vafai SB, Ong SE, et al. A mitochondrial protein compendium elucidates complex I disease biology. Cell. 2008;134(1):112–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Howell GR, Walton DO, King BL, Libby RT, John SW. Datgan, a reusable software system for facile interrogation and visualization of complex transcription profiling data. BMC Genomics. 2011;12:429.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Zhang B, Gaiteri C, Bodea LG, Wang Z, McElwee J, Podtelezhnikov AA, et al. Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease. Cell. 2013;153(3):707–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Peirce MJ, Brook M, Morrice N, Snelgrove R, Begum S, Lanfrancotti A, et al. Themis2/ICB1 is a signaling scaffold that selectively regulates macrophage toll-like receptor signaling and cytokine production. PLoS One. 2010;5(7):e11465.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  48. 48.

    Wang GL, Jiang BH, Rue EA, Semenza GL. Hypoxia-inducible factor 1 is a basic-helix-loop-helix-PAS heterodimer regulated by cellular O2 tension. Proc Natl Acad Sci U S A. 1995;92(12):5510–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Chidlow G, Wood JPM, Casson RJ. Investigations into hypoxia and oxidative stress at the optic nerve head in a rat model of Glaucoma. Front Neurosci. 2017;11:478.

  50. 50.

    Darville MI, Crepin KM, Hue L, Rousseau GG. 5′ flanking sequence and structure of a gene encoding rat 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase. Proc Natl Acad Sci U S A. 1989;86(17):6543–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Halestrap AP, Meredith D. The SLC16 gene family-from monocarboxylate transporters (MCTs) to aromatic amino acid transporters and beyond. Pflugers Arch. 2004;447(5):619–28.

    CAS  PubMed  Article  Google Scholar 

  52. 52.

    Lan X, Han X, Li Q, Yang QW, Wang J. Modulators of microglial activation and polarization after intracerebral haemorrhage. Nat Rev Neurol. 2017;13(7):420–33.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Ma W, Silverman SM, Zhao L, Villasmil R, Campos MM, Amaral J, et al. Absence of TGFβ signaling in retinal microglia induces retinal degeneration and exacerbates choroidal neovascularization. eLife. 2019;8:e42049.

    PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    Williams PA, Harder JM, Foxworth NE, Cardozo BH, Cochran KE, John SWM. Nicotinamide and WLDS act together to prevent Neurodegeneration in Glaucoma. Front Neurosci. 2017;11:232.

  55. 55.

    Kobuke K, Furukawa Y, Sugai M, Tanigaki K, Ohashi N, Matsumori A, et al. ESDN, a novel neuropilin-like membrane protein cloned from vascular cells with the longest secretory signal sequence among eukaryotes, is up-regulated after vascular injury. J Biol Chem. 2001;276(36):34105–14.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Nie L, Guo X, Esmailzadeh L, Zhang J, Asadi A, Collinge M, et al. Transmembrane protein ESDN promotes endothelial VEGF signaling and regulates angiogenesis. J Clin Invest. 2013;123(12):5082–97.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Li X, Jung JJ, Nie L, Razavian M, Zhang J, Samuel V, et al. The neuropilin-like protein ESDN regulates insulin signaling and sensitivity. Am J Physiol Heart Circ Physiol. 2016;310(9):H1184–93.

    PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Zhang Y, Chen K, Sloan SA, Bennett ML, Scholze AR, O'Keeffe S, et al. An RNA-sequencing transcriptome and splicing database of glia, neurons, and vascular cells of the cerebral cortex. J Neurosci. 2014;34(36):11929–47.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Johnson EC, Deppmeier LMH, Wentzien SKF, Hsu I, Morrison JC. Chronology of optic nerve head and retinal responses to elevated intraocular pressure. Invest Ophthalmol Vis Sci. 2000;41(2):431–42.

    CAS  PubMed  Google Scholar 

  60. 60.

    Burgoyne CF, Downs JC, Bellezza AJ, Suh JKF, Hart RT. The optic nerve head as a biomechanical structure: a new paradigm for understanding the role of IOP-related stress and strain in the pathophysiology of glaucomatous optic nerve head damage. Prog Retin Eye Res. 2005;24(1):39–73.

    PubMed  Article  Google Scholar 

  61. 61.

    Son JL, Soto I, Oglesby E, Lopez-Roca T, Pease ME, Quigley HA, et al. Glaucomatous optic nerve injury involves early astrocyte reactivity and late Oligodendrocyte loss. Glia. 2010;58(7):780–9.

    PubMed  Google Scholar 

  62. 62.

    Sato-Hashimoto M, Nozu T, Toriba R, Horikoshi A, Akaike M, Kawamoto K, et al. Microglial SIRPα regulates the emergence of CD11c+ microglia and demyelination damage in white matter. eLife. 2019;8:e42025.

    PubMed  PubMed Central  Article  Google Scholar 

  63. 63.

    Herz J, Johnson KR, McGavern DB. Therapeutic antiviral T cells noncytopathically clear persistently infected microglia after conversion into antigen-presenting cells. J Exp Med. 2015;212(8):1153–69.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. 64.

    Hickman SE, Kingery ND, Ohsumi TK, Borowsky ML, Wang LC, Means TK, et al. The microglial sensome revealed by direct RNA sequencing. Nat Neurosci. 2013;16(12):1896–905.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Lou N, Takano T, Pei Y, Xavier AL, Goldman SA, Nedergaard M. Purinergic receptor P2RY12-dependent microglial closure of the injured blood-brain barrier. Proc Natl Acad Sci U S A. 2016;113(4):1074–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Nguyen JV, Soto I, Kim KY, Bushong EA, Oglesby E, Valiente-Soriano FJ, et al. Myelination transition zone astrocytes are constitutively phagocytic and have synuclein dependent reactivity in glaucoma. Proc Natl Acad Sci U S A. 2011;108(3):1176–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  67. 67.

    Neumann H, Takahashi K. Essential role of the microglial triggering receptor expressed on myeloid cells-2 (TREM2) for central nervous tissue immune homeostasis. J Neuroimmunol. 2007;184(1–2):92–9.

    CAS  PubMed  Article  Google Scholar 

  68. 68.

    Klünemann HH, Ridha BH, Magy L, Wherrett JR, Hemelsoet DM, Keen RW, et al. The genetic causes of basal ganglia calcification, dementia, and bone cysts: DAP12 and TREM2. Neurology. 2005;64(9):1502–7.

    PubMed  Article  CAS  Google Scholar 

  69. 69.

    Thrash JC, Torbett BE, Carson MJ. Developmental regulation of TREM2 and DAP12 expression in the murine CNS: implications for Nasu-Hakola disease. Neurochem Res. 2009;34(1):38–45.

    CAS  PubMed  Article  Google Scholar 

  70. 70.

    Guerreiro R, Wojtas A, Bras J, Carrasquillo M, Rogaeva E, Majounie E, et al. TREM2 variants in Alzheimer's disease. N Engl J Med. 2013;368(2):117–27.

    CAS  PubMed  Article  Google Scholar 

  71. 71.

    Krasemann S, Madore C, Cialic R, Baufeld C, Calcagno N, El Fatimy R, et al. The TREM2-APOE Pathway Drives the Transcriptional Phenotype of Dysfunctional Microglia in Neurodegenerative Diseases. Immunity. 2017;47(3):566–81 e9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Keren-Shaul H, Spinrad A, Weiner A, Matcovitch-Natan O, Dvir-Szternfeld R, Ulland TK, et al. A Unique Microglia Type Associated with Restricting Development of Alzheimer's Disease. Cell. 2017;169(7):1276–90 e17.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  73. 73.

    Ulland TK, Song WM, Huang SC, Ulrich JD, Sergushichev A, Beatty WL, et al. TREM2 Maintains Microglial Metabolic Fitness in Alzheimer's Disease. Cell. 2017;170(4):649–63 e13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Sun D, Lye-Barthel M, Masland RH, Jakobs TC. The morphology and spatial arrangement of astrocytes in the optic nerve head of the mouse. J Comp Neurol. 2009;516(1):1–19.

    PubMed  PubMed Central  Article  Google Scholar 

  75. 75.

    Orihuela R, McPherson CA, Harry GJ. Microglial M1/M2 polarization and metabolic states. Br J Pharmacol. 2016;173(4):649–65.

    CAS  PubMed  Article  Google Scholar 

  76. 76.

    Brown GC. Regulation of mitochondrial respiration by nitric oxide inhibition of cytochrome c oxidase. Biochim Biophys Acta. 2001;1504(1):46–57.

    CAS  PubMed  Article  Google Scholar 

  77. 77.

    Vats D, Mukundan L, Odegaard JI, Zhang L, Smith KL, Morel CR, et al. Oxidative metabolism and PGC-1beta attenuate macrophage-mediated inflammation. Cell Metab. 2006;4(1):13–24.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  78. 78.

    De Biase LM, Schuebel KE, Fusfeld ZH, Jair K, Hawes IA, Cimbro R, et al. Local Cues Establish and Maintain Region-Specific Phenotypes of Basal Ganglia Microglia. Neuro. 2017;95(2):341–56 e6.

    Google Scholar 

  79. 79.

    Friedman BA, Srinivasan K, Ayalon G, Meilandt WJ, Lin H, Huntley MA, et al. Diverse brain myeloid expression profiles reveal distinct microglial activation states and aspects of Alzheimer's disease not evident in mouse models. Cell Rep. 2018;22(3):832–47.

    CAS  PubMed  Article  Google Scholar 

  80. 80.

    Tay TL, Sagar, Dautzenberg J, Grün D, Prinz M. Unique microglia recovery population revealed by single-cell RNAseq following neurodegeneration. Acta Neuropathol Commun. 2018;6(1):87.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  81. 81.

    Rodríguez MJ, Martínez-Moreno M, Ortega FJ, Mahy N. Targeting microglial K (ATP) channels to treat neurodegenerative diseases: a mitochondrial issue. Oxidative Med Cell Longev. 2013;2013:194546.

    Article  CAS  Google Scholar 

  82. 82.

    Flowers A, Bell-Temin H, Jalloh A, Stevens SM, Bickford PC. Proteomic anaysis of aged microglia: shifts in transcription, bioenergetics, and nutrient response. J Neuroinflammation. 2017;14(1):96.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references

Acknowledgments

The Authors would like to thank Flow Cytometry, Histology, and Computational Sciences Services at The Jackson Laboratory, Gareth Howell for critical reading of the manuscript, and Mimi de Vries and Amy Bell for organizing mouse colonies.

Funding

EY011721, the Barbra and Joseph Cohen Foundation, startup funds from Columbia University, and the Precision Medicine Initiative (SWMJ). Vetenskapsrådet 2018–02124 (PAW). Simon John is an Investigator of HHMI. Pete Williams is supported by the Karolinska Institutet in the form of a Board of Research Faculty Funded Career Position and by St. Erik Eye Hospital philanthropic donations. Open access funding provided by Karolinska Institutet.

Author information

Affiliations

Authors

Contributions

JRT – analysed data, wrote the manuscript; JMH – analysed data, wrote the manuscript; PAW – conceived, designed, performed, and analysed experiments, wrote the manuscript; SWMJ – conceived and oversaw the project, wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Pete A. Williams or Simon W. M. John.

Ethics declarations

Ethics approval

All breeding and experimental procedures were undertaken in accordance with the Association for Research for Vision and Ophthalmology Statement for the Use of Animals in Ophthalmic and Research. The Institutional Biosafety Committee (IBC) and the Animal Care and Use Committee (ACUC) at The Jackson Laboratory approved this study.

Consent for publication

N/A

Competing interests

The Authors report no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tribble, J.R., Harder, J.M., Williams, P.A. et al. Ocular hypertension suppresses homeostatic gene expression in optic nerve head microglia of DBA/2 J mice. Mol Brain 13, 81 (2020). https://doi.org/10.1186/s13041-020-00603-7

Download citation

Keywords

  • Glaucoma
  • Optic nerve head
  • Microglia
  • RNA-sequencing
  • DBA/2 J
  • Neuroinflammation
  • Mitochondria