MINI REVIEW article

Front. Cell. Neurosci., 20 April 2023

Sec. Non-Neuronal Cells

Volume 17 - 2023 | https://doi.org/10.3389/fncel.2023.1173200

Reactive astrogliosis in the era of single-cell transcriptomics

  • 1. Laboratory of Gene Expression, Institute of Biotechnology of the Czech Academy of Sciences, Vestec, Czechia

  • 2. Faculty of Science, Charles University, Prague, Czechia

  • 3. Department of Translational Neuroscience, University Medical Centre Utrecht Brain Centre, Utrecht University, Utrecht, Netherlands

  • 4. Laboratory of Astrocyte Biology and CNS Regeneration, Center for Brain Repair, Department of Clinical Neuroscience, Institute of Neuroscience and Physiology, Sahlgrenska Academy at the University of Gothenburg, Gothenburg, Sweden

  • 5. Florey Institute of Neuroscience and Mental Health, Parkville, VIC, Australia

  • 6. University of Newcastle, Newcastle, NSW, Australia

  • 7. TATAA Biocenter AB, Gothenburg, Sweden

  • 8. Department of Cellular Neurophysiology, Institute of Experimental Medicine of the Czech Academy of Sciences, Prague, Czechia

Abstract

Reactive astrogliosis is a reaction of astrocytes to disturbed homeostasis in the central nervous system (CNS), accompanied by changes in astrocyte numbers, morphology, and function. Reactive astrocytes are important in the onset and progression of many neuropathologies, such as neurotrauma, stroke, and neurodegenerative diseases. Single-cell transcriptomics has revealed remarkable heterogeneity of reactive astrocytes, indicating their multifaceted functions in a whole spectrum of neuropathologies, with important temporal and spatial resolution, both in the brain and in the spinal cord. Interestingly, transcriptomic signatures of reactive astrocytes partially overlap between neurological diseases, suggesting shared and unique gene expression patterns in response to individual neuropathologies. In the era of single-cell transcriptomics, the number of new datasets steeply increases, and they often benefit from comparisons and integration with previously published work. Here, we provide an overview of reactive astrocyte populations defined by single-cell or single-nucleus transcriptomics across multiple neuropathologies, attempting to facilitate the search for relevant reference points and to improve the interpretability of new datasets containing cells with signatures of reactive astrocytes.

1. Introduction

Astrocytes carry out functions essential for maintaining homeostasis in the CNS. They induce and control neuronal synapses, supply neurons with energy, and participate in neural plasticity mechanisms. Astrocytes also support the blood–brain barrier integrity and function and can react to cytokines circulating in the vascular system (Pekny and Pekna, ; Hasel and Liddelow, ; Lee et al., ). In disease or upon injury, astrocytes respond by reactive astrogliosis. This is characterized by increased expression of immune response genes and cytoskeletal and extracellular matrix (ECM) components, which is followed by morphological and functional changes (Wilhelmsson et al., ; Pekny et al., ; Escartin et al., ).

Reactive astrocytes are classically identified by hypertrophy of cellular processes and increased expression of the glial fibrillary acidic protein (GFAP) (Wilhelmsson et al., ; Hol and Pekny, ; Escartin et al., ). Although astrocytes show some common features in response to different pathological stimuli, their heterogeneity is far more complex and dynamic than previously thought (Stahlberg et al., ; Rusnakova et al., ; Orre et al., ; Pekny et al., ; Batiuk et al., ; Bayraktar et al., ; Gomes et al., ; Habib et al., ; Wheeler et al., ; Hasel et al., ; Sadick et al., ). We and others have recently cautioned against using a single or only a few selected marker genes to define astrocyte reactive states (Escartin et al., ). Instead, reactive populations of astrocytes should be defined considering multiple parameters, including gene expression, proteomics, morphology, and function. In addition, the astrocyte phenotype is determined by location in the CNS, disease pathogenesis and progression, co-morbidities, age, sex, and other sources of heterogeneity.

High-throughput single-cell and single-nucleus RNA sequencing (scRNA-seq/snRNA-seq), along with spatial transcriptomics, have been widely used to determine the diverse cellular composition of the CNS and to classify subsets of cells. These innovative technologies also allow for the determination of the diverse populations of cell types and cell states involved in the pathogenesis of brain disorders (Keren-Shaul et al., ; Habib et al., ; Kenigsbuch et al., ). The analysis of transcriptomic data often employs reference datasets, when annotating cellular clusters and uses data integration to increase resolution and identify concordances with previous studies (Stuart et al., ). The existing transcriptomic databases gather data obtained through various protocols from a wide range of specimens, allowing for the exploration of pre-processed data [e.g., PanglaoDB by Franzen et al. ()] or browsing through an easily accessible collection of references along with key metadata [e.g., single-cell studies database by Svensson et al. ()]. However, as the number of available datasets grows, it becomes important, albeit difficult, to search for relevant reference points.

