Branched EMG electrodes for stable and selective recording of single motor unit potentials in humans.

Branched surface EMG electrodes are bipolar electrodes with the hot signal pole referenced to two or more short-circuited leading-off surfaces. This technique provides stable recording of single motor unit potentials during real movements, up to maximal muscle contractions. The selective characteristic of branched electrodes is based on the same principles as the double differential detection system and spatial filtering technique proposed later. Equi-weight calculations to assess the selectivity of different electrode types and their position are used. The main advantage of branched electrodes, especially high stability, is achieved by the wire electrode version. The design, manufacture, implementation, and application of wire electrodes are discussed in detail. During recording of motor unit potentials, electrodes are positioned subcutaneously over the muscle fascia. This positioning maximizes electrode stability. Appropriate orientation of the electrode relative to the muscle architecture ensures adequate selectivity for single motor unit recordings. Branched electrodes require ordinary EMG equipment (two or even one amplifier).     

In vivo magnetic and electric recordings from nerve bundles and single motor units in mammalian skeletal muscle. Correlations with muscle force

Recent advances in the technology of recording magnetic fields associated with electric current flow in biological tissues have provided a means of examining action currents that is more direct and possibly more accurate than conventional electrical recording. Magnetic recordings are relatively insensitive to muscle movement, and, because the recording probes are not directly connected to the tissue, distortions of the data due to changes in the electrochemical interface between the probes and the tissue are eliminated. In vivo magnetic recordings of action currents of rat common peroneal nerve and extensor digitorum longus (EDL) muscle were obtained by a new magnetic probe and amplifier system that operates within the physiological temperature range. The magnetically recorded waveforms were compared with those obtained simultaneously by conventional, extracellular recording techniques. We used the amplitude of EDL twitch force (an index of stimulus strength) generated in response to graded stimulation of the common peroneal nerve to enable us to compare the amplitudes of magnetically recorded nerve and muscle compound action currents (NCACs and MCACs, respectively) with the amplitudes of electrically recorded nerve compound action potentials (NCAPs). High, positive correlations to stimulus strength were found for NCACs (r = 0.998), MCACs (r = 0.974), and NCAPs (r = 0.998). We also computed the correlations of EDL single motor unit twitch force with magnetically recorded single motor unit compound action currents (SMUCACs) and electrically recorded single motor unit compound action potentials (SMUCAPs) obtained with both a ring electrode and a straight wire serving as a point electrode. Only the SMUCACs had a relatively strong positive correlation (r = 0.768) with EDL twitch force. Correlations for ring and wire electrode-recorded SMUCAPs were 0.565 and -0.366, respectively. This study adds a relatively direct examination of action currents to the characterization of the normal biophysical properties of peripheral nerve, muscle, and muscle single motor units.     

Motor unit acceleration maps and interference mechanomyographic distribution

The study analyses the two-dimensional distribution of surface mechanomyographic (MMG) signal generated by the activation of single motor units located in three transverse positions in the tibialis anterior muscle. In 12 healthy volunteers, surface MMG signals were recorded from the tibialis anterior muscle with a 3x4 grid of accelerometers spaced by 20 and 30mm in the transverse and longitudinal direction. Three intramuscular electromyographic (EMG) signals were recorded with wire electrodes inserted 20-mm apart, between the first and second most proximal accelerometers of each column of the grid. The subject was asked to activate three different motor units (target motor units) in three contractions with visual feedback from each of the three intramuscular recordings (three locations). The MMG signals from the 12 accelerometers were averaged using the intramuscular single motor unit action potentials as trigger in order to obtain surface motor unit acceleration maps (MUAMs). The peak-to-peak value of the averaged MMG depended on motor unit location (P<0.001) and on the transverse position of the accelerometer in the grid (P<0.05). Moreover, MUAM amplitude depended on the interaction between motor unit location and transverse accelerometer position (P<0.05), demonstrating an influence of motor unit location on the generated MUAM. The observed dependency of MUAMs on motor unit location provides a quantitative analysis of the effect of the volume conductor on the recorded surface MMG signal.     

Spatial analysis of motor unit potentials of the rat medial gastrocnemius

