Abstract

Caenorhabditis elegans nematodes are broadly used to investigate the impact of environmental factors on animal physiology and behavior. Here, C. elegans with internalized paramagnetic nanoparticles were placed inside a magnetic field (MF) to explore its effects on locomotion. We hypothesized that internalized paramagnetic nanoparticles combined with external MF affect C. elegans’ locomotion machinery. To test our hypothesis, we used adult C. elegans fed on bacteria mixed with paramagnetic nanoparticles of 1 μm, 100, and 40 nm diameter. The presence of nanoparticles inside the worms’ body (alimentary canal, body muscle) was verified by fluorescent and electron microscopy. A custom-made software was used to track freely moving C. elegans in the absence or presence of MF, sequentially, for 200 + 200 s. We used established metrics to quantify locomotion-related parameters, including posture, motion, and path features. Key attributes of C. elegans locomotion (stay ratio, forward over backward motion, speed) were affected only in worms with internalized nanoparticles of 100 nm in the presence of MF (reduced speed, increased stay ratio, decreased forward/backward ratio), in contrast to untreated worms. Our work shows that internalized particles of specific properties affect C. elegans locomotion under MF. Hence, it contributes to clarifying the effects of MF and activated nanoparticles on C. elegans locomotion, thus fueling further research.

1. Introduction

The effects of magnetic field (MF) on living organisms have been a target of numerous research efforts, with their number increasing significantly during the last decades [13]. In addition to the interest the scientific community shows on the impact of alternating MF on cells [46], the effects of static MF have also gained attention [7, 8], often with respect to complicated cellular processes like apoptosis [9]. The least studied type of MF, with respect to living cells, is MF of spatial gradient, even though the theoretical framework for the possible impact of high-gradient MF of various sources on cell’s molecular components and function is known [10].

Model organisms have been a successful resource to study MF effects on various types of cells and tissues [1114]. Invertebrate models, like Drosophila melanogaster, have been used since the 1980s [1520]. The nematode C. elegans, an emblematic model organism to study the impact of environmental factors on behavior and physiology [2125], was used in MF-related work only recently [2630]. Excitingly, the first animal magnetosensory neurons were identified in C. elegans [26], which were found to orient to the earth’s MF during vertical burrowing migrations, reflating interest in the field [3134]. The presence of biogenic magnetite has also been reported in C. elegans [35].

Nanoparticles uptake by C. elegans worms has been a successful means to evaluate the toxicity of heavy metals and pollutants [3638], and the importance of C. elegans as a model system for in vivo nanoparticle assessment has been highlighted [33, 39]. Worms’ behavior [40] and locomotion [33, 41, 42] have been evaluated under the influence of internalized metal nanoparticles, yielding interesting results. In addition, magnetized nanoparticles have been used to activate ion channels in C. elegans through heating [43], and internalized nanoparticles were used to locally enhance MF in the worms’ body to study the subsequent impact on its metabolism [44].

C. elegans locomotion has been a major behavioral output used to investigate the impact of genetic background, environmental factors, and diverse treatments on the worm’s nervous system [4550]. Locomotion features have been characterized and quantified extensively and are being used as an indicator of C. elegans physiological status and health span [5154]. With respect to C. elegans’ locomotion in the presence of MF, besides the impact of the earth’s MF [26, 55], the literature is sparse. It is reported that nematodes are not affected by mild static MF [27] and that their locomotion can be manipulated by magnetically responsive structures [56], but other than that not much is known about the impact of MF on C. elegans’ locomotive machinery.

The need to clarify MF effects on animal physiology, and the increasing interest in the sensitivity of C. elegans to MF despite the current gaps in the literature, indicate that a detailed investigation of MF gradient effects on worms’ behavior can provide important insights. In this context, we hypothesize that internalized paramagnetic nanoparticles combined with external MF affect C. elegans’ locomotive behavior. To test our hypothesis, we use internalized paramagnetic nanoparticles to generate an effect inside the worms’ body, in the closest possible proximity with tissues and cells. We select locomotion as a quantifiable behavioral expression to determine the effect of MF combined with internalized nanoparticles because of its role in animal fitness and as a reliable readout of physiological changes. Our results demonstrate a response of C. elegans locomotion machinery to internalized paramagnetic nanoparticles in combination with MF gradients, and they pave the way for future studies seeking to clarify the participation of excitable cells, muscles, and potentially even neurons, to this uncharacterized behavior.

2. Results

2.1. Magnetic Field Gradient Characterization

Simulation results for the two electromagnets used in the study agree with experimental results provided by the manufacturer (Figure S1). A schematic of the MF around the electromagnets and the computational domain of the COMSOL model are presented in Figure 1(a). A view of the MF on the plane of the worm plate surface is presented in Figure 1(b). The direction of the MF between the two electromagnets is indicated by arrows. Contours of MF intensity confirm that the MF is stronger near the electromagnets, as expected.

We focused on the MF and the forces it generates on particles located on the experimental plate surface, where the worms’ locomotion takes place. The MF was almost one-dimensional on the plate surface and was strongest near the electromagnets (Figure 1(c)). There were nine components for the gradient of the MF. In Figure 1(d), the strongest component of the MF gradient is shown (Bxx), which is parallel to the x-axis. The magnitude of the gradient was the largest near the electromagnets.

Nanoparticles create secondary MFs and forces in the presence of an external MF. Details regarding the computation of forces created by nanoparticles can be found in the Supplementary Information. The MF flux values for three nanoparticles located along the x-axis and separately the y-axis were calculated in MATLAB using Equations (2) and (3) of Supplementary Information (Figure 2(a)2(d)). As expected, the MF was strongest close to the particles for both locations and decayed rapidly as the distance from the particles increased (Figure 2(e)). The force between particles located along the x-axis was attractive, while the force between the particles located along the y-axis was repulsive. The attractive forces between the particles allow them to form chain-like structures if they are not interrupted by the medium in which the particles are located [57, 58]. The magnetic moment of the external MF for the particles on the worm plate surface is depicted in Figure 2(f).

2.2. Confirmation of Nanoparticles Uptake and Particle Localization in C. elegans’ Body

Nanoparticles mixed with bacterial food (Escherichia coli OP50) were successfully internalized, as verified by microscopy methods and selected according to the properties of each particle group (Figure 3). The presence of 1 μm paramagnetic particles (Table 1) in the worm’s intestine and in the pharynx around the grinder area was verified by bright field microscopy. The particles appear as dark (copper-colored) aggregates accumulated in the alimentary canal (Figure 3(a), right panel), whereas the control animals’ intestine area appears transparent (Figure 3(a), left panel). Uptake of 100 nm fluorescent, paramagnetic particles (Table 1) was verified by fluorescent microscopy. The particles appear to accumulate along the intestine lumen and in the pharynx, as shown when filters for rhodamine, the fluorescent substance with which the particles were coated (see Table 1 for particles’ properties), were used (Figure 3(b)). Successful feeding on 40 nm paramagnetic particles (Table 1) was confirmed by scanning electron microscopy. When using the circular backscatter detector, 40 nm particles were visualized as white dots under the worm cuticle, in the broad area downstream of the pharynx and along the alimentary canal (Figure 3(c)). The deformation of the sample due to the process followed allowed only for estimating the approximate location of the particles. The white dots which represent the particles appear in different sizes, which might be attributed to particle aggregates or to the different depths at which the particles were located.

To investigate the particles’ localization in the worms’ body and to explore whether they passed the intestine barrier, we used transmission electron microscopy (TEM). TEM images show that particles can be found in the intestine (40 nm, Figure 4(d)) and the intestine lumen (100 nm, Figure 4(e)). Interestingly, particles’ aggregates were also detected inside muscle tissue, very close to the body wall (40 nm, Figure 4(f)). Therefore, the nanoparticles can be in close proximity to excitable cells.

2.3. Analysis of C. elegans Locomotion

We examined posture, motion, and path features [59]. The posture feature quantified was the number of bends (bend count), the motion features included are related to motion state (forward/backward ratio and stay ratio) and traveling speed (, , and ), and the path features were the path curvature and the range (see Section 5).

First, we explored whether there are differences in locomotion features due to the presence of nanoparticles alone in the absence of MF. We compared the three groups of particle-fed nematodes to untreated worms when the MF was OFF (Figure S4). No statistically significant difference was detected between untreated and particle-fed animals in the absence of MF (Figure S5).

Next, the particle-fed nematode groups were analyzed during the OFF and ON states of the MF. Results show that only worms fed with 100 nm diameter paramagnetic nanoparticles were affected as they moved freely in MF, compared to nematodes without internalized nanoparticles (Figure 5). More specifically, the number of bends (bend count) realized per worm was not affected in any of the nematode groups (Figure 5, superpanel (a)).

