# Elementary calcium release events in the skeletal muscle cells of the honey bee Apis mellifera

Claude Collet, Mercédès Charreton, Laszlo Szabo, Marianna Takacs, Laszlo

Csernoch, Peter Szentesi

► To cite this version:

Claude Collet, Mercédès Charreton, Laszlo Szabo, Marianna Takacs, Laszlo Csernoch, et al. . Elementary calcium release events in the skeletal muscle cells of the honey bee Apis mellifera. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-96028-w. hal-03334447

HAL Id: hal-03334447

https://hal.science/hal-03334447v1

Submitted on 3 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Elementary calcium release events in the skeletal muscle cells of the honey bee Apis mellifera

Claude Collet<fnref n="1" />✉, Mercédès Charreton<fnref n="1" />, Laszlo Szabo², Marianna Takacs³, Laszlo Csernoch⁴ & Peter Szentesi⁴

Calcium sparks are involved in major physiological and pathological processes in vertebrate muscles but have never been characterized in invertebrates. Here, dynamic confocal imaging on intact skeletal muscle cells isolated enzymatically from the adult honey bee legs allowed the first spatio-temporal characterization of subcellular calcium release events (CREs) in an insect species. The frequency of CREs, measured in x-y time lapse series, was higher than frequencies usually described in vertebrates. Honey bee CREs had a larger spatial spread at half maximum than their vertebrate counterparts and a slightly ellipsoidal shape, two characteristics that may be related to ultrastructural features specific to invertebrate cells. In line-scan experiments, the histogram of CREs' duration followed a bimodal distribution, supporting the existence of both sparks and embers. Unlike in vertebrates, embers and sparks had similar amplitudes, a difference that could be related to genomic differences and/or excitation-contraction coupling specificities in honey bee skeletal muscle fibres. The first characterization of CREs from an arthropod which shows strong genomic, ultrastructural and physiological differences with vertebrates may help in improving the research field of sparkology and more generally the knowledge in invertebrates cell Ca²⁺ homeostasis, eventually leading to a better understanding of their roles and regulations in muscles but also the myotoxicity of new insecticides targeting ryanodine receptors.

Bees are very efficient pollinator insects, and are thus crucial both for plant biodiversity maintenance and agricultural food production. Amongst bees, owing to their easy domestication, honey bees have long been used as a research model in various fields, ranging from cell physiology to inter-individual communication and sociability<fnref n="1" />,². Moreover, the high levels of mortality in honey bee's populations observed worldwide in the last few years prompted researchers to further study this model with the hope of understanding origins of this decline. Honeybee cells provide valuable models to study calcium signalling and intracellular calcium homeostasis in insects³⁻¹¹. Calcium is an essential intracellular second messenger which plays pivotal roles in neuromuscular cellular processes such as signal transduction and modulation in motor or sensory neurons, contraction in cardiac, smooth, and skeletal muscle cells. More generally, Ca²⁺ plays a major role in development/maturation, cell shape modulation and migration, cell division, cell-membrane repair, hormonal secretion, and gene expression in virtually all cell types. In the honey bee, macroscopic intracellular Ca²⁺ variations have been recorded with Ca²⁺-imaging techniques in neurons from central brain structures such as mushroom bodies³,⁴ or antennal lobes⁷,⁸. Recently, Ca²⁺ was also monitored in peripheral sensory neurons from the antennae, i.e. the 'nose' of honeybees¹⁰. Global Ca²⁺ fluctuations were also studied in photoreceptor cells from the honey bee compound eye¹¹. Finally, cytoplasmic Ca²⁺ transients (i.e. transient Ca²⁺ increases that can be elicited by a stimulation) and excitation-contraction coupling have been described in bee skeletal muscle cells as well, using fluorescent dyes and electrophysiological techniques⁵,⁶. Recently, studying intracellular Ca²⁺ signalling in bees and pollinators gained renewed interest because new classes of insecticides targeting several molecular players involved in Ca²⁺ homeostasis have been launched on the pesticides market and are now widely used in agriculture. Notably, most recent diamide type insecticides such as flubendiamide and chlorantraniliprole induce detrimental effects in

<footnote n="1">¹UR406 Bees and Environment, French National Institute For Agriculture, Food and Environment, INRAE, 84914 Avignon, France. ²Department of Electrical Engineering, Sapientia Hungarian University of Transylvania, Tîrgu Mureş, Romania. ³Department of Animal Science, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, Debrecen, Hungary. ⁴Department of Physiology, Medical Faculty, University of Debrecen, Debrecen, Hungary. ✉email: claude.collet@inrae.fr</footnote>
honey bees by anarchically mobilizing $Ca^{2+}$ from intracellular stores through $Ca^{2+}$ release channels in neurons from the antenna and in skeletal muscle cells from the leg⁹,¹⁰.

A major molecular player in $Ca^{2+}$ signalling is undoubtedly the ryanodine receptor. The molecular identity of this receptor was unravelled thanks to the botanical alkaloid ryanodine, present in root extracts of the plant species *Ryania*, which was used in the 50's as a natural insecticide (now withdrawn owing its toxicity towards environment). Ryanodine receptors are huge macromolecular complexes functioning as ion channels to provide a route for $Ca^{2+}$ to be massively transferred on demand from intracellular storage organelles (sarco and endoplas- mic reticulums) to the cytoplasmic compartment¹,¹³. Increasing the intracellular cytoplasmic $Ca^{2+}$ concentration ($[Ca^{2+}]_i$) from submicromolar to tenths of micromolar ranges thus allows $Ca^{2+}$ to play its multiple regulatory roles¹⁴,¹⁵. At rest, $[Ca^{2+}]_i$ is tightly regulated and maintained at around 100 nM to avoid the toxicity of sustained $Ca^{2+}$ increases and this low resting $[Ca^{2+}]_i$ provides both a wide dynamic range and a high signal-to-noise ratio. Optical and electrophysiological techniques evolved and allowed characterizing cytoplasmic $Ca^{2+}$ fluctuations with continuously increased acuity along with recent years. Nowadays, laser-scanning confocal microscopes allow detecting optically subcellular cytoplasmic $Ca^{2+}$ fluctuations in the µm spatial and in the ms temporal ranges with the help of fast and high-affinity $Ca^{2+}$ dyes ($Ca^{2+}$ indicators). These localized fluctuations are believed to originate from the release of $Ca^{2+}$ by one or a cluster of several ryanodine receptors and in striated muscles, they are considered as the unitary bricks of voltage-induced global $Ca^{2+}$ release¹⁶. Localized elevation in $[Ca^{2+}]_i$ can also fuse to initiate propagating macroscopic $Ca^{2+}$ waves and have specific cytosolic localizations such as the perinuclear space where they can modulate nuclear $Ca^{2+}$ independently of the bulk cytosolic $Ca^{2+}$¹⁷. So far, local fluctuations in $[Ca^{2+}]_i$ have been studied only in muscle cells from vertebrate species. First discovered in cardiac myocytes¹⁸, and originally termed $Ca^{2+}$ sparks, similar signalling events were identified later in smooth muscle cells, skeletal muscle fibres, neurons and oocytes. The role and regulation of mechanisms of these spatially localized fluctuations are however elusive and remain to be fully identified with a rather unique research field dedicated to their study: 'sparkology'²ⁱ,²⁷. This research field is essential as these signals are frequently altered in genetic diseases linked to mutations in the ryanodine receptor genes (ryanodinopathies). Even under physi- ological conditions, it soon turned out that localized $Ca^{2+}$ fluctuations had many temporal patterns, not only an archetypal steep increase (~5 ms) and a fast mono-exponential decrease, and new terms emerged to better describe this signals diversity e.g. embers, macrosparks, compound sparks, bursts, sparklets...¹⁹. A more general term is now used to depict these $[Ca^{2+}]_i$ fluctuations : elementary $Ca^{2+}$ release events (CREs).

While CREs have been identified in many tissues from several vertebrate species (rat, mouse, frog and fish), to our knowledge, they have never been identified in an invertebrate tissue. Here, using honey bee as a model, we describe for the first time elementary CREs in an insect species and compare their spatiotemporal characteristics in the light of interspecific physiological and ultrastructural cellular differences. Part of this work was presented to the Biophysical Society²⁰.

## Material and methods

**Isolation of leg muscle cells.** Muscle fibres were prepared from young adult domestic honey bees (*Apis mellifera*), as described earlier⁵, with slight modifications. Newborn bees were collected outdoor, directly on a frame of a hive and used during the following 4 days. Bees were chilled (~ 4 °C) until lethargic immobility and then killed by decapitation. Tibias from the third pair of legs were dissected in a Tyrode's solution without cal- cium (in mM: NaCl 140, KCl 5, MgCl₂ 2, HEPES 10, pH 7.2), cut longitudinally and pinned open. Tibias were bathed 20 min in a dissociation solution at 37 °C (Tyrode's solution without calcium containing four enzymes: collagenase type 1, trypsin, papain, and pronase at concentrations 0.5, 1.5, 1.5 and 1 mg/ml, respectively). After rinsed in Tyrode's solution without calcium and normal Tyrode (2 mM CaCl₂), muscle masses were triturated with a pipette and dispersed in the central well of a clean glass bottom dish (MatTek Corp., Ashland, MA, USA). Cell were loaded for 30 min with the cell permeant dye Fluo-8-AM (5 µM).