The spatial analysis of the potentials of single motor units of the rat medial gastrocnemius muscle evoked by stimulation of the fibres of split ventral roots was carried out with a bipolar electrode moving in the direction perpendicular to the longitudinal axis of the muscle fibres. During this movement of the electrode a variability was observed in the time of the biphasic potential from its maximum to minimum, and in the peak-to-peak amplitude of these potentials. The potentials recorded outside the territory of the motor unit had a lower amplitude in relation to the potentials from the territory of the unit. This made localization of the motor unit on the cross-section of the muscle possible. Differences in the duration of the potential from maximal to minimal amplitude (maximum-minimum amplitude time--M-MAT) of each investigated motor unit from successive recording sites reflected the number of fibres contributing to the action potential and the distance of the recording surface of the electrode from the zone of the motor end-plates of this motor unit. The greatest diameter of the territory of the observed motor units reached 2.5 mm.     

Detection of motor unit action potentials with surface electrodes: influence of electrode size and spacing

A model of the motor unit action potential was developed to investigate the amplitude and frequency spectrum contributions of motor units, located at various depths within muscle, to the surface detected electromyographic (EMG) signal. A dipole representation of the transmembrane current in a three-dimensional muscle volume was used to estimate detected individual muscle fiber action potentials. The effects of anisotropic muscle conductance, innervation zone location, propagation velocity, fiber length, electrode area, and electrode configuration were included in the fiber action potential model. A motor unit action potential was assumed to be the sum of the individual muscle fiber action potentials. A computational procedure, based on the notion of isopotential layers, was developed which substantially reduced the calculation time required to estimate motor unit action potentials. The simulations indicated that: 1) only those motor units with muscle fibers located within 10–12 mm of the electrodes would contribute significant signal energy to the surface EMG, 2) variation in surface area of electrodes has little effect on the detection depth of motor unit action potentials, 3) increased interelectrode spacing moderately increases detection depth, and 4) the frequency content of action potentials decreases steeply with increased electrode-motor unit territory distance.     

Electromyographic measures of muscle activation and changes in muscle architecture of human elbow flexors during fatiguing contractions.

The study compared changes in intramuscular and surface recordings of EMG amplitude with ultrasound measures of muscle architecture of the elbow flexors during a submaximal isometric contraction. Ten subjects performed a fatiguing contraction to task failure at 20% of maximal voluntary contraction force. EMG activity was recorded in biceps brachii, brachialis, and brachioradialis muscles using intramuscular and surface electrodes. The rates of increase in the amplitude of the surface EMG for the long and short heads of biceps brachii and brachioradialis were greater than those for the intramuscular recordings measured at different depths. The amplitude of the intramuscular recordings from three muscles increased at a similar rate (P = 0.13), as did the amplitude of the three surface recordings from two muscles (P = 0.83). The increases in brachialis thickness (27.7 +/- 5.7 to 30.9 +/- 3.5 mm; P < 0.05) and pennation angle (10.9 +/- 3.5 to 16.5 +/- 4.8 degrees ; P = 0.003) were not associated with the increase in intramuscular EMG amplitude (P > 0.58). The increase in brachioradialis thickness (22.8 +/- 4.8 to 25.5 +/- 3.4 mm; P = 0.0075) was associated with the increase in the amplitude for one of two intramuscular EMG signals (P = 0.007, r = 0.79). The time to failure was more strongly associated with the rate of increase in the amplitude of the surface EMG than that for the intramuscular EMG, which suggests that the surface measurement provides a more appropriate measure of the change in muscle activation during a fatiguing contraction.     

Comparison of force and EMG measures in normal and reinnervated tibialis anterior muscles of the rat.

The relationship between motor unit force and the recorded voltage produced by activated muscle unit fibres (electromyogram, EMG) was examined in normal and reinnervated rat tibialis anterior muscles. The number, cross-sectional area, and radial distance from the recording electrode of muscle fibres in a given unit, obtained directly from a sample of glycogen-depleted motor units, were analysed in relation to the magnitude of the EMG signal produced by that unit. EMG peak to peak amplitude and area varied as approximately the square root of twitch force in both normal and reinnervated units. Furthermore, the EMG amplitude increased approximately as the total cross-sectional area of the motor unit (number of muscle fibres x the average cross-sectional area of the fibres) and inversely with approximately the square root of the distance of fibres from the recording electrodes on the surface of the muscle.     

Influence of amplitude cancellation on the simulated surface electromyogram.

