Skip to main content
  • Research article
  • Open access
  • Published:

A novel application of motion analysis for detecting stress responses in embryos at different stages of development



Motion analysis is one of the tools available to biologists to extract biologically relevant information from image datasets and has been applied to a diverse range of organisms. The application of motion analysis during early development presents a challenge, as embryos often exhibit complex, subtle and diverse movement patterns. A method of motion analysis able to holistically quantify complex embryonic movements could be a powerful tool for fields such as toxicology and developmental biology to investigate whole organism stress responses. Here we assessed whether motion analysis could be used to distinguish the effects of stressors on three early developmental stages of each of three species: (i) the zebrafish Danio rerio (stages 19 h, 21.5 h and 33 h exposed to 1.5% ethanol and a salinity of 5); (ii) the African clawed toad Xenopus laevis (stages 24, 32 and 34 exposed to a salinity of 20); and iii) the pond snail Radix balthica (stages E3, E4, E6, E9 and E11 exposed to salinities of 5, 10 and 15). Image sequences were analysed using Sparse Optic Flow and the resultant frame-to-frame motion parameters were analysed using Discrete Fourier Transform to quantify the distribution of energy at different frequencies. This spectral frequency dataset was then used to construct a Bray-Curtis similarity matrix and differences in movement patterns between embryos in this matrix were tested for using ANOSIM.


Spectral frequency analysis of these motion parameters was able to distinguish stage-specific effects of environmental stressors in most cases, including Xenopus laevis at stages 24, 32 and 34 exposed to a salinity of 20, Danio rerio at 33 hpf exposed to 1.5% ethanol, and Radix balthica at stages E4, E9 and E11 exposed to salinities of 5, 10 and 15. This technique was better able to distinguish embryos exposed to stressors than analysis of manual quantification of movement and within species distinguished most of the developmental stages studied in the control treatments.


This innovative use of motion analysis incorporates data quantifying embryonic movements at a range of frequencies and so provides an holistic analysis of an embryo’s movement patterns. This technique has potential applications for quantifying embryonic responses to environmental stressors such as exposure to pharmaceuticals or pollutants, and also as an automated tool for developmental staging of embryos.


The development of imaging technology over the past two decades has enabled significant advances in biology and resulted in an exponential increase in the size and complexity of the data sets generated [1, 2]. However, this ability to acquire data at a rapid rate has also led to major data management challenges [3]. Consequently the ability to manage and analyse image sequences, particularly those recorded at a high frame rate, is now recognised as the major bottleneck in laboratories using high throughput imaging systems [4, 5]. This bottleneck is being addressed by the development of bioinformatics, an emerging field tasked with developing techniques to extract biological data from image datasets allowing the subsequent deletion of memory intensive images [6].

Among the many methods of automated image analysis available to biologists is motion analysis, a technique that has been used for quantifying the behavior of mammals [7], fish [8], sperm [9], larvae [10], algal spores [11, 12] and bacteria [13]. These applications measure parameters including distance travelled, speed and angular velocity and they therefore rely on an observable forward motion in the sample studied. The application of these motion analysis systems to quantify embryonic movement effectively is limited, as the movement patterns embryos exhibit are often complex, subtle and diverse (e.g. muscle flexing, spinning, heart beat, blood flow and tail flicking), which creates challenges in quantifying embryonic movements in a holistic way. Here we demonstrate a technique for quantifying motion patterns in embryos that can distinguish the effects of environmental stressors on embryos at different developmental stages.

We applied sparse optic flow [14], a machine vision motion analysis technique (Figure 1), to image sequences of two widely used vertebrate model species at three early developmental stages: (i) Danio rerio embryos at stages 19 h, 21.5 h and 33 h [15], and (ii) Xenopus laevis at stages 24, 32 and 34 [16]; and (iii) to embryos of an invertebrate, the freshwater pond snail Radix balthica, at stages E3, E4, E6, E9 and E11 [17]. At each developmental stage image sequences were recorded under control conditions and in embryos exposed to abiotic stressors, for which developmental effects have previously been observed [18-26]. Danio rerio embryos were exposed to 1.5% ethanol [18 - 21], and a salinity of 5 [22-24], Xenopus laevis were exposed to a salinity of 20 [25] and Radix balthica embryos were exposed to salinities of 5, 10 and 15 [26].

Figure 1
figure 1

Optic flow of an E6 stage Radix balthica embryo showing angular rotations (red and yellow) and centre of mass (blue line).

