Animals
The data presented in this paper were produced in two experiments performed in 2018 (n = 60 mussels) and 2019 (n = 30 mussels). Zebra mussels (Dreissena polymorpha) were collected by hand from the rivers Götzinger Ache (Bavaria, Germany) in March 2018 and Schinderbach (Bavaria, Germany) in July 2019. The mussels were fed ad libitum with shellfish diet in 2018 and with Chlorella vulgaris (SAG Number 211-19; Algae collection of the University of Goettingen, Goettingen, Germany) in 2019 before and during the experiments. The mussel size was measured before sacrificing according to [14] (mean length 23 ± 2.2 mm (mean ± standard deviation); width, 12 ± 1.5 mm; height 11 ± 1.2 mm).
All experiments were performed according to German animal protection regulations which do not require registration or approval of experiments using zebra mussels.
Exposure of mussels to radial extracorporeal shock waves
The mussels were exposed to rESWs produced with a Swiss DolorClast device (Electro Medical Systems, Nyon, Switzerland), using the radial handpiece and 6-mm applicator (Figs. 1, 2a). During the first/second experiment performed in 2018/2019, n = 10/n = 5 mussels each were randomly selected and exposed to 1000 rESWs each produced using an air pressure of the rESWT device of, respectively, 0 bar (sham exposure), 2.0, 2.5, 3.0, 3.5 or 4.0 bar.
For exposure to rESWs, the mussels were fixed under water in aquarium sand (diameter 2–3 mm; Dupla Marin Reef Ground; Dohse Aquaristik, Grafschaft-Gelsdorf, Germany) to disperse and, thus, minimize the reflection of rESWs (Fig. 1). Using a drill stand, the distance between the applicator tip and the mussels was set at 2.5 mm to prevent any mechanical destruction of the mussel valve through direct contact with the applicator tip. Accordingly, the energy flux density (EFD) at 3.0 and 4.0 bar air pressure that hit the mussels was approximately 0.08 mJ/mm2 and 0.11 mJ/mm2 (the EFD generated using the 6 mm applicator of the handpiece of the rESWT device shown in Fig. 1 is similar to the EFD generated using the 15 mm applicator of this device [15]; the decrease in the EFD is almost linear between a distance of 1 mm and 5 mm to the applicator [15]. At a distance of 1 mm and 5 mm to the applicator, the following EFDs were measured using the 15-mm applicator [11]: 0.1 mJ/mm2 and 0.04 mJ/mm2 when operated at 3.0 bar air pressure, and 0.14 mJ/mm2 and 0.06 mJ/mm2 when operated at 4.0 bar air pressure). The rESWs were applied at a frequency of 8 Hz.
Immediately after exposure to rESWs or sham exposure, the mussels were incubated in calcein solution (10 mg/l; Product Number: C0875-5G; Sigma-Aldrich, St. Louis, MO, USA) for 24 h. To this end, all mussels were placed in the same aquarium which contained six liters of calcein solution, with each group of mussels in a separate glass chamber (10 × 15 × 15 cm). The position of each glass chamber within the aquarium was selected randomly. Afterwards, the mussels were housed (using the same glass chambers and aquarium) in ventilated tap water for two weeks. Then, the mussels were euthanized in 70% ethanol, and the dissected valves were dehydrated in increasing concentrations of ethanol (70%, 80% and 90% for six days each, followed by 100% for 12 days).
After fixation, both valves of each mussel were degreased in xylene for six days, followed by incubation in methanol for six days. Then, the mussel valves were embedded in methyl methacrylate (Product Number: 800590; Sigma-Aldrich) according to [16]. Polymerization took 14 days. Afterward, the polymerized methyl methacrylate blocks containing the valves (one valve per block) were cut into 400-µm-thick sections along the longest axis of the embedded valve using a ring saw microtome (SP 1600; Leica, Wetzlar, Germany) (Fig. 2b, c). The sections were ground and polished using a 400 CS micro-grinder (EXAKT Advanced Technologies, Norderstedt, Germany). The final section thickness was approximately 200 µm, measured in the middle of each section using a digimatic micrometer (Mitutoyo, Kawasaki, Japan).
Measurements of fluorescence signal intensity
Images were taken using a fluorescence microscope (Olympus BX51WI; Olympus, Tokyo, Japan) using a 4 × UPlanSApo objective (numerical aperture = 0.16) (Olympus), Alexa Fluor 488 filter (49011; Chroma, Bellows Falls, VT, USA), grayscale EM CCD camera (Model C9100-02, 1000 × 1000 pixels; Hamamatsu Photonics, Hamamatsu City, Japan) and SOLA LED lamp (Lumencor, Beaverton, OR, USA). All images were taken with the Stereo Investigator software (64 bit, Version 11.07; MBF Bioscience, Williston, VT, USA) and saved as 8 bit TIF files (i.e., with gray values ranging from 0 to 255). Using pilot measurements, the camera was adjusted so that no image was overexposed (i.e., all gray values were smaller than 255). This resulted in the following camera settings: exposure time, 24 ms; sensitivity, 80; gamma, 1.0.
In line with our previous study [10], the strongest fluorescence signal was found over the hypostracum (Fig. 2d). Analysis of mussels after sham exposure indicated that the signal over the hypostracum was indeed caused by exposure to rESWs (Fig. 2e, f). Accordingly, measurements of fluorescence signal intensity were performed over the hypostracum, using the linear pixel plot function of the Stereo Investigator software (MBF Bioscience). Four measurement lines each (spanning 243 ± 79 µm representing 135 ± 44 pixels, depending on the curvature of the valve) were positioned over the hypostracum as shown in Fig. 2d, representing Regions A-D indicated in Fig. 2a. Region A was next to the umbo, Region D was next to the shell growth zone, and Regions B and C were in between. As in our previous study [10], the umbo itself was excluded from the analysis because of strong autofluorescence of the ligament.
Statistical analysis
For each group of mussels (i.e., each intensity of the rESWs) mean and standard deviation of side- and region-specific fluorescence signal intensities were calculated. Outliers were identified using the Tukey's fences method [17] (with k > 1.5 indicating an outlier) and removed (outlier values were most probably caused by the methodology used for generating the sections, in particular by grinding and polishing). The corresponding calculations were performed using GraphPad Prism (Version 9.2.0 for Windows; GraphPad Software, San Diego, CA, USA). Fifty-seven out of the 720 individual data (six groups of mussels × 15 mussels per group × two valves per muscle × four regions per valve) (7.9%) were identified as outliers. The absolute and relative numbers of valves with 0/1/2/3/4 outlier values in their respective group were 145/21/9/2/3 and 80.6%/11.7%/5.0%/1.1%/1.7%, respectively. After removal of outliers, there were at least 12 (out of 15 maximally possible) values available for each combination of energy, side and region.
Then, differences in mean fluorescence signal intensities were investigated using general linear model analysis, with energy (i.e., the intensity of the rESWs), side (left/exposed vs. right/unexposed) and region (regions A–D shown in Fig. 2a) as fixed factors and the averaged fluorescence signal intensities (one value each per mussel, side and region) as depending factor. Post hoc analyses (energy, region) were performed using Bonferroni's multiple comparison test. Calculations were performed using SPSS (Version 26.0.0.0; IBM, Armonk, NY, USA). P values smaller than 0.05 were considered statistically significant.