The purpose of the study was to quantify the influence of selected motor unit properties and patterns of activity on amplitude cancellation in the simulated surface electromyogram (EMG). The study involved computer simulations of a motor unit population with physiologically defined recruitment and rate coding characteristics that activated muscle fibers whose potentials were recorded on the skin over the muscle. Amplitude cancellation was quantified as the percent difference in signal amplitude when motor unit potentials were summed before and after rectification. The simulations involved varying the level of activation for the motor unit population, the recording configuration, the upper limit of motor unit recruitment, peak discharge rates, the amount of motor unit synchronization, muscle fiber length, the thickness of the subcutaneous tissue, and the motor unit properties that change with advancing age. The results confirmed a previous experimental report (Day SJ and Hulliger M, J Neurophysiol 86: 2144-2158, 2001) that amplitude cancellation in the surface EMG can reach 62% at maximal activation. A decrease in the range of amplitudes of the motor unit potentials, as can occur during fatiguing contractions, increased amplitude cancellation up to approximately 85%. Differences in the amount of amplitude cancellation were observed across all simulated conditions, and resulted in substantial changes in the absolute magnitude of the EMG signal. The most profound factors influencing amplitude cancellation were the number of active motor units and the duration of the action potentials. The effects of amplitude cancellation were minimal (<5%) when the EMG amplitude was normalized to maximal values, with the exception of variations in peak discharge rate and recruitment range, which resulted in differences up to 17% in the normalized EMG signal across conditions. These results indicate the amount of amplitude cancellation that can occur in various experimental conditions and its influence on absolute and relative measures of EMG amplitude.     

Measuring surface potential components necessary for transmembrane current computation using microfabricated arrays

This study was designed to test the feasibility of using microfabricated electrodes to record surface potentials with sufficiently fine spatial resolution to measure the potential gradients necessary for improved computation of transmembrane current density. To assess that feasibility, we recorded unipolar electrograms from perfused rabbit right ventricular free wall epicardium (n = 6) using electrode arrays that included 25-microm sensors fabricated onto a flexible substrate with 75-microm interelectrode spacing. Electrode spacing was therefore on the size scale of an individual myocyte. Signal conditioning adjacent to the sensors to control lead noise was achieved by routing traces from the electrodes to the back side of the substrate where buffer amplifiers were located. For comparison, recordings were also made using arrays built from chloridized silver wire electrodes of either 50-microm (fine wire) or 250-microm (coarse wire) diameters. Electrode separations were necessarily wider than with microfabricated arrays. Comparable signal-to-noise ratios (SNRs) of 21.2 +/- 2.2, 32.5 +/- 4.1, and 22.9 +/- 0.7 for electrograms recorded using microfabricated sensors (n = 78), fine wires (n = 78), and coarse wires (n = 78), respectively, were found. High SNRs were maintained in bipolar electrograms assembled using spatial combinations of the unipolar electrograms necessary for the potential gradient measurements and in second-difference electrograms assembled using spatial combinations of the bipolar electrograms necessary for surface Laplacian (SL) measurements. Simulations incorporating a bidomain representation of tissue structure and a two-dimensional network of guinea pig myocytes prescribed following the Luo and Rudy dynamic membrane equations were completed using 12.5-microm spatial resolution to assess contributions of electrode spacing to the potential gradient and SL measurements. In those simulations, increases in electrode separation from 12.5 to 75.0, 237.5, and 875.0 microm, which were separations comparable to the finest available with our microfabricated, fine wire, and coarse wire arrays, led to 10%, 42%, and 81% reductions in maximum potential gradients and 33%, 76%, and 96% reductions in peak-to-peak SLs. Maintenance of comparable SNRs for source electrograms was therefore important because microfabrication provides a highly attractive methods to achieve spatial resolutions necessary for improved computation of transmembrane current density.     

Analysis of motor units with high-density surface electromyography

Although the behaviour of individual motor units is classically studied with intramuscular EMG, recently developed techniques allow its analysis also from EMG recorded in multiple locations over the skin surface (high-density surface EMG). The analysis of motor units from the surface EMG is useful when the insertion of needles is not desirable or not possible. Moreover, surface EMG allows the measure of motor unit properties which are difficult to assess with invasive technology (e.g., muscle fiber conduction velocity or location of innervation zones) and may increase the number of detectable motor units with respect to selective intramuscular recordings. Although some limitations remain, both the discharge pattern and muscle fiber properties of individual motor units can currently be analyzed non-invasively. This review presents the conditions and methodologies which allow the investigation of motor units with surface EMG.     