As the embryos are shape-deformable it is sensible to encode optic flow as a summation of directed rotations. Optic flow was therefore calculated as the mean angles of positive (clockwise) and negative (anti-clockwise) rotation for each frame and, from these data, the centre of mass was calculated as X and Y polar coordinates (Figure 1). The patterns in these frame-by-frame movement parameters (Figure 2) were analysed using Discrete Fourier analysis, run using the R language for statistical computing [27], to quantify the distribution of energy of each of these parameters, within 18 to 30 frequency bins ranging between 600 - 0.1 s. A Bray-Curtis similarity matrix based on 72 variables (18 frequency bins for each of the parameters - mean angles of positive and negative rotation and the X and Y coordinates of the calculated centre of mass) was visualised using Multidimensional Scaling (Figure 3) and differences in movement patterns between developmental stages and treatment groups were tested statistically using ANOSIM [28] (Tables 1 and 2). This analysis revealed that, in all three species, for most of the developmental stages examined, individuals exposed to an environmental stressor were distinguishable from control organisms and also that many of the developmental stages in control treatments were distinguishable from each other. This technique could be applied to any species whose movement patterns can be visualized and could have applications within the fields of biotechnology, toxicology and developmental biology where the extraction of a holistic measure of embryonic movement patterns would be desirable.

Figure 2
figure 2

Frame by frame optic flow parameters of control and treatment embryos for each of the developmental stages of Danio rerio , Xenopus laevis and Radix balthica . Red - positive angle movements, blue - negative angle movements, green - X coordinate of centre of mass, yellow - Y coordinate of centre of mass. A time period (Danio rerio - 10 min, Xenopus laevis - 10 min, Radix balthica - 5 min) of a single embryo is shown, exhibiting typical movement patterns for each treatment response per developmental stage.

Figure 3
figure 3

Multidimensional scaling plots in two dimensions using Bray-Curtis similarity matrices. MDS plots using Bray-Curtis similarity matrices performed on logarithmically transformed spectral frequency data produced using Discrete Fourier Transform analysis of (i) negative angle (ii) positive angle (iii) centre of mass - rho and (iv) centre of mass - theta, frame-to-frame parameters for (a) Danio rerio, (b) Xenopus laevis and (c) Radix balthica under different environmental conditions and at different developmental stages. Bubble size represents the amount of tail flicks of Danio rerio and Xenopus laevis, the number of rotations of Radix balthica at stages E3, E4 and E6, and the number of complete gliding along the circumference of the egg capsule on stages E9 and E11.

Table 1 Comparison of analysis of treatment groups by both the automated motion analysis technique and observer quantification
Table 2 Significance levels, calculated using ANOSIM, of pairwise differences in developmental stages


Time series of the output parameters produced by optic flow (Figure 2) revealed patterns that accurately reflected the organism’s movement seen during manual observation of the same image sequence, e.g. tail flicks in zebrafish embryos resulted in peaks in the motion parameters recorded by the optic flow motion analysis. Subsequent analysis of the frame-to-frame motion analysis parameters was successful in distinguishing individuals exposed to environmental stress from those under control conditions within many of the developmental stages studied (Table 1), and to distinguish, within species, individuals under control conditions at different stages of development (Table 2).

Danio rerio

Danio rerio embryos at 33 hpf exposed to 1.5% ethanol were statistically distinguishable from those under control conditions (Table 1). All developmental stages of Danio rerio used here (19, 21.5, 33 hpf) exhibited increased tail flick frequency in response to exposure to 1.5% ethanol and this was visible in the optic flow motion parameter plots as increased peak frequencies (Figure 2). Embryos exposed to 1.5% ethanol also exhibited more vigorous tail flicking and this caused greater levels of movement within the egg than seen in embryos under control conditions, recorded by optic flow as changes to the calculated X and Y coordinates of the centre of mass. However, embryonic movement responses at stages 19 and 21.5 hpf were highly variable. Two embryos from both 19 and 21.5 hpf showed movement patterns typical of control embryos and one embryo at 19.5 hpf exposed to ethanol performed more tail flicking (182 tail flicks in 10 min) than all other embryos exposed to ethanol combined. This inter-individual variation perhaps goes someway to explaining why stress responses in these developmental stages were not distinguishable. Exposure to a salinity of 5 produced no discernible difference in movement responses compared with control embryos.

Danio rerio at 19, 21.5 and 33 h post fertilization (hpf) under control conditions were all distinguishable from each other via analysis of frame-to-frame parameters. Danio rerio at 19 h post fertilisation (hpf) performed a regular flicking at the end of the tail. This movement resulted in optic flow producing small peaks in the magnitude of clockwise and anti-clockwise movements and in the calculated X and Y coordinates of the embryo’s centre of mass. Embryos at 21.5 hpf performed a less regular tail flicking involving a greater length of the embryo’s tail and this was visible in the motion analysis data as peaks in clockwise and anticlockwise movements and also in the X and Y image coordinates of the embryo, with greater amplitude than those seen in 19 hpf embryos. The more vigorous embryonic movements caused occasional shifts in the orientation of the embryo within the egg and this resulted in periodic shifts in the X and Y image coordinates of the embryo. Danio rerio at 33 hpf exhibited the longest periods with no gross movement, but when embryos at this stage did perform tail flicks these were more vigorous than at 19 or 21.5 hpf and so this produced larger peaks in the clockwise and anticlockwise movements recorded. There was little change in the X and Y coordinates of the centre of mass in Danio rerio at 33 hpf due to the embryo having little space to move within the egg capsule.

