|Home | About | Journals | Submit | Contact Us | Français|
Genetically encoded calcium indicators (GECIs) are powerful tools for systems neuroscience. Recent efforts in protein engineering have significantly increased the performance of GECIs. The state-of-the art single-wavelength GECI, GCaMP3, has been deployed in a number of model organisms and can reliably detect three or more action potentials (APs) in short bursts in several systems in vivo. Through protein structure determination, targeted mutagenesis, high-throughput screening, and a battery of in vitro assays, we have increased the dynamic range of GCaMP3 by several-fold, creating a family of “GCaMP5” sensors. We tested GCaMP5s in several systems: cultured neurons and astrocytes, mouse retina, and in vivo in Caenorhabditis chemosensory neurons, Drosophila larval neuromuscular junction and adult antennal lobe, zebrafish retina and tectum, and mouse visual cortex. Signal-to-noise ratio was improved by at least 2–3-fold. In the visual cortex, two GCaMP5 variants detected twice as many visual stimulus-responsive cells as GCaMP3. By combining in vivo imaging with electrophysiology we show that GCaMP5 fluorescence provides a more reliable measure of neuronal activity than its predecessor GCaMP3. GCaMP5 allows more sensitive detection of neural activity in vivo and may find widespread applications for cellular imaging in general.
Calcium is a ubiquitous second messenger, playing an essential role in excitable cells and signal transduction. Calcium ions (Ca2+) enter neurons during action potential (AP) firing and synaptic input. AP firing and synaptic inputs can therefore be assessed, sometimes quantitatively, by measuring changes in intracellular [Ca2+](Yasuda et al., 2004). Genetically encoded calcium indicators (GECIs)(Mank and Griesbeck, 2008) and small molecule calcium-sensitive dyes (Cobbold and Rink, 1987) are both used to report [Ca2+] changes; but GECIs have the advantage that they enable chronic, non-invasive imaging of defined cells and compartments (Mao et al., 2008). State-of-the-art GECIs include the Förster resonance energy transfer (FRET) indicators D3cpVenus (D3cpV)(Palmer et al., 2006), TN-XXL (Mank et al., 2008), and YC3.60(Nagai et al., 2004), and the single-wavelength sensor GCaMP3(Tian et al., 2009). GCaMP3 is based on circularly permuted green fluorescent protein (cpGFP), calmodulin (CaM), and the Ca2+/CaM-binding “M13” peptide (M13pep). Several versions of the original GCaMP sensor (Nakai et al., 2001) have been published (Ohkura et al., 2005; Tallini et al., 2006; Akerboom et al., 2009). Recent versions include “GCaMP4.1” (Shindo et al., 2010), which was used to image Xenopus gastrulation, but no sequence information or comparison with other GECIs is published. “GCaMP-HS” consists of GCaMP2 with a subset of the “superfolder GFP” mutations (Pedelacq et al., 2006) and was used for imaging zebrafish motor neurons (Muto et al., 2011), but was also not compared to other sensors. The “G-GECO” sensors (Zhao et al., 2011) were created from GCaMP3 by random mutagenesis; they show ~2x greater fluorescence increase in purified protein (Ca2+-saturated vs. Ca2+-free) but were not tested in neurons. However, the sensors are dimmer than GCaMP3 in both the Ca2+-free and Ca2+-bound states, which can complicate imaging.
GCaMP3 has been used to detect activity in large neuronal populations in the motor cortex (Tian et al., 2009), barrel cortex (O’Connor et al., 2010), and hippocampus (Tank et al., 2010) of behaving mice. Long-term imaging of GCaMP3 has revealed learning-related circuit changes in vivo (Huber et al., 2012). GCaMP3 imaging has also been used to probe dendritic excitation in layer 5 dendrites in vivo (Xu et al., 2010; Mittmann et al., 2011), light-evoked responses in populations of neurons in mouse retina (Borghuis et al., 2011), zebrafish tectum (Del Bene et al., 2010), and walking Drosophila (Chiappe et al., 2010; Seelig et al., 2010), among others. However, GCaMP3 remains an imperfect GECI. Single action potentials are not reliably detected in vivo, and detection of active cells lags synthetic indicators (Tian et al., 2009). GCaMP5 was engineered from GCaMP3 using a combination of structure-guided design (Akerboom et al., 2009) and semi-rational library screening. We have enhanced sensitivity by increasing the dynamic range of the fluorescence response (ΔF/F = (F−F0)/F0)), the Ca2+-saturated brightness, and the Ca2+ affinity. We characterized ten new GCaMP5s in vitro and in vivo under a wide variety of conditions. Each sensor has improved properties relative to the parent sensor GCaMP3. The ideal GCaMP5 indicator for a given study can be selected from this set according to particular experimental requirements.
Site-directed mutagenesis of GCaMP3 was carried out using the method of Kunkel (Kunkel et al., 1991) or the QuikChange methodology (Agilent Technologies, USA). Single-stranded uracil-containing DNA template of pRSET-GCaMP3 for Kunkel mutagenesis was produced according to established protocols. Mutants were confirmed by DNA sequencing.
Libraries were generated using primers containing degenerate codons (NNS) following the site-directed mutagenesis protocols described above and transformed into E. coli XL1-Blue (Stratagene/Agilent, USA). The following day, colonies were scraped off plates, combined, and the plasmid library was isolated with Qiagen Miniprep kits (Qiagen, Germany) following the procedures provided, eluting in 100 μl water. 1 μl of the library was subsequently transformed into E. coli BL21(DE3) (EMD4 Biosciences, USA), plated on 244 × 244 mm square LB-agar plates containing 100 μg/ml ampicillin, and grown for 20 hours at 30 °C. Colonies were selected using a colony picker (QPix2XT, Genetix, UK) and grown in 800 μl ZYM-5052 medium (Studier, 2005) containing 100μg/ml ampicillin in 96 deep-well blocks for 48 hours at 30 °C, shaking vigorously at 700 rpm. 2 μl was taken from each well and mixed into a fresh deep-well block containing 800 μl LB medium + 100μg/ml ampicillin, grown overnight at 37 °C, pelleted and stored at 4 °C for sequence analysis. The E. coli BL21(DE3) cultures in the deep-well blocks containing the over-expressed mutants were pelleted by centrifugation (4000 × g, 20 min, 4 °C), frozen, thawed, resuspended in lysis buffer (20 mM TRIS, pH 8.0, 100 mM NaCl, 1 mg/ml lysozyme, 1.5 Kunitz units/ml DNAseI (Qiagen, Germany)) and subsequently incubated at 30 °C for 2–4 hours, shaking. Lysates were clarified by centrifugation (4000 × g, 30 min, 4 °C), and 100 μl was taken from each well into Greiner Bio-One black 96-well fluorescence plates (Greiner, Germany) in duplicate. To one plate 1 μl 100 mM CaCl2 was added (final Ca2+ concentration ~ 1 mM), to the other 1 μl 100 mM EGTA, pH 7.4. Fluorescence was measured in a Tecan Sapphire2 Spectrophotometer (Tecan Group Ltd., Switzerland), at 485 nm excitation and 510 nm emission, 5 nm slits, gain=90 V.
Cloning, expression in E. coli and purification was performed essentially as before (Akerboom et al., 2009). Briefly, for expression in BL21(DE3), pRSETa-GCaMP variants were transformed to BL21(DE3), and single colonies were grown for 3 days in ZYM-5052 media (Studier, 2005) at 25 °C, after which cells were pelleted by centrifugation. Cells were lysed by resuspending them in 4 × w/v lysis buffer 2 (20 mM TRIS.HCl pH 8.0, 100 mM NaCl), followed by a freeze-thaw cycle and subsequent cell-rupture by passing the cell suspension three times through a pre-cooled (4 °C) Avestin Emulsiflex-C5 (Avestin, Canada) and finally a 15 second sonication step at 30 mW amplitude on ice (Fisher Dismembrator Model 100 equipped with 3 mm tip). Lysate was clarified by centrifugation (30,000 × g, 4 °C, 45 min.), and cell-free extract (CFE) was incubated on a rotary incubator with 5% (v/v) Profinity IMAC Ni-NTA resin (Biorad, USA) at 4 °C for 4–16 hours. Resin was allowed to settle in 20 ml disposable columns (Biorad, USA), and was washed with 20 ml lysis buffer 2, followed by 10 ml wash buffer (20 mM TRIS.HCl pH 8.0, 100 mM NaCl, 10 mM imidazole). Proteins were eluted into elution buffer (20 mM TRIS.HCl, pH 8.0, 100 mM NaCl, 300 mM imidazole) and subsequently dialyzed extensively into lysis buffer 2 using Spectra/por membrane tubing (Spectrumlabs, USA), with a molecular weight cut-off (MWCO) of 10,000 Da, at 4 °C. Concentration and purity was determined using Agilent Protein 230 chips on a Agilent 2100 Bioanalyzer (Agilent Technologies, USA) as well as NuPAGE NOVEX BIS-TRIS SDS-PAGE gels (Invitrogen, USA). Proteins were concentrated when necessary using Millipore Centrifugal Filter Units of 10,000 Da MWCO (Millipore, USA).
5 μl purified GCaMP solution (~100 μM) was added to 100 μl pre-mixed, pH titrated triple-buffer (10 mM glycine, 10 mM sodium citrate, 10 mM TRIS.HCl) containing either 3 mM BAPTA or 1 mM CaCl2, in duplicate. Fluorescence was measured in 96-well Greiner Bio-One transparent fluorescence plates (Greiner, Germany) in a Tecan Sapphire2 Spectrophotometer (Tecan Group, Ltd. Switzerland), at 485 nm excitation and 510 nm emission, 5 nm slits, gain=90 V).
5 μl purified GCaMP solution (~100 μM) was added to 100 μl of either zero-free calcium buffer (10 mM EGTA, 100 mM KCl, 30 mM MOPS, pH 7.2) or 39 μM free calcium buffer (10 mM CaEGTA, 100 mM KCl, 30 mM MOPS, pH 7.2) (Invitrogen, USA). Excitation spectra were taken from 240 nm to 550 nm, emission was 600 nm. Emission spectra were taken from 470 to 700 nm, excitation wavelength was 430 nm (gain=80 V, 5 nm slit widths). Spectra were normalized to the calcium-saturated excitation and emission maxima, respectively. For absorbance measurements, 200 μl protein was dialyzed into 20 mM TRIS.HCl, pH 8.0, 100 mM NaCl, using Slide-a-Lyzer mini dialysis units with a 2,000 Da MWCO (Pierce, USA). The absorbance of 100 μl protein with 2 μl of either 100 mM CaCl2 or 100 mM EGTA added was measured from 240 to 700 nm, blanked against the chemically identical dialysis buffer containing 2 μl of 100 mM CaCl2 or 2 μl 100 mM EGTA. Spectra were normalized using the absorbance at 280 nm.
Calcium affinity assays were performed by mixing different volumes of the zero-free calcium buffer (10 mM EGTA, 100 mM KCl, 30 mM MOPS, pH 7.2) and 39 μM free-calcium buffer (10 mM CaEGTA in 100 mM KCl, 30 mM MOPS, pH 7.2) from the calcium calibration buffer kit (Invitrogen, USA) according to the manufacturer’s instructions. Measurements were performed by mixing 3 μl purified GCaMP (~100 μM) with 100 μl of different ratios of calcium buffer and 39 μM free-calcium buffer in 96 well transparent Greiner BIO ONE plates and measuring the fluorescence at 485 nm excitation and 510 nm emission, 5 nm slits, gain=90 V in duplicate as described earlier.
Photophysical properties for GCaMP2, GCaMP3, and the GCaMP5 variants A, D, and G were investigated in buffer solutions in the presence or absence of free Ca2+. Absorption and emission properties, including quantum yield, were determined using a UV/VIS spectrometer (Lambda 35, Perkin-Elmer) and a fluorescence spectrometer (LS-55, Perkin-Elmer). Two-photon-excited properties, including fluorescence excitation spectra, two-photon cross section, fluorescence decay time and emission spectroscopy, were conducted with laser pulses from an 80-MHz Ti: Sapphire laser (Chameleon Ultra II, Coherent). In all cases, near-IR laser pulses were focused into a solution containing the GCaMP proteins via a 60x, 1.2 NA water-immersion objective of an inverted epi-fluorescence microscope (IX81, Olympus), slightly overfilling the objective back aperture. Fluorescence generated by 2-photon excitation was collected by the same 60x objective used to excite the proteins, and after passing through a short-pass filter (FF01-720SP, Semrock) and a band pass filter (FF01-550/88, Semrock), was directed to the input face of a fiber-coupled avalanche photodiode (APD). Two different APDs were used, one optimized for speed for lifetime measurements (model PDF CCTB, Micro Photon Devices), and one optimized for low-noise for FCS and spectral measurements (SPCM_AQRH-14-FC, Perkin Elmer) coupled to a 100-μm-core multimode fiber (AFS105/125Y, Thorlabs). Output pulses generated by the fast-timing APD were fed to a TCSPC board (TimeHarp200, Picoquant). Output pulses from the low-noise APD were fed to an external autocorrelator (Flex03LQ-01, correlator.com) for spectra and FCS measurements. Emission spectra were recorded from an additional microscope port using a fiber-coupled 0.3-m spectrograph/CCD (model SP2358 with Pixis 256 CCD camera, Princeton Instruments).
Control of the amount of laser power delivered to the sample in the focal plane of the microscope was accomplished by rotating the linear polarization of the laser output using an achromatic -waveplate (AHWP05M-980, ThorLabs), mounted in a computer-controlled rotation stage, followed by passing this beam through an angle-fixed Glan-Laser polarizer (GL10-B, Thorlabs). We calibrated the system prior to measurements, where for each laser wavelength (in 10 nm increments), the laser power at the focus of the microscope objective was measured while the -wave plate was angle-adjusted, calibrating at each wavelength a specific power at the focus with a specific orientation of the waveplate. Laser wavelength and power adjustments, as well as data acquisition, were run under computer control. Data analysis was performed using Origin 8.0 software (OriginLab Corporation, USA), and FCS fits were performed using a custom software package (V. Iyer, unpublished).
Samples were prepared by dilution from stock solutions of purified proteins into buffers at either pH 7.25 or pH 9.5, in order to characterize the fluorescence properties under physiological pH and at elevated pH where most of the protein chromophores are in the deprotonated form in the presence of Ca2+. The pH 7.25 buffers (30 mM MOPS, 100 mM KCl) were either purchased as part of a Ca2+-calibration kit (C-3008MP, Invitrogen, USA), that contain either 10 mM K2EGTA or 10 mM CaEGTA, or prepared from MOPS and KCl stock solutions (Mediatech Inc., USA) and contain either 1 mM CaCl2 or 0.5 mM EGTA or 1 mM BAPTA. The pH 9.5 buffer (30 mM CHES, 100 mM KCl) was prepared from stock chemicals and was supplemented with either 1 mM CaCl2 or 0.5 mM EGTA. For measurements using very low amounts of protein (below 200 nM) such as FCS measurements and one-photon fluorescence measurements, we added 0.1 mg/ml BSA as a blocking agent to the buffer solutions to prevent loss of proteins from solution to the nearby glass or coverslip surfaces via absorption over the course of the measurements; this eliminated any decrease in fluorescence signal over the course of the measurements.
The concentration of chromophore-forming proteins was determined by two methods: alkali-denaturation (Ward, 2005), and two-photon-excited fluorescence correlation spectroscopy (FCS). For both methods, absorption spectra of either EGFP or GCaMP protein solutions were taken in the UV/VIS spectrometer at either pH 7.25 or pH 9.5. In the second method using FCS, GCaMP stock solutions were diluted several thousand-fold into pH 9.5 buffer (30 mM CHES, 100 mM KCl) containing 1 mM CaCl2, and excited at 960 nm in a fluorescence microscope over a range of laser powers. At each laser power, the mean fluorescence rate <F> and its fluctuations were recorded for 50–200 seconds, and the autocorrelation G(τ) of the fluorescence signal was computed. A fit to the autocorrelation function G(τ) based on a diffusion model determines the diffusion coefficient of the proteins, and the mean number of fluorophores in the excitation volume, given by <N> = 1/G(τ−>0) (Schwille et al., 1999). This measurement is repeated for a sample containing EGFP at known concentration (determined by alkali denaturation) and dilution in pH 9.5 buffer, providing a reference between a known concentration and a number of fluorophores in the excitation volume.
Quantum yield (QY) was determined for the fluorescent proteins in both pH 7.25 and pH 9.5 buffer in the presence of free Ca2+ using standard methods that measure the optical absorption and total fluorescence yield of samples at a fixed wavelength for both sample and a standard fluorescein with a QY of 0.93 in aqueous 0.1 N NaOH (Magde et. al.,- 2002) with approximately the same emission spectrum and emission peak.
The fluorescence lifetime was measured by time-correlated single-photon counting (TCSPC) using two-photon excitation at 960 nm in a fluorescence microscope setup, where detector pulses from the fast-timing APD and trigger signals from a PIN diode monitoring the laser pulse train, were fed to the TCSPC board. To achieve improved performance, a pulse-picker (Model 350-160, Conoptics, Inc.) was inserted in the beam to reduce the laser pulse frequency from 80 MHz to 20 MHz. The fluorescence lifetime of GCaMP samples was determined in either pH 9.5 buffer supplemented with CaCl2 or EGTA as described earlier, or pH 7.25 buffer (Invitrogen buffer with either 10 mM CaEGTA or 10 mM EGTA). The lifetime reference for the system was fluorescein (e−1 lifetime τ = 4.1+/−0.1 ns); in our setup the fluorescein decay was well fit to a single-exponential decay (4.0 nsec, χ2 =1.08). Measured fluorescent decays were fit to a single-exponential decay curve, or to a two-exponential decay, which improved the goodness-of-fit.
Two-photon excitation spectra and ΔF/F were measured for the GCaMPs in the presence and absence of free calcium at pH 7.25 and pH 9.5, respectively as described above, at 1 μM protein concentration. Two-photon spectra are taken with constant laser power delivered to the sample, although due to a wavelength-dependent pulse width of the femtosecond pulses, and changes in focal spot size (focused beam diameter scales as the excitation wavelength), the laser intensity varies gradually across the spectrum. We do not correct for this variation in intensity. Together with each run of GCaMP samples, a reference two-photon excitation spectrum of fluorescein was recorded, allowing us to determine the absolute two-photon cross-section of the GCaMPs using published cross-sections measured for these fluorophores (Xu and Webb, 1996; Drobizhev et al., 2011).
Fluorophores can be characterized by their specific brightness under two-photon excitation, measured in counts per second per molecule, at a specific laser intensity and wavelength. This is the average fluorescence rate detected per fluorescent molecule, and is measured using two-photon-excited FCS. This quantity reaches a maximum or peak value as the laser intensity is increased, beyond which the fluorescence rate per molecule decreases with higher intensity, due to photobleaching of the fluorophores in the volume of the focused laser beam. While the peak brightness will strongly depend on the molecular environment (in vitro versus intracellular/in vivo), this value can be used as a quantitative measure in comparing the photostability of different fluorophores.
To determine the peak brightness, GCaMPs were diluted to nominally 50 nM in pH 9.5 buffer containing either 1 mM CaCl2 or 0.5 mM EGTA. Solutions also contained 0.1 mg/ml BSA to prevent absorption of fluorophores to the nearby coverslip surface. As a control, EGFP at 50 nM was prepared and measured in the same buffer, without CaCl2 or EGTA. Measurements were taken for a series of laser powers (with power measured in the focal plane) at 940 nm, where fluorescence time course data was acquired for 50 – 200 seconds at each laser intensity. For FCS, we used the low-noise APD. The TTL output of the APD was fed to an autocorrelator and associated software to generate two quantities, the time-average fluorescence rate <F> and the measured autocorrelation G(τ) of the fluorescence data. FCS theory equates the quantity 1/G(0) to <N>, the average number of emitting molecules in the excitation volume. By acquiring <F> and G(t) simultaneously for each protein over a range of laser intensities, and fitting G(0) to determine <N>, we can define the two-photon brightness at each intensity as <F>/<N>, the effective detected fluorescence rate per emitting molecule at each intensity. This quantity has a maximum value or peak brightness, since as the intensity is raised, saturation and photobleaching begin to diminish the fluorescence rate.
Without using FCS, GCaMPs can be characterized by the power dependence of their (ΔF/F)max. For these measurements, GCaMPs were diluted to 0.5 μM in pH 7.25 buffer +/− Ca2+ and fluorescence recorded for laser excitation at 940 nm under increasing power over the range of 0.5 mW – 70 mW.
For expression in HEK293 cells, DNA was PCR amplified from the pRSETa-GCaMP construct, purified, and digested with BglII and NotI, and ligated into digested pEGFP-N1 (Takara-Bio, Japan), which had been digested with BglII and NotI to remove EGFP. Successful clones, with the EGFP gene replaced with the gene coding for GCaMP, were confirmed by sequence analysis, and plasmids were prepared using the endo-free plasmid Maxi kit (Qiagen, Germany). Equal amounts of plasmid for each GCaMP5 variant were transfected into HEK293 cells using the 96-well Nucleofector protocol in Amaxa plates (Lonzo, Germany), with each variant in 16 wells for 8 duplicate measurements. Cells were grown and incubated at 37 °C for 2 days, after which growth medium was aspirated and replaced with 100 μl pre-warmed (37 °C) cell buffer (1× TBS, 2 mM CaCl2). Whole-well fluorescence response traces evoked by addition of 100 μl acetylcholine in cell buffer (dilution series 10−2 –10−10 M) were recorded by the liquid handling 96-well Hamamatsu FDSS (Hamamatsu, Japan).
For cloning in pAAV, DNA was PCR amplified from pCMV-GCaMP clones, purified and digested with BamHI and HindIII (NEB, USA), gel-purified and ligated into pre-digested pAAV-humanSynapsin-1 vector. Ligations were transformed into E. coli Stbl2 (Invitrogen, USA) and plates were incubated for 2 days at 30 °C. Insertions were verified by DNA sequencing. Virus was produced by standard protocols.
GCaMP proteins were expressed in E. coli BL21(DE3) using the pRSETa plasmid and purified by immobilized metal ion affinity and size exclusion chromatography as described (Rodriguez Guilbe et al., 2008; Akerboom et al., 2009). Proteins were concentrated to 5–10 mg/ml for crystallization as described earlier. Crystallization was carried out at room temperature by sitting-drop vapor diffusion using commercially available sparse-matrix screens (Hampton Research, USA) by mixing 1.2 μl of protein solution with 1.2 μl of precipitant in 96-well plates. Individual proteins were crystallized under conditions listed in Table 2. Crystals were cryo-protected prior to X-ray data collection. X-ray data were reduced using Mosflm (Leslie, 1992)/Scala (Collaborative Computational Project, 1994) or D*Trek (Pflugrath, 1999) (Table 2). Each structure was solved by molecular replacement using Phaser (McCoy et al., 2007), or simple rigid body refinement starting from previously published GCaMP structures using the CCP4 package (Collaborative Computational Project, 1994). Iterative cycles of model building in Coot (Emsley and Cowtan, 2004) and refinement in Refmac/CCP4(Collaborative Computational Project, 1994) led to the final models described in Table 2.
To check the extent of in vitro dimerization of GCaMP5A, GCaMP5B and GCaMP5G compared to GCaMP2 and GCaMP3, an aliquot of each purified mutant was adjusted to 150 μM in lysis buffer 2, after which EGTA was added to 1 mM final concentration. This mixture was allowed to equilibrate at room temperature for 30 minutes, after which CaCl2 was added to 2.5 mM final concentration. 100 μl aliquots of each GCaMP were then injected onto a Superdex 200 10/300 GL column (GE Healthcare, USA) with 20 mM Tris, pH 8.0, 100 mM NaCl, 2 mM CaCl2 as the running buffer, and protein was eluted at a flow rate of 0.5 ml/min.
Primary cultures of hippocampal neurons were obtained from P0 rat pups by dissection, papain-based dissociation, and plating onto Matrigel-coated (BD Biosciences, USA) 24-well glass-bottom plates (MatTek, USA) and then cultured in DMEM/B27 medium (Invitrogen, USA). SIV-based lentiviral vectors containing hsyn1-GCaMP variant-IRES-nls-mCherry-WPRE-SV40 polyadenylation signal (hsyn1: human synapsin-1 promoter, IRES: internal ribosome entry site, nls: nuclear localization signal, WPRE: woodchuck hepatitis virus post-transcriptional regulatory element) constructs were produced through quadruple transfection of HEK293T cells in 10-cm plates. Harvested lentiviral particles were used to infect hippocampal neuronal cultures on the third day in vitro for 16 hours, and medium was replaced with DMEM/B27/4 microM AraC (Invitrogen & Sigma). On day 16–18 in vitro, infected neurons were stimulated using a custom-built 24-well multiplexed field stimulator with platinum wires and imaged using an Olympus IX81 motorized, inverted microscope (10X objective, 0.4 NA, Chroma ET-GFP or ET-TxRed filter sets) Prior Scientific H117 ProScan™ II motorized stage, Cairn Research optical feedback OptoLED illumination system and an EMCCD camera (Andor iXon+ 897, 34.8 frames per second). Field stimuli were delivered at 40V, 83 Hz, 1 millisecond pulses for the following trains: 1, 2, 3, 5, 10, 20, 40, 80, 160 field stimuli. The whole system was automated using Metamorph (Molecular Devices) and MATLAB (MathWorks) software. Imaging buffer included (in mM): 145 NaCl, 2.5 KCl, 10 glucose, 10 HEPES (pH 7.4), 2 CaCl2, 1 MgCl2, 0.01 3-(2-Carboxypiperazin-4-yl)-propyl-1-phosphonic acid (CPP, Tocris), 0.01 6-Cyano-7-nitroquinoxaline-2,3-dione (CNQX, Tocris), 0.01 gabazine (Tocris), and 1 α-methyl-4-carboxyphenylglycine (MCPG, Tocris). Images were processed and analyzed using custom software.
Methods were identical to those in ref (Shigetomi et al., 2010b). Briefly, we used an Olympus IX71 microscope equipped with an IXON DV887DCS EMCCD camera (Andor), epifluorescence condenser, control unit and Polychrome V monochromator (TILL photonics). The control of excitation and image acquisition was achieved using TILLVision software. We used an Olympus 60X 1.45 NA objective lens. Images were typically taken every one second. Exposure time and pixel binning were optimized to visualize fluorescence signals for each experiment (maximum binning was 4 × 4). Cultures were perfused with recording buffer (110 mM NaCl, 5.4 mM KCl, 1.8 mM CaCl2, 0.8 mM MgCl2, 10 mM D-glucose, 10 mM HEPES at pH 7.4 (adjusted with NaOH).
For imaging astrocyte activity in neuron/astrocyte co-cultures, Lck-GCaMP3 and Lck-GCaMP5G were sub-cloned into an AAV vector driven by the astrocyte-specific GFAP promoter (Brenner et al., 1994). Neuron-astrocyte co-culture was infected at day 3 in vitro and cells were images 10 days post-infection.
Methods were identical to (Borghuis et al., 2011).
The same experimental protocol was used as in previous GCaMP experiments (Tian et al., 2009). All recordings are from the AWCon sensory neuron. After 10 seconds of observation, odor ([IAA] was 10−4 M) was added for 5 minutes. At the last 10 seconds of odor addition, imaging was commenced for another minute to observe the neuron’s response to odor removal (the AWCon neuron activates in response to odor removal).
Imaging data was analyzed differently than our previous work (Tian et al., 2009). Expression in C. elegans from extra-chromosomal arrays can show considerable animal-to-animal variation in expression level. This complicates analysis of imaging results; we have attempted to improve our image analysis to better cope with differences in expression level. The previous analysis used the total fluorescence intensity of pixels above threshold to ensure only pixels from the neuron were quantified. However, fluctuations in fluorescence intensity across animals also produced large fluctuations in the number of pixels analyzed, which resulted in highly variable measurements between samples. To dampen this effect, the average intensity of pixels above threshold was used.
To allow imaging access to transgenic presynaptic neuromuscular junction (NMJ) boutons which expressed variants of GCAMP3, third instar Drosophila larvae were dissected using methods similar to those described previously (1976). Genetic constructs were prepared in the pJFRC7-20XUAS-IVS vector (Pfeiffer et al., 2010) and inserted in the VK00005 site (Bellen et al., 2006) to allow expression in motor neurons using the OK6-Gal4 driver (Goodman et al., 2002). The combination of OK6-Gal4, VK00005 and pJFRC7-20XUAS-IVS provided good labeling of Type 1b and 1s, but not Type II boutons in heterozygous animals (+/yellow1 white1118; +/OK6-Gal4; +/UAS-GCaMPxx). Actively wandering larvae were dissected in ice-cold Schneider’s insect medium (Sigma) pH shifted to 7.2 using NaOH. A 35 mm Petri dish previously 1/3 filled with Sylgard (Dow Corning, USA), was used to pin down the cuticle and body wall muscles. Utmost care was taken in the dissection to reduce the potential for muscle movement during imaging (ensuring central longitudinal muscles were not contacted with instruments at any time, cuticle stretching was sufficient to keep the preparation in place but not so much to cause spontaneous calcium release, choosing pin insertion locations that minimized potential damage to imaging region). Imaging commenced within 15 minutes of the segmental nerves being cut close to the ventral nerve cord (VNC). Before imaging, the Schneider’s insect medium was replaced with HL6 (Macleod et al., 2002) supplemented with 2 mM Ca2+, an osmolarity of ~340±10 mOsm, pH of ~7.2±0.05, and 7 mM L-glutamic acid added to reduce muscle movement (Macleod et al., 2004). During imaging, pH and temperature were monitored, with pH values ranging from ~7.4–7.6 (start to end) while temperature stayed relatively constant (~22±0.5°C).
Cut segmental nerves were drawn by suction into a heat-polished glass pipette, ~12 μm internal diameter (Macleod et al., 2002), which was connected to an ISO-Flex stimulator (A.M.P.I. Jerusalem, Israel) to deliver supra-threshold electrical pulses (~2.9 V, ~5 mA). After 3 s of baseline, each stimulus was 2 s in duration with individual 300 μs impulses delivered at a frequency of 1, 5, 10, 20, 40, 80 and 160 Hz. The inter-stimulus interval was ~20 seconds and each stimulus frequency was repeated five times, with a pause of ~35 seconds, before moving to a higher stimulus frequency. The stimulus order was kept constant as higher frequency stimuli were potentially deleterious and outside the normal range of firing frequencies (Macleod et al., 2010).
Wide-field Ca2+ imaging was performed on an upright Olympus BX61-WI microscope using a 40X water immersion objective (0.8 NA), a 2X lens extender (EX2C; Computar, Commack, NY) and an Andor EMCCD camera (Model DU897_BV, 512×512 pixels, 30 FPS, 100 EMgain, 5.2x pregain, −70°C, 32.9 ms exposure, Andor Technology, Belfast, UK) with 7% illumination from a DC powered X-CITE exacte light source (Lumen Dynamics, Ontario, Canada). The fluorescence filter set used was BrightLine 3035B (472/30, 495, 520/35; Semrock Rochester, NY). The imaging system was controlled by custom written journals in Metamorph (ver. 7.7.5, Molecular Devices, Sunnyvale, CA) while the stimulation and image timing was controlled (as a slave) by MATLAB (2010a, The MathWorks Inc, Natick, MA) using a USB-6343, X Series DAQ (National Instruments, Austin, TX) and custom written software, Ephus (Suter et al., 2010).
Fields of view (FOV) were exclusively from Muscle 13 due to its high native firing rate (Macleod et al., 2010), using predominantly segment A4, however segments A3 and A5 were used when 4 boutons (not including the tip bouton) were not simultaneously in focus or present. The left and right side of the animal were taken as separate samples to increase throughput of data collection, with the second imaging session approximately 20 minutes after the first. Once each FOV was selected and the corresponding nerve drawn into the glass pipette and connected to the stimulator, the Metamorph software (MM) controlled data acquisition.
An image-based auto-focusing routine (Adjust Focus) was used to find the optimal focus position using the Brenner algorithm based on nearest neighbor pixel values (Brenner et al., 1976). MM was configured to start the camera trigger and simultaneously start the stimulus timer in Ephus via a USB X Series DAQ. Using a Uniblitz shutter (Model VS25S1ZM1R3-24, Vincent Associates, Rochester, NY) set into the light path and controlled by either MM or Ephus, 11 dark frames were collected (shutter closed) first, followed by 394 frames with an open shutter. The image-based auto-focusing routine was repeated between each five replicate trials for each stimulus frequency acquired, to account for any small movements in the Z-axis.
Along with the image data, we collected analog signals for the temperature and pH of the bath, the current and voltage of the stimulus and the frame signals from the camera. Continuously monitoring temperature and pH enabled constant conditions across experiments for different calcium indicators. Recording the current, voltage and frame signals at 10 KHz enabled precise calculation of the timing of the calcium response.
MATLAB was used for all analysis. We performed background subtraction before calculating fractional changes in fluorescence from baseline values (ΔF/F0). During the 15-minute FOV acquisitions, we infrequently observed either slow XY drifts in the muscle position or actual muscle contractions, normally in neighboring muscles. This was despite 7 mM L-glutamic acid being used, which normally prevented muscle movement in acquisitions of <3 minutes (Macleod et al., 2004). In the rare cases that we observed noticeable changes in Z-position (<1:1000 trials), which were detected by changes in baseline fluorescence (F0), we discarded the data from the entire FOV. In order not to discard FOVs due to XY movement, we developed a custom MATLAB analysis routine that accounted for image movement between images in a stack (one stimulus, using cross correlation (Fienup et al., 2008)) and also between trials, by independently tracking each segmented ROI centroid, and moving the ROIs based on the centroid movement between trials. We only included data from each FOV if the ΔF/F0 continually increased with increasing stimulation and peaked at 80 Hz, as previously observed from intracellular recordings (Macleod et al., 2010). When this was not observed, presumably due to muscle damage or overstretching during the dissection, the data from the entire FOV was discarded.
For each GCAMP construct, 10 FOV that met the quality control parameters described above were used. These FOV were collected from 7 animals (i.e., more than one FOV collected from one animal). To describe the performance of an indicator, boutons within a FOV were averaged, replicate trials of a stimulus frequency were averaged, and the 10 FOV were averaged.
Flies were reared on standard cornmeal agar medium. We used the Gal4/UAS system (Brand et al., 1994) to direct the expression of the calcium sensors to PNs. GH146-Gal4 flies were a gift from L. Luo (Stanford University, Stanford, CA). All animals were adult females, 3–5 days after eclosion. Adult flies were dissected using previously described methods (Jayaraman and Laurent, 2007). Flies were anaesthetized in a vial on ice until movement stopped (<15 seconds) and then gently inserted into a hole in a piece of aluminum foil. Small drops of wax (55 °C) were used to suspend the fly in the hole, with the edge of foil defining a horizontal plane around the head and thorax, from the first antennal segment anteriorly to the scutellum posteriorly. The dorsal side of the foil was bathed in saline, while the ventral side (including antennae and maxillary palps) remained dry and accessible to odors. A window was cut in the dorsal head cuticle between the eyes, extending from the ocelli to the first antennal segment. Fat and air sacs dorsal and anterior to the brain were removed, but the perineural sheath was left intact. The proboscis was affixed with a small drop of wax to a strand of human hair to limit brain movement. Spontaneous leg movements were typically observed in this preparation for the duration of the recording (2–3 h). The saline composition used in all olfactory experiments was (in mM): 103 NaCl, 3 KCl, 5 N-tris (hydroxymethyl) methyl-2-aminoethane-sulfonic acid, 10 trehalose, 10 glucose, 26 NaHCO3, 1 NaH2PO4, 2.0 CaCl2, and 4 MgCl2, adjusted to 275 mOsm, pH 7.4.
Odors (different concentrations of octanol) were delivered using a custom-made odor-delivery system designed by Dmitry Rinberg, and a Teflon nozzle (entry diameter 1/8″) directed towards the antennae. Odors were delivered at different concentrations diluted in paraffin oil (Paraffin oil alone, 0.001%, 0.01%, 0.1%, 1.0% and 10%) in a constant stream of air (1 l/min) with an additional 10% dilution in air. For each concentration, five replicate deliveries were performed and the data averaged. Odor delivery times were measured using a mini-PID (Aurora Scientific Inc., Ontario, Canada). Odors were presented for 1s. All comparisons of sensor performance were made using experiments with identical odor presentation times. The results reported are based on data obtained from 5 GCaMP3-expressing flies (6 ALs) and 5 GCaMP5-expressing flies (6 ALs).
We imaged on a two-photon microscope using PrairieView software and an Olympus 40x, 0.8 NA LUMPlanFL/IR objective. A mode-locked Ti:Sapphire Chameleon Ultra II laser (Coherent, Santa Clara, CA) tuned to 930 nm was used for excitation. Fluorescence was collected using photomultiplier tubes (Hamamatsu, Japan) after bandpass filtering. Images were acquired in framescan mode (20 Hz) for a single plane of one antennal lobe. Fluorescence time series were then obtained by averaging across the spatial extent of the glomerulus in the frame. In all cases fluorescence changes were calculated relative to baseline fluorescence levels as determined by averaging over 2 seconds just before odor presentation.
mitfa−/− (nacre) zebrafish larvae (Lister et al., 1999) expressing Gal4 under the control of the pan-neuronal elavl3/HuC promoter (elavl3:Gal4) were injected at 1–8 cell embryonic stage with a solution containing 0.04% Phenol Red, 50 ng/μl transposase mRNA, and 25 ng/μl plasmid DNA with the coding sequence of GCaMP5A or GCaMP5G downstream of 14 upstream activation sequences (UASs) and flanked by two Tol2 sequences for stable genomic integration. An elavl3:Gal4, UAS:GCaMP3 stable line (Del Bene et al., 2010) was used for experiments with GCaMP3.
6 day post-fertilization (dpf) larvae were embedded in 2% low melting point agarose and imaged at 7dpf at 28°C. A custom-built two-photon microscope equipped with a mode-locked Ti:Sapphire Chameleon UltraII laser (Coherent) tuned to 920 nm and controlled by ScanImage v3.6 software (Pologruto et al., 2003) was used for acquiring image time series at 17Hz. Visual stimuli (moving vertical grey bars on black background) were generated with VisonEgg software (Straw, 2008) and presented with an 800X600 pixel organic light-emitting diode (OLED, eMagin, USA) with a green light filter to the left eye of the larvae. Imaging was performed in the right (contralateral) tectal neuropil. Image time series were x-y motion-corrected with a program written in MATLAB (Dombeck et al., 2007) and analyzed with ImageJ.
5 and 6dpf AB/nacre larval zebrafish expressing GCaMP2, GCaMP3 or GCaMP5G under the elavl3 promoter were paralyzed by immersing them in 1 mg/ml solution of bungarotoxin dissolved in E3 fish embryo water and were subsequently embedded in 2% low melting point agarose in a 35 mm Petri dish. They were placed in a custom 2-photon microscope and imaged using a Mai Tai HP Ti-Sapphire laser tuned to 950 nm. The visual stimulus used for the experiment consisted of a light dot (0.5 mm × 0.5 mm) projected, using an amber (590 nm) LED mounted into a miniature LCOS projector, onto an opal glass screen directly underneath the larvae. Stimulus light was filtered with a narrow bandpass filter, and each fish was run through one stimulus set with the laser off to detect stimulus bleed-through, which was always negligible. The dot appeared to the left or right of the larva and moved on a straight line at a speed of 3 mm/s until it disappeared on the opposite side. The larva was located halfway along the dot’s trajectory and perpendicular to it, with the point of closest approach of the dot being 0.5 mm rostral to the larva.
The experimental protocol consisted of 1 min dark, followed by a presentation every 30 s of the moving dot, alternating between left to right and right to left. There were ten such presentations (five in each direction). The experiment concluded with 1 min dark, and therefore lasted 7 min in total. Individual frames were captured at 138.32 ms per frame (7.23 Hz), using a quad-interlaced scan pattern that ensured that each cell was sampled evenly at 4 times this frame rate.
Data analysis: Movies were assessed for x-y drift during the experiment (usually < 1 pixel), and a sub-pixel translation correction was applied using MATLAB software (David Heeger, NYU). Neuronal somata were detected based on their dark nuclei. Mean images were smoothed with a Gaussian, and local minima were detected. These were classified as cell nuclei if the ratio of the brightness 3 pixels from the center was more than 3.5 the brightness 1 pixel from the center, i.e. they look like a bright ring around a dark centre, and they were sufficiently bright (>17,500 photons detected per experiment). Fluorescence was then averaged over a 7×7 pixel square. Baseline fluorescence (F) was defined as the average fluorescence in the 50 frames immediately preceding each left-right stimulus.
All procedures were carried out according to the UK Animals (Scientific Procedures) Act 1986 and approved by the UK Home Office. Fish were maintained on a 14:10 h light:dark cycle at a temperature of 28°C and bred naturally. Transient expression of GCaMP2 and GCaMP5 were driven by constructs containing the Ribeye-A promoter (Dreosti et al., 2009), which targets ribbon synapses of sensory neurons. Constructs at a concentration of 0.084 μg/μl were injected into embryos at the 1–4 cell stage. In order to inhibit melanin formation, larvae were kept in 200 μM 1-phenyl-2-thiourea (Sigma) from 28 h post-fertilization. Imaging was carried out as described previously (Dreosti et al., 2009). Briefly, whole zebrafish larvae (9–11 dpf) were immobilized in 2.5 % low melting point agarose on a glass coverslip. The temperature of the room was 23–24°C. Bipolar cell terminals were imaged in vivo using a custom-built 2-photon microscope equipped with a mode-locked Chameleon titanium–sapphire laser tuned to 915 nm (Coherent) with an Olympus LUMPlanFI 40x water immersion objective (N.A. 0.8). Emitted fluorescence was captured through both the objective and a sub-stage oil condenser, filtered through a HQ 520/60 nm 2P GFP emission filter (Chroma Technology) and detected by a set of GaAsP photodiodes (Hamamatsu, Japan). Scanning and image acquisition were controlled under ScanImage v.3.6 software (Pologruto et al., 2003). Movies were processed using the SARFIA suite of analysis routines (Dorostkar et al., 2010) running in Igor Pro 7 (Wavemetrics). These routines begin with the automated extraction of the fluorescence change in each terminal by defining ROIs using a filtering algorithm based on a Laplacian transform.
Light stimuli were delivered using an amber LED (590 nm, Phillips Luxeon, 350 mA, 3 V) filtered through a 600/10 nm BP filter (ThorLabs) and projected through a light guide onto the surface of the bath, very close to eye of the fish. The stimulus time-course was controlled in Igor Pro v. 4.01 (WaveMetrics) and time locked to image acquisition through ScanImage. The mean intensity of light stimuli was ~2 × 105 photons/μm2/s, which corresponds to a low photopic intensity.
Mice were anesthetized using isoflurane (3% for induction, 1.5 – 2% during surgery). A circular craniotomy (2 – 3 mm diameter) was placed above V1 (centered 2.7 mm lateral from lambda suture). OGB-1-AM (Invitrogen) was injected as previously described (Stosiek et al., 2003; Komiyama et al., 2010; Zariwala et al., 2011 ). A custom titanium head post was fixed to the skull using black dental cement (Contemporary Ortho-Jet). The craniotomy was covered with agarose (1.2 – 1.5%) and a round glass coverslip (Warner Instrument, 5 mm diameter, #1 thickness) was cemented to the skull to reduce motion of the exposed brain. For cell-attached recording, the exposed brain was covered with ~1 mm thick agarose without a coverslip. The animal was then transferred to the imaging setup, where it was placed on a warm blanket (37 C) and kept anesthetized using 0.5% isoflurane and sedated with chlorprothixene (20 – 40 μl at 0.33 mg/ml, intramuscular)(Niell and Stryker, 2008b).
Imaging was performed with a custom-built two-photon microscope (designs available at http://research.janelia.org/Svoboda/). The light source was a Mai-Tai 80 femtosecond pulsed laser (Spectra-Physics) running at 900 nm for GCaMPs and 800 nm for OGB-1. The objective was a 40X dipping lens (Olympus, 40X, 0.8 NA). Image acquisition was performed using ScanImage 3.7 (http://www.scanimage.org)(Pologruto et al., 2003). Images (512 × 250 pixels, 250 × 250 μm) were collected at 4 Hz.
In vivo cell-attached recording was performed using glass pipettes (~5–7 MΩ) filled with solution containing (125 mM NaCl, 5 mM KCl, 10 mM Glucose, 10 mM HEPES, 2 mM CaCl2, 2 mM MgSO4, 0.1 mM AlexaFluor 594). Signals were amplified using an AxoPatch 200B amplifier (Molecular Devices), filtered at 5 kHz, and digitized at 10 kHz. Images (32 × 20 pixels, 20 × 20 μm) were acquired at 50 Hz.
The moving grating stimuli were generated using the Psychophysics toolbox (Brainard, 1997; Pelli, 1997) in MATLAB. Each stimulus trial consisted of a 4s blank period (uniform gray at mean luminance) followed by a 4s drifting sinusoidal grating (0.05 cycles per degree, 1 Hz temporal frequency). The visual stimuli were synchronized to individual image frames using frame-start pulses provided by ScanImage 3.7. The gratings were presented through an LCD monitor (30 cm × 40 cm), placed 25 cm in front of the center of the right eye of the mouse. The monitor subtended an angle of ± 38° horizontally and −20° to +38° vertically around the eye of the mouse. For cell-attached recording, we used a smaller LCD monitor (12 cm × 16 cm) placed 10 cm in front of the right eye. Each stimulus trial consisted of a 4s blank period followed by a 2s drifting square wave grating (0.05 cycles per degree, 2 Hz temporal frequency).
Slow drifts in brain position in the X and Y directions were corrected using the TurboReg plug-in in ImageJ (Thevenaz et al., 1998). All remaining analyses were performed in MATLAB. Regions of interest (ROIs) corresponding to visually identifiable cell bodies were selected using a semi-automated algorithm. For GCaMP3, GCaMP5G and 5K, ring-shaped ROIs were placed at the cytosolic regions of the cells (excluding the nucleus; GCaMP expression is typically restricted to the cytoplasm (Tian et al., 2009)). For OGB-1, circular ROIs covering the whole soma were used. The fluorescence time course of each cell was measured by averaging all pixels within the ROI. The neuropil contamination was corrected using a published method (Kerlin et al., 2010); in short, the neuropil signal Fneuropil (t) surrounding each cell was measured by averaging the signal of all pixels within a 20 μm circular region from the cell center (excluding all selected cells). The true fluorescence signal of a cell body was estimated as
with r = 0.7 throughout the study. After neuropil correction, the ΔF/F0 of each trial was calculated as (F − F0)/F0, where F0 is the baseline fluorescence signal averaged over a 2 s period immediately before the start of grating stimulation. Visually responsive neurons were defined using ANOVA across blank and eight direction periods (p < 0.01)(Ohki et al., 2005). Of the responsive cells, orientation selective cells were defined by ANOVA across eight direction periods (p < 0.01).
The orientation selectivity index (OSI), tuning width, and direction selectivity index (DSI) were calculated for visually responsive cells. First, the preferred orientation (θpref) of the cell was determined as the angle of the grating that produced the strongest response. The orientation-tuning curve was constructed by measuring the mean ΔF/F0 over the 4 s stimulus period for each orientation. We then fitted the tuning curve as the sum of two Gaussian functions centered on θ pref and θ pref + π with equal width σ, different amplitudes A1 and A2 and a constant baseline B (Niell and Stryker, 2008b). The value of σ was required to be larger than 15° to reflect the limit of our stimulus set (45° separation) in resolving sharper tuning. The OSI was defined as
where Rpref and Rortho are the response amplitude at the preferred (θpref) and the orthogonal orientation ( ). Tuning width was defined as the half width at half-maximum (HWHM) of the fitted Gaussian ( ). Finally, DSI was calculated as
where Ropposite is the response in the opposite direction (θpref +π).
For simultaneous imaging and cell-attached recording, ring shaped ROIs were placed at the cytosolic regions of the cells. Neuropil signals were measured from the image region with the cell excluded. Neuropil compensation was performed as above (r=0.7) before calculating ΔF/F. For visual evoked response (Figure 9a–e), the baseline fluorescence (F0) was measured over a 1 s period immediately before the start of grating stimulation. For responses to few isolated AP (Figure 9e,f), we searched the spike trace for events with 1, 2, and 3 APs within a 200 ms windows, and no other APs during a 1s period before and a 0.5s period after the first AP. The baseline fluorescence (F0) was measured over a 0.2s period prior to the first AP. AP detection was quantified using template-matching with the average trace of 1, 2, and 3AP events as templates for detecting 1, 2, and 3 APs respectively. The 0AP traces (478 traces) were taken from the same fluorescence recordings during periods where no AP was detected for at least 2s. The decision variable was the projection of the fluorescence traces along the direction of the template vector. Detection efficiency was defined as the fraction of correctly detected events given a 5% false positive rate.
The single exponential model fit of the GCaMP5 signal (Figure 9h) was calculated by linear convolution of the detected spike point process with a single exponential kernel: . The amplitude A and the decay time constant τ were adjusted to minimize the mean square error of the fit. For the non-linear model, the output of the linear mode x(t)= s(t) · h(t) was passed to a nonlinear stage to generate the final output:
The four parameters A, τ, b, c were adjusted independently to minimize the mean square error.
Unless specified otherwise, all statistical methods were implemented in MATLAB or Microsoft Excel, using standard packages.
A small family of GCaMP5 variants was produced from the GCaMP3 scaffold by combining improvements generated from site-directed mutagenesis at the cpGFP/CaM proto-interface and targeted library screening at the M13pep-cpGFP and cpGFP-CaM linkers. We also tested the effect of mutations to the M13 peptide and near the third Ca2+-binding site of CaM (Figure 1a, Table 1). To facilitate engineering, we solved the crystal structure of several GCaMP2 and GCaMP3 variants (Table 2).
In our previous analysis of the GCaMP2 structure (Akerboom et al., 2009), we predicted that the chemical environment and solvent accessibility of the cpGFP chromophore would be modulated in the Ca2+-bound state by the inter-lobe linker of CaM (residues 377–383 of GCaMP3). Here, we tested a small number of mutations and found that Asp380Tyr increases the Ca2+-bound fluorescence and dynamic range (ΔF/F)max of GCaMP3 (1.5-fold and 1.4-fold; Table 1; GCaMP3-Asp380Tyr = “GCaMP5A”, GCaMP3-Thr302Leu, Arg303Pro and Asp380Tyr=“GCaMP5G”). Comparison of the crystal structures of GCaMP2, GCaMP5A and GCaMP5G confirm that Tyr380 is oriented towards the cpGFP chromophore; indeed this mutation reorganizes nearby side-chains and the water network in close proximity to the cpGFP chromophore resulting in an apparent decrease in solvent accessibility of the cpGFP barrel (Figure 1b). Consequently, the Asp380Tyr mutation raises the brightness of the calcium-bound state of GCaMP3 for both GCaMP5A and GCaMP5G; in addition, calcium affinity and cooperativity (Hill coefficient) are increased by ~25% for GCaMP5A (Table 1).
The structure of the M13pep-to-cpGFP linker (“linker 1”) is known to be critical to sensor function (Nakai et al., 2001). Crystal structure analysis of Ca2+-bound and -free GCaMP2(Wang et al., 2008; Akerboom et al., 2009) revealed that both amino acids of linker 1 (Leu60 and Glu61) make significant contacts to cpGFP and CaM and stabilize both the apo and Ca2+-bound conformations. Upon screening mutagenic libraries at the linker 1 positions (59–60), we found a number of variants that substantially increased the dynamic range of GCaMP3 in E. coli lysates (Table 1, Figure 1c). Mutation of linker 1 from Leu-Glu (GCaMP3) to His-Pro (GCaMP5B) or Gln-Pro (GCaMP5C) roughly doubles or triples the (ΔF/F)max of GCaMP3, respectively (Table 1, Figure 1c).
Mutation of the cpGFP-to-CaM linker (“linker 2”) has also been shown to affect sensor function (Souslova et al., 2007). Accordingly, we also screened targeted libraries at the linker 2 positions (302–303) and identified variants with increased dynamic range in E. coli lysates (Table 1, Figure 1c). Most prominent was Leu-Pro (GCaMP5D), which has a (ΔF/F)max roughly twice that of GCaMP3, but with decreased affinity (Table 1).
Presuming that the improvements to sensor function resulting from the individual mutation of the cpGFP/CaM interface and linker 1 (L1) and linker 2 (L2) might be additive, different combinations of these (in the GCaMP3 scaffold, see Figure 1a) named GCaMP5E-J, were characterized extensively in vitro (Table 1). Compared to GCaMP3, all these GCaMP5 variants have significantly higher ΔF/F (Table 1). The Arg380Tyr mutation not only increases the brightness of both the saturated (sat) and apo states in the GCaMP3 background (GCaMP5A), it also increases the brightness of both states for several of the other sensors with modified linkers, e.g. 5G/D and 5H/F. Linker mutations L1-QP (GCaMP5C), L1-HP (GCaMP5B) and L2-LP (GCaMP5D) decrease the fluorescence of the apo state, and their combination has a cumulative effect: GCaMP5E, GCaMP5F, GCaMP5H and GCaMP5I are all exceptionally dim in the apo state. This effect is strongest for GCaMP5F and GCaMP5H, both containing combinations of L1-QP and L2-LP, leading to a ~160-fold fluorescence increase upon calcium binding in vitro (Table 1).
A GCaMP5 sensor incorporating the CaM mutation Arg90Gly (Sorensen and Shea, 1996) was constructed in the background of GCaMP5A, to further increase Ca2+ affinity (GCaMP numbering Arg392Gly, GCaMP5K). While GCaMP5K, as intended, has a higher affinity (and Hill coefficient) for Ca2+ than GCaMP3 (Kd = 190 nM, n = 3.8), it comes at a cost of (ΔF/F)max (Table 1). Another gain in affinity was achieved by mutation of the M13 peptide, Ala52Val, in a “hydrophobic anchor” position (Hultschig et al., 2004) and similar to the “D2” M13 variant (Palmer et al., 2006b). This mutation, acquired serendipitously in GCaMP5D (resulting in GCaMP5L), increases calcium affinity ~2-fold while preserving (ΔF/F)max (Table 1).
We used an in depth biophysical characterization of this panel of variants to reduce the number taken forward for in vivo testing, and identified some with particularly useful characteristics. GCaMP5G has a significantly higher apo pKa (Table 1), indicating that the Ca2+-free state is more likely to have a protonated, dim chromophore at physiological pH. Apparent rates of fluorescence conversion were determined by stopped-flow mixing, beginning at 0 [Ca2+] and stepping to 250 nM - 10 μM. Initial rates (k1) of fluorescence after steps to 10 μM [Ca2+] were significantly faster for 5D and 5G than for GCaMP3 (Table 1), although rates after steps to 500 nM [Ca2+] were comparable. Peak brightness and quantum yield of GCaMP5s tested were comparable to GCaMP3 (Table 3). One- and two-photon excited brightness spectra (Mutze et al., 2012) of equimolar (1 μM) amounts GCaMP3 and GCaMP5G revealed that a larger fraction of molecules is in the fluorescent, de-protonated state for GCaMP5G compared to GCaMP3, and ΔF/F for GCaMP5G was over two-fold improved under 1- and 2-photon excitation (Figure 1d). The GCaMP5 variants retain the slight red shift of GCaMP3 relative to EGFP (Figure 1d).
We next characterized the GCaMP5 variants in HEK293 cells in response to acetylcholine (Ach)-induced Ca2+ mobilization (Tian et al., 2009). The ratio of peak to baseline fluorescence was increased for most variants, by up to ~17 fold for GCaMP5H (Table 1). This is the direct result of lower baseline fluorescence for all variants except GCaMP5A, in agreement with in vitro data (Table 1). GCaMP5s containing L1-QP (GCaMP5C, GCaMP5F, GCaMP5H) were not selected for further analysis in cultured neurons, due to the low baseline brightness of variants containing this linker. GCaMP5K was directly tested in primary rat hippocampal neurons due to its high affinity; the HEK293 cell testing was omitted for this variant.
A selection of eight GCaMP5 variants was expressed in primary rat hippocampal neurons by lentivirus-mediated gene transfer. After 16–18 days, infected neurons appeared healthy (Figure 2a); most GCaMP5s were dimmer than GCaMP3 at baseline, except GCaMP5K (“5K”), which was approximately twice as bright before stimulation (Table 1). Action potentials were elicited at 83 Hz via a custom-built extracellular field stimulator (1 AP per stimulus, data not shown) and imaged at 34.8 Hz with an EMCCD camera. In parallel experiments we incubated neurons with cell-permeable versions (acetoxymethyl ether, “AM” derivatives) of the small molecule calcium dyes Oregon Green BAPTA-1 (OGB-1) or fluo-4. Robust fluorescence responses were seen from 1–160 field stimuli (Figure 2b–d). We extracted fluorescence from individual neurons and computed background-subtracted fluorescence transients. Single field stimuli, corresponding to one action potential, evoked fluorescence transients that were detectable in single trials (Figure 2b). Peak ΔF/F (Figure 2c) and signal-to-noise ratio (SNR, Figure 2d) were significantly improved for most GCaMP5 variants compared to GCaMP3, over at least part of the stimulus regime (Figure 2c,d; ΔF/F and SNR for each variant including GCaMP3, OGB-1, Fluo-4 and G-GECO1.2).
GCaMP5G (“5G”), which showed the highest response at maximum stimulation in the cultured neuron screen, was selected for testing in the context of specific protein fusions. Variant 5G was fused to synaptophysin (Dreosti et al., 2009) creating “SyGCaMP5G” (targeted to the outside of synaptic vesicles), and to the Lck domain (Shigetomi et al., 2010b) (“Lck-GCaMP5G”; inside of the plasma membrane). SyGCaMP constructs were transfected into hippocampal neurons, and boutons were imaged (Figure 3a); the response of SyGCaMP5G (“Sy5G”) to small field stimuli was significantly higher compared to SyGCaMP2 and SyGCaMP3 (Figure 3b).
5G also proved superior to GCaMP3 for detection of Ca2+ hotspots in astrocytes. The membrane-targeted Lck-GCaMP5G (“Lck-5G”) detected spontaneous Ca2+ transients (“spotty calcium” signals)(Shigetomi et al., 2010a) and ATP-induced responses with ~2-fold greater response magnitude than Lck-GCaMP3 in astrocytes; the 5G-determined responses also had more apparent local structure, significantly increasing the resolution of spotty calcium signal detection (Figure 3c–e). Furthermore, the threshold of detection was lowered from ~10 field stimuli for Lck-GCaMP3 to ~2 field stimuli for Lck-5G, indicating that both dynamic range and sensitivity have been substantially improved for membrane-targeted 5G compared to GCaMP3 (Figure 3f). SNR was slightly improved for most of the stimulus range, although variability was higher due to lower baseline fluorescence (Figure 3f).
In parallel, membrane-targeted 5G and GCaMP3 were expressed in astrocytes (under control of the glial fibrillary acidic protein (GFAP) promoter) co-cultured with neurons. Fluorescence changes in astrocytes following neural activity evoked by electrical field stimulation were easily recorded using GFAP-5G (Figure 3g). These results demonstrate that the GCaMP5 sensors will be useful to study Ca2+ dynamics in a range of cell types, and that these indicators are suited for membrane targeting and protein fusions.
GCaMP5D (“5D”) displayed large (ΔF/F)max to 1–3 field stimuli in cultured neurons (Figure 2d), and was selected for further testing in in vitro mouse retina, as before (Borghuis et al., 2011). We introduced 5D into adeno-associated virus (AAV) under the neuron-specific hsyn1 promoter and infected mouse retinas in vivo (Methods)(Borghuis et al., 2011). Retinal ganglion cells (RGCs) were bright and appeared healthy, with nuclear-excluded fluorescence (Figure 4a-top). We recorded light-evoked responses from 5D-expressing RGCs in dissociated retinas (Figure 4a-bottom). RGCs showed fluorescence changes both to the onset of the infrared scan laser (910 nm; 20 mW after the objective)(Borghuis et al., 2011) and to a flash of visible light (420 nm) delivered with a light-emitting diode (LED) two seconds after scan onset (Figure 4b). Time-course, amplitude and polarity of the responses varied across the labeled population, consistent with the known diversity of RGC types (e.g., ON and OFF, brisk and sustained). The peak amplitude (ΔF/F)max of the evoked fluorescence response was 1.6 ± 1.3 (n = 55; mean ± s.e.m.; range 0.0 – 4.5) (Figure 4c). Of all 5D-expressing cells, 65% responded with (ΔF/F)max > 0.6 (mean 2.3 ± 1.0). For comparison, GCaMP3 and OGB-1 produced distributions of (ΔF/F)max with peaks around 0.2, with only ~30% of GCaMP3- and OGB-labeled RGCs showing (ΔF/F)max > 0.15(Borghuis et al., 2011) (Figure 4c). SNR for GCaMP5D was significantly improved compared to GCaMP3 (GCaMP3: SNR average=8.97 ± 1.08, n = 187; GCaMP5D: SNR average 13.43 ± 1.16; n = 64; p = 0.025) (Figure 4d).
To test the in vivo performance of GCaMP5 in worms, we selected GCaMP3, 5A and 5G for analysis. We monitored the activity of the Caenorhabditis elegans AWCon chemosensory neuron, which responds to odor presentation with graded hyperpolarization and to odor removal with graded depolarization (Tian et al., 2009) (Figure 5a). All GCaMPs were expressed as before; in all cases, no defect was observed in AWC-dependent spontaneous omega turning frequency (Tian et al., 2009) (Figure 5b). Individual worms were imaged in a microfluidic chamber during an odor addition-removal sequence with the odor isoamyl alcohol (IAA) (Figure 5c,d). All GCaMPs detected the known decrease in AWC calcium upon odor addition and the calcium increase upon odor removal (Tian et al., 2009). 5G displayed a lower ΔF/F compared to GCaMP3 for odor presentation; interestingly the ΔF/F was increased compared to GCaMP3 for odor removal for this variant. SNR for both odor addition and removal was comparable to GCaMP3 for this variant. GCaMP5A showed the strongest suppression in baseline fluorescence upon odor addition, and the largest increase after odor removal. In addition, the higher baseline fluorescence of 5A resulted in a three-fold improvement in SNR for odor presentation, and a two-fold increase in SNR for odor removal, over GCaMP3. The 5A fluorescence plateaued near peak, consistent with its higher affinity and cooperativity of Ca2+ binding (Table 1). Many neurons in C. elegans exhibit graded responses to stimuli, and 5A extends the observable range to calcium fluctuations too small for GCaMP3 to detect reliably. These results indicate that 5A may be an improved indicator for C. elegans neurons, with the choice between 5A and GCaMP3 depending on the specific application.
Both GCaMP3 and 5G were expressed in Drosophila melanogaster larvae using a motor-neuron promoter (OK6-Gal4). Type 1b neuromuscular junction (NMJ) boutons of 3rd instar larvae were imaged following activity evoked from electrically stimulating motor neuron axons using a suction electrode (Macleod et al., 2002) (Figure 6a). Larval NMJs showed robust fluorescence changes to evoked APs across a range of stimulation frequencies. 5G showed a 3-fold increase in (ΔF/F)max over GCaMP3, saturating at ~8. Single action potentials were clearly detectable in trial-averaged 5G responses but only rarely in single trial responses (Figure 6b). The shapes of the ΔF/Fmax and SNRmax curves were similar for both indicators, but responses with GCaMP5G were significantly greater than those with GCaMP3 across all stimuli frequencies (n = 6–10 NMJs/GCaMP, 7 larvae, for electrical stimuli at 1, 5, 10, 20, 40, 80 and 160 Hz, p-values at ΔF/Fmax were 0.0006, 0.01, 0.002, 4.8e-06, 6.6294e-10, 1.5453e-11, 1.3136e-12 respectively and p-values at SNRmax were 6.6174e-06, 0.006, 0.001, 4.6563e-06, 1.2206e-09, 9.6765e-10, 1.2233e-06 respectively).
In parallel we crossed UAS:GCaMP Drosophila flies with a Gal4 fly line (GH146-Gal4) expressing broadly in the olfactory projection neurons, and imaged calcium changes in the antennal lobe (AL) (Jayaraman and Laurent, 2007), in the DC1 glomerulus, in response to presentations of octanol (Figure 6c). For 1% octanol (the highest intensity stimulus we presented) we observed: 5G, 8.96 ± 3.02 (ΔF/F)max (range 6.11–14.08, n = 6 ALs); GCaMP3, 3.18 ± 0.90 (ΔF/F)max (range 2.02–4.72, n = 6 ALs), a 3-fold improvement (Figure 6d) in dynamic range. We presented a range of concentrations of octanol to obtain a tuning curve with GCaMP3 and 5G (Figure 6e). The shapes of the ΔF/Fmax and SNRmax curves were similar for both indicators, but responses with GCaMP5G were significantly greater than those with GCaMP3 across all concentrations except 0.01 % octanol (n = 6 ALs/GCaMP, 5 flies, for concentrations of 0.0001%, 0.001%, 0.01%, 0.1% and 1% octanol, p-values for ΔF/Fmax were 0.007, 0.005, 0.07, 0.0003, 0.0002, 0.0003 respectively and p-values for SNRmax were 0.002, 0.001, 0.11, 0.004, 0.002, 0.001 respectively).
We tested GCaMP2, GCaMP3, 5A and 5G in vivo in zebrafish (Danio rerio), at three different locations in the visual pathway (Nevin et al., 2010). An overview of the fish visual system for each of the 3 imaging locations is shown (Figure 7a,c,g). GCaMP2 and 5G were first expressed in sensory neuron ribbon synapses, under control of the Ribeye-A promoter. In Ribeye-A GCaMP fish (Figure 7a), 5G showed a ~2.5 fold larger (Δ:F/F)max compared to GCaMP2 in bipolar cell terminals in the inner plexiform layer following a visual stimulus (Dreosti et al., 2009) (Figure 7b), which might be related to Ca2+ spiking in some of these terminals (Dreosti et al., 2011). GCaMP3, 5A and 5G were next expressed pan-neuronally, using the elavl3 (HuC) promoter. In RGC axons and tectal dendrites (Figure 7c), 5A and 5G responses to visual stimuli were greater compared to GCaMP3 (Figure 7d) in both average and maximum peak response (Figure 7e). SNR was improved for both 5A and 5G. Interestingly, as for the worm odor addition/removal imaging, 5A outperformed GCaMP3 and 5G (Figure 7d–f). In tectal neuron somata (Figure 7g), 5G (ΔF/F)max was over 3-fold improved compared to GCaMP3 (Figure 7h,i). More importantly, the number of cells with detectable visual responses increased by 6-fold (35% of the total number of cells) with 5G compared to GCaMP3 (Figure 7h). SNR was 4-fold improved for 5G compared to GCaMP3 (Figure 7j).
L2/3 neurons in the mouse primary visual cortex (V1) show a broad distribution of spike rates in response to visual stimulation (0–20 Hz; median 4 Hz)(Niell and Stryker, 2008a). In contrast to other sensory areas in the rodent (O’Connor et al., 2010), the majority of V1 neurons can be driven to spike. L2/3 in V1 thus provides an ideal system to assay the sensitivity of different calcium indicators in vivo. A larger fraction of responding neurons in V1 indicates higher sensitivity for detecting action potentials.
We thus tested GCaMPs in V1 under similar experimental conditions used in published studies (Niell and Stryker, 2008a; Kerlin et al., 2010) (Zariwala et al., 2012). 5G and GCaMP5K (“5K”) were tested in vivo because of their superb SNR in cultured neurons (Figure 2e,,4).4). Both GCaMP5s and GCaMP3 were delivered by AAV-hsyn1 viral infection. For comparison, the synthetic indicator OGB-1-AM was bulk-loaded into V1 (Stosiek et al., 2003; Ohki et al., 2005). Three weeks after AAV infection, robust GCaMP fluorescence was observed in layer 2/3 neurons (5G shown in Figure 8a). We imaged many in vivo fields of view in mice ~3–4 weeks after infection with GCaMP5G or GCaMP3. Nuclear filling, a correlate of cytomorbidity (Tian et al., 2009), was seen in only a small fraction of cells and was comparable to GCaMP3 (Figure 8a,b).
To elicit neuronal activity, moving gratings were presented in eight orientations to the contralateral eye of lightly anesthetized mice (Niell and Stryker, 2008a; Kerlin et al., 2010). Two-photon imaging revealed visual stimulus-evoked GCaMP responses in subsets of layer 2/3 neurons (5G shown in Figure 8c,d 438/2041 cells in 5 mice; ANOVA, p <0.01). Many of these cells were orientation- (e.g. cells 1, 2; Figure 8c) or direction-selective (cell 3; Figure 8c), consistent with the known properties of these neurons (Mangini and Pearlman, 1980; Sohya et al., 2007; Niell and Stryker, 2008a; Kerlin et al., 2010; Zariwala et al., 2011). Of the responsive neurons, 56% (244/438) were orientation-selective (ANOVA across 8 orientation, p < 0.01), and 24% (105/438) were direction selective (DSI>0.5, Figure 8d). For 5G, the average orientation selectivity index (OSI) of visually responsive neurons was 0.75±0.24 (mean ± s.d); the direction selectivity index (DSI) averaged 0.32 ± 0.26. The orientation/direction selectivity was comparable to published reports based on electrophysiology (Niell and Stryker, 2008a), suggesting that expression of 5G does not compromise the tuning properties of mouse cortical neurons.
Identical experiments were performed with GCaMP3 (8 mice, 3910 cells) and OGB-1 (3 mice, 3606 cells). 5G and 5K allowed the detection of visual responses in more than twice as many neurons as GCaMP3, although sensitivity still lagged OGB-1 (Figure 8e). This improvement did not depend on compensation for the neuropil signal (Figure 8f; Methods)(Kerr et al., 2005). The fluorescence half-decay time after stimulus offset was not significantly different between 5G, 5K and GCaMP3; all three were significantly faster than OGB-1 (Figure 8g,i). The average 5G response amplitude (peak ΔF/F) at the preferred orientation was ~2-fold greater than GCaMP3 among low responders (50th–80th percentile of all neurons), mid responders (80th–97th percentile) and high responders (>97th percentile); for 5K the low and mid responders were ~2-fold greater, whereas the high responders were comparable to GCaMP3 (Figure 8h). This indicates both an improved sensitivity and dynamic range for both GCaMP5s compared to GCaMP3. Average GCaMP5G and OGB-1 OSI, DSI, and tuning width were almost indistinguishable (Figure 8j–l) and similar to published studies based on electrophysiology. These data show that long-term expression of GCaMPs does not perturb the synaptic circuits underlying orientation and direction tuning.
To characterize the relationship between spiking and GCaMP signals in vivo, we performed simultaneous cell-attached recording (Sato et al., 2007; Tian et al., 2009) and Ca2+ imaging during visual stimulation (Figure 9a) for GCaMP5K. The goal of these experiments was two-fold. First, we wanted to directly relate the fluorescence change of one of the most promising GECIs, GCaMP5K, to spiking activity. Second, in our in vivo experiments using GCaMP3 we have noticed large trial-to-trial fluctuations in the fluorescence signal (Huber et al., 2012). We thus tried to determine if the variability in fluorescence is explained by variability in neural activity coupled to the non-linear dynamics of the sensor, or if other noise sources need to be considered.
Visually evoked spikes were tightly correlated with increases in 5K fluorescence (Figure 9b). Single spike-induced fluorescence events were occasionally detected (e.g. Figure 9a, arrow). In simple cells, the phasic modulation of spike rate at the temporal frequency of the drifting grating (2 Hz in this set of experiments) was reflected in a step-like change in the fluorescence signal (Figure 9a). The orientation selectivity of the spiking responses was reflected in the 5K responses (Figure 9b).
The relationship between fluorescence change and spike rate within the stimulus period was supra-linear (Figure 9c–e). This supra-linear relationship was also evident for spontaneously occurring spikes (Figure 9g). The amplitude of 5K signal for single APs and bursts of two or three APs was 3.6±5.7%, 9.1±7.7% and 18.5±9.7% (mean±s.d.), respectively. For GCaMP3, the fluorescence response to 1–3 physiological APs was barely detectable in vivo (Tian et al., 2009). The half rise time and half decay time for single AP-induced signals were 28±5 s and 268±20s, respectively (mean±s.e.m.). We quantified the spike detection efficiency of 5K under our imaging conditions. The detection efficiency was 29.3% for 1 AP, 63.0% for 2 APs and 95.2% for 3 APs, at a 5% false positive rate (Figure 9f).
The supra-linear relationship between spiking and GCaMP5K fluorescence has to be taken into account in the interpretation of GCaMP-based calcium imaging. This supra-linearity provides superb SNR for imaging highly active neurons, however the supra-linear relationship also emphasizes differences in spike rate across trials, leading to a larger trial-to-trial variability. Indeed, trial-to-trial variability of GCaMP5K responses was larger compared to spikes (Figure 9i) and OGB-1 responses (Figure 9j) during repeated presentation of preferred stimuli.
This difference in variability is completely explained by a quantitative model relating 5K signals and spiking activity. For synthetic calcium indicators convolving spikes with a single exponential kernel accounts well for the signal of several synthetic calcium indicators (Yaksi and Friedrich, 2006; Greenberg and Kerr, 2009; Kerlin et al., 2010; Komiyama et al., 2010; Sato and Svoboda, 2010). When applied to GCaMP5K, however, a single exponential model underestimates the response during strong activity (arrows, Figure 9h), and over-estimates the response during weak activity (arrowheads, Figure 9h). Adding a simple supra-linearity after the linear convolution (see Methods) greatly improves the fit (Figure 9h, green curve). The nonlinear model reproduces the higher trial-to-trial variation of the GCaMP5K response compared to spikes (Figure 9i). Furthermore, removing the non-linear fit from the data, and replacing it with a linear fit, brought the trial-to-trial variability back to a similar level as estimated by spikes.
GCaMP3 has been widely used in diverse model organisms, facilitating a large number of new neuroscience applications. However, GCaMP3 has remained lacking in detection of sparse spiking activity. All of the GCaMP5s described showed improved dynamic range compared to GCaMP3. Several GCaMP5 indicators (5A, 5D, 5G and 5K) were further characterized in a wide variety of neuronal and astrocyte activity imaging preparations, and in several model organisms in vivo. Depending on the application, the user will have the option to pick the ideal GECI.
GCaMP5s were produced from the GCaMP3 scaffold by targeted library screening at the cpGFP/CaM proto-interface and the two inter-domain linkers. Functional analysis of GCaMP variants has largely validated the conjecture that crystal structures provide snapshots of conformational states that are directly relevant to the sensing mechanism. The primary utility of GCaMP crystal structures has been to delineate the protein domain proto-interfaces and the GFP chromophore environment, which guided targeted library mutagenesis, while also supporting some atomic-level predictions. The critical importance of linker length and composition to GCaMP function is consistent with recent results in other sensor classes (Horikawa et al., 2010; Alicea et al., 2011; Marvin et al., 2011). In GCaMP and other sensors based on circularly permuted fluorescent proteins, the linkers are in close proximity to the FP chromophore, and may both directly modulate fluorescence and contribute to the efficiency of larger-scale ligand-dependent domain rearrangements. We have consistently found proline to be enriched in the inter-domain linkers of high-SNR variants of cpGFP-based indicators (this work; (Alicea et al., 2011; Marvin et al., 2011)); the increased rigidity of the proline polypeptide backbone may limit inter-domain conformational sampling to states that differ more widely in fluorescence.
Some GCaMP5 variants show a greater than 150-fold increase in fluorescence upon calcium binding in vitro; this is roughly equivalent to the best small-molecule probes (but both have low baseline fluorescence, however). However, this ten-fold improvement did not translate into increased performance for these specific GCaMPs in more demanding preparations such as cultured neurons, underlining the importance of these medium-throughput cell-based assays in biosensor development. Performance in various in vivo and reduced preparations in a variety of cell types shows that GCaMP5A, GCaMP5G and GCaMP5K variants consistently outperform GCaMP3. This level of improvement was also consistent for GCaMP5G fused to synaptic vesicles and the plasma membrane. These results suggests that the improvements are “intrinsic” to the sensor, rather than due to “extrinsic” factors such as [Ca2+] levels in specific cells or sub-cellular locations, temperature, expression level, or vagaries of the particular system tested. Biophysical characterization showed that several factors contribute to this improvement: lower Ca2+-free fluorescence, higher Ca2+-bound fluorescence, and higher Ca2+ affinity.
GCaMP5G showed the largest responses to maximal stimulation when expressed in cultured neurons, and was therefore tested in all model organisms described. Comparative experiments with GCaMP5A and GCaMP5G in zebrafish and worm, and GCaMP5G and GCaMP5K in mouse, showed that all these GCaMP5s outperformed GCaMP3 in ΔF/F and/or SNR. In worm and zebrafish, GCaMP5A showed the largest (ΔF/F)max and SNR, indicating that for some preparations GCaMP5A might be the preferred version over GCaMP5G and GCaMP3. In mouse, GCaMP5G and GCaMP5K performed similarly well.
A recent report describes variants of GCaMP3, termed “G-GECOs” (Zhao et al., 2011), optimized by selecting for maximum ΔF/F in E. coli colonies expressing random mutagenic libraries. Each of the G-GECO variants described is significantly dimmer than GCaMP3 in both the Ca2+-free and -bound states, with a greater decrease in the former leading to a higher ΔF/F. Although G-GECO1 displayed improved KCl-evoked signal change in dissociated rat hippocampal neurons at maximum stimulation, an intrinsically dimmer sensor may complicate imaging in more complex preparations. Indeed, SNR of G-GECO1.0, G-GECO1.1 and G-GECO1.2 was lower than GCaMP3 over the complete range of field stimuli (Figure 2). In vitro (ΔF/F)max for the best G-GECOs are ~2-fold higher than GCaMP3, whereas some GCaMP5 variants show a 14-fold increase in the (ΔF/F)max of GCaMP3 in vitro. GCaMP5G, the most consistently high-performing variant across in vivo assays, outperforms G-GECOs in terms of (ΔF/F)max in vitro. The improved performance of the GCaMP5s versus the G-GECOs supports the strategy of structure-guided engineering as an efficient way to improve sensors, rather than random mutagenesis.
We compared the in vivo responses of mouse visual cortical neurons labeled with two of the best GCaMP5s (GCaMP5G and GCaMP5K), GCaMP3 and OGB-1, upon presentation of visual stimuli to the mouse. The GCaMP5 variants and OGB-1 showed similar fluorescence responses for most cells at the optimal stimulus orientation; GCaMP3 fluorescence responses were significantly lower. For all GCaMPs, after stimulation, fluorescence intensity returned to baseline level significantly faster compared to OGB-1. Annotation of the total fraction of visually responsive cells by GCaMP5 was two-fold improved over GCaMP, although still trailing OGB-1. This is in agreement with the zebrafish imaging, where GCaMP5G resulted in a larger fraction of visually responsive cells as well. These results are consistent with simultaneous imaging/cell-attached recordings; GCaMP5K detects single APs and bursts of 2–3 APs much better than GCaMP3 (Tian et al., 2009), although less well than OGB-1. Detection of larger bursts of activity is also much improved for GCaMP5s compared to GCaMP3 and OGB-1, increasing the effective dynamic range of imaging. A non-linear model of fluorescence dependence on spike rate accurately fits the in vivo data, and reduced trial-to-trial variability. Together these results show that GECIs are approaching small molecule indicators in terms of detection of sparse activity and neural activity quantification.
With its improved performance, GCaMP5 will directly enable more experiments in neuroscience and other fields of biology. Improvements in GECI transgene delivery will also contribute to increased utility. We have recently published a Cre-dependent GCaMP3 reporter mouse, allowing stable long-term expression in genetically defined neurons (Zariwala et al., 2012). Trans-synaptic delivery of calcium indicators is possible using rabies virus (Osakada et al., 2011). Zinc finger nucleases have facilitated chromosomal knock-ins in a variety of organisms, including those with few established genetic resources.
Although the GCaMP5s are currently the best GECIs for single-wavelength calcium monitoring, FRET-based sensors offer the advantage of easy donor/acceptor ratioing, primarily for motion artifact control. Several versions of the FRET-based GECI Yellow Cameleon, including YC2.6, YC3.6 (Nagai et al., 2004b), and the high-affinity YC-Nano (Horikawa et al., 2010), have been reported to detect sparse neural activity in various in vivo preparations (Grewe et al., 2010). Additionally, mutants of GCaMP3, GEX-GECO1 and GEM-GECO1(Zhao et al., 2011), offer excitation- and emission-based ratioing from a single FP chromophore. Variants of GCaMP5 harboring similar mutations, fusion to a second FP, or stoichiometric expression using viral 2A peptides may offer increased performance levels with a ratiometric output.
Although we have shown that the GCaMP5 variants constitute a significant improvement over G-GECO and GCaMP3, further GECI engineering remains. Detection of sparse spiking activity should be brought in line with the best small molecule indicators. Improvements in rise and decay kinetics are required to precisely monitor spike number and time. Long-term over-expression artifacts, such as the “cytomorbid” nuclear-filling phenotype, must be understood and eliminated, either through protein engineering or fine-tuned control of expression by promoter and enhancer adaptation. Improved GECIs, in combination with recent advances in light delivery and collection, fast scanning, image analysis, and behavioral paradigms, are setting the stage for chronic neural activity imaging to address fundamental questions in learning and memory, development, and the neural basis of behavior.
We would like to thank the Janelia Farm Shared Resources for critical support of this project: A. Arnold for imaging; A. Hu and B. Shields for histology; M. Ramirez, K. McGowan and the MBSR for molecular biology and virus production; H. White, D. Murphy and the CCSR for cell culture; Instrumentation Design & Fabrication for technical development and assistance; T. Laverty and the Fly Core for fly crosses and stock maintenance; R. Simmons, R. Rogall and H. Naz for media prep; and C. Kimmel and J. Cox for mouse support. We thank R. Behnam and A. Karpova for help in establishing lentiviral production. We also thank D. Rinberg and T. Tabachnik for olfactometer design and construction, and S. Hampel and J. Simpson for cloning. Use of the Advanced Photon Source at Argonne National Laboratory was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. Use of the Lilly Research Laboratories Collaborative Access Team (LRL-CAT) beamline at Sector 31 of the Advanced Photon Source was provided by Eli Lilly Company, which operates the facility. We thank Greg MacLeod for generous assistance in establishing the Drosophila NMJ preparation. H.B. and A.F. would like to thank J. Donovan, T. Thiele, and A. Tran for technical help, and C. Niell for helping with data analysis. A.F. has been supported by an EMBO long-term fellowship. Work in HB’s laboratory was funded by the NIH Nanomedicine Development Center “Optical Control of Biological Function”, work in SW’s laboratory was funded by NIH R01 NS045193.
Author contributionsL. L. L., J. A., J. S. M., L. T., E. R. S., design of project; J. A., J. S. M., E. R. S., L. T., L. L. L., sensor design; J. A. and J. S. M., E. coli screen; J. A. and L. T., HEK293 cell screen; J. A. and J. J. M., 2-photon spectroscopy; J. A. and E. R. S., X-ray crystallography; J. A., E. R. S., L. L. L., L. T., J. S. M., structure analysis; L. T., T-W. C., T. J. W., D. S. K., neuron screen; N. C. C., L. L., L.T., SyGCaMP & mGCaMP fusion screen; X. R. S. and S. S.-H. W., kinetics measurements; R. T. and B. S. K., Lck fusion astrocyte imaging; L. T., astrocyte/neuron imaging; T. J. W., Drosophila larval assay; T. J. W., A.A., V. J., adult Drosophila imaging; B. G. B., mouse retinal imaging; A. G. and C. B., C. elegans imaging; N. C. C., F. Esposti, L. L., zebrafish retinal bipolar cell imaging; M. B. O., R. P., F. Engert, zebrafish 2-photon optic tectum imaging; H. B., A. F., zebrafish tectal neuropil imaging; T-W. C., S. M., K. S., D. S. K., mouse in vivo imaging; T.-W. C., mouse in vivo electrophysiology; L. L. L., J. A., E. R. S., led the project; L. L. L., J. A., E. R. S., wrote the paper, with the help of the other authors.
All experiments were conducted according to protocols approved by the Institutional Animal Care & Use and Institutional Biosafety Committees of the Howard Hughes Medical Institute, Janelia Farm Research Campus, and of the corresponding committees at the other institutions.
Competing interests statement
The authors declare no competing financial interests.
Constructs and AAV viruses for GCaMP5G are available from AddGene (http://www.addgene.org) or the University of Pennsylvania Vector Core (http://www.med.upenn.edu/gtp/vector_core.shtml), respectively. The Janelia Farm GECI Project website (http://www.janelia.org/team-project/geci-project) has additional information about construct characterization and availability.