Neuromuscular adjustments that constrain submaximal EMG amplitude at task failure of sustained isometric contractions

The amplitude of the surface EMG does not reach the level achieved during a maximal voluntary contraction force at the end of a sustained, submaximal contraction, despite near-maximal levels of voluntary effort. The depression of EMG amplitude may be explained by several neural and muscular adjustments during fatiguing contractions, including decreased net neural drive to the muscle, changes in the shape of the motor unit action potentials, and EMG amplitude cancellation. The changes in these parameters for the entire motor unit pool, however, cannot be measured experimentally. The present study used a computational model to simulate the adjustments during sustained isometric contractions and thereby determine the relative importance of these factors in explaining the submaximal levels of EMG amplitude at task failure. The simulation results indicated that the amount of amplitude cancellation in the simulated EMG (～ 40%) exhibited a negligible change during the fatiguing contractions. Instead, the main determinant of the submaximal EMG amplitude at task failure was a decrease in muscle activation (number of muscle fiber action potentials), due to a reduction in the net synaptic input to motor neurons, with a lesser contribution from changes in the shape of the motor unit action potentials. Despite the association between the submaximal EMG amplitude and reduced muscle activation, the deficit in EMG amplitude at task failure was not consistently associated with the decrease in neural drive (number of motor unit action potentials) to the muscle. This indicates that the EMG amplitude cannot be used as an index of neural drive.     

Specificity of esophageal electrode recordings of posterior cricoarytenoid muscle activity

Esophageal electrodes have been used for recording the electromyographic (EMG) activity of the posterior cricoarytenoid muscle (PCA). To determine the specificity of this EMG technique, esophageal electrode recordings were compared with intramuscular recordings in eight anesthetized mongrel dogs. Intramuscular wire electrodes were placed in the right and left PCA, and the esophageal electrode was introduced through the nose or mouth and advanced into the upper esophagus. On direct visualization of the upper airway, the unshielded catheter electrode entered the esophagus on the right or left side. Cold block of the recurrent laryngeal nerve (RLN) ipsilateral to the esophageal electrode was associated with a marked decrease in recorded activity, whereas cold block of the contralateral RLN resulted only in a small reduction in activity. After supplemental doses of anesthesia were administered, bilateral RLN cold block essentially abolished the activity recorded with the intramuscular electrodes as well as that recorded with the esophageal electrode. Before supplemental doses of anesthesia were given, especially after vagotomy, the esophageal electrode, and in some cases the intramuscular electrodes, recorded phasic inspiratory activity not originating from the PCA. Therefore, one should be cautious in interpreting the activity recorded from esophageal electrodes as originating from the PCA, especially in conditions associated with increased respiratory efforts.     

Using two-dimensional spatial information in decomposition of surface EMG signals.

Recently, high-density surface EMG electrode grids and multi-channel amplifiers became available for non-invasive recording of human motor units (MUs). We present a way to decompose surface EMG signals into MU firing patterns, whereby we concentrate on the importance of two-dimensional spatial differences between the MU action potentials (MUAPs). Our method is exemplified with high-density EMG data from the vastus lateralis muscle of a single subject. Bipolar and Laplacian spatial filtering was applied to the monopolar raw signals. From the single recording in this subject six different simultaneously active MUs could be distinguished using the spatial differences between MUAPs in the direction perpendicular to the muscle fiber direction. After spike-triggered averaging, 125-channel two-dimensional MUAP templates were obtained. Template-matching allowed tracking of all MU firings. The impact of spatial information was measured by using subsets of the MUAP templates, either in parallel or perpendicular to the muscle fiber direction. The use of one-dimensional spatial information perpendicular to the muscle fiber direction was superior to the use of a linear array electrode in the longitudinal direction. However, to detect the firing events of the MUs with a high accuracy, as needed for instance for estimation of firing synchrony, two-dimensional information from the complete grid electrode appears essential.     

Electromyographic instantaneous amplitude and instantaneous mean power frequency patterns across a range of motion during a concentric isokinetic muscle action of the biceps brachii.