Xenopus laevis

Xenopus laevis exposed to a salinity of 20 in the three developmental stages used here were statistically distinguishable from control individuals (Table 1; Figure 3). Xenopus laevis at stage 24 exhibited an increase in tail flicking frequency when exposed to a salinity of 20, although there was significant inter-individual variation in this response ranging from 48 to 400 tail flicks, during the ten minute recording period. Despite this variation there was a clear distinction in the Multidimensional Scaling plot between control and salinity exposed embryos (Figure 3). Xenopus laevis at both stages 32 and 34 exhibited more frequent tail flicking when exposed to the salinity treatment. Tail flicking in stage 32 embryos exposed to salinity caused movement of the egg capsule within the well and this was observed in the motion analysis as changes to the calculated X and Y centre of mass parameters.

Xenopus laevis under control conditions at stage 34 was distinguishable from both stages 32 and 24, however stage 24 was not distinguishable from 32. Xenopus laevis at stage 24 performed only a very occasional tail flicking and this was a slow movement, relative to stages 32 and 34, leading to a gradual peak in either the clockwise or anti-clockwise magnitude of movement, depending on the direction of the tail flick. Embryos at stages 32 and 34 had greater regularity in the frequency of tail flicking than embryos at stage 24 and these tail flicks were often too rapid for optic flow to track continuously in the image sequence resulting in sharp peaks in the negative and positive movements calculated. In stage 32 and 34 individuals optic flow also identified and tracked blood flow within the organism and this produced a low level of continuous negative and positive movements which were not present in the motion parameter data for stage 24 embryos in which blood flow was not visible. Stage 34 individuals had hatched from the egg capsule but motion analysis output did not show any difference between them and stage 32 embryos, as the organism displayed the same body flicking behavior. The main difference in motion parameters between stages 32 and 34 was caused by a more pronounced heart beat and blood flow in stage 34, producing more continuous low level negative and positive movements (Figure 2).

Radix balthica

Embryos at all five developmental stages exposed to salinity treatments of 5, 10 or 15, with the exception of embryos at developmental stage E9 exposed to a salinity of 5, were statistically distinguishable from control embryos (Table 1; Figure 3). All developmental stages of Radix balthica used here exhibited a dose-dependent decreased frequency and magnitude of embryonic movements with increased salinity. Embryos from developmental stages E3, E4 and E6 exhibited fewer ciliary-driven oscillations than control embryos. Embryos at E3, which under control conditions undergo a “figure of eight” path within the egg, when exposed to salinity, rotated with little net movement within the egg. Radix balthica at stages E9 and E11 exposed to a salinity of 5 exhibited less gliding around the egg compared to control embryos, and this resulted in less positive and negative movement identified by optic flow as less change to the X and Y centre of mass motion parameters. When embryos at stages E9 and E11 were exposed to a salinity of either 10 or 15 they stopped gliding around the egg, using their cilia, resulting in only minor changes to X and Y centre of mass coordinates. However, despite the lack of overall movement in the egg, embryos in these treatments did continue to exhibit small scale movements of their muscular foot and shell which led to small peaks in the negative and positive movement parameters recorded by the motion analysis.

Radix balthica embryos under control conditions at the developmental stages; E3, E4, E6, E9 and E11 were all statistically distinguishable from each other with the exception of E9 from both E3 and E6 (Figure 3; Table 2). The movement patterns of R. balthica produced by optic flow were markedly different between early and late developmental stages due to the developmental transition from regular ciliary driven rotations to complex muscular movements more similar to those seen in adult snails (Figure 2). The manually observed rates of ciliary rotation in the three earliest stages (E3, E4 and E6) closely matched the cyclical patterns in negative and positive angles of rotation produced by optic flow.

The oscillatory pattern in clockwise and anti-clockwise movement became more variable as development progressed from ciliary rotation (E3, E4 and E6) to muscular movement (E9 and E11) producing marked differences in the frame-to-frame optic flow parameters. Optic flow was still able to accurately to track embryos exhibiting these more complex muscular movement behaviors, despite accurate manual quantification not being possible.