The forward/backward ratio (time the worm was in the forward/backward state over the total recording time) of Group 100 worms decreased when the MF was ON (Figure 5(b1)) and their stay ratio (the time the worm was in the paused state over the total recording time) increased (Figure 5(c1)), indicating an overall decrease of motility. Moreover, the speed of worms fed with 100 nm particles decreased when the MF was ON (Figure 5(d1)). There was no significant effect detected in the components of the and components (Figure 5, superpanels (e) and (f)). A trend for decreased is apparent, but not statistically significant, in all experiments of Group 100 (Figure 5(e1)).

None of the path features, namely path curvature and range, was found to be affected in a significant and consistent way in any of the groups tested (Figure 5, superpanels (f)–(h)).

3. Discussion

3.1. Internalized Nanoparticles Alone Do Not Affect C. elegans Locomotion

C. elegans nematodes are shown to be an excellent model for evaluating the impact of nanomaterials on animal physiology and behavior [60]. More specifically, metal nanoparticles of various types have been used to evaluate particle toxicity using C. elegans [42, 6165]. Particle coating and size, worm developmental stage, and duration of exposure have been shown to the affect translocation of particles in various tissues of the worm’s body [42, 61, 66]. Particles used in the present study are larger [42, 6163], and worms have been exposed to them for a shorter period (18–20 hr) than elsewhere [42, 67]. These differences may explain why most of the nanoparticles are found along the worms’ pharynx, upper intestine (Figure 3(a) and 3(b)), and lower intestine area (Figure 3(b)). The location of 40 nm particles in worms of Group 40 around the pharynx and grinder area (Figure 3(c)) is only approximate, as some deformation has been induced on the sample during preparation, and scanning electron microscopy (SEM) allows for detecting objects that are close to the body surface. Particles’ aggregates are also found in the intestine (Figure 4(d)), intestine lumen (Figure 4(e)), and in muscle tissue (Figure 4(f)).

Worms used in the present work are young adults and not larvae [42, 66, 67], which, as developing organisms, could be more vulnerable to toxic effects [68]. Moreover, particles used in the experiments presented here are larger [42, 6163] and made of different metals [41, 63], compared to the ones utilized in similar studies, which may result in different abilities to overcome the intestine barrier or translocate to other tissues, as well as to different toxicity per se.

3.2. Internalized Nanoparticles, Combined with MF, Have a Particle-Specific Effect on C. elegans Locomotion

Magnetotaxis in C. elegans has been demonstrated [26], with the participation of AFD neurons, the first to be identified as magnetosensory in the animal kingdom. It was suggested that endogenous magnetic material, previously reported in C. elegans, may be also involved [26, 35]. These findings sparked a vivid discussion in the C. elegans community [31, 6971]. In our experiments, the locomotive behavior of Group C worms, which did not contain any particles, was not affected by the externally applied MF (Figure 5).

At the same time, the presence of 100 nm-diameter internalized nanoparticles had an impact on specific locomotion features (Figure 5) when an MF was applied. Indeed, Group 100 worms move more slowly (Figure 5(b1)), spend less time moving, either forward or backward (Figure 5(e1)), and spend more time in the paused state (Figure 5(f1)). Therefore, the overall motility of these nematodes is decreased [33]. Modulation of the locomotive circuitry, including the switch from a more mobile to a less mobile state, is known to be mediated by dopaminergic neurons [72, 73], like the PDE mechanosensory neuron. Our findings set the stage for investigating one more role for dopamine as a potential regulator of behavioral plasticity in 100 nm-diameter particle-fed C. elegans inside MF.

It is possible that worms slowed down when they found themselves in a particular orientation inside the MF or when the internalized particles obtained a particular orientation with respect to the MF. As indicated in Figure 1, the properties of the MF change significantly in the direction of the MF. Therefore, any effect the MF may have on the particle-fed worms or on the internalized particles themselves could be changing when worms are moving along the MF.

There was no effect detected in Group 40 in any of the metrics examined. This can be due to the smaller MF or smaller gradient of MF of the particles in this group. Therefore, even though they are found in the intestine and close to the body wall muscles (Figures 4(d) and 4(f)), they are not potent enough to generate any detectable effect. Similarly, there was no impact detected in Group 1 nematodes. This suggests that it may be difficult for particles of this size to cross tissue barriers in the worms’ body and end up close to excitable cells. Indeed, although aggregates of 1 μm particles were spotted in the alimentary canal (Figure 3(a)), they were not detected at other locations by electron microscopy.

