Many soft-bodied invertebrates use a flexible, fluid-filled hydrostatic skeleton for burrowing. The aim of our study was to compare the scaling and morphology between surface-dwelling and burrowing earthworm ecotypes to explore the specializations of non-rigid musculoskeletal systems for burrowing locomotion. We compared the scaling of adult lumbricid earthworms across species and ecotypes to determine whether linear dimensions were significantly associated with ecotype. We also compared the ontogenetic scaling of a burrowing species, Lumbricus terrestris, and a surface-dwelling species, Eisenia fetida, using glycol methacrylate histology. We show that burrowing species are longer, thinner and have higher length-to-diameter ratios than non-burrowers, and that L. terrestris is thinner for any given body mass compared with E. fetida. We also found differences in the size of the musculature between the two species that may correlate with surface crawling or burrowing. Our results suggest that adaptations to burrowing for soft-bodied animals include a disproportionately thin body and strong longitudinal muscles.

Burrowing is a difficult form of locomotion because of the abrasive, heterogeneous and dense nature of many substrates. Despite the challenges, many vertebrates and invertebrates ranging from micrometers to meters in length burrow effectively in a variety of substrates. Their burrowing actions alter the soil environment and aid in nutrient recycling, air and water infiltration, and soil decompaction. Many invertebrate burrowers lack rigid skeletal elements, relying instead on a hydrostatic skeleton consisting of a liquid-filled internal cavity surrounded by a muscular body wall (Chapman, 1958; Kier, 2012). When the muscles in the body wall contract, the internal fluid is pressurized, allowing for skeletal support, muscle antagonism, skeletal leverage, locomotion and other skeletal functions (Chapman, 1950, 1958; Alexander, 1995). The hydrostatic skeleton can also accommodate deformation in the body due to muscle contraction. Earthworms, for example, possess a fully segmented hydrostatic skeleton with two predominant muscle orientations present in each segment – circumferential and longitudinal. Circumferential muscle contraction elongates the worm, allowing it to move forward and excavate a new burrow; the longitudinal muscles expand the worm laterally, allowing for anchorage and burrow consolidation (Trueman, 1975). In addition, the radial straining of the soil by the longitudinal muscles breaks up soil particles ahead of the worm, reducing the pressure required for axial elongation (Abdalla et al., 1969; Whalley and Dexter, 1994; Keudel and Schrader, 1999; Dorgan et al., 2008).

Soft-bodied burrowing invertebrates range in size from several hundred micrometers in length (e.g. nematodes) to several meters in length (e.g. earthworms) and burrow in a variety of terrestrial and marine environments. The effects of size on burrowing mechanics has not, however, been studied in detail (e.g. Piearce, 1983; Quillin, 2000; Che and Dorgan, 2010). In addition, the impact on subterranean organisms of anthropogenic changes in soil properties from chemicals and heavy machinery has been investigated previously, yet we do not know if there are size-dependent effects on burrowers (e.g. Ehlers, 1975; Roberts and Dorough, 1985; Chan and Barchia, 2007). This research may also provide insights important for the design of burrowing soft robots (e.g. Trimmer, 2008; Trivedi et al., 2008; Daltorio et al., 2013).

The physical characteristics of soil may impose size-dependent constraints on burrowers (Dorgan et al., 2008; Che and Dorgan, 2010; Kurth and Kier, 2014). For example, many soils exhibit strain hardening, in which the modulus of compression (stiffness) of the soil increases with increasing strain (Chen, 1975; Yong et al., 2012; Holtz et al., 2010). As an earthworm grows in cross-section, it must displace more soil radially as it burrows, which may result in an increase in the stiffness of the soil surrounding the burrow. Hatchling worms may avoid the strain hardening effect as a result of the relatively small volume of soil they must displace during burrowing. Adult earthworms, however, are often several times longer and wider than hatchlings and must create wider burrows in order to accommodate their bodies (Gerard, 1967; Quillin, 2000; Kurth and Kier, 2014). The formation of new burrows is common when conditions change or resources become scarce, forcing large worms to encounter and overcome strain hardening (Evans, 1947; Gerard, 1967). Thus, as a burrower grows there may be a selective advantage to becoming relatively thinner to mitigate the strain hardening effect (Piearce, 1983; Kurth and Kier, 2014).

In previous work, we showed that the burrowing earthworm Lumbricus terrestris becomes relatively thinner during growth and shows additional allometric changes in the musculature (Kurth and Kier, 2014). We hypothesized that these allometries may help to compensate for changes in strain hardening with growth. In order to examine this issue, we compared the linear dimensions of earthworms across ecotypes, as well as the ontogenetic scaling of a non-burrowing, surface-dwelling earthworm, Eisenia fetida. Not all earthworms burrow; there are three main ecotypes of earthworms that are largely differentiated by their burrowing patterns or lack thereof (Bouché, 1977). Surface-dwelling species like E. fetida are known as epigeic worms, which do not burrow and are instead found under leaf litter, in manure and under debris. There are also endogeic worms, which create ephemeral horizontal burrows in the upper 10–15 cm of soil and are geophagus (Edwards and Bohlen, 1996). Finally, there are anecic worms such as L. terrestris that build deep permanent/semi-permanent vertical burrows and feed on surface litter (Keudel and Schrader, 1999). We refer to these three ecotypes as surface dwellers, horizontal burrowers and vertical burrowers, respectively.

List of symbols and abbreviations

     
  • A

    muscle cross-sectional area

  •  
  • a

    scaling constant

  •  
  • b

    scaling exponent

  •  
  • D

    diameter

  •  
  • F

    force

  •  
  • L

    body length

  •  
  • L/D

    length to diameter ratio

  •  
  • M

    body mass

  •  
  • Mech adv

    mechanical advantage

We hypothesized that there would be both interspecific and ontogenetic scaling differences between earthworm ecotypes. We predicted that to mitigate strain hardening, the burrowing species would be thinner for any given body mass during development compared with surface dwellers, resulting in higher length-to-diameter ratios in the burrowing species. We also hypothesized that forces from the longitudinal musculature, which radially expand the worm during contraction, would be relatively larger and would develop more rapidly in the burrowers compared with the surface dwellers. These muscles are believed to be important in burrowing by anchoring the worm (with assistance from projections of hair-like setae), consolidating the burrow, relieving soil compaction ahead of the worm and pulling posterior segments into the burrow (Seymour, 1969; McKenzie and Dexter, 1988; Keudel and Schrader, 1999; Barnett et al., 2009). These muscles also move the bulk of soil during burrow formation and must generate sufficient force to overcome potential strain hardening effects in the soil (Barnett et al., 2009).

