Abstract
Background: Mobile phones may have harmful health effects and clinical examinations report ambiguous results of exposure concerning neurophysiological and cardiovascular actions. Materials and Methods: This study investigated heart rate asymmetry (HRA) and heart rate variability (HRV) parameters with 1:2 and 1:1 metronome-paced inspiration/expiration ratios during short-term 1,800MHz GSM cellular phone exposure in 20 healthy volunteers. Results: Significant HRA changes by Porta and Guzik indices were not found on exposure compared to sham exposure. Time-domain HRV parameters on exposure showed significant differences at 1:1 paced, but not at 1:2 paced breathing compared to sham exposure. A mild post-exposure effect was observed regarding root mean square of successive RR-differences. Conclusion: The findings reflect persisting acute effects of GSM handset emission on the autonomic nervous system. Exploring its influences on health status and survival needs further studies. Symmetrical breathing can be used as a sensitizing factor in other HRV/HRA analysis studies.
Evolution of species has taken place in the presence of natural ionizing- and non-ionizing electromagnetic radiation, including sunlight, cosmic radiation, natural radioactivity and atmospheric electromagnetic phenomena, among others (1). There is increasing concern regarding the potentially harmful biological effects of exposure to electromagnetic fields (EMF), especially on humans, from military, industrial and commercial wireless telecommunication systems. Mobile phones have become popular and indispensable in our everyday life. The first generation analogue network was replaced in 1991 by the digital Global System for Mobile Communication (GSM), which had more than 90% market share still in 2014 in spite of the introduction of the third-generation Universal Mobile Telecommunication System (UMTS) in 2001 and the fourth-generation Long-term Evolution (LTE) in 2009, with increasingly improved transmission capacity (2). All three generations of mobile communication standards are simultaneously present in Hungary in 2017 (3).
The mobile communication system consists of a cellular network of base and mobile stations (cellular phone or handset). The latter emits significantly less microwave power, however, its long-term close proximity to the head and torso results in more significant exposure in the general (not professional) population. The possible biological effects of microwave radiation are related to the energy absorption by living cells and their interactions at the molecular, cellular, tissue and organism level. Since the relative permeability of biological materials is around unity, independently of the exposure frequency, whereas their relative permittivity is 10-200, the magnetic rather than electric component of EMF may be responsible for any biological effects (4).
The direct biophysical mechanism of the interaction of EMF and living tissues is still not clear; however, there are several theories (1). Membranes play a crucial role in biology by ensuring compartmentalization and interaction of the inner and outer cellular spaces at the same time. The interplay between membrane potential and membrane conductivity regulates cell functions, and excitation of neurons and contractile cells among others. A direct membrane-polarizing effect of microwave radiation can be excluded since the average capacity of the cell membrane (1 μF/cm) is considered as a short circuit above ~1 MHz, hence there is no voltage drop and energy dissipation induced across the membrane itself (5). However, a low-frequency modulated radiofrequency field can act differently. Moreover, considering the resting membrane potential of −70mV results in an electric field gradient across the 5 nm-thick cell membrane of E≈14×106V/m, the external EMF used in communication are negligibly weak to be able to directly act on membrane potential (6). Furthermore, the quantum (ionizing) effect of microwave radiation at the molecular level can be disregarded considering the energy carried by photons with a frequency of 0.300-300GHz would be 1.24×10−6 to 1.24×10−3 eV. This amount of energy is well below the ionizing threshold of biological molecules (~13.6 eV) and even smaller than the 0.027eV energy of Brownian motion at 37°C (4). Today the most widely accepted and explored theory for interaction is based on the thermal effects of microwave radiation. The current safety limits are established based on this; however, there are several biological phenomena that cannot be explained simply by warming of tissues, and temperature changes cannot be detected at all behind some observed biological effects (1, 7). The physical basis for thermal effects is excitation by rotational resonance of dipole molecules induced by sinusoidal EMF, which is resisted by the molecular moment of inertia and friction with the surrounding particles, which increases the thermal energy of the object, resulting in its warming up. In the microwave range, water and similar small dipolar molecules are the target of resonant frequency, although conformational changes of proteins via their bound water or molecular arms can also be considered targets. Depending on the inhomogeneity of the tissues, local hot-spots can arise, and due to the temperature-dependent nature of relative permittivity, focal thermal breakdown can occur (4).
The present applicable standard of maximal transmission power of a hand unit is 1.0 W at GSM 1800/1900 and 2.0 W at GSM 850/900 (8). The European Union specific absorption rate (SAR) limit for the public is 2.0 W/kg averaged over 10 g of tissue, but in a few countries the limit is 1.6 W/kg averaged over 1 g of tissue (9).
Exposure to EMF may have negative health effects even below the power defined in safety standards. Investigations into possible adverse biological effect of cellular phones are focused on the central nervous system (CNS) and inner ear, because most of the microwave radiation is absorbed here due to the anatomical proximity to the signal source (10-15). Conclusions of human trials on the impact of EMF on the autonomic nervous system are controversial (16-19).
HRV analysis is a non-invasive assessment of the actions of the autonomic neuroendocrine system on the sinus node via the beat-to-beat fluctuations of the heart rate. This method is considered valuable in the prediction of progression and outcome in several diseases (20-22). Atlas et al. found no difference in HRV in 35 healthy volunteers breathing spontaneously during 10-minute exposure to EMF from a commercial mobile phone forced at 2 W output peak power at 900 MHz (23). Heart rate asymmetry (HRA) is a relatively novel marker of HRV, based on the Poincaré-plot. This asymmetry can be quantified as Guzik index (24) or Porta index (25), among others. HRA reflects the dynamics of respiratory arrhythmia, namely temporal asymmetry of pulse rate acceleration and deceleration as a consequence of inspiration/expiration period ratio, first published by our research group in 2012 (26): Change in inspiration/expiration ratio by double-paced breathing from 1:2 to 1:1 or 2:1 is followed by corresponding change of HRA parameters. HRA measures can be more sensitive markers of respiration–heart rate coupling than standard HRV parameters; additionally, paced breathing ‘stress-situation’ requiring auditory and cortical functions also acts as a sensitizing factor.
In the available literature, there is no publication on the effect of cellular phone exposure on HRA parameters and at different paced inspiration/expiration period ratios. The aim of present study was to examine the influence of the acute effects of pulsed microwave irradiation from a commercial cellular phone on HRV and HRA parameters during various double-paced breathing patterns in healthy volunteers in randomized double-blind repeated-measures crossover design.
Materials and Methods
The study was approved by the Regional Research Ethics Committee (approval number: 2013/4747) and was conducted in accordance with the Declaration of Helsinki and its later amendments. Twenty healthy volunteers (14 female and six male) were enrolled. Informed consent was obtained from all individual participants included in the study. The mean age was 25.2 (range 21 to 32) years. The majority of the volunteers did not smoke (17 altogether), two people smoked fewer than 10 cigarettes/day, and one person smoked more than 10 cigarettes/day. Sixteen volunteers consumed alcohol occasionally, while four did not consume alcohol at all. Regarding coffee consumption, the population was heterogeneous: six did not consume at all, seven occasionally and seven regularly drank coffee. The mean body mass index (BMI) of the group was within the normal range: 21.8±2.5 kg/m2. The participants were not allowed to consume alcohol, coffee or cigarettes for 4 hours before the test.
A commercial Nokia 6230i mobile telephone (Nokia Corporation, Helsinki, Finland) was used for radiofrequency exposure working on an 1,800 MHz GSM network (217 Hz pulse rate; 0.577 μs pulse width, actual emissions during the experiment were not measured) with SAR value of 0.70 W/kg at the head (27). The call to the mobile phone was initiated from a landline unit in a neighbouring room until automatic disconnection at about 1 minute and 25 sec, and then the number was redialled three more times over a 6-minute period. The participants had no information on the emission of the phone since it was in silent mode and was totally covered by an elastic bandage.
During each experimental stage, a 6-min-long single channel electrocardiogram (ECG) and respiration signal were recorded by a microcontroller-based, battery-powered, handheld data acquisition system (László Hejjel, Pécs, Hungary). An ultra-rapid thermistor probe was fixed under the volunteer's nose for monitoring respiration. The thermistor sensed the cooling and warming of the flowing air during inhalation and exhalation, respectively. The instrument also contained a metronome to trigger inspiration and expiration with two different frequency beeps. Standard 300-second tachograms (22) were extracted from the ECG records by ECGRdet v2.4 (László Hejjel), while the HRV analysis was carried out with Varian v2.2 (László Hejjel). The mean time between two R-waves on the ECG (RR interval; MeanRR), standard deviation of normal-to-normal RR-intervals (SDNN), coefficient of variation (CV%=SDNN/MeanRR), root mean square of successive RR–interval differences (RMSSD) in the time domain (21, 22); high-frequency (HF) integral in the range of 0.15-0.40 Hz (in normalized units), low-frequency (LF) component at 0.01-0.15 5Hz (in normalized units) and LF/HF were analysed by fast Fourier transformation in the frequency domain (21, 22); Porta index and Guzik index (24, 25) were also computed as HRA parameters.
Box and whisker diagrams of the time domain parameters. A: Mean time between two R-waves on the ECG (RR interval). B: Standard deviation of RR intervals (SDNN). C: Coefficient of variation (CV%). D: Root mean square of successive RR interval differences (RMSSD). Boxes: Lower and upper quartiles, horizontal line: median value, whiskers: minimum and maximum values, square: mean value. *Friedman test p-value, **post-hoc Wilcoxon's paired-sample test p-values, only significant differences are indicated. Insp: Only inspiration paced without exposure; 1:1 symmetrically paced inspiration/expiration; 1:2 asymmetrically paced inspiration/expiration; i: with irradiation; n: without irradiation.
The study was carried out at our Institution in 2014 in the late afternoon, in quiet circumstances. The volunteers lay in a comfortable supine position with head elevated at 30° 15 minutes prior to data acquisition to allow orthostatic adaptation. The ECG electrodes along with the thermistor probe were placed on the subject, and the mobile phone was fixed on their right ear with an elastic bandage imitating its position during use. The subjects were educated regarding metronome-controlled respiration during this adaptation period. Talking and movements were not allowed, since these can influence HRV parameters (28).
Five consecutive 360-second-long stages were recorded for each volunteer. A T=4.5s breathing cycle (0.22 Hz) was ensured by participants breathing according to the built-in beeper. Since previous studies have already demonstrated that single metronome-controlled versus spontaneous breathing does not differ as regards of HRV parameters, we did not investigate spontaneous respiration here in order to reduce the study duration and hence the stress on the participants (29, 30). During the first stage, only inhalation was according to the metronome. In the second and fourth stages, inhalation and exhalation were triggered with beeps of two different frequencies, in an asymmetrical manner at a 1:2 (physiological) inspiration:expiration ratio. During the third and fifth stages, symmetrical breathing was paced, i.e. 1:1 inhalation:exhalation ratio. The radiofrequency exposure from the mobile phone took place randomly either during the second and third, or during the fourth and fifth stages, blinded to both the volunteers and the person who analyzed the records. Valentini et al. stated in their meta-analysis that GSM-irradiation causes only transient cortical phenomena by EEG studies, therefore only the instantaneous effects of radiofrequency exposure were analysed (10). For this reason and in order to minimize the mental load of our investigation, only two 6-minute exposures were applied with 1:2 and 1:1 inspiration/expiration ratios in contrast to the frequently applied 30-minute microwave exposure. The entire measurement needed about 45 minutes per volunteer.
Box and whisker diagrams of the frequency domain parameters. A: High frequency power in normalized units (HFNU). B: Low frequency band power in normalized units (LFNU). C: Low frequency/high frequency power ratio (LF/HF). Boxes: Lower and upper quartiles, horizontal line: median value, whiskers: minimum and maximum values, square: mean value. By Friedman test p=0.107 for each frequency domain parameter, no significant changes. Insp: Only inspiration paced without exposure; 1:1 symmetrically paced inspiration/expiration; 1:2 asymmetrically paced inspiration/expiration; i: with irradiation; n: without irradiation.
Statistical analysis was prepared using the statistiXL v1.8 (statistiXL, Broadway Nedlands, Western Australia) plug-in in Microsoft Excel (Microsoft Corporation, Redmond, WA, USA). The figures were constructed with OriginPro v2017 (OriginLab Corporation, Northampton, MA, USA). Non-parametric repeated-measures ANOVA by Friedman was applied for primary assessment. Values of p≤0.05 were considered statistically significant. Wilcoxon paired-sample tests with the Holm-Bonferroni correction of significance level were applied post-hoc. In pairs, the following were compared: A: only inspiration controlled by metronome (Insp) with sham-exposed (i.e. not exposed) asymmetrically double-triggered breathing (1:2n); B: 1:2n pattern with the pattern of non-irradiated symmetrical (1:1n); C: irradiated asymmetrical (1:2i) with irradiated symmetrical breathing (1:1i); D: 1:2n with 1:2i; and E: 1:1n with 1:1i.
Results
On manual checking of the ECG and respiration recordings, no records were rejected due to poor signal quality, inappropriate breathing, arrhythmia or any other reason. Moreover, spectral analysis by fast Fourier transformation clearly showed the respiratory peak at around 0.22 Hz for each individual.
Time domain analysis. The mean RR-interval changed significantly according to Friedman's test. Wilcoxon's paired-sample test post-hoc showed a significant difference (p=0.002) between the Insp and 1:2 paced sham-exposed measurements even after Holm–Bonferroni correction. Further comparisons of the mean RR-interval showed no significant changes (Figure 1A).
The Friedman test of SDNN parameters was statistically significant (Figure 1B), post-hoc paired-sample assessment resulted in significant differences even after p-value correction between 1:2 and 1:1 breathing measurements under EMF exposure (p=0.009). The other comparisons regarding SDNN demonstrated no significant changes. Interestingly, the mean and median of SDNN were separated at 1:1 breathing with EMF irradiation. The CV% behaved similarly giving a statistically significant p=0.008 value on comparing 1:2i with 1:1i series by post-hoc Wilcoxon's paired-sample test (Figure 1C).
The RMSSD resulted in significant differences by Friedman test (p=0.022) and again, only 1:2i with 1:1i series gave significant differences by post-hoc Wilcoxon's paired-sample test (p=0.008, Figure 1D). The increase of the interquartile distance is evident here as well, representing a possible separation of the population during cellular phone EMF exposure at symmetrical but not asymmetrical breathing together.
Box and whisker diagrams of heart rate asymmetry (HRA) parameters. A: Porta index. B: Guzik index. Boxes: Lower and upper quartiles, horizontal line: median value, whiskers: minimum and maximum values, square: mean value. *Friedman test p-value, only significant differences are indicated; the examined post-hoc Wilcoxon's paired sample tests resulted in non-significant p-values after Holm-Bonferroni correction. Insp: Only inspiration paced without exposure; 1:1 symmetrically paced inspiration/expiration; 1:2 asymmetrically paced inspiration/expiration; i: with irradiation; n: without irradiation.
Frequency domain analysis. Among frequency domain parameters (LF, HF, LF/HF) no statistically significant differences were found by Friedman test (Figure 2). Due to the complementary nature of LF and HF power and the ranking in Friedman test, the p-values proved to be identical.
HRA analysis. The Friedman test resulted in statistically significant p-values for both Porta (p=0.016; Figure 3A) and Guzik (p=0.009; Figure 3B) indices. However, post-hoc Wilcoxon's paired-sample test of Porta index demonstrated only a tendency for differences between 1:2 and 1:1 breathing patterns in both comparisons sham-exposed and irradiated pairs (p=0.011 for 1:2n vs. 1:1n; p=0.024 at 1:2i for 1:1i) after Holm–Bonferroni correction of the significance level. Wilcoxon's paired-sample test for Guzik's index showed similar results: p=0.027 for 1:1n vs. 1:2n; p=0.011 for 1:1i vs. 1:2i. There were no significant differences in either index between sham-exposed and irradiated groups comparing the same breathing patterns: 1:1n vs. 1:1i and 1:2n vs. 1:2i (Figure 3). Both Pearson's correlation and Spearman's rank correlation confirmed a significantly strong positive correlation above 0.9 between Porta and Guzik indices with p<0.001 and p=1.022×10−18, respectively (Figure 4).
Pearson's correlation and Spearman‘s rank correlation tests between Porta and Guzik indices. The tests showed a significantly strong positive correlation. Pearson correlation: rP=0.932, p<0.0001; Spearman rank correlation: rS=0.948, p=1.022×10−18.
Recombination of the time domain measurements. In order to reject or confirm order-dependent effects, the five stages were regrouped in the order of data acquisition ignoring irradiation: 11 volunteers were exposed in the 2nd and 3rd stage, while nine were exposed in the 4th and 5th stage. Figure 5 demonstrates that there were no significant differences in SDNN among the groups with identical breathing patterns independently of the EMF exposure. Between the 1:1 and 1:2 groups, SDNN in the 2nd and 3rd stages differed, p=0.002 is considered a statistically significant difference even after Holm-Bonferroni correction.
Standard deviation (SDNN) changes according to the order of measurement and breathing pattern. Boxes: Lower and upper quartiles, horizontal line: median value, whiskers: minimum and maximum values, square: mean value. *Friedman test p-value, **post-hoc Wilcoxon's paired-sample test p-values, only significant differences are indicated. Insp: Only inspiration paced without exposure; 1:1 symmetrically paced inspiration/expiration; 1:2 asymmetrically paced inspiration/expiration; i: with irradiation; n: without irradiation.
Exploring possible post-exposure effects, the 2nd and 3rd irradiation (n=11) and 4th and 5th stage irradiation (n=9) groups were analyzed separately according to exposure versus sham-exposure, given that RMSSD reflects parasympathetic influence on heart rate. The 2nd and 3rd stage exposed groups gave the chance to examine post-exposure effects up to 12 minutes (during the 4th and 5th stages) following a 12-minute EMF exposure. During exposure (n=11+9) and before exposure (n=9), the 1:1 breathing resulted in greater median and mean RMSSD than with 1:2 breathing (Figure 6). Paradoxically, the post-exposure analysis (n=11) showed lower mean and median RMSSD with 1:1 breathing. Due to the small number of people in the split groups, we did not perform further statistical analysis. However, this may be the explanation for the significant differences according to breathing patterns only between the 2nd and 3rd but not between the 4th and 6th stages apparent in Figure 5.
Root mean square of successive differences (RMSSD) changes after separation of the first-exposed and second-exposed volunteers. Boxes: Lower and upper quartiles, horizontal line: median value, whiskers: minimum and maximum values, square: mean value. Horizontal axis with the breathing patterns, the serial number, and the presence of exposure. Insp – Only inspiration paced without exposure, 1:1 symmetrically paced inspiration/expiration; 1:2 asymmetrically paced inspiration/expiration; i: with irradiation; n: without irradiation, bottom line: stage serial number.
Discussion
Additionally to electric inhomogeneity and anisotropy of the human body, local EMF effects are influenced by the wavelength relative to the target size, penetration/damping and intracorporal reflections of the microwave besides passive and active thermoregulation processes, such as instantaneous regional blood flow. All these factors call for sophisticated numerical simulation studies, phantom model tests, careful extrapolation of experimental results, and further epidemiological investigations (4, 6). In the present study, the possible effects of pulsed, low-frequency-modulated EMF emitted by mobile phones on the autonomic modulation of the heart rate was investigated by HRV and HRA analysis at different paced breathing patterns. We supposed altered breathing–respiratory sinus arrhythmia coupling on exposure which could be detected by HRA analysis. We did not find statistically significant differences in HRA parameters due to exposure, and even the earlier published (26) reporting on 1:1 vs. 1:2 breathing did not reach statistically significant differences after Holm–Bonferroni correction of the significance level in spite of the similar study population and pace-adhered breathing verified by recordings. The powerful correlation of Guzik and Porta indices observed in this study also calls attention to their redundancy, since assuming a constant heart rate during the examination period, the number of heart beats and the inter-beat differences in the accelerating and decelerating set of points of the Poincaré plot are reciprocal, as detailed in our previous investigation (26).
On the other hand, we detected statistically significant elevation of the SDNN, CV% and RMSSD parameters during irradiation in symmetrical but not in asymmetrical breathing patterns. Thus 1:1 paced breathing seems to be a sensitizing factor in HRV analysis, at least concerning the effects of EMF. The higher inter-individual variation of HRV in the 1:1i group (Figure 1B-D) may be the result of the presence of some individuals with electromagnetic hypersensitivity, or due to the non-uniform power of exposure since actual emissions during the experiment were not measured; this is a limitation of the study. The power output of mobile phones is automatically set to the necessary minimum (31) in order to increase battery life and reduce emission. Individuals with electromagnetic hypersensitivity only comprise a minority of the population: they report symptoms that may be related to exposition to weak EMF (32). However, Hietanen and co-workers did not ascribe subjective complaints of adverse effects of mobile phone emission in a blind study on 20 healthy volunteers (33). In a meta-analysis (34) no evidence was found of a causal association of electrosensitivity and EMF. Eltiti and colleagues conducted a study similar in volume on EMF-sensitive and healthy volunteers: combined GSM900, GSM1800 and UMTS 2000MHz base station-like exposure did not affect subjective well-being, skin conductance or pulse rate (35). Another limitation of our study is the fixed 1:2 and 1:1 sequence of breathing patterns, resulting in longer (altogether 12 minutes) exposure or a pre-exposed state at 1:1 breathing, however, we considered the negative results or only a rapidly recovering transient effect of previous authors (10, 18, 23, 35, 36) regarding cellular phone exposition and HRV, hence we focused on the immediate effects of exposure on HRV and HRA parameters. However, randomization made it possible to examine post-exposure effects in 11 volunteers, where we found opposite changes of RMSSD with 1:2 versus 1:1 breathing patterns compared to the exposed (the same group n=11) or pre-exposed (other group n=9) cohort (Figure 6). This can reflect vagal withdrawal following 12-minute-long mobile phone EMF irradiation.
HRV analysis is a relatively frequently used non-invasive method to analyse the acute impact of mobile phone irradiation on the autonomous modulation of cardiac functions. In our previous study (23), we did not find any changes in standard HRV parameters during short-term 900 MHz GSM EMF exposure in healthy volunteers. Parazzini et al. concluded that short-term RF irradiation of healthy young persons by mobile phone does not result in statistically significant changes but suggested a weak interaction of time domain parameters and LF power (36). Other researchers came to similar conclusions regarding blood pressure and spectral HRV parameters in 120 volunteers (18). Upon habitual exposure to mobile phone radiation, the time domain HRV parameters were lower and the LF/HF ratio was higher, suggesting a shift to sympathetic activation compared to controls not using cellular phone at all (37).
The mean RR-interval at the inspiration-only controlled first run was significantly shorter (corresponding to higher pulse rate) compared to 1:2 paced sham-exposed measurements. The latter was the second stage only in about the half of the participants considering the double-blind randomized exposure. The 15 minutes of orthostatic adaptation prior to recording excludes orthostatic effects; a stress reaction caused by the beginning of the recording or starting the paced breathing might be one explanation.
Szmigielski and co-workers detected altered diurnal rhythm of blood pressure and heart rate in people regularly working exposed to EMF compared to not exposed controls with the same working pattern (16). The effect was additionally proportional to the electric field strength, ranging from 20 to 550 V/m at 0.738-1.503 MHz. Braune and colleagues used a remote-controlled 900 MHz GSM phone at the right hemisphere for 35 minutes in 10 healthy volunteers and found an increase in capillary vasoconstriction due to increased efferent sympathetic activity resulting in higher systolic and diastolic blood pressures but there were no changes in standard baroreflex tests (17). On the other hand, Barker et al. excluded any significant effect in mean arterial pressure and heart rate variability (HRV) parameters on short-term exposure by simulating GSM and TETRA handsets in 120 healthy subjects (18). Tahvanainen and co-workers did not find any change in arterial blood pressure or heart rate during or after 35 minutes of exposure at 900 or 1800 MHz at maximal allowed antenna power in a double-blind placebo-controlled crossover trial in 32 healthy volunteers (19).
Braune et al. suspected blood pressure elevation and increased efferent sympathetic activity due to cellular phone irradiation central effects since there were no changes in standard baroreflex tests (17). Our experimental setup demands more complex interaction between auditory (beep for pacing), cortical (recognition and differentiation of beep frequency and intentional breathing in and out), and autonomic centres (intentional breathing and reflex responses to intrathoracal and blood pressure changes), all of them localized in the vicinity of the microwave signal source. Valentini and co-workers concluded in their meta-analysis of the neurophysiological effects of mobile phone exposure that not only the cortex but also subcortical structures such as the hypothalamus and thalamus can be affected by pulsed microwave radiation since the alpha and sigma bands of EEG are enhanced during exposure (10). However, in our study we cannot exclude the role of baroreflexes following intrathoracal pressure changes due to breathing either at the receptor level or at processing in the brain stem. In addition, the carotid body is close to a mobile phone during use, so its function may also vary with exposure. We believe an interaction exists between the affected centres resulting in a delay of feedback which can increase HRV under 1:1 breathing and EMF exposure. This might be due to the excitation effect of cellular phone exposure on the brain (10), resulting in increased sympathetic activity in accordance with Braune et al. (17) and Ekici and co-workers (37), considering also the relatively longer inspiratory period during the 1:1 breathing compared to the 1:2 pattern. During inhalation, respiratory sinus acceleration is a result of vagal withdrawal and increased sympathetic flow. In our investigation, HF power (vagal effect) decreased and LF power (vagal and sympathetic interactions) increased upon exposure during 1:2 breathing (Figure 2), however, it did not reach statistical significance by the Friedman test. We suspect that microwave exposure acts mainly during inspiration.
Conclusion
To our knowledge, the present examination was the first investigating HRA and HRV parameters with different paced inspiration/expiration ratios during short-term GSM cellular phone exposure in a double-blind repeated-measures design crossover trial. Significant HRA changes were not found on exposure compared to sham exposure, however, time domain HRV parameters on exposure were significantly different at 1:1 paced but not at 1:2 paced breathing, considering symmetrical breathing as a sensitizing factor. We observed a mild post-exposure effect regarding RMSSD. This finding reflects persisting acute effects of GSM handset emission on the autonomic nervous system, although exploring its real effects on health status and on survival, as well as the elucidation of any underlying physiopathology, call for further higher-volume double-blind studies focused on respiration–circulation coupling. Our study reproduced the HRA changes accordingly to the metronome-controlled inspiration/expiration ratios. The strong correlation of Porta and Guzik indices was also demonstrated in the overall study raising their redundant nature. Symmetrical breathing (1:1 inspiration:expiration ratio) can be used as sensitizing factor or ’stress test’ in other HRV/HRA analysis studies.
Acknowledgements
The Authors express thanks to Levente Nyitrai, former student, for his assistance in data acquisition and analysis. The present scientific contribution is dedicated to the 650th anniversary of the foundation of the University of Pécs, Hungary.
Footnotes
This article is freely accessible online.
Conflicts of Interest
The Authors declare that they have no conflict of interest in regard to this study.
Funding
OriginPro v2017 software was funded by 3/2016 AOK-KA grant from the Medical Faculty of the University of Pécs, Hungary. Balázs Németh was supported by the ÚNKP-17-3-III New National Excellence Program of the Ministry of Human Capacities, Hungary.
Ethical Approval
All procedures performed in the study were in accordance with the ethical standards of the Regional Research Ethics Committee and with the 1964 Helsinki declaration and its later amendments.
- Received April 9, 2018.
- Revision received May 27, 2018.
- Accepted May 29, 2018.
- Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved