Genetically encoded calcium indicators (GECIs) produce unprecedentedly large signals that have enabled routine optical recording of single neuron activity in vivo in rodent brain. Genetically encoded voltage indicators (GEVIs) offer a more direct measure of neuronal electrical status, however the signal-to-noise characteristics and signal polarity of the probes developed to date have precluded routine use in vivo. We applied directed evolution to target modulable areas of the fluorescent protein in GEVI ArcLight to create the first GFP-based GEVI (Marina) that exhibits a ΔF/ΔV with a positive slope relationship. We found that only three rounds of site-directed mutagenesis produced a family of "brightening" GEVIs with voltage sensitivities comparable to that seen in the parent probe ArcLight. This shift in signal polarity is an essential first step to producing voltage indicators with signal-to-noise characteristics comparable to GECIs to support widespread use in vivo.
Genetically encoded calcium indicators (GECIs) produce unprecedentedly large signals that have enabled routine optical recording of single neuron activity in vivo in rodent brain. Genetically encoded voltage indicators (GEVIs) offer a more direct measure of neuronal electrical status, however the signal-to-noise characteristics and signal polarity of the probes developed to date have precluded routine use in vivo. We applied directed evolution to target modulable areas of the fluorescent protein in GEVI ArcLight to create the first GFP-based GEVI (Marina) that exhibits a ΔF/ΔV with a positive slope relationship. We found that only three rounds of site-directed mutagenesis produced a family of "brightening" GEVIs with voltage sensitivities comparable to that seen in the parent probe ArcLight. This shift in signal polarity is an essential first step to producing voltage indicators with signal-to-noise characteristics comparable to GECIs to support widespread use in vivo.
Entities:
Keywords:
ArcLight; Genetically encoded voltage indicators; biosensors; directed evolution; green fluorescent protein; mutagenesis
High sensitivity of
genetically encoded calcium indicators (e.g.,
GCaMPs) combined with high resolution two-photon microscopy have provided
unprecedented views of neuronal activity in vivo.[1] This methodology enables monitoring of the activity of
a large number of genetically defined neurons in the cortex of awake
rodents. However, intracellular calcium levels are a surrogate indicator
of the more relevant physiological measure, membrane potential. An
alternative to calcium indicators for directly reporting neuronal
activity are genetically encoded indicators of membrane potential.
However, while the first fluorescent-protein voltage indicator was
reported nearly two decades ago,[2] and many
additional constructs have been reported,[3−16] mainly due to low voltage sensitivity none have seen general adoption
for use in physiological studies.We developed ArcLight,[10] a genetically
encoded fluorescent voltage indicator based on fusion of the voltage
sensing domain (VSD) of the Ciona intestinalis voltage
sensitive phosphatase[17] and a mutated version
of the pHluorin fluorescent protein.[18] Due
to its relatively high sensitivity, ArcLight has been used to study
the electrical activity of Drosophila neurons[19−22] in situ, Caenorhabditis elegans neurons in situ,[23] human cardiomyocytes,[24,25] and mammalian neurons in culture[10,16,26] and in vivo.[27] The successful
use of ArcLight in various organisms is mainly due to its brightness
and good membrane localization. Unlike opsin-derived indicators,[9,13,14] ArcLight is based on the Aequorea victoria green fluorescent protein (GFP) which
has brightness that is orders of magnitude greater than opsins.[28] And while there are a number of GEVIs published
based on the fusion of a voltage sensing domain with GFP-like fluorescent
proteins[2−8,11,12,16,29−31] selective mutations in the ArcLight fluorescent protein have produced
a construct with unparalleled changes in absolute fluorescence output
with voltage. Nonetheless, relatively slow kinetics (activation time
τ > 9 ms) of the voltage-dependent fluorescence change in
ArcLight
significantly diminishes its performance in detecting action potentials
in neuronal cells and has negatively impacted its broader use in physiological
experiments.[16,32]The current multiphoton
methods of recording neuronal activity
in the mammalian cortex cannot follow the response kinetics of faster
indicators when recording the concurrent dynamics of large populations
of individual neurons in awake animals.[33,34] Modeling studies
clearly indicate that temporally exaggerated fluorescence responses
to action potentials will improve detection.[34,35] In the light of this, the most important limitations of ArcLight
and other voltage probes compared to GFP-based calcium indicators
(i.e., GCaMP6) are not response speeds but rather: (i) the absolute
(and relative) changes in their fluorescence output are far lower
than the state-of-the-art calcium indicators and (ii) the voltage
indicators with the largest signals combined with optimal brightness
(i.e., ArcLight, ASAP2f, and AceNeon) have a negative slope relationship
between voltage and fluorescence, unlike calcium indicators that have
positive slope relationships between calcium levels and fluorescence.[10,32,36]To improve the utility
of the GEVI ArcLight for in vivo recordings
of neuronal electrical activity, we sought to reverse the polarity
of the signal from a decrease in fluorescence intensity with depolarization,
to an increase. However, numerous studies on ArcLight failed to reveal
a clear mechanistic hypothesis for its robust voltage sensitivity.[37−39] The majority of ArcLight modifications targeting the voltage sensitive
domain resulted in several probes with improved kinetics but reduced
amplitude of the voltage-dependent response.[26,38,39] In addition, the vast majority of single
(non FRET) GFP-based voltage indicators, both published (FlaSh, VSFP3,
SPARC, Zahra1 and 2, ElectricPK, ArcLight, other species ArcLights,
Bongwoori, ASAP1, ASAP2f, ArcLightning, etc.) and unpublished (>500,
Platisa and Pieribone), produce fluorescence decreases in response
to depolarizations. To our knowledge, there is only one GFP-like voltage
indicator (FlicR1) that exhibits a depolarization-dependent increase
in fluorescence intensity substantial enough to be considered for
practical use.[40] Herein, to produce such
a probe, we undertook a directed evolutionary mutagenic strategy targeting
the fluorescent protein in ArcLight combined with automated screening.
We show that changes to only three residues in the fluorescent protein
(FP) of ArcLight resulted in a family of “brightening”
GEVIs, that exhibit comparable voltage sensitivity to ArcLight. The
most sensitive variant of these novel probes, the GEVI Marina is capable
of reporting action potentials and subthreshold events in spontaneously
active mammalian neurons in vitro.
Results and Discussion
The most advanced FP-based genetically encoded voltage indicators
(i.e., ArcLight, ASAP2f, Ace-mNeon) reduce their fluorescence intensity
with activation (i.e., depolarization).[10,32,36] As such, during in vivo recording, this bright resting
state fluorescence (F0) introduces a large
noise signal from surrounding cells that reduces the overall signal-to-noise
ratio (SNR) of the response compared to probes that exhibit low fluorescence
in the resting state. As neurons exist in a depolarized state far
more infrequently than in a resting (hyperpolarized) state, an indicator
that exhibits low resting-state fluorescence would yield improved
SNR. A key improvement of the later versions of GCaMPs (>GCaMP3)
is
their low fluorescence in neurons at rest. Action potentials cause
an influx of calcium and thereby a dramatic increase in the fluorescence
output of GCaMP6. The high ΔF/ΔV sensitivity of ArcLight is conferred by the presence of
specific residues in the fluorescent protein. We reasoned that modifications
to the fluorescent protein might result in a change in polarity of
the voltage-dependent signal.For the development of a GEVI
with a positive slope relationship
between membrane potential and fluorescence, we decided to use the
brighter version of ArcLight, super ecliptic ArcLight as a starting
point. We previously found that eGFP could replace the ArcLight fluorescent
protein (super ecliptic pHluorin) if four crucial mutations were made
(S389D, S444F, Q446T, and A469D).[37] Notably,
all of these residues populate the same surface of the β barrel
(beta sheets 7, 10 and 11) and have outward orientations. We began
by making saturated mutagenesis libraries of three of these crucial
residues (D389, F444 and T446) in the FP domain of ArcLight (Figures , S1, and 2). In addition, we made mutagenic libraries of six
other residues (R338, H390, Y442, T445, S447, and E464; Figures and S1), which are proximal to either the chromophore, or residues that
we found to be essential in producing the ArcLight effect.
Figure 1
Crucial residues
within the ArcLight fluorescent protein (super
ecliptic pHluorin) determine the polarity and size of the voltage-dependent
fluorescence response in genetically encoded voltage indicators. Two
different view angles of the 3D structure of eGFP (PDB identifier 1EMA) illustrating the
position of residues within super ecliptic pHluorin of ArcLight that
confer the inverted polarity and high sensitivity of the voltage indicator
Marina: D389A (blue), H390A (magenta) and Y442V (yellow). The crucial
residues for voltage sensitivity of ArcLight are D389 (blue), F444
(orange), T446 (light blue), and D469 (green).[37]
Figure 2
Directed evolution of GEVI ArcLight produces
an indicator with
reversed voltage-dependent response. Left side panels (A, C, and E):
for each of the saturated mutagenesis libraries, 96 colonies were
selected and tested for expression and voltage sensitivity. Screening
of each library was preformed two to four times on a custom-made,
semiautomated screening platform which allows for simultaneous field
stimulation and optical imaging in the 96-well dish format. We used
either an electrically active HEK293 cells or mouse cortical neurons
transiently expressing various mutants. Each dot on the graph represents
the highest value of depolarization-dependent fluorescence change
detected in each well. Right side-panels (B, D, and F): detailed study
of voltage sensitivity of selected indicators using combination of
whole-cell patch-clamp recordings and optical imaging in nonelectrically
active HEK293 cells. Excitation light intensity measured at the sample
plane was 18 mW/mm2. For each construct three to six cells
were tested. Data are plotted as mean ± SEM. (A) First round
of site-directed mutagenesis on residues proximal to the chromophore
(R338, H390, Y442, T445, S447, and E464) or residues which we found
to be essential for ArcLight voltage sensitivity (D389, F444, and
T446). F444X and T446X libraries were tested in neurons, all other
libraries were tested in spiking HEK293 cells. Some of the constructs
carrying mutations on residue D389 produce probes that respond to
membrane depolarization with modest increase in fluorescence intensity.
(B) Detailed study confirms an inverted signal polarity of several
D389 mutants compared to the parent ArcLight. The highest positive
depolarization-dependent signal was detected in an ArcLight D389A
mutant. (C) Second round of site-directed mutagenesis on residues
R338, H390, T445, S447, and E464 using ArcLight D389A as a template
identified several H390 mutants with increased voltage sensitivity
compared to parent indicator. (D) Voltage clamp recordings identified
mutant ArcLight A389 H390A as a high responder. (E) Additional mutagenic
library at residue Y442 of ArcLight A389 H390A produces variants with
further improved voltage sensitivity. (F) Voltage clamp analysis confirms
significant improvement in voltage sensitivity in two of Y442 mutants,
Y442L and Y442V. ArcLight A389 A390 Y442V was named Marina.
Crucial residues
within the ArcLight fluorescent protein (super
ecliptic pHluorin) determine the polarity and size of the voltage-dependent
fluorescence response in genetically encoded voltage indicators. Two
different view angles of the 3D structure of eGFP (PDB identifier 1EMA) illustrating the
position of residues within super ecliptic pHluorin of ArcLight that
confer the inverted polarity and high sensitivity of the voltage indicator
Marina: D389A (blue), H390A (magenta) and Y442V (yellow). The crucial
residues for voltage sensitivity of ArcLight are D389 (blue), F444
(orange), T446 (light blue), and D469 (green).[37]Directed evolution of GEVI ArcLight produces
an indicator with
reversed voltage-dependent response. Left side panels (A, C, and E):
for each of the saturated mutagenesis libraries, 96 colonies were
selected and tested for expression and voltage sensitivity. Screening
of each library was preformed two to four times on a custom-made,
semiautomated screening platform which allows for simultaneous field
stimulation and optical imaging in the 96-well dish format. We used
either an electrically active HEK293 cells or mouse cortical neurons
transiently expressing various mutants. Each dot on the graph represents
the highest value of depolarization-dependent fluorescence change
detected in each well. Right side-panels (B, D, and F): detailed study
of voltage sensitivity of selected indicators using combination of
whole-cell patch-clamp recordings and optical imaging in nonelectrically
active HEK293 cells. Excitation light intensity measured at the sample
plane was 18 mW/mm2. For each construct three to six cells
were tested. Data are plotted as mean ± SEM. (A) First round
of site-directed mutagenesis on residues proximal to the chromophore
(R338, H390, Y442, T445, S447, and E464) or residues which we found
to be essential for ArcLight voltage sensitivity (D389, F444, and
T446). F444X and T446X libraries were tested in neurons, all other
libraries were tested in spiking HEK293 cells. Some of the constructs
carrying mutations on residue D389 produce probes that respond to
membrane depolarization with modest increase in fluorescence intensity.
(B) Detailed study confirms an inverted signal polarity of several
D389 mutants compared to the parent ArcLight. The highest positive
depolarization-dependent signal was detected in an ArcLight D389A
mutant. (C) Second round of site-directed mutagenesis on residues
R338, H390, T445, S447, and E464 using ArcLight D389A as a template
identified several H390 mutants with increased voltage sensitivity
compared to parent indicator. (D) Voltage clamp recordings identified
mutant ArcLight A389 H390A as a high responder. (E) Additional mutagenic
library at residue Y442 of ArcLight A389 H390A produces variants with
further improved voltage sensitivity. (F) Voltage clamp analysis confirms
significant improvement in voltage sensitivity in two of Y442 mutants,
Y442L and Y442V. ArcLight A389 A390 Y442V was named Marina.Ninety six (96) Escherichia
coli colonies from
each site directed mutagenesis library were selected, and cDNA from
each colony was purified. Each of the 96 clones was tested by screening
two to four times in either primary neuronal cell cultures or in spontaneously
spiking HEK cell line[41] using a custom-made
automated screening platform (Vasan, Platisa, and Pieribone, personal
communication). Mutations to residues R338, H390, Y442, F444, T445,
T446, S447, and E464 produced signal sizes that were either unchanged
or smaller than those found in the wildtype ArcLight (Figure A). However, mutations to D389
caused either significantly reduced negative-going fluorescence changes
or in some cases, increase in fluorescence intensity (up to 7%ΔF/F using our standard stimulation protocol;
see Methods) with depolarizations in membrane
voltage. We sequenced all 96 D389 mutants and selected constructs
that represented all 19 different novel amino acids. Each mutant was
then tested via whole cell patch clamp electrophysiology combined
with microfluorimetry in nonelectrically active HEK293 cells. The
results of these screens can be seen in Figure B. The mutant D389A produced the largest
increase in fluorescence intensity upon membrane depolarizations (+7.1
± 0.5%ΔF/F for +100 mV
depolarization from a holding potential of −70 mV in expressing
HEK293 cells). The D389W and D389C substitutions also produced positive,
albeit smaller, responses to depolarizations (Figure B). While the D389A mutation in ArcLight
reversed the sign of the fluorescence response to depolarization from
negative to positive, the absolute change in fluorescence as well
as the ΔF/F was relatively
small.In an effort to increase voltage sensitivity, we next
performed
saturation mutagenesis of the same residues as described above (R338,
H390, Y442, T445, S447, and E464) but using ArcLight A389 as a template.
Again, mutations of R338, T445, and E464 all produced constructs with
positive going signals but these were largely similar or smaller in
magnitude than ArcLight A389 (Figure C and D). However, constructs with mutations at H390
and S447 produced signals with significantly larger response magnitudes
than ArcLight A389. We tested each of the 19 amino acid substitutions
at H390 using whole cell patch clamp electrophysiology and found that
mutants carrying either alanine or cysteine at residue H390 produced
response magnitudes that were 2.5 and 2 times larger than ArcLight
A389, respectively (17.8 ± 1.2%ΔF/F and 13.9 ± 1.2%ΔF/F for a + 100 mV depolarization tested in HEK293 cells).
The mutant ArcLight A389 A390 showed biexponential kinetics with time
constants of ∼3.6 ms and ∼30 ms for activation (fast
component accounting for 38% of response), and ∼12 ms and ∼47
ms for inactivation (fast component accounting for 68% of response)
(Figure S2). Our attempt to increase voltage
sensitivity of the ArcLight W389 mutant by making changes at position
390 failed to produce improved voltage sensitive variants. We did
not study this mutant further.The following round of site-directed
mutations was performed on
residue Y442 in ArcLight A389 A390 (Figure E) and resulted in several variants (Y442
L, V, A, M and C) that had increased depolarization dependent changes
in fluorescence intensity (Figure F). The two most promising constructs containing leucine
or valine at position 442. These mutants exhibited an ∼1.6
fold increase in signal magnitude compared to the parent ArcLight
A389 A390 (31.1 ± 4.2%ΔF/F and 29.2 ± 2.0%ΔF/F for a +100 mV depolarization in HEK293
cells, respectively). The ArcLight A389 A390 V442 mutant (named Marina)
has single exponential kinetics of activation with a time constant
of 29.2 ± 2.0 ms and biexponential kinetics of deactivation with
time constants 15.6 ± 1.5 ms for the fast component and 59.4
± 5.8 ms for the slow component (fast component accounting for
61% of response; Figure S2).The
detected reversal in polarity of the voltage-dependent response
of ArcLight seems to be dependent on the presence of amino acids with
hydrophobic (alanine, tryptophan, cysteine, and methionine) side chains
at position 389. The hydrophobic side chains of alanine and tryptophan
produced the largest effects. Hydrophilic side chains (i.e., aspartic
acid, glutamic acid, lysine, and arginine) at this position do not
produce probes with a positive slope relationship to voltage. Replacement
of polar tyrosine (Y) on another outward oriented residue 442 within
ArcLight FP with a range of hydrophobic amino acids (leucine, isoleucine,
valine, alanine, and methionine) resulted in variants with increased
voltage sensitivity. Again, introduction of hydrophilic residues (aspartic
acid or glutamic acid, data not shown) to this position diminished
voltage sensitivity.Position 442 appears to be more promiscuous
to residue size since
multiple amino acids (valine, leucine, and isoleucine) produced probes
with similar voltage sensitivities. We also demonstrated that an equally
significant improvement in voltage sensitivity (more than double)
of the GEVI Marina was achieved through the replacement of the hydrophilic
histidine (H) on position 390 with hydrophobic residues, either alanine
(A) or cysteine (C). The introduction of other larger hydrophobic
residues at position 390 caused smaller increases in voltage sensitivity.
Replacement of H390 with other polar, hydrophilic groups caused not
only decreases in signal magnitude, but some substitutions (arginine
and asparagine) produced a reversal of the probe response polarity.
Our results suggest that both the size and charge of the side chain
of residue 390 play a significant role in voltage sensitivity. An
increase in surface hydrophobicity of the FP supports the hypothesis
that the voltage sensitivity of single FP GEVIs might depend on interactions
between the FP surface and the plasma membrane.[29,37] The effects of changing residue 390 on voltage sensitivity are not
surprising given the large body of evidence describing different ways
in which H390 (H148 in avGFP) affects the fluorescent properties of
GFP, including altering proton transfer, and protein conformation
and flexibility.[42−46]Comparisons of the V/F curves
of the novel GEVI variants (Figures A, B and S3) showed that
all three mutations (D389A, H390A, and Y442V) introduced into the
parent GEVI ArcLight contributed to the increased sensitivity of the
probe to depolarizing events, while all three probes exhibited nearly
identical lack of sensitivity to hyperpolarizing events. Spontaneously
spiking HEK293 cells showed decreases in baseline fluorescence for
the GEVI Marina (180 ± 5 AU, n = 3045 cells)
compared to the parent ArcLight (266 ± 5 AU, n = 8987 cells). Mutational changes in Marina seem to lower the photobleaching
constants of the indicator compared to the parent probe ArcLight measured
in HEK293 cells (Figure C). Photobleaching curves for both probes were fitted to single exponential
curves resulting in time constants of 206.3 ± 20.8 s for Marina
and 132.3 ± 4.8 s for ArcLight.
Figure 3
Characterization of novel GEVI Marina.
In experiments shown in
A and B, HEK293 cells transiently expressing various GEVIs were simultaneously
electrically recorded using whole-cell patch clamp (300 ms 100 mV
depolarizing steps from −70 mV holding potential) and imaged
with a high-speed (1000 Hz) CCD camera. In all experiments, excitation
light intensity measured at the sample plane was 18 mW/mm2. (A) Depolarization changes in membrane potential cause a decrease
in fluorescence intensity in ArcLight (in blue) and an increase in
fluorescence intensity in ArcLight D389A (light green), ArcLight A389
H390A (green) and ArcLight A389 A390 Y442V-Marina (dark green). All
traces are single trials and were subtracted with a double exponential
fit for bleach correction. See also Figures S2 and S3 for unsubtracted traces. (B) V–F curve for all the indicators shown in A (4–5 cells
recorded for each indicator). For all cells, voltage steps of −50
to +100 mV were applied in increments of 10 mV from a resting potential
of −70 mV. Data are plotted as mean ± SEM. (C) GEVI Marina
exhibits lower rate of photobleaching compared to the parent probe
ArcLight. Electrically active HEK293 cells transiently expressing
GEVIs ArcLight (blue), ArcLight A389 A390 (light green) and ArcLight
A389 A390 V442 (dark green) were exposed to 60 cycles of intermittent
laser excitation (duty cycle of 5 s on/55 s off) resulting in 300
s total laser excitation illumination (18 mW/mm2). (AU
= arbitrary units.)
Characterization of novel GEVI Marina.
In experiments shown in
A and B, HEK293 cells transiently expressing various GEVIs were simultaneously
electrically recorded using whole-cell patch clamp (300 ms 100 mV
depolarizing steps from −70 mV holding potential) and imaged
with a high-speed (1000 Hz) CCD camera. In all experiments, excitation
light intensity measured at the sample plane was 18 mW/mm2. (A) Depolarization changes in membrane potential cause a decrease
in fluorescence intensity in ArcLight (in blue) and an increase in
fluorescence intensity in ArcLight D389A (light green), ArcLight A389
H390A (green) and ArcLight A389 A390 Y442V-Marina (dark green). All
traces are single trials and were subtracted with a double exponential
fit for bleach correction. See also Figures S2 and S3 for unsubtracted traces. (B) V–F curve for all the indicators shown in A (4–5 cells
recorded for each indicator). For all cells, voltage steps of −50
to +100 mV were applied in increments of 10 mV from a resting potential
of −70 mV. Data are plotted as mean ± SEM. (C) GEVI Marina
exhibits lower rate of photobleaching compared to the parent probe
ArcLight. Electrically active HEK293 cells transiently expressing
GEVIs ArcLight (blue), ArcLight A389 A390 (light green) and ArcLight
A389 A390 V442 (dark green) were exposed to 60 cycles of intermittent
laser excitation (duty cycle of 5 s on/55 s off) resulting in 300
s total laser excitation illumination (18 mW/mm2). (AU
= arbitrary units.)To examine the potential
to reveal electrical events in neuronal
cells, GEVI Marina was transiently expressed in mouse primary neuronal
culture in vitro. Representative mouse cortical neurons in vitro transiently
expressing GEVI Marina are shown in Figure . Simultaneous optical and electrical recordings
revealed that voltage sensitivity of Marina allows for detection of
membrane transients (single action potentials and subthreshold events)
from neuronal cell bodies and processes (Figure ). Spontaneous action potentials in Marina
expressing neurons could be reliably seen optically (Figure ; 1.4 ± 0.1%ΔF/F for action potentials of 37.4 ±
3.7 mV; average of 35 APs in 5 different cells) with a signal-to-noise
ratio (SNR) of 4.5 ± 0.2. Under identical experimental conditions,
the parent probe ArcLight resolves action potentials with an average
fluorescence change of 2.6 ± 0.7%ΔF/F for action potentials of 48.9 ± 8 mV; average of
21 APs in 3 different cells) with a signal-to-noise ratio (SNR) of
8.8 ± 1.5.
Figure 4
Spontaneous electrical activity in cultured mammalian
neurons detected
with GEVI Marina. (Top) GEVI Marina exhibiting membrane localization
in transiently expressing mouse cortical neurons in vitro. Image is
a Z-stack of nine 2 μm thick sections taken with a two-photon
microscope equipped with 20× objective. Bar scale: 20 μm.
(Bottom) Sample traces of single trial recordings of spontaneous action
potential bursts from three locations of the neuron shown. Fluorescence
recordings at 125 Hz (colored): % ΔF/F from the area of interest in the same color on the left;
soma electrode voltage recording (white). All optical traces have
double exponential subtraction to correct for FP bleach. An 80 ×
80 image (right) of a neuron with the regions of interest averaged
to produce the traces shown. Excitation light intensity on the sample
plane was 18 mW/mm2.
Figure 5
Simultaneous electrical and optical recordings of spontaneous action
potential activity from cortical neuronal cells in vitro expressing
GEVI Marina. (A and B) Single-trial fluorescence traces of activity
in neuronal cell bodies of two different cells. (C) Magnification
of trace marked in (B) with red borderline. Optical trace in green,
and electrical in black. Excitation light intensity on the sample
plane was 18 mW/mm2. All fluorescence traces were recorded
at 500 Hz using high speed CCD camera. All traces are bleach corrected.
Traces in (A) and (B) are unfiltered, trace in (C) is filtered using
low pass Kaisser-Bessel 30 filter (150 Hz cutoff).
Spontaneous electrical activity in cultured mammalian
neurons detected
with GEVI Marina. (Top) GEVI Marina exhibiting membrane localization
in transiently expressing mouse cortical neurons in vitro. Image is
a Z-stack of nine 2 μm thick sections taken with a two-photon
microscope equipped with 20× objective. Bar scale: 20 μm.
(Bottom) Sample traces of single trial recordings of spontaneous action
potential bursts from three locations of the neuron shown. Fluorescence
recordings at 125 Hz (colored): % ΔF/F from the area of interest in the same color on the left;
soma electrode voltage recording (white). All optical traces have
double exponential subtraction to correct for FP bleach. An 80 ×
80 image (right) of a neuron with the regions of interest averaged
to produce the traces shown. Excitation light intensity on the sample
plane was 18 mW/mm2.Simultaneous electrical and optical recordings of spontaneous action
potential activity from cortical neuronal cells in vitro expressing
GEVI Marina. (A and B) Single-trial fluorescence traces of activity
in neuronal cell bodies of two different cells. (C) Magnification
of trace marked in (B) with red borderline. Optical trace in green,
and electrical in black. Excitation light intensity on the sample
plane was 18 mW/mm2. All fluorescence traces were recorded
at 500 Hz using high speed CCD camera. All traces are bleach corrected.
Traces in (A) and (B) are unfiltered, trace in (C) is filtered using
low pass Kaisser-Bessel 30 filter (150 Hz cutoff).Voltage dependent dimerization between fluorescent
proteins of
different indicator molecules has been proposed as a mechanism of
voltage sensitivity of ArcLight and ArcLight-derived indicators.[37,47] Residues D389, F444, and T446 (S147, S202, and Q204 in avGFP, Figure S1) impart the high ΔF/ΔV sensitivity of ArcLight.[37] These residues are considered to be part of a larger group
of amino acids that form the dimer interface in GFP.[48,49] However, only substitutions A206K, L221K, and F223R, individually
or in combination, have been shown to disrupt dimer formation, rendering
GFP monomeric.[50] To test directly the possible
influence of the dimer interface forming residues on voltage sensitivity,
we made saturation libraries at these three positions (T448, L463,
and F465 in Marina; Figure ). None of these mutations produced indicators with increased
voltage sensitivity compared to parent probe. In addition, mutagenic
libraries of several other residues present in the dimerization interface
of GFP[49] N149, Y151, R168, S202, and Q204
(Q391, Y393, R410, F444, and T446 in Marina) also did not produce
a more sensitive indicator. Similar to what was seen for ArcLight
residues, F444, T446, and D469 play crucial roles in the voltage sensitivity
of Marina with many mutants diminishing or even reversing the voltage-dependent
signal. The fact that mutations to multiple residues that are present
in the dimerization interface disrupt voltage sensitivity of the probe
is compatible with the hypothesis that intermolecular dimerization
plays role in ArcLight probe function.
Figure 6
Effects of mutagenic
changes in dimerization interface of FP on
voltage-dependent response. (A and B) For each of the saturated mutagenesis
libraries, 96 colonies were selected and tested for expression and
voltage sensitivity. Screening of each library was preformed four
times on a custom-made, semiautomated screening platform which allows
for simultaneous field stimulation and optical imaging in the 96-well
dish format. Screening was performed in electrically active HEK239
cells transiently expressing various mutants. Each dot on the graph
represents the highest value of depolarization-dependent fluorescence
change detected in each well. (A) ArcLight A389 A390 Y442V-Marina
was used as a template for making site-directed mutagenesis libraries
on residues which we found to be essential for ArcLight voltage sensitivity
and/or involved in FP dimerization (Q391, R410, F444, T446, and D469).
(B) ArcLight A389 A390 Y442L was used as a template for making site-directed
mutagenesis libraries on residues involved in FP dimerization (Y393,
T448, L463, and F465).
Effects of mutagenic
changes in dimerization interface of FP on
voltage-dependent response. (A and B) For each of the saturated mutagenesis
libraries, 96 colonies were selected and tested for expression and
voltage sensitivity. Screening of each library was preformed four
times on a custom-made, semiautomated screening platform which allows
for simultaneous field stimulation and optical imaging in the 96-well
dish format. Screening was performed in electrically active HEK239
cells transiently expressing various mutants. Each dot on the graph
represents the highest value of depolarization-dependent fluorescence
change detected in each well. (A) ArcLight A389 A390 Y442V-Marina
was used as a template for making site-directed mutagenesis libraries
on residues which we found to be essential for ArcLight voltage sensitivity
and/or involved in FP dimerization (Q391, R410, F444, T446, and D469).
(B) ArcLight A389 A390 Y442L was used as a template for making site-directed
mutagenesis libraries on residues involved in FP dimerization (Y393,
T448, L463, and F465).Due to the multiple states of its chromophore, wildtype GFP
(wtGFP)
shows two absorption maxima: the major band absorbs at 398 nm and
the minor band absorbs at 475 nm, each of which have similar emissions
peaks (508 and 503 nm, respectively).[51,52] Excitation
causes the chromophore to shift from a neutral (absorption at 398
nm) to an anionic form (absorption at 475 nm) via an excited-state
proton transfer (ESPT) resulting in green emission. The ESPT reaction
involves multiple residues surrounding the chromophore. It has been
shown that mutations of critical residues can result in ESPT disruption
but also GFP rewiring with proton pathways relying on alternative
proton acceptors.[44] In wtGFP, mutation
of S65T causes the chromophore to exist primarily in the deprotonated
form resulting in a single excitation peak at ∼488 nm.[53,54]The original ArcLight (epArcLight), based on fusion of the
voltage
sensitive domain (VSD) of the Ciona voltage sensitive
phosphatase (CiVSP) and ecliptic pHluorin ArcLight) contains a serine
at the 65 residue of the GFP (as present in the original ecliptic
pHluorin residue 307 from the starting methionine in Ciona VSD).[10,39] The ecliptic pHluorin FP imparts dual excitation peaks (395 and
475 nm) in epArcLight, which resembles the spectrum of wild-type GFP[18] and GCaMP. The ArcLight in popular use, however,
is one in which the FP includes the S65T mutation (i.e., super ecliptic
pHluorin)[55] found in eGFP, which causes
a simplification in the excitation spectrum to a single peak centered
at ∼488 nm.[53] This version of ArcLight
was adopted because its resting brightness was higher.[10]To test the voltage sensitivity of different
forms of the GFP chromophores
(neutral vs anionic), we measured voltage sensitivity in both “negative”
and “positive” going variants of ArcLight using either
400 or 488 nm excitation light (Figures and S4). We examined
the voltage/fluorescence response of epArcLight and found that a depolarization-dependent
increase in fluorescence was present when ∼400 nm light was
used, which is a reversal of the response seen with ∼488 nm
excitation (Figure A). In contrast to epArcLight, super ecliptic pHluorin ArcLight shows
no depolarization-dependent increase in fluorescence when excited
at ∼400 nm (Figure B). In addition, adjustments of the linker length, which increased
voltage sensitivity of super ecliptic pHluorin ArcLight,[10] had no effect on the voltage sensitivity of
the probe at ∼400 nm (Figure C).
Figure 7
Effects of excitation light wavelength on voltage sensitivity
in
ArcLight-based GEVIs. The voltage sensitivity was tested using whole
cell patch clamp electrophysiology combined with optical recording
in nonelectrically active HEK293 cells transiently expressing various
GEVIs (see also Figure S4). Cells were
held at −70 mV holding potential and stepped to +30 mV for
300 ms. All traces were recorded at 1000 Hz using high speed CCD camera.
Traces were recorded using either ∼400 nm (traces in light
green) or ∼488 nm (traces in dark green) excitation light from
LED light source. Light intensity on the sample plane for both wavelengths
was adjusted to 8.5 mW/mm2. All traces are single trials
and nonfiltered with bleaching removed by subtracting a double exponential
fit. (A) Depolarization changes in membrane potential cause either
an increase or decrease in fluorescence intensity in Ecliptic pHluorin
ArcLight (S249) depending on excitation light wavelength used, ∼400
and ∼488 nm, respectively. Super Ecliptic pHluorin GFP-based
GEVIs show no depolarization dependent fluorescence change when ∼400
nm excitation light is used. (B) Super Ecliptic pHluorin ArcLight
(S249) and (C) Super Ecliptic pHluorin ArcLight (Q239). (D) ArcLight
(A389 A390), (E) ArcLight (S307 A389 A390) and Marina in (F) exhibit
significant depolarization dependent voltage sensitivity only when
illuminated with ∼488 nm excitation light.
Effects of excitation light wavelength on voltage sensitivity
in
ArcLight-based GEVIs. The voltage sensitivity was tested using whole
cell patch clamp electrophysiology combined with optical recording
in nonelectrically active HEK293 cells transiently expressing various
GEVIs (see also Figure S4). Cells were
held at −70 mV holding potential and stepped to +30 mV for
300 ms. All traces were recorded at 1000 Hz using high speed CCD camera.
Traces were recorded using either ∼400 nm (traces in light
green) or ∼488 nm (traces in dark green) excitation light from
LED light source. Light intensity on the sample plane for both wavelengths
was adjusted to 8.5 mW/mm2. All traces are single trials
and nonfiltered with bleaching removed by subtracting a double exponential
fit. (A) Depolarization changes in membrane potential cause either
an increase or decrease in fluorescence intensity in Ecliptic pHluorin
ArcLight (S249) depending on excitation light wavelength used, ∼400
and ∼488 nm, respectively. Super Ecliptic pHluorin GFP-based
GEVIs show no depolarization dependent fluorescence change when ∼400
nm excitation light is used. (B) Super Ecliptic pHluorin ArcLight
(S249) and (C) Super Ecliptic pHluorin ArcLight (Q239). (D) ArcLight
(A389 A390), (E) ArcLight (S307 A389 A390) and Marina in (F) exhibit
significant depolarization dependent voltage sensitivity only when
illuminated with ∼488 nm excitation light.In all “brightening” variants of ArcLight (ArcLight
A389 A390 and ArcLight S307 A389 A390), regardless of the amino acid
residue present at position 65 within the FP (either serine or threonine),
a voltage dependent fluorescence intensity change is present only
when ∼488 nm excitation light is used (Figure D–F). Therefore, it appears that for
the majority of pHluorin-based probes, changes in membrane potential
affect the deprotonated and not the protonated form of the fluorophore.
Possibly, the voltage change results in the translocation of the FP
in relation to an external fixed neighboring charge (possibly the
charged head groups of the membrane). This could alter the ionic environment
and thereby the fluorescence output of the fluorophore. However, these
voltage dependent movements might not affect the neutral form of the
fluorophore. The one exception to this is the ecliptic pHluorin version
of ArcLight, in which the voltage dependent translocation possibly
results in a shift of the fluorophore between the two states.The present studies indicate that molecular rearrangements in a
voltage sensing domain that accompany a change in membrane voltage
can result in either an increase or decrease in output fluorescence
from the attached fluorescent protein. We have found that mutations
to D389 (147 position in eGFP) in ArcLight produce an indicator with
a weak positive slope ΔF/ΔV relationship. In addition, we found that mutations of two more amino
acid residues in the FP, inward oriented 390 and outward oriented
442, result in a GEVI (Marina) with a large “positive”
response. The introduced mutations did not significantly alter the
time constants of the response compared to the parent probe ArcLight.
The most prominent Marina variant (ArcLight A389 A390 V442) is capable
of resolving action potentials and subthreshold events in spontaneously
active mammalian neurons. The reversal in signal polarity is an essential
first step to production of GEVIs with signal-to-noise characteristics
comparable to GCaMPs that will enable more widespread in vivo use.
Methods
Molecular Biology and Construction
of Mutagenesis Libraries
All constructs used in this study
were derived from genetically
encoded-voltage indicator ArcLight-Q239 (GenBank ID: AB183035, AY533296; Plasmid
#36856, Addgene, Cambridge, MA). In order to improve expression in
mammalian cells, the sequence of ArcLight Q239 was codon optimized
using mammalian codon preferences. In order to improve membrane trafficking
of the probe, we replaced residues 71–91 within CiVSD with
an amino acid sequence (KSRITSEGEYIPLDQIDINV)
from the Kir2.1 channel Golgi-to-plasma membrane trafficking signal.[56] In order to prevent ER aggregation, the Kir
2.1 channel ER export signal sequence (FCYENEV) was added to
C-terminus of ArcLight.[57] ArcLight Q239
was fused with a self-cleaving T2A peptide sequence (GSGEGRGSLLTCGDVEENPGP)
followed by nuclear localized mCherry. The restriction sites were
inserted between the promoter, voltage sensor domain, fluorescent
protein domain and nuclear localized tag protein (ApaI, KpnI, EcoRI, SphI, and NheI, respectively) to aid cloning. In order
to study the effects of excitation light wavelength on voltage sensitivity,
additional probes were used: epArcLight (Ecliptic pHluorin fused to
CiVSD at the S249 position) and ArcLight S249 (Super Ecliptic pHluorin
fused to CiVSD at the S249 position, Plasmid #36855, Addgene). For
neuronal expression of Marina, GEVI CMV promotor was replaced with
hSynapsin.Targeted amino acid mutations were introduced by
using the QuikChange II XL site-directed mutagenesis kit (Agilent
Technologies, Inc., Santa Clara, CA). For each site-directed mutagenesis
library, 96 E. coli colonies were selected and cDNA
purified. Each library’s complexity was confirmed by sequencing
10% of selected colonies. We found that a selection of 96 constructs
will contain all 20 possible amino acids at the specific residue,
and therefore does not require costly and time-consuming sequencing
to identify each mutation. Alternatively, in some cases site directed
mutagenesis libraries were produced using InFusion system (Clontech
Laboratories, Inc., Mountain View, CA). We found that with this system
a selection of 48 constructs will contain all 20 possible amino acids
at the specific residue. Following functional screening on a semiautomated
platform, the most promising constructs from each round were sequenced
and the best were used as templates in future rounds of library production.
Cell Culture
This study was carried out in strict accordance
with the recommendations in the Guide for the Care and Use of Laboratory
Animals of the National Institutes of Health. The protocol was approved
by the Pierce Animal Care and Use Committee. HEK 293 cells (AATC,
Manassas, VA) were kept in Dulbecco’s modified Eagle’s
medium (DMEM, high glucose; Invitrogen, Carlsbad, CA) supplemented
with 10% fetal bovine serum (FBS; Sigma-Aldrich, Saint Louis, MO).
For functional testing on the semiautomated screening platform, probes
were expressed in spontaneously spiking HEK cell line (kind gift from
Dr. Adam Cohen of Harvard University).[41] Cells were cultured in DMEM/F12, 10% FBS, 1% penicillin (100 U/mL),
streptomycin (100 μg/mL), Geneticin (500 μg/mL), and puromycin
(2 μg/mL) (Sigma-Aldrich). Cortical neurons were isolated from
E16 mouse embryos and maintained in Neurobasal medium with 0.5 mM
Glutamax-I and 1 mL of B-27 supplement (Invitrogen) per 50 mL of culture
medium. In all cases, cells were plated either on glass bottom 96-well
black dishes (Cellvis, Mountain View, CA) or on 12 mm coverslips (Carolina
Biological, Burlington, NC) coated with poly-d-lysine hydrobromide
(Sigma-Aldrich). Cultures were maintained in a humidified incubator
at 37 °C in a 5% CO2 environment. Transient transfection of constructs
was accomplished by using half of the manufacturer’s recommended
amount of DNA (0.1 μg per 96-well or 0.4 μg per 12 mm
coverslip in 24-well dish) and Lipofectamine2000 (0.25 μL per
96-well or 1 μL per 24-well; Invitrogen).
Functional
Testing of Potential Probes
A custom built
semiautomated screening platform was built around a Nikon Eclipse
Ti-E inverted microscope equipped with a Perfect Focus System and
fitted with a Nikon Plan Apo 20x 0.75 NA objective (Nikon, Tokyo,
Japan) and a pE-300white CoolLED LED illumination System (CoolLED,
Ltd., Andover, U.K.) light source. A custom built imaging chamber
was attached to a motorized Prior Proscan II stage (Prior Scientific,
Inc., Rockland, MA) and maintained constant temperature (37 °C)
and humidity of plates during experiments. All image trains were collected
using a Hamamatsu ORCA Flash 4.0 sCMOS camera (Hamamatsu, Japan) at
a frame rate of either 50 or 100 Hz. For ArcLight-based constructs
(GFP) we used a 472/30 nm excitation filter, 495 nm dichroic mirror
and 520/35 nm emission filter (Semrock, Rochester, NY). For nuclear
localized mCherry visualization we used a 560/40 nm excitation filter,
585 nm dichroic mirror and 630/75 nm emission filter (Chroma Technologies
Corp., Bellows Falls, VT). For field stimulation of each single well,
we used a custom-made field electrode and actuator system attached
to the roof of the imaging chamber. Electrical stimulation, imaging
workflow and signal detection was controlled with a custom application
written in LabView (National Instruments, Inc., Austin, TX). For cell
detection, 16 images of the well bottom were collected and analyzed
by a cell identification algorithm that quantified the fluorescence
intensity of nuclear localized mCherry. The four fields with the largest
number of cells were then selected by the program and used for analysis.
For signal detection, a series of images were collected at 100 fps
from the preselected fields of views using a GFP filter set. During
collection of each image series, a fixed pattern of electrical stimulation
was applied to the well. This routine was repeated until time-series
images have been recorded from all 96 wells. The region-of-interest
(ROI) mask for each cell within a field of view is created by a thresholding
algorithm applied to the mCherry channel. A “cell mask”
was generated by enlarging each of the mCherry nuclear features. A
mean intensity trace for each cell is created from the time-series
images recorded using the GFP filter set. The magnitude of the fluorescence
signal in each cell was quantified (ΔF/F), and various values of construct performance were evaluated
to rank the signal size of a particular construct. Individual cells
expressing the same construct showed large variation in fluorescence
response amplitudes. We found it necessary to screen each construct
in at least four separate plates to get an accurate picture of its
performance. We also found that the maximum response amplitude across
cells and wells was the most reliable indicator of particular constructs
performance in subsequent patch clamp experiments. The values for
resting light intensity of HEK cells expressing either GEVI ArcLight
or Marina were derived from mean intensity traces recorded using the
GFP filter set at the time points prior voltage stimulation. For each
cell, resting light intensity value is average of three frames.
Whole Cell Patch Clamp Electrophysiology
Microelectrode
recordings were performed in a perfused chamber with the bath temperature
kept at 33–35 °C by a temperature feedback controller
(Warner Instruments, Hamden, CT). Bath solution contained: 129 mM
NaCl, 4 mM KCl, 1 mM CaCl2, 1 mM MgCl2, 10 mM d-glucose,
and 10 mM HEPES, pH 7.4 and was adjusted to 310 mOsm with sucrose.
We used 3–5 MΩ glass patch pipet (capillary tubing with
1.5/0.75 mm OD/ID-World Precision Instruments, Sarasota, FL) that
were pulled on a P-97 Flaming/Brown type micropipet puller (Sutter
Instrument Company, Novato, C). For HEK293 cell whole cell patch recordings,
the pipet solution contained 125 mM K-gluconate, 8 mM NaCl, 0.6 mM
MgCl2, 0.1 mM CaCl2, 1 mM EGTA, 4 mM Mg-ATP, 0.4Na-GTP and 10 mM HEPES,
pH 7.2 and adjusted to ∼290 mOsm. Voltage-clamp recordings
in the whole-cell configuration were performed using MultiClamp 700B
amplifier (Molecular Devices, LLC, Sunnyvale, CA) with a holding potential
of −70 mV. Spontaneous activity of cultured cortical neurons
was recorded in current clamp mode without holding current injection.
The pipet solution for neuronal recordings contained: 130 mM K-gluconate,
10 mM KCl, 4 mM Mg-ATP, 0.3 mM Na2-GTP, 10 mM HEPES and 10 mM Na-phosphocreatin,
pH 7.3 and adjusted to ∼310 mOsm. Imaging was performed on
an Olympus BX61WI upright microscope using either a LUMFL 60×
N.A. 1.10 or a LUMPlan FL 40× N.A. 0.80 water immersion objective
(Olympus, Waltham, MA). In the experiments with laser light, all ArcLight-derived
constructs were excited with a 488 nm 50 mW laser (DL488-050, CrystaLaser,
Reno, NV), using a 495 nm dichroic mirror and 520/35 nm emission filter
(Semrock). Excitation light measured at the preparation was 18 mW/mm2, and was adjusted for each recording using a continuous circular
neutral density filter (ThorLabs, Inc., Newton, NJ) to the minimum
required to record optical signals. For experiments in which constructs
were excited with either 400 or 488 nm light, we used a SOLA Light
Engine LED light source (Lumencor, Inc., Beaverton, OR). Excitation
light measured at the preparation was 8.5 mW/mm2. For excitation
at 400 nm we used an ET405/40 nm excitation filter, T495lprx nm dichroic
mirror, and 525/50 nm emission filter (Chroma Technologies Corp.).
For experiments with 488 nm light excitation, we used a 469/35 nm
excitation filter, dichroic mirror 497 nm, and emission filter 525/39
nm (Semrock). The objective image was projected onto the 80 ×
80 pixel chip of a NeuroCCD-SM camera controlled by NeuroPlex software
(RedShirt Imaging, Decatur, GA). For image demagnification, we used
either an Optem zoom system A45731 0.13 or Optem C-to-C mount 25-70-54
0.383 (Qioptiq LINOS, Waltham, MA). Images of all cell types were
recorded at a frame rate of 125, 500, or 1000 Hz, and depicted optical
traces were simple spatial averages of intensity of all pixels within
the region of interest (ROI), with signals processed as previously
reported[10,58] with double-exponential fitting/subtractions
used to compensate for rapid and slow photobleaching. All traces are
single trials. SNRs were calculated as the maximum fluorescence response
to action potential (as identified by the patch-clamp recording) divided
by the SD of the baseline fluorescence. Initially, we used NeuroPlex
software (RedShirtImaging, GA) to view the image sequences and output
optical and electrophysiological recordings. The data were further
processed and statistically analyzed using Igor Pro 6 (Wavemetrics,
OR) and Numbers (Apple, Inc., CA).
Statistical Analysis
Results are expressed as means
± SEM, as indicated in figure legends. The statistical significance
was evaluated using Student’s t test using
IgorPro software (Wavemetrics, OR).
Authors: Deborah Stoner-Ma; Andrew A Jaye; Kate L Ronayne; Jérôme Nappa; Stephen R Meech; Peter J Tonge Journal: J Am Chem Soc Date: 2008-01-08 Impact factor: 15.419
Authors: Helen H Yang; François St-Pierre; Xulu Sun; Xiaozhe Ding; Michael Z Lin; Thomas R Clandinin Journal: Cell Date: 2016-06-02 Impact factor: 41.582
Authors: Divya Sitaraman; Yoshinori Aso; Xin Jin; Nan Chen; Mario Felix; Gerald M Rubin; Michael N Nitabach Journal: Curr Biol Date: 2015-10-08 Impact factor: 10.834
Authors: François St-Pierre; Jesse D Marshall; Ying Yang; Yiyang Gong; Mark J Schnitzer; Michael Z Lin Journal: Nat Neurosci Date: 2014-04-22 Impact factor: 24.884
Authors: Jordan S Leyton-Mange; Robert W Mills; Vincenzo S Macri; Min Young Jang; Faraz N Butte; Patrick T Ellinor; David J Milan Journal: Stem Cell Reports Date: 2014-02-06 Impact factor: 7.765
Authors: Yao-Hui Sun; Hillary K J Kao; Che-Wei Chang; Alexander Merleev; James L Overton; Dalyir Pretto; Sergey Yechikov; Emanual Maverakis; Nipavan Chiamvimonvat; James W Chan; Deborah K Lieu Journal: Stem Cells Date: 2019-09-30 Impact factor: 6.277
Authors: Masoud Sepehri Rad; Yunsook Choi; Lawrence B Cohen; Bradley J Baker; Sheng Zhong; Douglas A Storace; Oliver R Braubach Journal: Biophys J Date: 2017-11-01 Impact factor: 4.033
Authors: Vincent Villette; Mariya Chavarha; Ivan K Dimov; Jonathan Bradley; Lagnajeet Pradhan; Benjamin Mathieu; Stephen W Evans; Simon Chamberland; Dongqing Shi; Renzhi Yang; Benjamin B Kim; Annick Ayon; Abdelali Jalil; François St-Pierre; Mark J Schnitzer; Guoqiang Bi; Katalin Toth; Jun Ding; Stéphane Dieudonné; Michael Z Lin Journal: Cell Date: 2019-12-12 Impact factor: 41.582