E3 stage embryos are spherical and exhibited a “figure of eight” type movement within the egg capsule, which produced a sinusoidal pattern in the X and Y coordinates of its centre of mass (Figure 2). E4 stage embryos are “crescent shaped” and occupy approximately three times the volume of the E4 stage, which resulted in the coordinates of the centre of mass having a frequency close to that of its rotational rate. There was some long term movement in the X and Y coordinates of the calculated centre of mass caused by lateral, as opposed to spinning, movement within the egg capsule. E6 stage embryos occupy approximately one third of the egg capsule and are significantly more asymmetrical than the E4 stage and also exhibit variable rates of spinning which leads to a less uniform pattern. In E9 stage embryos, cilia-driven rotational behaviour is largely absent and replaced by prolonged periods of apparent inactivity, during which the embryos undergo muscular contractions and flexing of the foot and shell. The loss of rotational movement is evidenced by the lack of a clear rhythm in the centre of mass time series replaced instead by a mostly linear pattern. Periods during which no large scale movements occur are interrupted by occasional gliding around the egg capsule. E11 stage embryos are much larger than E9 stage embryos, the former almost completely filling the egg capsule, however they exhibit similar patterns of behaviour to E9 embryos. This increase in size explains the decrease in the amount of “noise” surrounding the X and Y centre of mass from the E9 to E11 stage as the embryo has very little room to move within the egg capsule in the E11 stage which likely restricts small scale movement.

Manual quantification of embryonic behaviour (tail flick frequency for Danio rerio and Xenopus laevis embryos and spinning frequency for Radix balthica) was found to be less effective at distinguishing individuals exposed to stress treatments than the motion, and subsequent frequency, analysis presented here. Manual quantification of movement parameters could not distinguish treatments that were not distinguishable using this motion analysis technique. Furthermore, this technique was successfully able to distinguish some treatments that were not distinguishable using manual analysis (Table 1).


The application of motion analysis that we present here to detect the effects of abiotic stressors on different stages of embryos of three species incorporates data quantifying embryonic movements at a range of frequencies and so provides a holistic analysis of embryonic motion. This type of analytical approach may be applicable to a wide range of species, at various developmental stages, as it doesn’t rely on one particular type of movement, but rather being able to incorporate a diverse range of embryonic movements into a single analysis. This potential benefit is demonstrated here through the ability of this technique to distinguish environment-induced changes in motion in early molluscan developmental stages undergoing ciliary rotation, late molluscan development, during which the heart beat and slow muscular contractions are the prominent movement, and vertebrate developmental stages during which movement is mostly in the form of tail flicking or, in later development, the beating heart and associated blood flow. This technique was also better able to distinguish embryos exposed to stressors than analysis of manual quantification of movement patterns (Table 1).

While we have shown that this technique can be applied to embryos of two of the most common vertebrate model species and an invertebrate embryo, there is no reason why, in principle, such analysis cannot be extended to species belonging to other phyla (e.g. crustaceans and annelids). The spectral frequency analysis used here to analyse the frame-to-frame motion parameters, could be modified to target particular frequencies and refined for use on different species or to increase efficacy at detecting a particular behavioural response to a stressor. Refinement of the spectral analysis to target frequencies relevant for particular movements would enable this technique to be used as an automated toxicological assay able to quantify complex behavioral responses to stressors.

Image analysis and automated pattern recognition have been used to relatively good effect in reliable and quick species identification procedures [29-32] and the determination of developmental stages [33, 34], but have not been used successfully to stage live animals in real time. The innovative use of motion analysis presented here was successfully able to distinguish many of the developmental stages studied and so demonstrates its potential application to an automated system for staging living, animal embryos. Although most biologists are probably accustomed to staging organisms prior to their use in an experiment the efficacy and reliability of developmental staging has long been questioned [35]. Automated image analysis is recognized as offering advantages over manual analysis; including greater consistency and accuracy and is less time consuming [32]. Therefore this approach could offer the benefits of automated image analysis for staging living embryos. The ability of this technique to distinguish stress induced changes and embryonic developmental stages also demonstrates its potential application to automate both, screening of short term stress responses and longer term effects on rates of development e.g. decreased developmental rate in long term exposure to a stressful environment.


Image analysis is a key tool in biology for the extraction of biologically meaningful data from image datasets. Here we present a novel application of motion analysis coupled with spectral frequency analysis for providing a holistic measure of movement patterns. This technique is able to distinguish environmentally-induced alterations to movement patterns, and also individuals at different stages of development, and therefore could have applications in disciplines such as ecotoxicology, developmental biology and biotechnology.



