Analysis of the Autonomic Regulation in a Case of Facioscapulohumeral Muscular Dystrophy after Ken Ware Treatment

This is a study on autonomic neuroscience. In a previous paper in [1], we studied a subject affected from facioscapulohumeral muscular dystrophy before and after Ken Ware treatment (NPT). Using the non linear methodology of the Generalized Mutual Information (GMI) analysis of Sensory Motor Rhythm, we produced detailed results evidencing that the mentioned NPT treatment involved a net improvement of the patient under his subjective psychological condition, and in particular, under the neurological and sensory motor profile. We quantified with accuracy the improvement that the subject realized during such treatment. Of course, previous studies of several authors have evidenced that muscular dystrophies are strongly linked to a profound ANS disfunction. Therefore, the aim of the present study was to analyze the ANS of the subject before and after the treatment. We performed analysis in time as well as in frequency domain and by using non linear methods. The basic result of the paper was that, according to our analysis, the subjects started with a serious ANS disfunction before the NPT treatment and that a net improvement was obtained after this therapy. All the examined parameters resulted strongly altered before the treatment and all they returned in the normal range after the NPT.


Introduction
When speaking about Heart Rate Variability (HRV), we have to start considering such basic foundations.Baroreceptor activity is received from the Nucleus Tractus Solitarius (NTS) that provides to transmit on one hand the signal to Nucleus Ambigous/Dorsal Motor Nucleus (NA/DMN) integrated by the hypothalamus action which generates finally parasympathetic heart rhythm modulation and on the other hand the signal to Rostral Ventrolateral Medulla oblongata (RVLM) that integrates hypothalamus contribution generating orthosympathetic heart rhythm modulation.A complete scheme of the complex arrangement may be found in [2] On the basis of this scheme, the HRV may be considered as an analysis that on one hand gives important information about the fitness of heart to respond to the autonomic nervous system (ANS) modulation, and on the other hand (and it is here the important feature relating the present paper) gives direct indication on the manner in which some central brain structures of dominant relevance are explicating their basic function when inducing heart rhythm modulation.In this sense, we may say that HRV dynamics is the "mirror image" of the psychophysiological mechanisms involved in some central brain regions.By this way, we conclude that HRV studies represent currently a non invasive marker not only to explore the cardiovascular system but also to indirectly investigate some important neurological involved functions.
To support this thesis, we have to quote here a recent neuroimaging study [3].This study was conducted by Julian F. Thayer, Fredrik Åhs, Mats Fredrikson, John J. Sollers and Tor D. Wager who evidenced the existing main relationship between Hear Rate (HR) and regional brain blood flow.These authors indicated in detail that the more significant involved brain regions were represented from the amygdala and from the ventromedial prefrontal cortex.In detail, these authors identified the medial pfc (MPFC) (right pregenual cingulate), the MPFC (right subgenual cingulated) and the SLEA (left sublenticular extended amigdala/ventral striatum).
Our conclusion is that such results evidence the importance of the HRV studies in exploring cognitive as well as emotive functions in clinical conditions of subjects affected from psychological disorders as well as psychiatric pathologies.This is not a restriction of course since, as we will attempt to evidence in the present paper, the heart rate variability analysis may be considered of importance also in the study of subjects affected from neuromuscular pathologies.
In 2003, we started to analyze HRV in different experimental conditions.And in the course of several years of studies we have also introduced novel methods [4]- [6] of analysis in order to obtain valuable quantifications of the VLF, LF, and HF bands where in particular the quantification obtained in the LF band is related mainly to the modulation due to the orthosympathetic activity while the HF band relates more directly parasympathetic activity and respiration.
The aim of the present study is to evaluate HRV and thus ANS disfunction in a subject affected from Facioscapulohumeral Muscular Dystrophy and to estimate its improvement after that the subject has been submitted to a treatment that we have called NPT that is due to Ken Ware and that we have explained in detail recently elsewhere [1].
We may reassume briefly the reason to investigate here the ANS disfunction before and after Ken Ware treatment.As previously outlined in [1], Duchenne (DMD) and Becker (BMD) muscular dystrophies are disorders due to complete or partial loss of dystrophin protein.The presence of persistent sinus tachycardia, atrial and ventricular arrhythmias, sweating and chills has suggested an ANS involvement in these disorders.Yotsukura et al. first described an impairment of ANS in 55 DMD, characterized by an increase in sympathetic activity and a reduction in parasympathetic output [7] [8].ANS impairment was confirmed also by studies that were developed by Lanza et al. [9] and by Inoue et al. [10] [11].In a study developed from Vita in 2001 [12], 20 BMD patients were investigated with a battery of six cardiovascular autonomic tests and power spectral analysis of HRV.Although 11 patients revealed abnormal findings at some cardiovascular tests, none of them had a definite autonomic damage, as indicated by two or more abnormal tests.This clinical observation suggested therefore that autonomic involvement did not represent a major finding in BMD.In 2006, the prognostic value of HRV for sudden death was evaluated in a population of 30 BMD patients and dilated cardiomyopathy compared with a control group of normal subjects [13].The important result was that an increment of sympathetic tone in BMD was observed, characterized by a lower HRV in the frequency and time domains and a higher mean heart rate than in the control group.
When distinguishing the case of the myotonic dystrophy (DM) in the myotonic dystrophy type 1 (DM1) and type 2 (DM2) with different clinical phenotypes, we find that in such conditions we have indication of ANS in-volvement such as heart rhythm conduction disturbances, orthostatic hypotension, disorders of sweating, abnormal gut movements, and dysfunction of genitourinary apparatus as studied in particular from [14].Very important results may be found in [10] [15]- [26].
Facio Scapulo Humeral Muscular Dystrophy (FSHD) is the third most frequent form of muscular dystrophy.In [27], HRV was evaluated in 55 FSHD patients and evidenced the presence of autonomic modifications, characterized by a slight increase in sympathetic activity and a progressive decrease in parasympathetic output, which became more evident with the progression of the disease.
Emery-Dreifuss Muscular Dystrophy (EDMD) is genetically heterogeneous since it may be caused by a mutation in the STA gene encoding emerin (X-linked inheritance) or as AD trait determined by mutations in the LMNA gene encodinglamin A/C.Scintigraphic studies [28] evidenced a diffuse and severe decrease in accumulation of MIBG in the heart, suggesting an abnormality in the cardiac sympathetic nerve terminals in the patient.
Miyoshi muscular dystrophy (MMD) is caused by mutation in the gene encoding the protein dysferlin.In 1994, Tomoda et al. [29] performed accurate autonomic studies in 2 MMD affected girls, including laser doppler flowmetry, study of the component analysis of the cardiographic R-R interval and sympathetic skin response (SSR).They demonstrated marked abnormalities such as sensitive vasoconstrictive response, a suppressed peak of low frequency components, and an absence of SSR, compared with healthy controls.
We retain to have delineated a sufficient picture of the importance of evaluation of HRV and ANS disfunction in muscular dystrophy.The present paper is a case report.We give the results on ANS renormalization in a severe case of muscular dystrophy after the treatment that is introduced by Ware and discussed and analyzed in detail by us elsewhere [1].