The purpose of this study was to examine the electromyographic (EMG) instantaneous amplitude (IA) and instantaneous mean power frequency (IMPF) patterns for the biceps brachii muscle across a range of motion during maximal and submaximal concentric isokinetic muscle actions of the forearm flexors. Ten adults (mean +/- SD age = 22.0 +/- 3.4 years) performed a maximal and a submaximal [20% peak torque (PT)] concentric isokinetic forearm flexion muscle action at a velocity of 30 degrees s(-1). The surface EMG signal was detected from the biceps brachii muscle with a bipolar electrode arrangement, and the EMG IA and IMPF versus time relationships were examined for each subject using first- and second-order polynomial regression models. The results indicated that there were no consistent patterns between subjects for EMG IA or IMPF with increases in torque across the range of motion. Some of the potential nonphysiological factors that could influence the amplitude and/or frequency contents of the surface EMG signal during a dynamic muscle action include movement of the muscle fibers and innervation zone beneath the skin surface, as well as changes in muscle fiber length and the thickness of the tissue layer between the muscle and the recording electrodes. These factors may affect the EMG IA and IMPF patterns differently for each subject, thereby increasing the difficulty of drawing any general conclusions regarding the motor control strategies that increase torque across a range of motion.     

Spatial distribution of surface action potentials generated by individual motor units in the human biceps brachii muscle

This study analyses the spatial distribution of individual motor unit potentials (MUPs) over the skin surface and the influence of motor unit depth and recording configuration on this distribution. Multichannel surface (13 × 5 electrode grid) and intramuscular (wire electrodes inserted with needles of lengths 15 and 25 mm) electromyographic (EMG) signals were concurrently recorded with monopolar derivations from the biceps brachii muscle of 10 healthy subjects during 60-s isometric contractions at 20% of the maximum torque. Multichannel monopolar MUPs of the target motor unit were obtained by spike-triggered averaging of the surface EMG. Amplitude and frequency characteristics of monopolar and bipolar MUPs were calculated for locations along the fibers’ direction (longitudinal), and along the direction perpendicular (transverse) to the fibers. In the longitudinal direction, monopolar and bipolar MUPs exhibited marked amplitude changes that extended for 16–32 mm and 16–24 mm over the innervation and tendon zones, respectively. The variation of monopolar and bipolar MUP characteristics was not symmetrical about the innervation zone. Motor unit depth had a considerable influence on the relative longitudinal variation of amplitude for monopolar MUPs, but not for bipolar MUPs. The transverse extension of bipolar MUPs ranged between 24 and 32 mm, whereas that of monopolar MUPs ranged between 72 and 96 mm. The mean power spectral frequency of surface MUPs was highly dependent on the transverse electrode location but not on depth. This study provides a basis for the interpretation of the contribution of individual motor units to the interference surface EMG signal.     

Phonomyogram and electromyogram relationships with isometric force reinvestigated in man

Contracting muscle generates sounds which can be recorded easily by means of a microphone. To determine if a phomomyogram (PMG) can be used to monitor muscle force, a comparison was made between simultaneous recordings of PMG and monopolar electromyogram (EMG) from the isometrically contracting biceps brachii muscle and the external flexion force. Locations of the monopolar electrode and microphone were identified in relation to the motor point. Whatever the recording site, PMG amplitude was proportional to EMG amplitude and both showed a quadratic relationship to muscle force. Changes in the PMG spectrum with force were similar to those in EMG, i.e. the mean power frequency increased up to about 30% maximal voluntary contraction and then reached a plateau. Despite a slightly higher variability, PMG was shown to be a valid index of muscular isometric force. At the same force, the amplitude of both PMG and EMG was lower in the prone than in the supine position of the hand. This result indicated a selective recording of biceps brachii muscle activity.     

Inter-operator agreement in decomposition of motor unit firings from high-density surface EMG.