**Measurement of $Ca^{2+}$ release events (CREs) in intact bee muscle fibres.** Fluo-8-AM loaded cells were kept in normal Tyrode's solution and observed with a laser scanning confocal microscope Zeiss LSM 5 LIVE (Zeiss, Oberkochen, Germany) equipped with a Plan-Neofluar 40× oil immersion objective (NA = 1.3). Fluo-8 was excited with the 488 nm wavelength of an argon laser and the emitted fluorescent light was meas- ured through a band-pass filter (505–570 nm) and digitized at 12 bit. First the fibres were examined by eye and the fibres showing spontaneous activity were recorded with series of 512×512 (XYT) images, captured at 1 Hz. The pixel size was 0.328 µm and the pixel dwell time was 1 ms. CRE detection and analysis were performed using methods and algorithms previously described²¹,²². In brief, the region of the fibre and the background was defined on the first image of each series. The averaged background fluorescence was subtracted from each pixel of all images in one series. $Ca^{2+}$ release events were detected by the stationary wavelet method. The filtering was made by soft thresholding wavelet detection. Finally the amplitude and full width at half maximum (FWHM) of the sparks were calculated. Two FWHM values were determined: perpendicular to (X) and parallel (Y) with the Z-lines. To detect the Z-lines the frequency spectrum in each line of all images was calculated using fast Fourier transform. To remove the Z-lines from the images, inverse FFT of the frequency components corresponding to sarcomeres was used. Bleaching correction was performed on all images by calculating the change in fluo- rescence intensity of each pixel in time. Finally images were normalized using stationary wavelet transform. The event frequency was determined by counting the number of events on all images from a series. Events that appeared on consecutive images at the same location were counted as a single event. To get the frequency the number of events was normalized to fibre area and duration of recording. To assess the amount of $Ca^{2+}$ released during a CRE, the Signal Mass (SM) was calculated as described earlier²³. To estimate FWHM in the z direction (i.e. the optical axis) in two-dimensional measurements, the average of the two FWHMs in the focal plane (x
<figure><img src="image_3.png" /><figcaption>**Figure 1.** General structure of bee muscle fibres. (A) Transmitted light for a honey bee tibia muscle fibre (mean sarcomere length 4,9 µm). (B) Confocal longitudinal section through the fluorescent central row of nuclei labelled with Hoechst stain (same fibre as in A). (C) Transmitted light for a Fluo-8 stained fibre (mean sarcomere length 4,4 µm). (D) Confocal fluorescence image of the fibre shown in C. Scale bars = 20 µm.</figcaption></figure>