The setae that project outward during longitudinal muscle contraction are necessary in both crawling and burrowing to prevent backslip of the animal, and are similarly arranged in a ‘lumbricine’ pattern in all lumbricid earthworms (Sims and Gerard, 1985). Therefore, we focused our research on differences in longitudinal musculature instead of seta morphology or arrangement, although exploration of this topic may be of interest in future work.

In contrast to the longitudinal muscles, we predicted that forces from the circumferential muscles, which cause thinning and elongation of the worm, would be larger in the surface dwellers. The circumferential muscles are particularly important in surface crawling, extending the worm forward during each peristaltic wave of contraction and aiding penetration into litter and under debris; in fact, the largest pressures exerted in surface-crawling earthworms occur during circumferential muscle contraction (Gray and Lissman, 1938; Arthur, 1965; Seymour, 1969).

We also found significant interspecific and ontogenetic differences in scaling, consistent with our hypotheses (Kurth and Kier, 2014). Our results demonstrate that many aspects of the hydrostatic skeleton of earthworms develop in different ways between species, reflecting the ecological context of the organism.

Scaling of functionally relevant morphological features

A variety of organisms, including L. terrestris, exhibit allometric growth, in which the relative proportions change with body size rather than remaining constant, as in isometric growth (Huxley and Tessier, 1936; Schmidt-Nielsen, 1997). Since the density of an animal typically does not change with size, the mass (M) is proportional to the volume (V). If an organism scales isometrically, linear dimensions such as length (L) or diameter (D) are predicted to scale to the animal's V1/3 and thus M1/3 and any area, such as surface area or muscle cross-sectional area, will scale as V2/3 and thus M2/3.

In earthworms, the size and scaling of morphological features (e.g. diameter, muscle cross-sectional area) can vary from segment to segment down the length of the body (Quillin, 2000; Kurth and Kier, 2014). If an earthworm exhibits isometric growth, these morphological features must scale isometrically across all segments. To account for potential scaling differences across segments, multiple segments were measured in this study.

Scaling of linear dimensions

In an isometrically scaling vermiform animal, the L/D ratio will not change with size. This is because both L and D are linear dimensions and should scale as M1/3 across segments. However, Kurth and Kier (2014) found allometry in the overall dimensions of L. terrestris, which changes the relative force and displacement of the musculature during growth (Kier and Smith, 1985; Vogel, 2013). An increase in the L/D ratio during growth, as is found in L. terrestris, increases the distance advantage (the ratio of distance output to distance input) for the circumferential muscles and increases the mechanical advantage (the ratio of force output to force input) for the longitudinal muscles (Kier and Smith, 1985; Vogel, 2013; Kurth and Kier, 2014). Since mechanical advantage and distance advantage are reciprocal, an increase in the L/D ratio decreases the mechanical advantage of the circumferential musculature and decreases the distance advantage of the longitudinal musculature.

We predict that burrowing species will have higher L/D ratios than surface dwellers because surface dwellers are not under selective pressure to minimize their diameters for burrowing. We also predict that this trend will be reflected ontogenetically, with L. terrestris having a higher L/D ratio and smaller diameter for a given body mass than E. fetida. A lower L/D means that surface dwellers will have lower mechanical advantage during longitudinal muscle contraction and higher mechanical advantage during circumferential muscle contraction for a given size than burrowers.

Scaling of muscle cross-sectional areas and force output

The scaling of muscle physiological cross-sectional area (A) determines how relative force production by the musculature changes with size, because force due to muscle contraction is proportional to cross-sectional area. If the circumferential and the longitudinal musculature scales isometrically, the cross-sectional area in each segment will be proportional to M2/3. The final force output the animal exerts on the environment, however, is a product both of the force generated by the muscles and the mechanical advantage produced by the skeleton itself:
(1)

where F is the force output to the environment, mech adv is the mechanical advantage and A is the muscle cross-sectional area (Kurth and Kier, 2014). We predict that the scaling of force output for E. fetida will be lower during longitudinal muscle contraction but higher during circumferential muscle contraction than L. terrestris.

Interspecific scaling of linear dimensions

Because L/D is dimensionless, we first compared this value across species without respect to body size. We found a significant relationship between L/D and ecotype across species and clades (independent contrasts, P<0.05; Fig. 1). Surface-dwelling worms generally had the lowest L/D of the three ecotypes, whereas vertical burrowers had the highest L/D ratios of the three ecotypes. Horizontal burrowers had moderate L/D, which were significantly higher than surface dwellers and significantly lower than vertical burrowers (P<0.05).

Fig. 1.

Simplified phylogenetic tree comparing L/D ratio and ecotype. The phylogeny is based on a tree built by Pérez-Losada et al. (2012). Text colors indicate ecotype. The numbers adjacent to the phylogeny indicate the L/D value for each species. Each L/D value is an average from three adult specimens per species. No Hormogaster elisae specimens were available for analysis, so it was only used to root the phylogenetic tree.

Fig. 1.

Simplified phylogenetic tree comparing L/D ratio and ecotype. The phylogeny is based on a tree built by Pérez-Losada et al. (2012). Text colors indicate ecotype. The numbers adjacent to the phylogeny indicate the L/D value for each species. Each L/D value is an average from three adult specimens per species. No Hormogaster elisae specimens were available for analysis, so it was only used to root the phylogenetic tree.

These differences in L/D ratios result from dissimilarities in the scaling of length and diameter between ecotypes (Fig. 2). We found that while both burrowing and surface-dwelling species increased in length with similar scaling exponents (b=0.410 and 0.401 for burrowers and surface dwellers, respectively), burrowing species were significantly longer for a given body volume than surface dwellers (a=0.737 and 0.686 for burrowers and surface dwellers, respectively). Burrowers and surface dwellers also increased in diameter at similar rates (b=0.295 and 0.300 for burrowers and surface dwellers, respectively), but burrowers were thinner for a given body volume than surface dwellers (a=−0.316 and −0.291 for burrowers and surface dwellers, respectively).

Fig. 2.