High-density surface EMG can be used to obtain a spatially selective representation of several motor unit action potentials. Recently, a decomposition of the signal into the underlying motor neuron firing patterns has been described. The reliability of the algorithm has not yet been tested. Eleven healthy subjects participated. High-density surface EMG was recorded from the vastus lateralis muscle during an isometric knee extension. Two independent operators analyzed the signals. After operator-supervised cluster analysis of spikes, motor unit action potential templates were constructed and an automatic template matching was performed. The decomposition was adjusted by hand. Agreement between operators was calculated for the number of coincident firings. Bland-Altman plots of peak-to-peak amplitude were constructed and limits of agreement were calculated. For completely decomposed motor unit action potential trains the between-operator agreement of firing events was very high. The peak-to-peak amplitude of monopolar motor unit action potentials was 115microV (SD 74microV). The agreement was within 3microV and independent of amplitude. With partial decomposition agreement within 26microV was achieved. For bipolarly derived motor unit action potentials the peak-to-peak amplitude was 54microV (SD 49microV), the agreement was within 3microV. Only for recordings obtained from a force level below 5% of the maximum voluntary contraction full decomposition was possible. It was concluded that when full decomposition is achieved, two independent operators are likely to arrive at nearly identical firing patterns.     

Consensus for experimental design in electromyography (CEDE) project: Single motor unit matrix

The analysis of single motor unit (SMU) activity provides the foundation from which information about the neural strategies underlying the control of muscle force can be identified, due to the one-to-one association between the action potentials generated by an alpha motor neuron and those received by the innervated muscle fibers. Such a powerful assessment has been conventionally performed with invasive electrodes (i.e., intramuscular electromyography (EMG)), however, recent advances in signal processing techniques have enabled the identification of single motor unit (SMU) activity in high-density surface electromyography (HDsEMG) recordings. This matrix, developed by the Consensus for Experimental Design in Electromyography (CEDE) project, provides recommendations for the recording and analysis of SMU activity with both invasive (needle and fine-wire EMG) and non-invasive (HDsEMG) SMU identification methods, summarizing their advantages and disadvantages when used during different testing conditions. Recommendations for the analysis and reporting of discharge rate and peripheral (i.e., muscle fiber conduction velocity) SMU properties are also provided. The results of the Delphi process to reach consensus are contained in an appendix. This matrix is intended to help researchers to collect, report, and interpret SMU data in the context of both research and clinical applications.     

Amplitude cancellation reduces the size of motor unit potentials averaged from the surface EMG.

The purpose of the study was to evaluate the influence of selected physiological parameters on amplitude cancellation in the simulated surface electromyogram (EMG) and the consequences for spike-triggered averages of motor unit potentials derived from the interference and rectified EMG signals. The surface EMG was simulated from prescribed recruitment and rate coding characteristics of a motor unit population. The potentials of the motor units were detected on the skin over a hand muscle with a bipolar electrode configuration. Averages derived from the EMG signal were generated using the discharge times for each of the 24 motor units with lowest recruitment thresholds from a population of 120 across three conditions: 1) excitation level; 2) motor unit conduction velocity; and 3) motor unit synchronization. The area of the surface-detected potential was compared with potentials averaged from the interference, rectified, and no-cancellation EMGs. The no-cancellation EMG comprised motor unit potentials that were rectified before they were summed, thereby preventing cancellation between the opposite phases of the potentials. The percent decrease in area of potentials extracted from the rectified EMG was linearly related to the amount of amplitude cancellation in the interference EMG signal, with the amount of cancellation influenced by variation in excitation level and motor unit conduction velocity. Motor unit synchronization increased potentials derived from both the rectified and interference EMG signals, although cancellation limited the increase in area for both potentials. These findings document the influence of amplitude cancellation on motor unit potentials averaged from the surface EMG and the consequences for using the procedure to characterize motor unit properties.     

Electromyogram recordings from freely moving animals

Electromyography can be used to record activity from sets of muscles in awake, freely moving animals using implanted intramuscular electrodes. As a tool, EMG has a wide range of applications ranging from inferring neural processes to analyzing movement. The amplitude of the rectified and filtered electromyogram (EMG) can be used as an indirect measure of muscle activity. Although it is often tempting to correlate the EMG with muscle force, the fact that force varies more with different activation strategies than with EMG estimates must be taken into account. The purpose of this article is to provide the researcher wishing to introduce the technique of recording EMGs from conscious animals using intramuscular electrodes with a step-by-step guide. It includes details on the manufacture of electromyograph electrodes, recording, and analysis considerations along with a section on solving common problems. For the sake of clarity, this article focuses on using the cat as a model and on the implantation of hindlimb muscles with intramuscular wire electrodes. However, the procedures can be adapted for use on other striated muscles and species.