and y directions) was used. Figure 1 and Supplementary Video 1 were made under the ImageJ software (v1.52a, Rasband, W.S. ImageJ, U.S. National Institutes of Health, Bethesda, Maryland, USA, https://imagej.nih.gov/ij/, 1997–2018).

Line-scan images (XT) were taken at 1 or 2 ms/line and the pixel size was 0.328 µm. In some cases 0.2 ms/line was used to study the temporal properties of CREs in high time resolution. The length of the line was variable and depended on the actual fibre width. The images (F[x,t]) were normalized to baseline fluorescence (F₀[x]). Line-scan images containing spontaneous calcium release events were analysed using an automatic event detection program²¹, which calculated the amplitude (ΔF/F₀), rise time, full width at half maximum (FWHM), and full duration at half maximum (FDHM) of the identified events. An event was identified as an ember if it had a clear plateau and was longer than 40 ms in time. Otherwise, it was considered as a spark. Events with an amplitude smaller than 0.1 ΔF/F₀ units were excluded from the analysis. All measurements were carried out at room temperature (22–23 °C).

**Depolarization-evoked Ca²⁺ transients measurements.** Ca²⁺ transients were elicited on CRE free (silent) intact fibres with field stimulation and recorded in line-scan mode. Two platinum electrodes were placed adjacent to the fibres and 50 ms long depolarizations with amplitude ranging between 40 and 80 V were used to evoke Ca²⁺ transients (S88 Stimulator, Grass Technologies, Warwick, RI, USA). Whole-cell Ca²⁺ transients were recorded in patch-clamp experiment as described earlier⁵.

**Nuclei staining.** Fibres were stained with Hoechst 33258 (5 µg/ml) in Tyrode during 30 min. After being loaded, fibres were rinsed twice in Tyrode's solution and nuclei were visualized under a Zeiss LSM710 confocal inverted microscope with a 63 × 1.42 NA objective (Zeiss, Oberkochen, Germany).
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Spontaneous elementary calcium release events in two dimensional (XY) experiments. (A) Sum of 30 background corrected images of a Fluo-8 loaded fibre representing 30 s. (B) A representative spark. Histograms of parameters of calcium sparks: amplitude (C), FWHM perpendicular (D) and parallel (E) to the fibre axis. Calibration on panel A is 50 µm. Calibrations on panel B are 5 µm.</figcaption></figure>

**Resting membrane potential measurement.** Single muscle fibres were superfused with normal Tyrode's at room temperature (24 °C). Membrane potentials were recorded using 3 M KCl filled sharp glass microelectrodes with tip resistances ranging between 30 and 40 MΩ only in cells with intact surface membrane and clear cross striations. The electrodes were connected to the input of a Multiclamp-700A amplifier (Axon Instruments, Foster City, CA, USA) under current-clamp conditions. Membrane potentials were digitized at 100 kHz using Digidata 1322 A/D card (Axon Instruments) under software control (pClamp 9.2, Axon Instruments, Foster City, CA, USA) and stored for later analysis.

**Statistical analysis.** Pooled data were expressed as mean ± standard error (SE) of the mean. The differences were assessed using one-way analysis of variance (ANOVA) and all pairwise multiple comparison procedures (Student–Newman–Keuls method). A p-value of less than 0.05 was considered statistically significant. Graphpad Prism 6.0 software (GraphPad software, San Diego, USA) was used for statistical analysis and non-linear regression.

**Results**

Honey bee intact skeletal muscle fibres isolated from the legs present a nice sarcomere striation under light microscope and the nuclei are located in the centre of the fibres⁵ instead of right underneath the sarcolemma as in other adult organisms like amphibians and mammals (Fig. 1A, B, Hoechst staining). The axial row of nuclei give them a so-called ‘tubular’ appearance in histological transverse sections. The sarcomere pattern in healthy muscle fibres used in this study showed a period of ~ 4 µm, indicating a relaxed status for the contractile machinery. Their membrane potential, measured with sharp microelectrodes (-58.0 ± 6.9 mV, n = 3) was consistent with measurements made in other insect species²⁴ and in enzymatically-isolated mammalian muscle fibres that are polarized around -60 mV²⁵⁻²⁷. An extensive t-tubular organization has been demonstrated earlier in these non-flight muscle fibres, as well as functional Ca²⁺ and K⁺ ionic channels and excitation-contraction coupling⁵,²⁸. Here, examining Fluo-8-loaded leg fibres in confocal fluorescence (Fig. 1C, D) revealed that around one fifth of fibres produced spontaneous calcium release events when bathed in a physiological extracellular solution (Tyrode's solution).

In spontaneously active fibres, localized calcium release events appeared randomly in time and space (Supplementary video 1, Fig. 2A). The events appearing in XYT image sequences showed classical spark shape (Fig. 2B). The analysis of a number (4898) of CREs revealed that their average amplitude is quite small (0.220 ± 0.001 ΔF/F₀), but some higher amplitude events (up to 0.6 ΔF/F₀) were observed as well (Fig. 2C). This amplitude can be compared to ~1–2 ΔF/F₀ global Ca²⁺ transients we observed during field stimulations (Supplementary Fig. S1) and previously measured under voltage-clamp in bee muscle cells loaded with Fluo-3⁵. The average frequency of CREs was high (1696 ± 367 Hz/mm², n = 25 fibres). The CREs were wide in space since both FWHMs (X and Y) were greater than 3 µm in average (3.73 ± 0.02, and 3.30 ± 0.02 µm, respectively). The probability histograms of FWHMs (Fig. 2D and E) show that some events had a broad size (> 10 µm). The shape of honey bee CREs in space were not completely circular. The average ratio of FWHMs (X/Y) were significantly greater than 1
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Sparks and embers on line-scan image. (A) Representative background corrected line-scan image measured with 2 ms/line scanning speed. The trace below the image was calculated by averaging 5 lines at the position marked by the arrow. Enlarged image of a spark (B) and an ember (C) with its spatial and time profile. The positions of the events are marked with yellow rectangular in panel A.</figcaption></figure>

<table><caption>**Table 1.** Characteristics of honey bee elementary calcium release events detected in line-scan (XT) experiments. Events with less than 0.1 amplitude were discarded. All events with greater than 40 ms duration were classified as ember. The numbers in parenthesis are the number of event analysed. ***Denotes significant differences from spark, p < 0.001. FWHM – full width at half maximum, FDHM – full duration at half maximum.</caption><thead><tr><th>(22 fibres)</th><th>Amplitude (ΔF/F<sub>0</sub>)</th><th>Rise time (ms)</th><th>FWHM (μm)</th><th>FDHM (ms)</th><th>Signal mass (μm<sup>3</sup>)</th><th>Frequency (mm<sup>-1</sup> s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>All events (3338)</td><td>0.185 ± 0.001</td><td>15.88 ± 0.20</td><td>2.82 ± 0.02</td><td>35.82 ± 0.34</td><td>7.41 ± 0.19</td><td>265.1 ± 15.0</td></tr><tr><td>Sparks (2244)</td><td>0.185 ± 0.001</td><td>10.45 ± 0.11</td><td>2.65 ± 0.02</td><td>25.52 ± 0.16</td><td>6.00 ± 0.20</td><td>174.0 ± 14.4</td></tr><tr><td>Embers (1094)</td><td>0.186 ± 0.002</td><td>27.02 ± 0.41***</td><td>3.18 ± 0.04***</td><td>56.95 ± 0.58***</td><td>10.29 ± 0.41***</td><td>63.1 ± 6.6***</td></tr></tbody></table>

(1.22 ± 0.01; p < 0.001). The signal mass (SM), calculated to estimate the amount of Ca²⁺ released during a CRE, was 14.52 ± 0.29 µm³ on average.

To study elementary CREs in more details, we performed line-scan experiments. Figure 3A shows a representative background-corrected line-scan image in pseudo colour. This fibre was very active, and several sparks (Fig. 3B) and long lasting embers (Fig. 3C) were observed in this record. Most CREs appeared randomly in space and time. However, there were some frequently firing sites which generated more than one event in the same place. The trace calculated from Fig. 3A at the position marked by an arrow, shows a representative example of this phenomenon. The time profile of the active calcium release unit demonstrates that the spontaneous event originating from the same place can be a spark, an ember, or a burst²⁹.

First, the spatiotemporal parameters of CREs were calculated supposing a homogeneous population (Table 1, Fig. 4). However, as it is clearly visible in Fig. 3 there are some long lasting events, called traditionally embers³⁰. The histogram of duration of all CREs was indeed well fitted by the sum of two Gaussian functions (Fig. 4D). The peak position of the first Gaussian was found at 25.5 ± 0.5 ms, while the second function peaked at 44.3 ± 3.9 ms. The intersection of the two Gaussian functions was at around 40 ms and this value was thus chosen as the threshold of duration to distinguish a spark from an ember (see Table 4, column entitled “Threshold duration of embers”). All CREs with a duration longer than 40 ms were categorized as embers.

Based on this classification we found that 67% of the events were sparks, with a high frequency (Table 1). The average amplitude of sparks was 0.185 ± 0.001 ΔF/F₀, but the histogram of amplitudes shows some events with higher amplitudes up to 0.4 ΔF/F₀ (Fig. 5A). The average FWHM of sparks was found to be 2.65 ± 0.02 µm,
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Parameters of calcium release events. Histogram of amplitudes (A), FWHMs (B), rise times (C), and durations (D). The number of calcium release events was 3338. The parameters of Gaussian fits in panel D are 0.109 and 0.076 for the amplitude and 25.54 and 44.26 s for the mean to spark (red) and embers (blue), respectively. Pink line represents the sum of the two (red and blue) Gaussian functions.</figcaption></figure>

interestingly we found some very wide sparks with a FWHM greater than 7 µm (Fig. 5B). The distribution of rise times clearly shows that the majority of sparks developed in less than 15 ms (Fig. 5C). The average duration was identical with the first peak position of the theoretical Gaussian function (25.52 ± 0.16 ms). The histogram of duration did not show the classical bell shaped curve since the upper limit was chosen arbitrarily (Fig. 5D).

We classified 33% of events as embers with significantly lower frequency than sparks (Table 1). Their average amplitude was identical to that of sparks (0.186 ± 0.002 ΔF/F₀; p > 0.6). The maximum of ember amplitude was 0.5 ΔF/F₀ (Fig. 6A). The average FWHM was significantly larger than for sparks (Table 1), likely because calcium had more time to diffuse during embers (Fig. 6B). The average rise time was also significantly longer than for sparks (Table 1) showing that the development of embers is a slower process (Fig. 6C). The average duration was somewhat bigger than the second peak position of the theoretical Gaussian function (56.95 ± 0.56 ms). We found some very long events and more than 5% of the embers were longer than 100 ms (Fig. 6D).

High speed confocal scanning was used to study in more details the time profile of CREs (Fig. 7A). This speed (200 µs/line) allowed us to explore the rising and the falling phases of the events. We could confirm our observation that the amplitude of sparks and embers are roughly the same in honey bee muscle. Furthermore if one compares the time profile of a spark (Fig. 7B) to that of an ember (Fig. 7C), it is clear that embers need more time to reach their peak.

## Discussion

In the honey bee leg skeletal muscle, excitation-contraction coupling relies on a calcium entry through plasma membrane voltage-gated calcium channels, which activates ryanodine receptors to release more Ca²⁺ from intracellular stores⁵,⁶. Whereas macroscopic calcium signalling has been studied in invertebrate species, subcellular calcium release events such as Ca²⁺ sparks have never been described in detail. Our work may help understanding
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Parameters of calcium sparks. Histogram of amplitudes (A), FWHMs (B), rise times (C), and durations (D). The number of sparks was 2244.</figcaption></figure>

the role and regulation of these events, since interspecific comparisons often lead unravelling cellular phenomena. Data collected in the present study were compared to all available data on CREs recorded in intact freshly isolated skeletal muscle fibres and in intact cell-cultured myotubes, gathered in Tables 2–4. In our comparison, ‘intact’ refers to the integrity of the plasma-membrane, without any attempt at manipulating the intracellular cytoplasmic composition and with physiological extracellular solutions bathing the cells. In honey bee fibres, both sparks (i.e. CREs of short duration) and embers which are CREs of long duration first characterized in adult frog muscle cells³⁰, were recorded. Similar short- or long-lasting events have been characterized, for instance in adult mammalian fibres subjected to surface membrane skinning³¹ or in intact foetal mammalians muscle fibres³². Embers were seen in intact mouse fibres challenged with depolarization³³. By comparison, embers are more difficult to detect spontaneously in frog muscle fibres unless pharmacological manipulations were applied¹⁶.

A major difference lies in the spatial spreading of CREs from adult bees, which are clearly wider than CREs from adult intact vertebrate muscle fibres (FWHM, Tables 2–4). First, it could be related to the fact that sarcomeres are almost twice as long in bee leg muscle fibres 4.2 ± 0.1 µm⁶, as compared as vertebrate ones (~2–2.5 µm). Secondly, insect myofibrils are radially oriented around a central row of nuclei, (an arrangement which is not seen in adult vertebrate muscle fibres), giving a so-called ‘tubular’ morphology to the fibres. Accordingly, insect tubular muscle fibres have relatively straight t-tubules arranged in radial symmetry from the surface membrane to central nuclei in order to convey action potentials throughout the fibre volume³⁴. This spatial configuration may lead to less restricted diffusion along long portions of straight T-tubules. The average ratio of FWHMs (X/Y) honey bee CREs was greater than 1 (1.23 ± 0.01), which indicates that ellipsoidal sparks diffuse slightly more freely in the X direction (along the long axis of the fibre) than in the Y direction (along T-tubules). Elongated sparks were detected in frog skeletal muscle as well²²,³⁵ and CREs non-circularity has been previously attributed to a ‘channelled diffusion’ due to tubes of mitochondria¹⁷, but this hypothesis remains to be explored in bee fibres. Thirdly, an interesting feature is that bee T-tubules are running at the junction of A and I bands, like in mammalian skeletal muscle fibres, so that two T-tubules run per sarcomere⁵, increasing the chance to recruit release sites for large CREs. Fourthly, electron microscopy and rotary shadowing techniques showed that RyR disposition in invertebrate muscle differ from their disposition in vertebrates³⁶, an arrangement that may favour wider CREs in honey bees. Finally, ultrastructural differences were also identified in the invertebrate
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Parameters of calcium embers. Histogram of amplitudes (A), FWHMs (B), rise times (C), and durations (D). The number of embers was 1094.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Events on line-scan images at fast scan speed (0.2 ms/line). (A) Normalized line-scan images showing sparks (top) and an ember (bottom). Enlarged image of a spark (B) and an ember (C) with its spatial and time profile. The positions of the events are marked with yellow rectangular in panel A.</figcaption></figure>

coupions, i.e. the functional arrangement of RyRs and voltage-gated calcium channels from the plasma membrane (DHPRs, for dihydropyridine receptors), and previous studies underlined that this architecture was reminiscent of the cardiac striated muscle³⁷. Accordingly, other ultrastructural similarities can be pointed out in between bee skeletal and vertebrate cardiac muscles, such as the presence of dyads in both tissues instead of the triads observed in vertebrate skeletal muscles⁵. Besides structural and ultrastructural differences, that clearly settle the
<table><caption>**Table 2.** Interspecific comparison of calcium sparks' parameters measured in two dimensional (XY) experiments made in intact cells. In this table studies are ordered according to i) the developmental stage (from adult to foetal stages) and ii) the species (Honey bee, Frog and Mouse). All cells were freshly isolated muscle fibres. *Chun values are from adult P1, P7 and P14 diaphragm (D) and EDL muscle. Foetal diaphragm is from E18, i.e. 3 days before birth. #Personal communication. Data marked with ~ is approximated from figures. All data depicted at room temperature. FDB - flexor digitorum brevis; EDL - extensor digitorum longus; D-diaphragm; AM - acetoxymethyl ester; FWHM - full width at half maximum; EDHM - equivalent diameter at half maximum.</caption><thead><tr><th>Type of intact cells</th><th>Organism</th><th>Age (days after birth)</th><th>Sarcomere (μm)</th><th>Dye</th><th>Ext. Sol</th><th>[Ca<sup>2+</sup>]<sub>ext</sub> (mM)</th><th>Frame rate (image/s)</th><th>Incidence (% cells)</th><th>Frequency (events/s/mm²)</th><th>Amplitude (ΔF/F<sub>0</sub>)</th><th>FWHM or EDHM (μm)</th><th>Signal mass (μm³)</th><th>Reference</th></tr></thead><tbody><tr><td>Adult tibia fibres</td><td>Honey bee</td><td>1-4</td><td>~ 4</td><td>Fluo-8-AM</td><td>Tyrode</td><td>2</td><td>1</td><td>&lt; 20%</td><td>1696 ± 367</td><td>0.220 ± 0.001</td><td>3.73 ± 0.02 (X)<br>3.30 ± 0.02 (Y)<br>3.30 ± 0.02 (Y)</td><td>14.52 ± 0.29</td><td>This study</td></tr><tr><td>Adult FDB fibres</td><td>Frog</td><td>-</td><td>~ 2</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>30 and 60</td><td>~ 20%*</td><td>1470 ± 300</td><td>0.188 ± 0.001</td><td>1.73 ± 0.01 (x)<br>1.51 ± 0.01 (y)</td><td>1.02 ± 0.01</td><td>Cserne-Szappanos 2020</td></tr><tr><td>Adult FDB fibres</td><td>Mouse</td><td>&gt; 90</td><td>-</td><td>Fluo-3-AM</td><td>Krebs</td><td>2</td><td>1.2</td><td>&lt; 5%</td><td>-</td><td>-</td><td>-</td><td>-</td><td>Conklin 1999</td></tr><tr><td>Adult fibres*</td><td>Mouse</td><td>1</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>1</td><td>-</td><td>~ 80 (EDL)<br>~ 55 (D)</td><td>~ 0.9 (EDL)<br>~ 0.9 (D)</td><td>1.56 ± 0.04 (EDL)<br>~ 1.6 (D)</td><td>2.17 ± 0.20 (EDL)<br>~ 2.0 (D)</td><td>Chun 2003*</td></tr><tr><td>Adult fibres*</td><td>Mouse</td><td>7</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>1</td><td>-</td><td>~ 5 (EDL)<br>~ 25 (D)</td><td>~ 0.9 (EDL)<br>~ 0.7 (D)</td><td>~ 1.25 (EDL)<br>~ 1.32 (D)</td><td>~ 1.2 (EDL)<br>~ 0.8 (D)</td><td>Chun 2003*</td></tr><tr><td>Adult fibres*</td><td>Mouse</td><td>14</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>1</td><td>-</td><td>~ 7 (EDL)<br>~ 5 (D)</td><td>~ 1 (EDL)<br>0.65 ± 0.05 (D)</td><td>1.10 ± 0.05 (EDL)<br>~ 1.3 (D)</td><td>0.66 ± 0.12 (EDL)<br>0.80 ± 0.25 (D)</td><td>Chun 2003*</td></tr><tr><td>Adult interossei fibres</td><td>Mouse</td><td>-</td><td>-</td><td>Fluo-4AM</td><td>Ringer</td><td>2</td><td>0.9</td><td>-</td><td>15.1 ± 1.2</td><td>-</td><td>-</td><td>-</td><td>Teichmann 2008</td></tr><tr><td>Foetal intercostal fibres</td><td>Mouse</td><td>3 days before birth</td><td>-</td><td>Fluo-3-AM</td><td>Krebs</td><td>2</td><td>1.2</td><td>&lt; 20%</td><td>-</td><td>-</td><td>-</td><td>-</td><td>Conklin 1999</td></tr><tr><td>Foetal diaphragm fibres</td><td>Mouse</td><td>3 days before birth</td><td>-</td><td>Fluo-4AM</td><td>Ringer</td><td>1.8</td><td>1</td><td>~ 20%</td><td>~ 80</td><td>1.24 ± 0.02</td><td>~ 1.67 (D)</td><td>3.79 ± 0.22 (D)</td><td>Chun 2003*</td></tr></tbody></table>

conditions for broad CREs in insect muscle fibres, genomic differences could favour it as well. Homology data on RyR sequences suggest that invertebrates and vertebrates show important genomic differences. First, insects have only one *ryr* gene whereas amphibians have two (*ryr* α also named *ryr* 1, and *ryr* β also named *ryr* 3) and mammalians have three<sup>38,39</sup>. Two isoforms are actually expressed in vertebrate skeletal muscles (*ryr* 1 and *ryr* 3 in mammalians and amphibians), since the *ryr* 2 gene is expressed in cardiac mammalian muscle only. After the 'honey bee genome sequencing consortium' made available the honey bee (*Apis mellifera*) genome sequence<sup>40</sup>, the predicted amino acid sequence for the *Apis mellifera* RyR protein (*AmRyR*) has been obtained from an annotated genomic sequence NC_037639.1 (Gene 408680 on chromosome LG2) using traditional gene prediction methods (GNOMON, supported by EST evidence). Like in other species, predicted *AmRyR* protein isoforms (X1 to X9) contain ~ 5000 amino acids. According to BLAST tools (NCBI database), the homology degree of the honey bee *AmRyR* with vertebrate RyRs is around 45%. Such a difference may well settle the basis for strong differences into the protein functions and regulations leading to morphological differences in bee CREs as compared to other species. In zebrafish skeletal muscle cells, two RyR isoforms have been shown to be essential to achieve normal sparks production<sup>41</sup>. The second frog RyR isoform may also be related to high spark production in this species<sup>42</sup>. In foetal mammalian muscles, the two isoforms of RyR (1 and 3) are expressed and it was suggested that they both contribute to generate sparks<sup>43</sup>. Conversely, at the mature adult stage, most intact mammalian skeletal fibres (which no longer express RyR3) require specific conditions to fire sparks and at a relatively low frequency. Our data show that a single *ryr* gene allows the production of a high CREs frequency with wide spatial spread in intact insect fibres.

Functional differences could also lead to the specific shape properties of bees CREs. In the honeybee muscle, as in other invertebrates, excitation-contraction coupling relies on i) Ca²⁺ entry through the plasmalemmal membrane via voltage-gated channels and ii) on a calcium induced calcium release (CICR) mechanism<sup>5,44</sup>. CICR is thus a crucial mechanism in the skeletal muscle from bees and other insects (which are devoid of TTX-sensitive voltage-gated channels) where voltage-gated calcium currents activate quickly to achieve a rapid action potential overshoot<sup>6,45,46</sup>. The role of CICR in excitation-contraction coupling implies that honey bee RyRs must be highly sensitive to Ca²⁺, a feature that favours broad CREs. CICR has been shown in frog muscle fibres as well
<table><caption>**Table 3.** Interspecific comparison of calcium sparks' and depolarization-evoked transients' parameters measured in line-scan (XT) experiments performed in intact cells. In this table studies are ordered according to i) the developmental stage (from adult to foetal/larval stages) and ii) the species (Honey bee, Frog, Mouse, Rat, Zebrafish). All preparations were freshly isolated muscle fibres except the last four where the cells were in cell culture. *Chun values are from P1, P7, and P14 EDL muscle and from foetal (3 days before birth) diaphragm, rise time was calculated between 10-90% of the peak amplitude. †sarcomere stretched to ~3 µm. §parameters of CREs originate from sparks and embers together. Data marked with ~ is approximated from figures. All data depicted at room temperature except Perni et al. 2015 (18 °C). Hollingworth et al. 2001 cooled the fibres to 17-20 °C. All studies were from intact cells bathed in a physiological solution (Ringer, Krebs or Tyrode), except for Perni et al. 2015 (caffeine added) and Shirokova et al. 1998 (150 TEA-MeSO<sub>3</sub>)<sup>57</sup>. FDB - flexor digitorum brevis; EDL - extensor digitorum longus; D- diaphragm; AM - acetoxyethyl ester; FWHM - full width at half maximum; FDHM - full duration at half maximum; FTHM - full time at half maximum.</caption><thead><tr><th>Type of intact cells</th><th>Organism</th><th>Age (days)</th><th>Sarcomere (µm)</th><th>Dye</th><th>Ext. Sol</th><th>(Ca<sup>2+</sup>)<sub>int</sub> (mM)</th><th>Speed of recording (line/ms)</th><th>Frequency (events/s/mm)</th><th>Sparks amplitude (ΔF/F<sub>0</sub>)</th><th>FWHM (µm)</th><th>FDHM or FTHM (ms)</th><th>Rise time (ms)</th><th>Signal mass (µm<sup>2</sup>)</th><th>Depolarization-induced macroscopic Ca<sup>2+</sup> transient (ΔF/F<sub>0</sub>)</th><th>Reference</th></tr></thead><tbody><tr><td>Adult tibia fibres</td><td>Honey bee</td><td>1-4</td><td>-4</td><td>Fluo-8-AM</td><td>Tyrode</td><td>2</td><td>1</td><td>174.0 ± 14.4</td><td>0.185 ± 0.001</td><td>2.65 ± 0.02</td><td>25.52 ± 0.16</td><td>10.45 ± 0.11</td><td>6.00 ± 0.20</td><td>1.5 ± 0.16</td><td>This study</td></tr><tr><td>Adult iliofibularis and semitendinosus fibres</td><td>Frog</td><td>-</td><td>-3†</td><td>Fluo-3 (injected)</td><td>Ringer</td><td>1.8</td><td>2.048</td><td>23 ± 9</td><td>0.92 ± 0.02</td><td>0.95 ± 0.04</td><td>5.66 ± 0.23</td><td>0-100% rise time 3.49 ± 0.15</td><td>1.28 ± 0.17</td><td>18.6 ± 0.9</td><td>Hollingworth 2001</td></tr><tr><td>Adult FDB fibres</td><td>Mouse</td><td>> 90</td><td>-</td><td>Fluo-3-AM</td><td>Krebs</td><td>2</td><td>2.05</td><td>&lt;20</td><td>0.9 ± 0.3</td><td>2.0 ± 0.5</td><td>40.3 ± 35</td><td>16.4 ± 14</td><td>-</td><td>-</td><td>Conklin 1999</td></tr><tr><td>Adult fibres</td><td>Mouse</td><td>1</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>2</td><td>-</td><td>0.7 ± 0.1 (EDL)<br>0.8 ± 0.1 (D)</td><td>1.5 ± 0.1 (EDL)<br>0.7 ± 0.2 (D)</td><td>13.2 ± 0.7 (EDL)<br>10.5 ± 1.1 (D)</td><td>8.9 ± 0.7 (EDL)<br>8.0 ± 1.2 (D)</td><td>-</td><td>-</td><td>Chun 2003*</td></tr><tr><td>Adult fibres</td><td>Mouse</td><td>7</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>2</td><td>-</td><td>0.6 ± 0.1 (EDL)<br>0.7 ± 0.1 (D)</td><td>1.7 ± 0.1 (EDL)<br>1.6 ± 0.2 (D)</td><td>12.6 ± 0.9 (EDL)<br>11.3 ± 1.1 (D)</td><td>7.3 ± 0.7 (EDL)<br>6.1 ± 0.6 (D)</td><td>-</td><td>-</td><td>Chun 2003*</td></tr><tr><td>Adult fibres</td><td>Mouse</td><td>14</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>2</td><td>-</td><td>1.5 ± 0.5 (EDL)<br>0.5 ± 0.1 (D)</td><td>1.6 ± 0.2 (EDL)<br>1.5 ± 0.1 (D)</td><td>13.5 ± 1.6 (EDL)<br>11.2 ± 0.9 (D)</td><td>9.3 ± 1.4 (EDL)<br>7.6 ± 1.3 (D)</td><td>-</td><td>-</td><td>Chun 2003*</td></tr><tr><td>Foetal fibres</td><td>Mouse</td><td>3 days before birth</td><td>-</td><td>Fluo-3-AM</td><td>Krebs</td><td>2</td><td>2.05</td><td>&lt;20</td><td>1.6 ± 0.6</td><td>3.5 ± 1.1</td><td>84.5 ± 46</td><td>42.8 ± 24</td><td>-</td><td>-</td><td>Conklin 1999</td></tr><tr><td>Foetal fibres</td><td>Mouse</td><td>3 days before birth</td><td>-</td><td>Fluo-4-AM</td><td>Ringer</td><td>1.8</td><td>2</td><td>-</td><td>0.9 ± 0.1</td><td>1.9 ± 0.1</td><td>13.9 ± 0.4</td><td>Rise(10-90%)<br>7.6 ± 0.3</td><td>-</td><td>-</td><td>Chun 2003*</td></tr><tr><td>Larval tail fibres</td><td>Zebrafish</td><td>3 days post-fertilization</td><td>-1.5</td><td>Fluo-4-AM</td><td>Ringer with 0.3 mM Caffeine</td><td>2</td><td>-</td><td>-44</td><td>-0.8</td><td>-1.3</td><td>-5.2</td><td>-4</td><td>-3</td><td>-</td><td>Perni 2015</td></tr><tr><td>Primary myotubes culture</td><td>Mouse</td><td>-</td><td>-</td><td>Fluo-3-AM</td><td>TEA-MeSO<sub>3</sub></td><td>2</td><td>2</td><td>-</td><td>0.94 ± 0.05</td><td>2.1 ± 0.3</td><td>11.8 ± 1.6</td><td>-</td><td>-</td><td>-</td><td>Shirokova 1998</td></tr><tr><td>Primary myotubes culture</td><td>Rat</td><td>-5 days in vitro</td><td>-</td><td>Fluo-4-AM</td><td>Tyrode</td><td>1.8</td><td>1.54</td><td>3.2 ± 1.9</td><td>0.49 ± 0.03</td><td>1.67 ± 0.16</td><td>39.7 ± 4.3</td><td>15.4 ± 2.2</td><td>-</td><td>0.62 ± 0.05</td><td>Fodor 2008</td></tr><tr><td>C2C12 cell line myotubes</td><td>Mouse</td><td>3-6 days in vitro</td><td>-</td><td>Fluo-3-AM</td><td>Ringer</td><td>10</td><td>0.02-007</td><td>0.3</td><td>-300 nM</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>Gyorke 1996</td></tr><tr><td>C2C12 cell line myotubes</td><td>Mouse</td><td>> 5 days in vitro</td><td>-</td><td>Fluo-4-AM</td><td>Tyrode</td><td>2</td><td>0.8</td><td>16</td><td>0.49 ± 0.01</td><td>3.8 ± 0.08</td><td>92.8 ± 8.7</td><td>55.9 ± 4.2</td><td>-</td><td>-2.5</td><td>Altomare 2013<sup>5</sup></td></tr></tbody></table>