Here, we present an overview of reactive astrocyte populations described in neurodegenerative diseases, acute injuries, and other pathological conditions using scRNA-seq and snRNA-seq technologies in rodents and humans. We summarize the key studies and the marker genes defining the populations of reactive astrocytes (Table 1, Supplementary Table 1) and identify the core of reactive astrocyte-associated genes that are repeatedly listed across datasets (Figure 1). We believe this compact resource will help to identify relevant reference datasets to study reactive astrocytes, and thereby contribute to a better understanding of astrocyte heterogeneity and function in disease.

Table 1

StudyClusterGene signatureNotes
Zamanian et al. () and Liddelow et al. ()Pan-reactiveLcn2, Steap4, S1pr3, Timp1, Hsbp1, Cxcl10, Cd44, Osmr, Cp, Serpina3n, Aspg, Vim, GfapReactive signatures as in Liddelow et al., , Figure 1a
A1H2.T23, Serping1, H2.D1, Ggta1, Iigp1, Gbp2, Fbln5, Ugt1a, Fkbp5, Psmb8, Srgn, Amigo2
A2Clcf1, Tgm1, Ptx3, S100a10, Sphk1, Cd109, Ptgs2, Emp1, Slc10a6, Tm4sf1, B3gnt5, Cd14
Leng et al. ()IRAS1FTH1, MT2A, CXCL1, CCL2, CXCL8, CCL20, CXCL6, CXCL3, C15orf48, MT1E, LGALS1, CSF3, SOD2, CXCL2, PTX3, SERPINE1, TMSB10, SAA1, IL32, SRGNTop 20 genes from Supplementary Table 4a, cluster 1, ordered based on avg_diff
IRAS2FTH1, TMSB10, CCL2, IL32, B2M, TMSB4X, LGALS1, SOD2, MT-ND3, HLA-B, NNMT, S100A11, PFN1, PRELID1, CALR, RPS27A, PSME2, CYP1B1, ICAM1, NFKBIATop 20 genes from Supplementary Table 4a, cluster 2, ordered based on avg_diff
Hasel et al. ()Cluster 4Lcn2, Ifitm3, Timp1, B2m, Vim, Gfap, Psmb8, Bst2Top genes enriched in cluster as in Figure 2c
Cluster 8Igtp, Gm4951, Ifit3, Iigp1, Irgm1, Tap1, Gbp2, Ifit1, Isg15, Cxcl10
Habib et al. ()Gfap-highMyoc, Gfap, Id3, Igfbp5, Aqp4, Serpinf1, Id1, Apoe, 6330403K07Rik, Slc38a1, Clu, Fbxo2, Vim, Ckb, Fabp7, Cd9, Ptchd2, Fxyd6, Unc5c, Scd2Top 20 genes from Supplementary Table 2, cluster 1 vs. 6, ordered by log2FC
DAARpl36-ps3, Gfap, Vim, Clu, Igfbp5, Id3, Gsn, Cd9, Fxyd1, Ctsb, Apoe, C4b, Kcnip4, Plce1, Cst3, Ggta1, Serpina3n, S100a6, Cnn3, Mt2Top 20 genes from Supplementary Table 2, cluster 1 vs. 4, ordered by log2FC
Morabito et al. ()GFAPhigh/CHI3L+CHI3L1, TPST1, RGS6, AC083837.2, ARHGEF3, ZFYVE28, ELL2, SAMD4A, GNA14, SLC44A3-AS1, IFI16, CPNE8, CD44, WWTR1, AC010655.4, ADAM12, ETV6, CLIC4, OSMR-AS1, PKNOX2Top 20 genes from Supplementary Data 1d, cluster ASC3, ordered by log2FC
Zhou et al. ()Astro0, Astro1PDE4DIP, PFKP, COL5A3, NCAN, SLC5A3Figure 5c, genes upregulated in AD clusters Astro0 and Astro1
Lee et al. ()A-C5aSerpina3n, Aqp4, Igfbp5, A2m, Gsn, Gfap, Cd44, Ctsb, Slc16a1, Ddr1, Serpinf1, C4b, Ctsd, Plat, Serping1, Timp2, Tagln3, Vim, Ccdc3, Mlc1Supplementary Table 8, Cohort II astrocyte
A-C5bCst3, Scrg1, Cd81, Reep5, Lxn, Tmco1, Cavin3, Nde1, Gstm5, Clu, Nudt10, Tmem50a, Gm44974, Tmem208, Cacng6, Gm12892, Gm38379, Hsd11b1, Ccdc13, Ccdc116Supplementary Table 8, Cohort II astrocyte
Grubman et al. ()a1LINGO1, HSPA1A, BCYRN1, BOK, DNAJB2, MT-CO2, MT-ND3, HSP90AA1, GFAP, MT-ND1, MT-ND2, CRYAB, HSPB1, MT-CYB, MT-ATP6, MT-ND4, MT-CO3, CHI3L1, DNAJB1, AEBP1Top 20 genes ordered by log2FC as at web interface: http://adsn.ddnetbio.com/, log2FC > 0.5, FDR < 0.05
a2XIST, MAOB, NEAT1, GOLGB1, GFAP, CD44, PLEKHA5, CH17-189H20.1, ID3, VCAN, PLCE1, HIF3A, DCLK1, RP11-507B12.2, MAP1B, GADD45G, TSHZ2, ADAMTSL3, GALNT15, RHPN1
Mathys et al. ()Ast1SLC26A3, MT1G, RASGEF1B, LINGO1, MT1M, MT1E, CSRP1, KCNJ10, VIM, C10orf54, MT2A, GJA1, GLUL, TMBIM6, CLU, WFS1, PBXIP1, PPP1R1B, SDC4, HSP90AB1Top 20 genes from Supplementary Table 6, ordered by log2FC
Lau et al. ()AD-up-regulatedPTOV1, KIAA0930, PRRX1, PLCG2, ACAP3, IRS2, RIMS1, GFAP, HSPA1B, PTGES3, TTTY14, LINC00609, VCAN, DLGAP1, LINGO1, MT-ND5, DNAJB1, HSPB1, DNAJB2, HMGB1, HSPH1, CRYAB, AEBP1, FP671120.1, NEAT1, MTRNR2L12, DHFR, HSPA1A, RASGEF1B, HSP90AA1, SLC26A3, FP236383.1Supplementary Figure 5d, top genes enriched in a1 and a6 populations
Leng et al. ()GFAPhighGFAP, CD44, TNC, HSPB1, HSP90AA1Figures 5c, d – expression of reactivity genes; markers of the cluster are not available
Sadick et al. ()Pan-astrocytic DEGsNEAT1, RANBP3L, PLCG2, SLC39A11, PFKP, SESN1, PLPP1, LRMDA, ABCA1, SLC25A18, RASSF8, LINC00511, WWOX, NPAS2, CPNE8Supplementary Table 6, overlap of DEGs in each cluster
Serrano-Pozo et al. ()astR1, astR2ACTB, APP, AQP4, CD44, CLU, CRYAB, ECE1, FGFR1, GFAP, HSPB1, HSPB8, LAMA1, MAOB, MAP1B, MAPT, MT1E, MT2A, PRDX6, S100B, VCANCommonly selected markers of reactive clusters from Figure 5b
Smajic et al. ()CD44highKCNE4, SHISA9, CP, EMP1, NEAT1, CHI3L1, TPD52L1, CACNB2, HPSE2, ATRNL1, HS6ST3, MYRIP, PLCXD3, SLCO3A1, ST6GAL1, MAN1C1, NFASC, CERS6, SYNPO2, AEBP1, DCLK1, CD44, GALNT15, S100A6, CRYAB, FTH1, CPAMD8, C8orf34, TTN, NRXN3, HSP90AA1, HSPA1A, HSPB1, UBCFigure 4F, astrocytic DEGs upregulated in PD overlapped with genes enriched in the astrogliosis trajectory toward CD44high
Al-Dalahmah et al. ()Cluster 1MT2A, MT1F, MT1E, MT1X, MT1G, MT1M, FTL, MT3, FAM171B, TPD52L1, EGLN3, CPE, GJB6, EPHX1, TKT, CLU, RGCC, NUPR1Table 2, selected markers
Cluster 2HSP90AA1, HSPB1, CRYAB, MTATP6P1, HSPA1A, UBC, CALM1, HSP90AB1, ATP1B1, GPM6B, CCK, VAMP2, MAP1A, NRGN, EEF1A1, EIF1, TMSB4X
Cluster 5MALAT1, GFAP, NEAT1, FP236383.2, SORBS1, PLP1, CEBPD, TSC22D1, PRRG3, SP100, HILPDA, COL6A3, RASD1, CHI3L1, SLC14A1, SAT1, HNRNPH1, PCDH9, PRKCA, BCL6, LINC00969, ANKRD36B, HIF3A, ITGB8, DTNA, HSPA1B, ETNPPL
Cluster 6ID2, S100B, FOS, DBI, HES1, JUNB, FABP7, H3F3B, ID3, PTGDS, MT1H, IFITM3
Wheeler et al. ()MAFGFyb, C1qc, Ccl4, Ly86, Lpcat2, Tyrobp, C1qa, Cyba, Rgs10, Ctsc, Cx3cr1, Bcl2a1b, Ccl2, Selplg, Hexb, Fcer1g, Cebpb, Ctsh, Mpeg1, Gm8995Top 20 genes from Supplementary Table 3, cluster 5, ordered by log2FC, padj < 0.05
Sanmarco et al. ()LAMP1+TRAIL+Fscn1, Egr1, Ier2, Tnfaip3, Il1a, Nfkbia, Btg2, Marcksl1, Gadd45b, Jun, Cst3, Cd83, Fos, Nfkbiz, Icam1, Ier5, Junb, Ppp1r15a, Tob2, SparcTop 20 genes from Supplementary Table 8, cluster 0, ordered by log2fc, padj < 0.05
Absinta et al. ()AIMSIGKC, VIM, APOE, CLU, CST3, ACTG1, FTL, CAPS, ACTB, CD81, ITM2C, SRP14, CEBPD, S100B, ZFP36L2, AGT, TUBA1A, SPARC, GAPDH, CKBTop 20 genes from Supplementary Table 7, cluster 6, ordered by log2FC
Reactive/stressed astrocytesNAMPT, PLOD2, HSPH1, HILPDA, TPST1, ELL2, TRIO, ARID5B, ABCA1, ZFYVE28, SERPINH1, GPR158, HSP90AA1, SAMD4A, CNN3, ACTN1, ARHGEF3, NOD1, FAM189A2, SPOCK1Top 20 genes from Supplementary Table 7, cluster 1, ordered by log2FC
RASGEF1B, MT-ND2, SLC26A3, MT-CO3, MT-ND1, MT-ND3, MT-CO2, MT-CO1, MT-ND4, MT-CYB, MT-ATP6, AC105402.3, MT-ND5, NEAT1, BNIP3L, LINC00511, SLC16A1-AS1, GNB1, GOLGA4, SGMS1Top 20 genes from Supplementary Table 7, cluster 5, ordered by log2FC
Shi et al. ()Cluster 3Apod, Ttr, Ier2, Igfbp5, Btg2, Cpe, Spp1, Txnip, Dnajb1, Cyr61, Klf2, Egr1, Gstm1, Mt1, Mt2, Nfkbia, Zfp36, Sox9, Junb, Cst3Top 20 genes from source data for Figure 10f, ordered by log2FC, padj < 0.05
Ma et al. ()AST12_CPlp1, Ptgds, Gfap, Mal, Mbp, Cldn11, Apod, Hsbp1, S100a10, MobpTop 10 markers in Figure 3i
AST24_CSpp1, Ccl14, Gfap, Cd14, Ccl12, Ccl3, Cxcl12, Lgals3, C5ar1, Ccl6Top 10 markers in Figure 4i
Zamboni et al. ()AC3Stat1, Bst2, Iigp1, Ifitm3, H2-T23, H2-D1Figure 3b, selected markers
Li et al. ()Cluster 4Aqp4, Gfap, S100a6, Vim, NesCluster signature genes from Figure 4b
Cluster 5S100a6, Vim, Nes, Tubb5, Mki67
Wang et al. ()Cluster 4Spp1, Cd44, Srbd1, Smarcd2, Slc16a3, Acsl4, Ccl12, Lrp1, PvrSupplementary Data Sheet 1, ordered by log2FC, padj < 0.05
Chancellor et al. ()Astrocyte2SLC14A1, HIF3A, PDK4, LINC01088, NPL, ZHX3, ADGRV1, PPP2R2B, KANSL1, CDH20, PARD3, MAPK10, MAGI2, PPM1E, CYP4F3, GRAMD2B, ECHDC2, MIR99AHG, RORA, ARL17BTop 20 genes from Supplementary File 3, Table 8, ordered by avg_diff
Astrocyte3SERPINA3, CD44, BCL6, GLIS3, IRS2, CHI3L1, CADPS, ANGPTL4, CLU, MAPK4, GLUL, VEGFA, ARHGEF3, CSNK2A1, MAN1C1, MT2A, EFCAB3, SIAH1, PFKP, KCTD21-AS1

