ABSTRACT
We analyzed midline kinematics and obtained electromyograms (EMGs) from the superficial red muscle at seven longitudinal positions in four largemouth bass swimming steadily at standardized speeds of 0.7, 1.2, 1.6, 2.0 and 2.4 lengths s-1. Analysis of variance was used to test for significant variation attributable to both speed and longitudinal position. EMGs propagated posteriorly were unilateral and alternated between the left and right sides. Despite the propagation of EMGs, all the red muscle along an entire side of the fish was simultaneously active for as much as one-quarter of the locomotor cycle. When expressed as a proportion of the locomotor cycle, EMG durations at a given site did not vary significantly with speed but did vary longitudinally, ranging from values of 0.45 cycles anteriorly to 0.35 cycles posteriorly. The amplitudes of lateral displacement and bending depended on longitudinal position and also increased by a maximum of approximately 50 % with increased swimming speed, whereas for all longitudinal positions the intensity of EMGs increased approximately fourfold with increased swimming speed. Electrical activity of red muscle did not correspond simply to the time of muscle shortening. Instead, the timing of EMG onset and offset relative to both lateral bending and displacement changed significantly with both longitudinal position and increased speed, such that the phase shifts between the EMGs and kinematic values were generally greatest for posterior sites at the fastest speeds. At a single longitudinal position, the phase shift between the EMG and maximal lateral bending could change by more than one-tenth of a cycle from the slowest to the fastest swimming speed. Phase lags per body segment of EMG onset and EMG offset did not vary significantly with either swimming speed or longitudinal position. EMG offset was propagated posteriorly faster than EMG onset, and both EMG onset and EMG offset were generally propagated faster than both lateral bending and displacement. Largemouth bass have a similar number of vertebrae to carp, and these two species also have a very similar pattern of muscle activation that differs substantially from that of the trout, which has nearly twice as many vertebrae.
Introduction
The rhythmic axial undulations of fishes swimming at steady speeds have been an exemplary experimental system for gaining insights into the neural control and muscle function that are responsible for locomotor movements. For example, fictive preparations of such diverse species as sharks (Bone, 1966), lamprey (Wallén and Williams, 1984) and goldfish (Fetcho and Svoboda, 1993) have shown that circuitry within the spinal cord is sufficient to generate a pattern of alternating unilateral activation of muscle that is posteriorly propagated and is characteristic of in vivo observations of muscle activity. In addition, recent application of the work loop method for studying in vitro preparations of fish muscle have yielded considerable insights into power output and how it is affected by differences in swimming speed and longitudinal location within the body of a fish (Altringham and Johnston, 1990; Rome and Swank, 1992; Rome et al. 1993; Johnson et al. 1994). Furthermore, the spatial segregation of different fiber types in the axial muscles of fish has facilitated these in vitro studies of muscle physiology as well as in vivo observations of muscle recruitment that are associated with speed (Bone et al. 1978; Johnston and Moon, 1980; Rome et al. 1985; Jayne and Lauder, 1994a).
Despite a large literature on the kinematics and muscle activity involved in the undulatory swimming of vertebrates, much remains to be learned about the interrelationships between morphological variation and the coordination of muscle activity. In contrast to the relatively simple posterior propagation of motor output exemplified by fictive preparations, the propagation of swimming movements of intact organisms appears likely to be a complex result of motor output, the passive and active mechanical properties of the organism and the resistive forces imposed by the environment (Blight, 1977; Jayne and Lauder, 1993; Long et al. 1994). The best-studied fish taxa such as lamprey, eel, trout and carp (Grillner and Kashin, 1976; Williams et al. 1989; van Leeuwen et al. 1990) include so much phylogenetic and morphological diversity that it is difficult to interpret the consequences of any one particular morphological difference. For example, compared with many subcarangiform swimmers, eels and lampreys have both a more elongate external shape and greater numbers of vertebrae (Scott and Crossman, 1973); therefore, external shape differences confound understanding the effects of vertebral number and vice versa.
To isolate the effects of external body shape on the relationship between muscle activity and undulatory swimming movements, we have been quantifying morphology and conducting a series of experiments on two closely related subcarangiform swimmers, the largemouth bass Micropterus salmoides and the bluegill sunfish Lepomis macrochirus (Jayne and Lauder, 1993, 1994a,b, 1995a,b). Both of these species belong to a well-defined monophyletic group of perciform teleost fishes (Lauder and Liem, 1983; Mabee, 1993) that has minimal variation in the numbers of vertebrae (Scott and Crossman, 1973), despite the fact that the bass conforms to a generalized fusiform shape (Webb, 1994) whereas the bluegill has a more specialized deep-bodied shape.
Because much of the previous electromyographic research with the subcarangiform swimming of fishes has been motivated primarily by interest in fiber type recruitment (Hudson, 1973; Bone et al. 1978; Johnston and Moon, 1980; Rome et al. 1985), electromyograms (EMGs) have generally been obtained for few (three or fewer) longitudinal locations, and simultaneous kinematic information has often been sparse (for notable exceptions, see Williams et al. 1989; van Leeuwen et al. 1990). Consequently, the longitudinal extent of muscle activation, the relationship between EMGs and movement, and the effects of longitudinal position on this relationship are still not well understood for most subcarangiform swimmers. Kinematic studies have commonly documented the effects of speed on the amplitude of lateral displacement, but data on the amplitude of lateral bending are less abundant. Because the strain of the axial musculature is proportional to lateral bending, this is a critical kinematic variable to quantify in order to understand in vivo muscle function (Rome et al. 1988). Thus, even for a single species of fish, it is presently uncommon to have integrated information on the muscle activity associated with both lateral displacement and bending at several longitudinal locations and on the effects of speed on these quantities. We studied the steady swimming of largemouth bass specifically to obtain such an integrated understanding of muscle activity and movement.
The goal of the present study is to analyze red muscle activity and axial kinematics during steady swimming of the largemouth bass, specifically with the aims of (1) determining muscle activity from a large number of longitudinal locations to facilitate inferring the whole-body pattern of muscle activation, (2) quantifying the amplitude and timing of both lateral displacement and bending, and relating these to the timing of muscle activity, and (3) determining the effect of locomotor speed on the interrelationships between kinematics and muscle activity.
Materials and methods
Experimental subjects and protocol
We obtained specimens of largemouth bass, Micropterus salmoides (Lacepède), from ponds in southern and central California. The bass were fed a maintenance diet of earthworms and goldfish and were held in captivity before experiments for less than 3 months. All fish were maintained individually in 40 or 80 l tanks at a temperature of 20±2 °C, which averaged the same value as the water temperature that was used during experiments (20±0.5 °C).
After standardizing methods during preliminary experiments with three individuals, we performed experiments with seven additional individuals from which the data were analyzed. Four of these seven individuals were chosen for complete statistical analysis with the goal of obtaining a maximally balanced statistical design for the analyses of variance (described below). Individual fish not used for the complete statistical analysis had missing data due either to electrodes being pulled out during recording sessions or to subsequent determination that the electrode tips were not in the superficial red muscle fibers.
For the four individuals analyzed, mean values (and ranges) of mass, standard length and total length (L) were 88 g (77–99 g), 16.7 cm (15.8–17.5 cm) and 19.8 cm (18.9–20.4 cm), respectively. Mean tail span approximated 5.1 cm. The numbers of precaudal vertebrae ranged from 14 to 16, and the numbers of vertebrae between the trunk and hypurals ranged from 15 to 16.
After implanting electrodes (see below), we placed the fish into the flow tank and allowed them to recover from anesthesia for 2–3 h before we elicited swimming at steady speeds of 0.7, 1.2, 1.6, 2.0 and 2.4 L s-1. The working section of the calibrated flow tank was 18 cm×18 cm×46 cm, with the long dimension being parallel to the flow direction. As in previous experiments with steady swimming (Jayne and Lauder, 1995a), we only analyzed results from trials where the variation in swimming speed among successive tail beats was less than 5 %. Fish commonly swam nearly in the center of the working section of the flow tank, and we only used trials for which the fish were at least 4 cm away from any wall of the tank and the surface of the water. The gap span ratio of the swimming fish (Webb, 1993) approximated 2.4 when fish were centered and 1.5 for sequences in which fish were closest to a wall of the flow tank. At the conclusion of each experiment, fish were killed with an overdose of anesthetic and fixed in formalin to allow post-mortem dissections and X-ray photograph.
Kinematics
Video tapes (synchronized with electromyographic recordings) were obtained simultaneously for both a lateral view and a back-lit ventral view (Jayne and Lauder, 1995a) of the swimming fish using a NAC HSV-400 high-speed video system operating at a rate of 200 images s-1. The lateral view was used primarily to ensure that the vertical position of the fish was approximately midway between the bottom of the tank and the surface of the water.
For each swimming speed, video images were digitized at equally spaced time intervals approximating 1/20 of the tail-beat cycle duration for a total time interval slightly greater than four cycles at each speed for each individual. Ventral-view images were analyzed by digitizing 25–30 points along each side of the fish using custom-made software. For each digitized image, additional software was used to reconstruct a smoothed outline and the midline of the fish. As described in greater detail in Jayne and Lauder (1995a), the midline was partitioned into lengths equaling the skull, individual vertebrae, the hypural bones and intervals one-fifth of the length of the caudal fin (Fig. 1). We obtained the anatomical dimensions necessary for partitioning the midline from X-ray photographs of each of the four fish used for statistical analysis (see Jayne and Lauder, 1995a).
Midline kinematic data were filtered with a low-pass (finite impulse response) filter with a cut-off frequency approximating 2–3 times the tail-beat frequency. From the reconstructed fish midlines, we determined the maximum amplitudes of the filtered values of lateral displacement (zmax) and the angle of lateral bending (βmax) between adjacent axial segments for seven longitudinal locations corresponding to the sites of the implanted electrodes as well as for reference points at the end of the vertebral column and near the trailing edge of the caudal fin (Fig. 1). Amplitudes of kinematic variables were measured as half of the difference between two successive maxima and minima during half of the tail-beat period; hence, during each cycle, values fluctuated between ±zmax and ±βmax. For both zmax and βmax within each tail-beat cycle, lag times were calculated relative to the posterior end of the vertebral column and the other longitudinal locations, and these times were converted to units of phase (z phase, β phase) by dividing by the duration of the tail-beat cycle. For each longitudinal position, we also calculated the phase shift (z-β shift) between zmax and βmax such that negative values indicate that zmax preceded βmax (as described in Jayne and Lauder, 1995a).
Electromyography
After anesthetizing the fish, we implanted a total of thirteen bipolar electrodes in superficial red muscle at the locations shown in Fig. 1 and Table 1. All thirteen channels were simultaneously used for recording during the experiments. Procedures for electrode construction and anesthesia are described in more detail in Jayne (1988) and Jayne and Lauder (1995b). We only analyzed results from electrodes that remained in position for the entire experimental period and after preservation. We used post-mortem dissections and X-ray photographs to confirm the exact location of each electrode relative to the axial skeleton.
Electromyograms (EMGs) were amplified 20 000 times using Grass model P511 K preamplifiers using high-and low-bandpass filter settings of 100 Hz and 3 kHz, respectively, and a 60 Hz notch filter. The analog EMGs were recorded with a TEAC XR-5000 FM data recorder using a tape speed of 9.5 cm s-1. A pulse generator provided coded output to both the NAC video system and the TEAC tape recorder, allowing synchronization of EMGs and video recordings. We converted the analog signal to a digital file using a sampling rate of 8 kHz (Jayne et al. 1990). We then filtered the digital EMGs using a finite impulse response filter that reduced the portion of the signal below 100 Hz to less than 10 % of its original amplitude.
To analyze the digitally filtered data, we used a custom-made computer program to determine the onset and offset times (to a resolution of ±1 ms) and the rectified integrated area of each EMG burst. Relative durations of EMG bursts were calculated by dividing the burst duration by the duration of the locomotor cycle, which was determined from lateral displacement data of the most posterior vertebra.
We used primarily four variables (see Table 2: on-β shift, off-β shift, on-z shift and off-z shift) to describe the timing of muscle activity relative to kinematic data. We subtracted the time of maximal convexity, βmax (towards the side of the fish containing the electrode) from the time of EMG onset, and the resulting difference was converted to a phase shift (on-β shift) by dividing by the cycle duration. In a similar fashion, we calculated the phase shift between EMG offset and the time at which a region was maximally concave towards the side of the electrode (off-β shift). Thus, values of zero for both on-β shift and off-β shift correspond to simple electrical activity during the shortening of contractile tissue. We also calculated the phase shift between the time of maximal lateral displacement (zmax) and EMG onset (on-z shift) and offset (off-z shift). For all lag time and phase shift variables, negative values indicate that the first event preceded the later event in the pair of events in a variable name. For example, a negative value of on-β shift indicates that EMG onset preceded the time of maximal lateral convexity.
Statistical analyses
We sought to obtain a complete set of results from seven longitudinal locations during each of five steady swimming speeds for all four individuals. However, because of the relatively small size of the superficial red muscle in Micropterus salmoides (Jayne and Lauder, 1995b; Johnson et al. 1994), we did not always obtain usable EMGs from every site shown in Fig. 1. Furthermore, a clear periodicity to the kinematic variables was not always evident since the amplitudes at site 1 (Fig. 1) at low speeds were occasionally near the limit of resolution of our movement data (approximately 0.25 ° of intervertebral flexion). Consequently, some descriptive statistics are included for site 1 when possible, but our primary statistical analysis (Table 2) was restricted to sites 2–7.
It is desirable to use statistical methods that account for individual variation, but for sites 2–7, we did not always have measurements for every individual at every speed. Thus, we analyzed our data with a two-way analysis of variance (ANOVA) with speed (five levels) and longitudinal location (six levels) as fixed and crossed factors. To minimize problems associated with pseudo-replication, we use only one mean value (N=3–4 tail beats) for each longitudinal location for each individual at each speed. If we obtained EMGs from both the left and right electrodes at a particular longitudinal location, we pooled left and right observations and used the resulting mean for our two-way ANOVAs. Because several ANOVAs were performed on the same series of results, a variety of methods should be considered for adjusting the criterion for deciding whether a test was statistically significant. We indicate where performing such a correction for multiple comparisons might alter the conclusions from our ANOVAs.
Results
Kinematics
General features of the kinematics of the Micropterus salmoides in this study conformed closely to those described previously for slightly larger individuals (Jayne and Lauder, 1995a). For the five standardized speeds, 0.7, 1.2, 1.6, 2.0 and 2.4 L s-1, mean (N=4) values ± S.D. of the duration of the tail-beat cycle were 519±109, 394±59, 319±35, 278±34 and 247±46 ms, respectively. The corresponding tail-beat frequencies (f) ranged from 1.93 to 4.05 Hz and increased linearly with specific swimming speed (U), as described by the regression f=1.26U+1.05 (r2=0.998, P<0.001).
Lateral displacement (z) and bending (β) (Fig. 2) were propagated posteriorly and varied with time in a sinusoidal fashion. For sites 2–7 (Fig. 1), which were used for the statistical analysis of EMGs, zmax and βmax increased significantly both with increased speed (most commonly from 0.7 to 1.6 L s-1) and with increasingly posterior location (Table 2; Fig. 3A,B). As described previously (Jayne and Lauder, 1995a), zmax and βmax were not in phase, and the values of z-β shift changed significantly with site position (Table 2). When the phase of zmax and βmax for sites 2–7 was calculated relative to the end of the vertebral column, z phase and β phase were not significantly affected by speed (Table 2).
Including kinematic data from the more extreme longitudinal locations of site 1 and the caudal fin (Fig. 1) in ANOVAs (results not shown in Table 2) of the kinematic phase and phase shift variables altered the conclusions regarding the significance of speed. For example, for values of z-β shift from site 1 through the caudal fin site, a two-way ANOVA revealed marginally significant effects of the speed (P=0.03) and of the speed×site (P=0.02) terms, whereas the effect of longitudinal position was highly significant (P<0.001), as had been the case for the results for sites 2–7 summarized in Table 2. When the phases of zmax and βmax were calculated relative to the caudal fin site rather than relative to the end of the vertebral column, the effect of longitudinal position was still highly significant (P<0.001) and the interaction terms were not significant; however, in this case, z phase (P=0.01) and β phase (P=0.002) decreased significantly with increased speed. Thus, analyses including site 1 and the caudal fin consistently showed more extensive effects of speed on kinematics than analyses confined to sites 2–7. However, because the primary goal of this study was to investigate congruent patterns of change in EMGs and kinematics with changing speed and longitudinal position, we concentrated the analysis on sites 2–7.
Electromyography
The fundamental pattern of red muscle activity consisted of posteriorly propagated EMGs that were unilateral and alternated between the left and right sides at a given longitudinal location (Fig. 4). The duration of EMGs (in ms) decreased significantly both with increased speed and with more posterior position within each swimming speed (Table 2). Relative EMG duration (in cycles) did not change significantly with swimming speed, but did decrease significantly with more posterior position (Table 2; Fig. 5). Using the relative EMG durations pooled for all available individuals (Fig. 5) at all five speeds, the mean values ± S.D. for sites 1–7 were 0.448±0.017, 0.451±0.016, 0.462±0.029, 0.432±0.010, 0.419±0.012, 0.403±0.022 and 0.335±0.043 cycles, respectively.
Noticeable variation did occur in the relative EMG durations among different tail beats within individual fish. For example, approximately 20 % of the observations for the relative EMG duration at sites 1–3 ranged from 0.5 to 0.53 cycles. An EMG duration exceeding 0.5 cycles implies some simultaneous activity of contralateral sites. However, when EMGs from a pair of contralateral electrodes were examined, there was never a consistent pattern of bilateral activity. In contrast to the more anterior sites, relative EMG durations at sites 4–7 never exceeded half a cycle. Furthermore, for the most posterior position (site 7), all values of relative duration (with the exception of one) were less than 0.39 cycles.
The relative intensity of red muscle activity increased significantly with increased swimming speed, and the pattern of change in relative intensity with speed was similar for all longitudinal positions (Table 2; Fig. 4). EMG intensity increased by approximately fourfold between swimming speeds of 0.7 and 2.4 L s-1 (Fig. 3C), whereas the amplitudes of kinematic variables only increased by as much as 50 % over the same speed range (Fig. 3A,B). EMG intensity also did not appear to be asymptotically approaching a maximum near 2.4 L s-1 (Fig. 3C). In contrast, the most conspicuous increase in zmax occurred between 0.7 and 1.6 L s-1 (Fig. 3A), and for βmax, most longitudinal positions also showed reduced rates of increase between 2 and 2.4 L s-1 (Fig. 3B).
Despite the fact that muscle activity propagated posteriorly, there were usually periods when the entire longitudinal extent of the red muscle on one side of the fish was simultaneously active (Fig. 6). The longitudinal extent of muscle activity can also be seen in Fig. 7 by observing the number of longitudinal positions for which EMGs (the horizontal bars) intersect a line (representing an instant in time) that is perpendicular to the x-axis. Using the mean values for experimentally observed muscle activity (Fig. 7C–E), the durations of simultaneous ipsilateral activity of sites 1–7 for speeds from 1.6 to 2.4 L s-1 were 20 ms (0.062 cycles), 50 ms (0.179 cycles) and 46 ms (0.247 cycles), respectively. The timing of muscle activity at site 1 could not be determined with confidence for 0.7 and 1.2L s-1 (Fig. 7A,B). Using data from site 2, extrapolated to obtain estimated values for site 1 at the two slowest swimming speeds (Fig. 7A,B), the durations of simultaneous ipsilateral activity for 0.7 and 1.2 L s-1 were 0 and 22 ms (0.055 cycles), respectively. Thus, the proportion of the locomotor cycle with simultaneous ipsilateral muscle activity tended to increase from site 1 to site 7 with increased swimming speed. Simultaneous activity from sites 1–7 on one side of the fish generally began close to the time when the trailing edge of the tail was maximally displaced to the contralateral side (see Fig. 6 at 0.27 cycles) and ended near or before the time when the tail was crossing the axis of forward progression.
The timing of muscle activity relative to kinematics is complicated by the fact that all phase-shift variables relating muscle activity to kinematics varied significantly with both speed and longitudinal position (Table 2: on-β shift, off-β shift, on-z shift and off-z shift). The time difference between EMGs and kinematic events generally showed a progressive posterior increase within each of the swimming speeds (Fig. 8).
If the timing of EMGs corresponded exactly to shortening of contractile tissue (as predicted from lateral bending), then both on-β shift and off-β shift would equal zero. As shown in Fig. 8A,B, anterior EMG onset was nearly synchronous with the time when a region was maximally convex towards the side with the recording electrode (on-β shift=0). Mean values of on-β shift were as large as –0.170 cycles for the most posterior sites at the fastest swimming speed (Fig. 8A), and faster swimming speeds generally had greater on-β shift values. If EMG duration were constant among different longitudinal locations, then one would expect that longitudinal variation in off-β shift would resemble that of on-β shift. However, the longitudinal variation in EMG duration is such that it exaggerates the progressive posterior increase in amplitudes of off-β shift compared with on-β shift. For example, at site 7, mean values of on-β shift ranged from –0.065 to –0.170 cycles, whereas off-β shift ranged from –0.239 to –0.317 cycles.
The phase-shift variables do not indicate whether there is a constant time lag between muscle activity and kinematic events because the cycle duration changes with swimming speed. Consequently, we also analyzed the lags (in ms) between EMGs and lateral bending and displacement. Two-way ANOVAs (not shown in Table 2) demonstrated that all four EMG–kinematic lag variables (on-β lag, off-β lag, on-z lag and off-z lag) had highly significant differences attributable to both swimming speed and longitudinal position, with a pattern of variation grossly similar to that observed for the phase-shift variables (compare Fig. 8A,B with Fig. 8C,D).
As shown in Fig. 8C, mean values of on-β lag ranged from –57 ms to 25 ms, and roughly 90 % of the values for sites 4–7 (0.572–0.774L) ranged from –45 ms to –25 ms. Values of off-β lag ranged from –156 ms to 8 ms, and approximately 90 % of the values for sites 4–7 were less than –50 ms (Fig. 8D). Interestingly, the greatest disparity of EMG phase relative to β occurred for the fastest speed (Fig. 8A,B), whereas the greatest lag times between EMGs and bending generally occurred for the slowest swimming speed (Fig. 8C,D). The phase shifts and time lags between EMGs and lateral displacement showed overall patterns of change associated with swimming speed and longitudinal position similar to those for the timing of EMGs relative to lateral bending (Fig. 7).
Longitudinal propagation of kinematics and EMGs
We were not always able to record EMGs simultaneously from adjacent ipsilateral sites, but we commonly did obtain EMGs from most of the longitudinal positions as a result of having electrodes in both the left and right sides. Thus, to allow pooling of EMG results from the left and the right sides, we added values of on-β shift to those of β phase to obtain the phase of EMG onset, and then, for each successive pair of adjacent longitudinal locations, we obtained the phase difference and divided it by the number of intervening body segments (calculated separately for each individual fish). The resulting intersegmental phase lags (in cycles per body segment) provide a useful indication of the time course of events, and the inverse indicates a standardized rate of propagation (also referred to as EMG wavelength) in terms of numbers of body segments per cycle. Similar calculations were performed for EMG offsets. Calculations of intersegmental phase lags for kinematic variables were simpler (used left-side data only) since there were no missing values.
Neither the intersegmental phase lags (IPL) for lateral displacement (z) nor those for lateral bending (β) changed significantly with swimming speed (Table 3). z IPL did have highly significant variation among longitudinal locations (Table 3), with the smallest values consistently occurring at intermediate longitudinal locations (Fig. 9A). The mean z IPL for all data pooled for sites 2–6 was 0.0195 cycles segment-1. Mean values (pooled across all speeds) of z IPL for sites 1–2 and for sites 6–7 were 0.0289 and 0.0257 cycles segment-1, respectively. Values of β IPL had marginally significant variation with longitudinal position (Table 3) which was no longer significant if a Bonferroni correction was made for multiple comparisons (P=0.003). Values of β IPL generally decreased posteriorly (Fig. 9B). Using pooled values from all swimming speeds, mean values of β IPL from sites 2–3 and sites 5–6 were 0.0255 and 0.0135 cycles segment-1, respectively.
The intersegmental phase lags of EMG onset and offset showed no highly significant (Table 3) or consistent pattern (Fig. 9C,D) of variation with either longitudinal position or swimming speed. Using data pooled for all swimming speeds from sites 2–7 (Fig. 9, N=25), mean ± S.E.M.values of IPL for EMG onset and offset were 0.0182±0.0020 and 0.0102±0.0018 cycles segment-1, respectively. Despite the rather high variability in the estimates of IPL for EMGs (Fig. 9C,D), onset IPL was significantly greater than offset IPL (paired t-test, t=4.18, d.f.=24, P<0.001). Hence, as indicated by the inverse of the mean values of IPL, EMG offsets propagated posteriorly faster than EMG onsets (98 versus 55 segments per cycle).
Similar paired comparisons can be made between values of IPL for EMGs and kinematic events. For the two slowest swimming speeds, smaller sample sizes of EMG intersegmental phase lags contributed to a large among-site variance, which obscured most trends for paired comparisons between EMG and kinematic intersegmental phase lags (Fig. 9). However, when comparisons were restricted to the more precise data from the three fastest speeds, some interesting differences in IPL became apparent. For example, for the three fastest swimming speeds, the values of IPL for β (Fig. 9B) were significantly larger than those for EMG onset (Fig. 9B,versusFig. 9C, paired t-test, t=4.62, d.f.=14, P<0.001). Thus, EMG onset propagated faster than lateral bending at homologous longitudinal locations. For this comparison, the disparity in values of IPL was greatest anteriorly and diminished posteriorly until, near site 6, the values of IPL were nearly identical. Similarly, for the three fastest speeds, the values of IPL for EMG offset (Fig. 9D) were smaller than those of β (Fig. 9B) and the difference between these two quantities increased anteriorly along the fish.
Discussion
Currently, only limited comparative data are available for the intersegmental phase lags of EMGs during the undulatory swimming of vertebrates. For a combination of intact and fictive preparations of anguilliform swimmers, EMG intersegmental phase lags generally approximate 0.01 cycles segment-1 and appear to be unaffected by either longitudinal position or cycle duration (Grillner and Kashin, 1976; Wallén and Williams, 1984). For fictive preparations of dogfish (Grillner, 1974, calculated from his Fig. 5A) and goldfish (Fetcho and Svoboda, 1993), EMG intersegmental phase lags approximated 0.004 and 0.021 cycles segment-1, respectively, and these values and did not vary with frequency of EMGs. Consequently, our study of Micropterus salmoides and all other studies of fish swimming currently support the generalization that EMG intersegmental phase lags are independent of swimming speed.
For the swimming of intact trout, Williams et al. (1989) found that the rate of EMG propagation depended on longitudinal position, with anterior locations having faster speeds of propagation; however, the propagation speeds of lateral bending appeared to be constant along the entire length of the fish. Williams et al. (1989) expressed rates of propagation in the trout in terms of percentage length rather than body segments; therefore, direct comparisons of intersegmental phase lags are not possible for the trout and Micropterus salmoides. Data from studies of trout and carp (Williams et al. 1989; van Leeuwen et al. 1990) suggest that the propagation of EMG offset is faster than that of EMG onset in both taxa (Fig. 10). Hence, bass, trout and carp, which are all subcarangiform swimmers, appear to have different propagation speeds of EMG onset and offset, whereas available data for diverse anguilliform swimmers appear to indicate that EMG onsets and offsets propagate with the same speed (Grillner and Kashin, 1976; Williams et al. 1989; Jayne, 1988).
We designed our experiments specifically to detect effects of speed and longitudinal location on the relationship of EMGs to movement. Consequently, we found effects of speed on certain variables which have previously been found to be constant with speed in addition to finding effects of speed on other variables which have not been previously quantified. For example, for the steady swimming of trout, Williams et al. (1989) found that speed did not affect the timing of EMGs relative to lateral bending, although this relationship did vary longitudinally. In contrast, for Micropterus salmoides, the timing of EMGs relative to lateral bending varied significantly both with speed and with longitudinal location (Table 2), regardless of whether values were expressed as lag times or as phase shifts (Fig. 8). In the caudal regions of Micropterus salmoides, values of on-β shift (Fig. 8A) changed by as much as 0.105 cycles between 0.7 and 2.4 L s-1, while values of on-β lag were generally between –30 and –40 ms (Fig. 8C). Values of off-β shift changed with speed similarly to those of on-β shift. At a given longitudinal location, values of off-β lag were larger than those of on-β lag, and most off-β lags ranged from –50 to –100 ms. For Micropterus salmoides, Johnson et al. (1994) found that the times to peak tension and half-relaxation for isometric twitch were 43 and 67 ms, respectively, for in vitro samples of red muscle removed near sites 6–7 (Fig. 1) and tested at 20 °C. Interestingly, the disparity between the times of EMG onset and the time of maximal convexity was nearly always less than the time to peak tension for the superficial red muscle. For other species of fish that swim using the subcarangiform mode, considerable longitudinal variation may occur for the time course of force development and relaxation. For superficial red muscle of the scup, Rome et al. (1993) found that the relaxation time could increase nearly twofold posteriorly. From the available data on the fast axial musculature of fish, it appears that the times of both force development and relaxation tend to increase posteriorly, but the magnitude of longitudinal changes in relaxation times generally exceeds that of the times required to develop tension (Wardle, 1985; Davies and Johnston, 1993; Altringham et al. 1993). Because these previously described patterns of longitudinal variation in contractile properties are remarkably congruent with the longitudinal patterns of variation which we found for values of on-β lag and off-β lag (Fig. 8C,D), it is tempting to conclude that the timing differences between force and lateral bending are minimal compared with differences in the timing of electrical activity and lateral bending. This suggestion may be further supported by observations of positive in vitro work loops for the superficial red muscle of Micropterus salmoides using strains and stimulation phases representative for a very posterior location (site 6–7) during steady swimming (Johnson et al. 1994). Because not all species of fish appear to show longitudinal variation in contractile properties, and both strain and stimulus have complex effects on force development and relaxation (Johnston et al. 1993), further experiments clarifying the extent of longitudinal variation in the contractile properties of Micropterus salmoides muscle would be very desirable.
The strain of the superficial red muscle of Micropterus salmoides can be estimated by converting values of βmax to a percentage length change, as described previously by Jayne and Lauder (1995a). Values of βmax in the present study were statistically indistinguishable (two-tailed paired t-test, t=0.72, d.f.=15, P=0.48) from those published previously for
Micropterus salmoides swimming steadily unencumbered by electrodes (present study Table 4 using sites 2, 3, 4 and 6 versusJayne and Lauder, 1995a, Table 4). Consequently, our methodology for obtaining EMGs appeared to have few obvious effects on the swimming behavior of Micropterus salmoides over the observed range of swimming speeds. As shown in Table 4, estimated strains of the superficial red muscle were between ±1.9 % and ±6.8 % depending on the exact longitudinal location and swimming speed. Besides agreeing well with our previous observations for Micropterus salmoides (1.7–6.6 %; Jayne and Lauder, 1994a), these estimates of superficial red muscle strain are very similar to values presently available for the steady swimming of a variety of other fish species (van Leeuwen et al. 1990; Rome and Swank, 1992; Rome et al. 1993).
Fairly detailed comparisons are now possible for the timing of red muscle activity and lateral bending for Micropterus salmoides, the trout (Williams et al. 1989) and the carp Cyprinus carpio (van Leeuwen et al. 1990), and this information is summarized in Fig. 10, which shows EMGs and kinematic data plotted on a time scale standardized to equal the duration of a single cycle. The largemouth bass appears to be more similar to the carp for two main reasons. First, in both bass and carp, EMG onsets occur before the time of maximal convexity, whereas in the anterior region of the trout the EMGs lag behind this kinematic event. Second, in both bass and carp, there is a noticeable portion of the locomotor cycle for which all of the red muscle was simultaneously activated on one side of the fish, whereas the greatest longitudinal extent of simultaneous ipsilateral muscle activity in the trout is always less than its body length. As discussed earlier, decreased swimming speed may decrease the longitudinal extent of muscle activity in Micropterus salmoides, and data shown in Fig. 10 are for a trout swimming more slowly than the bass. However, Williams et al. (1989) reported that this pattern of muscle activation was also characteristic of faster swimming speeds. In the light of the substantial effects of speed on EMGs and their relationships to kinematics in Micropterus salmoides, future comparisons will best be accomplished using standardized swimming speeds for different fish taxa.
Similar numbers of vertebrae are found in both bass (30–32) and carp (35–36), whereas trout (60–66) have nearly twice as many (Scott and Crossman, 1973). Consequently, similar numbers of vertebrae may be a key factor contributing to the similar muscular control of swimming movements in bass and carp. Bass, carp and trout may all be considered subcarangiform swimmers, and the overall shape of moderately sized individuals from all three species resembles the generalized fusiform condition (Webb, 1994). Because high numbers of vertebrae are often correlated with elongate shape, these two factors are often confounding when attempting to interpret the diversity of fish swimming modes and their underlying muscular basis. In the light of these encouraging results comparing bass with carp and trout, interpreting the diversity of morphology and swimming modes in fishes should be facilitated by additional studies that attempt to isolate the effects of additional morphological traits (such as variation in body depth) when numbers of vertebrae are similar (e.g. by comparing Micropterus salmoides and Lepomis macrochirus motor and kinematic patterns).
ACKNOWLEDGEMENTS
We are grateful to J. Seigel and the Section of Fishes at the Los Angeles County Museum for the X-ray photographs of specimens. H. Nguyen, J. Davis, A. Lozada and B. Malas laboriously digitized the video images. Support was provided by NSF grants BNS 8919497 to B.C.J. and G.V.L., and NSF BSR 9007994 to G.V.L. The high-speed video system was obtained under NSF BBS 8820664.