In parallel, we have been able to estimate the MF and the gradient of MF only for the 1 µm particles due to lack of available information on the magnetic properties of the 40 and 100 nm particles. Nonetheless, we can compare the MF and gradient of MF between the particles based on their size (see Supplementary Information). We found that smaller particles have a larger gradient of MF compared to larger particles in their proximity. The overall impact each particle type has on worms’ physiology depends on the magnitude and the gradient of the MF. Both depend on the material properties of the particles, which determine the magnetic moment. Our experimental observations show that the strongest effect among the three types of studied particles occurred in the case of 100 nm particles.

The particles’ coating (Table 1) was not expected to affect their magnetic behavior. It could, however, impact their interaction with cells. Since the magnetic and physical properties of the particles are the most influential regarding the secondary MF effects, we focused on the particle size for our data analysis. In addition, the thickness of the coating was small (a single monolayer of streptavidin for the 1 μm particles and a 2–3 nm thick layer of polymer for the 40 nm particles, according to the manufacturers), which means that the magnetic core could still affect cells and tissues close to the particle. The experimental procedure did not allow us to know the quantity of particles ingested by each individual worm, nor the exact location and the precise interactions of particle aggregates in each individual nematode. Therefore, we cannot draw conclusions based on the particles’ quantity or exact aggregate location inside each worm tested.

Ideally, results shown in Figures 1 and 2 should be combined to assess the synergistic properties of the external component of the MF (generated by the electromagnets) and of the secondary component of the MF (generated by the magnetized particles). However, the MF generated by the particles is very localized (in the microscale), as shown in Figure 2. Thus, we discuss below the potential effects of the MF induced by the electromagnets and of the secondary MFs separately and on their own respective scale. We also discuss the effect of particles’ aggregates behavior under the influence of an external MF.

3.3. Factors that Could Shape the Impact of Internalized Particles in the Presence of External MF

The forces that were created either by the external MF or the paramagnetic particles themselves were small (Figure 2(e)), and they were not strong enough to mechanically push the worms to move along their track. This suggests that some other mechanism is responsible for the detected changes in worms’ locomotion. The magnitude of the external and secondary MFs had the same order of magnitude, and they were small (Figures 1 and 2(a) and 2(c)). However, the gradient of the MF fields in the vicinity of the particles was substantially large (Figures 1(b) and 1(d)).

Effects were likely to be more pronounced where the external MF was stronger since that would result in stronger secondary MFs generated by the particles (until the magnetization of the particles becomes saturated). Therefore, the locomotion of worms crawling under stronger MF, namely near the electromagnets (Figure 1), was more likely to be affected. The spatial distribution of the MF can partly explain the variations observed in the experiments. In addition, the presence of a worm crawling into areas of higher magnetic flux, therefore of higher impact, is a random event. The variations in the values of the metrics examined might mirror that randomness. Lastly, since we had no direct control over the location where particles resided in the worms’ body, their presence in body locations where they could affect mechanosensitive ion channels could be another explanation for the variability that we observe in our experiments (indicative cases: Figure 5(b1), ON1; Figure 5(c1), ON1; Figure 5(c2), ON1; Figure 5(c3), ON1).

Zablotskii et al. [10] provide several examples in which the gradient of MF can affect cellular and subcellular mechanisms. The gradient of the secondary MF fields obtained from our simulations (up to 2 × 105 T/m for the 1 μm particles, Figures 2(b) and 2(d)) was well above the threshold that Zablotskii et al. [10] suggest may impact cells with mechanosensitive ion channels (103 T/m, see also [10]). The gradient of the secondary internal MF fields was also above the threshold the same authors pose for magnetically induced changes in gene expression; however, we consider such a possibility highly unlikely in our case due to the very short time the external MF was applied (∼90 s). Therefore, the secondary field generated by the paramagnetic particles upon application of external MF could lead to local gradients of MF inside the worm’s body, large enough to interfere with the functionality of excitable cells (e.g., body muscle cells, see also Figure 4(f)).

This could have happened by affecting the cells’ ion channels, provided that the particles were very close or even in contact with the cells’ membrane. Experimental data presented here cannot provide insight into exactly which cells might have been the target of the observed MF effect. However, TEM findings (Figure 4) showed that particles’ aggregates could be located in body wall muscles (Figure 4(f)). Therefore, the possibility that excitable, e.g., muscle cells, were affected by ingested particles in the presence of MF is considerable. In addition, it is possible that mechanosensitive neurons are excited as well, given the conceivable effect of MF on mechanosensitive ion channels, as described earlier [10]. These could be dopaminergic neurons, like PDE, ADE, or CEP [73, 74], or mechanosensitive neurons that innervate the body wall muscles, like PVD [75]. Moreover, the impact on the intestine (Figures 4(d) and 4(e)) could affect C. elegans physiology and locomotion dynamics, given the multiple roles of this complex tissue [76, 77]. Further experiments are needed to clarify the mechanism behind the observed changes.

The possibility of MF having a direct action on the magnetosensitive neurons described in C. elegans [26] cannot be excluded. The neurons reported to have a magnetosensory function are the AFD sensory neuron pair [26], which are also known to respond to temperature [78], CO2 [79], and moisture [80] gradients. Interestingly, AFD neurons can transmit both stimulatory and inhibitory signals to the interneuron AIY, operating as a behavioral regulator in a thermosensory circuit [81]. Whether AFD contributes in a similar way to the behavioral dynamics reported here could be explored in future experiments via calcium imaging of freely moving nematodes under the reported experimental conditions. In the present study, we did not have direct evidence of sensory or motor neurons being affected by MF, neither was that possibility explored.

When an MF is applied to a population of paramagnetic particles, the particles self-organize into arrays, columns, or chains, depending on the nature of the applied MF and the properties of the particle-containing medium [57, 8284]. In our experiments, when the MF is turned ON, it is likely that the internalized particles start moving as they organize into self-assembled structures. It is possible that this motion applies pressure on or stretches the surrounding tissue, resulting in a decrease in overall motility.

4. Conclusions

We demonstrate a particle-specific effect of internalized paramagnetic nanoparticles, in combination with externally applied MF, on the dynamics of C. elegans’ locomotion. Established locomotion metrics, i.e., speed and motion state, are significantly different between untreated worms and worms treated with 100 μm particles when moving inside an MF. In parallel, no significant difference is detected between untreated and particle-treated worms in the absence of an MF. Possible explanations on the mechanism that leads to the observed results are related to the effects of MF gradients on cells [10, 85] mediated by magnetic nanoparticles [86]. These are also related to the formation of paramagnetic particles’ aggregates inside MF [8284]. The exact mechanism that underlies the observed effect is not clear. Our experiments are the first to comparatively explore the effect of four different types of paramagnetic particles in young adults of an established nematode model system. Our findings pave the way for more targeted experiments on the sensitivity of animals to MF gradients. C. elegans nematodes, due to their excellent trackability, could play a key role in such an effort.

5. Materials and Methods

5.1. Magnetic Field Characterization

COMSOL Multiphysics (COMSOL, Burlington, MA) software was used to characterize the MF that is generated by the two electromagnets in the experimental setup (Figure S2). The data for the magnetic flux density of the electromagnets (available from the manufacturer) was used to calibrate the parameters of the electromagnets in COMSOL Multiphysics. The COMSOL Multiphysics model was also used to estimate the intensity of the external MF, the gradient of the external MF, and the forces that are applied on paramagnetic particles by the external MF. MATLAB was used to calculate the forces applied on the paramagnetic nanoparticles. More details are given in the Supplementary Information section.

5.2. Nanoparticles Internalization

We investigated the locomotion of four groups of young adult wild-type N2 C. elegans hermaphrodites fed on (see also Table 1):(i)plain bacterial food source E. coli OP50, control animals-Group C, n = 30,(ii)E. coli OP50 mixed with 1 μm-diameter paramagnetic particles (Dynabeads MyOne Streptavidin C1, Invitrogen, Thermo Fisher Scientific, USA), Group 1, n = 32,(iii)E. coli OP50 mixed with 100 nm-diameter paramagnetic particles (nanomag-CLD-red, Micromod Partikeltechnologie GmbH, Germany), Group 100, n = 33, and(iv)E. coli OP50 mixed with 40 nm-diameter paramagnetic particles (iron oxide nanocrystals, Ocean NanoTech, USA), Group 40, n = 25.

In all cases, particles were isolated from the initial suspension by brief centrifugation and were resuspended in OP50 in a final concentration of 0.5 mg/ml OP50-particle mix. Freshly made 60 mm standard nematode growth medium (NGM) plates were seeded with 100 μl of plain OP50 or OP50-particle mix. Plates were left to dry overnight at room temperature (RT), and ∼20 worms were transferred in them the next day. Nematodes were left to feed on the plain or enriched bacterial lawn for 18–20 hr at 20°C. Then, they were either prepared for microscopy or 12–15 of them were transferred to a fresh, unseeded 35 mm NGM plate for locomotion recording. In the second case, worms were left to acclimatize in the new plate for ∼15 min before recording. There were three reasons for transferring worms to a new, smaller plate. First, we wanted the worms to experience the effect of only internalized nanoparticles under MF and not of the remaining particles on the plate surface. Second, the presence of enriched bacterial lawn on the plate surface interfered with the tracking algorithm and could have affected the worms’ locomotion due to its viscosity. Third, by using 35 mm plates, we decreased the distance between the electromagnets and the worms (as shown in Figure 6) so that the worms experience a stronger external MF.