Gene signatures of selected populations of reactive astrocytes (mostly based on single-cell studies).

AIMS, astrocytes inflamed in multiple sclerosis; avg_diff, average difference; DAA, disease-associated astrocytes; DEGs, differentially expressed genes; FDR, false discovery rate; log2FC, log2 fold change; padj, adjusted p-value; PD, Parkinson's disease. The signatures have been presented in the original studies, and the source is indicated in the notes. If additional filtering was applied to prioritize the genes, the thresholds are noted.

Figure 1

2. Neuroinflammation

Zamanian et al. () provided a landmark transcriptomic study defining distinct states of reactive astrocytes emerging during neuroinflammation and ischemic stroke in mice injected with bacterial lipopolysaccharide (LPS) and in the middle cerebral artery occlusion (MCAO) model. Reactive astrocytes in both conditions upregulated intermediate filament proteins (known also as nanofilament proteins), ECM, and adhesion molecules, regulators of metal ion homeostasis, and components of immune response pathways. However, while LPS reactive astrocytes showed signs of an elevated immune response, MCAO astrocytes upregulated genes involved in metabolic activity and proliferation, presumably connected with post-stroke neuroplasticity responses.

Lipopolysaccharide and MCAO reactive astrocytes, occasionally still termed A1 and A2 (Liddelow et al., ; Escartin et al., ), have become perceived as detrimental/pro-inflammatory or beneficial/anti-inflammatory subtypes. A1-like astrocytes were found to be induced by the cytokine trio Il-1α, TNF, and C1q secreted by microglia and were identified in post-mortem samples from patients suffering from Alzheimer's disease (AD), Parkinson's disease (PD), amyotrophic lateral sclerosis (ALS), and multiple sclerosis (MS) (Liddelow et al., ; Guttenplan et al., ). Although such a binary division is now accepted as oversimplified and insufficient to encompass the real spectrum of reactive astrocyte heterogeneity (Escartin et al., ), these early experiments provided a useful point of reference for later studies.

