4.1 Muscle excitation
4.1.1 Visualization
To visualize sEMG patterns across channels, the median RMS and MNF muscle excitation across limbs and hand grasp configurations were calculated as participants attempted to perform each hand grasp (Fig. 4, Fig. 5, and Supplementary Figs S1-S18 online). For the RMS characteristic, the unaffected limb behaved as expected with only a subset of sEMG channels recording muscle excitation. For example, during wrist extension and flexion, the extensor carpi radialis longus and flexor carpi radialis muscles are activated separately for each wrist movement43. A similar phenomenon can be seen in Fig. 4 as exemplified by participant SHR-C and throughout participant data sets. Different patterns of sEMG activity were observed for the RMS characteristic across attempted hand grasp configurations of participants' affected limbs, thereby providing a visual indication of the extent to which children with UCBED could actuate their affected muscles. Interestingly, although not in all cases, when the participants were asked to attempt a hand grasp with their affected limb the majority of sEMG channels recorded muscle excitation above the relaxed state (see Fig. 4 and the RMS Supplementary Figures online). Additionally, it is important to note that for some participants the affected limb RMS characteristic across the sEMG channels demonstrated relatively large degrees of variability, e.g., participant SHR-B produced a large RMS spread across sEMG channels for most hand grasp configurations (see Supplementary Fig. S3 online). In contrast, participant SHR-D had very few visually distinct patterns of RMS muscle excitation with similar spread across grasps. Therefore, participants have shown (1) visually distinct patterns of muscle excitation across hand grasp configurations, (2) minimal reproducibility with large spreads, and (3) consistent but few patterns of muscle excitation, all of which play a role in the ability to decode motor intent.
Insert Fig. 4
For the MNF characteristic, patterns were illustrated in frequency-based muscle excitation for the unaffected and affected limbs among participants across sEMG channels. The MNF muscle excitation for both limbs showed parallels to their RMS counterparts. While the unaffected excitation differed across sEMG channels from grasp to grasp, the majority of the affected excitation was above the relaxed state, as illustrated by participant SHR-C Fig. 5 and the MNF Supplementary Figures online. Although there were variations in the degree of visual patterns for the RMS characteristics across grasps, the MNF characteristics revealed additional distinct patterns for most participants across both limbs (see MNF Supplementary Figures online). This indicates that separately and together RMS and MNF characteristics may provide sufficient information to differentiate between hand grasp configurations.
Insert Fig. 5
4.1.2 Grasp excitation consistency
Within-grasp consistency was calculated from Kendall’s Coefficient of Concordance, W, to quantify participants' ability to consistently perform patterns of RMS and MNF muscle excitation for each limb over the multiple trials of hand grasp configurations. Here consistencies ranged from W < 0.20, 0.20 ≤ W < 0.40, 0.40 ≤ W < 0.60, 0.60 ≤ W < 0.80, and W ≥ 0.80 with poor, minimal, weak, moderate, and strong consistency, respectively. In general, for the unaffected limb, measures of consistency were moderate to strong which is indicative of typical repetitive muscle actuation and few participants had poor to weak consistency. Specifically, for the RMS consistency, all participants except for SHR-B had at least 9 hand grasp configurations with moderate to strong consistency i.e., a W of greater than 0.60. Interestingly, SHR-B had the top two lowest RMS consistencies in the unaffected limb for the grasp-type cylindrical wrap with wrist rotated (CR) and cylindrical wrap (CW), with values of W = 0.25 and W = 0.27, respectively. SHR-B only had 4 out of 10 hand grasp configurations with moderate to strong consistencies. Similarly, participant SHR-I had the third lowest consistency of W = 0.30 for wrist flexion (WF). Alternatively, for the MNF unaffected limb within-grasp consistency, all participants except SHR-B had moderate to strong consistency values for at least 5 hand grasp configurations while SHR-B had only 2. The top three lowest consistencies were present in participants SHR-F (W = 0.16 for wrist extension (WE)), SHR-A (W = 0.21 for pulp pinch (PP)), and SHR-B (W = 0.24 for wrist rotation (WR)) in the poor to minimal consistency range. All hand grasp consistency values for the unaffected limb can be seen in the first column of Fig. 6.
When attempting to move their missing limb, participants were able to produce repeatable patterns of affected muscle excitation, although there was variability across participants. Repeatable patterns are an important facet for robust control of prostheses given that 6–9 common hand grasp movements can account for nearly 80% of activities in daily living25,44, indicating the efficacy of prosthesis use in this population. When evaluating the RMS within-grasp consistency in the affected limb, we found at least 5 of the 10 hand grasp configurations had moderate to strong consistency. Here, participant SHR-A had the lowest consistency value of W = 0.14 for the cylindrical wrap (CW). Participant SHR-B had the second and third lowest consistency values: W = 0.17 in the cylindrical wrap (CW) and W = 0.22 in the cylindrical wrap with wrist rotated (CR). Additionally, the MNF within-grasp consistency for the affected limb showed that all participants except SHR-B had moderate to strong consistency values for 4 or more hand grasp configurations. SHR-G had the lowest consistency value of W = 0.13 for index flexion (IF). Participant SHR-B had only 1 hand grasp, wrist flexion (WF), with a moderate consistency value, and had the second through the sixth lowest consistency values from poor to weak. All other consistency values for the affected limb can be seen in the second column of Fig. 6. Together, we see that 5 of 10 grasps for RMS and 4 of 10 grasps for MNF approached moderate to strong consistency (i.e., nearing the consistency for 6–9 hand grasp movements). In conclusion, this indicates that the participants were able to perform reproducible attempted hand grasp configurations in their affected limb to a similar degree as their unaffected limb.
Insert Fig. 6
4.2 Across grasp dissimilarity
To further understand if these hand grasps were dissimilar or distinguishable from one another, additional evaluation was performed given the emergence of muscle excitation patterns in previous visualizations across sEMG channels when participants attempted each hand grasp configuration. To determine whether sufficient information for distinguishable structure in hand grasp configurations was present, we performed the sdRDM analysis. As such, the analysis demonstrated whether pairwise distances between grasps for RMS and MNF measures of spread (i.e., interquartile range (IQR)) and amplitude were typical of muscle excitation movements. To note, typical muscle excitation has been observed to be the variability in sEMG signals across channels from task to task, or even within the same movement15. Furthermore, to qualitatively visualize the distinguishability across hand grasp configurations, multidimensional scaling (MDS) was used.
The sdRDM analysis of the median and IQR measures for RMS and MNF indicated that distinguishable structures of sEMG data were present when the participants with UCBED attempted to perform missing hand grasp configurations. This analysis also showed distinguishable structures of sEMG data were present in participants’ unaffected limbs, as is expected of typical muscle contractions. However, this was not the case for all participants e.g., the analysis did not provide sufficient evidence to suggest there was a distinguishable structure for the median RMS of SHR-B’s unaffected and affected limbs with p = 0.188 and p = 0.267, respectively. Additionally, sdRDM analysis of the affected limb did not provide sufficient evidence to suggest a distinguishable structure for the following participants: SHR-D (p = 0.106, RMS IQR), SHR-F (p = 0.118, MNF IQR), and SHR-I (p = 0.163, RMS IQR). All other sdRDM analyses for the affected and unaffected limb showed a distinguishable structure of the hand grasp configurations; all p-values for each limb across participants and characteristics can be found in Table 2.
Table 2
Split-data representational dissimilarity analysis to distinguish the structure of amplitude and spread of RMS and MNF characteristics.
Participants | Unaffected | Affected |
RMS EDI | MNF EDI | RMS EDI | MNF EDI |
Median | IQR | Median | IQR | Median | IQR | Median | IQR |
SHR-A | 0.552 p < 0.001 | 0.483 p < 0.001 | 0.629 p < 0.001 | 0.235 p = 0.005 | 0.319 p = 0.003 | 0.459 p = 0.006 | 0.519 p < 0.001 | 0.452 p < 0.001 |
SHR-B | 0.095‡ | 0.546 p < 0.001 | 0.252‡ | 0.728† | 0.047‡ | 0.256 p = 0.036 | 0.167‡ | 0.510† |
p = 0.188* | p = 0.018 | p < 0.001 | p = 0.267* | p = 0.023 | p < 0.001 |
SHR-C | 0.550 p < 0.001 | 0.472 p < 0.001 | 0.720 p < 0.001 | 0.377 p < 0.001 | 0.535 p < 0.001 | 0.294 p = 0.003 | 0.571 p < 0.001 | 0.423 p < 0.001 |
SHR-D | 0.455 p < 0.001 | 0.463 p < 0.001 | 0.667 p < 0.001 | 0.205‡ | 0.466 p < 0.001 | 0.114‡ | 0.665† | 0.199 p = 0.015 |
p = 0.020 | p = 0.106* | p < 0.001 |
SHR-E | 0.421 p < 0.001 | 0.453 p < 0.001 | 0.648 p < 0.001 | 0.235 p = 0.007 | 0.489 p < 0.001 | 0.273 p < 0.001 | 0.599 p < 0.001 | 0.212 p = 0.018 |
SHR-F | 0.554 p < 0.001 | 0.459 p < 0.001 | 0.438 p < 0.001 | 0.603 p < 0.001 | 0.474 p < 0.001 | 0.254 p = 0.004 | 0.590 p < 0.001 | 0.119‡ |
p = 0.118* |
SHR-G | 0.635† | 0.578† | 0.783† | 0.449 p < 0.001 | 0.629† | 0.519† | 0.529 p < 0.001 | 0.303 p = 0.007 |
p < 0.001 | p < 0.001 | p < 0.001 | p < 0.001 | p = 0.002 |
SHR-H | 0.457 p < 0.001 | 0.368‡ | 0.468 p < 0.001 | 0.214 p = 0.007 | 0.351 p < 0.001 | 0.212 p = 0.008 | 0.233 p = 0.005 | 0.364 p < 0.001 |
p < 0.001 |
SHR-I | 0.552 p < 0.001 | 0.497 p < 0.001 | 0.453 p < 0.001 | 0.366 p < 0.001 | 0.474 p < 0.001 | 0.116 | 0.309 p < 0.001 | 0.180 p = 0.021 |
p = 0.163* |
*Red p-values from the sdRDM analysis indicate a failure to reject the null hypothesis (i.e., there was not sufficient evidence to suggest a distinguishable structure of the hand grasp configurations), given a significance level of α = 0.05. |
†Green EDI values indicate the maximum within a column. |
‡Yellow EDI values indicate the minimum within a column. |
Further investigation of the sdRDM exemplar discriminability index (EDI) highlighted participants that may be good candidates for control of dexterous prostheses i.e., large amplitude-based (median) EDI values suggest increased distinguishability of prompted missing hand grasp configurations. Additionally, the large EDI values for IQR measures indicated two findings. First, typical muscle excitation, which is indicative of distinguishable movements and therefore potentially effective control of dexterous prostheses. Second, variability may be too large across sEMG signals which may indicate poor consistency, upon which the muscle excitation plots and measures of consistency should be further investigated. Here, maximum and minimum EDI values are highlighted within the RMS or MNF measures for the unaffected and affected limbs across participants. SHR-G had the majority of maximum EDI values for 3 of the 4 measures of the unaffected limb and 2 of the 4 measures for the affected limb. In contrast, SHR-B had the majority of minimum EDI values for 2 of the 4 measures of the unaffected and affected limb, respectively. EDI values across limbs and measures for RMS and MNF are highlighted for each participant in Table 2.
Insert Table 2
To qualitatively visualize the differences across hand grasp configurations, the complete data set across trials (unsplit) was used to create the RDM for each limb. As previously mentioned, nonmetric MDS was used on each limb’s RDM to reflect the higher dimensional distances across hand grasps in a three-dimensional subspace. The MDS illustrated distinguishable differences between the various hand grasp configurations with few grasps close and/or overlapping in the subspace as seen by the separation of points in the top panel of Fig. 7a and b. Figure 7 shows both the median and IQR of RMS and MNF characteristics with the MDS (top) and corresponding RDM (bottom). Additionally, all MDS and RDM plots for participants can be seen in Supplementary Figs S19-S27 online. These results reveal, through a qualitative visualization, that the measures analyzed provide information to distinguish attempted hand grasps by visual separation.
Insert Fig. 7
4.3 Differences across limbs
4.3.1 Consistency across limbs
To get a better understanding of the differences in the participants’ ability to attempt consistent hand grasp configurations with their affected limb, statistical comparisons were made to the consistency of their unaffected limb. We found RMS consistency in the affected limb to be statistically lower than the unaffected limb for participants SHR-A, C, D, E, and H. Alternatively for MNF, participants SHR-C, SHR-E, and SHR-I had statistically lower consistency in the affected limb when compared to the unaffected limb. In each of the preceding cases, the unaffected limb had an overall higher median than the affected limb. The consistency for RMS and MNF across participants with highlighted statistical differences are shown in the top and bottom panels of Fig. 8, respectively. The resulting consistencies in RMS and MNF across limbs (Fig. 8) indicated that some participants had difficulty reproducing hand grasp configurations to a similar degree as that of their unaffected limb. This inability to reproduce grasps, in turn, may hinder their potential to use multi-grasp prostheses.
Insert Fig. 8
4.3.2 Relatedness across limbs
To understand if the structure of the affected limb RDM was related to that of the unaffected limb, the non-parametric RDM label randomization test was used34. It was found that the majority of participants had some degree of shared information such that the pair-wise distances between hand grasps for their affected limb was related to that of their unaffected limb. There was relatedness with the exception of the following participants: SHR-F (p = 0.060, RMS IQR), SHR-B (p = 0.055, median MNF), and SHR-F (p = 0.054, median MNF). Additionally, Table 3 shows the non-permuted Kendall’s Tau-a (\({\tau }_{a}\)) correlations used for the RDM test with corresponding p-values. These results indicate that the majority of participants had related information across limbs for amplitude and spread of RMS and MNF characteristics.
Table 3
Relatedness of RDMs across limbs for amplitude and spread of RMS and MNF characteristics.
Participants | RDM Relatedness Across Limbs |
RMS Correlation | MNF Correlation |
Median | IQR | Median | IQR |
SHR-A | 0.219 p = 0.002 | 0.266 p = 0.001 | 0.234 p = 0.002 | 0.288 p < 0.001 |
SHR-B | 0.401 p < 0.001 | 0.182 p = 0.015 | 0.119 p = 0.055* | 0.381 p < 0.001 |
SHR-C | 0.255 p = 0.005 | 0.159 p = 0.017 | 0.318 p < 0.001 | 0.293 p < 0.001 |
SHR-D | 0.272 p < 0.001 | 0.215 p = 0.005 | 0.180 p = 0.008 | 0.364 p = 0.001 |
SHR-E | 0.465 p < 0.001 | 0.262 p < 0.001 | 0.298 p = 0.001 | 0.254 p < 0.001 |
SHR-F | 0.200 p = 0.005 | 0.116 p = 0.060* | 0.123 p = 0.054* | 0.200 p = 0.003 |
SHR-G | 0.439 p < 0.001 | 0.349 p < 0.001 | 0.204 p = 0.002 | 0.301 p < 0.001 |
SHR-H | 0.224 p = 0.014 | 0.158 p = 0.013 | 0.423 p < 0.001 | 0.200 p = 0.005 |
SHR-I | 0.179 p = 0.007 | 0.144 p = 0.030 | 0.202 p = 0.011 | 0.137 p = 0.023 |
*Red cells highlight values from the RDM label randomization test, which indicate there was not sufficient evidence for the relatedness across limbs (with a significance level of α = 0.05). |
Insert Table 3