A method for determination of muscle fiber diameter using single fiber potential (SFP) analysis

We have used computer simulation to study the relationship between the muscle fiber diameter and parameters: peak-to-peak amplitude and duration of the negative peak of the muscle fiber action potential. We found that the negative peak duration is useful in the determination of fiber diameter via the diameter dependence of conduction velocity. We have shown a direct link between the underlying physiology and the measurements characterizing single fiber potential. Using data from simulations, a graphical tool and an analytical method to estimate the muscle fiber diameter from the recorded action potential has been developed. The ability to quantify the fiber diameter can add significantly to the single fiber electromyography examination. It may help study of muscle fiber diameter variability and thus compliment the muscle biopsy studies.


Introduction
Single fiber EMG (SFEMG) is a powerful technique to study the pathophysiology of the motor unit. Two types of measurements are made using this technique: fiber density (FD) and jitter. The FD is useful to study the grouping of muscle fibers of a motor unit within its territory. In neuropathy, FD increases due to reinnervation. The ''jitter'' measurements are used to assess the efficacy of the neuromuscular transmission [9].
In myopathy, the FD may be increased slightly due to fiber splitting, regeneration of muscle fibers, etc. The jitter is usually normal. In this manner, SFEMG is not particularly useful in the diagnosis of myopathy. It is often used to rule out other diseases such as a neuropathy or a neuromuscular junction disease [9]. A primary change in myopathy is the increased variability of the muscle fiber diameter. This is assessed quantitatively using muscle biopsy studies [2]. There is no electrophysiologic technique to assess this characteristic of the motor unit.
The abnormalities of muscle fiber diameter are observed indirectly on electrodiagnostic studies. On routine needle EMG examination, increased variability gives motor unit potentials with polyphasic waveforms. Atrophy can give low amplitude (in what follows by amplitude we mean the peak-to-peak value of potential change, counted in mV) potentials [1,8,9].
The variability of muscle fiber diameter may be studied indirectly by investigating the muscle fiber conduction velocity. The muscle fibers are stimulated directly using an intramuscular needle, and their action potentials are recorded at a few millimeter distance. The latency of the potential is used to compute the velocity [10].
In principle, the shape of the single muscle fiber action potential also contains information about the muscle fiber diameter. In the so-called line source mode, the extracellular muscle fiber action potential V(t) is computed as [3]: where i(t) is the transmembrane current and u is a weight function. The transmembrane current is proportional to the square of the fiber diameter [3]. Thus, larger fibers will produce a larger amplitude potential. The amplitude is also affected by the distance of the fiber from the electrode [4].
Higher the distance, lower is the amplitude. Hence, amplitude alone cannot be used as a marker of the fiber diameter. The weight function is the potential recorded by the electrode as a unit current source propagates from the endplate to the tendon. The waveform of this function is thus dependent on the propagation velocity of muscle fiber, and hence the fiber diameter. The change in weight function waveform with radial distance and with fiber diameter has not been investigated systematically.
In this study, we have used computer simulation to study the relationship between the muscle fiber diameter and parameters of potential: peak-to-peak amplitude and duration of the negative peak of SFP. These relationships are used to develop a graphical and an analytical tool to estimate the muscle fiber diameter from the recorded action potential. Although our goal is similar, our method of analysis is quite different from that used by Rodriguez et al. [6,7]. They used a computer model to recursively obtain the best match for measured waveform. This is not always possible in clinical environment. We believe that the graphical method described in this study offers more simplicity and gives a better understanding of the action potential waveform.