Using single-cell transcriptomics, Leng et al. () explored the heterogeneity of astrocytes derived from human-induced pluripotent stem cells (hiPSCs) treated by the cytokine trio used by Liddelow et al. () and defined two clusters of neuroinflammation-induced reactive astrocytes. While the IL-1/IL-6-responsive astrocytes (IRAS1) were characterized by acute phase response genes, the TNF/IFN-responsive astrocytes (IRAS2) upregulated interferon signaling, indicating their specialized function. Using an integrative analysis, both signatures were identified in an MS mouse model (Wheeler et al., ), in hiPSC-derived astrocytes stimulated with cytokines (Barbar et al., ), and in a study by Hasel et al. ().

Hasel et al. () addressed the astrocyte response in LPS-induced neuroinflammation in mice. A neuroinflammation-related fast-responding cluster was found almost exclusively in the LPS samples. As IRAS2, it was induced by IL1α, TNF, and C1q cytokines, regulated by Stat1 and Stat2, and showed differential expression of genes involved in the interferon response, antigen processing, and antigen presentation genes. Noteworthy, these astrocytes were localized close to blood vessels and around ventricles, suggesting their region-specific function, which cannot be observed in cell cultures (Leng et al., ). The signature of this cluster was found to be expressed also in models of AD (Habib et al., ; Zhou et al., ), MS (Wheeler et al., ), and stab wound injury (Zamboni et al., ). Another interesting cluster localized in the white matter (WM), was present in control and treated mice, and upon LPS treatment activated the expression of inflammatory genes. Together, these studies indicated that, in response to neuroinflammation, astrocytes form functionally and regionally diverse reactive populations.