Interspecific differences in the scaling of linear dimensions. (A) Log-transformed graph comparing body length (L) with body volume between burrowing and surface dwelling adult lumbricid species. (B) Log-transformed graph comparing anterior diameter (Dant) with body volume between burrowing and surface-dwelling lumbricid species. The regressions shown in A and B were fitted to empirical data using OLS regression (solid line for burrowers, dashed line for surface dwellers) and the regression equations for both ecotypes are shown. Asterisk indicates a significant difference between species with the Bonferroni correction. N=13 surface-dwelling species; N=15 burrowing species.

Fig. 2.

Interspecific differences in the scaling of linear dimensions. (A) Log-transformed graph comparing body length (L) with body volume between burrowing and surface dwelling adult lumbricid species. (B) Log-transformed graph comparing anterior diameter (Dant) with body volume between burrowing and surface-dwelling lumbricid species. The regressions shown in A and B were fitted to empirical data using OLS regression (solid line for burrowers, dashed line for surface dwellers) and the regression equations for both ecotypes are shown. Asterisk indicates a significant difference between species with the Bonferroni correction. N=13 surface-dwelling species; N=15 burrowing species.

Ontogenetic scaling of linear dimensions

We found a significant difference between the scaling of L/D between the two species (Fig. 3; Table 1). While both L. terrestris and E. fetida grew disproportionately long (bLt=0.393, bEf=0.383) and disproportionately thin (bLt=0.290, 0.275, 0.277; bEf=0.293, 0.300, 0.308; for anterior, middle and posterior segments) at similar rates, E. fetida was always significantly wider at a given body mass than L. terrestris, as shown by the differences in log a, the y-intercept of the log-transformed graph (log aLt=0.630, 0.605, 0.550; log aEf=0.861, 0.883, 0.850 for anterior, middle and posterior segments, respectively) (Table 1; supplementary material Fig. S1). Because of these differences in diameter, E. fetida had a lower L/D for any given body mass compared with L. terrestris (log aLt=1.407, log aEf=1.202; averaged across segments) despite a similar increase in L/D with size for both species (bLt=0.114, bEf=0.087; averaged across segments) (supplementary material Fig. S2).

Fig. 3.

Ontogenetic scaling of linear dimensions. (A) Log-transformed graph comparing body length (L) with body mass (M) for L. terrestris and E. fetida. (B) Log-transformed graph comparing anterior diameter (Dant) with body mass for L. terrestris and E. fetida. The regressions shown in A and B were fitted to empirical data using OLS regression (solid line for L. terrestris, dashed line for E. fetida) and the regression equations for both species are shown. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Fig. 3.

Ontogenetic scaling of linear dimensions. (A) Log-transformed graph comparing body length (L) with body mass (M) for L. terrestris and E. fetida. (B) Log-transformed graph comparing anterior diameter (Dant) with body mass for L. terrestris and E. fetida. The regressions shown in A and B were fitted to empirical data using OLS regression (solid line for L. terrestris, dashed line for E. fetida) and the regression equations for both species are shown. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Table 1.

Scaling of linear dimensions in the vertical burrower L. terrestris and the surface dweller E. fetida

Scaling of linear dimensions in the vertical burrower L. terrestris and the surface dweller E. fetida
Scaling of linear dimensions in the vertical burrower L. terrestris and the surface dweller E. fetida

Ontogenetic scaling of muscle cross-sectional area

We found differences in muscle cross-sectional area between species (Fig. 4; Table 2). In the anterior segments, L. terrestris had larger longitudinal muscles for a given body mass than E. fetida (log aLt=0.511; log aEf=0.354) and its longitudinal muscles grew at faster rates than those of E. fetida (bLt=0.612; bEf=0.539), but these differences were not statistically significant (Fig. 4). Conversely, E. fetida had larger anterior circumferential muscles at a given body mass than L. terrestris (log aLt=−0.713; log aEf=−0.640), despite a faster growth of these muscles in L. terrestris (bLt=0.674; bEf=0.543). Muscles in the middle and posterior segments were more similar in the two species (Table 2; supplementary material Figs S3 and S4). The longitudinal muscles from the middle segments scaled similarly (bLt=0.541; bEf=0.552) and were similar in cross-sectional area at a given body mass (log aLt=0.375; log aEf=0.392), whereas the circumferential muscles grew at a faster rate in L. terrestris (bLt=0.800; bEf=0.627) but were larger at a given body mass in E. fetida (log aLt=−0.974; log aEf=−0.731). The posterior longitudinal segments showed the opposite scaling trend from the anterior segments; the longitudinal muscles of E. fetida increased in cross-sectional area at a faster rate (bLt=0.564; bEf=0.640) and were larger at a given body mass (log aLt=0.379; log aEf=0.437) though these differences were not statistically significant. The posterior circumferential muscles showed no significant difference in scaling exponents between the two species (bLt=0.792; bEf=0.743), although the circumferential muscle cross-sectional area of E. fetida was larger at a given body mass than that of L. terrestris (log aLt=−1.048; log aEf=−0.609).

Fig. 4.

Ontogenetic scaling of muscle cross-sectional areas. (A) Log-transformed graph comparing longitudinal muscle cross-sectional area in the anterior segments (Al,ant) to body mass for L. terrestris and E. fetida. (B) Log-transformed graph comparing circumferential muscle cross-sectional area in the anterior segments (Ac,ant) with body mass for L. terrestris and E. fetida in the anterior segments. The regressions shown in Fig. 1A,B were fitted to empirical data using OLS regression (solid line for L. terrestris, dashed line for E. fetida) and the regression equations for both species are shown. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Fig. 4.

Ontogenetic scaling of muscle cross-sectional areas. (A) Log-transformed graph comparing longitudinal muscle cross-sectional area in the anterior segments (Al,ant) to body mass for L. terrestris and E. fetida. (B) Log-transformed graph comparing circumferential muscle cross-sectional area in the anterior segments (Ac,ant) with body mass for L. terrestris and E. fetida in the anterior segments. The regressions shown in Fig. 1A,B were fitted to empirical data using OLS regression (solid line for L. terrestris, dashed line for E. fetida) and the regression equations for both species are shown. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Table 2.

Scaling of muscle cross-sectional area in the vertical burrower L. terrestris and the surface dweller E. fetida

Scaling of muscle cross-sectional area in the vertical burrower L. terrestris and the surface dweller E. fetida
Scaling of muscle cross-sectional area in the vertical burrower L. terrestris and the surface dweller E. fetida

Ontogenetic scaling of calculated mechanical advantage and force production

Because the L/D ratio increased in both E. fetida and L. terrestris, both had similar trends in the scaling of mechanical advantage (Fig. 5). We calculated increases in mechanical advantage during longitudinal muscle contraction for both species (bLt=0.104; bEf=0.078); L. terrestris had higher mechanical advantage for a given body mass than E. fetida (aLt=1.872; aEf=1.649). The mechanical advantage of the circumferential muscle decreased with growth in both species (bLt=−0.104; bEf=−0.078), but L. terrestris exhibited lower mechanical advantage at a given body mass (log aLt=−1.872; log aEf=−1.649).

Fig. 5.

Comparison of calculated mechanical advantage with body mass. Mechanical advantage was calculated by normalizing each worm's average L/D across segments with mass and calculating the reciprocal of distance advantage over 25% radial strain. (A) Mechanical advantage from longitudinal muscle contraction (mech advl) and (B) circular muscle contraction (mech advc) mechanical advantage as a function of earthworm body mass. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Fig. 5.

Comparison of calculated mechanical advantage with body mass. Mechanical advantage was calculated by normalizing each worm's average L/D across segments with mass and calculating the reciprocal of distance advantage over 25% radial strain. (A) Mechanical advantage from longitudinal muscle contraction (mech advl) and (B) circular muscle contraction (mech advc) mechanical advantage as a function of earthworm body mass. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

We also found significant differences in the scaling of calculated force production between the two species (Fig. 6; Table 3). In the anterior segments (Fig. 6), we found calculated force output during longitudinal muscle contraction at any given body mass was greater for L. terrestris than for E. fetida (log aLt=2.383; aEf=2.003). Calculated longitudinal muscle force production increased at a greater rate with mass in L. terrestris than E. fetida (log bLt=0.716; bEf=0.617), although this difference was not statistically significant with the Bonferroni correction. In the case of calculated circumferential muscle force production, however, E. fetida had a greater force output at a given body mass than did L. terrestris (log aLt=−2.584; log aEf=−2.288), but both species had similar growth rates (bLt=0.568; bEf=0.465).

Fig. 6.

Comparison of calculated force production with body mass. Force production was estimated for each worm using its mechanical advantage and muscle cross-sectional area. Mechanical advantage was calculated by normalizing each worm's L/D ratio with mass and calculating the reciprocal of distance advantage over 25% radial strain (A) Force production during anterior longitudinal muscle contraction (Fl,ant) and (B) force production during anterior circumferential muscle contraction (Fc,ant) as a function of earthworm body mass. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Fig. 6.

Comparison of calculated force production with body mass. Force production was estimated for each worm using its mechanical advantage and muscle cross-sectional area. Mechanical advantage was calculated by normalizing each worm's L/D ratio with mass and calculating the reciprocal of distance advantage over 25% radial strain (A) Force production during anterior longitudinal muscle contraction (Fl,ant) and (B) force production during anterior circumferential muscle contraction (Fc,ant) as a function of earthworm body mass. Asterisk indicates a significant difference between species with the Bonferroni correction. N=25 per species.

Table 3.

Scaling of calculated force production in the vertical burrower L. terrestris and the surface dweller E. fetida

Scaling of calculated force production in the vertical burrower L. terrestris and the surface dweller E. fetida
Scaling of calculated force production in the vertical burrower L. terrestris and the surface dweller E. fetida

We found that most of the differences in force production between E. fetida and L. terrestris were consistent across segments (Table 3; supplementary material Figs S5 and S6). Calculated longitudinal muscle force production in the middle and posterior segments was greater for a given mass in L. terrestris than E. fetida (log aLt=2.245; log aEf=2.041 in middle segments; log aLt=2.251; log aEf=2.086 in posterior segments), but these segments did not show significant inter-specific differences in the rates of longitudinal muscle force production with mass (bLt=0.649, bEf=0.630 in middle segments; log bLt=0.668, log bEf=0.717 in posterior segments). Circumferential muscle force production in the middle and posterior segments also exhibited similar trends to the anterior segments, with higher intercepts in E. fetida (log aLt=−2.838, log aEf=−2.380 in the middle segments; log aLt=−2.920, log aEf=−2.258 in the posterior segments) and similar scaling exponents between the two species (bLt=0.681; bEf=0.550 in the middle segments; bLt=0.688; bEf=0.665 in the posterior segments).

Previous studies found that the hydrostatic skeleton in L. terrestris scales allometrically, but the reasons for these growth patterns remain unclear (Quillin, 2000; Kurth and Kier, 2014). We hypothesized that one important factor may be compensation for increases in soil strain hardening as the animal becomes larger. We tested this hypothesis by comparing the hydrostatic skeleton across ecotypes in earthworms using interspecific and ontogenetic methods. Our results are consistent with the strain hardening hypothesis and suggest that a disproportionately thin diameter and large forces during longitudinal muscle contractions are key burrowing adaptations in soft-bodied animals.

Linear dimensions and mechanical advantage

We found that burrowing species across clades had higher L/D ratios than surface dwellers, consistent with previous research by Piearce (1983). These differences in L/D were reflected in both the interspecific and ontogenetic scaling of linear dimensions. The interspecific scaling comparison shows that both ecotypes grew disproportionately long and thin, but burrowing species were significantly longer and thinner than surface-dwelling species. Ontogenetically, both the burrowing L. terrestris and surface-dwelling E. fetida also grew disproportionately long and thin. At any given body mass, however, L. terrestris was significantly thinner than E. fetida. Since burrowers would experience greater selective pressures for thin bodies than surface dwellers in order to mitigate strain hardening underground, our interspecific and ontogenetic results are consistent with the strain hardening hypothesis.

Ontogenetic changes in mechanical advantage showed similar trends between species since both increased their L/D ratio during growth. For both species, mechanical advantage increased with body size for longitudinal muscle contraction and decreased with body size for circumferential muscle contraction. The magnitudes of mechanical advantage, however, differed between the two species because of differences in L/D ratios. Lumbricus terrestris had greater mechanical advantage during longitudinal muscle contraction, whereas E. fetida had greater mechanical advantage during circumferential muscle contraction. We believe that these differences in mechanical advantage highlight the relative importance of the longitudinal and circumferential muscles in burrowing and crawling, respectively, as discussed below.

