Abstract
Mitochondria play a variety of functional roles in cortical neurons, from metabolic support and neuroprotection to the release of cytokines that trigger apoptosis. In dendrites, mitochondrial structure is closely linked to their function, and fragmentation (fission) of the normally elongated mitochondria indicates loss of their function under pathological conditions, such as stroke and brain trauma. Using in vivo two-photon microscopy in mouse brain, we quantified mitochondrial fragmentation in a full spectrum of cortical injuries, ranging from severe to mild. Severe global ischemic injury was induced by bilateral common carotid artery occlusion, whereas severe focal stroke injury was induced by Rose Bengal photosensitization. The moderate and mild traumatic injury was inflicted by focal laser lesion and by mild photo-damage, respectively. Dendritic and mitochondrial structural changes were tracked longitudinally using transgenic mice expressing fluorescent proteins localized either in cytosol or in mitochondrial matrix. In response to severe injury, mitochondrial fragmentation developed in parallel with dendritic damage signified by dendritic beading. Reconstruction from serial section electron microscopy confirmed mitochondrial fragmentation. Unlike dendritic beading, fragmentation spread beyond the injury core in focal stroke and focal laser lesion models. In moderate and mild injury, mitochondrial fragmentation was reversible with full recovery of structural integrity after 1–2 weeks. The transient fragmentation observed in the mild photo-damage model was associated with changes in dendritic spine density without any signs of dendritic damage. Our findings indicate that alterations in neuronal mitochondria structure are very sensitive to the tissue damage and can be reversible in ischemic and traumatic injuries.
SIGNIFICANCE STATEMENT During ischemic stroke or brain trauma, mitochondria can either protect neurons by supplying ATP and adsorbing excessive Ca2+, or kill neurons by releasing proapoptotic factors. Mitochondrial function is tightly linked to their morphology: healthy mitochondria are thin and long; dysfunctional mitochondria are thick (swollen) and short (fragmented). To date, fragmentation of mitochondria was studied either in dissociated cultured neurons or in brain slices, but not in the intact living brain. Using real-time in vivo two-photon microscopy, we quantified mitochondrial fragmentation during acute pathological conditions that mimic severe, moderate, and mild brain injury. We demonstrated that alterations in neuronal mitochondria structural integrity can be reversible in traumatic and ischemic injuries, highlighting mitochondria as a potential target for therapeutic interventions.
Introduction
In addition to supplying cellular energy, neuronal mitochondria are involved in the biosynthesis of neurotransmitters and the regulation of intracellular Ca2+ signaling, acting as a “sink” and “store” of Ca2+, and thereby are important for synaptic communication and neuronal plasticity (Chan, 2006; Mattson et al., 2008). Traumatic and ischemic brain injuries are multifactorial disorders that lead to a spectrum of pathological changes, having a significant influence on the patients' health status. Although the primary impact differs, the mechanisms of injury progression can be similar to some extent (Kunz et al., 2010; Dirnagl, 2012). Interruption of blood flow eliminates the supply of oxygen and glucose and prevents the brain from generating ATP. Ion and water imbalances have long been considered as chief contributors to acute, irreversible brain injury together with a complex sequence of pathological events, including spreading depolarizations (SDs), excitotoxicity, calcium overload, generation of reactive oxygen species, and activation of cell death pathways (Lee et al., 1999; Hossmann, 2006; Hartings et al., 2016). All of these conditions involve mitochondria participation at some level.
Previously, the role of mitochondrial dysfunction has been addressed in vivo using endpoint measures of histological and biochemical markers (Matsumoto et al., 1999; Schinzel et al., 2005) or classic ultrastructural studies (Garcia et al., 1978; Solenski et al., 2002). Employment of in vitro preparations, such as primary cultured neurons, acute and organotypic brain slices have provided a mechanistic concept of mitochondria alteration during traumatic and ischemic injury-like events. Spreading depolarizations (Hartings et al., 2016), glutamate excitotoxicity (Rintoul et al., 2003; Abramov and Duchen, 2008; Brustovetsky et al., 2009), and/or ionic-oxidative stress lead to mitochondrial Ca2+ overload. This triggers the mitochondrial membrane potential collapse, generation of excessive reactive oxygen species, cessation of oxidative phosphorylation, and the opening of the mitochondrial permeability transition pore, which results in the eventual outer membrane rupture with release of cytochrome c and other apoptogenic proteins (Halestrap, 2006).
Many lines of evidence obtained in vitro suggest that mitochondrial functions are directly related to their morphology. Abnormal mitochondrial dynamics and fragmentation of the mitochondrial network into small spherical structures are considered hallmarks of mitochondrial injury (Rugarli and Langer, 2012; Kremneva et al., 2013; Otera et al., 2013). Such collapse of the mitochondrial network may be the result of either membrane permeability transition or disruption of mitochondria dynamics induced by traumatic or ischemic conditions. However, the relevance of these in vitro observations for in vivo conditions remains largely unknown.
Recent advances in optical techniques, such as two-photon laser scanning microscopy (Svoboda et al., 1997; Peron et al., 2015) and genetic animal models (Feng et al., 2000; Misgeld et al., 2007), allow repetitive real-time visualization of subcellular structures. Previously, we have developed an automatic approach, based on supervised statistical computer learning, for quantifying the degree of mitochondrial fragmentation in cortical neurons, and validated it using the irreversible global ischemia paradigm (Lihavainen et al., 2015). Here, we took advantage of this automatic analysis algorithm and performed a longitudinal in vivo two-photon imaging study of the neuronal mitochondrial fragmentation and recovery in response to mild, moderate, and severe injury using several animal models of acute neurological disorders.
Materials and Methods
Animals and surgical procedures
We used male and female adult mice ranging from 2 to 5 months. In total, 73 mice were used: 37 Thy1-mitoCFP mice (B6.Cg-Tg(Thy1-CFP/COX8A)C1Lich/J; RRID:IMSR_JAX:007940) (Misgeld et al., 2007), 10 Thy1-YFP(H) mice (B6.Cg-Tg(Thy1-YFP)HJrs/J; RRID:IMSR_JAX:003782), 6 Thy1-GFP(M) mice (B6.Cg-Tg(Thy1-EGFP)MJrc/J; RRID:IMSR_JAX:007919) (Feng et al., 2000), 5 dual Thy1-GFP(M)×Thy1-mitoCFP mice, and 15 C57BL/6J mice. The mice were bred and kept in group cages in the certified animal facilities of the University of Helsinki and Medical College of Georgia and provided with food and water ad libitum. Local authority (ELÄINKOELAUTAKUNTA-ELLA) approved the animal license (ESAVI/2857/04.10.03/2012) that included all procedures conducted at the University of Helsinki. All procedures conducted at Medical College of Georgia followed National Institutes of Health guidelines for the humane care and use of laboratory animals and underwent yearly review by the Animal Care and Use Committee at the Medical College of Georgia. Mice were allowed free access to room air and were anesthetized with intraperitoneal injection of mixture of ketamine (80 mg/kg body weight) and xylazine (10 mg/kg) during surgery and all longitudinal imaging experiments. In six terminal imaging experiments followed by transcardiac perfusion for tissue collection for electron microscopy (EM), mice were anesthetized with an intraperitoneal injection of urethane (1.5 mg/g). In these mice, the trachea was cannulated, and animals were ventilated mechanically with an SAR-830 ventilator (CWE). The depth of anesthesia, blood oxygen saturation level (>90%), and heart rate (450–650 beats/min) were monitored continuously with a MouseOx pulse oximeter (STARR Life Sciences) equipped with a mouse thigh sensor. A heating pad was used to maintain the core temperature of 37°C. Implantation of cranial window was performed as described previously (Holtmaat et al., 2009; Risher et al., 2010). Briefly, a subcutaneous injection of 0.1% lidocaine was administered to reduce local pain at the incision site, and then skin and connective tissue attached to the skull were removed. To form an ∼3 × 3 mm round cranial window, cranial bone over the somatosensory cortex was carefully removed. The brain was then covered with sterile cortex buffer containing the following (in mm): 125 NaCl, 5 KCl, 10 glucose, 10 HEPES, 2 CaCl2, and 2 MgSO4, pH 7.3, and a 5-mm-diameter #1.5 glass coverslip (Electron Microscopy Sciences) was placed over the window and sealed with a custom-made metal holder using dental cement. Mice recovered for 21–40 d before two-photon imaging. In 6 mice used for terminal imaging experiments followed by EM studies, an optical chamber was constructed by covering the intact dura with a thin layer of 1.5% agarose prepared in a cortex buffer. The chamber was left open to facilitate access with a glass recording microelectrode. The Ag/AgCl pellet ground electrode (A-M Systems) was installed under the skin above the nasal bone.
In vivo two-photon microscopy imaging of blood flow and neuronal mitochondria
Mice were imaged with a FV1000MPE two-photon microscope (Olympus) using the 25×/1.05 NA (numerical aperture) water-immersion objective. To achieve stable imaging and precise repositioning between imaging sessions, mice were placed under the microscope by attaching the metal holder to the custom-built frame fixed to the motorized stage (Prior Scientific). Fluorescence was excited with Mai Tai DeepSee titanium:sapphire femtosecond laser (SpectraPhysics) tuned to 860 nm. As measured by using the FWHM of a point spread function obtained with subresolution beads, the microscope has the resolution of 0.332 μm in the lateral dimension and 1.640 μm in the axial dimension. Emission light was collected by PMT detectors with a pinhole entirely opened using the following bandpass filters: 420–500, 515–560, and 590–650 nm. The laser power at the objective front focal plane (0.05–0.40 μJ) required for image acquisition was adjusted to the lowest possible level for each animal. Three-dimensional (3D) images were taken at 0.5–1 μm increments across a large scale 500 × 500 μm imaging field, yielding a nominal spatial resolution of 2.05 pixels/μm (12 bits/pixel, 2 μs pixel time) or across 169 × 169 μm field (6.06 pixels/μm, 12 bits/pixel, 2 μs pixel time) for higher resolution. A 0.1 ml bolus of 5% (w/v) 70 kDa dextran conjugated with Texas Red was injected into the tail vein, and blood flow was imaged in a repetitive line scan mode (at a rate of 1.5 kHz for 30 s) along the central axis of surface arterioles and venules. Moving red blood cells (RBCs) produce dark bands with a slope on the resulting space-time images (Kleinfeld et al., 1998). The slope was used to determine RBC velocity and flow direction.
Cortical injuries
Four different injury types, ranging from severe to mild, were inflicted in anesthetized mice, and the same volume of affected cortex was imaged before and after injury at several time points. At least two control volumes were selected at the distances between 0.5 and 2.0 mm away from the injury site. They were imaged before and after injury at several time points. To compare changes in mitochondrial morphology associated with each injury type, image sessions were performed during the acute (0–3 h), subchronic (2, 4, and 7 d), and chronic (2–3 weeks) phases of injury.
Global ischemia model was used to cause severe brain injury, followed by short-term imaging for 1–2 h. Transient global cerebral ischemia was induced by bilateral common carotid artery occlusion (BCCAO) during two-photon imaging (Murphy et al., 2008; Risher et al., 2012). Briefly, common carotid artery bifurcations were exposed after skin incision on the midline of the ventral neck and separated from the vagus nerves and connective tissue. One suture (3-0, polyester, Ethicon) was placed around each artery. For transient occlusion, sutures were tensioned on the microscope stage for 5 min (n = 5 mice) or 8 min (n = 7 mice) and then subsequently released for reperfusion.
Photothrombotic stroke model was used to cause severe focal brain injury that can be followed by longitudinal imaging for several weeks. Rose Bengal photosensitization (RBPS) was induced in the somatosensory cortex using the photothrombotic method (Watson et al., 1985; Labat-gest and Tomasi, 2013) of unilateral focal ischemic stroke with the following modification. Three or four large-scale z-stacks were acquired before RBPS. Then, a cortical region, 720 μm in diameter, near cranial window border was illuminated with a green epifluorescent light (535 ± 25 nm) through a 25×/1.05 NA objective (the average power through the objective was 3–10 mW) after the tail vein injection of the Rose Bengal solution (10 mg/kg, in PBS). Typically, photosensitization was performed for 1–5 min until clot formation was confirmed visually in the majority of the vessels. RBC velocity and blood flow direction were measured 20 min before RBPS and then every 30 min for 3 h. Ischemic core and remote areas were identified by these measurements.
A moderate brain injury was achieved with the focal laser lesion (FLL). FLL was applied at the depth of 80–110 μm from dura mater by targeting a spherical region of 50 μm in diameter with the laser beam and irradiating it for 0.42 ± 0.16 s with the 1.98 W of infrared laser energy at 800 nm, resulting in 5–20 μJ at the objective focal plane. To minimize bleeding, the FLL region was selected far away from large vessels. The FLL injury typically was confined to a roughly spherical area with a 60–530 pl volume surrounding the laser targeted region (as indicated by the immediate loss of CFP or YFP fluorescence and disruption of capillary network), which was surrounded by a perilesion site of 40–60 μm.
A mild photo-damage (MPD) was used to inflict a mild brain injury. MDP was achieved by zooming in with the factor of 3× to the center of a large scale image, resulting in 165 × 165 μm field and subsequent z-stack imaging of 100 optical sections with 1 μm spacing and increased laser energy at the 860 nm laser wavelength, yielding 0.5–1 μJ in the objective focal plane the imaging field. The site and degree of injury were strictly controlled; thus, the response of mitochondrial morphology to the MPD was highly reproducible.
Image processing and quantification
Because the relatively poor z-axial resolution of two-photon microscopy hindered proper 3D object analyses, mitochondria were analyzed in individual optical sections. With this approach, we accurately measured the “projected length” of mitochondria in the 2D plane of the optical section. Segmentation of individual mitochondria within each optical section was completed with custom MATLAB-based scripts (The MathWorks) (Lihavainen et al., 2012). Briefly, images from z-stacks were passed to the Difference-of-Gaussians (DOG) filter that removed low-frequency fluctuations in intensity. Mitochondria binary masks were generated from the resulting image by thresholding the pixel values. To remove false-positive objects and obtain better mitochondria size matching, a morphological erosion algorithm was applied. In addition, some connected components were discarded based on statistical properties associated with each locality computed from the original image. A set of morphological features was extracted from high resolution mitochondria images using National Institutes of Health ImageJ software (RRID: SCR_003070) to represent their morphology and subtypes (tubule and globule) classified by MATLAB-based MicroP program (Peng et al., 2011).
The resolution and signal-to-noise ratio of the large scale images limit accurate segmentation of individual mitochondria. To overcome these limitations, we used our recently developed supervised-learning-based analysis algorithm (Lihavainen et al., 2015) that estimates the degree of fragmentation of neuronal mitochondrial populations in a large region. The fragmentation score represents a continuous variable ranging from completely intact (score = 0) to completely fragmented mitochondria (score = 1). The method automatically extracts the texture and morphological features from an image patch (50 μm block to 101 pixels; see Fig. 2A), based on training of the classifier and regression model. For training, a set of images is visually examined and each patch is classified and scored. Then, the previously trained classifier detects the regions containing mitochondria, and the fragmentation score is obtained from the regression model. Here, we used the classifier and regression model trained earlier on the dataset from control Thy1-mitoCFP mice and mice exposed to cardiac arrest (Lihavainen et al., 2015).
Image stacks acquired from Thy1-YFP(H) and Thy1-MitoCFP×Thy1-EGFP(M) mice were preprocessed with ImageJ software. To analyze dendritic morphology in dual transgenic mice, CFP signal was subtracted from the GFP channel. All z-stacks were stabilized and aligned over time. To quantify structural changes in dendritic morphology after injury, the percentage of beaded dendrites was determined from maximal intensity projections, which were cropped into 25 μm blocks. Each block was manually assigned as beaded or intact based on criteria reported previously (Murphy et al., 2008; Sword et al., 2013; Steffensen et al., 2015). Data were double-checked with automated MATLAB-based BlebQuant program (Chen et al., 2011).
Automatic segmentation with Otsu's thresholding was applied to convert each optical section of labeled dendritic tuft to a binary image and to identify pixels corresponding to either dendritic structure or to the background. Density of dendritic tufts was calculated as a sum of foreground pixels per micrometer of lesion or perilesional area. The same calculation was done to compute density of mitochondria from binary image segmented as described above for “projected length” measurement.
To assess spine density, the same dendritic segments (20–45 μm) of the tuft that projected laterally were identified and cropped from z-stack at different time points. Obtained substacks were deconvolved in AutoQuant X3 (Media Cybernetics; RRID:SCR_002465). Dendritic spines were counted manually using 3D tools of NeuronStudio (Rodriguez et al., 2008), and spine density was calculated as the number of spines per micrometer of the dendritic segment length.
Experimental procedures used to obtain samples for EM analysis
Transient global cerebral ischemia was induced on the microscope stage by BCCAO as described above. The cortical slow DC potential and spontaneous EEG activity were recorded with a glass microelectrode (filled with 0.9% NaCl, 1–2 mΩ) inserted through the dura to the site of imaged dendrites within layer I of the somatosensory cortex. Signals were recorded with a MultiClamp 200B amplifier, filtered at 1 kHz, and digitized at 10 kHz with a Digidata 1322A interface board (Molecular Devices). Ischemia was verified by laser speckle imaging, and the onset of dendritic beading was monitored with two-photon imaging.
In vivo two-photon microscopy imaging of dendritic beading
Images were collected with infrared-optimized 40×/0.8 NA water-immersion objective (Carl Zeiss) using the Zeiss LSM 510 NLO META multiphoton system mounted on the motorized upright Axioscope 2FS microscope. The scan module was coupled directly with the Spectra-Physics titanium:sapphire broadband, mode-locked laser (Mai-Tai) tuned to 910 nm for two-photon excitation. The 3D time-lapse images were taken at 1 μm increments (with ∼20 sections per z-stack) across a 75 × 75 μm imaging field. After confirmation of the intact dendritic structure in sham-operated mice (n = 3) or dendritic beading after BCCAO (n = 3), mice were perfusion-fixed through the heart at 180 mmHg with mixed aldehydes (2% PFA, 2.5% glutaraldehyde, 2 mm CaCl2, 4 mm MgSO4, in 0.1 m sodium cacodylate, pH 7.4), and brain was processed for EM.
Electron microscopy and image analysis
A small piece of tissue containing somatosensory cortex was microdissected from the coronal section of the brain and subsequently processed with standard microwave-enhanced procedures through osmium, uranyl acetate, dehydration with a graded ethanol series, and embedding in Epon–Araldite resin (Kirov et al., 1999). Thin sections at 150 μm from the cortical surface were cut with a diamond knife on an EM UC6 ultramicrotome (Leica), collected on pioloform-coated copper Synaptek slot grids (Electron Microscopy Sciences), and stained with uranyl acetate and lead citrate. These protocols produced well-stained and readily identifiable neuronal processes (see Fig. 3D). Five to 13 fields were photographed randomly from each mouse (total 46 fields) at the JEOL 1230 transmission electron microscope (JEOL). Images were captured digitally at a magnification of 5000× using the UltraScan 4000 camera (Gatan). The images were randomized, coded, and analyzed blind as to condition. Each cross-sectioned dendrite containing cross-sectioned mitochondria was identified, and mitochondria cross-sectional area was measured with the freely available RECONSTRUCT software (Fiala, 2005) (RRID:SCR_002716). Each image provided ∼107 μm2 to evaluate mitochondria swelling. For serial section transmission electron microscopy (ssTEM), samples were obtained at 150 μm from the cortical surface from 2 sham-operated and 2 mice subjected to BCCAO. One series from each animal, containing 167–220 sections, each ∼53 nm thick, were cut and collected on pioloform-coated slot grids, and serial sections were photographed at 5000×. The 3D alignment, surfaced reconstructions, and analyses were completed using RECONSTRUCT.
Laser speckle imaging.
Two-dimensional (2D) maps of cerebral blood flow (CBF) with high spatiotemporal resolution were acquired by laser speckle imaging as described previously (Dunn et al., 2001; Sigler et al., 2008; Risher et al., 2010). Briefly, the cortical surface was illuminated through an anamorphic beam expander (Edmund Optics) by a 785 nm StockerYale laser (ProPhotonix) at an angle of ∼30° and imaged with a 4×/0.075 NA objective (Zeiss). Real-time speckle imaging was used with custom written LabView software (National Instruments; RRID:SCR_014325) (Yang et al., 2011) and modified for use with the Pantera 1M60 camera (Dalsa) using the XCLIB DLL library (EPIX). Laser speckle contrast was obtained by dividing the SD image by the mean of each raw image with a 5 × 5 pixel sliding window for immediate display. Fifty frames were saved individually as 32 bit images in LabView and averaged offline using ImageJ to obtain a single image for figures.
Statistics
Origin software (OriginLab) was used for statistical analyses and to plot graphs. A two-tailed paired t test and one-way ANOVA followed by Tukey's post hoc test were used to compare group means at different time points. One-sample and unpaired Student's t test was applied for analyses of dendritic beading percentage. Kolmogorov–Smirnov test and Mann–Whitney rank sum test were used for statistical analysis of cross-sectional mitochondria area in EM. Changes in RBC velocity were expressed as percentage of baseline values. χ2 test was used to analyze mitochondria subgroups in two-photon and in EM images. The linear regression analysis was applied to quantify the strength of the relationship between two variables. The sample size of each experimental group is given in Results. All data analyses were performed blind to the experimental conditions. Data are presented as mean ± SEM, except mitochondrial fragmentation score, which is given as median ± SEM. The significance criterion was set at p < 0.05.
Results
Structure of cortical neuronal mitochondria in vivo
We used in vivo two-photon microscopy and chronic cranial window preparation (Holtmaat et al., 2009; Paveliev et al., 2014) combined with Thy1-mitoCFP transgenic mouse strain (Misgeld et al., 2007) to study structure and distribution of mitochondria labeled with CFP fused to the mitochondria-targeted sequence from subunit VIII of human cytochrome c oxidase. It was previously shown that this mitochondrial targeting sequence is sufficient for ubiquitous and exclusive expression of fusion proteins in the mitochondria of transgenic mice (Shitara et al., 2001; Chandrasekaran et al., 2006; Misgeld et al., 2007). Volume reconstruction of somatosensory cortex showed dense CFP labeling from cortical layer V and II/III pyramidal neurons (Fig. 1B). Such a high number of labeled cells did not allow distinguishing and tracing mitochondria in individual neurons, compared with sparse labeling of dendrites in Thy1-YFP(H) mice (Fig. 1C).
As visualized by mitoCFP in neocortex of adult mice, in vivo neuronal mitochondria appeared continuous, branched, and varicose with no abnormalities in morphology, such as fragmentation, vacuolization, or aggregation (Fig. 1D). Here, we selected mitochondria in molecular layer I and granular layer II/III and analyzed their morphology using the automated segmentation method (Lihavainen et al., 2012). Based on “projected length” of mitochondria segmentation in the 2D plane of the optical section, neuronal mitochondria were classified into globules and tubules using the automatic object subtyping method (Peng et al., 2011). Mitochondria highly varied in size (Fig. 1E) with most dendritic mitochondria having projected lengths of <5 μm (average 2.8 ± 0.6 μm), whereas a small proportion of mitochondria were >6 μm long. Some mitochondria within cortical layer I were extremely long and formed large networks that extended well beyond 20 μm (Popov et al., 2005). Average distribution of mitochondrial projected lengths calculated from 9 animals revealed that the projected length decreased significantly as a function of depth from the pia (Fig. 1E). The average projected length of mitochondria was 38.2 ± 2.4% smaller at the depth of 180 μm below cortical surface compared with that at the 50 μm depth (p < 0.001, n = 9 mice, n = 800 mitochondria per mouse per depth). Similarly, the mitochondrial tubular-to-globular ratio changed significantly with depth (p < 0.01) (Fig. 1G). Tubular mitochondria were predominant in superficial cortical layers, whereas globular mitochondria were observed more frequently in the granular layer II/III.
This observation may have at least two underlying reasons. First, the projected length of mitochondria may decrease with depth because the full length of mitochondria becomes smaller on average in deeper cortical layers. Second, projected length may decrease simply because the angle between the mitochondria-bearing dendrites and the vertical axis becomes sharper with depth (Fig. 1C). Although we were not able to directly test the former hypothesis, we addressed the latter by analyzing the projected length of dendritic branches in the well-characterized Thy1-YFP(H) mice (Fig. 1F,H). Indeed, we observed very similar regularities in the dendrite projected lengths and in the depth dependence of object subtype ratios for the dendritic labeling as we did for mitochondrial labeling. Therefore, it is likely that the mitochondrial depth dependence in projected length and tubular-to-globular ratio are determined primarily by the angular position of dendritic tufts to the horizontal imaging plane. Thus, to minimize the bias, further mitochondrial analysis was conducted only in dendritic tufts within molecular layer I, where dendrites run almost parallel to the brain surface.
Rapid, reperfusion-dependent mitochondrial fragmentation in the BCCAO model of global ischemia
In the BCCAO model, abrupt reduction of blood flow results in a severe energy deficit and, in the absence of immediate reperfusion, leads to a permanent dendritic beading that marks irreversible neuronal injury (Murphy et al., 2008; Risher et al., 2012). In this respect, it resembles the ischemic core when compared with focal stroke. In our previous study, when global ischemia was induced by cardiac arrest (Lihavainen et al., 2015), the peak of neuronal mitochondrial fragmentation was observed within 5 min after heart failure. Therefore, it was desirable to select duration of occlusion that would result in severe but reversible fragmentation. In the preliminary set of experiments, 9 wild-type C57BL mice were injected intravenously with the 70 kDa Texas Red dextran and RBC velocity was monitored during BCCAO. Occlusion caused a sharp drop in RBC velocity to 6.4 ± 4.7% of baseline toward the end of the fifth minute. More prolonged interruption of blood flow resulted in severe vasoconstriction as evidenced by a decrease in arterioles diameter from 24.1 ± 3.7 μm at baseline to 6.8 ± 0.5 μm at 8 min of occlusion (p < 0.01, n = 16 arterioles in 9 mice). Based on these observed changes in blood flow and mitochondrial fragmentation, we selected 5 min as a minimal sufficient duration of common carotid artery occlusion. Therefore, to study reperfusion-dependent mitochondrial fragmentation, the first group of animals (n = 5) was subjected to a brief 5 min BCCAO and the second group (n = 5) was subjected to a more prolonged 8 min occlusion.
In both groups, mitochondrial fragmentation began to develop during the first 5 min of occlusion (Fig. 2B–D,F), and it was linearly correlated with the reduction of blood flow (r = −0.719, p < 0.001; Fig. 2E,F). Indeed, by the end of the 5 min occlusion, the mitochondrial fragmentation score increased from the baseline level of 0.27 ± 0.06 to 0.67 ± 0.06 (Fig. 2C,D,F). Reperfusion during the first 3 min after a brief 5 min BCCAO resulted in a rapid restoration of blood flow to 81.9 ± 7.7% of baseline with a complete recovery by 30 min (108.4 ± 7.29%) (Fig. 2E; black trace). Recovery of blood flow was correlated with recovery of mitochondria tubular structure (r = −0.346, p < 0.05; Fig. 2B–D,F). Indeed, the fragmentation score significantly decreased from 0.65 ± 0.05 at the beginning of reperfusion to 0.41 ± 0.06 at 30 min (p = 0.02; paired t test, n = 5 mice) (Fig. 2F, black trace). In contrast to the brief 5 min BCCAO group, reperfusion in the second group of mice subjected to 8 min occlusion was incomplete, with blood flow recovering to only 81.8 ± 8.3% of baseline at 30 min of reperfusion and with 91.8 ± 8.5% of recovery at 90 min (Fig. 2E, red trace). The 8 min BCCAO also resulted in a further increase in the mitochondrial fragmentation score to 0.84 ± 0.06 at the end of occlusion (Fig. 2F, red trace). Moreover, no significant recovery of mitochondrial tubular structure was revealed upon reperfusion with fragmentation score, decreasing to only 0.60 ± 0.05 at 90 min of reperfusion (p < 0.05 relative to baseline, paired t test, n = 3 mice), which suggests irreversible fragmentation (Fig. 2F, red trace). Together, these findings indicate that neuronal mitochondria in vivo undergo fragmentation in parallel with a loss of blood flow during onset of global cerebral ischemia, and the timing and efficiency of reperfusion are crucial in determining the degree of reversibility in mitochondrial structural changes.
Next, mitochondrial fragmentation in response to global cerebral ischemia seen with two-photon microscopy was verified with quantitative EM analyses. Three sham-operated control urethane-anesthetized Thy1-GFP(M) mice were imaged with two-photon microscopy and then, after confirmation of the intact dendritic structure, were perfusion fixed through the heart. In the other 3 mice, rapid global ischemia was induced on the microscope stage with DC potential and EEG monitoring (Fig. 3A) (Murphy et al., 2008; Risher et al., 2012). Occlusion was accompanied by reduction of blood flow (Fig. 3B) followed by a rapid SD onset at 6.1 ± 2.0 min after beginning of BCCAO (Fig. 3A). As expected, BCCAO-induced SD invariably beaded dendrites (Fig. 3C), and then mice were immediately perfusion fixed for EM within 12.1 ± 1.7 min of occlusion and without any reperfusion. Electron micrographs from neuropil of sham-operated control mice showed intact dendrites with structurally unimpaired mitochondria (Fig. 3D). On the contrary, dendrites from animals subjected to stroke revealed signs of swelling with loss of microtubules and severely swollen mitochondria. Indeed, single-section EM analyses shown that the average cross-section dendritic mitochondria area after BCCAO (0.223 ± 0.011 μm2, n = 409 mitochondria) was significantly larger than in sham-operated control mice (0.062 ± 0.002 μm2, n = 393 mitochondria) (p < 0.001, Mann–Whitney rank sum test; Fig. 3E,F). The largest mitochondria cross-sectional area in control mice was 0.2 μm2. In comparison, in mice subjected to BCCAO, 36.4% of mitochondria had cross-sectional area >0.2 μm2 (Fig. 3F), revealing mitochondrial swelling.
Mitochondrial fragmentation that was first revealed with two-photon microscopy (Fig. 2) was supported with ssTEM. Long filamentous mitochondria have been previously found in dendrites of hippocampal neurons (Popov et al., 2005). Accordingly, every reconstructed dendritic segment (8.67–14.19 μm in length, average 11.86 ± 0.70 μm, n = 7) of cortical neurons from two series obtained from two sham-operated mice had long mitochondrial filaments (Fig. 3G). Transverse serial sections of seven dendritic segments allow reconstruction of 12 mitochondria with a length of 1.79–13.81 μm (average 8.52 ± 1.34 μm; Fig. 3H). It is actually unclear how long these mitochondria were because only two of these mitochondria (2.54 and 4.38 μm in length) from these dendritic segments were fully confined within series containing 167 and 203 sections, with the rest extending beyond the first or the last section, or both. In contrast, 11 reconstructed dendritic segments (6.4–15.67 μm in length, average 11.27 ± 0.95 μm; Fig. 3G) from two series from two mice subjected to BCCAO had 47 short discrete mitochondria with a length of 0.15–4.74 μm (average 1.54 ± 0.16 μm; Fig. 3H). Some of these globular mitochondria were still interconnected by thin mitochondrial segments indicating ongoing fission events known to lead to the fragmentation of filamentous mitochondrial network, which facilitates mitophagy or can result in apoptosis after severe damage (Rugarli and Langer, 2012). Forty of these globular mitochondria (0.29–4.74 μm in length, average 1.53 ± 0.17 μm) from 11 dendritic segments were fully reconstructed within 170 and 220 serial sections, with only 7 mitochondria being cut off at the beginning or at the end of the series. Thus, ssTEM confirmed that global ischemia was associated with fragmentation and substantial ultrastructural changes in dendritic mitochondria, not a redistribution of CFP.
Longitudinal imaging of mitochondrial fragmentation in the RBPS model of focal stroke
Bilateral occlusion of common carotid arterys allows studying the acute phase of ischemia during severe metabolic compromise, but longitudinal studies of the recovery phase are difficult due to high mortality rates and great variability of this model. Therefore, we used the RBPS model of focal stroke, which allows better controlled longitudinal studies of mitochondrial fragmentation in the ischemic core and areas of the tissue located 1.5–2 mm away from the core (Fig. 4A). RBPS resulted in the permanent ischemia in the core area, where RBC velocity was reduced by 87.3 ± 12.4% and 94.1 ± 5.6% as measured at 30 min and 3 h, respectively (n = 16 blood vessels in 7 mice), whereas in the remote area blood flow remained unaffected (from 8.7 ± 4.6 mm/s at baseline to 9.3 ± 5.5 and to 8.4 ± 5.0 mm/s, at 30 min and 3 h, respectively, n = 19 blood vessels in 7 mice). In agreement with the data from the global ischemia model with a prolonged 8 min occlusion, severe mitochondrial fragmentation was observed in an injury core with a score that increased from 0.27 ± 0.07 at baseline to 0.77 ± 0.11 at 20 min after RBPS (p < 0.01, n = 7 mice; Figs. 4B,C, 5A). Mitochondria fragmentation continued to develop as revealed by an increased fragmentation score at 2 d later (0.93 ± 0.05, p < 0.001, n = 7 mice). Indeed, we observed no recovery of mitochondrial morphology over 3 weeks. In contrast, in the remote area, mitochondrial fragmentation score gradually increased, reaching significance at 3 h after RBPS compared with the baseline score (from 0.31 ± 0.07 to 0.53 ± 0.07, respectively; p < 0.05, n = 7 mice, paired t test; Figs. 4D,E, 5B). Mitochondrial fragmentation continued to rise, attaining a score of 0.79 ± 0.06 on day 2. However, recovery of mitochondrial structure (0.39 ± 0.07, p = 0.52, n = 4 animals) was observed 3 weeks after RBPS. Together, these data indicate that mitochondria are extremely sensitive to the loss of blood flow and can serve as an in vivo indicator of ischemic tissue.
Dendritic beading is an early morphological sign of acute neuronal injury associated with a variety of pathological conditions, including ischemia (Hori and Carpenter, 1994; Zhang et al., 2005; Risher et al., 2010). Therefore, we quantified mitochondrial fragmentation and dendritic beading concurrently during the first 3 h following RBPS in the core and in the remote area using dual transgenic Thy1-MitoCFP×Thy1-EGFP(M) mice. Similarly to Thy1-mitoCFP mice, mitochondrial fragmentation score in dual transgenic mice increased in the core from 0.34 ± 0.09 at baseline to 0.72 ± 0.11 at 20 min after RBPS (p < 0.01, paired t test, n = 5 mice). In contrast to changes in mitochondrial morphology, dendritic injury in the core evolved relatively slow. Under our experimental conditions, only 19.7 ± 10.9% of beaded dendrites were found at 20 min after RBPS (p = 0.15, n = 5 mice), but most of the dendrites (78.1 ± 7.3%) were beaded by 2 h (p < 0.001) (Fig. 5C,E). With blood flow largely preserved in the remote areas, it was plausible that dendritic structure in these areas will remain mostly unchanged. Accordingly, dendritic structure remained stable throughout 3 h of imaging following RBPS, with the percentage of beaded dendrites equal to 7.3 ± 5.5% and 20.5 ± 8.1% at 20 min and 3 h, respectively (Fig. 5D,F). Together, these data provide strong evidence that mitochondrial structure is very sensitive to the loss of blood flow in contrast to dendritic structure, which was more stable in our RBPS model of focal ischemia.
Reversible fragmentation of dendritic mitochondria in the FLL model
We observed mitochondria loss in the ischemic core of the focal stroke and irreversible changes in mitochondrial structure after prolonged global ischemia. The loss of mitochondria in focal cerebral ischemia might simply reflect the loss of neurons at the site of ischemic lesion. Subsequently, irreversible mitochondrial fragmentation observed at 90 min after BCCAO might also be attributed to permanent damage to neuronal somata after prolonged global ischemia. Thus, we decided to limit the damage to dendritic compartment to test whether mitochondrial fragmentation is reversible after moderate injury to dendritic tufts by FLL. FLL was produced by laser ablation in the absence of any photosensitizing agents, such as Rose Bengal. FLL resulted in the loss of CFP fluorescence and elimination of individual branches of the apical dendrites in a controlled volume (typically 60–530 pl) at the site of injury (Fig. 6A,B). Spontaneous dendritic regrowth into the lesion site occurred during 7–14 d, and it was correlated (r = −0.73, p < 0.02) with restoration of CFP fluorescence (Fig. 6C–G). Interestingly, although FLL induced a significant mitochondrial fragmentation at the perilesion site (Fig. 6C,D), there were no structural changes in perilesional dendrites (Fig. 6E). Indeed, mitochondrial fragmentation score increased from 0.21 ± 0.07 before FLL to 0.52 ± 0.08 at 30 min after FLL at the perilesion site (p < 0.05; n = 5 mice; Fig. 6F), whereas there were no changes in the percentage of the beaded dendrites (15.6 ± 9.7%, 30 min after FLL, p = 0.18, one-sample t test, n = 5 mice). Full recovery of mitochondrial structure in the perilesion perimeter of the damaged volume was observed during 1 week (Fig. 6F). Together, these data indicate that localized injury to dendritic compartment that spares neuronal somata leads to only a transient loss of dendrites and their mitochondria. These data further emphasize that mitochondria morphology is easily affected by brain trauma.
Reversible fragmentation of mitochondria in intact dendrites after the MPD
Next, we decided to establish a mild brain injury model in which we could examine dynamics of mitochondrial fragmentation in intact dendritic compartment. We chose to induce mild injury by means of a photo-damage because the degree of injury could be tightly controlled by laser intensity and exposure time. First, we found that neuronal mitochondria and dendrites exhibited different sensitivity to laser light intensity. At the infrared laser energy levels ranging from 5 to 20 μJ (as measured in the objective focal plane), both mitoCFP and dendritic YFP fluorescence was lost. However, at the lower laser energy level from 0.5 to 1 μJ, dendrites remained intact whereas mitochondria became clearly fragmented. Finally, at the laser energy ranging from 0.05 to 0.40 μJ, both mitochondria and dendrites remained intact. Therefore, we selected the energy range between 0.5 and 1 μJ to determine mitochondrial fragmentation and recovery in intact dendrites during injury inflicted by MPD.
MPD was induced in a small cortical volume (typically 165 × 165 × 100 μm). Within minutes after MPD induction, mitochondria were found to become fragmented in the photo-damaged region (Fig. 7A,B, white box). Indeed, mitochondrial fragmentation score increased from 0.23 ± 0.08 to 0.76 ± 0.09 in the first 30 min after injury (p < 0.01; n = 6 mice; Fig. 7C). The fragmentation score declined over the next several days, reaching preinjury level at 7 d after laser exposure (p = 0.29, paired t test, n = 6 mice). Notably, the fragmentation did not extend to the surrounding tissue (0.21 ± 0.05 before MPD vs 0.35 ± 0.06 at 30 min after MPD; p = 0.10; n = 6 mice; Fig. 7D). To define spatial distribution of fragmented mitochondria over time, we generated digital surfaces over heatmap pixels with a score >0.6. Volume containing fragmented mitochondria significantly shrank from 0.49 ± 0.04 μm3 at 30 min after MPD to 0.17 ± 0.05 μm3 on day 4 after injury (Fig. 7E). These findings suggest that re-fusion (and/or replacement) of fragmented mitochondria in vivo is a relatively slow process that may take several days.
As mentioned above, MPD did not trigger any pathological changes in dendritic structure, such as beading either acutely or over 2 weeks following MPD (Fig. 7F). However, MPD resulted in a significant decrease in dendritic spine density from 0.44 ± 0.07 μm before to 0.32 ± 0.05 μm at 2 d after MPD (p < 0.01, n = 5 mice; Fig. 7G). Spine density remained declined at day 4 (0.37 ± 0.02 μm; p < 0.01, n = 5) but recovered by the end of 1 week. This recovery of dendritic spine density was paralleled by recovery of tubular mitochondrial structure, suggesting a novel link between mitochondria integrity and spine turnover in the adult brain.
A brief summary of our findings is presented in Figure 8.
Discussion
These findings are the first to show that in vivo neuronal mitochondrial fragmentation occurs within minutes of injury onset and it is either localized in the case of mild or moderate damage, or spreads gradually in the case of severe injury. Mitochondrial structural changes are accompanied by dendritic damage during severe global ischemia and in the core of focal injury. Yet, during mild trauma or in the areas located far away from the injury core, dendrites responded differently from mitochondria, remaining intact while mitochondrial fragmentation occurred. Over time, mitochondria were able to regain their unfragmented morphology in dendrites that were spared from terminal injury. The present work provides the first direct in vivo demonstration of such heterogeneous sensitivity to the injury between dendrites and their mitochondrial organelle.
Mitochondria are very dynamic organelles that can move within dendrites to the site of active synapses (Li et al., 2004; Chang et al., 2006) and undergo fusion and fission to maintain optimal morphology to match specific bioenergetic demands (Chen et al., 2005; Parone et al., 2008). Fission and fusion of mitochondria are controlled by specific proteins. It was shown that transient cerebral ischemia causes upregulation of two major players, the fission dynamin-related protein 1 (Drp1) and the fusion protein optic atrophy 1 (Opa1), in the ischemic penumbra peaking at 2 d after stroke (Liu et al., 2012). Fusion and fission are necessary to retain healthy mitochondrial population through fusion, which is thought to maintain mitochondrial bioenergetics, and fission, which partitions injured mitochondria from healthy mitochondrial network (Twig et al., 2008; Mouli et al., 2009). Upregulation of Drp1 and Opa1 after stroke may indicate increased fission and fusion events in the attempt for mitophagy and cell survival. Here, using morphometric analysis of mitochondria, we were able to accurately evaluate dynamic changes in the mitochondrial tubular-to-globular ratio after cortical injuries. When damage was limited to dendritic compartments and neuronal somata were not affected, as in the core of focal laser-lesion and MPD models, or in the remote cortical areas away from the core of focal stroke model, we have shown unequivocal recovery of tubular mitochondrial structure. The slow time course of this recovery over several days to weeks suggests that it mainly occurred by replacing damaged mitochondria with new organelles formed in those neuronal compartments that were spared from the injury. Similar time course of mitochondrial structural recovery was observed in the core of focal laser lesion; notably, it co-occurred with spontaneous regrowth of dendrites. Indeed, accumulating evidence indicates that the endogenous neuroprotective response to cerebral ischemia induces mitochondrial biogenesis that creates new functional mitochondria (Yin et al., 2008; Anne Stetler et al., 2013). Interestingly, recovery of mitochondria after MPD was paralleled by recovery of dendritic spine density, suggesting that these two events are closely related as mitochondria are thought to dock in an activity-dependent manner near dendritic spines to support synaptic transmission (Li et al., 2004; Chang et al., 2006). Hence, it appears that postinjury stimulation of mitochondrial recycling and recovery may play a major role in neuronal survival and restoration of synaptic connectivity and could become a potential therapeutic target.
Conceivably, mitochondrial fragmentation observed in our injury models involves rapid fission-like events. Previously, fission of mitochondria was detected with live confocal microscopy imaging in cultured cortical neurons following oxygen-glucose deprivation (Wappler et al., 2013) and with EM following middle cerebral artery occlusion in mice (Barsoum et al., 2006). A recent histological study by Owens et al. (2015) confirmed these EM observations in brains of mice perfusion-fixed at several time points after common carotid artery occlusion-induced global ischemia. Interestingly, widespread fission of neuronal mitochondria was detected in all hippocampal subregions, but mitochondria were able to re-fuse in ischemia-resistant CA3 and dentate gyrus neurons, suggesting that increased fusion might be neuroprotective in these cells. Here, using ssTEM, we revealed mitochondrial fragmentation after BCCAO and applied real-time in vivo imaging to show that fragmentation developed rapidly within 5 min after occlusion. Moreover, once triggered by noxious stimuli, fragmentation continued to develop in the BCCAO model, as well as in all other injury models we studied. Yet, in contrast to prolonged 8 min BCCAO, mitochondria were able to re-fuse by 30 min after a brief 5 min occlusion. To the best of our knowledge, these are the first in vivo experiments to show such fast reversible perturbation in mitochondrial structure in response to ischemia-induced injury. It is conceivable that the balance between excessive mitochondrial fragmentation and re-fusion after ischemia may underlie transition from injury to repair in tissue at risk of permanent brain damage, but delayed neuronal death was not evaluated in our study, which warrants future experiments to address this issue.
Previous studies have shown that excessive mitochondrial fragmentation causes increased production of reactive oxygen species, cytochrome c release, and apoptosis (Youle and Karbowski, 2005; Knott et al., 2008; Anne Stetler et al., 2013). Lack of oxygen during ischemia or trauma inhibits mitochondrial oxidative phosphorylation, leading to mitochondrial depolarization and morphological disruption. Indeed, inhibition of mitochondrial oxidative metabolism with chemical treatment inevitably results in mitochondrial fragmentation (Benard et al., 2007; Liot et al., 2009). Mitochondrial depolarization (e.g., in the presence of electron transport chain uncoupling drugs) elicits mitochondrial fission (Ishihara et al., 2003; Meeusen et al., 2004). A fission event by itself can result in mitochondrial depolarization, and then depolarized mitochondria can be selectively disposed by autophagy (Twig et al., 2008). In this respect, widespread mitochondrial fragmentation detected in our models of acute neurological disorders could signify the cellular response to increased number of injured mitochondria that must be eliminated.
In their pioneering in vivo work, Liu and Murphy (2009) used two-photon imaging of the mitochondrial membrane potential (Δψm) with rhodamine 123 during global ischemia. In agreement with brain slice studies (Bahar et al., 2000), they found that Δψm depolarization coincided with SD and occurred within ∼90 s of BCCAO. The Δψm depolarization was cyclosporine A sensitive, indicating activation of the mitochondrial permeability transition pore, which plays an important role in delayed neuronal death. SD also triggers dendritic beading and accelerates neuronal injury (Kirov, 2014), but beading and mitochondrial depolarization during SD represent two mechanistically separable pathways as blocking Δψm depolarization in vivo with cyclosporine A did not affect dendritic beading (Liu and Murphy, 2009). Similarly, in cultured neurons, glutamate excitotoxicity-induced dendritic beading and mitochondrial dysfunction (depolarization and structural collapse into rounded/swollen structures) are independent processes, despite their simultaneous occurrence (Greenwood et al., 2007). Indeed, SD-induced dendritic beading results from water accumulation (Kirov et al., 2004; Andrew et al., 2007) and involves, at least in part, activation of select chloride-coupled neuronal cotransporters (Steffensen et al., 2015; Sword et al., 2016). Studies have shown that the collapse of mitochondrial structure is a consequence of calcium overload (Pivovarova et al., 2004) and water accumulation, which are both present during SD (Greenwood et al., 2007). Accordingly, using EM, we demonstrated that SD-beaded dendrites were swollen and contained watery cytoplasm with swollen mitochondria. This is in agreement with previous ultrastructural studies of cortical neurons that revealed swollen dendrites and mitochondria after only 5 or 10 min of global ischemia without reperfusion (Tomimoto and Yanagihara, 1992, 1994). A rupture of the outer mitochondrial membrane and release of apoptotic proteins may be one of the possible harmful consequences of such excessive mitochondrial swelling (Brustovetsky et al., 2002).
Gradual development of dendritic beading was observed over 2 h in the core of RBPS model. Although the cortical DC potential was not recorded, the slow degeneration of dendritic structure is consistent with the lack of SD because SD would rapidly bead dendrites and dramatically accelerate dendritic injury (Murphy et al., 2008; Risher et al., 2010). Without recurrent SDs, beading is gated by the degree of ischemia as a flowing capillary ∼80 μm away could maintain dendritic structure in the ischemic tissue (Zhang and Murphy, 2007), whereas arterioles can supply oxygen beyond >100 μm (Kasischke et al., 2011). It is conceivable that local residual blood flow in the irradiated core of our RBPS model would initially support dendrites, whereas gradual loss of dendritic structure would parallel a gradual loss of flowing capillaries. It should be noted that, without SD, severely metabolically compromised dendrites ultimately will bead due to the cation influx driven by Gibbs-Donnan forces and insufficient cation efflux caused by inhibition of ATP-dependent sodium pump (Dreier et al., 2013). Such mechanism is different from mechanism of rapid SD-induced dendritic beading (Steffensen et al., 2015; Sword et al., 2016). In contrast, rapid mitochondrial fragmentation in the core of RBPS model reflects high sensitivity of mitochondrial structure to the loss of blood flow and lack of oxygen compared with the dendrites.
In conclusion, we present in vivo evidence of rapid loss and relatively slow recovery of dendritic mitochondrial structural integrity in several models of acute neurological disorders. Because mitochondrial structural rearrangements can influence neuronal function and affect neural pathological outcomes, we expect that novel therapeutic approaches targeting the mitochondrial organelle will benefit from these in vivo experimental insights into injury-induced mitochondrial dynamics.
Footnotes
This work was supported by the Academy of Finland Grant SA-126321, FGSN, CIMO, and University of Helsinki to L.K., Emil Aaltonen Foundation to E.L., Academy of Finland Grants 295027 and 305342, and Jane and Aatos Erkko Foundation Grant 610536 to A.S.R., and by National Institutes of Health Grant NS083858 to S.A.K. We thank Libby Perry and Brendan Marshall (Electron Microscopy Core at the Medical College of Georgia) for assistance with electron microscopy; and Julia Kolikova (Neurotar Ltd) and Marina Tibeykina (In Vivo Microscopy Unit at the University of Helsinki) for assistance with animal surgeries.
E.P. and D.T. are part-time employees of Neurotar Ltd. L.K. is a co-owner of Neurotar Ltd. The remaining authors declare no competing financial interests.
- Correspondence should be addressed to either of the following: Dr. Leonard Khiroug, In Vivo Microscopy Unit, Neuroscience Center, PO Box 56 (Viikinkaari 4), University of Helsinki, FI-00014, Helsinki, Finland, leonard.khirug{at}helsinki.fi; or Sergei A. Kirov, Multiphoton Imaging Laboratory, Department of Neurosurgery and Brain and Behavior Discovery Institute, Medical College of Georgia at Augusta University, 1120 15th Street, CB-3706, Augusta, GA 30912-2630, skirov13{at}gmail.com