3. Alzheimer's disease

The experimental neuroinflammation model provided valuable insights into astrocyte heterogeneity and predicted potentially higher complexity in neurodegenerative diseases. This has been confirmed by extensive AD research, as discussed in the following paragraphs. For a comprehensive overview of sc-/snRNA-seq AD studies without focus on astrocytes, refer to Cuevas-Diaz Duran et al. () and Saura et al. (), or explore the scREAD database (Jiang et al., ).

In a key study defining astrocyte reactivity in AD, Habib et al. () analyzed the hippocampus of the 5xFAD mouse using snRNA-seq. The authors identified two clusters of astrocytes with increased expression of Gfap, one of which was more abundant in 5xFAD samples. Accordingly, these were termed disease-associated astrocytes (DAAs). The DAAs partially shared the expression profile with Gfap-high astrocytes that were also present in the control samples, but the DAAs specifically upregulated genes related to endocytosis, complement, aging, and amyloid metabolism. Furthermore, the presence of homeostatic and intermediate populations suggested dynamic transitions of individual astrocyte states. DAA-like populations were also found in 5xFAD cortices and aged brains of mice and humans (Mathys et al., ), indicating that the signature is not specific to a single AD mouse model or brain region. These results were extended by Morabito et al. (), who detected Gfap-high and DAA signatures enriched among nuclei from prefrontal cortices of AD patients and using chromatin accessibility data, they identified FOSL2 as a potential regulator of the DAAs. Moreover, they found a corresponding cluster of GFAPhigh/CHI3L+ astrocytes proportionally increased in AD.