Methods
The line source model described by Nandedkar and Stålberg [3] was used for simulating single muscle fiber action potentials. This model has been tested in simulations of normal and abnormal EMG signals, and shows a good concordance with experimental and clinical recordings [4,5].
The single fiber potential V recorded at the electrode at time t, coming from a single fiber, is given by a convolution of a weight function u(t) and current i(t) given by (1). The effect on the potential V of the distance of the electrode from the fiber is described by the weight function, which is defined as follows [3]: where r r is the radial conductance and K depends on the ratio of axial to radial conductance, r and z are the radial and axial distance of electrode from the current source. The axial distance z is related to time t by the speed of propagation of the potential along the fiber: z = vt, where the velocity of propagation v is linearly dependent on fiber diameter. The action potentials were simulated for various combinations of muscle fiber diameter (25-90 lm) and radial distance (50-500 lm). The range of fiber diameters corresponds to the data from muscle biopsy [2]. In normal muscle the mean value is 50 ± 5 lm and in various neuromuscular disorders it is in range from 25 to 90 lm. The range of radial distances was determined so that the largest fibers would not be closer to the electrode than their radius (50 lm). The electrode records mainly from muscle fibers that are within 400 lm of the recording surface [4]. Hence the maximum radial distance was limited to 500 lm. The velocity of propagation of the current is determined by the diameter of the fiber according to the following formula [3]: v ¼ 2:2 þ 0:05ðd À 25Þ, where the fiber diameter d is in lm and the propagation velocity is in ms -1 .
For the simulation we have used our own software which is based on the formulation of the simulation method presented by Nandedkar and Stålberg [3]. In this method the potential is a result of the convolution of the weight function and source current given by (1) with the source current described by the formula derived by Nandedkar and Stålberg [3]. Using our software we are able to model the potential from one or more fibers and determine its parameters. In the modeling we have assumed the same parameters for conductance as well as for the source current model as given in [3]. The model has been used by us to examine the properties of SFP in various neuromuscular disorders (Nandedkar et al. [5]).
The amplitude was measured from maximum positive to maximum negative peaks (Fig. 1). The zero-crossing following the initial downward and subsequent upward peaks were identified. The time difference between them is the duration of the negative peak ( Fig. 1)-henceforth it will be denoted by t z .
The change in the amplitude and duration with radial distance and muscle fiber diameter was investigated. In the experimental recordings, the peak duration and amplitude are measured while the muscle fiber diameter and radial distance are the ''unknown'' variables. Hence we may write: By fixing t z , we obtain parametric curve with coordinates: where the symbol j t z means that it is a parametric dependence of d and r on a for a fixed t z . Similarly by fixing a, we may obtain parametric dependence of r and d on t z : We have approximated the dependencies between d, r, a, and t z as follows: the negative peak duration was approximated by the following two bi-quadratic polynomials: where both F i¼1;2 are expressed as: where e i;j ¼ c i;j;1 þ xðc i;j;2 þ xc i;j;3 Þ ð 10Þ In our approach, the F i are quadratic functions of y which is either fiber diameter d for i ¼ 1, or fiber to electrode distance r for i ¼ 2. The coefficients of these quadratic polynomials are by themselves given by quadratic Eq. (10) with x equal to the log 10 ðaÞ. The coefficients c i;j;k are determined by fitting Eqs. (7) and (8) to results of simulation calculations.
We have found experimentally that the form of dependence given by Eqs. (7) and (8) is better suited to least squares bi-quadratic approximation than Eqs. (3) and (4).
In order to calculate the diameter, x ¼ log 10 ðaÞ has to be calculated from the amplitude of the potential and then used in Eq. (10) to calculate the coefficients e 1;1 . . .e 1;3 . Using these coefficients the quadratic Eq. (9) may be solved for y i.e., for the diameter. Since Eq. (9) is quadratic in y it may be solved directly. Similarly by calculating coefficients e 2;1 . . .e 2;3 Eq. (9) may be solved to obtain radius. The same Eq. (9) is used to derive both fiber diameter and radial distance. The quantity that is determined depends on whether one uses the set of coefficients c 1;j;k -for the determination of d, or the coefficients given in the second column in Table 1-for the determination of r from a given a and t z .