Danio rerio (obtained from a stock population at Plymouth University, UK) were cultured in glass beakers (vol. = 200 ml), containing aerated, deionized water at 27°C. Periodically water in the glass beakers was changed at which time dead embryos were removed. Embryos at 19, 21.5 and 33 h post fertilization (hpf) [15] were placed into 20 ml glass vials and exposed to either: 1.5% ethanol (1.5% EtOH (purity = 99.99%, Fisher, UK) in deionized water) for 30 minutes, a salinity of 5 (Instant Ocean® Salt in deionized water) for 60 minutes or to a control treatment of deionized water. Glass vials were sealed with microfilm to prevent losses by volatilization and maintained at 27°C. Embryos were transferred to 96 multiwell plates for image acquisition.

Xenopus laevis (obtained from Bristol University, UK) were cultured in glass beakers (vol. = 200 ml) containing aerated, deionized water at 23°C. Twenty minutes prior to image sequence acquisition embryos at stages 24, 32 and 34 [16] were transferred to 20 ml glass vials and exposed to either: a salinity of 20 or a control of deionized water and maintained at 23°C. Embryos were transferred to 96-multiwell plates prior to image acquisition.

Radix balthica embryos (obtained from a stock population at Plymouth University) at developmental stages E3, E4, E6, E9 and E11 were placed individually in multiwell plates (96 wells, volume = 0.3 ml/well) containing a control of Artificial Pond Water (APW) [36] or a salinity of 5, 10 or 15 (Instant Ocean® salt in APW). Embryos were exposed to treatments 10 minutes prior to image sequence acquisition.

Imaging setup

Images from a 4 MP monochrome camera (Allied Vision Technology Pike 421B) operating with 1600 × 1200 pixels at 15 Hz for the Danio rerio and Xenopus laevis image sequence acquisition and with 1000 × 1000 pixels at 7.5 Hz for Radix balthica. The camera was connected to zooming optic lens systems (Danio rerio - Keyence VHZ100R; Xenopus laevis - Keyence VHZ20R; Radix balthica - Zoom 70 XL Optical System), which was inverted beneath a Prior Optiscan motorized XY stage, using a purpose-built frame, and controlled using Micromanager software package [37]. Ten minute long image sequences were recorded for six individuals from each of: Danio rerio at 19 h, 21.5 h and 33 h stages from a deionized water control, 1.5% ethanol and a salinity of 5 [15]; Xenopus laevis at stages 24, 32 and 34 in a control treatment of deionized water and a salinity of 5 [16]; Radix balthica at stages E3, E4, E6, E9, E11 in a control treatment of Artificial Pond Water [36] and salinities of 5, 10 and 15 [17] (Figure 3).

Motion analysis

Off-line, the image sequence was analysed for optical flow using the OpenCV toolkit by applying the sparse optic flow algorithm [14]. This temporal tracking technique first extracts corner features from each image. It then attempts to register image features from frame to frame, and assigns a velocity variance for each. A feature with high variance is probably not being tracked accurately and may be rejected. All remaining features are candidates for measuring an embryo’s motion. To ensure optic flow was accurately tracking embryonic movement all image sequences were watched whilst the analysis was running.

Four measures of motion were extracted for each frame (by comparison with previous frames). These are positive and negative angles of rotation and rho-theta angular changes (in polar coordinates) to the centre of mass for the features tracked within the embryo (see Figure 1 for example). Note the yellow and red vectors are negative and positive rotations (feasible because the embryos are gelatinous and flexible) and the blue vector arrow is the frame to frame change to the centre of mass.

The method is not specific to the embryo and can be misled by the motion of other targets, such as the embryo’s egg capsule. However, this is normally infrequent. Frame-to-frame measurements from the motion analysis were saved to a CSV format file.

Spectral frequency analysis