Zhou et al. () noticed inter-species differences in astrocyte changes between the 5xFAD mouse model and the prefrontal cortex of AD patients. While only a limited transcriptional change was reported in mice, several astrocyte clusters differed in proportions between human controls and AD samples. Of note, different genotypes of the microglial receptor TREM2, a known AD risk factor (Wolfe et al., ), were found to not affect reactive astrogliosis. Another analysis of the hippocampus of the AD mouse models of amyloidosis (PS2APP) and amyloidosis combined with tauopathy (TauPS2APP) revealed the presence of reactive astrocytes (Lee et al., ). Their comparison across conditions showed comparable transcriptional changes not only in PS2APP and TauPS2APP but also in a Trem2 knockout, indicating that astrocytes respond similarly in different models of AD and in a TREM2-independent fashion.

Focusing on AD patient samples, three early studies investigated the entorhinal cortex—the site of early pathological changes (Grubman et al., ) and the prefrontal cortex (Mathys et al., ; Lau et al., ) affected later in the disease progression. Interestingly, populations of AD-associated astrocytes identified in individual studies showed overlapping gene signatures, illustrated by common upregulation of heat shock proteins, LINGO1 or RASGEF1B. These results highlighted conserved transcriptomic changes in different brain regions affected in the early and late AD pathogenesis. This has been confirmed by an independent study comparing the entorhinal cortex and the superior frontal gyrus within the prefrontal cortex of AD patients at different stages of the pathology, identifying GFAPhigh astrocytes in comparable abundance across regions and stages (Leng et al., ).

Using an improved enrichment strategy (based on LHX2+NeuN sorting), Sadick et al. () increased the astrocyte proportion from ≈ 15% of nuclei in previous studies (estimated by the authors) to over 50%, allowing for an in-depth analysis of AD prefrontal cortices. Interestingly, no cluster was identified as disease-related; instead, each showed individual changes in AD samples. The global transcriptomic alteration in all astrocytes was represented by pan-astrocytic differentially expressed genes. Furthermore, multiple clusters upregulated metallothioneins, heat shock proteins, and other previously reported reactivity-related genes. Thus, reactive astrocytes might not be limited to a single population, but rather originate in multiple homeostatic populations.

A similar conclusion was reached in a comprehensive study by Serrano-Pozo et al. (), who analyzed five brain regions consecutively affected during AD progression. The authors identified two reactive astrocyte clusters characterized by a reactive signature. Importantly, analyses of multiple intermediate clusters and trajectories supported the homeostatic origin of reactive astrocytes, as proposed earlier by Habib et al. () and Sadick et al. ().

4. Other proteinopathies

Alzheimer's disease, PD, and Huntington's disease (HD) belong to neurodegenerative proteinopathies. PD and HD have also been analyzed with single-cell transcriptomics (Ma and Lim, ; Malla et al., ), but only a few studies have focused on astrocytes. Other proteinopathies with an important or causative role of astrocytes remain to be investigated, e.g., Alexander disease (Pajares et al., ).

In the midbrain of PD patients, Smajic et al. () identified CD44high astrocytes, representing the terminal activated state in astrogliosis trajectory. Astrocytes originating in PD patients were enriched at the end of this activation trajectory, suggesting their role in PD-associated neuroinflammation. Al-Dalahmah et al. () analyzed the cingulate cortex of patients suffering from HD and identified three HD clusters of astrocytes, expressing variable levels of metallothioneins and GFAP, possibly transitioning one into another. Interestingly, similarly to Habib et al. (), they also identified GFAP high astrocytes in controls, implying the existence of astrocyte heterogeneity even in healthy individuals. Collectively, these PD and HD studies together with the AD studies show the presence of similar reactive astrocyte populations across multiple disorders.

5. Multiple sclerosis