(with a crucial role attributed to the second RyR isoform) but it plays a null or minor role in rat fibres<sup>47</sup>. Despite a demonstrated CICR in frog, CREs from frog fibres are narrower than CREs from mammalian fibres, both in intact cells (Table 3) or in permeabilized fibres<sup>48</sup>. In the two vertebrate species, CREs width is narrower than in bees, which is clearly in agreement with a more important role for CICR in insects. In half of adult bee leg muscle fibres, a low-voltage-activated (LVA) current has been identified in addition to the nifedipine-sensitive high-voltage activated (HVA) current<sup>5,49</sup>. Both HVA- and LVA- calcium currents exists in mammalian fibres as well, but the density of the later gradually decreases postnatally to eventually disappear after 2-3 weeks of age<sup>50</sup>. Whereas DHPRs responsible for the HVA- (L-type) current are involved in triggering and/or repressing sparks in immature and adult mammalian muscle fibres<sup>27,32</sup>, pharmacological proofs indicate that the LVA- (T-type) calcium channel doesn't modulate sparks frequency<sup>32</sup>. A role for the LVA- current in the genesis or regulation of honey bee CREs thus remains to be studied.

Amplitude of honey bee sparks measured in *XYT* experiments appeared smaller than in almost all other intact preparations (Table 2, 3). This parameter is however highly subjected to the experimental set-up performances, imaging analysis and thresholding strategies, as illustrated in *XT* experiments from frog, where amplitude can vary by a factor 5 (Table 3). Namely, using the traditional threshold-based sparks detection method, the former study of sparks in intact frog fibres described an average amplitude of ~1 ΔF/F<sub>0</sub><sup>23</sup>, whereas a recent publication using the wavelet detection method described an average amplitude of ~0.2 ΔF/F<sub>0</sub><sup>22</sup>. The difference may be in part attributed to the more efficient detection of CREs of low amplitude with the wavelet method, but other factors may be involved because the maximal amplitudes detected in the later publication reached ~0.85 ΔF/F<sub>0</sub><sup>22</sup>. Interestingly, in intact muscle fibres from a mouse model of myopathy (*mdx*), a similar analysis strategy (à trous implementation of the two-dimensional undecimated discrete wavelet transform') led to the detection of CREs up to 1.3 ΔF/F<sub>0</sub>, with a half cumulative frequency value (CF<sub>50</sub>) of 0.53<sup>51</sup>. The apparently low amplitude of honey
<table><caption>**Table 4.** Interspecific comparison of calcium embers' parameters measured in line-scan (XT) experiments performed in intact cells. In this table studies are ordered according to i) the developmental stage and ii) the species (Honey bee, Frog, Mouse, Rat). Preparations were freshly isolated muscle fibres or cell culture. *Chun values are from EDL and diaphragm muscles. Rise time was calculated between 10–90% of the peak amplitude. †The embers frequency was given as a percentage of total CRE in these studies. Data marked with ~ is approximated from figures. All data depicted at room temperature. FDB – flexor digitorum brevis; AM – acetoxymethyl ester; FWHM – full width at half maximum.</caption><thead><tr><th>Type of intact cells</th><th>Organism</th><th>Age (days)</th><th>Sarcomere (μm)</th><th>Dye</th><th>Ext. Sol</th><th>[Ca<sup>2+</sup>]<sub>ext</sub> (mM)</th><th>Speed of recording (line/ms)</th><th>Threshold for duration of embers (ms)</th><th>Frequency (events/ s/ mm)</th><th>Amplitude (ΔF/F<sub>0</sub>)</th><th>FWHM (μm)</th><th>Duration (ms)</th><th>Rise time (ms)</th><th>Signal mass (μm<sup>2</sup>)</th><th>Reference</th></tr></thead><tbody><tr><td>Adult tibia fibres</td><td>Honey bee</td><td>1–4</td><td>~ 4</td><td>Fluo-8 AM</td><td>Tyrode</td><td>2</td><td>1</td><td>40</td><td>63.1 ± 6.6</td><td>0.186 ± 0.002</td><td>3.18 ± 0.04</td><td>56.95 ± 0.58</td><td>27.02 ± 0.41</td><td>10.29 ± 0.41</td><td>This study</td></tr><tr><td>Adult FDB fibres</td><td>Mouse</td><td>1–14</td><td>–</td><td>Fluo-4 AM</td><td>Ringer</td><td>1.8</td><td>2</td><td>30</td><td>~ 4% of CREs †</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td><td>Chun 2003*</td></tr><tr><td>Primary myotubes culture</td><td>Rat</td><td>~ 5 days in vitro</td><td>–</td><td>Fluo-4 AM</td><td>Tyrode</td><td>1.8</td><td>1.54</td><td>60</td><td>2.4 ± 1.4</td><td>0.27 ± 0.05</td><td>–</td><td>229 ± 33</td><td>–</td><td>–</td><td>Fodor 2008</td></tr><tr><td>C2C12 cell line myotubes</td><td>Mouse</td><td>> 5 days in vitro</td><td>–</td><td>Fluo-4 AM</td><td>Tyrode</td><td>2</td><td>0.8</td><td>50</td><td>53.4% of CREs †</td><td>–</td><td>–</td><td>–</td><td>–</td><td>–</td><td>Altomare 2013</td></tr></tbody></table>

bee sparks as compared to canonical sparks may thus be more related to species specificities and/or experimental conditions than to the detection method we used. Another striking difference between bees and vertebrates is that in bees, embers had almost the same amplitude as sparks, whereas embers are usually smaller than sparks in mammals, as shown in non-intact preparations¹⁶,³¹. This finding implies the idea that in bee muscle during a CRE, the same number of RyRs open and only the duration of their opening vary stochastically. The amplitude of sparks from honey bee is ~ eightfold smaller than the amplitude of voltage-gated global $Ca^{2+}$ transients (Table 3), a result consistent with the idea that sparks are the unitary building blocks of macroscopic calcium release. This ratio is in the same range in intact frog fibres (in *Rana pipiens* fibres, sparks are ~ 20 smaller than voltage-activated macroscopic transients²³). In intact mouse primary myotubes and in the C2C12 mouse cell line, sparks are 1.4 and ~ fivefold smaller than macroscopic transients, respectively⁵²,⁵³. To the best of our knowledge, no study relates both sparks amplitude and macroscopic transients in intact non voltage-clamped mouse fibres. A ratio of ~ 6 can however be calculated from macroscopic $Ca^{2+}$ transients recorded under similar conditions as those depicted in Table 3 for fibres of the mouse FDB muscle⁵⁵.

Honey bee sparks had similar kinetics properties as other species. Most features of honey bee embers were reminiscent to other species ones as well. Previous studies classified the embers using their durations (Table 4). Although there was no clear border between the duration of sparks and embers in bee muscle (see Fig. 4D), the distribution of the probability histogram of bee CREs duration was nevertheless successfully fitted by a bimodal Gaussian function, according to a previously described methodology⁵⁵ and the value 40 ms was retained as the threshold for ember detection. Importantly, the amplitude could not be used as the second separator parameter since CREs in bee fibres are more homogeneous in amplitude than in mammalian muscles where the amplitude of embers is much smaller than that of the sparks.

To our knowledge, few studies depict incidence of sparks in intact skeletal muscle (i.e. the percentage of fibres where such events can be detected). CREs were seen in < 20% of adult honey bee intact fibres. By comparison, the incidence of spontaneous CREs was less than 5% of intact adult mouse fibres whereas it was more elevated in embryonic fibres⁴³. A comparison of sparks frequency collected in XT line-scan experiments in intact cells shows that in honey bee, frequency is four fold to 600 fold higher than in any other mature or immature cells from frog, mouse, rat or zebrafish (Table 3, XT). In *XYT* images series from adult bees' fibres, the frequency is one or two orders of magnitude higher as compared to intact cells from adult or immature mammalian cells (20–300 fold higher, Table 2). For instance, in mouse cells the frequency of CREs in *XYT* has been described to be between 5 and 80 Hz/mm² only (Table 2, *XYT*). Interestingly, the *XYT* frequency observed in bees is almost identical to what was reported in intact adult frog muscle recently²². Again, these differences in CREs frequency among species may thus be a consequence of the relative incidence that CICR has in excitation-contraction coupling or to ultrastructural differences discussed before. It is important to point out that frequency similarities between bees and frogs are not only related to the analysis strategy, because a similar procedure used in mouse fibres led to a 100 fold lower frequency⁵¹. Finally, the proportion of embers as compared to the total number of CREs in *XT* experiments shows that honey bee embers are much more frequent than in adult or foetal mouse cells (33% in honey bee vs 4% in mouse) and either higher or identical in myotubes from rat and mouse respectively (Table 4).

In honey bee fibres, CREs appear most of the time in a stochastic manner and rarely at a given triad. However, a repetitive mode of activation was observed in specific hot spots. This second mode, previously termed *rep-mode*, has been identified when voltage-clamped vertebrate fibres (cut fibres mounted in a Vaseline gap chamber) were stimulated by a low-amplitude depolarization⁵⁶. This activity mode can also happen in chronically depolarized frog fibre right after a brief repolarization protocol. Such an activity was seen in intact adult mouse fibres as well⁴³ and in embryonic and postnatal mouse fibres, where 29% of all events location had reoccurring $Ca^{2+}$ release events³². A pharmacological manipulation with caffeine is also increasing the appearance of the repetitive CREs in intact frog fibres²². However, whereas CREs in *rep-mode* were short in duration in frog fibres, bee hot spots led to CREs with various kinetics ranging from sparks to burst and embers.
In addition to provide the first interspecific comparison of sparks between vertebrates and invertebrates, the present study sets the basis for a future comparison of sparks in other invertebrates, where these subcellular localized Ca²⁺ signals have not been described yet. We have overcome the difficulty of isolating single fibres from insects (firmly attached to the insect exoskeleton through tonofilaments), and this methodological advancement will allow an in-depth characterization of sparks in insects in general. Our methods could in particular be transposed to common insect models used either in fundamental research (fruit flies, cockroaches, grasshoppers) or models of interest in agronomy (aphids, moths, bugs, etc.). This study would be of particular interest to characterize the interspecific difference in susceptibility to new insecticides targeting RyRs, the phthalic acid diamides and anthranilic diamides. Diamides are highly efficient to control a variety of pests. They are widely used to kill a number of insect species, including lepidopterans (moths), hemipterans (aphids), dipterans (flies) and so on. Their major mode of action is to interfere with normal function of ryanodine receptors. They have been designed to ensure a favourable vertebrate selectivity ratio, that is to say: being more toxic to invertebrates than to vertebrates. However, insects exposed to these compounds, whether they are considered either as pests or beneficial insects, show strong symptoms of intoxication (from sublethal effects to mortality). For instance, we have shown that sublethal doses of chlorantraniliprole were highly deleterious to honey bees that encounter long lasting locomotor deficits lasting for several days, leading to high susceptibility to other threats⁹. Higher doses kill the bees quickly. Insect ryanodine receptors are quite divergent from those of vertebrates (~45% a.a. identity), but RyRs from insect species from all orders show a strong molecular proximity (~ 80%). The honeybee RyR (XP_006569098.1) shows 83, 80, 79 and 78% amino acid identity with RyRs from the beetle *Tribolium castaneum* (NP_001308588.1), the moth *Plutella xylostella* (XP_037970554.1), the aphid *Myzus Persicae* (XP_022160123.1) and the fly *Drosophila melanogaster* (NP_476991.1), respectively. Although fruit flies *D. melanogaster* are considered as models, the intracellular homeostasis in insect tubular muscle fibres has so far not been studied thoroughly owing to the difficulty of isolating muscle fibres from the cuticle, which is a prerequisite for calcium imaging. The honey bee fibres provide a suitable model in this context and the dissociation procedure we developed here to study sparks will be helpful in many instances.

## Conclusion

In addition to the characterization of localized Ca²⁺-release events of the honey bee, a detailed comparison of their major specificities as compared with their vertebrate counterparts has been carried out. The successful isolation of insect muscle fibres and the thorough characterization of subcellular Ca²⁺ signalling in the honey bee muscle paves the way for further in vitro studies, which are clearly needed to elucidate the origins of the established differential toxicity of the RyR-targeting diamides insecticides between vertebrates and invertebrates. Furthermore, while a so-called 'vertebrate selectivity toxicological ratio' already exists for diamides, no protective selectivity ratio has been established yet for useful insects. Pest insects and non-pest insects share a high degree of homology in their RyR and a comparative study of Ca²⁺ sparks among insects may help to better characterize the toxicity of diamides, which constitute a threat to useful insects such as pollinators.

## Data availability

The datasets generated and analysed during the current study are available on reasonable request.

Received: 19 March 2021; Accepted: 4 August 2021

Published online: 18 August 2021

## References

1. Galizia, G., Eisenhardt, D. & Giurfa, M. *Honeybee Neurobiology and Behavior* (Springer, 2012).

2. Winston, M. L. *The Biology of the Honey Bee*. (Harvard University Press, 1991).

3. Bicker, G. Transmitter-induced calcium signalling in cultured neurons of the insect brain. *J. Neurosci. Methods* **69**, 33–41 (1996).

4. Szyszka, P., Galkin, A. & Menzel, R. Associative and non-associative plasticity in kenyon cells of the honeybee mushroom body. *Frontiers in Systems Neuroscience* **2**. https://doi.org/10.3389/neuro.06.003.2008 (2008).

5. Collet, C. Excitation-contraction coupling in skeletal muscle fibers from adult domestic honeybee. *Pflugers Arch* **458**, 601–612. https://doi.org/10.1007/s00424-009-0642-6 (2009).

6. Collet, C. & Belzunces, L. Excitable properties of adult skeletal muscle fibres from the honeybee Apis mellifera. *J. Exp. Biol.* **210**, 454–464. https://doi.org/10.1242/jeb.02667 (2007).

7. Deisig, N., Giurfa, M., Lachnit, H. & Sandoz, J. C. Neural representation of olfactory mixtures in the honeybee antennal lobe. *Eur. J. Neurosci.* **24**, 1161–1174 (2006).

8. Galizia, C. G., Kuttner, A., Joerges, J. & Menzel, R. Odour representation in honeybee olfactory glomeruli shows slow temporal dynamics: an optical recording study using a voltage-sensitive dye. *J. Insect Physiol.* **46**, 877–886 (2000).

9. Kadala, A., Charreton, M., Charnet, P. & Collet, C. Honey bees long-lasting locomotor deficits after exposure to the diamide chlorantraniliprole are accompanied by brain and muscular calcium channels alterations. *Sci. Rep.* **9**, 2153. https://doi.org/10.1038/s41598-019-39193-3 (2019).

10. Kadala, A., Charreton, M. & Collet, C. Flubendiamide, the first phthalic acid diamide insecticide, impairs neuronal calcium signaling in the honey bee's antennae. *J. Insect. Physiol.* **104**.086. https://doi.org/10.1016/j.jinsphys.2020.104086 (2020).

11. Walz, B., Baumann, O., Zimmermann, B. & Ciriacy-Wantrup, E. V. Caffeine- and ryanodine-sensitive Ca²⁺-induced Ca²⁺ release from the endoplasmic reticulum in honeybee photoreceptors. *J. Gen. Physiol.* **105**, 537–567 (1995).

12. Meissner, G. & Lu, X. Dihydropyridine receptor-ryanodine receptor interactions in skeletal muscle excitation-contraction coupling. *Biosci. Rep.* **15**, 399–408. https://doi.org/10.1007/bf01788371 (1995).

13. Collet, C. & Jacquemond, V. Sustained release of calcium elicited by membrane depolarization in ryanodine-injected mouse skeletal muscle fibers. *Biophys. J.* **82**, 1509–1523. https://doi.org/10.1016/S0006-3495(02)75504-5 (2002).

14. Carafoli, E. & Krebs, J. Why calcium? How calcium became the best communicator. *J. Biol. Chem.* **291**, 20849–20857. https://doi.org/10.1074/jbc.R116.735894 (2016).

15. Collet, C. in *Toxicologie chez l'abeille domestique* (eds C. Collet, Sandoz J.C., & P. Charnet) (CNRS Editions, in press, 2021).
16. Csernoch, L. Sparks and embers of skeletal muscle: the exciting events of contractile activation. *Pflügers Arch. Eur. J. Physiol.* **454**, 869–878. https://doi.org/10.1007/s00424-007-0244-0 (2007).

17. Cheng, H. & Lederer, W. J. Calcium sparks. *Physiol. Rev.* **88**, 1491–1545. https://doi.org/10.1152/physrev.00030.2007 (2008).

18. Cheng, H., Lederer, W. & Cannell, M. Calcium sparks: elementary events underlying excitation-contraction coupling in heart muscle. *Science* **262**, 740–744. https://doi.org/10.1126/science.8235594 (1993).

19. Niggli, E. & Shirokova, N. A guide to sparkology: the taxonomy of elementary cellular Ca2+ signaling events. *Cell Calcium* **42**, 379–387. https://doi.org/10.1016/j.ceca.2007.02.010 (2007).

20. Collet, C. et al. Subcellular calcium events and calcium waves in leg skeletal muscle fibers isolated from the honey bee Apis mel-
lifera. *Biophys. J.* **114**, 289a. https://doi.org/10.1016/j.bpj.2017.11.1655 (2018).

21. Szabó, L. Z., Vincze, J., Csernoch, L. & Szentesi, P. Improved spark and ember detection using stationary wavelet transforms. *J. Theor. Biol.* **264**, 1279–1292. https://doi.org/10.1016/j.jtbi.2010.04.005 (2010).

22. Cserne Szappanos, H. et al. High Time Resolution Analysis of Voltage-Dependent and Voltage-Independent Calcium Sparks in Frog Skeletal Muscle Fibers. *Frontiers in Physiology* **11**, https://doi.org/10.3389/fphys.2020.599822 (2020).

23. Hollingworth, S., Peet, J., Chandler, W. K. & Baylor, S. M. Calcium sparks in intact skeletal muscle fibers of the frog. *J. Gen. Physiol.* **118**, 653–678. https://doi.org/10.1085/jgp.118.6.653 (2001).

24. Wareham, A. C., Duncan, C. J. & Bowler, K. The resting potential of cockroach muscle membrane. *Comp. Biochem. Physiol. A Physiol.* **48**, 765–797. https://doi.org/10.1016/0300-9629(74)90619-7 (1974).

25. Bekoff, A. & Betz, W. Properties of isolated adult rat muscle fibres maintained in tissue culture. *J. Physiol.* **271**, 537–547. https://
doi.org/10.1113/jphysiol.1977.sp012013 (1977).

26. Allard, B., Bernengo, J. C., Rougier, O. & Jacquemond, V. Intracellular Ca2+ changes and Ca2+-activated K+ channel activation induced by acetylcholine at the endplate of mouse skeletal muscle fibres. *J. Physiol.* **494**, 337–349. https://doi.org/10.1113/jphysiol.1996.sp021496 (1996).

27. Pickering, J. D., White, E., Duke, A. M. & Steele, D. S. DHPR activation underlies SR Ca2+ release induced by osmotic stress in isolated rat skeletal muscle fibers. *J. Gen. Physiol.* **133**, 511–524. https://doi.org/10.1085/jgp.200910191 (2009).

28. Baumann, O. Distribution of ryanodine receptor Ca2+ channels in insect photoreceptor cells. *J. Compar. Neurol.* **421**, 347–361. https://doi.org/10.1002/(SICI)1096-9861(20000605)421:3%3c347::AID-CNE4%3e3.0.CO;2-Z (2000).

29. Parker, I. & Gil Wier, W. Variability in frequency and characteristics of Ca2+ sparks at different release sites in rat ventricular myocytes. *J. Physiol.* **505**, 337–344. https://doi.org/10.1111/j.1469-7793.1997.337bb.x (1997).

30. González, A. et al. The spark and its ember: separately gated local components of Ca2+ release in skeletal muscle. *J. Gen. Physiol.* **115**, 139–158. https://doi.org/10.1085/jgp.115.2.139 (2000).

31. Kirsch, W. G., Uttenweiler, D. & Fink, R. H. Spark- and ember-like elementary Ca2+ release events in skinned fibres of adult mammalian skeletal muscle. *J. Physiol.* **537**, 379–389. https://doi.org/10.1111/j.1469-7793.2001.00379.x (2001).

32. Chun, L. G., Ward, C. W. & Schneider, M. F. Ca2+ sparks are initiated by Ca2+ entry in embryonic mouse skeletal muscle and decrease in frequency postnatally. *Am. J. Physiol. Cell Physiol.* **285**, C686–697. https://doi.org/10.1152/apjcell.00072.2003 (2003).

33. Csernoch, L., Pouvreau, S., Ronjat, M. & Jacquemond, V. Voltage-activated elementary calcium release events in isolated mouse skeletal muscle fibers. *J. Membr. Biol.* **226**, 43. https://doi.org/10.1007/s00232-008-9138-0 (2008).

34. Gilai, A. Electromechanical coupling in tubular muscle fibers. II. Resistance and capacitance of one transverse tubule. *J. Gen. Physiol.* **67**, 343–367. https://doi.org/10.1085/jgp.67.3.343 (1976).

35. Brum, G., Gonzalez, A., Rengifo, J., Shirokova, N. & Rios, E. Fast imaging in two dimensions resolves extensive sources of Ca2+ sparks in frog skeletal muscle. *J. Physiol.* **528**, 419–433. https://doi.org/10.1111/j.1469-7793.2000.00419.x (2000).

36. Loesser, K. E., Castellani, L. & Franzini-Armstrong, C. Dispositions of junctional feet in muscles of invertebrates. *J. Muscle Res. Cell Motil.* **13**, 161–173. https://doi.org/10.1007/BF01874153 (1992).

37. Takekura, H. & Franzini-Armstrong, C. The structure of Ca2+ release units in arthropod body muscle indicates an indirect mechanism for excitation-contraction coupling. *Biophys. J.* **83**, 2742–2753. https://doi.org/10.1016/S0006-3495(02)75284-3 (2002).

38. Xu, X., Bhat, M. B., Nishi, M., Takeshima, H. & Ma, J. Molecular Cloning of cDNA encoding a drosophila ryanodine receptor and functional studies of the carboxyl-terminal calcium release channel. *Biophys. J.* **78**, 1270–1281. https://doi.org/10.1016/S0006-3495(00)76683-5 (2000).

39. Troczka, B. J., Richardson, E., Homem, R. A. & Davies, T. G. E. An analysis of variability in genome organisation of intracellular calcium release channels across insect orders. *Gene* **670**, 70–86. https://doi.org/10.1016/j.gene.2018.05.075 (2018).

40. Weinstock, G. M., Robinson, G. E. & Consortium, T. H. G. S. Insights into social insects from the genome of the honeybee Apis mellifera. *Nature* **443**, 931–949. https://doi.org/10.1038/nature05260 (2006).

41. Perni, S. et al. Structural and functional properties of ryanodine receptor type 3 in zebrafish tail muscle. *J. Gen. Physiol.* **145**, 253. https://doi.org/10.1085/jgp.20141130302112015c (2015).

42. Felder, E. & Franzini-Armstrong, C. Type 3 ryanodine receptors of skeletal muscle are segregated in a parajunctional position. *Proc. Natl. Acad. Sci. USA* **99**, 1695–1700. https://doi.org/10.1073/pnas.032657599 (2002).

43. Conklin, M. W., Barone, V., Sorrentino, V. & Coronado, R. Contribution of ryanodine receptor type 3 to Ca(2+) sparks in embryonic mouse skeletal muscle. *Biophys J.* **77**, 1394–1403. https://doi.org/10.1016/S0006-3495(99)76988-2 (1999).

44. Gyorke, S. & Palade, P. Calcium-induced calcium release in crayfish skeletal muscle. *J. Physiol.* **457**, 195–210. https://doi.org/10.1113/jphysiol.1992.sp019373 (1992).

45. Ashcroft, F. M. & Stanfield, P. R. Calcium and potassium currents in muscle fibres of an insect (*Carausius morosus*). *J. Physiol.* **323**, 93–115 (1982).

46. Gielow, M. L., Gu, G. G. & Singh, S. Resolution and pharmacological analysis of the voltage-dependent calcium channels of Dros-
ophila larval muscles. *J. Neurosci.* **15**, 6085–6093 (1995).

47. Rios, E. & Zhou, J. Control of dual isoforms of Ca2+ release channels in muscle. *Biol. Res.* **37**, 583–591 (2004).

48. Zhou, J. et al. Ca2+ Sparks and embers of mammalian muscle. Properties of the sources. *J. Gen. Physiol.* **122**, 95–114. https://doi.
org/10.1085/jgp.200308796 (2003).

49. Rousset, M. et al. Honeybee locomotion is impaired by Am-CaV3 low voltage-activated Ca2+ channel antagonist. *Sci. Rep.* **7**, 41782. https://doi.org/10.1038/srep41782 (2017).

50. Beam, K. G. & Knudson, C. M. Effect of postnatal development on calcium currents and slow charge movement in mammalian skeletal muscle. *J. Gen. Physiol.* **91**, 799–815. https://doi.org/10.1085/jgp.91.6.799 (1988).

51. Teichmann, M. D. H. et al. Inhibitory control over Ca2+ Sparks via mechanosensitive channels is disrupted in dystrophin deficient muscle but restored by mini-dystrophin expression. *PLoS ONE* **3**, e3644. https://doi.org/10.1371/journal.pone.0003644 (2008).

52. Fodor, J. et al. Altered expression of triadin 95 causes parallel changes in localized Ca2+ release events and global Ca2+ signals in skeletal muscle cells in culture. *J. Physiol.* **586**, 5803–5818. https://doi.org/10.1113/jphysiol.2008.160457 (2008).

53. Altomare, C. et al. Altered functional differentiation of mesoangioblasts in a genetic myopathy. *J. Cell Mol. Med.* **17**, 419–428. https://doi.org/10.1111/jcmm.12023 (2013).

54. Caputo, C. & Bolaños, P. Fluo-3 signals associated with potassium contractures in single amphibian muscle fibres. *J. Physiol.* **481**, 119–128. https://doi.org/10.1113/jphysiol.1994.sp020423 (1994).

55. Rueda, A., Fernández-Velasco, M., Benitah, J.-P. & Gómez, A. M. Abnormal Ca2+ Spark/STOC coupling in cerebral artery smooth muscle cells of obese Type 2 diabetic mice. *PLoS ONE* **8**, e53321. https://doi.org/10.1371/journal.pone.0053321 (2013).
## Acknowledgements

The authors thank Ms. R. Ori for her excellent technical assistance, Dr A. Aouacheria for observations of nuclei under a confocal microscope and Dr C. Dienes for his help in the resting membrane potential measurements. The research was supported by Grants from the French National Research Agency (ANR-13-BSV7-0010 and ANR-20-CE34-0017), the Ministry of Europe and Foreign Affairs (PHC Balaton 2019 #41908SC), the Hungarian Research Fund (2018-2.1.13-TÉT-FR-2018-00005) and by the Thematic Excellence Programme (TKP2020-IKA-04) which has been implemented with the support provided from the National Research, Development and Innovation Fund of Hungary, financed under the 2020-4.1.1-TKP2020 funding scheme.

## Author contributions

Conceived and designed the experiments: C.C., P.S. Performed the experiments: C.C., P.S. Analysed the data: C.C., P.S. Contributed reagents/materials/analysis tools: L.S., M.T., M.C. Wrote the paper: C.C., P.S., L.C.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-96028-w.

**Correspondence** and requests for materials should be addressed to C.C.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021

<footnote n="56">56. Klein, M. G., Lacampagne, A. & Schneider, M. F. A repetitive mode of activation of discrete Ca2+ release events (Ca2+ sparks) in frog skeletal muscle fibres. *J. Physiol.* **515**(Pt 2), 391–411. https://doi.org/10.1111/j.1469-7793.1999.391ac.x (1999).</footnote>

<footnote n="57">57. Shirokova, N., Garcia, J. & Ríos, E. Local calcium release in mammalian skeletal muscle. *J. Physiol.* **512**, 377–384. https://doi.org/10.1111/j.1469-7793.1998.377be.x (1998).</footnote>