We found it surprising that, for E. fetida, the mechanical advantage during circumferential muscle contraction decreased with growth, given the importance of circumferential muscles in surface crawling (Gray and Lissman, 1938; Seymour, 1969). As we discuss below, however, an increase in circumferential cross-sectional area appears to compensate for the loss of mechanical advantage so that the circumferential muscles in E. fetida are significantly larger than those in L. terrestris.

Differences in calculated force production

The segments measured in E. fetida are estimated to produce significantly higher circumferential forces and significantly lower longitudinal muscle forces along the length of the body when compared with similar segments in L. terrestris. These differences agree with previous research suggesting that circumferential muscles are of great importance for crawling while the longitudinal muscles are essential for burrowing (Chapman, 1950; Seymour, 1969). Powerful circumferential muscle forces would permit surface-dwelling worms to squeeze in between rocks, litter and debris on the surface, allowing these worms access to new habitats. Conversely, robust longitudinal muscle forces would allow burrowing earthworms to overcome strain hardening in soil by exerting sufficient force to laterally displace soil, expand the burrow walls, break up soil particles ahead of the burrow, anchor the worm (with assistance from the setae) and pull posterior segments into the burrow (Seymour, 1969; McKenzie and Dexter, 1988; Keudel and Schrader, 1999; Barnett et al., 2009).

Scaling similarities

Although our results showed significant differences in the magnitude of musculoskeletal dimensions and calculated forces (i.e. different intercepts) between surface dwellers and burrowers, it is unclear why both burrowers and surface dwellers exhibit scaling similarities (i.e. similar scaling exponents). For example, both burrowing and surface-dwelling ecotypes grow disproportionately long and thin and are predicted to exhibit similar increases in circumferential and longitudinal muscle forces with size. These shared scaling trends may be the result of ecological, physiological or functional similarities between the species.

For instance, both ecotypes may growth disproportionately thin because the relative surface area for gas exchange would be enhanced in larger individuals. Since the burrowing earthworms are more likely to encounter hypoxic regions than surface dwellers, there may be increased selection pressure for a high L/D ratio in burrowing species.

Similar increases in the rates of force production with size may result from the shared functions of these muscles across ecotypes. The circumferential muscles in all earthworms must grow sufficiently powerful to push the animal forward and excavate through debris or soil. The longitudinal muscles in all species must provide sufficient forces to anchor the earthworm, prevent backslip, pull posterior segments forward and push away constrictive soil or debris.

Interspecific measurements and phylogenetic reconstruction

We used sexually mature earthworm specimens preserved in 70–95% ethanol in the collections of the Smithsonian Institution Museum of Natural History (Washington, DC). A phylogeny is available of species in the Lumbricidae family (Pérez-Losada et al., 2012), so we focused our analysis on genera from this family to avoid pseudoreplication (Felsenstein, 1985). We further narrowed the study by only comparing lumbricid species whose ecotypes are well documented, for a total of 29 species studied (Bouché, 1977; Sims and Gerard, 1985; Edwards and Bohlen, 1996). For each species, we chose museum specimens with the most recent fixation dates to minimize tissue distortion from fixation.

Most worms appeared elongate in fixation, and we selected only these elongate worms for measurement to ensure consistent resting lengths across specimens. To measure each species' linear dimensions, we first pulled each worm by the anterior end along the bench surface in order to straighten the body and ensure similar elongation. We used calipers to measure the length and anterior diameter of three adult specimens per species (for a total of 87 worms measured) and calculated an average length and diameter, which was then used to calculate the length-to-diameter ratio. We were also able to compare the interspecific scaling of linear dimensions in burrowing and surface-dwelling ecotypes. Because many specimens we measured had been dissected and were missing inner organs, we used body volume as a proxy for body mass. No Hormogaster elisae specimens were available for analysis, so it was only used to root the phylogenetic tree as discussed below.

We used TreeGraph2™ (Stöver and Müller, 2010) to construct a simplified phylogeny based on a lumbricid earthworm phylogeny by Pérez-Losada et al. (2012). Pérez-Losada et al. (2012) used molecular data from multiple specimens of each species, which resulted in significant variation in branch length and branch placement between specimens within a species. The authors attributed this variation to the sampling of cryptic species. Because we do not know which specimens were misidentified, we simplified the phylogeny by placing each species in the clade where most specimens per species appeared. Because of the high variation and uncertainty in branch length, we also made all branch lengths equal in our simplified tree. Although this reduces our statistical power, the reduction is relatively minor and tends to produce only false negative results (Grafen, 1989; Martins and Garland, 1991; Swenson, 2009).

Eisenia fetida collection and maintenance

Eisenia fetida earthworms were supplied by Uncle Jim's Worm Farm (Spring Grove, PA, USA) as well as raised from hatchlings bred in a colony maintained in the laboratory. Adult worms (∼0.1–0.7 g) were purchased, raised from purchased juveniles and colony hatchlings. Hatchlings were raised from cocoons deposited by adults bred in the laboratory colony. All worms were housed in plastic bins filled with moist peat moss (Inouye et al., 2006) at 15°C (Presley et al., 1996) and fed dried infant oatmeal (Ownby et al., 2005).

Histology and morphometrics

The measurements and calculations follow those described in Kurth and Kier (2014) for L. terrestris in order to allow consistent comparisons between E. fetida and L. terrestris. Each E. fetida worm was anesthetized in a 10% ethanol solution in distilled water (v/v) until quiescent, patted dry and weighed. The length was obtained after dragging the worm by the anterior end along the lab bench to straighten the body and extend the segments to a consistent resting length. Because E. fetida does not add segments with growth, we measured the length of the entire body (supplementary material Fig. S7). The worm was then killed and three blocks of tissue containing 20 segments each were removed (segments 1–20, 21–40 and 41–60, numbering from anterior) to account for morphological differences across segments. We focused on segments in the anterior half of the worm since it is of greatest importance in locomotion (Yapp, 1956).

The tissue blocks were fixed in 10% formalin in distilled water (v/v) for 24–48 h. After fixation, the blocks were further dissected for embedding (segments 9–14, 29–34 and 49–54). We refer to segments 9–14 as ‘anterior’, segments 29–34 as ‘middle’ and segments 49–54 as ‘posterior’. The anterior, middle and posterior segments were then cut so that both transverse and sagittal sections could be obtained from each location.

The tissue blocks were partially dehydrated in 95% ethanol and embedded in glycol methacrylate plastic (Technovit 7100, Heraeus Kulzer GmbH, Wehrheim, Germany). Sections of 3–7 μm thickness were cut with a glass knife. We used a Picrosirius/Fast Green stain in order to differentiate muscle from connective tissue (López-De León and Rojkind, 1985). We stained the slides at 60°C for 1–2 h then rinsed the slides in distilled water, dried them and mounted coverslips. We used Sigma Scan (Systat Software, San Jose, CA, USA) to make morphological measurements on micrographs. Longitudinal muscle cross-sectional area (Al) and diameter (D) were measured from transverse sections, whereas circumferential muscle cross-sectional area (Ac) was measured from sagittal sections (Fig. 7). The earthworms prepared in this way were flattened slightly and thus had an elliptical cross-section. To determine an equivalent diameter of a circular cylinder we measured the major and minor axes, calculated the area of the ellipse and then calculated the diameter of a circle of the same area.

Fig. 7.

Photomicrographs of 0.1 g Eisenia fetida and 0.1 g Lumbricus terrestris specimens. All sections were 7 μm in thickness and were stained with Picrosirius/Fast Green. (A) Transverse section of L. terrestris showing the cross-sectional area of the longitudinal musculature. (B) Transverse section of E. fetida showing the cross-sectional area of the longitudinal musculature. (C) Parasagittal section of L. terrestris showing the cross-sectional area of the circumferential musculature. (D) Parasagittal section of E. fetida showing the cross-sectional area of the circumferential musculature. LM, longitudinal muscle; CM, circumferential muscle.

Fig. 7.

Photomicrographs of 0.1 g Eisenia fetida and 0.1 g Lumbricus terrestris specimens. All sections were 7 μm in thickness and were stained with Picrosirius/Fast Green. (A) Transverse section of L. terrestris showing the cross-sectional area of the longitudinal musculature. (B) Transverse section of E. fetida showing the cross-sectional area of the longitudinal musculature. (C) Parasagittal section of L. terrestris showing the cross-sectional area of the circumferential musculature. (D) Parasagittal section of E. fetida showing the cross-sectional area of the circumferential musculature. LM, longitudinal muscle; CM, circumferential muscle.

Calculation of mechanical advantage and force output

The scaling of mechanical advantage was calculated using the average L/D ratio across segments for each worm. The average L/D ratio was observed to change as a function of size and thus the mechanical advantage of the musculature changes during growth (supplementary material Fig. S2; Kier and Smith, 1985; Kurth and Kier, 2014). Since the mechanical advantage is the reciprocal of the distance advantage, we calculated the mechanical advantage of the circumferential musculature as the absolute value of the decrease in body diameter (D) during circumferential muscle contraction divided by the resulting increase in body length (L), as a function of the L/D ratio, for a 25% decrease body in diameter. Kinematic data for E. fetida are unavailable but a 25% change in body diameter has been empirically recorded from kinematic data in L. terrestris during crawling (Quillin, 1999). Since L. terrestris and E. fetida are closely related phylogenetically, we believe that a 25% change in diameter is a reasonable assumption for both species. Likewise, the mechanical advantage of the longitudinal muscle was calculated as the absolute value of the decrease in body length of the worm divided by the resulting increase in body diameter, as a function of the L/D ratio:
(2)

These calculations thus provided estimates of the mechanical advantage of both the longitudinal (l) and circumferential (c) musculature as a function of size.

Force production was calculated in each worm as the product of mechanical advantage and muscle cross-sectional area in both the circumferential and longitudinal muscles. We made the assumption of constant stress in the muscles with ontogeny, though this assumption has not been empirically tested in obliquely striated muscle.

Statistical analysis

We used R statistical software for both the phylogenetic and ontogenetic analyses (R Development Core Team, 2014). We used the ape package in R (Paradis et al., 2004) to perform independent contrasts on the phylogeny for statistical analysis. Independent contrasts allowed us to test for correlations between ecotype and L/D ratio while avoiding pseudoreplication (Felsenstein, 1985). We treated ecotype as a continuous variable to allow for transitional/intermediate ecotypes in ancestral nodes.

We also used linear regression on log transformed interspecific and ontogenetic scaling data. We fitted both sets of scaling data to the power function y=aMb, where y represents the morphological trait of interest, a is the scaling constant, M is body mass and b is the scaling exponent. Log transforming these data allowed us to perform regression analyses, as b becomes the slope of the line and log a becomes the intercept.

We used the caper function (Orme et al., 2012) in R (R Development Core Team, 2014) to perform phylogenetically corrected regression on the interspecific scaling data. We pooled horizontal and vertical burrowers together for this analysis because only three vertical burrowing species were measured, and all three were similar in body size. To test for differences in slope and intercept between burrowing and surface-dwelling ecotypes, we performed an ANCOVA analysis on the phylogenetically corrected regression data. Although there may be error in the x-variable (i.e. volume) that is not accounted for in a standard ANCOVA, ANCOVAs using model II regression techniques are not well developed or commonly used (Sokal and Rohlf, 1994). Thus, standard ANCOVAs are still commonly used in scaling studies (e.g. Niven and Scharlemann, 2005; Davies and Moyes, 2007; Snelling et al., 2011; Xu et al., 2014).

We used the lmodel2 package (Legendre, 2011) in R to perform ordinary least squares (OLS) and reduced major axis regression (RMA) on the ontogenetic scaling data. In our ontogenetic analysis, the symbols bLt and log aLt denote the slopes and intercepts of L. terrestris, while the symbols bEf and log aEf denote the slopes and intercepts of E. fetida. To determine differences in slope and intercept between the two species, we used a standard ANCOVA. We also compared RMA scaling exponent bLt and constant aLt for L. terrestris against the corresponding 95% confidence intervals for E. fetida (Heins et al., 2004). Our data generally showed high coefficients of determination (R2>0.85) and both OLS regression and RMA regression fit similar scaling exponents in our analysis and were consistent in distinguishing significant scaling differences between species. Because of the similarity and agreement between the models, only the ANCOVA and OLS results for both species are reported to remain consistent with the statistical reporting from the interspecific scaling study. To account for multiple comparisons, a Bonferroni correction was used on the P-value outputs from the ANCOVAs. Most P-values remained significant.

We thank Brina M. Montoya for her insights into soil properties, Tony Purdue for help with histology, and Karen Osborn for access to the annelid collections at the Smithsonian Institution.

Author contributions