For each group tested, experiments were run over three different experimental days. Therefore, on each experimental day, we processed 10–14 worms for a specific group. These 10–14 worms were treated simultaneously (on the same plate), and each one of them is considered a biological replicate.

5.3. Fluorescent, Scanning, and Transmission Electron Microscopy (SEM and TEM)
5.3.1. Fluorescent Microscopy

Worms were transferred to an unseeded NGM plate and were washed with 0.5 ml of 1× phosphate-buffered saline (PBS). Next, they were transferred to a glass slide, where they were anesthetized on fresh agar pads [87] using 10 mM NaN3 [88]. Samples were imaged using a BX51WI Olympus fluorescent microscope (Olympus, Tokyo, Japan) coupled with an ORCA-flash4.0 camera (Hamamatsu Photonics, Hamamatsu City, Japan).

5.3.2. Scanning Electron Microscopy

Samples were prepared as described previously [87, 89], with modifications and dissection omitted. Briefly, worms were transferred to an unseeded NGM plate and were washed with 0.5 ml of 1× PBS. Next, they were transferred to a glass coverslip and were anesthetized using 10 mM NaN3 [88]. Samples were imaged using FEI Helios 650 nanolab SEM/FIB (FEI, Thermo Fisher Scientific, Waltham, MA).

5.3.3. Transmission Electron Microscopy

Samples were prepared based on the literature [90, 91], with modifications. Briefly, tissues were fixed in pre-warmed 3.2% paraformaldehyde, 0.2% glutaraldehyde in 0.1 M sodium cacodylate buffer for 1 hr at RT and incubated overnight at 4°C. Next, they were rinsed 3 × 10 min with 0.1 M sodium cacodylate buffer, and post-fix in 2% OsO4 in 0.1 M sodium cacodylate buffer for 1 hr at RT. Another rinse 3 × 10 min with 0.1 M sodium cacodylate buffer followed, and then worms were embedded in resin mold in histogel. Samples were dehydrated for 15 min each in 50%, 70%, 90%, 95%, and finally, two changes of 100% ethanol, cleared in two 15 min changes of propylene oxide and infiltrated in propylene oxide:epon (Embed812), as follows: (a) 3:1, 1 hr, (b) 1:1, 1 hr, (c) 1:3, 1 hr, (d) full strength, 2 hr or overnight, two changes. Finally, samples were embedded in beam capsules in full epon and polymerized at 60°C for 24 hr. Samples were imaged using a JEM-1400-plus transmission electron microscope (JEOL, Peabody, MA).

5.4. Worm Recording and Tracking
5.4.1. Recording

A 35 mm plate containing worms of a specific group was placed between the two electromagnets, as shown in Figure 6 so that the plate surface and, therefore, the worms were positioned close to the center of the electromagnets. First, a 200 s movie (1 frame/s) was recorded in the absence of MF (OFF state), and immediately after, a second 200 s movie (1 frame/s) was recorded with the MF on (ON state), using QCapture Pro software (QImaging, Surrey, Canada) and a Micropublisher3.3 RTV camera (QImaging, Surrey, Canada), mounted on an Olympus SZ61 microscope (Olympus, Tokyo, Japan). We recorded over 200 s intervals because we were interested in detecting the transient effect of MF on locomotion dynamics. The two electromagnets used were a 4.0″ Dia. Electromagnet, 12 VDC, and a 3.5″ Dia. Electromagnet, 12 VDC, both from APW Company, Rockaway, NJ. Electromagnets were operated at 1.67 and 3 A, respectively, as indicated by the manufacturer, using a 1762 DC power supply (BK Precision, Yorba Linda, CA). By using a noncontact infrared thermometer (Omega Engineering, Norwalk, CT) we verified that the plate surface temperature remained constant throughout the recording period.

5.4.2. Tracking