Frame-to-frame motion analysis data were analysed for spectral content using the Discrete Fourier Transform (DFT) [38]. This reveals cyclic frequency of motion information for each embryo. This is completed using the R language for statistical computing [27]. Data are collected over a 10 min period and so frequency analysis can reveal repetitive oscillations from 3.75 Hz to 10 min for Radix balthica (recorded at 7.5 fps) and 7.5 Hz to 10 min for Danio rerio and Xenopus laevis (recorded at 15 fps). These are calculated for each embryo image sequence. Specific wavelength bins were further analysed from the DFT. For Danio rerio and Xenopus laevis the following frequency bins were used:- [1]: 598 - 298 s; [2]: 200 - 149 s; [3]: 120 - 99 s; [4]: 85 - 54 s; [5]: 49 - 33 s; [6]: 29 - 20.5 s; [7]: 21 - 17 s; [8]: 16 - 15 s; [9]: 14 - 13 s; [10]: 12 - 11 s; [11]: 10 - 5 s; [12]: 4.9 - 3.5 s; [13]: 3.3 - 2.5 s; [14]: 2.4 - 1.7 s; [15]: 1.6 - 1 s; [16]: 0.9 - 0.33 s; [17]: 0.3 - 0.25 s; [18]: 0.24 - 0.1 s. Analysis of Radix balthica’s frame-to-frame motion analysis output was performed using 18 frequency bins for developmental stage comparisons - [1]: 598 - 299 s; [2]: 291 - 145 s; [3]: 120 - 100 s; [4]: 85 - 75 s; [5]: 66 - 60 s; [6]: 54 - 50 s; [7]: 46 - 43 s; [8]: 40 - 37 s; [9]: 35 - 30 s; [10]: 29 - 24 s; [11]: 20 - 17 s; [12]: 15 - 12.5 s; [13]: 12.2 - 10 s; [14]: 9.5 - 5 s; [15]: 2 - 1.1 s; [16]: 1 - 0.68 s; [17]: 0.66 - 0.5 s; [18]: 0.3 - 0.26 s, and 30 frequency bins for environmental stressor analysis - [1]: 300 - 150 s; [2]: 100 - 75 s; [3]: 60 - 50 s; [4]: 43 - 37.5 s; [5]: 33.5 - 30 s; [6]: 27.5 - 25 s; [7]: 23 - 21.5 s; [8]: 20 - 19 s; [9]: 17.5 - 16.5 s; [10]: 16 - 14.5 s; [11]: 13.5 - 12.5 s; [12]: 12 - 11 s; [13]: 10.5 - 9.5 s; [14]: 9 - 6.5 s; [15]: 6.3 - 5 s; [16]: 4.9 - 3.5 s; [17]: 3.3 - 2.5 s; [18]: 2.4 - 2 s; [19]: 1.9 - 1.5 s; [20]: 1.4 - 1 s; [21]: 0.9 - 0.8 s; [22]: 0.7 - 0.65 s; [23]: 0.6 - 0.5 s; [24]: 0.49 - 0.4 s; [25]: 0.39 - 0.34 s; [26]: 0.33 - 0.28 s; [27]: 0.27 - 0.22 s; [28]: 0.21 - 0.18 s; [29]: 0.17 - 0.15 s; [30]: 0.14 - 0.13 s. Radix balthica during its development performs movement types ranging from ciliary driven rotation to muscular driven crawling. As a result, to discern treatments more frequency bins were required for analysis of Radix balthica than for Danio rerio or Xenopus laevis, as movement differences fell within rather than between frequency bins. These extra frequency bins enhanced the frequency discrimination and therefore the detection of subtle changes to these complex movement behaviours.