Although single-cell transcriptomic studies of MS have focused mainly on neurons, microglia, and oligodendrocytes (Falcao et al., ; Jakel et al., ; Schirmer et al., ), astrocyte heterogeneity has not remained unnoticed [reviewed in Lo et al. ()]. Using a model of experimental autoimmune encephalomyelitis (EAE), Wheeler et al. () identified a population of astrocytes expressing pro-inflammatory genes. A similar population was also found in active WM lesions in MS patients. Due to high levels of the transcriptional regulator MAFG, this population was termed MAFG-driven astrocytes. Integration revealed the presence of MAFG astrocytes in other MS datasets, regardless of the CNS region analyzed (Lake et al., ; Jakel et al., ; Schirmer et al., ). Soon after, Hasel et al. () and Leng et al. () detected the signature of the neuroinflammatory populations of astrocytes in the dataset produced by Wheeler et al. (), showing a shared response of mouse and human astrocytes in experimental neuroinflammation and MS.

A follow-up study identified the LAMP1+TRAIL+ astrocyte population characterized by anti-inflammatory and pro-apoptotic properties, which was depleted in EAE samples (Sanmarco et al., ), and was proposed to act protectively. The impaired balance of MAFG and LAMP1+TRAIL+ astrocytes in EAE/MS may lead to outnumbering of the protective astrocytes by the pro-inflammatory population, likely exaggerating the inflammatory environment in the affected tissue.

Considering the spatial aspect of MS, Absinta et al. () identified reactive astrocytes enriched in multiple areas around and within MS lesions, expressing heat shock proteins, cell death regulators, and mitochondrial genes. Another population of astrocytes inflamed in MS upregulated GFAP, APOE, VIM, and several A1 markers. It was found almost exclusively at the edges of chronic active lesions, showing that distinct reactive astrocytic populations arise in MS depending on lesion localization and its severity.

6. Neurotrauma and stroke

Neurotrauma and stroke represent another group of CNS disorders triggering substantial cellular heterogeneity that has been investigated by single-cell transcriptomics (Wilhelmsson et al., ; Moulson et al., ; Cao et al., ).

Shi et al. () reported varying phagocytic properties of astrocytes in mouse models of hemorrhagic and ischemic strokes. Reactive astrocytes characterized by upregulated synapse pruning and lysosome-related processes were shown to be responsible for this difference. Although they were present in both conditions, their proportion decreased from 20% in ischemic to <2% in hemorrhagic mice, suggesting model-specific patterns in astrocyte reactivity. Another transcriptomic analysis of the early phases of ischemic stroke in mice identified a cluster of astrocytes enriched in MCAO samples (Ma et al., ). Interestingly, their transcriptional signature at 12 h and 24 h after stroke differed, transitioning from a state with upregulated oxidative phosphorylation and cell junction genes to a state defined by an increased immune and inflammatory response, implicating a dynamic nature of astrocyte reactivity in the early phases after the injury.

Zamboni et al. () investigated the neurogenic potential of astrocytes in mouse somatosensory cortex subjected to stab wound injury. They identified a small cluster of reactive astrocytes expressing A1 and immune response genes. Importantly, trajectory analysis revealed the separation of the reactive and neurogenic branches, showing that they represent two independent and non-sequential states.

Two distinct clusters of reactive astrocytes were identified in a mouse model of spinal cord injury (SCI) (Li et al., ). Both were defined by an increased expression of Gfap, Vim, and Nes, but one cluster also expressed Mki67 at high levels, indicating that some reactive astrocytes are proliferative, peaking 3 days after the injury. The SCI in a rat model induced a reactive astrocyte population expressing immune response-related genes (Wang et al., ), such as Cd44, which has been associated with reactive astrocytes in neurodegeneration (Grubman et al., ; Lee et al., ; Leng et al., ; Smajic et al., ), and Spp1, which characterizes ischemia-induced astrocytes (Shi et al., ; Ma et al., ). Collectively, these represent examples of shared reactive signatures in different pathophysiological conditions and species.

The heterogeneity of reactive astrocytes has been investigated also in chronic traumatic encephalopathy (CTE), a tauopathy induced by repetitive head impacts (Chancellor et al., ). The authors identified two astrocyte populations enriched in CTE. While one population was characterized by genes associated with a dysfunctional metabolism, the other expressed genes associated with neuroinflammation and reactivity. In addition to CTE, traumatic brain injury can increase the risk of neurodegenerative diseases or induce post-traumatic epilepsy, with astrocytes playing a pivotal role in these secondary processes (Graham and Sharp, ; Small et al., ). The exact mechanisms remain to be elucidated, but the identification of astrocyte genes responsible for the onset of post-traumatic complications might largely contribute to the discovery of biomarkers for severity prediction, their modulation, and mitigation of deleterious consequences.

7. Inter-species comparison