Every movie was imported to MATLAB (MathWorks, Natick, MA) for postprocessing. Each worm was tracked individually. To this end, we developed a custom tracking code in MATLAB (Figure S3). In the first step, all frames were used to construct the movie background, which consisted of all the objects that did not move for long periods of time during the entire recording. Then, each frame was subtracted from the background to extract the foreground, which consisted of all moving objects. Next, the user was prompted with the initial frame of the movie, of which the background had already been subtracted, to select the worm to be tracked by the software. After the user selected the worm, the code created a small examining frame around it and excluded the targeted worm from the rest of the movie frame. Then, the cropped figure was converted to a binary image. After the binary image was enhanced, the shape of the binary object, i.e., the worm, and its global position were stored. Next, the code proceeded to the next movie frame and used the extracted global location of the worm as the center of the small examining frame. The small examining frame must be large enough to capture the motion of the worm in two successive movie frames. Since there was more than one worm freely moving in each experiment, there were occasions in which more than one object were included in the small examining frame. For such occasions, the user was prompted by the code to manually indicate again the worm to be tracked. This way, the worm that was initially selected to be tracked was always encapsulated by the examining frame. The code continued the tracking process until the last frame of the movie was processed, and it stored the shape of the worm and its global location for each frame. Once finished, the user runs the code again to track another worm.

The tracking code can be found at https://github.com/mirzakhalili/MultiWormTracker.

5.5. Locomotion Analysis

The following features of C. elegans morphology and experimental setup properties were used for the quantification of C. elegans locomotion parameters.

5.5.1. Morphology Features

(1)Length: the worm length was defined as the chain-code pixel length of worm skeleton, which was converted into mm.(2)Centroid: the worm density was assumed to be constant throughout its body, so the centroid of mass was the same as geometric centroid. Since the swing of the head or tail (first or last 1/12 chain-code length part of the worm) can significantly influence centroid determination, they were ignored when computing the centroid.

5.5.2. Setup

(1)Coordinates system: the x-axis was set along the direction of the MF, between the two electromagnets, and z-axis was normal to the plate, pointing upward (Figure 6). Thus, by applying the right-hand rule, the coordinates system was established. Since we did not identify head/tail orientation for the worms, the coordinates system was important for the detection of directionality.(2)Unit conversion: any feature regarding length was derived first in pixels. With a known length recorded with the same experimental setup, the conversion between pixels and microns was determined.

Locomotion-related parameters of interest were divided into three categories: posture features, motion features, and path features, as described extensively by Yemini et al. [59], with minor modifications. A brief description of the examined features follows below.

5.5.3. Posture Feature

(1)Bend count: this metric () corresponds to the number of bends along a single worm. First, the supplementary angles (Figure S4) were computed along the worm skeleton. Next, a Gaussian filter over each 1/12 of the chain-code length of the skeleton was applied to the supplementary angles to smooth out any high-frequency changes and was then normalized. The filter had a constant proportional to the reciprocal of the standard deviation, α = 2.5. By checking the sequence of supplementary angles, the bend count was incremented whenever the angle reached 0° or changed sign. The check started from the first 1/12 segment to the last 1/12 segment to ignore small bends near the tail and the head.

5.5.4. Motion Features

(1)Motion state: worm’s motion state can be divided into two types, the forward/backward state and the paused state. The worm was considered to be in the forward/backward state when its instantaneous speed was greater or equal to 5% of its mean length per second, and it was considered in the paused state when the instantaneous speed was less than 5% of its mean length per second. Therefore, the ratio of the time the worm was in the forward/backward state over the total recording time, namely the , and the ratio of the time the worm was in the paused state over the total recording time, namely the stay ratio , were calculated.(2)Speed: velocity is defined as the signed difference between a single worm’s centroids of two sequential frames in the coordinate over the time gap between two frames (1 s) and can be projected on two orthogonal axes x and y in the plane of the plate (Figure 6). The absolute value of velocity and its components gave , , and , respectively.

5.5.5. Path Features

(1)Path curvature: this metric is defined as the angle, in radians, of the worm’s path divided by the distance traveled in microns. Three successive frames were used to approximate the start, middle, and end of the worm’s instantaneous path curvature. The angle was measured by the difference in tangent angles between the second to last frame centroid and the first to the second frame centroid. Then, the path curvature was obtained by dividing the angle by the distance between the first and last centroid.(2)Range: range is defined as the distance between the worm’s centroid and the centroid of the worm’s path in each frame.

5.6. Statistical Analyses

Locomotion features were analyzed using parametric tests since the Anderson–Darling normality test p-value was larger than 0.05 for all samples, thus confirming the normality null hypothesis that the data are not inconsistent with a Gaussian distribution.

For each group or nematodes, each group being fed on a particular type of nanoparticles, as described earlier, tested worms were tracked in three independent experiments conducted on three experimental days (OFF1/ON1, OFF2/ON2, and OFF3/ON3). Therefore, for Group 100, animals fed on E. coli OP50 mixed with 100 nm-diameter iron core paramagnetic particles, n = 33 (individual experiments: n1 = 13, n2 = 11, n3 = 9); for Group 1, animals fed on E. coli OP50 mixed with 1 μm-diameter paramagnetic particles, n = 32 (individual experiments: n1 = 14, n2 = 10, n3 = 8); for Group 40, animals fed on E. coli OP50 mixed with 40 nm-diameter iron core paramagnetic particles, n = 25 (individual experiments: n1 = 10, n2 = 7, n3 = 8); and for Group C, namely control animals, fed on plain food source E. coli OP50, n = 30 (n1 = 11, n2 = 10, n3 = 9).

For each group, we compared results from the ON state and OFF state for each experiment separately, and for all three experiments combined (Figure 5). Comparisons were made by applying Student’s t-test (unpaired for individual experiments, paired for combined data), and differences were considered statistically significant when . We considered a locomotion feature to be affected by internalized paramagnetic particles in the presence of MF only when there was a statistically significant difference detected on at least two experiments and when comparing combined data.

To find out whether there are differences in locomotion features due to the presence of nanoparticles alone (in the absence of a MF), we compared the three groups of particle-fed nematodes to untreated worms, when the MF is OFF. Comparisons were made by one-way analysis of variance (ANOVA), graphs and ANOVA results are presented in Figures S4 and S5. No statistically significant difference was detected between untreated and particle-fed animals, in the absence of MF.

All analyses were performed in GraphPad Prism 9.0 (GraphPad Software, San Diego, CA, USA).

To design the experiment, we run a statistical power analysis, using opensource software (Figure S6). We prepared the sample size used in the experiments based on this estimation. We slightly increased the sample size number to ensure that we have enough worms, since some of them might be injured or lost during the process.

Data Availability

All data supporting the results are available upon reasonable request; code is available on GitHub as described in the manuscript.

Disclosure

The manuscript was already published as a poster [92] based on the link https://www.biorxiv.org/content/10.1101/248369v2. Ehsan Mirzakhalili is presently at Bioengineering, University of Pennsylvania, Philadelphia, PA. Yang Zhang is presently at Robotics System Development, Carnegie Mellon University, Pittsburgh, PA.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

BE and EG conceived the idea; EG, EM, and BE designed the experiments; EG, EM, and YZ ran experiments; EM created tracking algorithm and ran simulations; YZ ran tracking algorithm, processed, and analyzed recordings; EG, YZ, and EM collected and analyzed data; EG and EM wrote the paper, with input from YZ. All authors reviewed and edited the manuscript and gave final approval for publication.

Acknowledgments

We thank Nikos Chronis, Kenn Oldham, and Jinhong Qu for the use of selected equipment and Syeda Maisa for helping with preliminary worm videos. Scanning electron microscopy was performed at the Electron Microbeam Analysis Laboratory (EMAL) with support from the University of Michigan College of Engineering; we thank John Mansfield, Kai Sun, and Haiping Sun for the training. Transmission electron microscopy was performed at the Microscopy & Image Analysis Laboratory (MIL) at the University of Michigan Medical School; we are grateful to Jeff Harrison and Pennelope Blakely for their help and guidance. We thank Surojit Sural for his input on power analysis and to Hannah Seidel and Elisa Frankel for feedback. We are particularly grateful to Hong Zhan for his help with TEM images interpretation. This work was supported by the Division of Civil, Mechanical and Manufacturing Innovation, National Sciences Foundation (award #1334908 to B.E.), and the University of Michigan Office of Research (grant #U055203 to E.G.). EG is the recipient of NIH–NIA award K01-AG057833. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Supplementary Materials

Figure S1: Calibration of the parameters in COMSOL Multiphysics simulations to match the available data for magnet 1 (a) and magnet 2 (b) that are used in the experiments. Figure S2: All the components of the gradient of the magnetic field (MF) for the particles in the vertical (along the y-axis) and horizontal (along the x-axis) configuration. Figure S3: Locomotion analysis process. Figure S4: Comparison of particle-fed nematodes with untreated nematodes, for all locomotion features, in the absence of MF (OFF). Each graph presents combined data for all three experiments for each worm group for one locomotion feature. Figure S5: One-way ANOVA for the comparisons between all three groups of particle-fed animals against the untreated (control) group for the OFF state. OFF state for each of the three groups was compared against the control for each locomotion feature (see also Figure S4). Figure S6: Power analysis for locomotion dynamics experiments. (Supplementary Materials)