Materials, Methods and the NPT Treatment
As it was outlined by us in detail in [1] the case is of a subject 60-year-old-male suffering for more than 30 years from the rare form of Muscular Dystrophy, Facioscapulohumeral Muscular Dystrophy (FSHD).(FSHD) is a genetic muscle disorder in which the muscles of the face, shoulder blades and upper arms are among the most affected.EEG, ECG, EMG were all recorded simultaneously using BioRadio 150 technology.EEG was recorded in channel 8 with electrodes placed at C3-F3 of the left hemisphere.ECG was recorded in Channel 7 with electrodes placed at V1 and V2.An appropriate software by BioPac sytem enabled us to obtain the relative R-R intervals for HRV analysis.The patient was recorded for 5 minutes at rest before therapy in order to perform quantitative estimations before and after the treatment.The exercise equipment used is specialized exercise equipment supplied to the International NeuroPhysics Functional Performance Institute by "Synergy Fitness (http://synergyfitness.com.au/gallery/15/Synergy_4_Pin_Loaded) and was described in detail min [1].In substance this treatment that we have called the Ken Ware NeuroPhysics Treatment (NPT) does not involve physical manipulation of the treated individual.Usually early phase equipment is the leg press (which may be substituted with the seated leg curl), lateral pull down, and pectoral exercise, the pec dec.Each exercise machine loaded with very light weight, mostly representative of the equipment's minimum lifting loads.Most movements have only two degrees of freedom (df) (up-down or forward-back).Only one machine used in the early phase of treatment has 4 df, namely the pec dec.The initial treatment in our case, involved four models of variable resistance exercise equipment, which were selected in consideration to the level of ability of the patient.No stretching or warm-up exercises are permitted to preserve a patient's system normal resting state, innocent of perturbations until the (NPT) starts.The initial assessment session(s) evidences the patient's chronic system condition and how it interacts with and evaluates its environments in general.The specific activities and equipment used at the beginning of treatment are expanded after initial sessions' effects area assessed and systematically modulated as necessary, through the unique applications and qualities of NPT.In the next treatment phase, additional exercises that are slightly more complex to perform, but with still very lightly weighted, are added to the patients program, again depending on ability or disability.Some of these pieces of equipment have many (df); the handles are attached to cables that can move in any direction.The patients are at this time (re) introduced to socalled "cardio" exercises (this nuisance term promoted by the exercise industry un-holistically implies the cardiovascular system is separate, requiring its own form of exercise).In the reinforcement and maintenance phase of the program, additional equipment with prescribed repetitions and weights is introduced.Individuals in trea-t ment perform the following activities with trained coaching support).In conclusion, we may identify six phases of the treatment.
The first phase is devoted the subject to learn to relax physically, mentally, and emotionally.The second phase is finalized to adjust position for proper contact points and symmetrical balance on equipment.
The third phase is devoted to maintain up right posture.
In the fourth and fifth phases, the subject closes the eyes and performs ultra-slow movements.Finally, the sixth phase is devoted to form and to operationalize the intention to distribute effort and energy throughout the body (pec dec).For further detail the reader is addressed to ref. [1].

