Advertisement
Research Article| Volume 14, ISSUE 1, P183-191, January 2021

Download started.

Ok

Low intensity repetitive magnetic stimulation reduces expression of genes related to inflammation and calcium signalling in cultured mouse cortical astrocytes

  • Darren Clarke
    Correspondence
    Corresponding author. Experimental and Regenerative Neuroscience, School of Biological Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia.
    Affiliations
    Experimental and Regenerative Neuroscience, School of Biological Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia

    Perron Institute for Neurological and Translational Science, Nedlands, WA, 6009, Australia
    Search for articles by this author
  • Jamie Beros
    Affiliations
    Experimental and Regenerative Neuroscience, School of Biological Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia

    Perron Institute for Neurological and Translational Science, Nedlands, WA, 6009, Australia
    Search for articles by this author
  • Kristyn A. Bates
    Affiliations
    Experimental and Regenerative Neuroscience, School of Biological Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia
    Search for articles by this author
  • Alan R. Harvey
    Affiliations
    Perron Institute for Neurological and Translational Science, Nedlands, WA, 6009, Australia

    School of Human Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia
    Search for articles by this author
  • Alexander D. Tang
    Affiliations
    Experimental and Regenerative Neuroscience, School of Biological Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia

    Perron Institute for Neurological and Translational Science, Nedlands, WA, 6009, Australia
    Search for articles by this author
  • Jennifer Rodger
    Affiliations
    Experimental and Regenerative Neuroscience, School of Biological Sciences, The University of Western Australia, Nedlands, WA, 6009, Australia

    Perron Institute for Neurological and Translational Science, Nedlands, WA, 6009, Australia
    Search for articles by this author
