An efficient and simple methodology for computer-based analysis of nematode swimming behavior in liquid is described. The method requires little to no investment for C. elegans laboratories. The hardware used is standard, and the computer software for the behavioral analysis (CeleST) is an open source one.
Dissecting the neuronal and neuromuscular circuits that regulate behavior remains a major challenge in biology. The nematode Caenorhabditis elegans has proven to be an invaluable model organism in helping to tackle this challenge, from inspiring technological approaches, building the human brain connectome, to actually shedding light on the specific molecular drivers of basic functional patterns. The bulk of the behavioral studies in C. elegans have been performed on solid substrates. In liquid, animals exhibit behavioral patterns that include movement at a range of speeds in 3D, as well as partial body movements, such as a posterior curl without anterior shape change, which introduce new challenges for quantitation. The steps of a simple procedure, and use of a software that enables high-resolution analysis of C. elegans swim behavior, are presented here. The software, named CeleST, uses a specialized computer program that tracks multiple animals simultaneously and provides novel measures of C. elegans locomotion in liquid (swimming). The measures are mostly grounded in animal posture and based on mathematics used in computer vision and pattern recognition, without computational requirements for threshold cut-offs. The software tool can be used to both assess overall swimming prowess in hundreds of animals from combined small batch trials and to reveal novel phenotypes even in well-characterized genetic mutants. The preparation of specimens for analysis with CeleST is simple and low-tech, enabling wide adaptation by the scientific community. Use of the computational approach described here should therefore contribute to the greater understanding of behavior and behavioral circuits in the C. elegans model.
Defining how genetics, epigenetics, experience, and environment influence behavior is a major challenge in modern biology. Simple, genetically amenable models that can be computationally tracked and measured can make critical contributions toward this goal. The nematode Caenorhabditis elegans is one such model. The purpose of this article is to demonstrate how C. elegans swimming locomotion can be tracked and quantitated to extract information on eight novel features with biological relevance.
C. elegans has a relatively short lifespan of about 2 – 3 w and a reproductive period of approximately 4 d at 20 °C1,2. In standard laboratory cultures, this microscopic nematode is grown on Petri plates containing Nematode Growth Media (NGM) that are spread with a bacterial food source3,4. WT N2 animals move actively in elegant sinusoidal waves on agar-filled plates; they change movement rates when roaming (food seeking), dwelling (browsing), or recovering from a meal (inactive satiety quiescence)5. Impairment6 or age7-12 can also dramatically alter movement on plates.
Genes that modulate behavior, healthspan or aging, can be functionally characterized by analyzing C. elegans movement patterns or locomotion. One approach to measure healthspan is to classify aging adults into three categories (e.g., A, B, and C) according to their locomotion on plates, with vigorous physical activity being class A and paralysis being class C7,8,13. Although such qualitative sorting is useful to reveal differences in physical fitness, the categories are broad without clear boundaries and thus their scoring is subject to experimenter bias.
A growing number of trackers have elevated the objectivity, sophistication, and precision of the analysis of C. elegans movement on solid media14-18. C. elegans locomotion on plates is mostly restricted to the plane in which the animal makes contact with the solid surface of the media. (C. elegans can also induce exploration by lifting the head away from the rest of the body that is still in contact with a solid surface, positioning the body in multiple planes. However, this behavior is unusual.) When placed in liquid, this nematode initiates an undulatory movement, or swimming, that features more extensive dimensional movement, a greater range of speed and depth of motion, and an increased incidence with age of anterior moving differently from posterior as compared to animals on solid surfaces. As a quick analysis of physical fitness and response to new environment, an experimenter can place an individual animal into a liquid drop and score its frequency of body bends under the dissecting scope. Video recording can facilitate present and future scorings of the swimming vigor of the animal. However, the manual approach limits the number of features that can be scored, and is fully constrained to scoring one animal at a time.
Locomotion in liquid has been less explored than locomotion on solid media. Indeed, there are few software options that are robust and easy to implement in the laboratory to measure locomotion in liquid19-24. The software CeleST (C. elegans Swim Test) offers simplicity of use and mathematical-based measures that deliver data (curvature scores) directly relevant to the nature of the motion8; (detailed description of features and advantages in Restif et al.8). Moreover, computational analysis enables the elucidation of phenotypic features that are impossible for the human eye to score. Here, data that exemplify the resolution of this analysis approach are presented and an easy to implement protocol to record swim trials for subsequent CeleST analysis is described.
1. Nematode Growth and Handling
2. Set up for Swim Analysis
3. Animal Preparation for Analysis of Physical Fitness in Liquid
4. Analysis of Physical Fitness in Liquid — Video Recording of Swimming
5. Analysis of Physical Fitness in Liquid — Video Treatment
6. Analysis of Physical Fitness in Liquid — Measurement
7. Analysis of Physical Fitness in Liquid — Output of Results
Note: The software can report on multiple measures of the swim motion, which cover from overt to very subtle behavioral traits not readily scored by eye (Videos 3 and 4). Here the focus is on 8 parameters that typically show a good dynamic range: Wave initiation rate, Body wave number, Asymmetry, Stretch, Curling, Travel speed, Brush stroke, and Activity index.
By analyzing locomotion in liquid (swimming), phenotypes that are not readily apparent on solid media (crawling) can be elucidated. To quantitate swimming locomotion we developed specific software that measures ten novel parameters of swimming behavior8. The eight most useful of these parameters are described in detail in Table 1. These parameters are named Wave initiation rate, Body wave number, Asymmetry, Stretch, Curling, Travel speed, Brush stroke, and Activity index. Studies exemplifying the power of the software have defined the functional decline of hundreds of aging adults with WT, behavioral or longevity mutant backgrounds8, and have analyzed the well-studied longevity mutants age-1(hx546) and daƒ-16(mgDƒ50), which harbor mutations that disrupt the normal insulin signaling pathway. The gene age-1 encodes for a phosphatidylinositide 3-kinase (PIK3) catalytic subunit, and, when it harbors the mutation hx546, causes lifespan extension and stress resistance27-29. The gene daƒ-16 encodes for a forkhead box O (FOXO) transcription factor that shortens lifespan and impairs the stress response when deleted30-33.
Certain parameters of swimming such as Wave initiation rate, Travel speed, Brush stroke and Activity index gradually declined with age even in favorable genetic backgrounds (Figure 1). In line with current knowledge, long-lived age-1(hx546) mutants showed a more vigorous physical performance than WT at advanced and extremely old ages. Also as anticipated, short-lived daƒ-16(mgDƒ50) mutants displayed compromised performance, especially at extremely old ages. Remarkably, it was only under the scrutiny of the CeleST computer vision and mathematical algorithm package that the superior swim performance of age-1(hx546) mutants was detectable at the onset of adulthood. The fact that age-1(hx546) results in enhanced physical performance at young adult life suggests that this mutation affects normal development and/or young adult phenotype in a way not previously appreciated (Figure 1).
Body wave number, Asymmetry, Stretch, and Curling parameters trended up with age in WT and aging mutant adults (Figure 2). Interestingly, the resolution level of the software revealed finer behavioral traits like the sustained symmetry of age-1(hx546) mutants throughout their lifespan and the inability of extreme old daƒ-16(mgDƒ50) mutants to stretch and curl up to the extent that same-age WT and age-1(hx546) adults do.
In addition to the overall inevitable loss of physical performance due to age, each individual adult displays a unique progression pattern through the aging process, even when genetics and environment are virtually homogeneous7. (By controlling genetics and environment, the possible confounding effects of these factors are minimized, unveiling the significant contribution of stochasticity to age-related degeneration.) A synchronized C. elegans population of similar genetic background kept in a controlled environment still contains a mix of different classes of individuals according to their aging peculiarities. Although all start as healthy adults, some rapidly lose their physical fitness (bad agers, class C) while others maintain vigor for longer period of time (graceful agers, class A). Bad agers thus appear to have a considerably shorter healthspan than graceful agers.
As further detailed in our study8, graceful agers maintained youthful physical fitness as observed by comparison with the swim profile of much younger adults (Figures 3, 4 and 5). This sustained fitness is comparable to the physical performance of long-lived age-1(hx546) mutants at post-reproductive age (D 11) (Figures 1 and 2). On the contrary, bad agers dramatically lost much of their physical capacity soon after reproduction, performing at levels similar to those of extreme old and progeric daƒ-16(mgDƒ50) adults (Figures 1 – 4). These similarities can be drawn by gross comparison, however particular signatures are readily noticeable at a closer look. For instance, although there is some correlation between the extent of both stretch and curling in extreme old wild types and aging mutants (Figure 2), this relation is not observed in bad agers (Figures 4 and 5), which show higher propensity to curl up but not to stretch in the samples tested. The software we present thus adds dimension to the analysis of physical fitness or locomotory capacity by providing the tools for more sophisticated studies that were not tractable before. In summary, CeleST provides comprehensive readouts in the form of eight novel measures highlighted here, which define the behavioral fingerprint of specific genetic, epigenetic, and environmental backgrounds, enabling the identification of unique and common parameter patterns that can be the signatures of specific conditions (environmental, pharmacological, nutritional), biological processes, or organism states such as healthspan.
Figure 1: CeleST Software Reports on Wave Initiation Rate (A), Activity Index (B), Brush Stroke (C) and Travel Speed (D) for WT, age-1(hx546), and daƒ-16(mgDƒ50) Adults on D 4 (young adulthood), 11 (post-reproduction) and 20 (extreme old). '#' in y axis means 'number'. WTs are colored in grey, age-1 in green and daƒ-16 in red. Error bars are the standard error of the mean (SEM). Same-age WT and aging mutants were compared for statistical significance using one-way ANOVA followed by Dunnett's multiple comparison test. **, p = 0.001 – <0.01; ***, p = 0.0001 – <0.001. n = 62 in each data point from four independent trials. Note that here, and for Figure 2, each individual 30 s video is made with 4 animals, and for each trial we score a total of 16 animals from 4 swim videos, this is done for 4 biological replicates for each data point shown. Please click here to view a larger version of this figure.
Figure 2: Software Reports on Body Wave Number (A), Asymmetry (B), Stretch (C), and Curling (D) for WT, age-1(hx546) and daƒ-16(mgDƒ50) Adults on D 4 (young adulthood), 11 (post-reproduction) and 20 (extreme old). '#' in y axis means 'number'. WT are colored in grey, age-1 in green, and daf-16 in red. Error bars are the standard error of the mean (SEM). Same-age WT and aging mutants were compared for statistical significance using one-way ANOVA followed by Dunnett's multiple comparison test. *, p = 0.01 – <0.05; **, p = 0.001 – <0.01; ***, p = 0.0001 – <0.001. n = 62 in each data point from four independent, 30 s swim trials. Please click here to view a larger version of this figure.
Figure 3: Software Reports on Wave Initiation Rate (A), Activity Index (B), Brush Stroke (C), and Travel Speed (D) for young WT adults (D 4), and Same-age Graceful and Bad Agers (D10 and 11). '#' in y axis means 'number'. Young WTs are colored in grey, class A graceful agers in green, and class C bad agers in red. Error bars are the standard error of the mean (SEM). Class A graceful agers and class C bad agers were compared to D 4 young adults using one-way ANOVA followed by Dunnett's multiple comparison test. ****, p <0.0001. n = 27 in each data point from two independent, 30 s swim trials. Graph is slightly modified from Restif et al. (2014)8, which was published under the Creative Commons Attribution (CC BY) license http://creativecommons.org/licenses/by/4.0/. Please click here to view a larger version of this figure.
Figure 4: Software Reports on Body Wave Number (A), Asymmetry (B), and Curling (C) for Young WT Adults (D 4), and Same-age graceful and Bad Agers (D 10 and 11). '#' in y axis means 'number'. Young WTs are colored in grey, class A graceful agers in green and class C bad agers in red. Error bars are the standard error of the mean (SEM). Class A graceful agers and class C bad agers were compared to D 4 young adults using one-way ANOVA followed by Dunnett's multiple comparison test. **, p = 0.001 – <0.01; ****, p <0.0001; n/a, non applicable since only one animal out of the total sample size curled. n = 27 in each data point from two independent, 30 s swim trials. Graph is slightly modified from Restif et al. (2014)8, which was published under the Creative Commons Attribution (CC BY) license http:/creativecommons.org/licenses/by/4.0/. Please click here to view a larger version of this figure.
Figure 5: Software Report on Stretch for Young WT adults (D 4), and Same-age Graceful and Bad Agers (D 10 and 11). Young WTs are colored in grey, class A graceful agers in green and class C bad agers in red. Error bars are the standard error of the mean (SEM). Class A graceful agers and class C bad agers were compared to D 4 young adults using one-way ANOVA followed by Dunnett's multiple comparison test. n = 27 in each data point from two independent trials. Graph is slightly modified from Restif et al. (2014)8, which was published under the Creative Commons Attribution (CC BY) license http:/creativecommons.org/licenses/by/4.0/. Please click here to view a larger version of this figure.
Video 1: Swimming of a Representative group of C. elegans adults. Please click here to view this video. (Right-click to download.)
Video 2: CeleST Software Computation of individual Curvature Maps of the Swim Performances of the tested Animals. Curvature maps are computed in the background; they do not appear in the software interface with the user. Please click here to view this video. (Right-click to download.)
Video 3: Software calculation of swim Measures based on Individual Curvature Maps. Please click here to view this video. (Right-click to download.)
Video 4: Software Calculation of Swim Measures that do not Rely on Curvature Maps. Please click here to view this video. (Right-click to download.)
The use of C. elegans as a model system continues to increase due to its genetic malleability, experimental tractability and anatomy annotated to meticulous detail. For instance, the neuronal structure and connectivity of the C. elegans hermaphrodite is clearly mapped out34-36, greatly facilitating investigations of specific neuronal circuits that control particular behaviors. 302 neurons constitute the adult hermaphrodite's nervous system, which process a broad range of sensory inputs into basic behavioral outputs like locomotion. The structure of the more complex male nervous system has also been described37, enabling sex-specific circuitry to be addressed.37
C. elegans behavior has been extensively studied on standard culture plates containing solid media. Since WT C. elegans moves in predictable sinusoidal waves on agar-filled plates, deviations from the overall pattern can be detected by eye and scored manually. This approach, however, is subject to the experimenter's criterion and is labor intensive. Hardware and software tools designed to track and measure C. elegans locomotion on solid media remove the subjectivity bias and allow for large-scale studies, permitting more sophisticated biological questions to be addressed. The recent behavioral database created by the Schafer lab16 is an excellent example of the analytical extension and depth that has been achieved with a computational system for locomotion on plates.
When WT C. elegans is placed in liquid, the animal quickly adapts its motion to the new environment, initiating a swim. C. elegans swimming utilizes a greater range of motion than crawling and can be more irregular8. Software like CeleST is intended to fill the gap for detailed analysis of C. elegans behavior in liquid, permitting quantitation of motion-associated parameters that are not readily measured by the un-aided eye, or that can be accomplished more rapidly than manual scoring. In 8 hard h an individual could process up to 200 videos, ~1,000 records per day.
The software defines swimming assessing parameters that serve as a comprehensive fingerprint of physical fitness and behavior. In addition to enriching understanding of the complex facets of C. elegans behavior in liquid and their underlying molecular pathways, this software can be used to explore multiple aspects of biology including pharmacological responses, aging, and distinct behavior. Presented here, the overview of the quantified changes that occur in the physical performance of C. elegans adults as they age is one example of such application of the software (for a more detailed account, see Restif et al.8). In the context of aging, some measured parameters declined while others increased or did not consistently change in wild type. Trends were confirmed to a large extent by the computational profile of longevity mutants, and the relative profiles of graceful and bad ager cohorts of same-age populations kept in uniform environmental conditions. The high resolution of the software can also reveal subtle phenotypes previously unknown in extensively characterized mutants (e.g., age-1(hx546) videodan Figure 1).
There are a few particularly critical steps of the protocol described. Maintaining a constant temperature environment between the swim environment and the strain plate culture is important for swim reproducibility, so experimenters are highly encouraged to go to great pains to avoid random temperature changes. Swim media should be at the same temperature as the plates. Likewise, careful attention to the drop size for the swim will help ensure reproducibility. Finally, it is prudent to think in advance about the offloading of the large video files that accumulate. Processing images on a site apart from the video capture computer is recommended.
The use of the software presented here for swim analysis has certain limitations. First, although the programs can simultaneously track multiple animals, if more than five animals are jointly analyzed, there is an increased chance that the animals will swim across each other in the video images. When the program cannot unambiguously determine which animal was which, it censors those data frames. Although this program feature ensures that data for individual animals are of high quality, it limits the throughput. Second, the images should be fairly clean, that is free of dust, smudges and glare from lights, as associated signals can also confound the image analysis. As noted in protocol section 2.1.1, a very low-tech investment that can greatly aid image capture by eliminating complications with fluctuations in environmental lighting is to cover the stage area with a dark cloth that does not allow ambient light to reach the stage. Third, the program is optimized for animals at adult stage. Young larvae swim very fast and have small bodies, which increases program error. Fourth, some of the software utilizes MATLAB, and when there are version upgrades and/or operating system upgrades, some program links may be disrupted. Currently, the software is optimized for use on MATLAB 2015b and Mac OS version 10.10, but we expect soon to post a software version that is more robust against such changes. Finally, the video data files can become large fast, and require storage space to be allotted.
In summary, presented here is a simple methodology that can be easily implemented by any laboratory without much investment to create videos of C. elegans swimming for CeleST analysis. Features of the software package include extensive automation from tracking through analysis, simultaneous multi-animal tracking, and use of mathematical bases (i.e., curvature measures) for quantitation of most locomotion parameters. The software is open source, with code and demos publicly available as detailed in Restif et al.8. Although the program features advanced computer vision analysis for tracking, other published tracking systems (e.g., Greenblum et al., 201438) are compatible with the parameter analysis of the software presented here. Future improvements will be directed toward converting the software into a more robust package that does not restrict use to the specific versions of the operating systems mentioned above (also indicated in the Table of Materials).
The authors have nothing to disclose.
CeleST development was supported by NIH grants R21AG027513 and U01AG045864. Data and some short video representations are adapted from Restif et al. (2014)8, which were published under the Creative Commons Attribution (CC BY) license http://creativecommons.org/licenses/by/4.0/. We thank Ricardo Laranjeiro for manuscript help.
REAGENT | |||
N2 | Caenorhabditis Genetics Center (CGC) | C. elegans wild type (ancestral). | |
OP50 Escherichia coli | Caenorhabditis Genetics Center (CGC) | Biosafety Level: BSL-1. | |
OP50-1 Escherichia Coli | Caenorhabditis Genetics Center (CGC) | Streptomycin resistant strain of OP50. Biosafety Level: BSL-1. | |
Streptomycin sulfate salt | Sigma-Aldrich | S6501 | |
Printed Microscope Slides | Thermo Fisher Scientific | Gold Seal Fluorescent Antibody Microslides: 3032-002 | have two etched 10mm diameter circles delineated by white ceramic ink |
Nematode Growth Medium (NGM) | For 1L: 17g Agar, 3g NaCl, 2.5g Peptone, 1mL 1M CaCl2, 1mL 5mg/mL Cholesterol in ethanol, 25mL 1M KPO4 buffer, 1mL 1M MgSO4, H2O to 1 L. Sterilize by autoclaving. | Stiernagle, T. Maintenance of C. elegans. WormBook, 1-11, doi:10.1895/wormbook.1.101.1 (2006) | |
M9 buffer | For 1L: 3g KH2PO4, 6g Na2HPO4, 5g NaCl, 1mL 1 M MgSO4, H2O to 1 L. Sterilize by autoclaving. | Stiernagle, T. Maintenance of C. elegans. WormBook, 1-11, doi:10.1895/wormbook.1.101.1 (2006) | |
EQUIPMENT | |||
CeleST | Driscoll Lab, Rutgers University | C. elegans Swim Test | Open Source, see http://celestmod.github.io/CeleST/ and http://celest.mbb.rutgers.edu/ |
MATLAB | www.mathworks.com/downloads | MatLab version 2015b (best) | The CeleST version demonstrated here has best functionality with Mac OS 10.10 and MatLab 2015b. MATLAB 2015B introduced changes to how MATLAB handled graphics, including a new coding convention and syntax. These changes resolved an issue that couldn't be resolved elegantly (primarily because the internals of MATLAB really needed the major graphics overhaul implemented in MATLAB 2015B). For this reason, CeleST should always be run on MATLAB 2015B or later versions. However for users without access to MATLAB version 2015B or later (or MATLAB at all), we have created a CeleST program that doesn't need MATLAB on the computer at al. An installer is downloaded by the prospective user and then it installs itself onto the computer through a couple prompts like most programs. |
Mac OS | www.apple.com | Version 10.10 | Currently, CeleST has been ported to the major operating systems (Windows, Mac, and Linux). The current code can be run on any of the operating systems and there are versions for each operating system that don't even require users to have MATLAB to use CeleST (this version requires a large download). The Windows version has been tested the least and is most prone to bugs as such. Linux has been moderately tested. And Mac has been and continues to be tested extensively (primarily because it's the operating system in our lab). |
Stereomicroscope | Zeiss | Stemi 2000-C | |
Transmitted Light Base | Diagnostic Instruments | TLB 3.1 | |
Digital CCD Camera | QImaging | Rolera-XR Mono Fast 1394 (ROL-XR-F-M-12) | |
Digital Video Recording Software | Norpix | Streampix Version 3.17.2 |