The data resulting from the DFT comprises 72 parameters (18 frequency bins for each of positive and negative angles of rotation and rho-theta angular changes, in polar coordinates, to the centre of mass) for each individual. These data were transformed (Log X + 1) and for each species a Bray-Curtis similarity matrix was calculated. This matrix was used to generate Multidimensional Scaling (MDS) plots and ANOSIM [28] was used to test for the degree of similarity within groups of embryos from particular developmental stages (Primer-E).


  1. O'Donoghue SI, Gavin A, Gehlenborg N, Goodsell DS, Hériché J, Nielsen CB, North C, Olson AJ, Procter JB, Shattuck DW, Walter T, Wong B: Visualizing biological data—now and in the future. Nat Methods 2010, 7: S2-S4. 10.1038/nmeth.f.301

    Article  PubMed  Google Scholar 

  2. Walter T, Shattuck DW, Baldock R, Bastin ME, Carpenter AE, Duce S, Ellenberg J, Fraser A, Hamilton N, Pieper S, Ragan MA, Schneider JE, Tomancak P, Hériché J: Visualization of image data from cells to organisms. Nat Methods 2010, 7: S26-S41. 10.1038/nmeth.1431

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Swedlow JR, Goldberg IG, Eliceiri KW: Bioimage informatics for experimental biology. Annu Rev Biophys 2009, 38: 327-346. 10.1146/annurev.biophys.050708.133641

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Eils R, Athale C: Computational imaging in cell biology. J Cell Biol 2003, 161: 477-481. 10.1083/jcb.200302097

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Zhou X, Wong STC: Informatics challenges of high-throughput microscopy. IEEE Signal Proc Mag 2006, 23: 63-72.

    Article  Google Scholar 

  6. Peng H: Bioimage informatics: a new area of engineering biology. Bioinformatics 2008, 24: 1827-1836. 10.1093/bioinformatics/btn346

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Spink AJ, Tegelenbosch RAJ, Buma MOS, Noldus L: The EthoVision video tracking system - a tool for behavioral phenotyping of transgenic mice. Physiol Behav 2001,73(5):731-744. 10.1016/S0031-9384(01)00530-3

    Article  CAS  PubMed  Google Scholar 

  8. Williams LR, Wong K, Stewart A, Suciu C, Gaikwad S, Wu N, DiLeo J, Grossman L, Cachat J, Hart P, Kalueff AV: Behavioral and physiological effects of RDX on adult zebrafish. Comp Biochem Phys C 2012, 155: 33-38.

    CAS  Google Scholar 

  9. Liu QH, Li J, Xiao ZZ, Ding FH, Yu DD, Xu XZ: Use of computer-assisted sperm analysis (CASA) to evaluate the quality of cryopreserved sperm in red seabream (Pagrus major). Aquaculture 2007, 263: 20-25. 10.1016/j.aquaculture.2006.11.017

    Article  Google Scholar 

  10. Amsler MO, Amsler CD, Rittschof D, Becerro MA, McClintock JB: The use of computer-assisted motion analysis for quantitative studies of the behaviour of barnacle (Balanus amphitrite) larvae. Mar Freshw Behav Physiol 2006,39(4):259-268. 10.1080/10236240600980640

    Article  Google Scholar 

  11. Greer SP, Iken K, McClintock JB, Amsler CD: Bioassay-guided fractionation of antifouling compounds using computer-assisted motion analysis of brown algal spore swimming. Biofouling 2006,22(2):125-132. 10.1080/08927010600602082

    Article  CAS  PubMed  Google Scholar 

  12. Iken K, Greer SP, Amsler CD, McClintock JB: A new antifouling bioassay monitoring brown algal spore swimming behaviour in the presence of echinoderm extracts. Biofouling 2003,19(5):327-334. 10.1080/08927010310001612045

    Article  PubMed  Google Scholar 

  13. Lopez-de-victoria G, Zimmerfaust RK, Lovell CR: Computer-assisted video motion analysis - a powerful technique for investigating motility and chemotaxis. J Microbiol Methods 1995,23(3):329-341. 10.1016/0167-7012(95)00047-O

    Article  Google Scholar 

  14. Lucas BD, Kanade T: An iterative image registration technique with an application to stereo vision. In Proceedings of the 7th International Joint Conference on Artificial Intelligence (IJCAI). San Francisco, USA; 1981:121-130.

    Google Scholar 

  15. Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF: Stages of embryonic development of the zebrafish. Dev Dynam 1995,203(3):253-310. 10.1002/aja.1002030302

    Article  CAS  Google Scholar 

  16. Nieuwkoop PD, Faber J: Normal Table of Xenopus laevis (Daudin):A systematical and chronological survey of the development from the fertilized egg till the end of metamorphosis. Amsterdam: North-Holland Pub Co; 1956.

    Google Scholar 

  17. Cumin R: Normentafel zur Organogenese von Limnaea Stagnalis (Gastropoda, Pulmonata) mit besonderer Berücksichtigung der Mitteldarmdrüse. Rev Suisse Zool 1972, 79: 709-774.

    Google Scholar 

  18. Ali S, Champagne DL, Alia A, Richardson MK: Large-scale analysis of acute ethanol exposure in zebrafish development: a critical time window and resilience. PLoS One 2011,6(5):e20037. 10.1371/journal.pone.0020037

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Bilotta J, Saszik S, Givin CM, Hardesty HR, Sutherland SE: Effects of embryonic exposure to ethanol on zebrafish visual function. Neurotoxicol Teratol 2002,24(6):759-766. 10.1016/S0892-0362(02)00319-7

    Article  CAS  PubMed  Google Scholar 

  20. Lockwood B, Bjerke S, Kobayashi K, Guo S: Acute effects of alcohol on larval zebrafish: a genetic system for large-scale screening. Pharmacol Biochem Behav 2004,77(3):647-654. 10.1016/j.pbb.2004.01.003

    Article  CAS  PubMed  Google Scholar 

  21. MacPhail RC, Brooks J, Hunter DL, Padnos B, Irons TD, Padilla S: Locomotion in larval zebrafish: influence of time of day, lighting and ethanol. Neurotoxicology 2009,30(1):52-58. 10.1016/j.neuro.2008.09.011

    Article  CAS  PubMed  Google Scholar 

  22. Sawant MS, Zhang S, Li L: Effect of salinity on development of zebrafish, Brachydanio rerio. Curr Sci 2001,81(10):1347-1350.

    Google Scholar 

  23. Boisen AMZ, Amstrup J, Novak I, Grosell M: Sodium and chloride transport in zebrafish soft water and hard water acclimated (Danio rerio). Biochim Biophys Acta-Biomembr 2003,1618(2):207-218. 10.1016/j.bbamem.2003.08.016

    Article  CAS  Google Scholar 

  24. Craig PM, Wood CM, McClelland GB: Gill membrane remodeling with soft-water acclimation in zebrafish (Danio rerio). Physiol Genomics 2007,30(1):53-60. 10.1152/physiolgenomics.00195.2006

    Article  CAS  PubMed  Google Scholar 

  25. Karraker NE, Arrigoni J, Dudgeon D: Effects of increased salinity and an introduced predator on lowland amphibians in Southern China: Species identity matters. Biol Conserv 2010,143(5):1079-1086. 10.1016/j.biocon.2010.01.020

    Article  Google Scholar 

  26. Tills O, Spicer JI, Rundle SD: Salinity-induced heterokairy in an upper-estuarine population of the snail Radix balthica (Mollusca: Pulmonata). Aquat Biol 2010,9(1):95-105.

    Article  Google Scholar 

  27. R Development Core Team: R: A language and environment for statistical computing.

  28. Clarke KR: Non-parametric multivariate analyses of changes in community structure. Aust J Ecol 1993, 18: 117-143. 10.1111/j.1442-9993.1993.tb00438.x

    Article  Google Scholar 

  29. Yu DS, Kokko EG, Barron JR, Schaalje GB, Gowen BE: Identification of ichneumonid wasps using image analysis of wings. Syst Entomol 1992, 17: 389-395. 10.1111/j.1365-3113.1992.tb00558.x

    Article  Google Scholar 

  30. Dörge T, Carstensen JM, Frisvad JC: Direct identification of pure Penicilium species using image anaysis. J Microbiol Methods 2000,41(2):121-133. 10.1016/S0167-7012(00)00142-1

    Article  PubMed  Google Scholar 

  31. Schroder S, Wittmann D, Drescher W, Roth V, Steinhage V, Cremers AB: The New Key to Bees: Automated Identification by Image Analysis of Wings. In Pollinating Bees - The Conservation Link Between Agriculture and Nature. Edited by: Kevan P, Imperatriz Fonseca VL. Brasilia: Ministry of Environment; 2002:209-216.

    Google Scholar 

  32. MacLeod N, Benfield M, Culverhouse P: Time to automate identification. Nature 2010, 467: 154-155. 10.1038/467154a

    Article  CAS  PubMed  Google Scholar 

  33. Sundbland LG, Geladi P, Dunberg A, Sundberg B: The use of image analysis and automation for mitotic index in apical conifer meristems. J Exp Bot 1998, 49: 1749-1756.

    Article  Google Scholar 

  34. Kellogg RA, Chebira A, Goyal A, Cuadra PA, Zappe SF, Minden JS, Kovačević J: Towards an image analysis toolbox for high-throughput Drosophila embryo RNAi screens. In 4th Institute of Electrical and Electronics Engineers (IEEE) International Symposium on Biomedical Imaging: From Nano to Macro. Arlington, USA; 2007:288-291.

    Google Scholar 

  35. Hopwood N: A history of normal plates, tables and stages in vertebrate embryology. Int J Dev Biol 2007, 51: 1-26. 10.1387/ijdb.062189nh

    Article  PubMed Central  PubMed  Google Scholar 

  36. Rundle SD, Spicer JI, Coleman RA, Vosper J, Soane J: Environmental calcium modified induced defenses in a freshwater gastropod. Proc Biol Sci 2004, 271: 67-70. 10.1098/rsbl.2003.0106

    Article  Google Scholar 

  37. Edelstein A, Amodaj N, Hoover K, Vale R, Stuurman N: Computer control of microscopes using μManager. Curr Protoc Mol Biol 2010, 92: 14.20.1-14.20.17.

    Google Scholar 

  38. Press WH, Flannery BP, Teukolsky SA, Vetterling WT: Fourier Transform of Discretely Sampled Data. In The Art of Scientific Computing. 2nd edition. Cambridge: Cambridge University Press; 1989:494-498.

    Google Scholar 

Download references


We thank Steven Soffe (University of Bristol) for supplying us with Xenopus laevis embryos, Ted Henry and Kath Sloman for supplying us with Danio rerio embryos and Ann Torr for her culture of Radix balthica embryos.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Oliver Tills.

Additional information

Competing interests

The authors have a patent pending, entitled: Method and system for determining characteristics of an embryo and uses thereof (UK Patent application number: GB 1016616.3; Date of filing 2/10/2010).

This work was funded by the Marine Institute at Plymouth University through the HEIF4 programme.

Authors’ contributions

TB and OT carried out experimental work and bioimaging. PC carried out programming required for motion analysis and frequency analysis. TB, OT, SDR and JIS participated in the data analysis. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Tills, O., Bitterli, T., Culverhouse, P. et al. A novel application of motion analysis for detecting stress responses in embryos at different stages of development. BMC Bioinformatics 14, 37 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: