Introduction

In vivo fluorescence optical imaging is a powerful method to characterize biological processes in live animal tissues. Yet, the spatial precision at which these processes can be studied is constrained by the fact that lateral resolutions in conventional optical imaging are diffraction limited to ~200 nm, a length scale well above the nanometre-range interactions of most biomolecules. In recent years, single-molecule (SM) fluorescence microscopy techniques have provided means to study the mobility and the function of biomolecules beyond this resolution limit and with a precision of a few nanometres in cultured cells1. Extending these techniques to live animal imaging could open new avenues to examine the nanoscale behaviours of signaling molecules under homeostatic control within live tissues and during various developmental or pathogenic stages. Despite previous efforts to track individual proteins and study protein–protein interactions in embryos and early larval stages of zebrafish2,3,4 and C. elegans nematodes5,6, SM fluorescence imaging in intact adult animals remains highly challenging. In vivo, tissue autofluorescence often results in poor signal-to-noise ratio when attempting to detect individual fluorophores. This issue can be addressed using high-contrast imaging techniques such as highly inclined and laminated optical sheet excitation (HILO)7 or selective plane illumination microscopy8,9; however, advanced labelling methods are still required for site-specific SM imaging. In particular, maintaining a nanomolar SM detection regimen at subcellular locations in a live animal necessitates site-directed and stochiometric labelling or detection of a few biomolecules, despite their endogenous expression at levels generally well above those required for SM studies.

SM detection regimes can be achieved using traditional fluorescent protein (FP) fusions expressed in vivo; however, it generally requires extensive and potentially toxic photobleaching or alteration of protein expressions6. While this can be avoided using fusions to photoactivable FPs10, the photoactivation of a limited number of molecules at specific organelles within live animal tissues is complicated by the fact that FP mature before the fusions reach their final cellular destination. An alternative is to label biomolecules post-translationally and in a controllable manner by intravital injection of synthetic fluorescence tags. However, most labelling strategies demand micromolar concentrations of probes11, which generally results in high levels of nonspecific binding to tissues that cannot be washed-off. Moreover, residual non-targeted probes circulating in the body of injected animals further hinder SM imaging.

Recently, we introduced complementation-activated light microscopy (CALM), a SM technique that allows nanometre-precision imaging of individual biomolecules in living cells, independently of their expression levels, at micromolar concentration of fluorescent probes, and without the need to wash excess probes12. In CALM, proteins of interest are fused to a dark, non-mature split-GFP (green fluorescent protein) that lacks a 16 amino-acid β-strand. They are then fluorescently activated by irreversible binding and complementation with synthetic peptide sequences encoding the missing GFP β-strand (M3 peptides). These small synthetic M3 peptides permit site-specific imaging and tracking of individual fusion proteins by direct detection of activated split-GFP (dCALM) or by single pair Förster resonance energy transfer with fluorescently labelled M3 peptides (CALM-spFRET)12. Here we applied CALM to tissue-specific imaging and SM tracking of CD4 and voltage-dependent Ca2+ channel (VDCC) membrane proteins in live adult C. elegans.

In C. elegans, VDCCs are located in the membrane of excitable cells where they mediate the influx of extracellular Ca2+ in response to depolarizing signals. They play an essential role for electrochemical coupling in motor neurons by triggering neurotransmitter release13 and for excitation–contraction coupling in body-wall muscle cells where they elicit action potentials14. A dysregulation of VDCC activity in dystrophin-mutant C. elegans has recently been proposed to be a key factor for the genesis of muscular dystrophy15; however, the channels’ location and membrane dynamics have not been studied in living animals. Here we imaged and tracked individual VDCC at the sarcolemma of muscle cells and within neuromuscular synapses of normal and dystrophin-mutant C. elegans, using CALM. We demonstrate that dystrophin is a load-bearing element and an adaptive tension transducer involved in modulating the nanoscale confinement of VDCC at the sarcolemma in response to changes in muscle tonus. Our results illustrate that SM imaging and tracking by CALM allows a phenotyping of diseases such as muscular dystrophy with subresolution precision directly in intact animals.

Results

Live tissue fluorescence imaging by CALM in C. elegans

To explore the applications of dCALM and CALM-spFRET for live tissue imaging, we first used a C. elegans strain expressing a transmembrane CD4-split-GFP fusion under the muscle-specific myo-3 promoter16 (Fig. 1a). Split-GFP complementary M3 peptides were microinjected at high concentration in the fluid-filled pseudocoelomic cavity of live worms to allow their distribution to all organs and the activation of CD4-split-GFP on muscle cells (Fig. 1b and Supplementary Movie 1). Injected M3 peptides rapidly triggered the specific GFP activation of CD4-split-GFP on body-wall and vulval muscles (Fig. 1c and Supplementary Fig. 1). When worms were injected with M3 peptides conjugated to an Alexa 647 fluorophore (M3-A647)-specific CALM-spFRET imaging of CD4-split-GFP could also be performed at 680 nm using only 488 nm excitation (Fig. 1d). The large-excitation Stokes shift afforded by CALM-spFRET prevented a direct excitation of residual M3-A647 nonspecifically bound to other tissues and allowed muscle-specific near-infrared imaging of CD4-split-GFP in living worms (Fig. 1d). FRET signals on muscle cells were exclusively triggered by 1:1 stochiometric and irreversible binding of M3-A647 to CD4-split-GFP, as confirmed by in vitro fluorescence correlation spectroscopy (FCS) with recombinant split-GFP at picomolar concentrations (Supplementary Fig. 2). Neither complementation nor FRET was observed in transgenic worms injected with unconjugated A647 (Fig. 1d).

Figure 1: Ensemble CALM imaging of CD4-split-GFP in C. elegans.
figure 1

(a) GFP immunofluorescence staining of C. elegans expressing CD4-split-GFP under a muscle-specific promoter. CD4-split-GFP is detected in body-wall muscle somas, muscle arms (arrowheads, inset) and in vulval muscles (arrow). Scale bars: 150 and 25 μm (inset). (b) Schematic of dCALM and CALM-spFRET by microinjection of synthetic M3 peptides in the pseudoceolomic cavity of worms and activation of transmembrane CD4-split-GFP at muscle cell membranes. Adapted with permission from ref. 62 and from Z. F. Altun and D. H. Hall (wormatlas.org). (c) Live worm dCALM confocal imaging of activated CD4-split-GFP at the surface of body-wall muscles. Scale bar: 100 μm. (d) CALM-spFRET confocal imaging of CD4-split-GFP in live C. elegans. Single excitation at 488 nm allows muscle-specific near-infrared imaging of CD4-split-GFP (top) by eliminating interferences from nonspecifically bound peptides detected by dual 488 and 647 nm excitation (middle, arrowheads). No complementation is observed when injecting worms with the unconjugated A647 dye (bottom). Scale bars: 25 μm.

Single-molecule tracking of CD4-split-GFP on muscle cells

In CALM, the simple titration of complementary M3 peptides in extracellular fluids permits a controlled SM imaging of split-GFP fusion proteins specifically at the outer plasma membrane of cells because the peptides do not activate intracellular split-GFP if they are not injected directly in the cytoplasm12. To control the activation rate of membrane CD4-split-GFP and achieve continuous SM detection, we injected worms with low concentrations of M3 peptides. SM dCALM imaging was performed at the interface between the basal lamina and the sarcolemma of body-wall muscles using 488 nm HILO excitation. Under these conditions, activated CD4-split-GFP appearing as individual GFP diffraction-limited spots (Fig. 2a,b and Supplementary Movie 2) could be tracked by two-dimensional (2D) Gaussian fitting of their point-spread function with a mean localization accuracy of 32±8 nm (Fig. 2a,e and Supplementary Fig. 3). Worms were additionally imaged by SM CALM-spFRET after microinjection of M3-A647 peptides. At the sarcolemma, individual complemented CD4-split-GFPs were identified as diffraction-limited and colocalized diffusing spots that emitted in both GFP and A647 detection channels when excited at 488 nm (Fig. 2c). Single bleaching steps and blinking events where both GFP and A647 fluorescence signals disappeared (Fig. 2c) indicated an effective intramolecular spFRET between CD4-split-GFP and M3-A647 and confirmed our ability to specifically target, activate, detect and track individual transmembrane CD4-split-GFP directly in C. elegans.

Figure 2: Single-molecule tracking of CD4-split-GFP on muscle cells by dCALM and CALM-spFRET.
figure 2

(a) Individual activated CD4-split-GFP trajectories at the plasma membrane of a dorsal muscle cell. Scale bar: 20 μm. (b) Single frame image from the white square in a showing a 2D and 3D representation of the GFP signal for two individual activated CD4-split-GFP (arrowheads). Scale bar: 2 μm. (c) CALM-spFRET imaging of an individual CD4-split-GFP activated by an M3-A647 peptide (yellow arrow head). Individual time frames show the overlay of GFP and A647 channels during 488 nm excitation only. Blinking (1.2 s) and single-step photobleaching (1.4 s) are observed. Diffraction-limited spots are intentionally expanded to facilitate visualization. Scale bar: 1 μm. (d) Analysis of CD4-split-GFP diffusion by PDSD for time lag Δt=100 ms. The distribution is best described by a three-parameter fit indicative of three regimes of diffusion for CD4-split-GFP at the sarcolemma. (e) Representative example of individual CD4-split-GFP trajectories. (f) Diffusion coefficient analysis of the three diffusion regimes for 6193 CD4-split-GFP on 21 muscle cells from five worms. The high-mobility regime is best fitted with a free diffusion model at short times (red fit), while both low-mobility regimes show confined diffusion (green and blue fits). Inset: Close up on low-mobility regime 3. Error bars represent the s.e. of each mean.

When analysing the diffusion of activated CD4-split-GFP (6,193 trajectories, 21 muscle cells, five worms) by probability distribution of the squared displacements (PDSD)17, three diffusion regimes were identified (Fig. 2d). A first high-mobility regime (10%) with a diffusion coefficient DCD4−1=1.14 × 10−1±5 × 10−3μm2s−1 showed free diffusion for a length scale of up to (Fig. 2f). Two additional low-mobility regimes with diffusion coefficients DCD4−2=2.34 × 10−2±3.4 × 10−3μm2s−1 (38%) and DCD4−3=3.31 × 10−3±8.0 × 10−4μm2s−1 (52%) were each confined to membrane microdomains 420±20 and 110±4 nm in diameter, respectively (Fig. 2f). The diffusion coefficients and confinement sizes of these three regimes are consistent with the generally heterogeneous and restricted diffusion of other transmembrane proteins in cultured muscle cells18,19; however, the largely confined and slow diffusive behaviours of CD4-split-GFP in vivo likely stem from extracellular matrix effects that are absent in cultured cells and from tracking at lower temperatures. Importantly, these results demonstrate that the nanoscale diffusion of individual membrane proteins and their interactions with plasma membrane microdomains can be quantitatively studied at the SM level in intact animals using CALM.

CALM imaging of plasma membrane L-type VDCC in C. elegans

To image VDCC by CALM, we fused split-GFP to the N terminus of the α2δ VDCC subunit (UNC-36) that associates with both neuronal (UNC-2)20 and muscular (EGL-19)21 α1 subunits to form active membrane Ca2+ channels (Fig. 3a). unc-36 null allele (e251) worms were rescued by reintroducing UNC-36-split-GFP under the control of the endogenous unc-36 promoter via single copy insertion using MosSCI protocols22. Mutant worms were also rescued with UNC-36-GFP to compare CALM imaging with conventional GFP imaging of UNC-36. The resulting integrated worms had fully rescued VDCC activities. The phenotypic impairment of locomotion for UNC-36 mutation was completely restored in thrashing assays and the conducting activities of VDCC were identical to that of wild-type channels (Fig. 3b,c and Supplementary Fig. 4).

Figure 3: Imaging of plasma membrane L-type VDCCs in live C. elegans by genetic engineering of UNC-36 fusion proteins.
figure 3