Both authors were involved in the planning of the experiments. J.A.K. collected and analyzed the data, and drafted the initial manuscript. W.M.K. revised the manuscript and assisted with data analysis.

Funding

This study was funded by the National Science Foundation [Grant IOS-0951067 to W.M.K.].

Abdalla
,
A. M.
,
Hettiaratchi
,
D. R. P.
and
Reece
,
A. R.
(
1969
).
The mechanics of root growth in granular media
.
J. Agric. Eng. Res.
14
,
236
-
248
.
Alexander
,
R. McN
. (
1995
).
Hydraulic mechanisms in locomotion
. In
Body Cavities: Function and Phylogeny (Selected Symposia and Monographs UZI, 8)
(ed.
G.
Lanzavecchia
,
R.
Valvassori
and
M. D.
Candia Carnevali
), pp.
187
-
198
.
Modena
:
Mucchi
.
Arthur
,
D. R.
(
1965
).
Form and function in the interpretation of feeding in lumbricid worms
.
View. Biol.
4
,
204
-
251
.
Barnett
,
C. M.
,
Bengough
,
A. G.
and
McKenzie
,
B. M.
(
2009
).
Quantitative image analysis of earthworm-mediated soil displacement
.
Biol. Fertil. Soils
45
,
821
-
828
.
Bouché
,
M. B.
(
1977
).
Strategies lombriciennes
.
Ecol. Bull.
25
,
122
-
132
.
Chan
,
K. Y.
and
Barchia
,
I.
(
2007
).
Soil compaction controls the abundance, biomass and distribution of earthworms in a single dairy farm in south-eastern Australia
.
Soil Till. Res.
94
,
75
-
82
.
Chapman
,
G.
(
1950
).
Of the movement of worms
.
J. Exp. Biol.
27
,
29
-
39
.
Chapman
,
G.
(
1958
).
The hydrostatic skeleton in the invertebrates
.
Biol. Rev.
33
,
338
-
371
.
Che
,
J.
and
Dorgan
,
K. M.
(
2010
).
It's tough to be small: dependence of burrowing kinematics on body size
.
J. Exp. Biol.
213
,
1241
-
1250
.
Chen
,
W. F.
(
1975
).
Limit Analysis and Soil Plasticity, Developments in Geotechnical Engineering
, Vol.
7
.
Amsterdam
:
Elsevier
.
Daltorio
,
K. A.
,
Boxerbaum
,
A. S.
,
Horchler
,
A. D.
,
Shaw
,
K. M.
,
Chiel
,
H. J.
and
Quinn
,
R. D.
(
2013
).
Efficient worm-like locomotion: slip and control of soft-bodied peristaltic robots
.
Bioinspir. Biomim.
8
,
035003
.
Davies
,
R.
and
Moyes
,
C. D.
(
2007
).
Allometric scaling in centrarchid fish: origins of intra- and inter-specific variation in oxidative and glycolytic enzyme levels in muscle
.
J. Exp. Biol.
210
,
3798
-
3804
.
Dorgan
,
K. M.
,
Arwade
,
S. R.
and
Jumars
,
P. A.
(
2008
).
Worms as wedges: effects of sediment mechanics on burrowing behavior
.
J. Mar. Res.
66
,
219
-
254
.
Edwards
,
C. A.
and
Bohlen
,
P. J.
(
1996
).
Biology and Ecology of Earthworms
, Vol.
3
, 3rd edn.
London
:
Chapman and Hall
.
Ehlers
,
W.
(
1975
).
Observations on earthworm channels and infiltration on tilled and untilled loess soil
.
Soil Sci.
119
,
242
-
249
.
Evans
,
A. C.
(
1947
).
LVII.—A method of studying the burrowing activities of earthworms
.
Ann. Mag. Nat. Hist: Ser. 11.
14
,
643
-
650
.
Felsenstein
,
J.
(
1985
).
Phylogenies and the comparative method
.
Am. Nat.
125
,
1
-
15
.
Gerard
,
B. M.
(
1967
).
Factors affecting earthworms in pastures
.
J. Anim. Ecol.
36
,
235
-
252
.
Grafen
,
A.
(
1989
).
The phylogenetic regression
.
Philos. Trans. R. Soc. B Biol. Sci.
326
,
119
-
157
.
Gray
,
J.
and
Lissman
,
H. W.
(
1938
).
Studies in animal locomotion. VII. Locomotory reflexes in the earthworm
.
J. Exp. Biol.
15
,
506
-
517
.
Heins
,
D. C.
,
Baker
,
J. A.
and
Guill
,
J. M.
(
2004
).
Seasonal and interannual components of intrapopulation variation in clutch size and egg size of a darter
.
Ecol. Freshwat. Fish
13
,
258
-
265
.
Holtz
,
R. D.
,
Kovacs
,
W. D.
and
Sheahan
,
T. C.
(
2010
).
Introduction to Geotechnical Engineering
, pp.
864
.
Upper Saddle River, NJ
:
Prentice Hall
.
Huxley
,
J. S.
and
Tessier
,
G.
(
1936
).
Terminology of relative growth
.
Nature
137
,
780
-
781
.
Inouye
,
L. S.
,
Jones
,
R. P.
and
Bednar
,
A. J.
(
2006
).
Tungsten effects on survival, growth, and reproduction in the earthworm, Eisenia fetida
.
Environ. Toxicol. Chem.
25
,
763
-
768
.
Keudel
,
M.
and
Schrader
,
S.
(
1999
).
Axial and radial pressure exerted by earthworms of different ecological groups
.
Biol. Fertil. Soils
29
,
262
-
269
.
Kier
,
W. M.
(
2012
).
The diversity of hydrostatic skeletons
.
J. Exp. Biol.
215
,
1247
-
1257
.
Kier
,
W. M.
and
Smith
,
K. K.
(
1985
).
Tongues, tentacles and trunks: the biomechanics of movement in muscular-hydrostats
.
Zool. J. Linn. Soc.
83
,
307
-
324
.
Kurth
,
J. A.
and
Kier
,
W. M.
(
2014
).
Scaling of the hydrostatic skeleton in the earthworm Lumbricus terrestris
.
J. Exp. Biol.
217
,
1860
-
1867
.
Legendre
,
P.
(
2011
).
lmodel2: Model II Regression. R package version 1.7-0. http://CRAN.R-project.org/package=lmodel2
.
López-De León
,
A.
and
Rojkind
,
M.
(
1985
).
A simple micromethod for collagen and total protein determination in formalin-fixed paraffin-embedded sections
.
J. Histochem. Cytochem.
33
,
737
-
743
.
Martins
,
E. P.
and
Garland
,
T.
(
1991
).
Phylogenetic analyses of the correlated evolution of continuous characters: a simulation study
.
Evolution
45
,
534
-
557
.
McKenzie
,
B. M.
and
Dexter
,
A. R.
(
1988
).
Radial pressures generated by the earthworm Aporrectodea rosea
.
Biol. Fertil. Soils
5
,
328
-
332
.
Niven
,
J. E.
and
Scharlemann
,
J. P. W.
(
2005
).
Do insect metabolic rates at rest and during flight scale with body mass?
Biol. Lett.
1
,
346
-
349
.
Orme
,
D.
,
Freckleton
,
R.
,
Thomas
,
G.
,
Petzoldt
,
T.
,
Fritz
,
S.
,
Isaac
,
N.
and
Pearse
,
W.
(
2012
).
CAPER: comparative analyses of phylogenetics and evolution in R. R package version 0.5. http://CRAN.R-project.org/package=caper
.
Ownby
,
D. R.
,
Galvan
,
K. A.
and
Lydy
,
M. J.
(
2005
).
Lead and zinc bioavailability to Eisenia fetida after phosphorus amendment to repository soils
.
Environ. Pollut.
136
,
315
-
321
.
Paradis
,
E.
,
Claude
,
J.
and
Strimmer
,
K.
(
2004
).
APE: analyses of phylogenetics and evolution in R language
.
Bioinformatics
20
,
289
-
290
.
Pérez-Losada
,
M.
,
Bloch
,
R.
,
Breinholt
,
J. W.
,
Pfenninger
,
M.
and
Domínguez
,
J.
(
2012
).
Taxonomic assessment of Lumbricidae (Oligochaeta) earthworm genera using DNA barcodes
.
Eur. J. Soil Biol.
48
,
41
-
47
.
Piearce
,
T. G.
(
1983
).
Functional morphology of lumbricid earthworms, with special reference to locomotion
.
J. Nat. Hist.
17
,
95
-
111
.
Presley
,
M. L.
,
McElroy
,
T. C.
and
Diehl
,
W. J.
(
1996
).
Soil moisture and temperature interact to affect growth, survivorship, fecundity, and fitness in the earthworm Eisenia fetida
.
Comp. Biochem. Physiol. A Physiol.
114
,
319
-
326
.
Quillin
,
K. J.
(
1999
).
Kinematic scaling of locomotion by hydrostatic animals: ontogeny of peristaltic crawling by the earthworm Lumbricus terrestris
.
J. Exp. Biol.
202
,
661
-
674
.
Quillin
,
K. J.
(
2000
).
Ontogenetic scaling of burrowing forces in the earthworm Lumbricus terrestris
.
J. Exp. Biol.
203
,
2757
-
2770
.
R Development Core Team
(
2014
).
R: A Language and Environment for Statistical Computing
.
Vienna
,
Austria
:
R Foundation for Statistical Computing
. .
Roberts
,
B. L.
and
Dorough
,
H. W.
(
1985
).
Hazards of chemicals to earthworms
.
Environ. Toxicol. Chem.
4
,
307
-
323
.
Schmidt-Nielsen
,
K.
(
1997
).
Animal Physiology: Adaptation and Environment
, 5th edn.
Cambridge
:
Cambridge University Press
.
Seymour
,
M. K.
(
1969
).
Locomotion and coelomic pressure in Lumbricus terrestris L
.
J. Exp. Biol.
51
,
47
-
58
.
Sims
,
R. W.
and
Gerard
,
B. M.
(
1985
).
Earthworms: Keys and Notes for the Identification and Study of the Species
, Vol.
31
.
London
:
Brill Archive
.
Snelling
,
E. P.
,
Seymour
,
R. S.
,
Matthews
,
P. G. D.
,
Runciman
,
S.
and
White
,
C. R.
(
2011
).
Scaling of resting and maximum hopping metabolic rate throughout the life cycle of the locust Locusta migratoria
.
J. Exp. Biol.
214
,
3218
-
3224
.
Sokal
,
R. R.
and
Rohlf
,
F. J.
(
1994
).
Biometry: The Principles and Practice of Statistics in Biological Research
, 3rd edn.
New York
:
W.H. Freeman and Company
.
Stöver
,
B. C.
and
Müller
,
K. F.
(
2010
).
TreeGraph 2: combining and visualizing evidence from different phylogenetic analyses
.
BMC Bioinformatics
11
,
7
.
Swenson
,
N. G.
(
2009
).
Phylogenetic resolution and quantifying the phylogenetic diversity and dispersion of communities
.
PLoS ONE
4
,
e4390
.
Trimmer
,
B. A.
(
2008
).
New challenges in biorobotics: incorporating soft tissue into control systems
.
Appl. Bionics Biomech.
5
,
119
-
126
.
Trivedi
,
D.
,
Rahn
,
C. D.
,
Kier
,
W. M.
and
Walker
,
I. D.
(
2008
).
Soft robotics: biological inspiration, state of the art, and future research
.
Appl. Bionics Biomech.
5
,
99
-
117
.
Trueman
,
E. R.
(
1975
).
The Locomotion of Soft-Bodied Animals
.
New York
:
American Elsevier Publishing Company, Inc
.
Vogel
,
S.
(
2013
).
Comparative Biomechanics: Life's Physical World
, 2nd edn.
Princeton, NJ
:
Princeton University Press
.
Whalley
,
W. R.
and
Dexter
,
A. R.
(
1994
).
Root development and earthworm movement in relation to soil strength and structure
.
Arch. Agron. Soil Sci.
38
,
1
-
40
.
Xu
,
L.
,
Snelling
,
E. P.
and
Seymour
,
R. S.
(
2014
).
Burrowing energetics of the Giant Burrowing Cockroach Macropanesthia rhinoceros: an allometric study
.
J. Insect Physiol.
70
,
81
-
87
.
Yapp
,
W. B.
(
1956
).
Locomotion of worms
.
Nature
177
,
614
-
615
.
Yong
,
R. N.
,
Nakano
,
M.
and
Pusch
,
R.
(
2012
).
Environmental Soil Properties and Behavior
.
Boca Raton, FL
:
CRC Press Llc
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information