Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 rsos.royalsocietypublishing.org Research Cite this article: Vernooij CA, Reynolds RF, Lakie M. 2016 Physiological tremor reveals how thixotropy adapts skeletal muscle for posture and movement. R. Soc. open sci. 3: 160065. http://dx.doi.org/10.1098/rsos.160065 Received: 29 January 2016 Accepted: 5 April 2016 Subject Category: Biology (whole organism) Subject Areas: physiology/behaviour/neuroscience Keywords: physiological tremor, mechanical resonance, thixotropy, posture, electromyography, muscle Author for correspondence: Carlijn A. Vernooij e-mail: [email protected] Physiological tremor reveals how thixotropy adapts skeletal muscle for posture and movement Carlijn A. Vernooij1,2 , Raymond F. Reynolds1 and Martin Lakie1 1 School of Sport, Exercise and Rehabilitation Sciences, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK 2 Institut des Sciences du Mouvement E.J. Marey (UMR 7287), Aix-Marseille Université and CNRS, 163 Avenue de Luminy, CP 910, Marseille 13009, France People and animals can move freely, but they must also be able to stay still. How do skeletal muscles economically produce both movement and posture? Humans are well known to have motor units with relatively homogeneous mechanical properties. Thixotropic muscle properties can provide a solution by providing a temporary stiffening of all skeletal muscles in postural conditions. This stiffening is alleviated almost instantly when muscles start to move. In this paper, we probe this behaviour. We monitor both the neural input to a muscle, measured here as extensor muscle electromyography (EMG), and its output, measured as tremor (finger acceleration). Both signals were analysed continuously as the subject made smooth transitions between posture and movement. The results showed that there were marked changes in tremor which systematically increased in size and decreased in frequency as the subject moved faster. By contrast, the EMG changed little and reflected muscle force requirement rather than movement speed. The altered tremor reflects naturally occurring thixotropic changes in muscle behaviour. Our results suggest that physiological tremor provides useful and hitherto unrecognized insights into skeletal muscle’s role in posture and movement. 1. Background Human motor activity consists of periods of immobility, enlivened by periods of movement. How are these very different postural and dynamic roles addressed by the neuromuscular system? In many species, there is a division of skeletal muscle into tonic and phasic types with profound mechanical differences, for example amphibians [1], insects [2], birds [3] and reptiles [4]. 2016 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited. Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 2.1. Experimental volunteers The experiment was carried out on 15 healthy, right-handed volunteers (23.7 ± 9.9 years old, three female) who gave written informed consent, and was undertaken in accordance with the declaration of Helsinki. Permission was obtained from the ethics committee of the University of Birmingham. None of the volunteers suffered from known neurological or muscular disorders. Participants were asked to refrain from exercising and not to consume alcohol or caffeine 24 h before participation. ................................................ 2. Material and methods 2 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 Often the tonic types are arranged as bundles which parallel the phasic fibres, so that the same muscle can act in different ways depending on which bundle is activated. However, mammals, including humans, do not have specific muscle types for tonic or phasic behaviour. Instead, they draw on a common pool of skeletal muscle motor units, invariably recruited in order of increasing motor neuron size, to satisfy both posture and movement [5]. There are certainly mechanical differences between the initially recruited slow motor units and later recruited faster motor units. One study on human extensor hallucis brevis muscles showed contraction times with an extreme range of 35–98 ms [6]. These differences are actually rather slight, considering that posture may endure for periods of many seconds or minutes, whereas movements may be over in a fraction of a second. Thus, they seem more suited to generating movements of different contraction speeds rather than for the fundamentally opposed roles of movement and posture. However, motor unit size and type are not the only determinants of muscle behaviour. There are profound mechanical differences between muscles when they are lengthening or shortening, and when they are static (or very nearly so). This difference is caused by muscular thixotropy, which can be defined as ‘the dependence of muscle stiffness on the history of length changes’. Thixotropic changes greatly alter the muscle’s mechanical response to artificial stimulation and imposed movement [7–11]. While at rest, muscles become stiff and resistant to lengthening or shortening. If there is a change in muscle length which exceeds the minor elasticity inherent in cross-bridges, this ‘stiction’ falls to a low level, re-establishing itself rapidly as posture returns (the phenomenon was noted as early as 1929 by DennyBrown [12]; he called it stationary rigidity). Thixotropic stiffening will ensure stability at rest, but will also permit unimpeded movement, thus enabling skeletal muscle to fulfil completely different postural and dynamic roles. The implication is clear; muscle can be expected to act very differently when stationary and when moving. Such changes are easily demonstrated by in vitro experiments [13]. However, the changes in muscle behaviour in vivo are relatively inscrutable and have not been much studied. In a few recent studies, we have explored the altered relationship between the input to a muscle, measured as electromyography (EMG), and its output, measured as limb tremor (acceleration), as a consequence of muscle thixotropy (hand: [10,14] and finger: [11,15]). With movement, the peak in the tremor spectrum shifts to a lower frequency, whereas the shape of the EMG spectrum changes only slightly. The explanation is that the stiffness of the muscle changes (decreases), but its inertial load does not, thus altering the tuning of the muscle/limb oscillator and the way that it responds to the buffeting of motor unit firing. However, these studies separately examined posture and movement, and did not track the mechanisms underlying the transition between a period of rest and a period of movement. In this paper, we describe the effect of transitions between posture and movement on muscle properties. Subjects used their middle finger to track a target that evolved from a static position into a very slow vertical sinusoidal movement before returning to a stationary position. We recorded surface EMG from the finger extensor muscle because it represents the neural drive to move the finger to overcome gravitational and elastic forces. We simultaneously recorded finger tremor as an indication of muscle output. We hypothesized the following: (i) tremor frequency will alter as a limb alternates between posture and movement; (ii) the frequency modulation can be parsimoniously explained by mechanical changes in muscle tissue and (iii) EMG magnitude will reflect the load but its frequency composition will not relate to tremor frequency. Our results show that human physiological tremor, sometimes seen as little more than a curiosity, can illuminate the way in which muscle carries out its postural and dynamic duties and provide fundamental insights into the nature of posture and movement. Some of these findings have been briefly reported elsewhere [16,17]. Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 average performance 0 target –15 representative performance 0 20 40 60 time (s) 80 100 120 speed (deg s–1) speed of the target and finger 3.0 average performance 1.5 target 0 0 20 40 60 time (s) 80 100 representative performance 120 Figure 1. The target (thick black line), a representative performance (thin black line) and the average performance of all subjects/all trials (dashed) in the 120 s tracking task. Positive values represent an upward movement of the target and finger (extension). Note, for clarity, the traces have been vertically displaced so they do not superimpose. 2.2. Apparatus The participants sat in a comfortable chair with the right forearm pronated and supported by a plastic curved rest. The hand as well as the index and ring finger were securely taped to a horizontal U-shaped aluminium support with a gap for the middle finger adjusted to fit individual subjects. This allowed for unhindered and isolated flexion–extension of the middle finger around metacarpophalangeal joint 3. A light duraluminium splint underneath the middle finger prevented movement at interphalangeal joints. The arm and the hand rest were individually connected to a heavy steel table by magnetic supports, so they could be optimally positioned. A miniature three-axis accelerometer (model SCA3000, Active Robots, UK, 12.7 × 20.32 mm) was attached above the nail bed of the middle finger to measure its vertical acceleration. A retroreflective laser rangefinder (YP11MGV80, Wenglor Sensoric, Germany) was pointed at a white plastic reflective surface (approx. 2 × 3 cm) placed on top of the accelerometer to record vertical finger position. A computer screen approximately 1.5 m in front of the subject displayed finger position as a white cross and a computer-controlled target in the form of a red ball. The target was a stereotyped positional waveform which included transitions between posture and movement as shown in figure 1. Specifically, the target was initially static at an individually determined comfortable neutral (neither extended nor flexed) position of the finger. After 10 s, the target moved into a vertically orientated sinusoid with a frequency linearly rising from 0 to 0.05 Hz over 50 s (i.e. a chirp signal). The target then decelerated as a mirror image of this chirp and ended stationary in the neutral position for 10 s. The entire sequence occupied 120 s. The peak-to-peak amplitude of the waveform was approximately 30° of finger movement (15° extension (up) and 15° flexion (down)), which was well within the limits of the subjects’ comfortable range of motion. Participants were asked to track the target by keeping the white cross and red ball on the computer screen aligned. It was emphasized to the subjects that this was not a test of accuracy or precision and that they should remain as relaxed as possible. The maximum angular velocity required to track was 3 degree s−1 . Subjectively, this seemed a very slow movement, and this stereotyped task was easy to perform. Each subject repeated the task 10 times and refrained from moving their hand or fingers 10 s before the start of every trial. Surface EMG was recorded from the belly of the extensor digitorum communis muscle (m. EDC) with a Bagnoli system (Delsys Inc., USA). ................................................ +15 3 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 position (deg) position of the target and finger amplitude (×10–5 mV) 10 10 10 8 8 8 8 6 6 6 6 4 4 4 10 20 30 40 50 frequency (Hz) 10 20 30 40 50 frequency (Hz) 4 10 20 30 40 50 frequency (Hz) 5s 4 10 20 30 40 50 frequency (Hz) 5s 5s 5s amplitude (×10–3 m s–2) acceleration 120 s 3.5 3.0 2.5 2.0 1.5 1.0 0.5 10 20 30 40 50 frequency (Hz) 3.5 3.0 2.5 2.0 1.5 1.0 0.5 10 20 30 40 50 frequency (Hz) 3.5 3.0 2.5 2.0 1.5 1.0 0.5 10 20 30 40 50 frequency (Hz) 3.5 3.0 2.5 2.0 1.5 1.0 0.5 10 20 30 40 50 frequency (Hz) Figure 2. Frequency spectra of EMG and acceleration (tremor) at different phases of target tracking by a representative subject. Top: EMG frequency spectra. Middle: finger position. Bottom: acceleration frequency spectra. From left to right, spectra are shown for 5-s snapshots of static posture, maximum upward finger velocity, maximum downward finger velocity and static posture. The tremor spectra are relatively robust, whereas the EMG is much more variable. There is a considerable modulation of tremor frequency and size between posture and maximum velocity snapshots while there is no obvious direct correspondence with the related EMG spectra. 2.3. Data analysis EMG was amplified ×1000. EMG, acceleration and positional data were sampled at 1000 Hz and digitized by an MC 6026 PCI card. Further analysis was carried out offline using custom MATLAB scripts (MathWorks Matlab 2011a, USA). EMG was band-pass filtered (35–200 Hz, fourth-order Butterworth zero-phase-lag filter) and rectified. Acceleration signals were high-pass filtered (0.1 Hz, fourth-order Butterworth dual filter) to correct for artefactual modulation within the range of target frequencies. For one representative subject, the amplitude spectra of the finger acceleration and EMG were obtained by NEUROSPEC software (v. 2.0, 2008) in order to display four 5 s snapshots at different phases of the 120 s target tracking (figure 2); static (second 5–10), maximum upward finger velocity (second 52.5–57.5), maximum downward finger velocity (second 62.5–67.5) and static (second 115–120). Position, EMG and acceleration were then down-sampled to 100 Hz. To capture the non-stationary features of these signals (frequency and power) over the tracking movements, we calculated the wavelet transformation of EMG and acceleration. Wavelets depict the frequency components of a signal over time as if it were a frequency spectrum over a sliding time-window, therefore enabling us to characterize the frequency and power of EMG and acceleration associated with the different movement speeds and positions of the finger. The mean wavelet power was calculated using a continuous wavelet transform (CWT) ([10,18]; for a tutorial, see [19]). The CWT scales and translates a mother wavelet shape at each infinitesimal time step of a continuous signal. The wavelet coefficient C of signal s over time t is calculated √ as: Ca,b = ∫ s(t)(1/ a)Ψ ∗ ((t − b)/a) dt, where ψ is the mother wavelet, a is the frequency scale of ψ, b is the position of ψ and * is the complex conjugation. We used a complex Morlet mother wavelet with a 2 bandwidth f b of 1 Hz and a central frequency f c of 1.5 Hz, so that Ψ (x) = (1/ π fb ) e2tπfc x e−x /fb . The wavelet was scaled to correspond to a frequency range of 5–37.5 Hz. Wavelet analysis allowed us to examine the peak power and frequency of acceleration and EMG over the 120 s duration of the tracking trial. Note that ‘power’ here is relative, as the wavelet analysis transforms the data. Absolute velocity ................................................ 10 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 EMG Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 (b) 0.4 1.5 0.2 0 0 20 40 60 time (s) 80 100 0.6 0.4 y = 0.51x + 0.15 r2 = 0.87 0 1.5 3.0 speed (deg s–1) 0.2 120 3.0 10 1.5 14 18 22 0 postural frequency range 0 20 40 60 time (s) 80 100 120 tremor frequency (Hz) 6 speed (deg s–1) (d) reversed tremor frequency (Hz) (c) 22 y = –13.65*x0.28 + 20.64 r2 = –0.92 18 14 10 6 0 1.5 3.0 speed (deg s–1) Figure 3. (a) Finger acceleration (tremor) power (arbitrary units, solid line ± s.d.) as depicted by wavelet analysis is shown alongside movement speed (dashed line). (b) Tremor power plotted as a function of movement speed. A linear regression provides a satisfactory fit, although it may misrepresent the relationship when movement is slow. (c) Tremor frequency (solid line ± s.d.) as depicted by wavelet analysis and movement speed (dashed line). The postural tremor frequency range is also indicated (details in text). Note that the scale for frequency runs from 6 to 22 Hz and is reversed to emphasize the striking correspondence with movement speed. (d) Tremor frequency plotted as a function of movement speed. A power function provides an excellent fit. 0.02 0 0.01 0 20 40 60 80 100 –15 120 EMG frequency (Hz) (c) 0.03 0.02 y = 4.67e–4x + 0.01 r2 = 0.93 0.01 –15 0 +15 (d) 15 14 13 12 11 0 postural frequency range 20 40 60 time (s) 80 100 120 EMG frequency (Hz) EMG power (arb. units) +15 EMG power (arb. unis) (b) 0.03 position (deg) (a) 15 14 13 12 11 y = 0.03x + 13.10 r2 = 0.34 –15 0 +15 position (deg) Figure 4. (a) EMG power (solid line ± s.d.) as depicted by wavelet analysis and movement speed (dashed line). (b) EMG power plotted as a function of finger position. A linear regression provides an excellent fit. (c) EMG peak frequency as depicted by wavelet analysis. The postural frequency range is also shown (details in text). Note that the modulation of frequency is different when the finger is rising and when it is falling. Also note that the EMG peak frequency does not correspond to the tremor frequency as seen in figure 3. (d) The relationship of EMG peak frequency to finger position. The regression line shows that EMG frequency is on average 0.81 Hz greater when the finger moves from −15 to +15°, but the relationship is weak and the coefficient of determinance (r2 ) is low. (speed) of the finger was obtained by differentiating the positional signal (second-order Savitzky–Golay filter, time window = 151 data points). 2.4. Statistics We wanted to know whether, and at which time point during the transition between posture and movement, the frequency of acceleration and its neural drive would significantly change. Therefore, the ................................................ 0.6 5 0.8 tremor power (arb. units) 3.0 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 tremor power (arb. units) 0.8 speed (deg s–1) (a) Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 3.1. Tracking task performance Figure 1 shows the position of the target waveform with the average tracking performance and a representative tracking performance. The associated speed of target and finger are also plotted. It is clear from this figure that the two were, on average, very well matched. This was to be expected for a slow, repetitive and predictable target movement [20]. Small differences between the target and the response are clearly visible in the position and speed records for representative single trials. Figure 2 shows the frequency spectra for EMG (top) and tremor acceleration (bottom) of one representative subject at four phases of the tracking; static, maximum upward velocity, maximum downward velocity and static, respectively. The corresponding finger position is shown also (middle). Whereas there is a large modification of the tremor spectra in both size and frequency between posture (left and right panels) and movement (middle panels), the differences in EMG are much smaller between the different tracking phases. 3.2. Relationship between acceleration (tremor) and finger movement Figure 3 illustrates the evolution of wavelet-derived tremor (acceleration) power (figure 3a) and frequency (figure 3c) over time during the tracking task (average ± s.d.). The s.d. is very small and hard to detect, depicting the small variance between subjects and trials. There is a striking resemblance between the two acceleration variables and the superimposed average finger speed. Acceleration exhibits an increase in power and decrease in frequency almost immediately after movement onset, and covers a wide range of frequencies. When finger speed exceeds 0.15 degree s−1 the frequency drops below the postural frequency bands (figure 3c; CI = 16.6–20.6 Hz) shown in grey. At the fastest speed (3 degree s−1 ), acceleration frequency reduces to 7.8 Hz. When the finger speed drops towards zero for short periods, the tremor frequency just re-attains the postural frequency band. When the ongoing movement is faster, it fails to re-attain the postural frequency band. Only when the finger is stationary for some time towards the end of the experiment does the frequency definitely re-enter the postural band. The right-hand panels of figure 3 display the relationship between finger speed, and acceleration power (figure 3b) and frequency (figure 3d) averaged over 0.5 s bins. There is a positive linear relationship between power and speed (average individual r2 = 0.64 or r2 = 0.87 for mean signals) and negative power-relationship between frequency and speed (average individual r2 = −0.41 or r2 = −0.92 for mean signals). 3.3. Relationship between EMG (neural input) and finger movement There is a stark contrast between acceleration and EMG, both in terms of the overall range of frequencies covered and the way that they respond to movement. Figure 4a shows EMG wavelet power almost exactly superimposed upon finger position, and the two signals are strongly linearly related (figure 4b; average individual r2 = 0.69 or r2 = 0.93 for mean signals). Not surprisingly, more EMG in the extensor muscle is recorded when the finger moves to a correspondingly higher (more extended) position. Figure 4c shows that, in complete distinction to peak acceleration frequency described above, EMG frequency during movement remains mostly within the postural frequency range (CI = 12.5–15.1 Hz). This relatively small range covers the middle frequencies of those covered by the acceleration (7.8–20.6 Hz). Relatively, steep drops in EMG frequency occur when the finger moves downwards, but only once is there a short period where the frequency significantly (i.e. average ± CI) decreases below the postural range. Overall, there is a slight increase in EMG frequency with higher finger positions (figure 4d; average individual r2 = 0.06 or r2 = 0.34) for mean signals. ................................................ 3. Results 6 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 95% confidence interval (CI) was calculated for the peak frequency of acceleration and EMG over time. For both variables, the frequencies covered by the CI over the last 5 s of the two static periods (at the start and end of the movement) were designated the ‘postural’ frequency range (figures 3c and 4c). We then measured the time when the limits of the CI during the tracking movement would exceed this range, of which the first and last time points were used to signify a borderline finger movement that signifies a change in frequency. Wavelet data as well as finger position and finger velocity were averaged over trials and subjects and finally smoothed (moving average with time constant = 0.4 s). Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 EMG frequency (Hz) 0 13 12 11 y = 53.13x + 12.42 r2 = 0.13 0.01 0.02 EMG power (arb. units) EMG power (arb. units) (c) 0.03 –15 deg 3 m s–2 2.5 0.02 2.0 1.5 1.0 0.01 0.2 0.4 0.6 0.8 acceleration power (arb. units) 0.5 0 m s–2 18 3.0 m s–2 y = 3.17x–0.64 + 4.45 r2 = –0.78 2.5 2.0 1.5 14 1.0 10 0.5 6 0.2 0.4 0.6 0.8 acceleration power (arb. units) (d ) line of unity 15 0 m s–2 3.0 m s–2 2.5 14 2.0 13 1.5 1.0 12 11 y = 0.07x + 12.35 0.5 r2 = 0.13 0 m s–2 6 10 14 18 22 acceleration frequency (Hz) Figure 5. In (a), the greyscale represents finger position. In (b–d), it represents finger speed. Black always represents the lowest values. (a) EMG peak frequency as a function of EMG power. The highest finger positions are associated with the greatest EMG power and the highest peak frequencies. (b) Acceleration frequency as a function of acceleration power. The highest finger speeds are strongly associated with the greatest acceleration power and the lowest frequencies. (c) EMG power as a function of acceleration power. The highest finger speeds are associated with the greatest acceleration power, but there is no strong relationship between speed and EMG power or acceleration power and EMG power. (d) The relationship between EMG peak frequency and acceleration frequency. The highest finger speeds are associated with the lowest frequency of both acceleration and EMG but the two variables alter by very different amounts as indicated by the line of unity (dashed line). 3.4. Power–frequency relationship within EMG and acceleration The average relationship between power and frequency for EMG and for acceleration is shown in figure 5. In each panel, finger kinematics are represented by the use of greyscale (finger position in a; finger speed in b–d), where a fade to black implies lower values. Figure 5a displays a significant positive relationship between EMG power and frequency. In general, as EMG power increases, so does its frequency (r2 = 0.13, p < 0.01), but there is considerable variability. The superimposed greyscale captures the much stronger positive relationship between EMG power and finger position (also shown in figure 4b). Figure 5b shows a significant negative relationship between power and frequency of finger acceleration. The largest drop in frequency occurs when the finger begins to move, i.e. at low finger speeds where the acceleration is small (denoted by darker black dots). A power-curve fit gave a strong correlation of r2 = −0.78 (p < 0.01). As any acceleration has to be a consequence of some neural input, whether broad-band or frequencyspecific input, we also studied the relationship between EMG and acceleration by correlating their power and frequency. Figure 5c shows an expected tendency for more acceleration with increased EMG power, but the strength of this correlation is weak (r2 = 0.18, p < 0.01) and the relationship is complicated; acceleration power can be small when EMG is both high and low. This corresponds to the situation when the finger is stationary at its top or bottom position, respectively. The greyscale shows the previously described positive correlation between finger speed and acceleration power (seen in figure 3b). Figure 5d shows no strong modulation of acceleration frequency with EMG frequency. Although a significant positive linear relationship does exist (solid line; r2 = 0.13, p < 0.01), this line is very different from the line of unity (dashed line) indicating that there is no direct causal relationship. Again, the strong negative relationship between acceleration frequency and finger speed is obvious in the greyscale. 7 ................................................ 14 22 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 15 acceleration frequency (Hz) (b) 15 deg EMG frequency (Hz) (a) Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 4. Discussion The power and frequency of physiological finger tremor (measured as acceleration) are strongly modulated by the speed of finger movement, with faster movement being associated with progressively larger and slower oscillations. The change is obvious in the initial stages of movement; acceleration frequency drops substantially and power increases as soon as movement begins (figure 3). This is not true for the associated EMG power which correlated nearly linearly with finger position (figure 4a,b). Throughout the 120 s trial, the EMG frequency stays within a relatively narrow range (12.5–15.1 Hz), which is situated neither at the high nor low end of the range of acceleration frequencies (7.8–20.6 Hz; figures 3 and 4). There is no simple relationship between the frequency of the acceleration (which alters a lot) and the frequency of the EMG (which alters only a little). This suggests that finger tremor under normal, healthy circumstances might not necessarily be generated by a specific neural input. This gives rise to an important question: are the characteristics of finger tremor determined by the mechanical properties of the limb? The strong correlation of extensor EMG power with finger position (figure 4) was entirely expected because it reflected the increased muscle force required to raise the finger. The activity of the main flexor muscle is not reported here as it hardly exceeded the background level and any activation present was erratic. The maximal velocity (which was attained slowly) was only 3 degree s−1 , which is extremely slow. Therefore, for the speeds that we studied the braking forces were almost certainly passive and mechanical. The downward force was mainly applied by elastic forces in the flexor muscles owing to passive stretch or, more likely, gravity. Some modulation of EMG frequency was also anticipated. For several reasons, there is a generalized increase in EMG frequency as force is increased (reviewed by [23]), for instance, because EMG will be modulated according to rate coding. EMG frequency drops abruptly as the finger commences downwards movement (flexion) and rises more slowly and progressively as the finger is extended (figure 4c). This is probably owing to the different efficiency of the extensor muscle as it acts concentrically and eccentrically. In addition, some EMG frequency modulation could have been caused by measurement artefacts owing to a change in muscle fibre orientation relative to the EMG electrode [23]. Independent of its cause, and more importantly here, there were no corresponding hysteretic effects in tremor acceleration frequency (figure 3c). This divergent behaviour is additional simple evidence against the idea that physiological tremor frequency is mainly a reflection of the underlying frequency of the EMG. 4.2. The resonant nature of physiological tremor In general, physiological finger tremor is described as being composed of a number of acceleration frequency components of distinct origins. Often a combination is proposed of two centrally determined modes of oscillation (a postural tremor generator and a kinetic tremor generator) or a combination of mechanical limb properties and neural oscillations of central origin [24,25] or reflex origin [26,27]. However, our previous research suggested that all frequencies reported in finger tremor could be economically explained by resonance [11,15,28] which differs in frequency depending on the state of the musculature. When the finger and muscle do not move, muscle stiffness is high, thereby generating ................................................ 4.1. Acceleration and EMG changes with movement rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 Here, we examined the behaviour of physiological finger tremor during transitions between posture and movement. This is important, because tremor during (slow) movement has been claimed to be fundamentally different from postural tremor [21]. However, although the relationship of the phase of tremor with rapid voluntary movement has been studied before, the way in which postural tremor at rest evolves into the kinetic tremor of movement, and vice versa, has not been characterized. Our results show that there are predictable (large) alterations in the power and frequency of tremor as posture becomes movement. These alterations are not shown in the concomitantly recorded EMG which drives the muscle. There has been a tendency in the literature to discuss ‘10 Hz tremor’ as though it provides ‘a window into the operation of the nervous system’ [22]. What this study demonstrates is that systematic and quite large alterations in tremor frequency can be generated by movement, and it is not easy to explain these changes as being produced by a neurogenic oscillator. On the other hand, it is quite easy to attribute them to mechanical causes. We discuss whether alterations in the stiffness of the muscles, which occur during movement, produce the altered tremor power and frequency and whether physiological tremor can provide practical and previously unrecognized insights into the stiffness of the musculature. 8 Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 The finding that finger tremor is modulated by the speed of finger movement ties in neatly with a mechanical resonance origin. It is known that extrafusal and intrafusal muscular stiffness is dependent on the recent history of movement [8,29]. When muscle moves, its stiffness reduces greatly, for example by a factor of approximately 15 in the human calf muscles [30]. Although still currently debated in the literature, conceivably the underlying mechanism consists of the elastic characteristic of crossbridges and the number of actin–myosin attachments [8,13,31,32] with the additional contribution of unfolding of gap–filament proteins like titin [33,34]. Whatever the mechanism, during posture, muscle stiffness is high. During movement, the postural short-range elastic stiffness (SREC) is transformed into a smaller, approximately constant, frictional resistance [31]. If a steady posture is subsequently maintained, stiffness will slowly regenerate over the following seconds. This phenomenon, known as thixotropy, makes muscle stiffness very strongly dependent on length changes. With movement, there is substantially decreased muscle stiffness and thus a decreased resonant frequency. The muscle length, resulting from a different finger position, is not the determining factor, because thixotropic stiffening occurs at all muscle lengths [13]. In a reduced preparation (single muscle or single muscle fibre), the transition from SREC to frictional behaviour is very abrupt [13]. In a limb which is controlled by different synergistic and antagonistic muscles and very many muscle fibres, the transition will be expected to be less abrupt and this is what we report here. However, it is clear that the greatest reduction in frequency occurs for very small movements (figure 3d). This behaviour is entirely consistent with a reduction in stiffness once a very small range of movement is exceeded. Hill [31] described the range of the SREC in vitro as 0.2% of muscle length. Several studies have shown that in vivo the required muscular movement to exceed the SREC is very little [8,35]. The tiny movements associated with ‘posture’, i.e. the tremor itself, are accommodated within the SREC range. In contrast, even with very slow movements, the SREC is continually exceeded, and the frequency remains lower than the postural range. This is exactly what we find here (figure 3c). 4.4. The implications for physiological tremor recording These results emphasize the importance of controlling very carefully the conditions under which tremor is recorded. Often, twin peaks in the postural tremor spectrum are discussed in the literature, e.g. [36,37], but this does not necessarily imply they occur simultaneously as is always assumed. When holding the finger extended the muscle stiffness and resonant frequency will be high. However, an occasional small postural adjustment will elicit a temporary drop in muscle stiffness, which will generate a temporary low resonant frequency. A Fourier transform displays the average frequency spectrum of a tremor recording, typically lasting 30–60 s. If small postural adjustments are interspersed with posture, both high and low frequencies will be present in the tremor FFT, both produced by the same mechanism. If the Fourier transform of postural tremor is taken over a few seconds only (as in figure 2), the low frequency peak may be absent. 4.5. What is posture, what is movement? Hill [31] and others have shown very clearly that the threshold of the SREC is one of position. The SREC will re-exert itself at any muscle length if the muscle is allowed to rest there for sufficient time. Rest, presumably, is not absolute, but is in reality a period where velocity is very low. There is no such thing as a fixed position, because there is always minor movement, positional drift and tremor. Figure 3 shows that the confidence limits of postural tremor frequency are exceeded when movement speed reaches a speed of 0.15 degree s−1 . Although this is by definition a movement, its speed is so slow that it is not greatly above the threshold for detection if it was to be passively applied [38]. The postural frequency range is only confidently re-entered when the finger is stationary for some time (110 s onwards; figure 3c). This strongly suggests posture and movement are quantitatively, rather than qualitatively, different. ................................................ 4.3. Thixotropic stiffening of muscle raises resonant frequency 9 rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 a high-frequency component greater than 20 Hz. On the other hand, when muscle moves its stiffness is low and generates a frequency component approximately 8 Hz. (This concept is further discussed below.) This study adds to that research by showing clearly that acceleration power and frequency are modulated by the speed of movement. The smooth modulation and lack of substantial change in the EMG suggest there is a single peripheral determinant of frequency that adjusts its properties with movement. Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 4.6. Is muscle thixotropy the basis of posture and movement? 10 References 1. Sommerkamp H. 1928 Das substrat der Dauerverkfirzung am Froschmuskel (Physiologische und pharmakologische Sonderstellung bestimmer Muskelfasern. Arch. Exp. Path. Pharmak. 128, 99–115. (doi:10.1007/BF01863244) 2. Cochrane DG, Elder HY, Usherwood PNR. 1972 Physiology and ultrastructure of phasic and tonic skeletal muscle fibres in the locust, Schistocerca gregaria. J. Cell Sci. 10, 419–441. 3. Ginsborg BL. 1960 Some properties of avian skeletal muscle fibres with multiple neuromuscular junctions. J. Physiol. 154, 581–598. (doi:10.1113/ jphysiol.1960.sp006597) 4. Gleeson TT, Putnam RW, Bennett AF 1980 Histochemical, enzymatic, and contractile properties of skeletal muscle fibers in the lizard Dipsosaurus dorsalis. J. Exp. Zool. 214, 293–302. (doi:10.1002/jez.1402140307) 5. Bawa PNS, Jones KE, Stein RB. 2014 Assessment of size ordered recruitment. Front. Hum. Neurosci. 8, 532. (doi:10.3389/fnhum.2014.00532) 6. Sica REP, McComas AJ. 1971 Fast and slow twitch units in a human muscle. J. Neurol. Neurosurg. Psychiatry 34, 113–120. (doi:10.1136/jnnp.34.2.113) 7. Lakie M, Robson LG. 1988 Thixotropic changes in human muscle stiffness and the effects of fatigue. Q. J. Exp. Physiol. 73, 487–500. (doi:10.1113/expphysiol. 1988.sp003169) 8. Proske U, Morgan DL, Gregory JE. 1993 Thixotropy in skeletal muscle and in muscle spindles: a review. Prog. Neurobiol. 41, 705–721. (doi:10.1016/03010082(93)90032-N) 9. Axelson HW, Hagbarth KE. 2001 Human motor control consequences of thixotropic changes in muscular short-range stiffness. J. Physiol. 535, 279–288. (doi:10.1111/j.1469-7793.2001. 00279.x) 10. Reynolds RF, Lakie M. 2010 Post-movement changes in the frequency and amplitude of physiological tremor despite unchanged neural output. J. Neurophysiol. 104, 2020–2023. (doi:10.1152/jn.00513.2010) 11. Vernooij CA, Lakie M, Reynolds RF. 2014 The complete frequency spectrum of physiological tremor can be recreated by broad-band mechanical or electrical drive. J. Neurophysiol. 113, 647–656. (doi:10.1152/jn.00519. 2014) 12. Denny-Brown D. 1929 On the nature of postural reflexes. Proc. R. Soc. Lond. B 104, 252–301. (doi:10.1098/rspb.1929.0010) 13. Campbell KS, Lakie M. 1998 A cross-bridge mechanism can explain the thixotropic short-range elastic component of relaxed frog skeletal muscle. J. Physiol. 510, 941–962. (doi:10.1111/j.1469-7793. 1998.941bj.x) 14. Lakie M, Vernooij CA, Osborne TM, Reynolds RF. 2012 The resonant component of human physiological hand tremor is altered by slow voluntary movements. J. Physiol. 590, 2471–2483. (doi:10.1113/jphysiol.2011.226449) 15. Vernooij CA, Reynolds RF, Lakie M. 2013 A dominant role for mechanical resonance in physiological finger tremor revealed by selective minimisation of voluntary drive and movement. J. Neurophysiol. 109, 2317–2326. (doi:10.1152/jn.00926.2012) 16. Vernooij CA, Lakie M, Reynolds RF. 2013 Physiological finger tremor size reflects altered mechanical properties of muscle resulting from changes in neural control. In Abstracts of Progress in Motor Control IX, 13–16 July 2013, Montréal, CA, pp. 73–74. ................................................ Ethics. Ethical approval was obtained from the ethics committee of the University of Birmingham. Volunteers gave written informed consent, and the experiment was undertaken in accordance with the Declaration of Helsinki. None of the volunteers suffered from known neurological or muscular disorders. Participants were asked to refrain from exercising and not to consume alcohol or caffeine 24 h before participation. Data accessibility. All raw data are available through the Dryad repository: http://dx.doi.org/10.5061/dryad.r1j0s. Authors’ contributions. C.V., R.R. and M.L. contributed substantially to conception and design, or acquisition of data, or analysis and interpretation of data; C.V., R.R. and M.L. drafted the article, or revised it critically for important intellectual content; and C.V., R.R. and M.L. have read and approved the manuscript before submission. Competing interests. We have no competing interests. Funding. No external funding is received for this study. Acknowledgements. Thanks to Ken Dawkins and Steve Allen for technical support. rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 The question of the distinctiveness of movement and posture is an old one, with some authorities, but not all, suggesting that specialized neural pathways are involved [39]. Our results do not answer that question. However, it is undeniable that in man the ‘final common pathway’ for both movement and posture is a single set of skeletal muscles which must satisfy the very different requirements of both roles. Muscle stiffness increases greatly as muscle lengthening or shortening ceases, and our results suggest that tremor provides a window into this natural behaviour. They show that as muscle transitions from lengthening or shortening to posture, tremor frequency increases, reflecting an increase in muscle stiffness and resonant frequency. This thixotropic stiffening of muscle is due to the reformation of actin–myosin attachments in combination with non-cross-bridge mechanisms. It is likely that this has a stabilizing influence, making it easy to minimize speed for long periods of time. By definition, this is posture. Man-made servo systems for controlling manufacturing processes, elevator systems, valves and the like, very commonly incorporate an electromechanical brake, so that a static position can be economically maintained. Thixotropic muscle stiffening, which permits the economical and effective control of posture and movement by a common actuator, may be nature’s version of the brake. Meanwhile, with little movement, thixotropic loosening induces very large reductions in muscle stiffness (up to 15-fold; [35]), allowing unimpeded and rapid motion. An example where the functional impact of thixotropy is clear is when raising a gun to shoot a target. This task generates a conflict between movement and stability. The marksman’s muscles are not very stiff when the gun is raised, to enable a large angular momentum. His muscles then get very stiff when he takes aim to maintain joint orientation. The transition in muscle stiffness solves the conflicting requirements of flexibility and speed during fast action, and accuracy and stability at rest. Thus, skeletal muscle is optimized for not only posture, but also for movement. Downloaded from http://rsos.royalsocietypublishing.org/ on May 5, 2017 32. Altman D, Minozzo FC, Rassier DE. 2015 Thixotropy and rheopexy of muscle fibers probed using sinusoidal oscillations. PLoS ONE 10, e0121726. (doi:10.1371/journal.pone.0121726) 33. Kellermayer MSZ et al. 2008 Muscle thixotropy: more than just cross-bridges? Response to comment by Campbell and Lakie. Biophys. J. 94, 329–330. (doi:10.1529/biophysj.107. 122309) 34. Powers K, Schappacher-Tilp G, Jinha A, Leonard TR, Nishikawa K, Herzog W. 2014 Titin force is enhanced in actively stretched skeletal muscle. J. Exp. Biol. 217, 3629–3636. (doi:10.1242/jeb.105361) 35. Loram ID, Maganaris CN, Lakie M. 2007 The passive, human calf muscles in relation to standing: the short range stiffness lies in the contractile component. J. Physiol. 584, 677–692. (doi:10.1113/ jphysiol.2007.140053) 36. Daneault J-F, Carignan B, Duval C. 2011 Finger tremor can be voluntarily reduced during a tracking task. Brain Res. 1370, 164–174. (doi:10.1016/j. brainres.2010.11.047) 37. Stiles RN, Randall JE. 1967 Mechanical factors in human tremor frequency. J. Appl. Physiol. 23, 324–330. 38. Fitzpatrick RC, McCloskey DI. 1994 Proprioceptive, visual and vestibular thresholds for the perception of sway during standing in humans. J. Physiol. 478, 173–186. (doi:10.1113/jphysiol.1994.sp020240) 39. Kurtzer I, Herter TM, Scott SH. 2005 Random change in cortical load representation suggests distinct control of posture and movement. Nat. Neurosci. 8, 498–504. (doi:10.1038/nn1420) 11 ................................................ 25. Vaillancourt DE, Newell KM. 2000 Amplitude changes in the 8–12, 20–25, and 40 Hz oscillations in finger tremor. Clin. Neurophysiol. 111, 1792–1801. (doi:10.1016/S1388-2457(00)00378-3) 26. Christakos CN, Papadimitriou NA, Erimaki S. 2006 Parallel neuronal mechanisms underlying physiological force tremor in steady muscle contractions of humans. J. Neurophysiol. 95, 53–66. (doi:10.1152/jn.00051.2005) 27. Hagbarth KE, Young RR. 1979 Participation of the stretch reflex in human physiological tremor. Brain 102, 509–526. (doi:10.1093/brain/102. 3.509) 28. Lakie M, Vernooij CA, Osler CJ, Stevenson AT, Scott JPR, Reynolds RF. 2015 Increased gravitational force reveals the mechanical, resonance nature of physiological tremor. J. Physiol. 593, 4411–4422. (doi:10.1113/JP270464) 29. Proske U, Gandevia SC. 2012 The proprioceptive senses: their roles in signaling body shape, body position and movement, and muscle force. Physiol. Rev. 92, 1651–1697. (doi:10.1152/physrev.00048. 2011) 30. Loram ID, Lakie M. 2002 Direct measurement of human ankle stiffness during quiet standing: the intrinsic mechanical stiffness is insufficient for stability. J. Physiol. 545, 1041–1053. (doi:10.1113/ jphysiol.2002.025049) 31. Hill DK. 1968 Tension due to interaction between the sliding filaments in resting striated muscle. The effect of stimulation. J. Physiol. 199, 637–684. (doi:10.1113/jphysiol.1968.sp 008672) rsos.royalsocietypublishing.org R. Soc. open sci. 3: 160065 17. Vernooij CA, Reynolds RF, Lakie M. 2013 Physiological finger tremor size mirrors speed of finger movement due to muscle thixotropy. In Proc. of the Physiology Society 7 : 37th Congress of IUPS 2013, 22–26 July 2013, Birmingham, UK, pp. PCD252. 18. Gilbertson T, Lalo E, Doyle L, Di Lazzaro V, Cioni B, Brown P. 2005 Existing motor state is favored at the expense of new movement during 13-35 Hz oscillatory synchrony in the human corticospinal system. J. Neurosci. 25, 7771–7779. (doi:10.1523/ JNEUROSCI.1762-05.2005) 19. Samar VJ, Bopardikar A, Rao R, Swartz K. 1999 Wavelet analysis of neuroelectric waveforms: a conceptual tutorial. Brain Lang. 66, 7–60. (doi:10.1006/brln.1998.2024) 20. Poulton EC. 1974 Tracking skill and manual control. New York, NY: Academic Press. 21. Vallbo ÅB, Wessberg J. 1993 Organization of motor output in slow finger movements in man. J. Physiol. 469, 673–691. (doi:10.1113/jphysiol.1993. sp019837) 22. McAuley JH, Marsden CD. 2000 Physiological and pathological tremors and rhythmic central motor control. Brain 123, 1545–1567. (doi:10.1093/brain/ 123.8.1545) 23. Farina D, Merletti R, Enoka RM. 2004 The extraction of neural strategies from the surface EMG. J. Appl. Physiol. 96, 1486–1495. (doi:10.1152/japplphysiol. 01070.2003) 24. Bye RT, Neilson PD. 2010 The BUMP model of response planning: intermittent predictive control accounts for 10 Hz physiological tremor. Hum. Mov. Sci. 29, 713–736. (doi:10.1016/j.humov.2010.01.006)