(a) Schematic of split-GFP or full GFP fusion to UNC-36 (α2δ subunit) and its association with the α1 and β subunits to form active transmembrane VDCC. (b) C. elegans thrashing assays showing impaired motility in unc-36(e251) null allele mutant worms and functional rescue upon UNC-36-GFP and UNC-36-split-GFP gene insertion. The mean number of body bending per min (±s.d.) was determined for n=20 worms for all conditions (Welch‘s t-test, ***P<0.001). (c) Representative traces of inward VDCC currents recorded from body muscle cells in response to depolarizing pulses from −60 mV to −20, 0, +20 and +40 mV. Compared with wild-type N2 worms (n=8), the increase in activation midpoint potential (E0.5±s.e.m.) and time-to-peak kinetics of inward currents (±s.e.m.) observed for unc-36 (e251) worms (n=7) are restored to wild-type levels in UNC-36-GFP (n=7) and UNC-36-split-GFP worms (n=8; *E0.5: ANOVA test P=0.0141, Dunnett’s post test P<0.05 for unc-36(e251), P>0.05 for others; *time-to-peak: Kruskal–Wallis test P=0.0026, Dunn’s post test P<0.05 for unc-36(e251), P>0.05 for others). (d) Ensemble confocal fluorescence imaging of UNC-36-GFP and alexa-594 anti-GFP antibody in live C. elegans. UNC-36-GFP is expressed in muscle cells and motor neurons at the nerve cords. Anti-GFP antibodies target only the edge of muscle cells and the nerve cords. Scale bars: 5 μm. (e) Ensemble dCALM confocal imaging of UNC-36-split-GFP in live C. elegans expressing the acetylcholine receptor subunit UNC-29-tagRFP. Membrane UNC-36-split-GFP is detected at the sarcolemma of muscle cell bodies (left) and within neuromuscular synapses (red) along the nerve cords (arrowheads, right). Scale bars: 5 μm. (f) Immunofluorescence detection of GFP and of the beta-1 integrin PAT-3 showing UNC-36-split-GFP localized in the A-band region of sarcomeres between rows of dense bodies (DBs). Scale bar: 5 μm.

In live worms, UNC-36-GFP localized in muscle cells and motor neurons. Muscular UNC-36-GFPs were distributed at the boundary of body-wall muscles and as a punctuated pattern at the sarcolemma, while neuronal UNC-36-GFPs were distributed at the nerve cords (Fig. 3d). To assess the effective plasma membrane translocation of UNC-36 in these tissues, live worms were microinjected with fluorescently labelled anti-GFP antibodies. Anti-GFP staining was observed at the boundary of body-wall muscle cells and as punctuates along the nerve cords but not at the sarcolemma (Fig. 3d). The lack of sarcolemma staining was likely because of the inability of large antibodies to access crowded extracellular spaces at the interface between the epidermal and body-wall muscle cells. Indeed, when we injected much smaller M3 peptides and performed dCALM imaging, GFP fluorescence from membrane-activated UNC-36-split-GFP was specifically detected at the nerve cords and the muscle sarcolemma (Fig. 3e). Activated GFP signals were weaker than those of UNC-36-GFP suggesting a low, yet systematic, plasma membrane translocation of muscular and neuronal UNC-36. Activated UNC-36-split-GFP was localized at the membrane of neuromuscular synapses based on colocalization with the fluorescent acetylcholine receptor subunit UNC-29-tagRFP (Fig. 3e). This position is consistent with the expected location of VDCC in presynaptic active zones of motor neurons20. UNC-36-split-GFP was also detected at the sarcolemma, consistently with its interaction with the muscular EGL-19 α1 subunit of VDCC21. It localized in the A-band region, between rows of dense bodies, the structures that anchor sarcomeric actin filaments to the sarcolemma (Fig. 3f).

These results demonstrate that CALM affords specific in vivo targeting and imaging of membrane Ca2+ channels within confined extracellular spaces at the basal lamina–muscle interface and at neuromuscular synapses in C. elegans tissues.

Tracking of individual VDCC using CALM in C. elegans

We then used SM dCALM and HILO excitation to track individual VDCC in live worms. At the sarcolemma, activated VDCC appeared as individual diffraction-limited fluorescent spots during imaging (Supplementary Movie 3) and were distributed as punctuates on muscle cells (Fig. 4a), in a pattern similar to that observed using confocal microscopy. For each muscle cell, hundreds of individual VDCC could be tracked with a precision better than 30 nm (Supplementary Fig. 3). Muscle-wide trajectory maps revealed a largely confined diffusion of the channels at the sarcolemma (Fig. 4a). Individual VDCC tracked along the nerve cords were also seen diffusing within neuromuscular synapses when we overlapped UNC-36-split-GFP trajectory maps with UNC-29-tagRFP images for worms co-expressing both fusion proteins (Fig. 4b). The single emitter nature of both muscular and neuronal VDCC was confirmed using fluorescence intensity time trace analyses that displayed sudden activation, single-step photobleaching and blinking events typical of SM detection by dCALM12 (Fig. 4c). No activation of UNC-36-split-GFP was observed when imaging non-injected integrated worms (Supplementary Movie 4).

Figure 4: Single-molecule dCALM imaging and tracking of individual muscular and neuronal VDCCs in C. elegans.
figure 4

(a) Sum intensity image of a single-molecule dCALM movie (left) and corresponding trajectory map of individual VDCC (right) on a body-wall muscle cell. Trajectories with at least three steps are shown. Scale bar: 10 μm. Inset: representative diffusion trajectories of single VDCC from the sum and tracking images (white squares). Scale bar: 100 nm. (b) Overlay image of individual VDCC trajectories with UNC-29-tagRFP-labelled motor neuron synapses (red). Individual channels diffuse at the nerve cord (arrowheads) within synapses (inset). Scale bars: 2 μm and 200 nm (inset). (c) Background-corrected fluorescence intensity time traces of individual muscular or neuronal VDCC.

Muscle contraction modulates VDCC diffusion in nanodomains

Taking advantage of the nanometre-precision imaging capabilities of dCALM, we compared the nanoscale distribution and diffusion of individual VDCC in resting body-wall muscles of worms immobilized with polystyrene beads and in muscles undergoing sustained contraction for worms anaesthetized with the cholinergic agonist levamisole.

In resting muscles, two populations of slow and fast diffusing VDCC undergoing confined or short-scale free diffusion were identified at the sarcolemma (10,481 trajectories, 39 muscle cells, six worms, Fig. 5a,b). Ensemble PDSD analyses revealed that most channels (81%) diffused slowly (Dvdcc−slow=1.22 × 10−3±1.2 × 10−4μm2s−1) in membrane nanodomains 82±2 nm in diameter (Fig. 5c,e and Supplementary Table 1). The second VDCC population (19%) diffused 20-fold faster (Dvdcc−fast=2.22 × 10−2±1 × 10−3μm2s−1) and underwent free diffusion for a length scale of 230±8 nm (Fig. 5c,e and Supplementary Table 1), after which membrane obstacles induced a subdiffusive behaviour. Both populations were independently distributed at the sarcolemma (Supplementary Fig. 5), did not change diffusion and diffused more slowly than CD4-split-GFP, indicating that we imaged two behaviours of stably assembled VDCC rather than a dynamic association of the glycosylphosphatidylinositol-anchored UNC-36 (ref. 23) with the channels’ α1 subunit.