Several single-cell transcriptomic studies indicate shared signatures of reactive astrocytes across neurological diseases, e.g., MS, AD, and neurotrauma (Hasel et al., ; Leng et al., ). Moreover, overlapping signatures were found in human and mouse reactive astrocytes by Habib et al. () and Leng et al. (), while others also noticed inter-species differences (Zhou et al., ).

To characterize shared and distinct aspects of astrocyte reactivity across species, Li et al. () compared the responses of mouse and human astrocytes to various pathological stimuli. While they found partial conservation of astrocyte gene expression, they also reported significant differences in the expression of genes involved in mitochondrial metabolism and defense response. Human astrocytes were susceptible to oxidative stress damage, which in mice was mitigated by efficient detoxification and mitochondrial resilience. Furthermore, in hypoxia, mouse astrocytes promoted neuronal growth and repair, representing an advantage in regeneration capacity. Human astrocytes, on the other hand, upregulated inflammatory response and antigen presentation upon stimulation with a TNFα and double-stranded RNA. Further comparison with AD and MS suggested that the inflammatory response of astrocytes is a common denominator in multiple neuropathologies.

While mouse models provide important insight into the biology of reactive astrocytes, inter-species differences must be considered when extrapolating mouse data to human neurological diseases.

8. Conclusion

We provided an overview of single-cell transcriptomic studies characterizing reactive astrocytes in various CNS diseases to facilitate the interpretation of new datasets and to improve our molecular understanding of reactive astrogliosis. There is a shared “reactive signature” in multiple pathologies, suggesting the existence of a “pan-reactive” expression program, which is partly conserved across species. Integration and meta-analyses are needed to understand the similarities and differences and to assess their translational potential. Genes expressed in reactive astrocytes were also found in other glial cells, including microglia and oligodendrocytes (Keren-Shaul et al., ; Wilhelmsson et al., ; Lee et al., ; Kaya et al., ; Kenigsbuch et al., ), implicating common activation mechanisms also across cell types. Future studies focusing on cellular populations will characterize context-specific formation and functions of reactive astrocytes and describe their interplay with other cell types. Large-scale analyses of single-cell datasets could then identify targets for interventions, facilitate the search for suitable drug candidates (Xu et al., ), and improve the treatment of acute and chronic CNS diseases.

Statements

Author contributions

ZM and LV wrote the manuscript and prepared the figures. EH, MP, and MK reviewed and edited the complete manuscript. All authors have read and agreed on the final version of the manuscript.

Funding

LV, ZM, and MK have been supported by the Ministry of Education, Youth, and Sports, under the frame of EJP RD, the European Joint Programme on Rare Diseases: CZ.1.05/1.1.00/02.0109, RVO 86652036, the European Union's Horizon 2020 research and innovation program under the EJP RD COFUND- EJP N° 825575, Czech Science Foundation GACR 23-05327S and 23-06269S, and project LX22NPO5107 (MEYS): Financed by EU—Next Generation EU. EH has been supported by grants from Alzheimer Nederland (WE.03-2017-4), NWO Gravitation BRAINSCAPES (463002004), ZonMW Memorabel (733050816), ZonMW TOP (91217035), EJP RD ALEXANDER (JTC2019/ZonMW 463002004), and la Caixa Foundation LCF/PR/HR21/52410002 (Astromad). MP has been supported by the grants from Swedish Research Council (2017-02255, 2019-00284, and 2020-01148), the ALF (724421 and 965939), W. and M. Lundgren's Foundation, The Swedish Brain Foundation (FO2018-0252 and FO2021-0082), The Swedish Stroke Foundation, Hagströmer's Foundation Millennium, T. Söderberg's Foundations, la Caixa Foundation LCF/PR/HR21/52410002 (Astromad), and P. Eriksson's, E. Jacobson's, R., and U. Amlöv's Foundations.

Conflict of interest

MK is employed by TATAA Biocenter AB.

The remaining 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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

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

References

Summary

Keywords

reactive astrogliosis, astrocytes, cell populations, single-cell RNA-seq, neuroinflammation, neurodegeneration, CNS diseases

Citation

Matusova Z, Hol EM, Pekny M, Kubista M and Valihrach L (2023) Reactive astrogliosis in the era of single-cell transcriptomics. Front. Cell. Neurosci. 17:1173200. doi: 10.3389/fncel.2023.1173200

Received

24 February 2023

Accepted

27 March 2023

Published

20 April 2023

Volume

17 - 2023

Edited by

Miren Revuelta, University of the Basque Country, Spain

Reviewed by

Brandon Peter Lucke-Wold, University of Florida, United States

Updates

Copyright

*Correspondence: Lukas Valihrach

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics