Abstract
Repetitive transcranial magnetic stimulation (rTMS) is increasingly used as a treatment for neurological and psychiatric disorders. Although the induced field is focused on a target region during rTMS, adjacent areas also receive stimulation at a lower intensity and the contribution of this perifocal stimulation to network-wide effects is poorly defined. Here, we examined low-intensity rTMS (LI-rTMS)-induced changes on a model neural network using the visual systems of normal (C57Bl/6J wild-type, n = 22) and ephrin-A2A5−/− (n = 22) mice, the latter possessing visuotopic anomalies. Mice were treated with LI-rTMS or sham (handling control) daily for 14 d, then fluorojade and fluororuby were injected into visual cortex. The distribution of dorsal LGN (dLGN) neurons and corticotectal terminal zones (TZs) was mapped and disorder defined by comparing their actual location with that predicted by injection sites. In the afferent geniculocortical projection, LI-rTMS decreased the abnormally high dispersion of retrogradely labeled neurons in the dLGN of ephrin-A2A5−/− mice, indicating geniculocortical map refinement. In the corticotectal efferents, LI-rTMS improved topography of the most abnormal TZs in ephrin-A2A5−/− mice without altering topographically normal TZs. To investigate a possible molecular mechanism for LI-rTMS-induced structural plasticity, we measured brain derived neurotrophic factor (BDNF) in the visual cortex and superior colliculus after single and multiple stimulations. BDNF was upregulated after a single stimulation for all groups, but only sustained in the superior colliculus of ephrin-A2A5−/− mice. Our results show that LI-rTMS upregulates BDNF, promoting a plastic environment conducive to beneficial reorganization of abnormal cortical circuits, information that has important implications for clinical rTMS.
Introduction
Activity-dependent mechanisms are necessary for maintaining neural circuits throughout life (Butz et al., 2009) and promoting plasticity during learning and memory (Martin et al., 2000). However, the ability to reorganize the structure and function of neural circuits is largely restricted to a critical period of development. Therefore, adults with neurological and psychiatric disorders, many of which are associated with neural circuit and plasticity abnormalities (Thickbroom and Mastaglia, 2009; Arnsten and Rubia, 2012), have limited capacity for repair.
Noninvasive brain stimulation techniques such as repetitive transcranial magnetic stimulation (rTMS) are thought to facilitate neural plasticity and demonstrate clinical promise in treating many neurological and psychiatric disorders (Wassermann and Zimmermann, 2012). rTMS uses trains of magnetic pulses to induce electric currents in the brain that modulate cortical activity (Pascual-Leone, 2006) and can produce persistent excitability and behavioral modifications (Maeda et al., 2000; Valero-Cabré et al., 2005). Although these changes are generally attributed to induced neuronal spiking within the stimulated region (Zangen et al., 2005), associated activation of the wider network also contributes to rTMS-induced functional plasticity (Lee et al., 2003; Rounis et al., 2006). This is consistent with rTMS modulating metabolism (Valero-Cabré et al., 2005; Valero-Cabré et al., 2007) and gene expression (Aydin-Abidin et al., 2008) in regions downstream of the target area.
However, the induced field is not restricted to the target region of rTMS: perifocal areas also receive stimulation, albeit at lower intensities, below the neuron firing threshold (Wagner et al., 2009). Although the functional importance of this perifocal subthreshold stimulation is unclear, it may contribute to network-wide effects. Support for this hypothesis comes from human studies which show that low-intensity rTMS (LI-rTMS) modifies cortical excitability (Cook et al., 2004; Capone et al., 2009), alters pain perception (Robertson et al., 2010), and alleviates depression symptoms (Martiny et al., 2010). Moreover, animal and in vitro studies demonstrate that low-intensity magnetic stimulation alters calcium signaling (Pessina et al., 2001), gene expression (Mattsson and Simkó, 2012), neurotransmitter synthesis (Sieroń et al., 2004), and the structure and function of a single neural pathway (Rodger et al., 2012). Together, these studies suggest that LI-rTMS, like suprathreshold stimulation (Valero-Cabré et al., 2005; Valero-Cabré et al., 2008; Reithler et al., 2011), may have network-wide effects.
In the present study, we evaluate to what extent LI-rTMS induces structural reorganization in afferent and efferent pathways within a network. We use a well-defined model, the mouse visual system, which has topographic projections between multiple centers. We examined LI-rTMS effects on the afferent geniculocortical [dorsal LGN (dLGN) to primary visual cortex (V1)] and the efferent corticotectal (V1 to superior colliculus) projections in both normal (wild-type) and abnormal (ephrin-A2A5−/−) adult mice, which have ectopic terminals in addition to normally located projections (Haustead et al., 2008; Wilks et al., 2010). To examine a possible mechanism underlying structural reorganization, we measured BDNF expression in the visual cortex and superior colliculus of wild-type and ephrin-A2A5−/− mice after a single stimulation or after multiple LI-rTMS sessions.
Materials and Methods
Subjects and housing.
Ephrin-A2A5−/− mice were originally a generous gift from David Feldheim and carry a homozygous-null mutation of the ephrin-A2 and A5 genes (Feldheim et al., 2000). Ephrin-A2A5−/− mice were bred from heterozygous parents and were backcrossed for >10 generations on a C57Bl/6J background. Mice were genotyped at weaning, as described previously (Feldheim et al., 2000). Wild-type mice of the same genetic background were purchased from the Animal Research Centre of Murdoch University (Perth, Western Australia). Mice were bred in the Biomedical Research Facility of The University of Western Australia and housed in single-sex cages on a 12 h light/dark cycle with food and water provided ad libitum. Anatomical and BDNF data were obtained from separate cohorts because of tissue-processing differences and to avoid the potential confounder of minor injury from the tracer injection altering BDNF expression (Hicks et al., 1997). Mice in both cohorts were mature young adults (2–7 months old; Grillo et al., 2013) and we confirmed there was no significant age effect on any measure (all p > 0.05). All procedures in this study were in accordance with National Institutes of Health guidelines and were approved by The University of Western Australia Animal Ethics Committee.
LI-rTMS administration.
For anatomical experiments, mice received either LI-rTMS or sham (no stimulation, handling control) for 10 min daily for 14 d. Magnetic pulses were generated using an electromagnetic pulse generator (Global Energy Medicine) connected to mains power by a standard 9 V AC/DC adaptor. The pulse generator was wired to a 3.5 mm jack (Jaycar) for attachment of a custom-built circle coil suitable for use on mice (Fig. 1A,B). The coil (copper wire, 0.25 mm diameter, 16Ω, 300 windings; Jaycar) had an 8 mm outer diameter and 6 mm inner diameter and was filled with a steel bolt to increase field penetration (Epstein and Davey, 2002). The coil was designed to achieve a similar coil-to-brain size ratio as used in humans, enabling focal stimulation of the head (Weissman et al., 1992). However, the magnetic field produced by our custom-built coils obviously differs from standardized devices currently used in human TMS/rTMS research (Lang et al., 2006; Ziemann et al., 2008; Lefaucheur, 2009).
Using a Hall device to measure magnetic field strength (Rodger et al., 2012), the coil generated a 12 mT magnetic field at its base, a programmed voltage input pulse width of 220 μs, and measured field rise time of 230 μs, which did not vary within the 10 min LI-rTMS session or between sessions. The field at the visual cortex was estimated to be ∼8 mT (Fig. 1C). This intensity is within the range known to affect cellular biochemistry (Morgado-Valle et al., 1998; Arias-Carrión et al., 2004), yet is below motor threshold, obviating the need for anesthetic or restraint, with their associated changes to neuronal excitability and circulating stress hormones (Gersner et al., 2011).
We delivered a high-frequency complex pattern of stimulation comprising 59.9 ms trains of 20 pulses, with trains repeated at 6.67 Hz for 1 min, 10.01 Hz for 8 min, and 6.25 Hz for the final minute (Rodger et al., 2012). This pattern was developed from electrobiomimetic principles to replicate endogenous patterns of electrical fields around activated nerves during exercise (Martiny et al., 2010) and is based on the patent PCT/AU2007/000454 of Global Energy Medicine. This pattern and intensity is consistent with the concept of a subthreshold “priming” stimulation at one frequency modulating cortical responsiveness to subsequent stimulation, as previously demonstrated using different combinations of stimuli (Iyer et al., 2003; Hamada et al., 2008).
To ensure that mice remained relatively still without restraint during LI-rTMS, they were habituated to handling and the coil for 3 d before the study, as described previously (Rodger et al., 2012). During stimulation, the coil was held between the ears, allowing reproducible targeting of visual cortices bilaterally (Fig. 1B). It was held above, but not in contact with, the head to ensure that any potential vibration of the coil could not confound the sham-stimulation procedure, in which the stimulator was switched off. After 10 min of stimulation, coil temperature was 26.5°C at the coil surface and 25.5°C 1 mm from the thermometer. Because this is within the range of room temperature (22–26°C), confounding effects of heat are unlikely. We also verified that the active coil did not emit audible or ultrasonic sound. The sound output emitted by the coil was measured in a sound proofed room using a 1/2-inch Bruel and Kjaer Type 4134 condensor microphone and a Hewlett Packard 3561A Dynamic Signal Analyzer. The microphone was positioned approximately the same distance from the coil as the animal's ear. Absolute sound intensity was calibrated using a Bruel and Kjaer Type 4231 sound calibrator. The maximum acoustic energy emitted was at 1.3 kHz and 37 dB SPL. At 8 kHz, the sound intensity had fallen to 10 dB SPL, which was only 7 dB above background. Above 10 kHz, there was no detectable acoustic energy emitted above background, which averaged 9 dB SPL from 10 to 20 kHz (measured with the LI-rTMS device turned off). This sound intensity (<7 dB in the frequency range detectable by mice) is below the auditory threshold (30 dB) for wild-type (Fernandez et al., 2010) and ephrin-A2A5−/− mice (Yates et al., 2014). Therefore, consistent with the low stimulation intensity and lack of auditory or heat/tactile confounders, mice did not display any gross motor, visuomotor, or behavioral response to either the stimulation or sham procedure.
Anatomical tracing of cortical projections.
On the 15th day after beginning sham or LI-rTMS, mice were anesthetized with xylazine (10 mg/kg) and ketamine (70 mg/kg). Mice were placed in a stereotaxic frame and a small part of the skull and dura were removed to expose the left primary visual cortex (V1), which was identified based on stereotaxic coordinates (Mouse Brain Atlas, Franklin and Paxinos, 2008) and confirmed visually using a landmark surface branch of the middle cerebral artery (Dorr et al., 2007). Injections into V1 were between 0.05 mm caudal and 1.5 mm rostral to lambda and 1.5 to 3 mm from the midline. Each mouse received two injections of fluorescent tracers, biotinylated dextran amine (BDA; 10,000 MW; Life Technologies): Alexa Fluor 488 (fluorojade) for the lateral injection and Alexa Fluor 555 (fluororuby) for the medial injection. Tracers were pressure injected using a nanojet injector fitted with a micropipette attached to a micromanipulator. The micropipette was lowered to 400 μm from the cortex surface to target layer 5 neurons, which project to the superior colliculus, and layer 4 geniculocortical axon terminals. Six 50 nL pulses of tracer were injected at each site.
Four days after cortical injections, mice were killed by intraperitoneal injection of 0.1 ml of sodium pentobarbitone (Lethabarb; Virbac) and transcardially perfused with 80 ml of 4% paraformaldehyde in 0.2 m phosphate buffer. The brain was dissected and postfixed for 24 h at 4°C before cryoprotection in graded concentrations of sucrose (15%, 30%) in PBS, which were changed after 24 h.
Brains were cryosectioned coronally (40 μm) and three parallel series of sections were collected onto gelatin-subbed slides. One of the series was used to visualize fluorescent tracers. Sections were dried on slides for 5 min before immersion in PBS for 2 min to remove optimal cutting temperature (OCT) embedding medium, dried for a further 5 min, and coverslipped with Dakomount mounting medium (Dako). Sections were examined using a Nikon fluorescent microscope and photographed using a cooled digital camera (QuantiFIRE; Optronics). The adjacent second series was stored at −80°C and selected slides were processed for fluorescent microscopy to confirm terminal zones (TZs) when cases in series 1 were ambiguous. Images of fluorescent sections were analyzed using Adobe Photoshop (CS6 version 13.1).
The third series was Nissl stained to measure dLGN volume and to confirm the location of injection sites within V1. Slides were air-dried, immersed in distilled water for 2 min (to remove OCT), heated to 50–60°C in ethanol with 5% glacial acetic acid solution, and then cooled to room temperature (2 h). Sections were rehydrated in distilled water (2 min), stained in prewarmed (55°C) cresyl violet acetate (0.5%; Proscitech) for 8 min, and then differentiated in ethanol with 5% glacial acetic acid for up to 13 min before dehydration in ethanol (3 × 3 min) and xylene (2 × 2 min), and mounted with Entellan (Merck). Slides with Nissl-stained sections were scanned for image analysis using an Aperio ScanScope and analyzed using Aperio ImageScope software (version 10.2.2.2352).
Analyses of anatomical tracing.
All analyses of anatomical tracing were performed blinded to stimulation condition. It was not possible to blind to genotype as well, because the presence or absence of additional TZs made the genotype apparent. We confirmed that injections were at an appropriate depth in V1 by examining fluorescent sections cross-checked against adjacent Nissl-stained sections and regions defined in the mouse brain atlas (Franklin and Paxinos, 2008). We tested for group differences in injection size by thresholding images to remove background (Haustead et al., 2008) and measuring the labeled area for each injection site in sections from a single series. To obtain injection volume (in cubic micrometers), area measures were summed and then multiplied by section thickness (40 μm) and the number of series. We applied the same method to obtain volumes of corticotectal TZs and dLGN measures.
We determined topography of the corticotectal projection by reconstructing TZ locations in the superior colliculus using every third section. Multiple TZs were scored if TZs were separated by at least one section in the examined series (i.e., >120 μm in the rostrocaudal axis) or visually distinct from an adjacent TZ (mediolateral axis; Wilks et al., 2010). There is an established relationship between cortical injection site on the mediolateral axis and TZ location in superior collicular rostrocaudal axis, with lateral injections resulting in rostral labeling and medial injections resulting in caudal labeling (Triplett et al., 2009; Wilks et al., 2010). Accordingly, we examined topography by analyzing injection site location mediolaterally as a function of TZ location rostrocaudally. Injection site location mediolaterally was defined as the distance from approximate injection center-point to the midline as a percentage of the total cortical hemisphere width (measured using Photoshop). TZ location was defined as the distance from the midpoint (rostrocaudally) of the TZ to the caudal extreme of the superior colliculus as a percentage of the total superior colliculus length in the rostrocaudal axis measured by multiplying number of sections by inter-section distance.
In a separate analysis, we distinguished between topographically normal and abnormal TZs within each mouse, using relative position in the azimuth axis (degrees) of a Cartesian plane for the medial (red) to lateral (green) injection sites and corresponding red-green TZs pairs. We restricted this analysis to mice in which both injections successfully labeled TZs. The difference from the predicted relative location was determined for each TZ pair using residuals from linear regression constrained to a perfect-fit (intercept = 0, slope = 1), with injection azimuth as the predictor and TZ azimuth as the outcome variable. Within each mouse, TZ pairs were ranked to determine the least disordered (closest to predicted location) and the most disordered (farthest from predicted location). We obtained TZ volumes for each distinct TZ (in cubic micrometers) as described for injection volume. Volumes were summed for the same TZ pairs used in the assessment of topographical disorder. Because all wild-type mice had only a single TZ per injection, data were included in analyses for both least and most disordered TZs. Two ephrin-A2A5−/− LI-rTMS mice had only a single TZ per injection and were included in both the least and most disordered TZ analyses because comparison between TZ pairs was not possible.
To assess LI-rTMS effects on the geniculocortical projection, we quantified the spatial distribution of retrogradely labeled neurons in the dLGN measured as cell dispersion and labeled cluster size. Cell dispersion was measured as the area defined by straight lines joining the outermost labeled neurons (convex-hull) in each section, converted to volume (in cubic micrometers). To control for dLGN size differences between mice, convex-hull volume was normalized to total dLGN volume obtained by measuring dLGN area in each Nissl-stained section (one series) and multiplying by inter-section distance. We also measured the mean area of the main labeled cluster from three representative sections for each animal. We averaged the area from three sections within the middle range of labeled dLGN sections instead of total cluster volume because labeling was too sparsely distributed at the rostral and caudal extremes, so cluster boundaries could not be identified accurately. The number of retrogradely labeled neurons was estimated in one series (one of every three sections in total), in accordance with stereological principles (West, 1993).
For all dLGN measures, preliminary analyses confirmed no significant differences in neuronal dLGN labeling between medial and lateral injections (data not shown). Therefore, data describing retrogradely labeled neurons from medial and lateral injections were averaged within mice rather than being analyzed separately. Two subjects were deemed outliers due to scores >3 SDs from the mean and were excluded.
BDNF expression.
Using another cohort of mice of the same age range, we evaluated a potential mechanism underlying LI-rTMS by measuring BDNF in the visual cortex and, to confirm previous results (Rodger et al., 2012), the superior colliculus. BDNF expression in the dLGN could not be examined because of the impossibility of precise microdissection of this region. Axons of the optic tract traverse the external surface of the dLGN and are known to transport BDNF (Quigley et al., 2000; Spalding et al., 2002); therefore, any results obtained could not be attributed exclusively to the dLGN. To compare the acute and cumulative effects of LI-rTMS, BDNF levels were measured in mice killed 2 and 24 h after a single stimulation session and 24 h after the last session of 14 d stimulation (n = 8). Mice were habituated to handling and the LI-rTMS coil for 4 d before commencing stimulation or sham. All other aspects of LI-rTMS administration were identical to that used in the anatomical experiments. All biochemical procedures were performed blinded to stimulation condition and genotype.
After terminal anesthesia of 0.2 ml of sodium pentobarbitone injected intraperitoneally, visual cortex and superior colliculus samples were rapidly excised and frozen. Samples were stored at −80°C until homogenized in 1 ml of ice-cold lysis buffer (100 mm PIPES, pH 7; 500 mm NaCl; 0.2% Triton X-100; 0.1% NaN3; 2 mm EDTA; EDTA-free mini protease inhibitors; Roche Diagnostics; Szapacs et al., 2004) and centrifuged (3320 × g for 1 h at 4°C). Protein content was determined using the BCA method (BCA protein assay kit; Pierce). The supernatant was analyzed by ELISA for BDNF according to the manufacturer's instructions (ChemiKine BDNF Sandwich ELISA; Millipore Bioscience Research Reagents). Optical absorbance measures (450 nm) were performed automatically using a spectrophotometer plate reader with EnSpire software (PerkinElmer). BDNF concentration was then calculated as a percentage of total protein for each sample. For comparison between groups, data were normalized to percentage of wild-type-sham within each region.
Statistical analyses.
The number of corticotectal TZs per injection was analyzed with Fisher's exact test. Topography was assessed further with linear regression. All other anatomical tracing data were analyzed using two-factor ANOVAs with between-subjects factors: genotype (wild-type, ephrin-A2A5−/−) and stimulation condition (sham, LI-rTMS). Follow-up post hoc tests with Sidak correction were constrained to comparisons between sham and LI-rTMS within each genotype. Because BDNF expression may not be independent between the visual cortex and superior colliculus within each mouse, group differences combined across the two regions were first analyzed using multivariate ANOVA (MANOVA) to protect against elevated risk of type one error, and then followed up with ANOVAs and Sidak-corrected post hoc tests performed separately for each region. The relationship between BDNF concentration in the visual cortex and superior colliculus was examined using Pearson's correlations and linear regression. The level for significance was set at p < 0.05. Statistical analyses were performed using GraphPad Prism version 6 and SPSS version 20 software.
Results
To assess LI-rTMS effects on cortical circuit reorganization, we examined topography of V1 projections and geniculate afferents using anterograde (corticotectal) and retrograde (geniculocortical) tracing (injection outcomes summarized in Table 1). Mean injection volumes did not differ significantly between groups (ANOVA, all p > 0.05), allowing direct comparison of measurements between experimental groups.
LI-rTMS improves topography of abnormal corticotectal projections in ephrin-A2A5−/−mice
First, we examined the effect of LI-rTMS on cortical efferent projections by analyzing the number of TZs labeled in the corticotectal projection of sham- and LI-rTMS-treated wild-type and ephrin-A2A5−/− mice. We confirmed that successful V1 injections resulted in a single TZ in wild-type mice and two TZs in the majority of ephrin-A2A5−/− mice (Table 1), indicative of abnormal topography (Fig. 2; Wilks et al., 2010). However, LI-rTMS did not reduce the number of ectopic TZs in ephrin-A5A5−/− mice (Fisher's exact test, p = 0.99).
To further investigate LI-rTMS effects on corticotectal topography, we mapped the location of corticotectal TZs within the superior colliculus relative to injection location using linear regression analysis (Fig. 3). Consistent with previous results (Triplett et al., 2009; Wilks et al., 2010), wild-types showed a strong, positive relationship between injection site (percentage medial-lateral) and TZ location (percentage rostral-caudal; F(1,34) = 9.31, p = 0.004; Fig. 3A). However, in ephrin-A2A5−/− mice, when all successful injections and distinct TZs were included in the analysis, this relationship was weak and nonsignificant in both the sham and LI-rTMS groups, suggesting abnormal topography (Fig. 3B). In either genotype, sham and LI-rTMS groups did not separately account for significantly more variance in the regression model than a global best-fit line, showing that there was no significant difference between LI-rTMS and sham (wild-types: F(2,32) = 0.23, p = 0.79; ephrin-A2A5−/−: F(2,54) = 0.25, p = 0.78).
Because normal and abnormal projections are present in ephrin-A2A5−/− mice and LI-rTMS may affect them differentially, we used regression analysis to identify the least and most disordered TZ in each mouse based on distance from location predicted by injection sites (see Materials and Methods). For the least disordered TZs, there were no significant differences between genotypes or stimulation conditions in mean distance from predicted location (ANOVA; all p > 0.05; Fig. 4A). In contrast, the most disordered TZs in ephrin-A2A5−/− mice differed significantly from wild-types in distance from predicted location (F(1,26) = 17.25, p < 0.001). However, there was no significant effect of stimulation condition (F(1,26) = 2.64, p = 0.12) and no significant interaction (F(1,26) = 3.17, p = 0.09). Follow-up tests were conducted despite the nonsignificant main effect of stimulation because, given previous findings in the retinotectal projection (Rodger et al., 2012), we hypothesized that LI-rTMS would affect ephrin-A2A5−/− mice, but not wild-types. Follow-up tests with Sidak correction showed that LI-rTMS significantly decreased disorder compared with sham in the most disordered TZs of ephrin-A2A5−/− mice (p = 0.04), but did not affect topography in wild-type mice (p = 0.99). In summary, although LI-rTMS did not decrease the number of ectopic corticotectal TZs, LI-rTMS improved topographic disorder by shifting the most abnormal corticotectal terminals toward a more topographically appropriate location.
All LI-rTMS groups showed slightly greater mean TZ volumes than sham (Fig. 4B), however, this was not significant (all p > 0.05). Consistent with previous reports (Wilks et al., 2010), wild-types had significantly greater mean TZ volumes than the most disordered ephrin-A2A5−/− TZs (F(1,26) = 7.97, p = 0.01), but not the least disordered TZs (F(1,26) = 2.14, p = 0.16).
LI-rTMS refines abnormal geniculocortical topography in ephrin-A2A5−/− mice
In addition to the structural plasticity of cortical efferent projections, we also evaluated LI-rTMS effects on V1 afferent terminals by quantifying the distribution and number of retrogradely labeled dLGN neurons (Fig. 5).
Consistent with previous reports (Cang et al., 2008), our convex-hull measurements showed that labeled neurons were scattered over a significantly higher percentage of the dLGN in ephrin-A2A5−/− mice compared with wild-types, representing greater dispersion (F(1,38) = 30.90, p < 0.001; Fig. 5E). In addition, LI-rTMS significantly reduced the abnormally high dispersion in ephrin-A2A5−/− mice (p < 0.001), but had no significant effect in wild-types (p = 0.95). Moreover, the mean area occupied by the main cluster, excluding scattered labeled neurons, remained constant, with no significant differences between groups (ANOVA; all p > 0.05; Fig. 5F). Therefore, our data suggest that LI-rTMS selectively affected topographically inappropriate geniculocortical afferents.
Further supporting LI-rTMS-induced reorganization of abnormal geniculocortical afferents in ephrin-A2A5−/− mice, there were significant effects of stimulation condition and genotype on the number of labeled dLGN neurons (stimulation condition: F(1,38) = 6.01, p = 0.01; genotype: F(1,38) = 9.27, p = 0.004; Fig. 5G). Although post hoc tests were nonsignificant, there was a strong trend for LI-rTMS to reduce the number of labeled neurons compared with sham-treated ephrin-A2A5−/− mice (Sidak correction: p = 0.07). This trend was not observed for wild-types (p = 0.40).
LI-rTMS increases BDNF concentration
To investigate a possible molecular mechanism for LI-rTMS-induced structural reorganization, we compared BDNF concentrations in the visual cortex and superior colliculus between sham- and LI-rTMS-treated mice. To distinguish between acute and cumulative effects, we measured BDNF after a single stimulation session (2 and 24 h after stimulation) and 14 daily sessions. Consistent with previous studies (Kuczewski et al., 2008), BDNF levels were in picogram amounts. Analysis was conducted relative to total protein in the sample and relative to wild-type sham controls (see Materials and Methods).
Because BDNF is anterogradely transported from visual cortex to superior colliculus (Conner et al., 1997), its expression in each region cannot be considered independent. Therefore, we used MANOVA (Pillai's trace) to analyze group differences combined across the two regions, which showed a significant effect of LI-rTMS time point on BDNF concentrations in the visual cortex and superior colliculus (V = 0.96, F(6,112) = 17.16, p < 0.001). Subsequent two-factor ANOVAs, conducted separately for the visual cortex and superior colliculus, showed that BDNF concentration differed significantly between genotypes (visual cortex: F(1,56) = 7.87, p = 0.007; superior colliculus: F(1,56) = 31.75, p < 0.0001) and LI-rTMS time points (visual cortex: F(3,56) = 20.79, p < 0.001; superior colliculus: F(3,56) = 20.68, p < 0.0001). There was also a significant interaction in both regions (visual cortex: F(3,56) = 3.69, p = 0.02; superior colliculus: F(3,56) = 10.49, p < 0.0001), indicating that wild-type and ephrin-A2A5−/− mice responded differently at specific LI-rTMS time points. In wild-types, LI-rTMS acutely increased BDNF levels in both the visual cortex (Fig. 6A) and the superior colliculus (Fig. 6B) compared with sham (2 h, p < 0.001; 24 h, p < 0.01), but these increases were not sustained after 14 d of stimulation (p = 0.99). In contrast, in ephrin-A2A5−/− mice, BDNF expression increased significantly in the visual cortex only at 24 h after a single stimulation session (p = 0.02) and was not sustained after 14 daily stimulation sessions (p = 0.19), whereas in the superior colliculus of ephrin-A2A5−/− mice, BDNF expression was significantly greater compared with sham for all LI-rTMS time points (all p <0.0001; Fig. 6B).
Finally we tested whether the raised collicular BDNF might arise from cortical synthesis and transport (i.e., BDNF produced in visual cortex may be transported to the superior colliculus, resulting in lower levels in visual cortex and correspondingly higher levels in the superior colliculus) through correlation analyses between visual cortex and superior colliculus BDNF concentrations at each time point (Fig. 6C–F). There were no significant correlations (all p > 0.05) and most correlations were weak, suggesting that transport between regions is unlikely to contribute to BDNF concentration as strongly as local production and/or transport from other regions (e.g., retina; Conner et al., 1997; Quigley et al., 2000; Spalding et al., 2002).
Discussion
We used both normal and abnormal mouse visual systems to investigate the effects of low-intensity magnetic stimulation on neural circuits. We present the first direct evidence that LI-rTMS induces structural reorganization across multiple relays within the visual network (Fig. 7). Our data suggest that differential BDNF signaling in visual cortex and superior colliculus may underlie LI-rTMS reorganization of abnormal neural circuits. Furthermore, the low intensity of our stimulation suggests that induced spiking activity may not be a prerequisite for magnetic stimulation effects. Together, our results suggest that LI-rTMS, and possibly rTMS clinically, may reinitiate some aspects of the developmental critical period that are permissive for neuroanatomical reorganization.
LI-rTMS affects multiple relays of the visual pathway: parallels with development
Human studies demonstrate that rTMS effects extend beyond the targeted region and affect interconnected pathways (Reithler et al., 2011). Consistent with this, we observed reorganization after LI-rTMS in three interconnected visual pathways: the retinotectal (Rodger et al., 2012), corticotectal, and geniculocortical projections. Moreover, the changes that we observed in these three circuits are reminiscent of visual system development (McLaughlin and O'Leary, 2005). Therefore, our finding that LI-rTMS induces structural reorganization in adult mice that are beyond the classical critical period for visual system map formation (Inoue et al., 1992) suggests that LI-rTMS may reactivate developmental processes (Cline, 2003; McLaughlin et al., 2003; McLaughlin and O'Leary, 2005).
A likely mechanism for such plasticity is the upregulation of BDNF induced by magnetic stimulation (present study; Gersner et al., 2011; Rodger et al., 2012). During development, closure of the critical period is associated with downregulation of BDNF (Huang et al., 1999). Furthermore, increased BDNF within the mature visual system, whether induced pharmacologically or behaviorally through environmental enrichment, reactivates critical period-like plasticity and reorganization of visual cortical afferents (Mandolesi et al., 2005; Sale et al., 2007; Maya Vetencourt et al., 2008; Baroncelli et al., 2010). Although it is unclear how strongly BDNF contributes to behavioral improvements in the visual system or other brain regions, dark-reared mice displayed normal visual cortical topography after injection with BDNF, suggesting that BDNF alone may be sufficient to reorganize visual cortical circuits during development (Gianfranceschi et al., 2003). However, work in adult systems suggests that, although BDNF may be a prerequisite for plastic change (Berardi et al., 1999), functionally useful reorganization may also require coordinated neural activity that is achievable with behavioral training (Sale et al., 2007; Sale et al., 2009). Therefore, LI-rTMS, in conjunction with behavioral training, could be a useful adjuvant therapy for rehabilitation.
In addition, reactivation of developmental processes could also facilitate interactions between retinotectal and corticotectal maps. In the developing visual system, retinotectal topography is established first within the superior colliculus, and the later developing corticotectal map partly depends on the retinotectal map for its organization (Triplett et al., 2009). It is therefore possible that LI-rTMS-induced reorganization within the retinotectal projection (Rodger et al., 2012) may contribute to, or even drive, the shift of corticotectal terminals. Therefore, developmental mechanisms, reminiscent of the critical period, may underlie LI-rTMS-induced visual pathway reorganization.
LI-rTMS selectively affects abnormal projections
Within the visual network, reorganization by LI-rTMS is not uniform; our data show that LI-rTMS preferentially induces plasticity in abnormal cortical circuits. This finding is consistent with animal studies (Valero-Cabré et al., 2008; Rodger et al., 2012; Sykes et al., 2013) and human studies showing null or transient stimulation effects in healthy participants, but longer-term effects in clinical populations (Siebner et al., 2004; Capone et al., 2009; Martiny et al., 2010). Disrupted homeostatic plasticity in abnormal systems may interact with BDNF signaling (Desai et al., 1999), enabling LI-rTMS to affect abnormal circuits selectively. Subthreshold changes to membrane potential, such as those induced by LI-rTMS, upregulate BDNF expression (Zafra et al., 1990; Singh and Miller, 2005; Pell et al., 2011) and may contribute to the initial increase in BDNF for all groups. However, in normal systems, such excitation also triggers homeostatic mechanisms (e.g., GABA-A receptor activation) that downregulate BDNF (Berninger et al., 1995), consistent with wild-types showing no change in BDNF after 14 d of stimulation. In contrast, in ephrin-A2A5−/− mice, competition between ectopic and appropriately located terminals is likely to perturb homeostatic plasticity (Cao et al., 2007). Therefore, in the superior colliculus, where cortical efferent ectopic terminals are retained, albeit more accurately located, the stimulus for upregulation of BDNF is sustained. In contrast, in the visual cortex, geniculocortical errors in ephrin-A2A5−/− mice are less severe than corticotectal abnormalities (Wilks et al., 2010) and were effectively corrected by LI-rTMS; therefore, the transient BDNF increase in the visual cortex, which returned to baseline by 14 d, may reflect normal homeostatic downregulation of BDNF (Huang et al., 1999).
BDNF also modulates growth cone response to axonal guidance cues, including semaphorins and Eph/ephrins (Fitzgerald et al., 2008; Marler et al., 2008; Marler et al., 2010; Naska et al., 2010); therefore, LI-rTMS-induced upregulation of BDNF may promote appropriate circuit reorganization by increasing sensitivity to environmental cues. Our findings linking LI-rTMS-induced structural reorganization to BDNF expression are compatible with long-term behavioral modifications in healthy (Valero-Cabré et al., 2008) and parietally injured (Afifi et al., 2013) cats. In addition, our findings suggest a mechanism to explain the clinical efficacy of low-intensity magnetic stimulation, as well as rTMS, in disorders associated with plasticity and BDNF abnormalities (Thickbroom and Mastaglia, 2009; Autry and Monteggia, 2012; low intensity: Prato et al., 2005; Martiny et al., 2010; rTMS: Wassermann and Zimmermann, 2012).
Projection-specific effects of LI-rTMS may reflect multiple mechanisms of action
In addition to acting primarily on the abnormal visual circuits, a key result of our study is that LI-rTMS corrects topographical errors differently within specific components of the visual circuit, suggesting different mechanisms. A possible explanation is that BDNF has different effects on axonal terminals depending whether signaling occurs via TrkB or p75NTR/Troy (Schecterson and Bothwell, 2010), suggesting that the differential expression of these receptors throughout the visual circuit may determine whether abnormal projections are moved or eliminated. In contrast to our previous work in the retinotectal pathway of ephrin-A2A5−/− mice, where LI-rTMS removed ectopic terminals (Rodger et al., 2012), here, we show that abnormal corticotectal projections are retained but moved to more topographically appropriate locations. Although LI-rTMS upregulated BDNF in the superior colliculus of ephrin-A2A5−/− mice, where both projections converge (Sefton et al., 2004), the different outcomes in these pathways may be due to regional differences in neurotrophin receptor expression. In the retinotectal projection, high p75NTR/Troy in retinal ganglion cells (Ugolini et al., 1995; Park et al., 2005), together with low TrkB expression (Rickman and Brecha, 1995), would promote ectopic terminal elimination (Singh et al., 2008; Sherrard et al., 2009). In contrast, corticotectal axons coexpress TrkB and p75NTR (McQuillen et al., 2002), allowing the higher affinity of BDNF for TrkB to promote reorganization, rather than removal, of ectopic corticotectal terminals (Bibel et al., 1999; Majdan and Miller, 1999). We also show plasticity in the geniculocortical projection of ephrin-A2A5−/− mice after transient upregulation of BDNF. Geniculocortical neurons express TrkB without p75NTR (Silver and Stryker, 2001; McQuillen et al., 2002); withdrawal of TrkB signaling after a period of transient activation results in axon retraction (Luikart and Parada, 2006; Sherrard et al., 2009). This suggests that the reduced dLGN cell dispersion we observed after LI-rTMS is due to removal, rather than reorganization, of aberrant geniculocortical terminals. Our present results do not allow causal relationships to be established between BDNF signaling and LI-rTMS-induced reorganization across the different relays of the visual pathway; future research could address this interesting issue through lesion studies and analysis of behavioral performance.
Projection-specific effects may also be due to differences in magnetic field intensities delivered across the brain (Cohen et al., 1990; Deng et al., 2013). Visual brain regions received stimulation synchronously but at different intensities: V1 and superior colliculus are in close proximity in the mouse brain (Figs. 1, 7) and received similar intensities (6–8 mT), but the dLGN received ∼2 mT and the retina <1 mT (Rodger et al., 2012). Given that human studies show different immediate excitability changes depending on stimulation intensity (Kähkönen et al., 2005; Batsikadze et al., 2013), these differences in field strength across the visual pathway could contribute to the different patterns of reorganization after LI-rTMS: terminal elimination within the weakly stimulated retinotectal (Rodger et al., 2012) and geniculocortical projections, but reorganization within the more strongly stimulated corticotectal projection.
Potential mechanisms underlying effects of low-intensity magnetic stimulation
Our LI-rTMS stimulation was subthreshold, estimated to be 0.5 V/m (Martiny et al., 2010), which is several orders of magnitude below the ∼50 V/m that depolarizes neurons (Volz et al., 2013). Therefore, our data show for the first time that circuit reorganization can occur in the absence of induced spiking activity, which is thought to be the primary mechanism underlying rTMS network effects (Thickbroom, 2007; Benali et al., 2011; Volz et al., 2013). Such network changes in the absence of induced neuronal firing is consistent with the frequency-specific changes to EEG waves reported for very-low-intensity magnetic stimulation (Di Lazzaro et al., 2013). Although mechanisms remain poorly understood, modeling studies indicate that stimulation at intensities used here can alter ion transfer across cell membranes, causing subthreshold changes to membrane potential (Poo, 1981). These subthreshold changes can have lasting effects on neuron activity by opening ion channels, modifying excitability and calcium signaling (Yu et al., 2010; Grehl et al., 2011). In addition, subthreshold presynaptic inputs modulate postsynaptic activity (Alle and Geiger, 2006; Shu et al., 2006). Therefore, transsynaptic modulation may underlie the plastic reorganization that we observed across multiple relays within the visual pathway after LI-rTMS.
Our data suggest mechanisms that may underpin the therapeutic effects of very low-intensity fields described in human studies (Capone et al., 2009; Martiny et al., 2010; Di Lazzaro et al., 2013). In addition, because we show that low-intensity pulsed magnetic fields of similar intensity to those delivered outside of the focal area of stimulation in human rTMS (Wagner et al., 2009) affects neuronal circuits, our study may have implications for understanding mechanisms underlying rTMS clinical effects. However, because our device delivers very different stimulation parameters compared with those typically used in human rTMS/TMS (Wagner et al., 2009; Pell et al., 2011), our data should not be extrapolated directly to a clinical context.
Conclusions
In summary, we show that multiple sessions of low-intensity, complex waveform magnetic stimulation induced structural reorganization, improving topographical accuracy in abnormal cortical circuits. Our results suggest that LI-rTMS may prime neuronal circuits via upregulation of BDNF, maintaining a plastic environment conducive to beneficial reorganization. Therefore, this study reveals the power of applying magnetic field stimulation, even at low intensities, to abnormal or dysfunctional neural circuits to correct structural and functional errors.
Footnotes
This work was supported by the National Health and Medical Research Council (NHMRC) Australia (project grant), the Australian Research Council (linkage grant), and the Neurotrauma Research Program (project grant to J.R., funded through the Road Trauma Trust Account, but does not reflect views or recommendations of the Road Safety Council). K.M. receives doctoral scholarship support from The University of Western Australia and an Australian Postgraduate Award. J.R. is an NHMRC Senior Research Fellow. We thank Andrew Garrett, Mark Hinder, and Gary Thickbroom for useful comments on the manuscript; Don Robertson for advice and assistance with coil sound measurements; Joao Paulo Coimbra for histology advice; and Marissa Penrose-Menz for technical assistance.
The authors declare no competing financial interests.
- Correspondence should be addressed to Jennifer Rodger, Experimental and Regenerative Neuroscience M317, School of Animal Biology, University of Western Australia, 35 Stirling Highway, Crawley, Western Australia 6009. jennifer.rodger{at}uwa.edu.au