Figure 5: Diffusion analysis of individual VDCC at the sarcolemma of normal C. elegans worms expressing UNC-36-split-GFP.
figure 5

(a) Distribution of diffusion coefficients for 10,481 VDCC trajectories (39 muscle cells, six worms) determined from individual MSD analysis and showing that two populations of slow and fast diffusing channels co-exist at the sarcolemma of resting muscles. (b) Example of individual VDCC trajectories and the evolution of their MSD over time. VDCC undergo either confined or short-scale free diffusion as determined using 2D diffusion models that best fit the MSD (red). Scale bars: 200 nm. (c) Ensemble probability distribution of the squared displacement (PDSD) analysis of VDCC diffusive behaviours (10,481 VDCC, 39 muscle cells, six worms) in resting muscles. (d) Ensemble PDSD analysis of VDCC diffusive behaviours (36,325 VDCC, 62 muscle cells, 18 worms) in muscles under sustained contraction with levamisole. (e) Diffusion coefficients (±s.d.) for both fast and slow VDCC populations in resting and contracted muscles (F-test, *P<0.05, ***P<0.001). (f) VDCC nanodomain nearest-neighbour distances at the sarcolemma of resting and contracted muscles. The central squares and bars represent the mean of the distribution and its median, respectively. The box length represents the interquartile range and the error bars the s.d. of the mean (Wilcoxon sign-rank test, **P<0.01 compared with resting muscles).

Additional Ripley’s K function and nearest-neighbour distance analyses of the clustering and dispersion pattern of VDCC nanodomains at the sarcolemma indicated that they were not clustered but significantly more dispersed than expected for complete spatial randomness, with a mean nearest-neighbour distance of 769±339 nm (Fig. 5f and Supplementary Fig. 6). This distance, below the 1-μm average width of sarcomeres in adult C. elegans24 and above the width of one-half A-band (~350 nm), suggests that individual sarcomere contains two VDCC nanodomains. Because of the stochastic nature of CALM, we could not accurately determine the exact number of VDCC per nanodomain. However, when we injected worms with high concentrations of M3 peptides to activate as many VDCCs as possible at the sarcolemma, a count of single photobleaching steps in individual nanodomains and additional quantitative analyses of their fluorescence intensity distribution indicated that some nanodomains comprise at least 10 VDCC (Supplementary Fig. 7).

In muscles undergoing sustained contraction induced by levamisol (36,325 trajectories, 62 muscle cells, 18 worms), slow diffusing VDCC (91%) were nearly immobile with a 13-fold decrease in lateral diffusion coefficient compared with resting muscles (DLeva vdcc−slow=9.4 × 10−5±4 × 10−6μm2s−1, P<0.001, Fig. 5d,e). No confinement was detected, indicative of a strong reduction of the nanodomains’ size or the near-immobilization of VDCC within intact nanodomains. The diffusion of the fast VDCC population (9%) was moderately reduced by 33% compared with resting muscles (DLeva vdcc−fast=1.49 × 10−2±7 × 10−4μm2s−1, P<0.05); however, the channels’ free diffusion distance was unchanged at 224±9 nm (Fig. 5d,e and Supplementary Table 1). Interestingly, the nearest-neighbour distance between nanodomains significantly decreased to 546±217 nm compared with resting muscles (P<0.01, Fig. 5f and Supplementary Fig. 6). Transient reductions in nanodomain interdistance were also observed during sporadic muscle contractions of immobilized worms (Supplementary Movie 5) suggesting that, under sustained or active contractions, the sarcolemma undergoes a deformation that brings stable VDCC nanodomains closer to each other in a process reminiscent of sarcolemma festooning: a ballooning of the membrane between dense bodies when muscles are shortened to less than their resting length25.

Thus, the nanoscale mobility and the distribution of VDCC are modulated by the tonic state of muscles. Our observations that a large majority of VDCC are constrained in two nanodomains on each sarcomere indicate that VDCC action potentials and extracellular Ca2+ influx are spatially controlled in C. elegans body-wall muscles and that they are further regulated, on the nanometre scale, by a near-immobilization of the channels during sustained cholinergic contractions. The nanoscale confinement of most VDCCs and the effect of muscle contraction on their diffusion suggest that this dominant population of Ca2+ channels is primarily involved in electromechanical coupling at the sarcolemma.

Dystrophin impacts the nanoscale diffusion of sarcolemma VDCC

At the sarcolemma, the interaction of dystrophin with phospholipids, γ-actin and the dystrophin-glycoprotein complex26,27 could favour its role as a tension-sensing scaffold and a load-bearing apparatus capable of transducing muscle tension to membrane-embedded ion channels28,29. To investigate the role of dystrophin as a potential modulator of VDCC lateral mobility, we imaged dystrophin-deficient C. elegans dys-1(cx18) mutant worms expressing UNC-36-split-GFP (Fig. 6a) by SM dCALM.

Figure 6: Diffusion studies of individual VDCC at the sarcolemma of dys-1(cx18) mutants expressing UNC-36-split-GFP.
figure 6

(a) Immunofluorescence staining of GFP and dystrophin in normal and dys-1(cx18) worms expressing UNC-36-split-GFP. Scale bars: 5 μm. (b) Ensemble PDSD analysis of VDCC diffusive behaviours (17,273 VDCC, 32 muscle cells, five worms) in resting dys-1 mutant muscles. (c) Ensemble PDSD analysis of VDCC diffusive behaviours (40,068 VDCC, 33 muscle cells, six worms) in contracted dys-1 mutant muscles with levamisole. (d) Diffusion coefficients (±s.d.) for both VDCC populations in resting and contracted muscles of normal and dys-1 mutant worms (F-test, *P<0.05, ***P<0.001, NS: nonsignificant). (e) VDCC nanodomain nearest-neighbour distances at the sarcolemma of resting and contracted muscles for normal and dys-1 mutant worms. Error bars represent the s.d. of the mean (Wilcoxon sign-rank test, **P<0.01, NS: nonsignificant).

At the macroscale there was no obvious difference in the sarcolemma distribution of VDCC between dys-1 mutant and normal worms (Supplementary Movie 6). However, clear changes in sarcolemma scaffolding and fluidity were observed when we tracked individual channels. In resting dys-1 mutant muscles (17,273 trajectories, 32 muscle cells, five worms), the diffusion coefficient of slow VDCC (85%) was not significantly different from that of normal worms (Ddys-1 VDCC−slow=1.36 × 10−3±2.0 × 10−4μm2s−1, P=0.12, Fig. 6b,d); however, their confinement size was reduced to 58±2 nm in diameter (Supplementary Table 1), indicating that dystrophin is required to maintain the confining area of VDCC nanodomains but not their fluidity. The fast VDCC population (15%) diffused twofold faster than in normal worms (Ddys-1 VDCC−fast=4.54 × 10−2±1.2 × 10−3μm2s−1, P<0.05, Fig. 6b,d), an increased lateral mobility similar to that observed for other membrane proteins in muscle cells lacking dystrophin30. As dystrophin stiffens membranes that mimic the sarcolemma31, we attribute this apparent higher sarcolemma fluidity to a loss of membrane stiffness in dystrophin-mutant muscles. The free diffusion distance of fast VDCC was increased to 303±10 nm compared with resting normal muscles (Supplementary Table 1), suggesting that dystrophin is also involved in maintaining sarcolemma scaffolds outside VDCC nanodomains in body-wall muscles. VDCC nanodomains remained well dispersed but with a nearest-neighbour distance reduced to 606±273 nm compared with normal worms (P<0.01, Fig. 6e and Supplementary Fig. 6). This enhanced festooning of the sarcolemma and its reduced stiffness are both consistent with dystrophin being a load-bearing element capable of minimizing membrane strain by damping transverse cytoplasmic forces that push towards the sarcolemma under resting muscular tension. The smaller confinement area but unchanged diffusion coefficient of nanodomain-associated VDCC further suggests that, in normal resting C. elegans muscles, dystrophin does not bind directly to the channels but participates in modulating the nanodomain confinement size and the nanoscale spatial distribution of VDCC.

In dys-1 mutant worms anaesthetized with levamisole (40,068 trajectories, 33 muscle cells, six worms), the slow VDCC population (92%), which was nearly immobile in normal contracted muscles, displayed diffusive behaviours similar to that observed in resting muscles lacking dystrophin, with a confined diffusion within nanodomains 52±2 nm in diameter (DLevadys-1VDCC−slow=9.8 × 10−4±1.2 × 10−4μm2s−1, P=0.06,Fig. 6c,d and Supplementary Table 1). The fast VDCC population (8%) was less mobile by 35% compared with resting dys-1 mutant muscles (DLeva dys-1 VDCC−fast=2.90 × 10−2±5 × 10−4μm2s−1, P<0.05) and diffused freely over a nonsignificantly changed length scale of 281±6 nm (Fig. 6c,d and Supplementary Table 1). The nearest-neighbour distance between VDCC nanodomains was similar to that of resting dystrophin-mutant muscles at 628±257 nm (Fig. 6e and Supplementary Fig. 6). This apparent lack of additional sarcolemma festooning and the unchanged diffusive behaviour of nanodomain-associated VDCC in dys-1 mutant muscles under sustained cholinergic stimulation indicate that, in C. elegans, a functional dystrophin is required for an effective modulation of VDCC nanoscale dynamics at the sarcolemma during contractions.

Thus, dystrophin is a key scaffolding protein of the sarcolemma where it plays a dual role as a load-bearing apparatus and a tension transducer that has an impact on the nanoscale diffusion of membrane-embedded VDCC in response to changes in muscle tones. By showing that dystrophin is involved in modulating the confinement size of VDCC that normally generate action potentials and extracellular Ca2+ influx in muscle cells, our SM dCALM data suggest that it participates, directly or indirectly, in electromechanical coupling at the membrane of body-wall muscles in live C. elegans.

Discussion

Using far-field fluorescence microscopy and CALM, we achieved highly specific, SM sensitivity and nanometre accuracy imaging of membrane proteins within confined tissue regions in intact living C. elegans worms. We showed that microinjected CALM complementary peptides diffusing throughout the body of animals can activate fluorescent signals in engineered muscular and neuronal tissues by targeting and binding irreversibly to cell surface split-GFP fusion proteins. By tuning the concentration of injected peptides, ensemble or continuous SM CALM imaging could be performed by direct detection of activated split-GFP or by spFRET with targeted near-infrared fluorophores regardless of protein expression levels and of peptide nonspecific binding. Compared with SM imaging with conventional FP6, CALM provides low-background and membrane-specific SM detection without a need for photobleaching or alteration of protein expression and thus permits SM studies at physiologically relevant expression levels directly in live animals. While we primarly focused on membrane protein imaging, similar SM CALM detection of intracellular biomolecules may be achieved by co-expressing split-FP fusions with genetically encoded M3 peptide sequences. Chemically or light-inducible transcriptional elements might be used to control the expression of complementary peptides and maintain an SM detection regime. In principle, multiplex in vivo CALM imaging can also be realized with additional split-CFP and split-YFP fusions and by varying the fluorophores employed for CALM-spFRET.

In C. elegans, we tracked thousands of individual transmembrane proteins and channels with subresolution precision and quantitatively studied their diffusion in membrane micro- and nanodomains directly in live adult animals. In contrast to ensemble and diffusion-limited measurements by fluorescence recovery after photobleaching or FCS, single-particle tracking by CALM provides detailed analyses of protein locations and subpopulation dynamics at the nanoscale even for slow diffusing and confined membrane biomolecules. While our technique can be implemented in any animal species amenable to fluorescence imaging, C. elegans provides many advantages for in vivo SM studies. Its optical transparency, simple anatomy, small size and easy propagation facilitate sample preparation and peptide injection for SM CALM. These advantages, combined with the means to generate knock-in strains32,33 and the availability of large bioresources of C. elegans mutants open the possibilities of correlating the nanoscale behaviour of individual biomolecules imaged by CALM with phenotypic changes occurring at various developmental or disease stages, as illustrated by our study of VDCC dynamics in normal and dystrophin-mutant adult worms.

In this study, we reveal that dystrophin is a load-bearing scafolding protein of C. elegans body-wall muscles that not only influences membrane stiffness but also tranduces muscle tension to membrane-embeded VDCC. We show that, in response to changing muscle tones, dystrophin participates in modulating the nanoscale diffusion of two sarcolemmal populations of VDCC: a majority VDCC fraction undergoing slow and confined diffusion in ordered membrane nanodomains, and a minority fraction of channels diffusing faster and more randomly in other parts of the sarcolemma. As observed for comparable Ca2+ channels in cardiac myocytes34, the different diffusion characteristics and locations of these two VDCC populations suggest that they might play distinct roles for extracellular Ca2+ signalling on body-wall muscles. The fact that muscle contractions and the lack of functional dystrophin primarily influence nanodomain-associated VDCC implies that this subpopulation of Ca2+ channels and dystrophin itself are likely involved in electromechanical coupling in C. elegans, a process that requires precise spatial and temporal controls of Ca2+ influx in animal muscular tissues. The smaller, faster and more randomly distributed diffusing population of VDCC might regulate other Ca2+-dependent cellular processes in muscle cells.

The exact nature of VDCC nanodomains remains to be determined. While caveolae have been proposed to be sites of VDCC confinement in mammalian muscle cells35,36, there is as yet no evidence of caveolae formation at the sarcolemma of C. elegans, despite the expression of two caveolin isoforms37,38. Rather than caveolae, flat caveolin scaffolds39 or membrane rafts40 are potential nanodomains that could confine VDCC and induce their slow diffusion. Many C. elegans glycosylphosphatidylinositol-anchored proteins, such as UNC-36, co-fraction with caveolin in detergent-resistant membranes41 and in mammalian cells, α2δ subunits partition with other VDCC subunits in lipid rafts42. In cells, these membrane domains can slow down the lateral diffusion of proteins by a 10- to 20-fold ratio18,43. Our observations that the diffusion of VDCC inside or outside nanodomains differ by a comparable ratio suggest that they might be confined by similar raft-like domains at the sarcolemma.

We speculate that the confined diffusion of VDCC in sarcolemmal nanodomains is important for normal muscle functions in C. elegans. The confinement of VDCC and their near-immobilization during contractions in normal body-wall muscles might amplify small extracellular Ca2+ flux from action potentials into high local Ca2+ concentrations by nanodomain Ca2+ boosting44. This nanometre-scale localization of VDCC Ca2+ loads could reduce cross-talks between intracellular Ca2+ sources and increase the opening probability of Ca2+-dependent ryanodine/UNC-68 receptors in the closely apposed sarcoplasmic reticulum45 to potentiate muscle contractions46,47. The impaired modulation of VDCC confinement observed here in dystrophin-mutant muscles might additionally have an impact on nanodomain coupling between VDCC and SLO-1 current rectifying BK channels48,49 which, when defective, increases the excitability of muscle cells15 and induces the hyper-contracted phenotype of dys-1(cx18) worms50. Overall, it appears that dystrophin plays a key role for the nanoscale modulation of VDCC dynamics at the sarcolemma in C. elegans body-wall muscle cells. Such a function was not anticipated based on conventional imaging but might participate to the impairment of Ca2+ homeostasis and of normal muscle function in Duchenne muscular dystrophy51,52.

In conclusion, SM CALM imaging opens new avenues to explore the basic principles of homeostatic controls and the molecular basis of diseases at the nanometre scale in intact animal models by providing subresolution localization of individual membrane proteins and unprecedented insights into their diffusive behaviours within live tissues.

Methods

C. elegans strains

C. elegans strains were maintained on NGM agar plates at 20 °C on Escherichia coli OP50 lawns. N2 Bristol wild-type worms, CD4-split-GFP (kyEx1940[Pmyo-3::SL2::GFP1-10-CD4; Podr-1::SL2::DsRed]) worms a kind gift of Dr Cornelia Bargmann, strain EG4322 (ttTi5605II; unc-119(ed9)III) and unc-36(e251) worms from the Caenorhabditis Genetics Center, dys-1(cx18) worms a kind gift from Dr Kathrin Gieseler, UNC-36-GFP (krSi3[UNC-36::GFP]; e251) worms, UNC-36-split-GFP (krSi6[UNC-36::GFP 1-10]; e251) worms, UNC-36-split-GFP dystrophin-mutant (krSi6[UNC-36::GFP 1-10]; e251; cx18) worms and UNC-36-split-GFP/UNC-29-tagRFP (krSi6[UNC-36::GFP 1-10]; kr208(unc-29::tagRFP); e251) worms were used in this study.

Plasmid construction and C. elegans genetics

UNC-36-split-GFP was obtained by NheI and SacI insertion of the split-GFP-coding sequence (GFP 1–10) after the signal peptide sequence in unc-36 cDNA20 by PCR fusion. A synthetic C. elegans intron was also introduced upstream of GFP 1–10 to improve transgene expression. This transgene sequence was inserted in a pCFJ151 vector backbone (Addgene) to create the targeting vector pHZ043 and generate integrated animals by direct insertion on chromosome II using MosSCI protocols22. An EGFP sequence was fused to UNC-36 cDNA using the same approach, resulting in targeting vector pHZ025.

For MosSCI experiments, a mixture of plasmids including the targeting vectors pHZ043 or pHZ025 at 50 ng μl−1, the Mos transposase expression vector pJL43.1 at 50 ng μl−1 and the transformation markers Pmyo-3::GFP and Prab-3::GFP each at 10 ng μl−1 were co-injected in worm strain EG4322. MosSCI insertion progenies were screened by restoration of locomotion phenotypes and absence of transformation markers. MosSCI insertion was verified using PCR. MosSCI worms were then crossed with unc-36 null allele (e251) mutant worms to generate UNC-36-split-GFP and UNC-36-GFP worms. UNC-36-split-GFP worms were further crossed with a worm strain expressing the postsynaptic acetylcholine receptor fused to a red fluorescent protein (UNC-29-tagRFP) to generate worms dually expressing UNC-36-split-GFP/UNC-29-tagRFP. UNC-36-split-GFP worms were also crossed with dystrophin-deficient dys-1(cx18) mutant worms to generate dystrophin-mutant worms expressing UNC-36-split-GFP. All crosses were verified using PCR.

C. elegans phenotyping

Effective rescue of VDCC activity by integrated expression of UNC-36-GFP and UNC-36-split-GFP was verified in thrashing assays where 20 adult worms were transferred separately in a drop of M9 medium at 20 °C and thrashes were counted in 1-min intervals after 2 min recovery. Data are reported as the mean number of body bends per min±s.d. of the mean. Statistically significant differences are estimated by a Welch two sample t-test.

Electrophysiology

Following microdissection of C. elegans, electrophysiological recordings of membrane currents and potentials were performed in the whole-cell configuration using an RK-400 patch-clamp amplifier (Bio-Logic)53. Acquisition and command voltage were controlled using the pClamp9 software driving a 1322A Digidata (Molecular Devices). Data were analysed and graphed using Origin software (OriginLab). The resistance of recording pipettes ranged between 2.5 and 2.8 MΩ. Recordings were performed after 1 min dialysis only on cells exhibiting resistances above 800 MΩ. Capacitance and resistance were not compensated.

The bath solution contained 140 mM Tetraethylammonium chloride (TEACl), 6 mM CaCl2, 5 mM MgCl2, 3 mM 4-Aminopyridine (4-AP), 10 mM Hepes and sucrose at 337 mosm l−1 (pH 7.2). The pipette solution contained 140 mM CsCl, 5 mM TEACl, 5 mM EGTA, 5 mM Hepes, 4 mM MgATP and sucrose at 328 mosm l−1 (pH 7.2). The membrane potential was hold at −60 mV, and currents were recorded by applying 200-ms voltage steps from −70 to +70 mV with 10-mV increments. Leak currents were subtracted from all recordings. Current–voltage relationships for N2 wild-type (n=8), unc-36(e251) (n=7), UNC-36-GFP (n=7) and UNC-36-split-GFP (n=8) worms were established by measuring the currents at the peak of the current and were fitted with the following equation:

where I(v) is the measured density of current, v the test voltage pulse, Gmax the maximum conductance, Vrev the reversal potential, V0.5 the half-activation voltage and k the steepness factor. Potentials of half-activation, E0.5 were obtained by fitting, for each muscle cell, the current–voltage relationships of the inward currents measured at the peak of the currents by equation (1). Time-to-peak kinetic was measured for the maximal current traces (+20 mV). The mean E0.5 and the mean time-to-peak values are reported with error bars representing the s.e.m.

All chemicals were obtained from Sigma-Aldrich and all experiments were performed at room temperature.

Immunohistochemistry

Mix-stage worms were freeze-fractured and fixed in methanol/acetone according to freeze fracture protocols54; however, only adult worms were imaged. For immunodetection of split-GFP fusion proteins, primary rabbit anti-GFP antibodies (Invitrogen A11122, dilution 1/500) and secondary polyclonal goat anti-rabbit Alexa-488-conjugated antibodies (Invitrogen A11008, dilution 1/800) were used. For immunostaining of dense bodies, mouse MH25 antibodies against PAT-3 β-integrin55 (dilution 1/200) and secondary polyclonal goat anti-mouse antibodies conjugated to Cy3 (Invitrogen A10521, dilution 1/500) were used. For immunostaining of dystrophin, mouse monoclonal antibodies against C. elegans dystrophin (BioTem 1H7-1B3-2C7-3F5, dilution 1/200), a kind gift from Dr Kathrin Gieseler, and secondary polyclonal goat anti-mouse antibodies conjugated to Cy3 (dilution 1/500) were used.

Worm microinjection

Worms were placed on dried 2% agarose injection pads in a drop of mineral oil (Sigma, H8898) and were forced to adhere to the pads by gentle brushing of the animal body. For ensemble dCALM and CALM-spFRET imaging, purified synthetic complementary M3-biotin and M3-A647 peptides12 were diluted in PBS at concentrations of 5 mM and 200 μM, respectively, before being microinjected into the worms pseudoceolomic cavity in the head region of the animals using a microinjector (Femtojet, Eppendorf) and FemtotipII microinjection capillaries. Worms were microinjected in three to five steps of 0.3 s with a few hundreds of femtolitres of the peptide solutions. Worms were allowed to recover in an M9 drop on NGM plates at 15 °C for 1 h and were then transferred to 20 °C for at least 3 h before imaging. For SM dCALM and SM CALM-spFRET of CD4-split-GFP, microinjections were performed with M3-biotin and M3-A647 peptide concentrations reduced to 200 and 20 μM, respectively, and worms’ recovery times were limited to 30–60 min before imaging. For SM dCALM of VDCC, worms were injected with higher M3-biotin concentrations than for SM dCALM of CD4-split-GFP and longer recovery times were employed to compensate for the low endogenous membrane expression of UNC-36-split-GFP. M3-biotin peptides at 1 mM were microinjected and worms were allowed to recover for 3–12 h at 15 °C before imaging. For in vivo membrane staining of UNC-36-GFP, anti-GFP Alexa-594 antibodies (Invitrogen A-21312) were diluted 200-fold in injection buffer (20 mM K2HPO4, 3 mM K+ citrate, 2% PEG 3000, pH=7.4) and were injected into the pseudocoelomic cavity56 in the head region before recovery in an M9 drop on a NGM plate at 20 °C for 3–6 h and imaging.

Ensemble microscopy imaging

Ensemble GFP, dCALM and CALM-spFRET imaging of live worms were performed on a confocal Leica TCS SP5AOBS microscope or a Leica DM5000B microscope equipped with a CSU10 spinning disc confocal scanner system (Yokogawa). Appropriate laser excitation wavelengths and emission filtering were used to image GFP and activated split-GFP fusions (488 nm excitation), anti-GFP antibodies or UNC-29-tagRFP (561 nm excitation) and M3-A647 peptides (633 nm excitation). FRET from complemented CD4-split-GFP to M3-A647 was detected by single laser excitation at 488 nm and simultaneous detection in the GFP and the A647 emission channels. Live worms were anaesthetized with 0.2% sodium azide or 50 mM muscimol before being mounted between 2% agarose pads and microscope coverslips for imaging. Freeze-fractured and immunostained worms were imaged on the same microscopes, except for immunostained CD4-split-GFP worms that were imaged by wide-field epifluorescence.

Single-molecule microscopy imaging