Results
The results of simulation studies are shown in Figs. 2, 3, 4, 5, and 6. The dependence of the weight function on parameters is shown in Fig. 2. The weight function is broader when the radial distance between the muscle fiber and the recording electrode is increased (Fig. 2a). For the same distance, the smaller fibers have a broader weight function due to slower conduction velocity (Fig. 2b). The simulated waveforms also show differences in the amplitude and durations. When their negative peaks are aligned and their amplitudes are normalized, larger fibers have a shorter negative peak duration (Fig. 3a). With identical radial distance, larger fibers have higher amplitude and shorter duration (Figs. 3b, 4). The amplitude changes much Fig. 1 Definitions of the SFP parameters. The amplitude was measured from the maximum positive (a) to maximum negative peak (b). The time difference between zero-crossings (c) and (d) is the duration of the negative peak t z Table 1 The values of the coefficients c 1;j;k and c 2;j;k used in Eqs. (9) and (10 more with radial distance than does the negative peak duration (Fig. 5). The amplitude and negative peak duration dependencies are shown on a single graph (Fig. 6). The solid curves represent the dependency of ðd; rÞj t z for a range of t z values (0.4-1.4 ms) and dashed curves represent the dependency of ðd; rÞj a for a range of amplitude values (50 lV-5 mV). A plot of these dependencies (Fig. 6) allows one to estimate the fiber diameter and radial distance knowing the amplitude and negative peak duration.
For example, if a potential with amplitude of 0.7 mV and a negative peak duration of 0.6 ms is recorded, we could determine the distance and diameter from the intersection of the curve for constant amplitude = 0.7 mV (solid line) and constant duration = 0.6 ms (dashed line). The point of intersection gives the fiber diameter (roughly 65 lm) and the radial distance (150 lm).
While in principle the graphical tool, provided that the curves for constant amplitude or negative peak duration are dense enough, may be used to determine fiber diameter and  Fig. 3 a SFP for fiber of diameter 30 and 90 lm located at a distance of 100 lm from the electrode. Both potentials have been scaled so that the maximum = 1. Potential due to larger fiber has been shifted in time (by ?5.9 ms) so that the maxima coincide. The increase of the width of potential with the decrease of diameter is clearly seen. b SFPs for fibers with diameter of 30, 45, 60, 75, and 90 lm located at a distance of 200 lm from electrode. It is seen how the amplitude of the SFP increases with the increase of fiber diameter and at the same time the potentials from larger fibers arrive earlier at the electrode Fig. 4 The negative peak duration (t z ) dotted line decreases with the increase in fiber diameter (at fixed fiber to electrode distance). The solid line is the dependence of amplitude on fiber diameter. With the increase of diameter the amplitude increases. The amplitude is measured in mV, the negative peak duration is measured in ms distance from needle we have approximated the data with analytical formulae which makes them more suitable for use. The graph, however, may be used to verify that the curves are smooth, nearly linear and that the two variables a and t z give rise to two families of curves that are well separated (Fig. 6).
The graph may be readily used to compare the diameters for two or more SFPs. If two SFPs have the same amplitude, so on the graph we are moving along one of the constant amplitude curves, then the one that has shorter negative peak duration is due to larger fiber. This is because the negative peak duration depends on conduction velocity which in turn is proportional to fiber diameter.
Hence the larger the fiber, the shorter the negative peak duration.
In order to explain the dependence of fiber diameter on amplitude for a constant negative peak duration one has to note that the amplitude depends primarily on two factorsthe distance from the electrode and fiber diameter (see Fig. 3b). For a constant t z , the amplitude changes mainly due to the change in fiber to electrode distance. The increase of radius, with the decreasing amplitude, leads to the broadening of the weight function. In order to keep the t z constant the fiber diameter has to increase with the increase of radius. Hence for two SFPs of the same negative peak duration the one for which the amplitude is smaller will be due to a larger fiber.
In order to calculate the values of the coefficients c 1;j;k and c 2;j;k we have simulated several hundreds of SFP with fiber diameter ranging from 25 to 90 lm and radial distance ranging from 50 to 500 lm calculating for each SFP its peak-to-peak amplitude and negative peak duration. These data were then approximated using formulae (9) and (10) and the coefficients have been determined by standard linear least squares method. These coefficients are given in Table 1. The mean error of the determination of t z from Eq. (7) is 5 Â 10 À3 , and for Eq. (8) it is 3 Â 10 À3 therefore we have rounded the coefficients to 4 significant places.
Using this method to approximate the data it is found that the root mean square error of the determination of diameter is 2 lm and for radial distance it is 6 lm. The maximum error for d is 8 lm, and for r it is 50 lm. The maximum errors occur at the ends of interpolation regions for r [ 310 lm, d\ 35 lm, or d [ 125 lm. Hence the approximations given by Eqs. (7) and (8) enable to determine d and r from a and t z with satisfactory accuracy.
For example, using these equations for a simulated SFP with fiber diameter d ¼ 55 lm, fiber to electrode distance r ¼ 80 lm for which a ¼ 1:357 mV, hence log 10 a ð Þ ¼ 0:13258 and t z ¼ 0:593 it is obtained that d ¼ 55:7 lm and r ¼ 80:9 lm. In this case the difference between the parameters used to model the SFP and their values obtained from calculation are less than 1 lm.
In order to estimate the effect of measurement errors of a and t z on the result we may assume a measurement error of AE20 lV for amplitude. For t z we assume an accuracy of AE0:04 ms which corresponds to sampling frequency of 25 kHz. It is found that with the assumed value for the error in a the change in the derived diameter is less than 0.3 lm. By changing t z by the assumed value the d changes by 8 lm. It can be easily verified using the coefficients of the fit that the diameter is more sensitive to measurement errors in t z than to errors in a. Therefore, to be able to determine fiber diameter and fiber to electrode distance a good quality of the SFP is required with sufficient Fig. 5 The dependence of t z (dotted line) and amplitude (solid line) on the electrode distance from the fiber for a fiber of diameter of 55 lm and electrode located at distances from 50 to 300 lm. The amplitude is measured in mV, the negative peak duration is measured in ms. The amplitude decreases with distance and the negative peak duration increases t z (ms) a (mV) Fig. 6 The dependence of fiber diameter (d) on fiber to electrode distance (r) for a fixed amplitude (a) of the potential or a fixed duration ðt z Þ. Solid curves represent constant duration and dashed curves represent constant amplitude. The values of amplitude are given along the top and right margins of the graph. The values of constant duration are given in the lower part of the graph resolution in time. It is also for this reason that the graph tool can be only used for approximate or qualitative estimates as it would have to contain tens of constant t z curves to provide comparable accuracy to the presented formulae.

Discussion
From the simulation studies it has been found that the negative peak duration of SFP is a quantity which is useful in the determination of fiber properties and in particular of fiber diameter. Rodriguez and co-workers [6,7] have shown that this quantity is related to the width of the source current generating the potential. We have found that t z is related to fiber diameter. This relation is easily understandable in view of the above findings because the duration is proportional to the ratio of the width (*1.5 mm) of the negative peak of the second derivative of intracellular potential [3] and source propagation velocity. The velocity of propagation is dependent on fiber diameter; therefore the larger the diameter, the less time it will take for the feature to pass the electrode, hence the negative peak duration will be shorter. Thus, it turns that t z is one of the basic quantities characterizing the SFP and it gives a direct link to the underlying physiology.
We have derived a set of analytical formulae using which it is possible to determine fiber diameter (and fiber to electrode distance) from the measurement of the SFPs amplitude and negative peak duration. The analysis of the sensitivity of the determination of fiber diameter to the errors in the measurement of amplitude and negative peak duration shows that a sampling rate of 25 kHz is required so that the error in the determined fiber diameter is less than 8 lm. The errors in the determination of amplitude have much smaller impact on the error of the determination of d. In order to be able to determine the fiber diameter and fiber to electrode distance, the noise level in the recorded SFP has to be low (\20 lV).
The ability to quantify the fiber diameter can add significantly to the SFEMG examination. It may help one to study muscle fiber diameter variability in many different muscles and thus compliment the muscle biopsy studies.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.