The Results
As previously outlined, we performed HRV analysis starting with an ECG recorded at 960 Hz.We used linear and non linear methods.As previously discussed, studies in literature have evidenced relationship between lowfrequency power (Power Spectral Density, PSD) of heart rate variability and muscle sympathetic nerve activity (MSNA).We studied HRV before, during the treatment, and days after the treatment.
By using Fourier transform (FFT or DFT), one obtains balancing inspection of the ANS modulation in HRV in the three VLF, LF and HF bands.The VLF Band relates hormonal activity, thermoregulation and possibly the rennin-angiotensin system, the LF band relates mainly the orthosympathetic activity and, finally, the HF band relates instead the parasympathetic activity in the frequency domain.
Let us explore the DFT results evaluating the ANS condition of the subject at rest before starting the treatment.We identify this at rest section as Lat Pull1bt.Let us give a look at the tachogram and to the DFT spectrum as they are reported in Figure 1.
Figure 1 evidences that we are in presence of a subject with very exaggerated beat-to-beat variations and presence of an arrhythmia to be classified.In the PSD spectrum peaks are represented in the three bands all having abnormal amplitudes in the three bands respectively.
Let us give a look at the results of the quantitative analysis.The bpm (beats per minute) of the subject resulted to be 68.7 bpm where we consider a normal range included between 59 and 96 bpm.Soon after we examined the variance as it results from the reconstructed tachogram starting with the given R-R time series.Usually, we attribute a great relevance to such index since it informs us directly on the variability induced on the heart rhythm by the ANS modulation and on the availability of the sinus node to respond to ANS activity.The value resulted of 135.5 mseconds and we adopt as normal range the interval 20 -60 mseconds in the age range of 50 -70 years old.This result indicated an exaggerated variability in heart rhythm induced from some kind of arrhythmia.The toal Power (PSD) resulted to be 18,397 msec 2 /Hz with natural logarithmic value of 9.81.Usually, we consider as normal range the interval 6.0 -8.1.Therefore, we concluded for the presence of an exaggerated modulation total value exceedingly the normal prefixed value.The VLF band gave as result the value of 3207 msec 2 /Hz with Ln value of 8.07 and the normal range is considered usually to be between 4.5 and 6.4.The obtained value resulted exaggerated exceeding the established normal range.
The LF band gave us a value of 7565 msec 2 /Hz with an Ln value of 8.93 while the normal range is included between 5.0 and 6.9.Finally the HF band gave us a result of 6956 msec 2 /Hz with an Ln value of 8.84 in the normal range included between 4.5 and 6.7.Also all the linear time indexes resulted exaggerated out of the normal range.They are reported in Figure 2.
The context previously delineated configures a net ANS disfunction but in order to deep the problem we need to add further investigations.This was the reason to apply to the data the CZF method that we introduced in literature as explained in detail in [4]- [6].The CZF method enables us to develop a detailed analysis of the hearth rhythm variability in the frequency domain evaluating in the detail the variability induced from orthosympathetic and parasympathetic modulating activities respectively.The results are reported in the Figure 3 and Figure 4.
The case of the CZF application with absolute value quantifies variability between each pair of subsequent value taking the result in mseconds while the case without absolute value quantifies such variability but accounting if the subsequent time beat value is greater or les of the previous one.The first estimates the value of the beat-to beat difference.The second approach estimates in addition if such interval increases (acceleration of beat to beat interval) or decreases.
In red we have time dynamics and on the right we Fourier analysis of the heart rhythm variability.
The important result that we obtain with this kind of analysis is that it evidences that in the LF band we have a strong reduction of variability and this is to say that we have a strong dominance of sympathetic activity while the HF band that of course characterizes the ANS modulation of the parasympathetic activity remains in the normal range.Therefore, the results appear to be unquestionable.The ANS heart rhythm modulation of this   subject is affected from an hyper orthosympathetic activation.This is in agreement with the results of other authors as we discussed in detail in the introduction.
At this point it is instructive to examine also the histograms computed one time for the heart rate (bpm) and one time for the R-R intervals.We give the results in Figure 5 and Figure 6 respectively.The analysis HR histograms also reveals a situation of the subjects affected by a profound arrhythmia with HR values violating a normal Gaussian distribution.Mean, Median and Mode are out of the normal requirement in a normal distribution and we observe substantially a distribution peaked about three central values, a peak in HR about 76 (bpm), another peak about 128 (bpm) and a final distribution about 180 (bpm) The same evaluation is obtained obviously when we consider the distribution of the R-R values as reported in Figure 6.We have substantially a plurimodal distribution that is always indicative of a serious happening arrhythmia in the examined subject.
To complete the analysis we may apply the Poincaré method.We give the results in Figure 7.This is a well know representation in bi-dimensional phase space.Also its interpretation, obtained on the basis of the distribution of the points and of the estimated values for SD1, SD2, ellipse area and centroid, is well known.As confirmation also of the results previously obtained, this is a clear case of atrial fibrillation.
In detail, the examination of the Poincaré plot given in Figure 7 in conjunction with the results previously given, enables us to classify the kind of arrhythmia of the subject.In the plot RR n /RR n+1 we identify two locations out of the ellipse, one located in the range 310 -496 ms and the other in the range 748 -868 ms.The Sd1 value is reduced to about 51 ms or 9.17 n.u.(orthosympathetic activity) against the 284.63 ms or 46.66 n.u. for SD2 with SD1/SD2 = 0.1966 strongly distant from a condition of acceptable balance parasympathetic/orhosympathetic modulation.The Ellipse Area is strongly out of a normal range.This picture enables us to conclude that the subject is suffering of a persistent arrhythmia to be classified as atrial fibrillation.
We may now conclude our examination by using non linear analysis, estimating in detail Approximate Entropy (ApEn), Sample Entropy (SnEp), and the Detrended Fluctuation Analysis (DFA).They are given in Table 1.
The results are given in Table 1.It is confirmed the profound condition of great cardiovascular risk of the subject and his serious ANS disfunction.This condition is clearly evidenced when considering the values of ApEn and SampEn that are so distant from a valuable value near to 1. Fortunately long as well as short range correlations as estimated by the DFA still remain in the range of an acceptable value.
We may now pass to examine the question of the treatment of the subject connecting in this manner the great  The subject has now 72 bpm.The most important thing is that the variability is now of 30 mseconds that is in the normal range.The total power (PSD) is now 902 msec 2 /Hz (Ln = 6.8) (against the value at rest of 18,397 msec 2 /Hz) and is now in the normal range.The VLF band gives a value of 464 msec 2 /Hz (Ln = 6.13) and is in the normal range, the LF band is 313 msec 2 /Hz (Ln = 5.74) and is now in the normal range, the HF band is 94 (ln = 4.54) and still is in the normal range.In Figure 8, we report the tachogram and the Fourier spectrum and both the behaviours re-enter now in the normal regime.
Let us start examining the results obtained by Fourier Transformation and given in Figure 7.
Looking at the Figure 8, we may appreciate under the neurological profile the valuable improvement induced from the NPT treatment.The ANS modulation has totally recovered.Also the values of the time linear indexes now result in their normal range as shown in Figure 9. Also the values of the CZF method result to be totally in the normal range.It is very interesting to examine the histograms and the Poicaré Plots.The histograms are given in Figure 10 and Figure 11, HR and R-R intervals respectively.
The histograms follow now the expected normal distribution that we observe in a subject with missing disfunctions.
Evaluating all such results, we may conclude that the ANS has totally recovered its modulating function.
Let us look at the Poicaré Plot in Figure 12.
Looking at Figure 12, we observe that the subject has totally recovered the characteristic distribution at cigar of normal subjects also if still remain some extra-distribution episodes and the form of the cigar is not so largely extended as in a health normal subject.Sd1 and SD2 values are however in the normal range of values as well as the values of the ellipse area and centroid.In detail, the SD1 and the SD2 values, relating the ANS modulation, enter now, as said, in the range of normal values with a correct balancing parasympathetic-orthosympathetic action, that in fact results to be Sd1/SD2 = 0.848 (against the value of 0.193 at rest before of the treatment) with SD1 = 43.94ms and 5.22 n.u. and SD2 = 51.79ms and 6.15 n.u.
Finally, from a neurological point of view, the most interesting result is given by using non linear indexes.Let us report the results relating the Approximate Entropy (ApEn), the Sample Entropy (SnEp) and the Detrended Fluctuation Analysis (DFA).They are given in Table 2.
All the results are now in the normal range.