SM dCALM and SM CALM-spFRET of CD4-split-GFP in live worms were performed on an inverted IX70 Olympus microscope equipped with custom-built total internal reflection (TIR) and highly incline and laminated optical (HILO) elements, a × 100 1.45 NA objective, a DV2 dual-view system (Photometrics), a QuantEM:512SC EMCCD (Photometrics) and appropriated filters12. SM dCALM imaging of VDCC and simultaneous dual-colour imaging of UNC-29-tagRFP was performed on an inverted Nikon Eclipse Ti-E microscope equipped with TIR and HILO optics, 488, 561 and 647 nm fibre-coupled excitation lasers, a × 100 1.49 NA objective, a two-camera imaging splitter (Andor) and two iXon EMCCD cameras (Andor). A multiband pass ZET405/488/561/647x excitation filter (Chroma), a quad-band ZT 405/488/561/647 dichroic mirror (Chroma), an emission splitting FF560-FDi01 dichroic mirror (Semrock) and two emission filters at 525/50 nm (Semrock) and 600/50 nm (Chroma) for GFP and tagRFP, respectively, were used.

HILO imaging of worms was achieved by moving circularly polarized laser excitation beams away from the optical axis at the back-focal plan of the objectives and at a position slightly off from the location normally required to obtain a critical TIR angle. For both optical set-ups, the size of pixels in images was determined by imaging a micrometre reticle and the alignment of emission channels was performed by imaging 40-nm diameter TransFluoSphere beads (488/685 nm, Invitrogen).

Non-anaesthetized worms (resting muscles) were deposited on 2% agarose pads in a drop of M9 containing 0.1 μm polystyrene beads at 2.5% per volume (Polyscience) and were mounted with a microscope coverslip57. Tracking was only performed in-between sporadic muscle contractions. Worms anaesthetized with levamisole (contracted muscles) were first deposited in a drop of 1 mM levamisole (Sigma) for 10–15 min before being mounted in levamisole between 2% agarose pads and coverslips. In each worm, multiple muscle cells were imaged at room temperature by continuous laser excitation and with camera acquisition rates at 100 ms per frame for CD4-split-GFP and 80 ms per frame for VDCC.

Tracking and diffusion analyses

SM tracking and diffusion analyses were performed using SlimFast, a single-particle detection and tracking software written in Matlab and based on multiple-target tracing algorithms58, which was generously provided by Christian Ritcher and Jacob Piehler. We used parallel Matlab on CPU clusters to localize, track and perform diffusion analysis on tens of thousands of individual split-GFP fusions in multiple worms. SM tracking in dCALM movies was performed by 2D Gaussian fitting of the point-spread function for each activated CD4-split-GFP or VDCC in each frame. The mean localization precision±s.d. for all individual fusion protein detected in CALM movies was determined by the method of Thompson et al.59 Trajectories were built by linking individual localized positions from one frame to the other, taking into account blinking statistics and local particle densities. Only the trajectories with at least three-step sizes were kept for ensemble and individual mean square displacement (MSD) analyses and for probability distribution of the displacement squared (PDSD) analyses.

For diffusion analyses, MSD and PDSD curves were computed17,43. For PDSD, computations were performed on the first ten observation time lags, a number sufficient to detect the confinement of CD4-split-GFP or VDCC on C. elegans muscles under our imaging conditions. In brief, for the first ten time lags t, each Pr2 curve was fitted with the general model:

where the fitting coefficient ri2(t) and αi(t) are the square displacement and the fraction corresponding to i numbers of diffusion regimes at each time lag t, respectively. Pr2 distribution for CD4-split-GFP was best fitted with i=3 diffusion regimes, while Pr2 distribution for VDCC were fitted with i=2 regimes.

Error bars in MSD curves correspond to the s.e.m. Error bars for each ri2 in ri2(t) curves were determined using , where N represents the number of independent data points used to build each probability distribution functions17,60.

Diffusion coefficients were obtained by fitting MSD or ri2(t) curves using:

(i) an anomalous diffusion model:

where α is the anomalous exponent and Dα is the time-dependent diffusion coefficient;

(ii) A free Brownian diffusion model with measurement error:

where σ is the position error and D is the diffusion coefficient;

or (iii) a circularly confined diffusion model with measurement error:

where R is the confinement radius, σ is the position error, D is the diffusion coefficient, A1=0.99 and A2=0.85 (ref. 43)43.

Fitting was performed using Origin software (OriginLab). When diffusion was found to be anomalous over the first ten time lags on ri2(t) curves, the linear part of the curves were fitted with the free Brownian diffusion model to assess short-scale diffusion coefficients and free diffusion distances. To build the distribution of diffusion coefficients for individual VDCC, the free Brownian diffusion model was fitted over the first three points of each individual MSD. All the diffusion coefficients D are reported in micrometre squared per second±s.d. of the fit value. Significant differences between diffusion coefficients were determined using F-tests.

Intensity time trace analyses on individual VDCC were performed by integrating the fluorescence intensity within a 3 × 3 pixel region centred on the point-spread function of each activated VDCC and for the whole length of acquired dCALM movies, including periods with no GFP signal before activation and after photobleaching. Background correction was performed by subtracting the fluorescence intensity at each time point from a similar 3 × 3 pixel region in the immediate vicinity of each individual VDCC. The slight differences in intensity between time traces are because of moderate variations in HILO excitation field uniformity when imaging different muscle cells or nerve cords in worms.

Distribution pattern analyses of VDCC nanodomains

Ripley’s K function analyses and nearest-neighbour distance measurements were performed with the software SpPack61 after 2D Gaussian fitting of the nanodomain positions in sum intensity projection images of SM dCALM movies.

Ripley’s K analyses with edge corrections were performed on muscle cells for multiple regions of interest totaling 650–2,000 μm2 of the sarcolemma membrane in multiple worms for each condition. 95% confidence intervals for complete spatial randomness were obtained by averaging Monte Carlo simulations consisting of 500 replicates in multiple regions of interest with areas and nanodomain seeding numbers similar to that of the worm data.

For nearest-neighbour distance measurements, 1,100–2,500 nanodomains were analysed for each condition and data are reported as the mean nearest-neighbour distance in nanometre±s.e.m. A non-parametric Wilcoxon sign-rank statistical test was performed to assess significant differences between data sets.

Additional information

How to cite this article: Zhan, H. et al. In vivo single-molecule imaging identifies altered dynamics of calcium channels in dystrophin-mutant C. elegans. Nat. Commun. 5:4974 doi: 10.1038/ncomms5974 (2014).