Open AccessPublished:December 23, 2020DOI:https://doi.org/10.1016/j.brs.2020.12.007

      Highlights

      • 1 Hz and 10 Hz low intensity rMS induced reductions in astrocyte gene expression in vitro, 5 h post-stimulation.
      • Modulated genes were primarily related to inflammation and calcium signalling pathways.
      • In a subset of genes modulated by LI-rMS, protein levels were also modulated at 5 and 24 h post-stimulation.

      Abstract

      Repetitive transcranial magnetic stimulation (rTMS) is a form of non-invasive brain stimulation frequently used to induce neuroplasticity in the brain. Even at low intensities, rTMS has been shown to modulate aspects of neuronal plasticity such as motor learning and structural reorganisation of neural tissue. However, the impact of low intensity rTMS on glial cells such as astrocytes remains largely unknown. This study investigated changes in RNA (qPCR array: 125 selected genes) and protein levels (immunofluorescence) in cultured mouse astrocytes following a single session of low intensity repetitive magnetic stimulation (LI-rMS – 18 mT). Purified neonatal cortical astrocyte cultures were stimulated with either 1Hz (600 pulses), 10Hz (600 or 6000 pulses) or sham (0 pulses) LI-rMS, followed by RNA extraction at 5 h post-stimulation, or fixation at either 5 or 24-h post-stimulation. LI-rMS resulted in a two-to-four-fold downregulation of mRNA transcripts related to calcium signalling (Stim1 and Orai3), inflammatory molecules (Icam1) and neural plasticity (Ncam1). 10Hz reduced expression of Stim1, Orai3, Kcnmb4, and Ncam1 mRNA, whereas 1Hz reduced expression of Icam1 mRNA and signalling-related genes. Protein levels followed a similar pattern for 10Hz rMS, with a significant reduction of STIM1, ORAI3, KCNMB4, and NCAM1 protein compared to sham, but 1Hz increased STIM1 and ORAI3 protein levels relative to sham. These findings demonstrate the ability of 1Hz and 10Hz LI-rMS to modulate specific aspects of astrocytic phenotype, potentially contributing to the known effects of low intensity rTMS on excitability and neuroplasticity.

      Keywords

      Introduction

      With the increasing use of repetitive transcranial magnetic stimulation (rTMS) for treating a wide range of neurological and neuropsychiatric conditions, there is a growing need to better understand the effects of electromagnetic stimulation on the brain. The forms of rTMS most commonly used in the clinic generate action potentials in neurons and elicit synaptic plasticity in the form of long-term potentiation (LTP) and long-term depression (LTD) [
      • Lenz M.
      • Müller-Dahlhaus F.
      • Vlachos A.
      Cellular and molecular mechanisms of rTMS-induced neural plasticity.
      ,
      • Tang A.D.
      • Thickbroom G.
      • Rodger J.
      Repetitive transcranial magnetic stimulation of the brain: mechanisms from animal and experimental models.
      ,
      • Thickbroom G.W.
      Transcranial magnetic stimulation and synaptic plasticity: experimental framework and human models.
      ]. Interestingly, rTMS at stimulation intensities below action potential threshold (low intensity, LI) also alters neuroplasticity in humans [
      • Martiny K.
      • Lunde M.
      • Bech P.
      Transcranial low voltage pulsed electromagnetic fields in patients with treatment-resistant depression.
      ,
      • Robertson J.A.
      • Théberge J.
      • Weller J.
      • Drost D.J.
      • Prato F.S.
      • Thomas A.W.
      Low-frequency pulsed electromagnetic field exposure can alter neuroprocessing in humans.
      ] and rodents [
      • Makowiecki K.
      • Harvey A.R.
      • Sherrard R.M.
      • Rodger J.
      Low-intensity repetitive transcranial magnetic stimulation improves abnormal visual cortical circuit topography and upregulates BDNF in mice.
      ,
      • Rodger J.
      • Mo C.
      • Wilks T.
      • Dunlop S.A.
      • Sherrard R.M.
      Transcranial pulsed magnetic field stimulation facilitates reorganization of abnormal neural circuits and corrects behavioral deficits without disrupting normal connectivity.
      ,
      • Tang A.D.
      • Bennett W.
      • Hadrill C.
      • Collins J.
      • Fulopova B.
      • Wills K.
      • et al.
      Low intensity repetitive transcranial magnetic stimulation modulates skilled motor learning in adult mice.
      ], suggesting that supra-threshold stimulation (i.e. induced action potential firing) may not be necessary to induce plasticity. LI-rTMS (<120 mT) alters many types of neural function, from molecular to network-level effects, that are similar to those observed at high intensity. In vitro, LI-repetitive magnetic stimulation (rMS, no cranium) at high frequency facilitates excitability of pyramidal neurons in mouse cortical slices [
      • Tang A.D.
      • Hong I.
      • Boddington L.J.
      • Garrett A.R.
      • Etherington S.
      • Reynolds J.N.
      • et al.
      Low-intensity repetitive magnetic stimulation lowers action potential threshold and increases spike firing in layer 5 pyramidal neurons in vitro.
      ] and modulates calcium release in cultured cortical neurons [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ]. In vivo, LI-rTMS increases neurotransmitter levels in rats [
      • Seewoo B.J.
      • Feindel K.W.
      • Etherington S.J.
      • Rodger J.
      Frequency-specific effects of low-intensity rTMS can persist for up to 2 weeks post-stimulation: a longitudinal rs-fMRI/MRS study in rats.
      ], and alters motor learning in mice [
      • Tang A.D.
      • Bennett W.
      • Hadrill C.
      • Collins J.
      • Fulopova B.
      • Wills K.
      • et al.
      Low intensity repetitive transcranial magnetic stimulation modulates skilled motor learning in adult mice.
      ] similar to the effects of HI-rTMS in humans [
      • Stöckel T.
      • Summers J.J.
      • Hinder M.R.
      Reversed effects of intermittent theta burst stimulation following motor training that vary as a function of training-induced changes in corticospinal excitability.
      ,
      • Teo J.T.H.
      • Swayne O.B.C.
      • Cheeran B.
      • Greenwood R.J.
      • Rothwell J.C.
      Human theta burst stimulation enhances subsequent motor learning and increases performance variability.
      ]. Because LI-rTMS effects are elicited independently of action potential induction and involve universal cellular machinery such as calcium signalling, these effects may be mediated by glia as well as neurons.
      Astrocytes are important modulators of neuroplasticity, primarily through the release of gliotransmitters at the synapse [
      • Kang J.
      • Jiang L.
      • Goldman S.A.
      • Nedergaard M.
      Astrocyte-mediated potentiation of inhibitory synaptic transmission.
      ,
      • Papouin T.
      • Dunphy J.
      • Tolman M.
      • Foley J.C.
      • Haydon P.G.
      Astrocytic control of synaptic function.
      ,
      • Pascual O.
      • Casper K.B.
      • Kubera C.
      • Zhang J.
      • Revilla-Sanchez R.
      • Sul J.-Y.
      • et al.
      Astrocytic purinergic signaling coordinates synaptic networks.
      ,
      • Perea G.
      • Araque A.
      Astrocytes potentiate transmitter release at single hippocampal synapses.
      ,
      • Serrano A.
      • Haddjeri N.
      • Lacaille J.-C.
      • Robitaille R.
      GABAergic network activation of glial cells underlies hippocampal heterosynaptic depression.
      ]. However, the contributions of astrocytes and other glial cells to the effects of non-invasive brain stimulation have only recently been considered [
      • Cullen C.L.
      • Young K.M.
      How does transcranial magnetic stimulation influence glial cells in the central nervous system?.
      ,
      • Gellner A.-K.
      • Reis J.
      • Fritsch B.
      Glia: a neglected player in non-invasive direct current brain stimulation.
      ]. In vitro, 1Hz (but not 10Hz, cTBS or BHFS) LI-rMS (18 mT) increases intracellular calcium release in the cytoplasm of purified cultures of mouse cortical astrocytes [
      • Clarke D.
      • Penrose M.A.
      • Penstone T.
      • Fuller-Carter P.I.
      • Hool L.C.
      • Harvey A.R.
      • et al.
      Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
      ]. Similarly, various frequencies of HI-rTMS and LI-rTMS alter the expression of GFAP (a commonly used cytoskeletal marker of astrocytic reactivity) and the number of GFAP-positive astrocytes in numerous in vivo rodent neurotrauma and neurodegeneration models [
      • Cacace F.
      • Mineo D.
      • Viscomi M.T.
      • Latagliata E.C.
      • Mancini M.
      • Sasso V.
      • et al.
      Intermittent theta-burst stimulation rescues dopamine-dependent corticostriatal synaptic plasticity and motor behavior in experimental parkinsonism: possible role of glial activity.
      ,
      • Clarke D.
      • Penrose M.A.
      • Harvey A.R.
      • Rodger J.
      • Bates K.A.
      Low intensity rTMS has sex-dependent effects on the local response of glia following a penetrating cortical stab injury.
      ,
      • Rauš S.
      • Selaković V.
      • Manojlović-Stojanoski M.
      • Radenović L.
      • Prolić Z.
      • Janać B.
      Response of hippocampal neurons and glial cells to alternating magnetic field in gerbils submitted to global cerebral ischemia.
      ,
      • Sasso V.
      • Bisicchia E.
      • Latini L.
      • Ghiglieri V.
      • Cacace F.
      • Carola V.
      • et al.
      Repetitive transcranial magnetic stimulation reduces remote apoptotic cell death and inflammation after focal brain injury.
      ] and in vitro cell cultures [
      • Chan P.
      • Eng L.F.
      • Ling Lee Y.
      • Lin V.W.
      Effects of pulsed magnetic stimulation on GFAP levels in cultured astrocytes.
      ]. HI-rMS has been recently shown to affect astrocyte inflammatory gene expression in an in vitro ischemic model [
      • Hong Y.
      • Liu Q.
      • Peng M.
      • Bai M.
      • Li J.
      • Sun R.
      • et al.
      High-frequency repetitive transcranial magnetic stimulation improves functional recovery by inhibiting neurotoxic polarization of astrocytes in ischemic rats.
      ], and reduce the expression of proteins linked to depression in cultured astrocytes [
      • Peng Z.w.
      • Xue F.
      • Zhou C.h.
      • Zhang R.g.
      • Wang Y.
      • Liu L.
      • et al.
      Repetitive transcranial magnetic stimulation inhibits Sirt1/MAO-A signaling in the prefrontal cortex in a rat model of depression and cortex-derived astrocytes.
      ]. However, almost nothing is known about the mRNA and protein expression changes that occur following LI-rTMS and how they might influence astrocytic phenotype.
      Here, the effect of LI-rMS on gene expression was examined in purified cortical astrocyte cultures generated from neonatal mouse pups. Astrocyte cultures were stimulated with low (1Hz) and high (10Hz) frequencies at a low intensity (18 mT) and RNA was extracted 5 h following stimulation to examine 125 genes of interest relating broadly to four main astrocytic biological functions: calcium signalling, inflammation and injury, plasticity, and morphological/cytoskeletal composition. We further studied if gene expression changes resulted in protein expression changes in a selection of regulated genes at 5 and 24 h post-stimulation.

      Methods

      Tissue culture

      Experimental procedures were conducted with ethics approval from The University of Western Australia Animal Ethics Committee (RA/3/100/1462) under the Australian Code of Practice for the Use of Animals in Research. Primary cortical cultures were prepared from euthanased (>160 mg/kg sodium pentobarbitone, i.p.) 1 day postnatal male and female C57Bl6J mice (n = 31 over 6 culture runs), obtained from the Animal Resources Centre, Western Australia. Cultures were prepared as previously shown [
      • Clarke D.
      • Penrose M.A.
      • Penstone T.
      • Fuller-Carter P.I.
      • Hool L.C.
      • Harvey A.R.
      • et al.
      Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
      ] and are described in detail in the supplementary materials. Cells were grown in T75 flasks until confluent (day in vitro 11, DIV11), and then detached (0.25% trypsin, Gibco) and replated at a density of 0.04 × 106 onto 13 mm diameter glass coverslips in 24-well plates. To ensure an enriched astrocyte culture, coverslips were assessed for the proportions of various cell types. Cells at DIV14 were consistently >97% positive for GFAP across runs, with the remainder of cells immunohistochemically and morphologically identified as neurons, microglia, or fibroblasts (data not shown).

      Repetitive magnetic stimulation

      On DIV14, cells were stimulated with LI-rMS or sham within the cell culture incubator [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ]. Custom-built round coils (inside diameter 8 mm, outside diameter 17 mm, thickness 9.5 mm) were attached to the bottom of the 24-well plates to stimulate the cells. The coils were powered by a 12V magnetic pulse generator (Global Energy Medicine). The applied frequencies were 1Hz (600 pulses, 10 min, ‘10m1Hz’), 10Hz (600 pulses, 1 min, ‘1m10Hz’; or 6000 pulses 10 min, ‘10m10Hz’), or sham (10 min, generator switched off, no pulses delivered). 1Hz LI-rMS was the primary focus of this study because of its previous effects on calcium levels in astrocytes in vitro [
      • Clarke D.
      • Penrose M.A.
      • Penstone T.
      • Fuller-Carter P.I.
      • Hool L.C.
      • Harvey A.R.
      • et al.
      Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
      ] and number of GFAP-positive astrocytes after injury in vivo [
      • Clarke D.
      • Penrose M.A.
      • Harvey A.R.
      • Rodger J.
      • Bates K.A.
      Low intensity rTMS has sex-dependent effects on the local response of glia following a penetrating cortical stab injury.
      ]. Although 10Hz did not elicit calcium changes in our previous work, applying 10Hz offered us the opportunity to pulse match with 1Hz. In addition, 10 Hz delivered at high intensity has been recently shown to alter astrocyte gene expression in an in vitro ischemic model [
      • Hong Y.
      • Liu Q.
      • Peng M.
      • Bai M.
      • Li J.
      • Sun R.
      • et al.
      High-frequency repetitive transcranial magnetic stimulation improves functional recovery by inhibiting neurotoxic polarization of astrocytes in ischemic rats.
      ] and reduce SIRT1 and MAO-A protein levels and AMPKα phosphorylation in vitro [
      • Peng Z.w.
      • Xue F.
      • Zhou C.h.
      • Zhang R.g.
      • Wang Y.
      • Liu L.
      • et al.
      Repetitive transcranial magnetic stimulation inhibits Sirt1/MAO-A signaling in the prefrontal cortex in a rat model of depression and cortex-derived astrocytes.
      ]. The field intensity at the surface of the coverslip above the centre of the coil was measured with a Hall probe to be 18 mT, and each frequency delivered a monophasic pulse with a rise time of 320μs, equating to a dB/dt of approximately 56T/s. The coils do not generate heat or vibration above background levels during stimulation [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ].

      Analysis of gene expression

      Changes in gene expression were examined via RT-qPCR for each frequency to investigate downstream changes relating to calcium signalling, inflammation and injury, plasticity, and morphological/cytoskeletal composition. RNA was collected from 3 biological replicates, each pooled from 45 coverslips (total 135), for each stimulation group, 5 h after the end of stimulation. The 5 h timepoint was chosen for methodological consistency with a previous experiment conducted in our laboratory investigating mRNA changes in primary neuronal cultures following LI-rMS [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ]. RNA was extracted and purified from the astrocytes on coverslips using a miRNeasy kit (Qiagen), with DNase digestion (RNase-free DNase Set, Qiagen). RNA purity and quantity were measured (NanoDrop) prior to cDNA synthesis (RT2 First Strand Kit, Qiagen) from 1 μg of RNA for each biological replicate within each stimulation group. cDNA for each sample was added to custom RT2 Profiler PCR Arrays and amplified on a Bio-Rad CFX384. The custom array consisted of 125 genes known to be influential in astrocyte biology and function and focused on 4 main biological processes (see supplementary information). The amplification protocol was 10 min at 95 °C followed by 60 cycles of 95 °C for 15s and 60 °C for 1 min. Following this, a melt curve was obtained by increasing the temperature from 65 °C to 95 °C in 1 °C increments. An identical fluorescence threshold value for the exponential phase of amplification was used for each run, and genes were excluded if they passed the threshold after cycle 40 or had abnormal melt curves. The genes Brca1, Cga, Crh, FoxA2, Grm1, and Hnf1a were excluded from analysis as they reached threshold after cycle 40.

      Protein expression

      Changes in protein expression for selected genes were examined via immunofluorescent staining. Coverslips were stained for protein products of genes that were significantly downregulated compared to sham stimulation: antibodies to STIM1 (cat PA5-82455, ThermoFisher), ORAI3 (cat. MA5-15778, ThermoFisher), KCNMB4 (cat. PA5-77611) and CD56 (additional protein name of NCAM1, cat. MA5-11563, ThermoFisher). Five coverslips per stimulation group replicate were analysed for protein intensity at 5h and 24 h following LI-rMS. See supplementary material for further detail on staining procedure and image analysis. For consistency, the coverslips for all groups were stained in one session for each labelling combination.

      Statistical analysis

      mRNA and protein data were examined for outliers and normality with SPSS (v21, IBM). RT-qPCR results were normalized on Microsoft Excel and analysed with SPSS. The geometric mean of genes Gapdh, Sdha, and Sod2 were used to normalize genes for analysis and the normalized mean expression levels (log2(2-ΔΔCt)) for each stimulation group were used to determine expression differences between stimulation groups [
      • Pfaffl M.W.
      A new mathematical model for relative quantification in real-time RT–PCR.
      ,
      • Pfaffl M.W.
      • Tichopad A.
      • Prgomet C.
      • Neuvians T.P.
      Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper – Excel-based tool using pair-wise correlations.
      ]. For each gene, a fold change value was used to compare expression levels between groups with a student’s t-test. Fold changes (log2) greater than ±1 with a p value < 0.05 were considered statistically significant. All significantly regulated genes above the fold change cut-off were further examined using gene set enrichment analysis via Webgestalt [
      • Liao Y.
      • Wang J.
      • Jaehnig E.J.
      • Shi Z.
      • Zhang B.
      WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs.
      ]. Genes were ranked by ± fold change ∗ –log10(p value), and a false discovery rate (FDR) threshold of 0.05 was applied. Protein expression was analysed with a generalised linear mixed model (log inverse gaussian distribution) with the number of replicates, coverslips and cell density applied as covariates. Stimulation group, time, and coverslip location were treated as fixed effects. P values were considered significant <0.05, with Sidak corrections for post-hoc analysis and multiple comparisons. Astrocyte protein level changes were graphed with Prism (GraphPad). Correlation (two-tailed) was assessed using Spearman’s rho, with p value significance determined as < 0.05.

      Results

      rMS stimulation decreased expression of select genes

      Among the 125 genes examined, a total of 21 genes were identified with significantly regulated expression levels compared to sham 5 h following a single stimulation session (Fig. 1, summarised in Table 1). All statistically significant genes above the fold change cut-off were downregulated compared to sham stimulation; there were no significantly upregulated transcripts following LI-rMS. All frequencies induced significant changes in gene expression, with the most consistently affected pathway across frequencies related to inflammation and injury; all stimulation protocols resulted in reduced expression of at least 3 inflammatory-related genes. There was considerable overlap of genes that were modulated by both 10Hz groups (i.e. 1m10Hz and 10m10Hz), but there were also genes that were modulated by only one of the 10Hz groups (Table 1 and Fig. 1). The two 10 min protocols (1Hz and 10Hz) shared one downregulated gene in common, the intermediate filament keratin (Krt19), and there were no genes in common that were significantly regulated following 600 pulses at either 1Hz or 10Hz. These findings indicate the effects of LI-rMS on astrocyte mRNA expression are different for 1Hz and 10Hz, and are influenced by the number of pulses delivered. The expression changes for all 125 genes are indicated in Supplementary Table 1. The analysis of enriched biological processes and molecular functions across the significantly regulated genes for each stimulation protocol is detailed in the supplementary material. The main protein activity and cellular compartmentalisation of the significantly regulated genes is indicated in Fig. 1.
      Fig. 1
      Fig. 1Downregulation of gene expression in astrocytes 5 h following different LI-rMS protocols. A Heatmap of the 21 significantly downregulated genes compared to sham stimulation. Colour values represent log2 fold change, asterisk indicates p value < 0.05. B Venn diagram illustrating the frequencyspecific changes in gene expression, with each gene coloured to a major biofunction of its protein product in astrocytes. C Representation of the protein compartmentalisation of the significantly regulated genes following rMS. Product is colour coded by the frequency protocol(s) that influenced it, and protein shape denotes functional role. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)
      Table 1Significantly regulated genes in astrocyte cultures 5 h following a single session of LI-rMS. Each row is the expression of a gene significantly regulated by at least one LI-rMS protocol. mRNA details and gene function are provided. The function denotes prominent functional classification of gene product in astrocytes: C cytoskeletal, I inflammation and injury, P plasticity, S signalling. The frequency of pulses and duration of stimulation are displayed, and the log2 fold change and associated p value in comparison to sham-stimulated coverslips are presented for each gene. Significant log2 fold change in noted in bold.
      1m 10Hz10m 10Hz10m 1Hz
      REFSEQ MRNAGene symbol and nameFunctionFCpFCpFCp
      NM_007498Atf3, Activating transcription factor 3I−0.170.91−1.140.07−1.120.02
      NM_007529Bcan, BrevicanI, P0.570.63−2.060.010.730.21
      NM_019952Clcf1, Cardiotrophin-like cytokineI0.600.26−1.470.01−0.100.91
      factor 1
      NM_009140Cxcl2, C-X-C motif chemokine ligand 2I−2.210.010.780.57−0.890.09
      NM_010516Ccn1, Cellular communication network factor 1I, P−2.470.04−0.890.38−0.450.68
      NM_013642Dusp1, Dual specificity phosphatase 1I−1.700.04−0.420.27−0.230.75
      NM_010332Ednra, Endothelin receptor type AI, S0.250.69−1.120.27−1.740.04
      NM_007913Egr1, Early growth response 1I−1.430.01−0.870.11−0.890.30
      NM_010442Hmox1, Heme oxygenase 1I−0.780.55−1.220.04−0.650.18
      NM_010493Icam1, Intercellular adhesion molecule 1I−0.470.40−0.920.11−1.430.01
      NM_010517Igfbp4, Insulin like growth factor-binding protein 4S, C−1.400.03−0.900.12−0.990.14
      NM_021452Kcnmb4, Potassium calcium-activated channel subfamily M regulatory beta subunit 4S−1.560.04−2.410.04−0.220.64
      NM_008471Krt19, Keratin 19C−0.490.29−1.060.03−1.180.04
      NM_010875Ncam1, Neural cell adhesion molecule 1I, P−1.560.02−1.510.03−0.250.51
      NM_198424Orai3, ORAI calcium release-activated calcium modulator 3S−1.640.03−1.740.030.740.23
      NM_133200P2ry14, Purinergic receptor P2Y14S−0.710.21−0.760.24−1.430.04
      NM_001081549Rcan1, Regulator of calcineurin 1I, S−1.690.030.360.600.330.78
      NM_010101S1pr3, Sphingosine-1-phosphate receptor 3S−1.510.01−1.590.01−0.920.10
      NM_016769Smad3, SMAD family member 3I−1.320.01−1.250.02−0.230.47
      NM_011486Stat3, Signal transducer and activator of transcription 3I−1.600.04−1.060.06−0.850.17
      NM_009287Stim1, Stromal interaction molecule 1S−1.640.01−1.430.01−0.400.41

      Protein level changes in astrocytes after LI-rMS

      To determine whether changes in gene expression were reflected in changes in protein levels in astrocytes, 4 of the 21 downregulated genes were selected for validation by immunohistochemistry on cultured astrocytes (Figs. 2 and 3). These 4 genes were chosen based on the relevance of their functions to known changes in astrocyte biology following LI-rMS in vitro [
      • Clarke D.
      • Penrose M.A.
      • Penstone T.
      • Fuller-Carter P.I.
      • Hool L.C.
      • Harvey A.R.
      • et al.
      Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
      ] and in vivo [
      • Clarke D.
      • Penrose M.A.
      • Harvey A.R.
      • Rodger J.
      • Bates K.A.
      Low intensity rTMS has sex-dependent effects on the local response of glia following a penetrating cortical stab injury.
      ,
      • Monai H.
      • Ohkura M.
      • Tanaka M.
      • Oe Y.
      • Konno A.
      • Hirai H.
      • et al.
      Calcium imaging reveals glial involvement in transcranial direct current stimulation-induced plasticity in mouse brain.
      ]. The genes Stim1, Orai3, Kcnmb4, and Ncam1 had significantly downregulated gene transcripts in both 10Hz protocols but no significant regulation in the 1Hz group compared to sham. In the astrocytic store-operated calcium entry system relating to endoplasmic reticulum calcium stores, STIM1 protein functions as a detector of depleted calcium stores, and ORAI3 is a calcium re-entry channel [
      • Kwon J.
      • An H.
      • Sa M.
      • Won J.
      • Shin J.I.
      • Lee C.J.
      Orai1 and Orai3 in combination with Stim1 mediate the majority of store-operated calcium entry in astrocytes.
      ,
      • Moreno C.
      • Sampieri A.
      • Vivas O.
      • Peña-Segura C.
      • Vaca L.
      STIM1 and Orai1 mediate thrombin-induced Ca2+ influx in rat cortical astrocytes.
      ]. KCNMB4 is a beta subunit that modulates calcium sensitivity in voltage- and calcium-activated potassium (BK) channels [
      • Gebremedhin D.
      • Yamaura K.
      • Zhang C.
      • Bylund J.
      • Koehler R.C.
      • Harder D.R.
      Metabotropic glutamate receptor activation enhances the activities of two types of Ca2+-activated K+ channels in rat hippocampal astrocytes.
      ], and NCAM1/CD56 protein is involved in cell adhesion, proliferation and plasticity activities in astrocytes [
      • Jucker M.
      • Mondadori C.
      • Mohajeri H.
      • Bartsch U.
      • Schachner M.
      Transient upregulation of NCAM mRNA in astrocytes in response to entorhinal cortex lesions and ischemia.
      ,
      • Krushel L.A.
      • Tai M.-H.
      • Cunningham B.A.
      • Edelman G.M.
      • Crossin K.L.
      Neural cell adhesion molecule (N-CAM) domains and intracellular signaling pathways involved in the inhibition of astrocyte proliferation.
      ]. Preliminary studies confirmed that these proteins were expressed in astrocyte cultures, as shown in Supplementary Figure 1 for expression of KCNMB4 and NCAM1 in a GFAP-positive astrocyte. In addition, the intensity levels of the transcription factor EGR1 and the cell adhesion protein ICAM1 were explored but not pursued further as they exhibited a weak signal to noise ratio in immunofluorescent staining.
      Fig. 2
      Fig. 2Changes in protein intensity of KCNMB4 and NCAM1 at 5 h and 24 h after LI-rMS. A Representative images of KCNMB4 (green) and NCAM1 (red) co-labelled with Hoechst (blue) at 24 h after stimulation for each protocol. B,C Violin plots displaying the effect of stimulation on mean pixel intensity of KCNMB4 (B) and NCAM1 (C). Dotted lines indicate median and interquartile ranges and black bars between groups represent statistically significant differences between groups (p < 0.05). (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)
      Fig. 3
      Fig. 3Changes in protein intensity of STIM1 and ORAI3 at 5 h and 24 h after LI-rMS. A Representative images of STIM1 (green) and ORAI3 (red) co-labelled with Hoechst (blue) at 24 h after stimulation for each protocol. B,C Violin plots displaying the effect of stimulation on mean pixel intensity of STIM1 (B) and ORAI3 (C).Dotted lines indicate median and interquartile ranges and black bars between groups represent statistically significant differences between groups (p < 0.05). (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)
      For KCNMB4 protein immunolabel intensity (Fig. 2A), statistical analysis indicated significant differences dependent on astrocyte density (F(1,268) = 44.65, p < 0.01), location on the coverslip (F(1,268) = 18.30, p < 0.01), stimulation group (F(3,268) = 6.04, p < 0.01), and time post-stimulation (F(1,268) = 6.12, p = 0.014), but interactions between frequency and other independent variables were not significant (p > 0.05). Post-hoc analysis revealed no significant difference between stimulated cells and sham at 5h post-stimulation, but a reduced intensity of 10m10 Hz at 24h post-stimulation compared to sham (−39.4%, p < 0.001, Fig. 2B).
      For NCAM1 protein immunolabel intensity (Fig. 2A), there were significant differences dependent on astrocyte density (F(1,274) = 58.19, p < 0.01), stimulation group (F(3,260) = 12.30, p < 0.01), time post-stimulation (F(1,274) = 49.07, p < 0.01), and LI-rMS frequency∗time post-stimulation (F(3,274) = 7.05, p < 0.01), but not location on the coverslip (F(1, 274) = 0.53, p = 0.47). Specifically, mean NCAM1 immunolabel intensity was significantly lower following 10m1Hz compared to sham at 5h following stimulation (−42.1%, p < 0.001). By 24h, there were significantly reduced mean intensities for NCAM1 protein compared to sham for 1m10Hz (−38.3%, p < 0.001), 10m10Hz (−46.8%, p < 0.001), and 10m1Hz (−38.2%, p < 0.001, Fig. 2C).
      STIM1 immunolabel intensity (Fig. 3A) showed significant differences dependent on astrocyte density (F(1,253) = 20.11, p < 0.01), location of cells on the coverslip (F(1,253) = 4.88, p = 0.028), stimulation group (F(3,240) = 38.56, p < 0.01), time post-stimulation (F(1,257) = 29.05, p < 0.01), and frequency∗time (F(3,256) = 21.11, p < 0.01). Within the 5h timepoint, both 10m10Hz and 10m1Hz protocols led to significantly increased protein levels compared to sham (10m10Hz + 29.8%, p < 0.001; 10m1Hz +102.7%, p < 0.001). However by 24h, immunolabel intensity was significantly reduced after 1m10Hz (−35.3%, p < 0.001) and 10m10Hz protocols (−45.5%, p < 0.001) but significantly increased in the 1Hz group (+36.3%, p < 0.001, Fig. 3B).
      Similar effects were found for ORAI3 immunolabel intensity (Fig. 3A), where significant differences were dependent on location of cells on the coverslip (F(1,287) = 4.20, p = 0.041), stimulation group (F(3,231) = 30.25, p < 0.01), time post-stimulation (F(1,293) = 4.55, p = 0.034), and frequency∗time (F(3,290) = 7.27, p < 0.01), but not on astrocyte density (F(1,287) = 0.29, p = 0.86). Specifically, immunolabel intensity following 10m1Hz was significantly higher compared to sham at both time points (5h + 79.8%, 24h + 55.4%, p < 0.05 for both timepoints), and the two 10Hz protocols showed reduced intensity compared to sham at 24h (1m10Hz −41.6%, 1m10Hz −48.7%, p < 0.001 for both protocols, Fig. 3C).
      Overall, the density of astrocytes and/or their location on the coverslip were significant factors in the statistical analysis for each of the four proteins. Correlation analysis was performed to examine the relationship between protein intensity and astrocyte density, and protein intensity and magnetic field strength at the imaging locations within the coverslip (Supplementary Figure 2). Protein labelling intensity was positively correlated with astrocyte density. Although the correlations did reach significance for magnetic field and protein labelling intensity, the weak correlation coefficient suggest that differences in magnetic field intensity at the imaging locations had minimal effect on protein levels (Supplementary Tables 2–4).

      Discussion

      Many clinical and experimental rTMS studies have found changes to neuroplasticity [
      • Vlachos A.
      • Müller-Dahlhaus F.
      • Rosskopp J.
      • Lenz M.
      • Ziemann U.
      • Deller T.
      Repetitive magnetic stimulation induces functional and structural plasticity of excitatory postsynapses in mouse organotypic hippocampal slice cultures.
      ,
      • Ziemann U.
      • Hallett M.
      • Cohen L.G.
      Mechanisms of deafferentation-induced plasticity in human motor cortex.
      ], even at intensities below action potential threshold [
      • Robertson J.A.
      • Théberge J.
      • Weller J.
      • Drost D.J.
      • Prato F.S.
      • Thomas A.W.
      Low-frequency pulsed electromagnetic field exposure can alter neuroprocessing in humans.
      ,
      • Rodger J.
      • Mo C.
      • Wilks T.
      • Dunlop S.A.
      • Sherrard R.M.
      Transcranial pulsed magnetic field stimulation facilitates reorganization of abnormal neural circuits and corrects behavioral deficits without disrupting normal connectivity.
      ], suggesting classically non-excitable cells such as astrocytes may play a role in mediating LI-rTMS effects. This study shows that LI-rMS at 1 and 10 Hz results in significant reduction in the expression of a subset of genes in purified astrocyte cultures. The regulated genes are involved in various functional processes of astrocyte biology, primarily calcium signalling and inflammation. Importantly, the reduction in gene expression was associated with reduced expression of protein related to calcium signalling and cell adhesion. Taken together, these results suggest that LI-rMS has multiple and significant effects on astrocyte biology.

      Calcium signalling a key target of rMS in astrocytes

      Calcium signalling is emerging as a key mediator of rTMS effects [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ,
      • Lenz M.
      • Platschek S.
      • Priesemann V.
      • Becker D.
      • Willems L.M.
      • Ziemann U.
      • et al.
      Repetitive magnetic stimulation induces plasticity of excitatory postsynapses on proximal dendrites of cultured mouse CA1 pyramidal neurons.
      ,
      • Lenz M.
      • Galanis C.
      • Müller-Dahlhaus F.
      • Opitz A.
      • Wierenga C.J.
      • Szabó G.
      • et al.
      Repetitive magnetic stimulation induces plasticity of inhibitory synapses.
      ] and our results confirm downstream changes in the expression of genes related to regulation of intracellular calcium levels in astrocytes. In particular, 1Hz increased expression of STIM1 and ORAI3, proteins that are involved in sensing and restoring astrocytic endoplasmic reticulum calcium stores following depletion [
      • Kwon J.
      • An H.
      • Sa M.
      • Won J.
      • Shin J.I.
      • Lee C.J.
      Orai1 and Orai3 in combination with Stim1 mediate the majority of store-operated calcium entry in astrocytes.
      ,
      • Moreno C.
      • Sampieri A.
      • Vivas O.
      • Peña-Segura C.
      • Vaca L.
      STIM1 and Orai1 mediate thrombin-induced Ca2+ influx in rat cortical astrocytes.
      ]. The finding is consistent with our previous work showing that 1Hz LI-rMS increased astrocytic cytoplasmic calcium levels during stimulation in vitro [
      • Clarke D.
      • Penrose M.A.
      • Penstone T.
      • Fuller-Carter P.I.
      • Hool L.C.
      • Harvey A.R.
      • et al.
      Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
      ]. This suggests that the changes to the expression of calcium-related genes and proteins may be a consequence of acutely increased calcium levels during and/or after stimulation. However, we cannot rule out that LI-rMS might have affected gene transcriptional or translational machinery, perhaps by upregulating immediate early genes such as the transcription factor c-Fos [
      • Grehl S.
      • Martina D.
      • Goyenvalle C.
      • Deng Z.-D.
      • Rodger J.
      • Sherrard R.M.
      In vitro magnetic stimulation: a simple stimulation device to deliver defined low intensity electromagnetic fields.
      ], and/or by regulating the expression of microRNA [
      • Consales C.
      • Merla C.
      • Marino C.
      • Benassi B.
      The epigenetic component of the brain response to electromagnetic stimulation in Parkinson’s Disease patients: a literature overview.
      ] or chromatin organisation [
      • Etiévant A.
      • Manta S.
      • Latapy C.
      • Magno L.A.V.
      • Fecteau S.
      • Beaulieu J.-M.
      Repetitive transcranial magnetic stimulation induces long-lasting changes in protein expression and histone acetylation.
      ].
      Interestingly, the changes in gene expression reported here do not suggest changes to the genes directly involved in calcium release; our previous work showed that the increase in calcium during rMS is from intracellular stores because the increase persisted in calcium free media, but not when calcium stores in the endoplasmic reticulum were depleted by thapsigargin [
      • Clarke D.
      • Penrose M.A.
      • Penstone T.
      • Fuller-Carter P.I.
      • Hool L.C.
      • Harvey A.R.
      • et al.
      Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
      ]. The key mediator of calcium release from astrocyte intracellular stores is the protein inositol triphosphate receptor type 2 (IP3R2) [
      • Hua X.
      • Malarkey E.B.
      • Sunjara V.
      • Rosenwald S.E.
      • Li W.h.
      • Parpura V.
      Ca2+-dependent glutamate release involves two classes of endoplasmic reticulum Ca2+ stores in astrocytes.
      ], which is also necessary for increased astrocytic calcium release during transcranial direct current stimulation in vivo [
      • Monai H.
      • Ohkura M.
      • Tanaka M.
      • Oe Y.
      • Konno A.
      • Hirai H.
      • et al.
      Calcium imaging reveals glial involvement in transcranial direct current stimulation-induced plasticity in mouse brain.
      ]. However, in the present study, there was no significant transcriptional changes in Itpr2, the gene coding for IP3R2, suggesting that a single session of rMS may not cause long-term changes in calcium release properties, but rather impacts primarily on signalling downstream of calcium release.

      The timing of rTMS-induced gene expression

      In this study, changes to mRNA expression were examined at 5 h post LI-rMS, consistent with our previous study examining LI-rMS effects on mRNA expression of cultured neurons [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ]. However, there was no overlap of significantly regulated genes (primarily survival and apoptotic-related functions in neurons). Other studies investigating the effect of HI-rTMS on gene expression in vivo have reported changes as acutely as 45 min post-stimulation for the immediate early gene c-fos in multiple brain regions [
      • Ji R.-R.
      • Schlaepfer T.E.
      • Aizenman C.D.
      • Epstein C.M.
      • Qiu D.
      • Huang J.C.
      • et al.
      Repetitive transcranial magnetic stimulation activates specific regions in rat brain.
      ], and changes to excitatory and inhibitory neurotransmitter transporter and endoplasmic reticulum stress-related genes at 1 h following stimulation, which, for some genes, persisted at 4, 12, and 24 h post stimulation [
      • Ikeda T.
      • Kobayashi S.
      • Morimoto C.
      Effects of repetitive transcranial magnetic stimulation on ER stress-related genes and glutamate, γ-aminobutyric acid and glycine transporter genes in mouse brain.
      ]. Interestingly, Ikeda et al. found that the expression of certain genes varied above and below baseline in vivo and in PC12 cells in vitro during the 24 h following the end of HI-rTMS [
      • Ikeda T.
      • Kobayashi S.
      • Morimoto C.
      Effects of repetitive transcranial magnetic stimulation on ER stress-related genes and glutamate, γ-aminobutyric acid and glycine transporter genes in mouse brain.
      ]. These fluctuations over time likely reflect a dynamic response to the immediate gene expression changes induced by stimulation, as has been shown in cultured inhibitory and excitatory neurons following KCl depolarisation [
      • Spiegel I.
      • Mardinly Alan R.
      • Gabel Harrison W.
      • Bazinet Jeremy E.
      • Couch Cameron H.
      • Tzeng Christopher P.
      • et al.
      Npas4 regulates excitatory-inhibitory balance within neural circuits through cell-type-specific gene programs.
      ]. Recently, it was reported that the plasticity-related genes Ntrk2 and Mapk9 were significantly upregulated at 24 but not 6 h following HI-rMS in vitro in human SH-SY5Y cells with a neural phenotype [
      • Thomson A.C.
      • Kenis G.
      • Tielens S.
      • de Graaf T.A.
      • Schuhmann T.
      • Rutten B.P.F.
      • et al.
      Transcranial magnetic stimulation-induced plasticity mechanisms: TMS-related gene expression and morphology changes in a human neuron-like cell model.
      ]. These studies highlight the complex temporal relationship of neuronal gene expression in response to rTMS/rMS, which the current experiment now extends to astrocyte gene expression and protein levels. For example, there were changes to Kcnmb4 gene expression for 1m10Hz, but no change to protein levels, and significantly different protein levels of STIM1 and ORAI3 despite no significant changes in mRNA levels. This inconsistency between gene expression and protein levels after rMS is interesting and may reflect differences in protein stability and turn-over [
      • Trippe J.
      • Mix A.
      • Aydin-Abidin S.
      • Funke K.
      • Benali A.
      Theta burst and conventional low-frequency rTMS differentially affect GABAergic neurotransmission in the rat cortex.
      ] and/or indirect, downstream changes to protein synthesis. An extended time course of both mRNA and protein would be interesting to conduct to elucidate chronic effects of rMS on astrocyte gene expression and protein levels.

      Inflammation and rMS

      Amongst the four functional domains, injury and inflammatory genes were the most abundantly downregulated (10/21 genes). Changes to inflammatory pathways are likely to be associated with calcium signalling. Exposure to pro-inflammatory molecules decreases calcium signalling and the production of G protein-coupled receptors that regulate calcium signalling in astrocytes [
      • Hamby M.E.
      • Coppola G.
      • Ao Y.
      • Geschwind D.H.
      • Khakh B.S.
      • Sofroniew M.V.
      Inflammatory mediators alter the astrocyte transcriptome and calcium signaling elicited by multiple G-protein-coupled receptors.
      ,
      • Morita M.
      • Higuchi C.
      • Moto T.
      • Kozuka N.
      • Susuki J.
      • Itofusa R.
      • et al.
      Dual regulation of calcium oscillation in astrocytes by growth factors and pro-inflammatory cytokines via the mitogen-activated protein kinase cascade.
      ]. Additionally, increased calcium signalling has been reported in injury models that induce astrocyte reactivity [
      • Shigetomi E.
      • Saito K.
      • Sano F.
      • Koizumi S.
      Aberrant calcium signals in reactive astrocytes: a key process in neurological disorders.
      ]. Our finding of concomitantly altered inflammatory and calcium related genes following LI-rMS strengthens the relationship between these functions, but causality remains unclear. rTMS has previously been shown to reduce inflammation in stroke and injury animal models [
      • Sasso V.
      • Bisicchia E.
      • Latini L.
      • Ghiglieri V.
      • Cacace F.
      • Carola V.
      • et al.
      Repetitive transcranial magnetic stimulation reduces remote apoptotic cell death and inflammation after focal brain injury.
      ,
      • Hong Y.
      • Liu Q.
      • Peng M.
      • Bai M.
      • Li J.
      • Sun R.
      • et al.
      High-frequency repetitive transcranial magnetic stimulation improves functional recovery by inhibiting neurotoxic polarization of astrocytes in ischemic rats.
      ,
      • Caglayan A.B.
      • Beker M.C.
      • Caglayan B.
      • Yalcin E.
      • Caglayan A.
      • Yulug B.
      • et al.
      Acute and post-acute neuromodulation induces stroke recovery by promoting survival signaling, neurogenesis, and pyramidal tract plasticity.
      ], however modulation of inflammation in intact or normal models has not been observed. Although we did not aim to induce inflammation, the preparation of astrocyte cultures is a significant insult to the cells and is reported to stimulate inflammatory pathways [
      • Lange S.C.
      • Bak L.K.
      • Waagepetersen H.S.
      • Schousboe A.
      • Norenberg M.D.
      Primary cultures of astrocytes: their value in understanding astrocytes in health and disease.
      ]. Consistent with this, growing astrocytes in media-containing serum induces a gene expression profile that overlaps with profiles described for models of astrocyte reactivity [
      • Cahoy J.D.
      • Emery B.
      • Kaushal A.
      • Foo L.C.
      • Zamanian J.L.
      • Christopherson K.S.
      • et al.
      A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function.
      ,
      • Foo Lynette C.
      • Allen Nicola J.
      • Bushong Eric A.
      • Ventura P.B.
      • Chung W.-S.
      • Zhou L.
      • et al.
      Development of a method for the purification and culture of rodent astrocytes.
      ]. Our results show LI-rMS downregulates inflammatory-related genes, but given the artificial context of these cells in vitro, it is difficult to draw conclusions about possible anti-inflammatory effects of LI-rTMS in vivo. Nonetheless, this study supports a growing body of evidence that rTMS may reduce inflammation [
      • Sasso V.
      • Bisicchia E.
      • Latini L.
      • Ghiglieri V.
      • Cacace F.
      • Carola V.
      • et al.
      Repetitive transcranial magnetic stimulation reduces remote apoptotic cell death and inflammation after focal brain injury.
      ,
      • Hong Y.
      • Liu Q.
      • Peng M.
      • Bai M.
      • Li J.
      • Sun R.
      • et al.
      High-frequency repetitive transcranial magnetic stimulation improves functional recovery by inhibiting neurotoxic polarization of astrocytes in ischemic rats.
      ,
      • Caglayan A.B.
      • Beker M.C.
      • Caglayan B.
      • Yalcin E.
      • Caglayan A.
      • Yulug B.
      • et al.
      Acute and post-acute neuromodulation induces stroke recovery by promoting survival signaling, neurogenesis, and pyramidal tract plasticity.
      ].

      Significance of changes for astrocyte function

      In the rTMS literature, protocols are generally designated a binary category such as “facilitatory” and “inhibitory”, which refer to their expected effects on neuronal excitability. However, astrocytes do not communicate through action potential propagation, but instead use calcium signalling to transmit and process information that modulates neuronal activity and/or synaptic plasticity [
      • Scemes E.
      • Giaume C.
      Astrocyte calcium waves: what they are and what they do.
      ,
      • Singh A.
      • Abraham W.C.
      Astrocytes and synaptic plasticity in health and disease.
      ]. Furthermore, astrocytes are involved in other modulatory activities that do not directly alter neuron excitability but rather influence their inflammatory and injury response. Accordingly, even though pathway analysis did not yield any significant results, most likely due to the very wide range of disparate targets included in the array, the downregulated genes suggest reductions in pro-inflammatory activity, calcium signalling, cytoskeletal structure, plasticity, and cell adhesion properties. While there is clear delineation between the regulation of genes by each protocol, the main biological functions of the regulated genes are similarly spread across the protocols, complicating the distinction of potentially altered functional processes. We have previously demonstrated that rMS altered the expression of genes relating to apoptosis and cell survival in neuronal cultures [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ]. Therefore, the frequency-specific effects of rMS on neurons and astrocytes are likely to be broader than just facilitatory or inhibitory to brain activity and may include anti-inflammatory and pro-survival effects. Bidirectional neuron-glia interactions are firmly established [
      • Fields R.D.
      • Burnstock G.
      Purinergic signalling in neuron–glia interactions.
      ,
      • Theodosis D.T.
      • Poulain D.A.
      • Oliet S.H.R.
      Activity-dependent structural and functional plasticity of astrocyte-neuron interactions.
      ], and future in vitro studies using mixed cultures, potentially including a significant proportion of other cell types such as endothelial cells and immune cells, will contribute to understanding and optimising therapeutic effects of rTMS in conditions that involve inflammation such as stroke, depression and neurodegenerative conditions.

      The influence of astrocyte density and magnetic field on protein changes

      There was a non-uniform spatial distribution of astrocytes on coverslips, and immunohistochemical analysis revealed that astrocytes in regions of greater density were more intensely stained for KCNMB4, NCAM1, and STIM1 proteins. While such altered expression may reflect increased adhesion, gap junction coupling and intercellular communication via ion flow [
      • Theis M.
      • Giaume C.
      Connexin-based intercellular communication and astrocyte heterogeneity.
      ], there was no evidence that astrocytic density significantly influenced the observed LI-rMS effects on protein levels for 1Hz and 10Hz frequencies. An additional variable is the distribution of induced magnetic and electric fields is not uniform across the coverslip. The magnetic field is strongest under the centre of the coil, but the electric field is greatest underneath the coil windings [
      • Deng Z.-D.
      • Lisanby S.H.
      • Peterchev A.V.
      Electric field depth–focality tradeoff in transcranial magnetic stimulation: simulation comparison of 50 coil designs.
      ,
      • Wassermann E.
      • Epstein C.
      • Ziemann U.
      • Walsh V.
      • Paus T.
      • Lisanby S.
      Oxford handbook of transcranial stimulation.
      ]. In this study, the greatest magnetic field was near the centre of coverslip, whereas the greatest induced electric field was approximately halfway between the centre and edges of the coverslip (Supplementary Figure 2). There was a weak correlation between the location of astrocytes on the coverslip and the intensity of KCNMB4 and ORAI3. However, LI-rMS did not affect this relationship. This is unexpected, because of the differential intensity of the induced electric field across the coverslip. The reductions in protein levels at the centre of the coverslip might be a consequence of intercellular signalling propagating through the astrocyte network [
      • Venance L.
      • Piomelli D.
      • Glowinski J.
      • Glaume C.
      Inhibition by anandamide of gap junctions and intercellular calcium signalling in striatal astrocytes.
      ] from regions that were directly affected by the stronger electric field. In future studies, shielding part of the coverslip with mu-metal or applying gap junction blockers could address the role of astrocyte communication in propagating LI-rMS effects [
      • Dufor T.
      • Grehl S.
      • Tang A.D.
      • Doulazmi M.
      • Traoré M.
      • Debray N.
      • et al.
      Neural circuit repair by low-intensity magnetic stimulation requires cellular magnetoreceptors and specific stimulation patterns.
      ,
      • Stout C.E.
      • Costantin J.L.
      • Naus C.C.G.
      • Charles A.C.
      Intercellular calcium signaling in astrocytes via ATP release through connexin hemichannels.
      ].

      Frequency and pulse number

      Both frequency- and pulse number-dependent changes have been observed in neuronal gene and protein expression following rTMS in vivo and in vitro [
      • Grehl S.
      • Viola H.M.
      • Fuller-Carter P.I.
      • Carter K.W.
      • Dunlop S.A.
      • Hool L.C.
      • et al.
      Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
      ,
      • Trippe J.
      • Mix A.
      • Aydin-Abidin S.
      • Funke K.
      • Benali A.
      Theta burst and conventional low-frequency rTMS differentially affect GABAergic neurotransmission in the rat cortex.
      ,
      • Aydin-Abidin S.
      • Trippe J.
      • Funke K.
      • Eysel U.
      • Benali A.
      High- and low-frequency repetitive transcranial magnetic stimulation differentially activates c-Fos and zif268 protein expression in the rat brain.
      ,
      • Ljubisavljevic M.R.
      • Javid A.
      • Oommen J.
      • Parekh K.
      • Nagelkerke N.
      • Shehab S.
      • et al.
      The effects of different repetitive transcranial magnetic stimulation (rTMS) protocols on cortical gene expression in a rat model of cerebral ischemic-reperfusion injury.
      ,
      • Volz L.J.
      • Benali A.
      • Mix A.
      • Neubacher U.
      • Funke K.
      Dose-dependence of changes in cortical protein expression induced with repeated transcranial magnetic theta-burst stimulation in the rat.
      ], and astrocyte activities such as gliotransmission and calcium signalling are selective to frequency of neuronal activity [
      • Wang X.
      • Lou N.
      • Xu Q.
      • Tian G.-F.
      • Peng W.G.
      • Han X.
      • et al.
      Astrocytic Ca2+ signaling evoked by sensory stimulation in vivo.
      ,
      • Covelo A.
      • Araque A.
      Neuronal activity determines distinct gliotransmitter release from a single astrocyte.
      ]. In this study, there was an overlap of 6 genes that were regulated by stimulation delivered at 10Hz for 600 and 6000 pulses. However, there was no overlap in gene expression between pulse matched stimulation at 1Hz and 10Hz, and the majority (4/6) of significantly regulated calcium-related genes were following 10Hz LI-rMS regardless of duration. Moreover, differences in gene expression between 1m10Hz and 10m10Hz suggests an influence of duration (dose), which mainly affected inflammatory-related genes, with less pulses counterintuitively resulting in more genes being downregulated. These results suggest that frequency is a significant determinant of LI-rMS effects on astrocytes. However, further characterisation of common frequencies with duration and pulse numbers controls is required to determine the frequency- and duration-specific (dose) effects of stimulation on astrocytes. It is important to note that dose-response in brain stimulation and pharmacological contexts differ, with the latter considering concepts such as threshold, slope, and maximal asymptote of the treatment [
      • Blumenthal D.K.
      ].

      Conclusion

      LI-rTMS has the ability to induce neuroplasticity, with a wide range of effects on neuronal plasticity. However, the effect of rTMS on astrocytes is not well understood. Our data shows that a single stimulation session of 1Hz or 10Hz can reduce the expression of genes and proteins in purified astrocyte cell cultures. The functions of the affected genes and proteins primarily relate to calcium signalling and inflammation, and suggest that LI-rTMS may have therapeutic effects that contribute to LI-rTMS induced neuroplasticity.

      Author declaration form

      We wish to confirm that there are no known conflicts of interest associated with this publication and there has been no significant financial support for this work that could have influenced its outcome.
      We confirm that the manuscript has been read and approved by all named authors and that there are no other persons who satisfied the criteria for authorship but are not listed. We further confirm that the order of authors listed in the manuscript has been approved by all of us.
      We confirm that we have given due consideration to the protection of intellectual property associated with this work and that there are no impediments to publication, including the timing of publication, with respect to intellectual property. In so doing we confirm that we have followed the regulations of our institutions concerning intellectual property.
      We further confirm that any aspect of the work covered in this manuscript that involved experimental animals has been conducted with the ethical approval of all relevant bodies and that such approvals are acknowledged within the manuscript.

      CRediT authorship contribution statement

      Darren Clarke: Conceptualization, Methodology, Formal analysis, Investigation, Writing - original draft, Writing - review & editing, Visualization. Jamie Beros: Investigation. Kristyn A. Bates: Conceptualization, Methodology, Resources, Writing - review & editing, Supervision, Funding acquisition. Alan R. Harvey: Conceptualization, Methodology, Resources, Writing - review & editing, Supervision, Funding acquisition. Alexander D. Tang: Methodology, Software, Formal analysis, Writing - review & editing, Supervision. Jennifer Rodger: Conceptualization, Methodology, Resources, Writing - review & editing, Supervision, Project administration, Funding acquisition.

      Declaration of competing interest

      The authors declare no competing interest.

      Acknowledgements

      This work was supported by the Neurotrauma Research Program.

      Appendix A. Supplementary data

      References

        • Lenz M.
        • Müller-Dahlhaus F.
        • Vlachos A.
        Cellular and molecular mechanisms of rTMS-induced neural plasticity.
        in: Platz T. Therapeutic rTMS in neurology: principles, evidence, and Practice recommendations. Springer International Publishing, Cham2016: 11-22
        • Tang A.D.
        • Thickbroom G.
        • Rodger J.
        Repetitive transcranial magnetic stimulation of the brain: mechanisms from animal and experimental models.
        Neuroscientist. 2015; 23: 82-94
        • Thickbroom G.W.
        Transcranial magnetic stimulation and synaptic plasticity: experimental framework and human models.
        Exp Brain Res. 2007; 180: 583-593
        • Martiny K.
        • Lunde M.
        • Bech P.
        Transcranial low voltage pulsed electromagnetic fields in patients with treatment-resistant depression.
        Biol Psychiatr. 2010; 68: 163-169
        • Robertson J.A.
        • Théberge J.
        • Weller J.
        • Drost D.J.
        • Prato F.S.
        • Thomas A.W.
        Low-frequency pulsed electromagnetic field exposure can alter neuroprocessing in humans.
        J R Soc Interface. 2010; 7: 467-473
        • Makowiecki K.
        • Harvey A.R.
        • Sherrard R.M.
        • Rodger J.
        Low-intensity repetitive transcranial magnetic stimulation improves abnormal visual cortical circuit topography and upregulates BDNF in mice.
        J Neurosci. 2014; 34: 10780-10792
        • Rodger J.
        • Mo C.
        • Wilks T.
        • Dunlop S.A.
        • Sherrard R.M.
        Transcranial pulsed magnetic field stimulation facilitates reorganization of abnormal neural circuits and corrects behavioral deficits without disrupting normal connectivity.
        Faseb J. 2012; 26: 1593-1606
        • Tang A.D.
        • Bennett W.
        • Hadrill C.
        • Collins J.
        • Fulopova B.
        • Wills K.
        • et al.
        Low intensity repetitive transcranial magnetic stimulation modulates skilled motor learning in adult mice.
        Sci Rep. 2018; 8: 4016
        • Tang A.D.
        • Hong I.
        • Boddington L.J.
        • Garrett A.R.
        • Etherington S.
        • Reynolds J.N.
        • et al.
        Low-intensity repetitive magnetic stimulation lowers action potential threshold and increases spike firing in layer 5 pyramidal neurons in vitro.
        Neuroscience. 2016; 335: 64-71
        • Grehl S.
        • Viola H.M.
        • Fuller-Carter P.I.
        • Carter K.W.
        • Dunlop S.A.
        • Hool L.C.
        • et al.
        Cellular and molecular changes to cortical neurons following low intensity repetitive magnetic stimulation at different frequencies.
        Brain Stimulation. 2015; 8: 114-123
        • Seewoo B.J.
        • Feindel K.W.
        • Etherington S.J.
        • Rodger J.
        Frequency-specific effects of low-intensity rTMS can persist for up to 2 weeks post-stimulation: a longitudinal rs-fMRI/MRS study in rats.
        Brain Stimulation. 2019; 12: 1526-1536
        • Stöckel T.
        • Summers J.J.
        • Hinder M.R.
        Reversed effects of intermittent theta burst stimulation following motor training that vary as a function of training-induced changes in corticospinal excitability.
        Neural Plast. 2015; 5 (2015)
        • Teo J.T.H.
        • Swayne O.B.C.
        • Cheeran B.
        • Greenwood R.J.
        • Rothwell J.C.
        Human theta burst stimulation enhances subsequent motor learning and increases performance variability.
        Cerebr Cortex. 2010; 21: 1627-1638
        • Kang J.
        • Jiang L.
        • Goldman S.A.
        • Nedergaard M.
        Astrocyte-mediated potentiation of inhibitory synaptic transmission.
        Nat Neurosci. 1998; 1: 683
        • Papouin T.
        • Dunphy J.
        • Tolman M.
        • Foley J.C.
        • Haydon P.G.
        Astrocytic control of synaptic function.
        Phil Trans Biol Sci. 2017; 37220160154
        • Pascual O.
        • Casper K.B.
        • Kubera C.
        • Zhang J.
        • Revilla-Sanchez R.
        • Sul J.-Y.
        • et al.
        Astrocytic purinergic signaling coordinates synaptic networks.
        Science. 2005; 310: 113-116
        • Perea G.
        • Araque A.
        Astrocytes potentiate transmitter release at single hippocampal synapses.
        Science. 2007; 317: 1083-1086
        • Serrano A.
        • Haddjeri N.
        • Lacaille J.-C.
        • Robitaille R.
        GABAergic network activation of glial cells underlies hippocampal heterosynaptic depression.
        J Neurosci. 2006; 26: 5370-5382
        • Cullen C.L.
        • Young K.M.
        How does transcranial magnetic stimulation influence glial cells in the central nervous system?.
        Front Neural Circ. 2016; 10
        • Gellner A.-K.
        • Reis J.
        • Fritsch B.
        Glia: a neglected player in non-invasive direct current brain stimulation.
        Front Cell Neurosci. 2016; 10
        • Clarke D.
        • Penrose M.A.
        • Penstone T.
        • Fuller-Carter P.I.
        • Hool L.C.
        • Harvey A.R.
        • et al.
        Frequency-specific effects of repetitive magnetic stimulation on primary astrocyte cultures.
        Restor Neurol Neurosci. 2017; 35: 557-569
        • Cacace F.
        • Mineo D.
        • Viscomi M.T.
        • Latagliata E.C.
        • Mancini M.
        • Sasso V.
        • et al.
        Intermittent theta-burst stimulation rescues dopamine-dependent corticostriatal synaptic plasticity and motor behavior in experimental parkinsonism: possible role of glial activity.
        Mov Disord. 2017; 32: 1035-1046
        • Clarke D.
        • Penrose M.A.
        • Harvey A.R.
        • Rodger J.
        • Bates K.A.
        Low intensity rTMS has sex-dependent effects on the local response of glia following a penetrating cortical stab injury.
        Exp Neurol. 2017; 295: 233-242
        • Rauš S.
        • Selaković V.
        • Manojlović-Stojanoski M.
        • Radenović L.
        • Prolić Z.
        • Janać B.
        Response of hippocampal neurons and glial cells to alternating magnetic field in gerbils submitted to global cerebral ischemia.
        Neurotox Res. 2013; 23: 79-91
        • Sasso V.
        • Bisicchia E.
        • Latini L.
        • Ghiglieri V.
        • Cacace F.
        • Carola V.
        • et al.
        Repetitive transcranial magnetic stimulation reduces remote apoptotic cell death and inflammation after focal brain injury.
        J Neuroinflammation. 2016; 13: 1-6
        • Chan P.
        • Eng L.F.
        • Ling Lee Y.
        • Lin V.W.
        Effects of pulsed magnetic stimulation on GFAP levels in cultured astrocytes.
        J Neurosci Res. 1999; 55: 238-244
        • Hong Y.
        • Liu Q.
        • Peng M.
        • Bai M.
        • Li J.
        • Sun R.
        • et al.
        High-frequency repetitive transcranial magnetic stimulation improves functional recovery by inhibiting neurotoxic polarization of astrocytes in ischemic rats.
        J Neuroinflammation. 2020; 17: 150
        • Peng Z.w.
        • Xue F.
        • Zhou C.h.
        • Zhang R.g.
        • Wang Y.
        • Liu L.
        • et al.
        Repetitive transcranial magnetic stimulation inhibits Sirt1/MAO-A signaling in the prefrontal cortex in a rat model of depression and cortex-derived astrocytes.
        Mol Cell Biochem. 2018; 442: 59-72
        • Pfaffl M.W.
        A new mathematical model for relative quantification in real-time RT–PCR.
        Nucleic Acids Res. 2001; 29: e45-e
        • Pfaffl M.W.
        • Tichopad A.
        • Prgomet C.
        • Neuvians T.P.
        Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper – Excel-based tool using pair-wise correlations.
        Biotechnol Lett. 2004; 26: 509-515
        • Liao Y.
        • Wang J.
        • Jaehnig E.J.
        • Shi Z.
        • Zhang B.
        WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs.
        Nucleic Acids Res. 2019; 47: W199-W205
        • Monai H.
        • Ohkura M.
        • Tanaka M.
        • Oe Y.
        • Konno A.
        • Hirai H.
        • et al.
        Calcium imaging reveals glial involvement in transcranial direct current stimulation-induced plasticity in mouse brain.
        Nat Commun. 2016; 7
        • Kwon J.
        • An H.
        • Sa M.
        • Won J.
        • Shin J.I.
        • Lee C.J.
        Orai1 and Orai3 in combination with Stim1 mediate the majority of store-operated calcium entry in astrocytes.
        Exp Neurobiol. 2017; 26: 42-54
        • Moreno C.
        • Sampieri A.
        • Vivas O.
        • Peña-Segura C.
        • Vaca L.
        STIM1 and Orai1 mediate thrombin-induced Ca2+ influx in rat cortical astrocytes.
        Cell Calcium. 2012; 52: 457-467
        • Gebremedhin D.
        • Yamaura K.
        • Zhang C.
        • Bylund J.
        • Koehler R.C.
        • Harder D.R.
        Metabotropic glutamate receptor activation enhances the activities of two types of Ca2+-activated K+ channels in rat hippocampal astrocytes.
        J Neurosci. 2003; 23: 1678-1687
        • Jucker M.
        • Mondadori C.
        • Mohajeri H.
        • Bartsch U.
        • Schachner M.
        Transient upregulation of NCAM mRNA in astrocytes in response to entorhinal cortex lesions and ischemia.
        Mol Brain Res. 1995; 28: 149-156
        • Krushel L.A.
        • Tai M.-H.
        • Cunningham B.A.
        • Edelman G.M.
        • Crossin K.L.
        Neural cell adhesion molecule (N-CAM) domains and intracellular signaling pathways involved in the inhibition of astrocyte proliferation.
        Proc Natl Acad Sci Unit States Am. 1998; 95: 2592-2596
        • Vlachos A.
        • Müller-Dahlhaus F.
        • Rosskopp J.
        • Lenz M.
        • Ziemann U.
        • Deller T.
        Repetitive magnetic stimulation induces functional and structural plasticity of excitatory postsynapses in mouse organotypic hippocampal slice cultures.
        J Neurosci. 2012; 32: 17514-17523
        • Ziemann U.
        • Hallett M.
        • Cohen L.G.
        Mechanisms of deafferentation-induced plasticity in human motor cortex.
        J Neurosci. 1998; 18: 7000-7007
        • Lenz M.
        • Platschek S.
        • Priesemann V.
        • Becker D.
        • Willems L.M.
        • Ziemann U.
        • et al.
        Repetitive magnetic stimulation induces plasticity of excitatory postsynapses on proximal dendrites of cultured mouse CA1 pyramidal neurons.
        Brain Struct Funct. 2015; 220: 3323-3337
        • Lenz M.
        • Galanis C.
        • Müller-Dahlhaus F.
        • Opitz A.
        • Wierenga C.J.
        • Szabó G.
        • et al.
        Repetitive magnetic stimulation induces plasticity of inhibitory synapses.
        Nat Commun. 2016; 7: 10020
        • Grehl S.
        • Martina D.
        • Goyenvalle C.
        • Deng Z.-D.
        • Rodger J.
        • Sherrard R.M.
        In vitro magnetic stimulation: a simple stimulation device to deliver defined low intensity electromagnetic fields.
        Front Neural Circ. 2016; 10
        • Consales C.
        • Merla C.
        • Marino C.
        • Benassi B.
        The epigenetic component of the brain response to electromagnetic stimulation in Parkinson’s Disease patients: a literature overview.
        Bioelectromagnetics. 2018; 39: 3-14
        • Etiévant A.
        • Manta S.
        • Latapy C.
        • Magno L.A.V.
        • Fecteau S.
        • Beaulieu J.-M.
        Repetitive transcranial magnetic stimulation induces long-lasting changes in protein expression and histone acetylation.
        Sci Rep. 2015; 5: 16873
        • Hua X.
        • Malarkey E.B.
        • Sunjara V.
        • Rosenwald S.E.
        • Li W.h.
        • Parpura V.
        Ca2+-dependent glutamate release involves two classes of endoplasmic reticulum Ca2+ stores in astrocytes.
        J Neurosci Res. 2004; 76: 86-97
        • Ji R.-R.
        • Schlaepfer T.E.
        • Aizenman C.D.
        • Epstein C.M.
        • Qiu D.
        • Huang J.C.
        • et al.
        Repetitive transcranial magnetic stimulation activates specific regions in rat brain.
        Proc Natl Acad Sci Unit States Am. 1998; 95: 15635-15640
        • Ikeda T.
        • Kobayashi S.
        • Morimoto C.
        Effects of repetitive transcranial magnetic stimulation on ER stress-related genes and glutamate, γ-aminobutyric acid and glycine transporter genes in mouse brain.
        Biochemistry and Biophysics Reports. 2019; 17: 10-16
        • Spiegel I.
        • Mardinly Alan R.
        • Gabel Harrison W.
        • Bazinet Jeremy E.
        • Couch Cameron H.
        • Tzeng Christopher P.
        • et al.
        Npas4 regulates excitatory-inhibitory balance within neural circuits through cell-type-specific gene programs.
        Cell. 2014; 157: 1216-1229
        • Thomson A.C.
        • Kenis G.
        • Tielens S.
        • de Graaf T.A.
        • Schuhmann T.
        • Rutten B.P.F.
        • et al.
        Transcranial magnetic stimulation-induced plasticity mechanisms: TMS-related gene expression and morphology changes in a human neuron-like cell model.
        Front Mol Neurosci. 2020; 13
        • Trippe J.
        • Mix A.
        • Aydin-Abidin S.
        • Funke K.
        • Benali A.
        Theta burst and conventional low-frequency rTMS differentially affect GABAergic neurotransmission in the rat cortex.
        Exp Brain Res. 2009; 199: 411
        • Hamby M.E.
        • Coppola G.
        • Ao Y.
        • Geschwind D.H.
        • Khakh B.S.
        • Sofroniew M.V.
        Inflammatory mediators alter the astrocyte transcriptome and calcium signaling elicited by multiple G-protein-coupled receptors.
        J Neurosci. 2012; 32: 14489-14510
        • Morita M.
        • Higuchi C.
        • Moto T.
        • Kozuka N.
        • Susuki J.
        • Itofusa R.
        • et al.
        Dual regulation of calcium oscillation in astrocytes by growth factors and pro-inflammatory cytokines via the mitogen-activated protein kinase cascade.
        J Neurosci. 2003; 23: 10944-10952
        • Shigetomi E.
        • Saito K.
        • Sano F.
        • Koizumi S.
        Aberrant calcium signals in reactive astrocytes: a key process in neurological disorders.
        Int J Mol Sci. 2019; 20: 996
        • Caglayan A.B.
        • Beker M.C.
        • Caglayan B.
        • Yalcin E.
        • Caglayan A.
        • Yulug B.
        • et al.
        Acute and post-acute neuromodulation induces stroke recovery by promoting survival signaling, neurogenesis, and pyramidal tract plasticity.
        Front Cell Neurosci. 2019; 13
        • Lange S.C.
        • Bak L.K.
        • Waagepetersen H.S.
        • Schousboe A.
        • Norenberg M.D.
        Primary cultures of astrocytes: their value in understanding astrocytes in health and disease.
        Neurochem Res. 2012; 37: 2569-2588
        • Cahoy J.D.
        • Emery B.
        • Kaushal A.
        • Foo L.C.
        • Zamanian J.L.
        • Christopherson K.S.
        • et al.
        A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function.
        J Neurosci. 2008; 28: 264-278
        • Foo Lynette C.
        • Allen Nicola J.
        • Bushong Eric A.
        • Ventura P.B.
        • Chung W.-S.
        • Zhou L.
        • et al.
        Development of a method for the purification and culture of rodent astrocytes.
        Neuron. 2011; 71: 799-811
        • Scemes E.
        • Giaume C.
        Astrocyte calcium waves: what they are and what they do.
        Glia. 2006; 54: 716-725
        • Singh A.
        • Abraham W.C.
        Astrocytes and synaptic plasticity in health and disease.
        Exp Brain Res. 2017; 235: 1645-1655
        • Fields R.D.
        • Burnstock G.
        Purinergic signalling in neuron–glia interactions.
        Nat Rev Neurosci. 2006; 7: 423-436
        • Theodosis D.T.
        • Poulain D.A.
        • Oliet S.H.R.
        Activity-dependent structural and functional plasticity of astrocyte-neuron interactions.
        Physiol Rev. 2008; 88: 983-1008
        • Theis M.
        • Giaume C.
        Connexin-based intercellular communication and astrocyte heterogeneity.
        Brain Res. 2012; 1487: 88-98
        • Deng Z.-D.
        • Lisanby S.H.
        • Peterchev A.V.
        Electric field depth–focality tradeoff in transcranial magnetic stimulation: simulation comparison of 50 coil designs.
        Brain Stimulation. 2013; 6: 1-13
        • Wassermann E.
        • Epstein C.
        • Ziemann U.
        • Walsh V.
        • Paus T.
        • Lisanby S.
        Oxford handbook of transcranial stimulation.
        Oxford University Press, 2008
        • Venance L.
        • Piomelli D.
        • Glowinski J.
        • Glaume C.
        Inhibition by anandamide of gap junctions and intercellular calcium signalling in striatal astrocytes.
        Nature. 1995; 376: 590-594
        • Dufor T.
        • Grehl S.
        • Tang A.D.
        • Doulazmi M.
        • Traoré M.
        • Debray N.
        • et al.
        Neural circuit repair by low-intensity magnetic stimulation requires cellular magnetoreceptors and specific stimulation patterns.
        Science Advances. 2019; 5eaav9847
        • Stout C.E.
        • Costantin J.L.
        • Naus C.C.G.
        • Charles A.C.
        Intercellular calcium signaling in astrocytes via ATP release through connexin hemichannels.
        J Biol Chem. 2002; 277: 10482-10488
        • Aydin-Abidin S.
        • Trippe J.
        • Funke K.
        • Eysel U.
        • Benali A.
        High- and low-frequency repetitive transcranial magnetic stimulation differentially activates c-Fos and zif268 protein expression in the rat brain.
        Exp Brain Res. 2008; 188: 249-261
        • Ljubisavljevic M.R.
        • Javid A.
        • Oommen J.
        • Parekh K.
        • Nagelkerke N.
        • Shehab S.
        • et al.
        The effects of different repetitive transcranial magnetic stimulation (rTMS) protocols on cortical gene expression in a rat model of cerebral ischemic-reperfusion injury.
        PloS One. 2015; 10e0139892
        • Volz L.J.
        • Benali A.
        • Mix A.
        • Neubacher U.
        • Funke K.
        Dose-dependence of changes in cortical protein expression induced with repeated transcranial magnetic theta-burst stimulation in the rat.
        Brain Stimulation. 2013; 6: 598-606
        • Wang X.
        • Lou N.
        • Xu Q.
        • Tian G.-F.
        • Peng W.G.
        • Han X.
        • et al.
        Astrocytic Ca2+ signaling evoked by sensory stimulation in vivo.
        Nat Neurosci. 2006; 9: 816
        • Covelo A.
        • Araque A.
        Neuronal activity determines distinct gliotransmitter release from a single astrocyte.
        eLife. 2018; 7e32237
        • Blumenthal D.K.
        B L.L. H-D R. K B.C. Goodman & gilman’s: the pharmacological basis of therapeutics. 13 ed. McGraw-Hill, New York2018