Conclusions
This is a case report of ANS analysis in a subject affected by muscular dystrophy.The HRV was analyzed by using linear and non linear analysis.The analysis performed initially on the subject at rest and before of the treatment by the linear time indexes and traditional FFT resulted out of the range.The use of the CZF method, histograms and Poincaré plot evidenced that the subject started with a strong ortho-sympathetic against parasympathetic unbalance and atrial fibrillation.After the treatment, the patient evidenced that all the standard HRV indexes re-entered in the normal range of values in consideration of his starting arrhythmia, and that in particular balancing orthosympathetic-para-sympathetic correct modulation was recovered.
Under the neurological profile, we have to outline here the profound correlation existing in autonomic regulation and muscular dystrophy.In particular, the use of standard methodologies in HRV analysis, as well as the use of the CZF method and of non linear estimations, enables us to estimate and to correlate the ANS dysfunction with the severity of the muscular pathology.In fact, the results on EMG analysis, published elsewhere [1] and      relating this subject, evidenced that ha had a so serious initial compromission that of course resulted to improve under the profile of the clinical and pathophysiological conditions after the treatment.

Figure 2 .
Figure 2. Linear time indexes of the subject at rest.

Figure 3 .
Figure 3.Time series of heart rhythm variability and values of variability (msec/Hz) in the case with absolute value (for details see [4]-[6]).The CZF method.

Figure 4 .
Figure 4. Time series of heart rhythm variability and values of variability (msec/Hz) in the case without absolute value (for details see [4]-[6]).The CZF method.

Figure 5 .
Figure 5. Histograms of HR for the subject at rest.

Figure 6 .
Figure 6.Histograms of R-R intervals for the subject at rest.

Figure 7 .
Figure 7. Poincaré Plot relating the subject at rest.

Figure 9 .
Figure 9. Values of Linear time indexes in HRV analysis after the treatment.

Figure 10 .
Figure 10.HR-Histogram distribution of the HR values after NPT treatment.

Figure 11 .
Figure 11.R-R Histogram distribution of the RR values after NPT treatment.

Figure 12 .
Figure 12.Poincaré Plot of the subject after NPT treatment.

Table 1 .
Estimation of approximate entropy, of entropy, and detrended fluctuation analysis in the subject at rest.