- 1Biomedical Engineering & Aerospace Neuroscience, Laboratory of Medical Physics, Faculty of Health Sciences, School of Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
- 2Greek Aerospace Medical Association and Space Research, Thessaloniki, Greece
- 3Aeromedical Center of Thessaloniki, Thessaloniki, Greece
- 4Laboratory of Aerospace and Rehabilitation Applications “Joan Vernikos” Arogi Rehabilitation Center, Thessaloniki, Greece
- 5Thirdage llc, Culpeper, VA, United States
Purpose: Increasing the level of gravity passively on a centrifuge, should be equal to or even more beneficial not only to astronauts living in a microgravity environment but also to patients confined to bed. Gravity therapy (GT) may have beneficial effects on numerous conditions, such as immobility due to neuromuscular disorders, balance disorders, stroke, sports injuries. However, the appropriate configuration for administering the Gz load remains to be determined.
Methods: To address these issues, we studied graded G-loads from 0.5 to 2.0g in 24 young healthy, male and female participants, trained on a short arm human centrifuge (SAHC) combined with mild activity exercise within 40–59% MHR, provided by an onboard bicycle ergometer. Hemodynamic parameters, as cardiac output (CO), stroke volume (SV), mean arterial pressure (MAP), systolic blood pressure (SBP), diastolic blood pressure (DBP), and heart rate (HR) were analyzed, as well as blood gas analysis. A one-way repeated measures ANOVA and pairwise comparisons were conducted with a level of significance p < 0.05.
Results: Significant changes in heart rate variability (HRV) and its spectral components (Class, Fmax, and VHF) were found in all g loads when compared to standing (p < 0.001), except in 1.7 and 2.0g. There were significant changes in CO, cardiac index (CI), and cardiac power (CP) (p < 0.001), and in MAP (p = 0.003) at different artificial gravity (AG) levels. Dose-response curves were determined based on statistically significant changes in cardiovascular parameters, as well as in identifying the optimal G level for training, as well as the optimal G level for training. There were statistically significant gender differences in Cardiac Output/CO (p = 0.002) and Cardiac Power/CP (p = 0.016) during the AG training as compared to standing. More specifically, these cardiovascular parameters were significantly higher for male than female participants. Also, there was a statistically significant (p = 0.022) gender by experimental condition interaction, since the high-frequency parameter of the heart rate variability was attenuated during AG training as compared to standing but only for the female participants (p = 0.004).
Conclusion: The comprehensive cardiovascular evaluation of the response to a range of graded AG loads, as compared to standing, in male and female subjects provides the dose-response framework that enables us to explore and validate the usefulness of the centrifuge as a medical device. It further allows its use in precisely selecting personalized gravity therapy (GT) as needed for treatment or rehabilitation of individuals confined to bed.
Introduction
Life evolved and adapted within the gravitational field of Earth and gravity is critical to our existence, since we have to rely on Earth’s gravity, as a fundamental reference, in structuring musculoskeletal support and in the organizing of bodily fluids. Along the years, we evolved specialized motion-sensing receptors in our inner ears acting like biological guidance systems. Also, the cardiovascular system has evolved to operate in Earth’s gravity while standing, sitting or lying down. Years of space research in weightlessness and analog environments led to the conclusion that physiological health and performance on Earth is significantly dependent on gravity (Vernikos et al., 1996; Vernikos, 1996, 1997, 2008, 2017; Blaber et al., 2013b; Tanaka et al., 2016). When gravity is not used efficiently (due to lifestyle, inactivity or bedrest) to maintain the level of health that is appropriate to living in Earth’s 1G environment, it leads to negative consequences in all human physiological systems (Sandler and Vernikos, 1986; Rittweger and Frings-Meuthen, 2013; McDonnell et al., 2019). Orthostatic intolerance (OI), an inability to regulate blood pressure (BP) on assuming an upright stance (Robertson, 1999; Lambert and Lambert, 2014), is commonly experienced by astronauts upon their return to Earth after both short and long duration spaceflight (Buckey et al., 1996; Vernikos, 1996; Lee et al., 2015), as well as in people who resume posture after being confined to bed for experimental, clinical or pathological reasons (Goswami, 2017; Goswami et al., 2017).
We now know from Space research, how human physiology is influenced by the lack of gravity loads (Gz) either acutely or over extended periods of time, with detrimental consequences on all organ systems of the body (Goswami, 2017; Goswami et al., 2017). Sitting and inactivity on Earth are also associated with an elevated risk of mortality from all causes, including cardiovascular disease (Katzmarzyk et al., 2009).
In the upright posture, gravity acts along the head to feet column and pulls fluids from the head to the feet along the Gz axis, leading to a BP decrease at brain and heart level and a significant increase at foot level (Furst, 2020). Removal of hydrostatic gradients, occurring in microgravity conditions and analog environments like bedrest and immobilization, or in physically impaired individuals, triggers adaptations in human physiology to compensate for these changes (Vernikos et al., 1996; Rittweger and Frings-Meuthen, 2013; McDonnell et al., 2019). In the supine position, gravity is applied across the short, chest-to-back fluid column (Gx) with equal BP at the levels of brain, heart and feet.
Standing up after bedrest, challenges the cardiovascular system due to the changes in fluid pressure and redistribution of fluid volume to the feet, resulting in involuntary reflex adjustments needed to get the blood back up to the head (Convertino et al., 1990). Such studies are of fundamental importance to physiologists, helping them understand how the human body adapts to both microgravity and bed rest as well as identify the underlying mechanisms involved in standing up after bedrest without fainting.
Space physiologists initially focused on adaptation processes to develop strategies (countermeasures) to mitigate the detrimental effects in human physiology resulting from the lack of gravity in space. Countermeasures tested included nutrition, vibration, or electrical stimulation. Though artificial gravity (AG) was also proposed, it was never tested systematically in space due to logistical and cost issues. (Goswami et al., 2008; Blaber et al., 2013a; Hargens et al., 2013). Studies using bedrest (BR) or head down tilt bedrest (HDTBR) on the ground to induce microgravity-related changes (LeBlanc et al., 1988; Jeong et al., 2013; Norsk, 2014) also indicated that these countermeasures were only partially effective. It should be noted that these countermeasures were administered once a day (Evans et al., 2004; Stenger et al., 2007, 2012; Hargens et al., 2013). However, when administered multiple times throughout the day, either standing alone or upright exercise, proved to be the most effective countermeasure (Vernikos et al., 1996). Based on these observations, the intermittent administration of Gz with or without exercise using a short-arm human centrifuge (SAHC) was proposed as the most promising countermeasure for mitigating the physiological multisystem deconditioning in space or bed rest (Vernikos, 1997; Evans et al., 2004; Frett et al., 2014; Clément et al., 2015).
The absence of hydrostatic pressure in 0g conditions, causes fluid shifts among compliance vessels and changes in the cardiovascular system, in flow resistance and in diastolic cardiac performance (Peterson et al., 2002), eliciting compensatory neurohormonal responses, such as those from carotid and atrial mechanoreceptors, vascular smooth muscle, stretch receptors, and several hormones. To counteract the effect of gravity on Earth, humans have developed neuronal and humoral adjustments in order to maintain a steady mean arterial pressure (MAP) in case of sudden change in posture or sudden acceleration. We considered that increasing acceleration and Gz while supine to mimic standing conditions would be beneficial not only to astronauts in space but also to bedrest patients or patients with immobility due to pathological causes such as neuromuscular diseases, osteoporosis, OI, geriatric individuals, and vestibular disorders. Daily posture change such as intermittent standing has been shown to counteract short term deconditioning (Convertino et al., 1990; Duvivier et al., 2018). There is a system-specific benefit of daily standing up in Earth’s 1G and there is a preventive value of exposure to daily short-period passive +1Gz vertical acceleration, or activity in +1Gz (Vernikos et al., 1996). Intermittent AG may effectively “train” cardiovascular or neuromotor mechanisms responsible for maintaining orthostatic integrity. However, the appropriate configuration including optimal g-level, angular velocity, gravity gradient, and exercise protocol have not been systematically established.
Goswami et al. (2015) used two G loads, – 1 or 2g at foot level – and monitored continuous BP and heart rate (HR). They concluded that the response of the cardiovascular and cerebrovascular systems during centrifugation with 2-g at the feet, was analogous to the orthostatic challenge exerted by standing up in natural gravity (1Gz).
In order to evaluate optimal G load, compared to standing, the present study established dose-response curves of significant cardiovascular indices. Hemodynamic parameters along the sequential levels of 0.5, 0.7, 1, 1.2, 1.5, 1.7, and 2.0g were integrated to construct dose response curves in order to identify the optimal G load that would be equivalent to standing and beneficial for rehabilitation and training. We simultaneously recorded non-invasively cardiac output (CO), stroke volume (SV), MAP, systolic blood pressure (SBP), diastolic blood pressure (DBP) and heart rate (HR), as well as blood gas analysis measurements of partial pressure of oxygen (pO2), partial pressure of carbon dioxide (pCO2), oxygen content (O–), carbon dioxide content (CO2–) and oxygen saturation (O2sat%). We specifically selected CO, SV, and MAP, since studies of extended bed rest in humans have revealed significant reduction in cardiac SV and CO when activity and posture change were restricted (Saltin et al., 1968). This data is designed to serve as a guide in identifying AG rehabilitation regimes for those following inactivity or bedrest, various mobility disorders or to mitigate detrimental effects of microgravity in astronauts.
This study proposes a robust analysis framework. More specifically, it includes a variety of G loads (seven separate G loads) as well as a large number of outcome measures. The combination of these two novelty factors will facilitate the use of AG over a range of pathologically induced immobility conditions.
Materials and Methods
Participants
Twenty-eight healthy volunteers (15 males and 13 females), aged 20–45 years. were recruited by the Aeromedical Center of Thessaloniki, Greece, to participate in this study to evaluate the physiological responses to AG. All participants were physically active, healthy, non-smokers, with no medical history of chronic disease. None of the participants had any history of cardiovascular, metabolic, or neurological diseases, and None of them were professional or elite athletes. Female participants had to test negative on a pregnancy test before enrollment in the study. There were non-physiological cardiovascular features for two participants and missing values for another two participants. So, they excluded from further analysis, resulting thus in twenty-four participants. Their demographics characteristics are reported in Table 1.
Exposure to whole body centrifugation is demanding and the number of enrolled participants is often limited. Therefore, we used a power calculation using G-Power 3.1 (Faul et al., 2013) for an error probability (α) of 0.05, a power (1 − β) of 0.80 and an average effect size (d) of 0.5 (Figure 1A). In order to achieve statistical significance for AG levels, this analysis yielded an estimated total sample size of 28 participants in the case of non-parametric analysis (Wilcoxon signed-rank test) and 27 for parametric testing (paired t-test). Our study finally enrolled 24 healthy participants (12 females). This population is compatible with previous studies in which cognitive function and hemodynamic parameters were assessed during varying gravitational loading (Stenger et al., 2007; Duda et al., 2012; Masatli et al., 2018; Winter et al., 2018; White et al., 2019; Laing et al., 2020) (N = 12–28).
Figure 1. (A) Visualization of the analysis performed for estimating the proper sample size for the proposed statistical analysis. The power analysis indicated that the proper sample size differed between parametric (28) and nonparametric tests (27). (B) Visualization of the short-arm human centrifuge with two rotating beds. Each bed is equipped with either aerobic or resistance training. There is also the feasibility of artificial gravity (AG) training with no exercise at all in the case of patients suffering from severe mobility difficulties. Cardiovascular biomarkers are monitored in real time, providing thus indices of the exercise intensity and alerts in case of critical situations. Multi-modal neurophysiological analysis during the exercise is being performed through a portable Nihon Kohden EEG device with 32 electrodes mapping electroencephalography (EEG), electrocardiography (ECG), electromyography (EMG) and electro-oculogrammic (EOG) activity. (C) Schematic of the +Gz intermittent training protocol. Each phase is labeled accordingly, indicating the corresponding +Gz load over time in minutes.
Participant Selection Criteria
Inclusion criteria for participation were the following: 18–45 years old, absence of neurological or psychiatric disorder, absence of vertigo, nausea or chronic pain, no history of syncope, adequate level of physical fitness for the ability to work out aerobically for half an hour. Participants with a height greater than 2 m, elite athletes, with chronic use of substances (drugs or alcohol), with recent (within 6 months) surgery, current arrhythmias, severe migraines, pregnancy, epilepsy, cholelithiasis or kidney stones, dehydration, recent wounds from surgery, recent fractures, acute inflammation or pain, newly inserted metal pins or plates and newly implanted stents were excluded from the study.
After extensive baseline physiological data collection and passing a comprehensive medical examination, they joined the Laboratory of Aerospace Applications and Rehabilitation “Joan Vernikos” at AROGI Rehabilitation Center that houses the SAHC, where the study was conducted.
Participants were fully briefed about the protocol, the sensations they might experience and possible adverse consequences of the acceleration (nausea, dizziness, etc.). Prior to the first study day, each participant was again familiarized with all aspects of the study, including the equipment and personnel involved. They underwent an orientation session on the centrifuge before the start of the study.
At any time, participants had the option to press a “panic button” if they wished to abort the test. Direct verbal contact was maintained between the study physician and each participant. All participants gave their written consent after being informed of any potential risks. They were also asked to complete participation questionnaires for eligibility criteria. The written consent forms were securely stored at the Aeromedical Center. The study was reviewed and approved by the Bioethical Committee of the School of Medicine of the Aristotle University of Thessaloniki (179/19.03.2020).
Study Protocol
An intermittent centrifugation protocol was selected to assess cardiovascular responses including maintenance of baroreflex and parasympathetic activity (Clément and Bukley, 2007; Clément et al., 2015). The study protocol was registered in the ClinicalTrials.gov (Identifier: NCT04369976).
Participants were advised to abstain from caffeine during the test day, avoid alcohol, and any medications during the preceding 12 h, avoid eating 2 h before the test and avoid exercising heavily 24 h before and after the test.
Participants joined the lab at a prearranged time, approximately the same time every study day – around noon for all – and biological material (blood, urine, and saliva) was taken before and after the end of each study. The participants were connected to non-invasive monitors. Clear Sight monitor1 for hemodynamic measurements (CO, MAP, SV, SBP, DBP, and HR), Cnoga Tensor Tip MTX for blood gas analysis (pO2, pCO2, PH, O –, CO2–, and O2sat%) and Nihon device for ECG and HRV measurements in different body positions such as standing, supine, 0.5, 0.7, 1, 1.2, 1.5, 1.7, and 2g (Figure 1C – schematic of the Gz intermittent training protocol).
We combined centrifugation on the SAHC with mild intensity exercise based on % maximum heart rate/maximum heart rate (MHR), due to recent studies suggesting that intermittent centrifugation combined with exercise has greater beneficial effects on the cardiovascular system, avoiding OI, and decreasing the possibility of pre-syncope during acceleration due to blood flow rushing towards the feet and contracting calf muscles to increase venous return (Yang et al., 2007; Bonjour et al., 2010; Yang et al., 2010; Diaz-Artiles et al., 2018). The symptoms considered as OI included nausea, sweating, confusion, dizziness, a narrowing of the visual field or the presence of hemodynamic decompensation criteria, such as new-onset ECG abnormalities, an abrupt drop in MAP of >20 mm Hg or a critical narrowing of the pulse pressure.
We adhered to international standards to ensure the safe installation and operation of the device in accordance with the international Directive 2006/42 EC and National Directive FEK 97/A/25-06-2010. Declaration of conformity (CE) marking was affixed in 2019, showing compliance with all relevant Directives and the Essential Health and Safety Requirements (EHSR) detailed within them. Additionally, the device was approved by the Calibration and Quality Services in 2020 (CQS) with calibration certificate number 399 × 0920. The construction of the device was funded by the Aeromedical Center of Thessaloniki, Greece. A safety harness was placed on the participants’ body and they were strapped with their head close to the center of the 2 m-radius centrifuge.
All signals were verified for quality and integrity, and safety procedures were reviewed. The participants’ eyes were covered with an eye cover and lights were off to remove visual cues during rotation. They were asked to avoid moving their heads during centrifugation to avoid Coriolis effect.
Participants remained in the standing position for 5 min. They were then asked to lie down on the centrifuge for 5 min for baseline data collection. They were exposed for a total of 35-min centrifugation intermittently administered for 5 min at each g load, followed by a 6 min “wash-out period” (pause). Gz included 0.5g (approximately 16 rpm), 0.7g (approximately 19 rpm), 1g (approximately 23 rpm), 1.2g (approximately 25 rpm), 1.5g (approximately 28 rpm), 1.7g (approximately 30 rpm), and 2.0g (approximately 32 rpm). The rpms were adjusted according to the participant’s height. The participants spun relatively fast and then stayed at the desired G level for 5 min. During each G load the mean value of the last 60 s was considered for evaluation. There were 7 × 5-min sessions of centrifugation and 7 × 6-min wash out periods. The total time of the procedure, including 5 min sitting and 5 min lying down and time to fasten seat belts, was 1 h and 40 min. After the 1.7 and 2.0g phase, there was a gradual reduction back to +1Gz and to level zero Gz.
For the hemodynamic measurements and monitoring, non-invasive medical devices were used. The Clear Sight beat to beat non-invasive cardiovascular monitor was attached to each participant’s finger for continuous measurements. The participants were asked to hold their hand at heart level and not move them for the entire centrifugation protocol due to the height reference that those devices have, since they are supposed to be connected at heart level to counteract of hydrostatic changes that occur during centrifugation. If signs of OI were present, the +Gz level was reduced or the protocol was terminated. Fortunately, no such phenomenon occurred and all participants were able to successfully complete the experimental protocol. Simultaneously, another noninvasive device, a Cnoga Tensor Tip MTX device was attached to the subject’s finger for blood gas analysis as well as oxygen saturation measurements.
The purpose of implementing a graded +Gz protocol was to reduce the occurrence of OI due to rapid inductions of high loads. Therefore, the first Gz loads from 0.5 to 1g (mild hyper-gravity) were used to reduce any sudden occurrence of an OI event and to prime the baroreflex system for a moderate gravitational stress stimulus. The g loads of 1.5, 1.7, and 2g induced the desired cardiovascular reactions. The wash-out period allowed for a return to baseline conditions before starting the following centrifugation round.
Systolic blood pressure mmHg, DBP mmHg, MAP mmHg, CO L/min, SV mmHg, and heart rate (HR), as well as partial pressure of oxygen (pO2) mmHg, partial pressure of carbon dioxide (pCO2) mmHg, oxygen content (O−−), carbon dioxide content (CO2−−), and oxygen saturation (O2sat%) were monitored continuously and non-invasively with finger plethysmography as the Clear sight (Edwards Life Sciences), the Cnoga Medical Tensor Tip (MTX device) and the Nihon device for continuous EEG/ECG monitoring, but only ECG is considered here. Total peripheral resistance (TPR) was calculated by dividing MAP/CO. Cardiac power (CP) was calculated as {MAP × CO}/451, where MAP = [(systolic blood pressure - diastolic blood pressure)/3] + diastolic blood pressure. To control for significant baseline anthropometric gender differences, CO, and SV were converted to cardiac index (CI L/min/m2) and stroke volume index (SVI ml/m2).
Additionally, HRV analysis was performed. The pre-processing for HRV consisted of normalization of the ECG time-series in the range of 0–1 in order to facilitate the robust detection of R peaks. The R peak detection was performed through a custom-based script in MATLAB, which estimated peaks greater than a given threshold (≥0.9), while the minimum distance between two consecutive R peaks was set at 0.6 s. The algorithm was previously validated in a semi-automatic sleep detection algorithm (Chriskos et al., 2018).
For the HRV Analysis, since the R–R timeseries is not uniformly sampled across time and does not share the same sampling frequency as the ECG, traditional methods like the Fourier transform are not robust in describing the signal properties. Therefore, we selected to use the Lomb–Scargle periodogram (LSP) in order to extract spectral HRV features (Ruf, 1999). This analysis resulted in the estimation of the relative energy contribution of the following frequency bands: (a) ultra low frequency (ULF) < 0.0033 Hz, (b) very low frequency (VLF) in the range 0.0033–0.04 Hz, (c) low frequency (LF) in the range 0.04–0.15 Hz, (d) high frequency (HF) in the range 0.15–0.4 Hz, and (e) very high frequency (VHF) > 0.4 Hz. The algorithm also estimates the frequency with the maximum power as well as the class in which it belongs. The latter is a scalar number in the range 1–5 (ULF: 1, VLF: 2, LF: 3, HF: 4, and VHF: 5). We also estimated the ratio LF/HF.
Statistical Analysis
The last 60 s from each phase were averaged and used for the statistical analysis of the hemodynamic variables, in order to highlight the maximal cardiovascular response (Habazetti et al., 2016), and as steady state without the interference of acute changes in +Gz transitions and with minimum variability (Diaz Artiles et al., 2016).
The goal of the analysis was to record a pattern of each cardiovascular variable in different AG loads and more specifically to check the correlations of those variables among standing body position and different g loads exerted at the feet of the individual.
Statistical tests were performed using SPSS (IBM SPSS Statistics 25). In this analysis, the 8 AG levels (standing_5 min, 0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g) were merged to five [standing_5 min, ultra-low intensity (0.5g), low intensity (0.7 and 1.0g), medium intensity (1.2 and 1.5g), and high intensity (1.7 and 2.0g) or to 3 (standing_5 min, 1.7 and 2.0g).
The five AG intensity levels were tested for normality using the Kolmogorov–Smirnov test and for homoscedasticity using a Levene’s test. A one-way repeated-measures ANOVA was conducted with AG levels (5 or 3) as the within-subjects factor on biosignal features which served as the dependent variables. In addition, pairwise comparisons were calculated using the Bonferroni post hoc procedure. In all cases, significance was taken at the a = 0.05 level. Biosignal features are presented as the mean ± standard deviation (Table 2).
Table 2. Description of the biomarkers that showed a statistically significant main effect of gender when comparing the AG training levels with the standing condition.
When resulting data violated an assumption of one-way repeated-measures ANOVA, a non-parametric Friedman’s ANOVA was used for comparison of AG levels (5 or 3) on bio signal features. Furthermore, non-parametric post hoc procedures conducted using Wilcoxon signed-rank tests, but correcting for the number of comparisons (the Bonferroni correction).
For HRV statistical tests were performed using SPSS (IBM SPSS Statistics 25). Spectral properties (Class, Fmax, VLF, LF, HF, VH, and LF/HF) at each AG intensity level (standing_5 min, 0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g) were tested for normality using the Kolmogorov–Smirnov test and for homoscedasticity using a Levene’s test. A one-way repeated-measures ANOVA was conducted with AG levels as the within-subjects factor on spectral properties which served as the dependent variables. In addition, pairwise comparisons were calculated using the Bonferroni post hoc procedure. In all cases, significance was taken at the a = 0.05 level. When resulting data violated an assumption of one-way repeated-measures ANOVA, a non-parametric Friedman’s ANOVA was used for comparison of AG levels on spectral properties. Furthermore, non-parametric post hoc procedures conducted using Wilcoxon signed-rank tests, but correcting for the number of comparisons (the Bonferroni correction).
Additionally, correlation analysis was performed and significant correlations were found for MAP and CO. The data collected from the participants in nine time-instances (standing pre, supine, 0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g) were checked to determine whether they followed a normal distribution. In some cases, this hypothesis was rejected. Therefore, Spearman correlation was performed in all datasets through the Python 2.7 environment with the use of Pandas library, since it is ideal for the evaluation of monotonic relationships between values, whether they are linear or not (Lehman et al., 2005). The correlation matrix was constructed, where ρ values of relationships with no statistical significance (p > 0.05) were set equal to 0. The main diagonal matrix the matrix was set to 0 as well.
Short Arm Human Centrifuge
The short-arm human centrifuge (SAHC) was proposed as a promising countermeasure and training tool within the context of space exploration (Evans et al., 2004; Frett et al., 2014; Clément et al., 2015). The new constructed SAHC (Patent 1009812/13/09/2019, EP3795131), (Figure 1B) is located at the Laboratory of Aerospace and Rehabilitation Applications “Joan Vernikos” at the premises of AROGI Rehabilitation Center, with a diameter of 4 m, with 4-bed capacity and is powered by electricity with a capability of a maximum +3.5 Gz-load at the feet, equipped with a detachable cycle ergometer and resistance pedals with elastic bands (Figure 1B). It was constructed according to both national and international regulations and appropriate certification for function and safety, suitable for patients with restricted mobility disabilities and the capability of gathering physiological data during centrifugation.
Before the study initiation, the height and weight of each participant were measured and more importantly, so was the distance from the central axis to the feet, in order to calculate the g-load at the feet. The head of the participant was placed close to the center of centrifuge rotation, while the feet were directed outwards from the center of rotation in order to allow linear distribution of the g-load in a head-to-toe axis. The g-load at the feet depended on the distance from the axis and the rotational speed of the centrifuge. The g load at the head level was close to zero. Although this does not simulate the typical hydrostatic difference created by standing up in natural gravity, it does generate similar physical stressors in the direction of the feet. Participants successfully completed activity at a voluntary level based on 40–60% of the MHR, provided by an onboard bicycle Ergometer during centrifugation, which they tolerated well, including during the spin-up and spin-down phases of each run.
Exercise Activity
Although the activity was at a voluntary level provided by the onboard bicycle ergometer, it was controlled and directed not to exceed or decrease the pace according to their personalized target heart rate zone, calculated based on both % MHR and Karvonen formula corresponding to 40–59% of MHR and 30–49% Karvonen corresponding to 5–8 METS (DeVries and Housh, 1994, p. 297; Wilmore and Costill, 1994, p. 524; Powers and Howley, 1997, p. 292). Karvonen formula was reasonably accurate method for estimating exercise intensity as demonstrated by Davis and Convertino (1975). Both methods seemed to be similar. The subjects were monitored real-time via a Bluetooth Polar H10 device, which is an accurate heart rate sensor, and the intensity of the exercise was kept at light intensity pace. The subject was followed during the whole process and was advised accordingly to keep exercise pace within the predetermined limits. The light intensity exercise was verified by HRV recorded on an ECG device placed on the subject during centrifugation. The verification procedure employed the first step of the HRV analysis described in the latter part of the above-mentioned section “Study Protocol.” More specifically, the identification of the R peaks resulted in the normalized heart rate which computed as the average number of R peaks per minute. Then the ratio of the heart rate divided by (220-age) was used to estimate whether the exercise training could be regarded as light intensity.
Results
Motion Sickness and Comfort Results
No soreness other than normal exercise fatigue was reported. None of the participants reported an overall motion sickness rating higher than 1 in a 0–10 scale and only one reported slight symptom (scale 1) due to momentarily moving his head to the opposite side. The spin-down process caused transient lightheadedness, due to the sensation of moving counter clock-wise during slowing down to stop spinning (Duda et al., 2012), but lasted only several seconds After that, the participants were allowed to remove their eye cover for the rest of the waiting time and during the wash-out periods.
Gender Analysis
A 2 (gender: male vs female) × 2 (condition: standing vs AG training) between subjects’ ANOVA was conducted to study cardiovascular variables differences (CO, CI, CP, MAP, SBP, DBP, HR, Class, Fmax, VHF, VLF, LF, HF, and LF/HF) between male and female participants. There was a statistically significant main effect of gender in CO and CP biomarkers (see Table 2). A statistically significant gender × condition interaction was observed only in the HF biomarker (see Table 3). The AG training was set as the mean biomarker value for the seven (0.5, 0.7, 1.0, 1.2, 1.5, 1.7, 2.0g) AG training loads.
Table 3. Description of the biomarkers that showed a statistically significant gender by experimental condition interaction when comparing the AG training levels with the standing condition.
In order to further investigate this interaction, we conducted Bonferroni post hoc comparisons. It indicated a statistically significant effect only for females (p = 0.004). More specifically, females showed decreased HF values during the AG training in comparison with the standing condition (see Table 3 for more details).
We further investigated which specific AG task was responsible for this decrease. So, we performed non-parametric analysis for the female participants comparing the standing condition with each AG load separately. The Friedman test showed a marginally significant effect of the AG level on HF on female participants (χ2(7) = 12.410, p = 0.088). Wilcoxon tests were used to follow-up these findings.
Since, we performed multiple comparisons among the standing condition and the seven AG loads, a Bonferroni correction was applied and so all effects are reported at a p = 0.007 (0.05/7) level of significance. It appeared that there were only marginal statistical differences, since none of them survived the multiple comparisons’ correction: (a) from standing to 0.5g (Z = −2.490, p = 0.013), (b) from standing to 1.0g (Z = −2.223, p = 0.026), (c) from standing to 1.7g (Z = −2.401, p = 0.016), and (d) from standing to 2.0g (Z = −2.223, p = 0.026).
Cardiovascular and Blood Gas Analysis
Cardiovascular and blood gas analysis variables were examined in graded g loads (0.5, 0.7, 1.0, 1.2, 1.5, 1.7 and 2.0g) and mean values and standard deviations of all measured values are depicted in Table 4.
Table 4. Mean (M) ± Standard Deviations (SD) of 15 biosignal features for the five AG levels, where the high intensity level due to its significance compared to standing, it was presented separately as 1.7 and 2.0g.
A one-way repeated-measures ANOVA was conducted with five intensity levels (standing 5 min, ultra-low intensity, low intensity, medium intensity, and high intensity) as the within-subjects factor on HR, pO2, pCO2, and which served as the dependent variables. Assumptions of normality using the Kolmogorov–Smirnov test and homogeneity of variance for homoscedasticity using a Levene’s test were met (see Figure 2).
Figure 2. HR for each of the 24 participants in the standing position (red dotted line) and at each artificial gravity (AG) load (denoted with a diamond of a different color: blue, 0.5g; green, 0.7g; red, 1.0g; magenta, 1.2g; cyan, 1.5g; yellow, 1.7g; black, 2.0g).
Heart Rate Variability
Investigating the temporal variation between each heartbeat, HRV revealed a statistically significant main effect of AG levels for Class , Fmax and VHF (Table 5). Performing Bonferroni post hoc procedure, pairwise comparisons indicated a trend to significant mean difference for Class between:
(i) Standing 5 min and 1.5g (mean difference = −0.378, p = 0.035, CI (95%) −0.740, −0.015) and
(ii) Fmax between standing 5 min and 0.5g (mean difference = −0.118, p = 0.026, CI (95%) −0.228, −0.008),
(iii) Standing 5 min and 0.7g
(iv) Standing 5 min and 1.2g
(v) Standing 5 min and 1.5g
(vi) No significant difference was revealed between standing and 1.7 and 2.0g for class and Fmax.
Table 5. Visualization of the spectral heart rate variability (HRV) features that yielded statistically significant results due to the artificial gravity (AG) training.
A significant mean difference of AG levels for VHF was revealed between:
(i) Standing 5 min and 0.7g
(ii) A total of 0.7 and 2.0g
suggesting the activation of the parasympathetic system at 0.7g initially and then at 2.0g, indicating the effectiveness of the centrifuge as training device.
For the pairwise comparisons, all the mean differences were significant at the 0.05 level, adjusted for multiple comparisons: Bonferroni. These results and their statistical significance are described in Table 6.
Table 6. Description of the heart-rate related exercise intensity ratio that showed a statistically significant a main effect of condition when comparing the AG training levels between them.
The pairwise comparisons of the maximum HRV frequency (Fmax) demonstrated a statistically significant increase between Standing when compared to 0.5g (p = 0.026), 0.7g (p = 0.005), 1.2g (p = 0.024), and 1.5g (p = 0.011), but no statistically significant difference between 1.7 and 2.0g, indicating that those loads were closer to standing.
Concerning VLF, LF, HF, and LF/HF, as assumptions of one-way repeated-measures ANOVA were violated, a non-parametric Friedman’s ANOVA was used for comparison of the eight AG levels on the abovementioned spectral properties. Furthermore, non-parametric post hoc procedures conducted using Wilcoxon signed-rank tests, but correcting for the number of comparisons by accepting something as significant only if its significance is less than a/number of comparisons (the Bonferroni correction). In this analysis, due to 8 AG levels (standing plus 7 AG levels), we have 28 comparisons, therefore, rather than use 0.05 as our critical level of significance, we had use 0.05/28 = 0.002. The Friedman test showed significant effect of AG levels on LF (χ2(7) = 15.833, p = 0.027) and LF/HF (χ2(7) = 17.045, p = 0.017) and marginal main effect for VLF (χ2(7) = 13.833, p = 0.054), and HF (χ2(7) = 5.909, p = 0.550). These results and their statistical significance are described in the Supplementary Material.
As shown in Figure 3, both standing condition and 2.0g AG induce greater sympathetic (LF) and lesser parasympathetic (HF) response. The VHF ratio is also lower for these two conditions. Wilcoxon tests were used to follow-up these findings. A Bonferroni correction was applied and so all effects are reported at a = 0.002 level of significance. It appeared that: (i) LF significantly increased from 0.7 to 2.0g (Z = −3.376, p = 0.001) and (ii) LF/HF significantly increased from 0.7 to 2.0g (Z = −3.406, p = 0.001), indicating the effectiveness of the centrifuge as training device.
Figure 3. Relative contribution (%) of the most significant HRV features between standing and all AG loads. The low frequency (LF) ranges is denoted with a blue bar (0.04–0.15 Hz), the high frequency (HF) is denoted with a green bar (0.15–0.40 Hz), and the very high frequency is denoted with a green bar (>0.40 Hz).
A repeated measure of ANOVA with a Greenhouse–Geisser correction and condition as within subjects’ factor (0.5, 0.7, 1, 1.2, 1.5, 1.7, and 2g) was conducted to study heart-rate related exercise intensity ratio differences among the aforementioned seven AG levels. There was a statistically significant main effect of condition [F(3.913,120) = 9.331, p < 0.001] (see Figure 4 and Table 6 for more details and post hoc comparisons).
Figure 4. Visualization of the Heart Rate Intensity ratio (%) for each artificial gravity (AG) condition. This is a cardiovascular index of the exercise intensity level. It is estimated as the heart rate (beats per minute) during each training phase divided by the maximum heart rate (MHR) intensity. The MHR for each participant is defined as 220-participant’s age. Ratio values ranging from 40 to 59% are indicative of light exercise.
Blood Gas Analysis
No significant main effect was revealed comparing the different g loads for pH, , pCO2, , and O2sat%. Comparing standing with 1.7 and 2.0g with the Friedman test showed significant effect of AG load on pCO2 (χ2(2) = 7.342, p = 0.025) and a marginal trend toward a significant main effect on ()(χ2(2) = 5.948, p = 0.051). Wilcoxon tests were used to follow-up these findings. A Bonferroni correction was applied and so all effects are reported at a 0.016 (0.05/3) level of significance. It appeared that () significantly decreased from 1.7 to 2.0g (Z = −2.433, p = 0.015). In addition, if a Bonferroni correction is not applied so the level of significance will be at 0.05, then, we observe that pCO2 significantly decreased: (i) from standing to 2.0g (Z = −2.034, p = 0.042) and ii) from 1.7 to 2.0g (Z = −2.171, p = 0.030). So, by increasing the AG level of 1.7–2.0g we propose that the respiration rate increased and consequently the mean value of pCO2 dropped. Note a substantial decrease of pCO2 at 1.2g where the baroreceptors were triggered.
Cardiac Output, Cardiac Index and Cardiac Power Output
The Friedman test showed significant effect of the AG load of CO (χ2(4) = 19.713, p = 0.001). Wilcoxon tests were used to follow-up these findings. It appeared that:
(i) CO significantly decreased: (a) from standing to ultra-low intensity (Z = −3.651, p = 0.000), (b) from standing to low intensity (Z = −3.346, p = 0.001), (c) from standing to medium intensity (Z = −3.258, p = 0.001) and (d) from standing to high intensity (Z = −2.844, p = 0.004) (Figure 5),
(ii) CI significantly increased from standing to ultra-low intensity (Z = −3.011, p = 0.003) and significantly decreased: (a) from standing to low intensity (Z = −3.285, p = 0.001) and from (b) standing to high intensity (Z = −2.857, p = 0.004),
(iii) Cardiac power significantly decreased: (a) from standing to ultra-low intensity (Z = −3.143, p = 0.002), (b) from standing to low intensity (Z = −3.157, p = 0.002), and (c) DBP significantly increased from low intensity to high intensity (Z = −2.937, p = 0.003).
Figure 5. Cardiac output (CO) for each of the 24 participants in the standing position (green) and at each artificial gravity (AG) load (denoted with a boxplot of a different color: red, 0.5g; blue, 0.7g; yellow, 1.0g; black, 1.2g; pink, 1.5g; magenta, 1.7g; gray, 2.0g).
The Friedman test showed a significant effect of AG level on CO (χ2(2) = 15.432, p = 0.000), CI (χ2(2) = 15.432, p = 0.000), CP (χ2(2) = 8.374, p = 0.015). Wilcoxon tests were used to follow-up these findings. A Bonferroni correction was applied and so all effects are reported at a 0.016 (0.05/3) level of significance. It appeared that:
(i) CO significantly decreased from standing to 1.7g (Z = −3.318, p = 0.001),
(ii) CI significantly decreased from standing to 1.7g (Z = −3.429, p = 0.001),
(iii) CP significantly decreased from standing to 1.7g (Z = −2.571, p = 0.010).
None or only minor changes were seen in the mean SV level during the exercise periods. Upon the end of exercise, a striking but transient increase in CO occurred, resulting from an increase in SV concomitant with a maintained HR (Figure 5). Converting CO to CI (CO/BSA), similar findings were revealed with an increase towards the higher g loads.
Mean Arterial Pressure
The Friedman test showed significant effect of AG level on MAP (χ2(4) = 16.201, p = 0.003). A second Friedman’s ANOVA was used to compare standing with 1.7g and 2.0g and there was no statistically significant difference. Both 1.7 and 2.0g loads of Art G result in a non-significant change in MAP with respect to standing (Figure 6).
Figure 6. Mean arterial pressure for each of the 24 participants in the standing position (red dotted line) and at each artificial gravity (AG) load (denoted with a diamond of a different color: blue, 0.5g; green, 0.7g; red, 1.0g; magenta, 1.2g; cyan, 1.5g; yellow, 1.7g; black, 2.0g).
Systolic Blood Pressure and Diastolic Blood Pressure
The Friedman test showed significant effect of the five AG level on SBP (χ2(4) = 9.595, p = 0.048). Wilcoxon tests were used to follow-up these findings. A Bonferroni correction was applied and so all effects are reported at a 0.005 level of significance. SBP marginally significant decrease from high intensity to low intensity (Z = −2.647, p = 0.008). A second Friedman’s ANOVA test comparing the standing with 1.7 and 2.0g did not show any statistical significance for SBP and DBP, which means that 1.7 and 2.0g results in a non-significant change in BP with respect to standing. SBP followed the MAP and leveled off after 1.7g (Figure 7).
Figure 7. Systolic blood pressure for each of the 24 participants in the standing position (red dotted line) and at each artificial gravity (AG) load (denoted with a diamond of a different color: blue, 0.5g; green, 0.7g; red, 1.0g; magenta, 1.2g; cyan, 1.5g; yellow, 1.7g; black, 2.0g).
Correlation Results
From the correlation matrix (Figure 8) of the hemodynamic parameters of CO and MAP standing significantly correlated with 1.7 and 2.0g. This enhances the significance of those parameters to construct dose response curves and also determine the optimal g load for training similar to standing condition along the Gz axis.
Figure 8. Correlation Matrices for MAP (left picture) and CO (right picture). ρ values with no statistical significance (p > 0.05) were set to zero (blue). Furthermore, the main diagonal of the matrix was set to zero (blue) as well. Dark red in the correlation matrix represents the higher ρs (0.7 for MAP and 0.8 for CO), blue represents the zeros, while the rest of the colors all values between.
Dose Response Curves
Clinical monitoring of the significant cardiovascular variables during centrifugation at graded Gz from 0.5 to 2.0g enabled the construction of dose response curves. These were constructed from measurements of MAP, CO and HR (Figure 9) in all 24 young healthy participants (12 male and 12 females) studied at seven different loads of AG 0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g. Furthermore, the median effective dose at which 50% of subjects exhibit the specified biomarker effect is calculated (ED50), which constitutes a useful guide for determining the appropriate G load for training in conditions that are similar to those of standing along the Gz axis (Figure 9). In addition, ED50 enables the selection of the most useful and safer G load, avoiding undermining morbidity, particularly in those induced by inactivity in bedrest, by various mobility disorders or to mitigate detrimental effects of microgravity in astronauts. Comparing CO, MAP, and HR under seven AG loads (0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g), we found increased levels of HR, CO, MAP, and in high loads of AG. The ED50 was about 1.08g for MAP, 1.13g for CO, and 0.91g for HR, respectively.
Figure 9. Visualization of the dose-response curves for mean aretrial pressure (MAP), cardiac output (CO), and heart rate (HR) in seven AG loads (0.5, 0.7, 1, 1.2, 1.5, 1.7, and 2g). The MAP-ED50 is about 1.08g, corresponding to a mean MAP of 85.09. The CO-ED50 is about 1.13g, corresponding to a mean CO of 6.17. The HR_ED50 is about 0.91g, corresponding to a mean HR of 75.88.
Discussion
The present study provided insight into determining the effective Gz suitable for training or rehabilitation, focusing mainly on the cardiovascular system. In order to come to this conclusion, data from 24 healthy volunteers were collected in seven different loads of AG (0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g). The data included SBP, DBP, MAP, CO, SV, HR, pO2, pCO2, O2–, CO2–, and O2sat%. The main findings involved the detection of g load equivalent to standing for hemodynamic variables and the construction of dose response curves for the most significant ones. The study was designed to begin to systematically build the framework for the study of gravity physiology, its threshold requirements, as well as its potential uses in training or rehabilitation therapy.
Standing on Earth is crucial for the maintenance of all physiological systems. According to the present study, standing equivalent benefits, assumed to be 1Gz, in microgravity analog environments or in the case of immobilized patients, can be achieved by the centripetal acceleration in levels higher than 1.5g in order to maintain cardiovascular function. For some cardiovascular variables, even centrifugation at 2g was not enough to reach comparable standing values, while others needed at least 1.7g.
In this study the spectral features of heart-rate variability analysis were investigated through the LSP, which demonstrated significant effects for a variety of features such as Class, Fmax, VHF, LF, and LF/HF. Based on the study of Migeotte et al. (2009), who underlined the need for further investigation on the changes in graded g load and the gravity gradient in order to find the appropriate g load for SAHC countermeasure, we aimed to examine whether the SAHC training could induce cardiac responses analogous to the standing position. HRV, which has been established during gravity transitions (Sasaki et al., 1999; Seps et al., 2002) was one of the tools used among the cardiovascular indices, to investigate the control and modulation of cardiovascular function through the autonomic nervous system during graded Gz (Malik, 1998). Rowell reported that the changes of HRV in the standing condition are mainly due to the sub-diaphragmatic venous pooling (Rowell, 1993), resulting in a decrease of venous return to the right atrium and the CO in standing and supine position during centrifugation. The higher LF component in the standing position suggest not only higher vagal tone, but also an increment in sympathetic tone, leading to sympathovagal balance decrease due to these changes of vagal activity. The pairwise comparison analysis for most of the features, identified significant differences between the standing position when compared to low and medium intensity level centrifugation (0.5, 0.7, 1.2, and 1.5g). Non-parametric analysis for sympathetic (LF and LF/HF) indices showed significant differences among 2.0 and 0.7g. These results indicate that AG training altered the HRV spectral properties. More specifically, the main effect of Fmax demonstrates that AG training induced a statistically significant shift of the frequency content among the standing and the low and medium AG intensity (0.5/0.7/1.2/1.5g). However, no statistically significant difference was found between standing when compared to the 1.7 and 2.0g. Figure 3 shows that both the standing condition and the 2.0g AG induces greater sympathetic (LF) and lower parasympathetic (HF) contribution. This implies that these specific AG levels are similar to the standing position. Moreover, analysis of the autonomic markers of the sympathetic activity (LF and the ratio of LF/HF) showed that there was a statistically significant increase at 2.0g (p = 0.011), when compared to the 0.7g. This finding implies that the 2.0g induces shorter and more uniform (less varied) R–R intervals due to increased sympathetic activity and simultaneous parasympathetic withdrawal. Also, it demonstrates a similarity with what happens during standing, where BP is regulated by low afferent discharge of the baroreceptor in the aortic arch and carotid sinus, leading to increased HR due to vagal withdrawal and sympathetic nerve activation (Smith et al., 1994; Heesch, 1999; Verma et al., 2018). AG training at 2.0g causes increased physiological stress and elevated fatigue level in comparison to the 0.7g (Makivić et al., 2013). Similarly, there were statistically significant changes among the 0.7g condition when compared to the standing (p = 0.019). Migeotte et al. (2009) reported that the change of HRV in hypergravity are due to venous pooling reduced filling of the right atrium and autonomic adaptation. The major point to be made is that VLF, as the sympathetic nervous system was activated, increased with the gradual increase of g load and was followed by a rapid increase of MAP by constricting arterioles, followed by increased heart rate and SV (Figure 2 and Table 4). Additionally, the higher muscle metabolism causes vasodilatation and therefore peripheral resistance decreases.
In our healthy population cohort, there was a statistically significant increase of VHF relative energy contribution, which denotes activation of the parasympathetic nervous system and the baroreceptor reflex responses to changes in arterial pressure, which resulted from the transition from 0.5 to 0.7g (p = 0.019) CI (95%) and then from 1.7 to 2g (p = 0.012) CI (95%) with a simultaneous decrease in heart rate and CO. These cardiovascular findings are in line with previous studies (Wehrwein and Joyner, 2013) indicating a position change (supine to upright). More specifically, it involves both increased atrial pressure and afferent firing regarding the carotid sinus nerve. These alterations attenuated sympathetic and increased parasympathetic activity. The opposite pattern occurs in cases of decreased arterial pressure.
The HRV findings validate our hypotheses that (a) both standing and 2.0g AG training share common HRV spectral patterns and (b) induce a much greater sympathetic tone than the ultra-low (0.7g) AG intensity. These results extend preliminary knowledge obtained by recent studies which indicated that 2.0g centrifugation shares more similarities with the Standing condition than the 1.0g level (Verma et al., 2018). Moreover, standing seems to induce a shift in the maximum HRV frequency in comparison to all AG levels other than those greater than 1.0g. The higher LF component within the standing position suggest not only higher vagal tone, but also an increment in sympathetic tone. Due to these changes in vagus nerve activity, sympathetic balance will decrease.
Statistical analysis revealed a significant main effect of AG load on the cardiovascular variables of HR, MAP, CO, CI, CP, and SBP. A significant effect of AG intensity level for heart rate was observed compared to standing with all g loads except high intensity load of 1.7 and 2.0g. Furthermore, pairwise comparisons among standing and 1.7 and 2g showed no statistically significant difference between 1.7g and standing, while it was marginally significantly different at 2g. SBP (Figure 7) and DBP dropped in 0.5g and increased significantly at 0.7g and higher g levels. The application of AG introduces a new stress condition that the cardiovascular system needs to overcome in order to assure the appropriate amount of blood flow to all parts of the body. This additional stress is reflected on the BP drop that occurs when the participant is first exposed to AG. Our findings (Table 5) indicate that the increased centrifugation intensity induced greater sympathetic activity through subsequent increases in heart rate and peripheral resistance. The increased g-load raises the need for elevated CO and vasoconstriction in order to facilitate both venous and BP return to normal level at the baroreceptors. This is a physiological compensatory mechanism against high leg pressure and subsequent increased orthostatic challenge. An alternative of this mechanism may be the inclusion of the cycling exercise which acts as a mechanical pump assisting venous return by reducing the blood pooling in the legs and the orthostatic challenge. This was firstly demonstrated a decade ago by pioneering studies (Clément and Bukley, 2007; Kaderka et al., 2010). We adopted a similar design in our study in order to enhance the beneficial role of the proposed intervention. Consequently, CO also increases to meet the new metabolic demands imposed by exercise. The higher AG combined to workload, result in higher CO, which significantly changed at all g intensity levels. Comparing standing with 1.7 and 2g there was a statistically significant difference with 1.7g but not with 2.0g, suggesting that 2.0g load was more comparable to standing (Figure 5). The mean CO value in the standing position (7.21) is higher than any mean value obtained during centrifugation. The lowest CO mean value during centrifugation was obtained at 0.7g (6.07) and the highest at 2.0g (6.46) (Figure 5). There is no significant difference between minimum CO values at any art G load. At 1.7 and 2.0g, maximum CO was restored to that of standing (Figure 5). The same applies to CI and CP.
When comparing MAP, SBP, and DBP between standing, 1.7 and 2g, there were no statistically significant differences, suggesting that MAP may indicate similar effects to standing in both 1.7 and 2.0g loads. From the correlation matrix for CO, the data reveal a strong correlation of standing with 1.7g (0.77) and 0.5g with lying down (0.92). Additionally, covariance was calculated at a positive value. Therefore, CO is a good indirect marker, indicative of a beneficial effect of g load on the human body (Figure 5).
We also observed a main effect of gender on CO and CP and a gender by condition interaction regarding the parasympathetic index of HF of the HRV activity. Although, we should be cautious due to the small sample size and the case of marginally statistically significant differences when correcting the significance level at p = 0.007 to avoid Type-1 errors, these gender differences may be due to the parasympathetic withdrawal response during the AG training when compared to the standing condition. Our hypothesis, which should be further tested by subsequent studies is that due to vagal parasympathetic withdrawal, HR was increased abruptly at the onset of exercise and progressively continued to increase during exercise because of a combination of sympathetic neural and humoral drive, and quickly returned to baseline after cessation of exercise because of withdrawal of neural sympathetic discharge. Women responded differently with increased parasympathetic withdrawal as reported in previous studies (Evans et al., 2001). The difference was attributed to the autonomic function, as well as other factors may contribute, such as myocardial structure, lower SV and hormonal status. The gender difference in orthostatic stress has been demonstrated in other studies as well (Hart and Charkoudian, 2014), with women responding with increased HR while men responded with an increase in vascular resistance (Shoemaker et al., 2001), or increase in HR, DBP, and SVRI (Masatli et al., 2018). Additionally, women presented greater baroreflex sensitivity than men (Drudi and Grenon, 2014). From research in Space, it is known that women are more prone to OI when returning from Space (Wenner et al., 2013) and a need for more effective countermeasure for cardiovascular deconditioning was always a priority for both genders of the human research roadmap (Vernikos et al., 2016) with particular interest in cardiovascular responses of women since this has not been thoroughly documented during intermittent AG exposure so far. Although the number of participants in this study was limited, one of the findings of this study was that the +Gz gradient required in order to elicit significant increases in cardiovascular responses was different and need further investigation to control for significant baseline anthropometric gender differences.
We should also highlight here the wide-range, dose response curves constructed from healthy individuals. These add to the knowledge obtained so far from previous studies where selective g levels of 1.0 and 2.0g were studied only without exercise (Goswami et al., 2015; Masatli et al., 2018; Verma et al., 2018), or dose response curves were generated by regression models since selective Gz of 1 G, and 1.4 G at the feet were used combined with different intensity exercise levels (Diaz-Artiles et al., 2018). To fill in this gap and better understand the relationship between gravitational dose and physiological response (Clément, 2017), dose response curves were constructed based on all intermediate stages tested from 0.5 to 2g in both men and women combined with mild intensity exercise well controlled by HR and HRV. The novelty of this work was the range of gravity levels (0.5, 0.7, 1.0, 1.2, 1.5, 1.7, and 2.0g) studied. The use of standing as the ground reference of Gz on Earth and the comprehensive database of cardiovascular variables (CO, CI, CP, MAP, SBP, DBP, HR, Class, Fmax, VHF, VLF, LF, HF, and LF/HF) analyzed, led to the construction of the dose-response curves. Although the analysis focused on healthy men and women, it is the authors’ goal and future plan to extend the application of the SAHC and continue to build and apply the data base in other age groups of healthy participants, as well as in the rehabilitation of patients suffering from various diseases and conditions, that restrict mobility or confine them to bed.
In summary, the present study was designed initially to introduce and validate a new SAHC device in Greece, to construct a reference framework required to base choices and decisions in research, in identifying optimal G load requirements for patients during AG therapy, for physical therapy and rehabilitation or for use of optimal AG as a spaceflight countermeasures. Mobility is the single crucial element underlying healthy living and healthy aging, as mobility on Earth is gravity-dependent.
Our encouraging findings should be considered in the context of several limitations. First of all, the use of Clearsight monitor to record non-invasive hemodynamics, which offer an estimate measure of CO based on internal models. These devices have a Height Reference that is supposed to be connected at heart level to counteract hydrostatic change (for example, when the participants move their hand where the cuffs are connected). We have decided to not use this correction, as other researchers (Goswami et al., 2015) did before, but instead, we asked all participants to keep their hand at heart level and not move for the entire centrifugation protocol. Another limitation is the Exercise protocol where although we measured light exercise according to % MHR and Carvonen formula and verified by continuous ECG measurement and HRV, we did not measure exactly the force exerted on the ergometer. The number of participants was also limited, especially for female participants, borderline to achieve significance, although sample size was similar to what was used in previous studies that examined mixed sex population during central hypovolemia testing (Evans et al., 2015, 2018; Goswami et al., 2019). Finally, we should also note that the proposed protocol could not reproduce normal standing condition since most people spend about 16 h/day in upright posture.
Conclusion
We introduced a new Short Arm Human Centrifuge in Greece with the first experiment to be conducted, demonstrating good operations and feasibility of scientific experiments, where subjects tolerated well the protocol without any undesired effects. Our main findings approached an answer to the key question of the analogous to standing level of g. We also identified the cardiovascular markers and their dose response curves, providing the answer to the appropriate configuration for individual training and rehabilitation due to pathological conditions which might be proved to be the novel therapeutic physiological approach for people with a broad spectrum of disabilities. Consequently, AG, a spinoff from Space technologies supporting humans in Space, may also support humans who suffer from Gz deprivation on Earth due to lifestyle or medical conditions. Overall, Intermittent AG protocol may effectively “train” cardiovascular hemodynamics or other parameters as neuromotor mechanisms responsible for maintaining orthostatic integrity in young male and female participants with the potential application to space and rehabilitation on Earth.
Data Availability Statement
The cardiovascular features or the raw data supporting theconclusions of this article will be made available by the corresponding author upon request.
Ethics Statement
The studies involving human participants were reviewed and approved by the Bioethics Committee of the School of Medicine of the Aristotle University of Thessaloniki, Greece. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
CK-P, JV, CAF, and PDB conceived and designed the research. CK-P, CN, and CAF conducted the experiments. CAF, CEP., and SG, analyzed the data. CK-P, CAF, CEP, and SG wrote the manuscript. DK, LB, CK-P, JV, and CAF revised the manuscript. PDB and JV supervised the study. All authors contributed to the article and approved the submitted version.
Funding
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. However, parts of this work were supported by the LLM Care (www.llmcare.gr/en) self-funded initiative that emerged as thenot-for-profit business exploitation of the Long-Lasting Memories (LLM Project) originally funded by the ICT-CIP-PSP (Information and Communication Technologies-Competitiveness and Innovation Framework Programme-Policy Support Programme) Program of the European Commission.
Conflict of Interest
JV is the director and founder of the ThirdAge company.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We appreciate the participation of our participants and the contribution of our medical students Aliki Karkala, Anna Nikolaidou, Agisilaos Krachtis, Asimoula Kavada, and Stavros Moschonas.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.644661/full#supplementary-material
Footnotes
References
Blaber, A. P., Hinghofer-Szalkay, H., and Goswami, N. (2013a). Blood Volume Redistribution During Hypovolemia. Aviat. Space Environ. Med. 84, 59–64. doi: 10.3357/ASEM.3424.2013
Blaber, A. P., Zuj, K. A., and Goswami, N. (2013b). Cerebrovascular autoregulation: lessons learned from spaceflight research. Eur. J. Appl. Physiol. 113, 1909–1917. doi: 10.1007/s00421-012-2539-x
Bonjour, J., Capelli, C., Antonutto, G., Calza, S., Tam, E., Linnarsson, D., et al. (2010). Determinants of oxygen consumption during exercise on cycle ergometer: the effects of gravity acceleration respiratory. Physiol. Neurobiol. 171, 128–134. doi: 10.1016/j.resp.2010.02.013
Buckey, J. C., Lane, L. D., Levine, B. D., Watenpaugh, D. E., Wright, S. J., Moore, W. E., et al. (1996). Orthostatic intolerance after spaceflight. J. Appl. Physiol. 81, 7–18.
Chriskos, P., Frantzidis, C. A., Gkivogkli, P. T., Bamidis, P. D., and Kourtidou-Papadeli, C. (2018). Achieving accurate automatic sleep staging on manually pre-processed EEG data through synchronization feature extraction and graph metrics. Front. Hum. Neurosci. 12:110. doi: 10.3389/fnhum.2018.00110
Clément, G. R., Bukley, A. P., and Paloski, W. H. (2015). Artificial gravity as a countermeasure for mitigating physiological deconditioning during longduration space missions. Front. Syst. Neurosci. 9:92.
Convertino, V. A., Doerr, D. F., Eckberg, D. L., Fritsch, J. M., and Vernikos-Danellis, J. (1990). Head-down bed rest impairs vagal baroreflex responses and provokes orthostatic hypotension. J. Appl. Physiol. 68, 1458–1464. doi: 10.1152/jappl.1990.68.4.1458
DeVries, H. A., and Housh, T. (1994). Physiology of Exercise for Physical Education, Athletics, and Exercisescience. Madison, WI: Brown and Benchmark.
Davis, J. A., and Convertino, V. A. (1975). A comparison of heart rate methods for predicting endurance training intensity. Med. Sci. Sports 7, 295–298.
Diaz Artiles, A., Heldt, T., and Young, L. R. (2016). Effects of artificial gravity on the cardio vascular system: computational approach. Acta Astronaut. 126, 395–410. doi: 10.1016/j.actaastro.2016.05.005
Diaz-Artiles, A., Heldt, T., and Young, L. R. (2018). Short-Term cardiovascular response to short-radius centrifugation with and without ergometer exercise. Front. Physiol. 9:1492. doi: 10.3389/fphys.2018.01492
Drudi, L., and Grenon, S. M. (2014). Women’s health in spaceflight. Aviat. Space Enviorn. Med. 85, 645–652. doi: 10.3357/ASEM.3889.2014
Duda, K. R., Jarchow, T., and Young, L. R. (2012). Squat exercise biomechanics during short-radius centrifugation. Aviat. Space Environ. Med. 83, 102–110. doi: 10.3357/asem.2334.2012
Duvivier, B. M. F. M., Bolijn, J. E., Koster, A., Schalkwijk, C. G., Savelberg, H. H., and Schaper, N. C. (2018). Reducing sitting time versus adding exercise: differential effects on biomarkers of endothelial dysfunction and metabolic risk. Sci. Rep. 8:8657. doi: 10.1038/s41598-018-26616-w
Evans, J. M., Ribeiro, L. C., Moore, F. B., Wang, S., Zhang, Q., Kostas, V., et al. (2015). Hypovolemic men and women regulate blood pressure differently following exposure to artificial gravity. Eur. J. Appl. Physiol. 115, 2631–2640. doi: 10.1007/s00421-015-3261-2
Evans, J. M., Knapp, C. F., and Goswami, N. (2018). Artificial gravity as a countermeasure to the cardiovascular deconditioning of spaceflight: gender perspectives. Front. Physiol. 9:716. doi: 10.3389/fphys.2018.00716
Evans, J. M., Stenger, M. B., Moore, F. B., Hinghofer-Szalkay, H., Rössler, A., Patwardhan, A. R., et al. (2004). Centrifuge training increases presyncopal orthostatic tolerance in ambulatory men. Aviat. Space Environ. Med. 75, 850–858.
Evans, J. M., Ziegler, M. G., Patwardhan, A. R., Ott, J. B., Kim, C. S., Leonelli, F. M., et al. (2001). Gender differences in autonomic cardiovascular regulation: spectral, hormonal, and hemodynamic indexes. J. Appl. Physiol. 91, 2611–2618. doi: 10.1152/jappl.2001.91.6.2611
Faul, F., Erdfelder, E., Buchner, A., and Lang, A.-G. (2013). G∗Power Version 3.1.7 [Computer Software]. Kiel: Uiversität Kiel.
Frett, T., Mayrhofer, M., Schwandtner, J., Anken, R., and Petrat, G. (2014). An innovative short arm centrifuge for future studies on the effects of artificial gravity on the human body. Microgravity Sci. Technol. 26, 249–255. doi: 10.1007/s12217-014-9386-9
Furst, B. (ed). (2020). “The effect of gravity and upright posture on circulation,” in The Heart and Circulation, (Cham: Springer), doi: 10.1007/978-3-030-25062-1_24 Available online at: https://link.springer.com/book/10.1007/978-1-4471-5277-4
Goswami, N. (2017). Falls and fall-prevention in older persons: geriatrics meets spaceflight! Front. Physiol. 8:603. doi: 10.3389/fphys.2017.00603
Goswami, N., Blaber, A. P., Hinghofer-Szalkay, H., and Montani, J. P. (2017). Orthostatic intolerance in older persons: etiology and countermeasures. Front. Physiol. 8:803. doi: 10.3389/fphys.2017.00803
Goswami, N., Blaber, AP., Hinghofer-Szalkay, H., and Convertino, VA. (2019). Lower body negative pressure: physiological effects, applications, and implementation. Physiol. Rev. 99, 807–851. doi: 10.1152/physrev.00006.2018
Goswami, N., Evans, J., Schneider, S., von der Wiesche, M., Mulder, E., Rössler, A., et al. (2015). Effects of individualized centrifugation training on orthostatic tolerance in men and women. PLoS One 10:e0125780. doi: 10.1371/journal.pone.0125780
Goswami, N., Loeppky, J. A., and Hinghofer-Szalkay, H. (2008). LBNP: past protocols and technical considerations for experimental design. Aviat. Space Environ.Med. 79, 459–471. doi: 10.3357/ASEM.2161.2008
Habazetti, H., Stahn, A., Nitsche, A., Nordine, M., Pries, A. R., Gunga, H. C., et al. (2016). Microvascular responses to (hyper-)gravitational stress by short-arm human centrifuge: arteriolar vasoconstriction and venous pooling. Eur. J. Appl. Physiol. 116, 57–65. doi: 10.1007/s00421-015-3241-6
Hargens, A. R., Bhattacharya, R., and Schneider, S. M. (2013). Space physiology VI: exercise, artificial gravity, and countermeasure development for prolonged space flight. Eur. J. Appl. Physiol. 113, 2183–2192. doi: 10.1007/s00421-0122523-5
Hart, E. C., and Charkoudian, N. (2014). Sympathetic neural regulation of blood pressure: influences of sex and aging. Physiology 29, 8–15. doi: 10.1152/physiol.00031.2013
Heesch, C. M. (1999). Reflexes that control cardiovascular function. Am. J. Physiol. 277, S234–S243. doi: 10.1152/advances.1999.277.6.S234
Jeong, S. M., Hwang, G. S., Kim, S. O., Levine, B. D., and Zhang, R. (2013). Dynamic cerebral autoregulation after bed rest: effects of volume loading and exercise countermeasures. J. Appl. Physiol. 116, 24–31. doi: 10.1152/japplphysiol.00710.2013
Kaderka, J., Young, L. R., and Paloski, W. H. (2010). A critical benefit analysis of artificial gravity as a microgravity countermeasure. Acta Astronaut. 67, 1090–1102. doi: 10.1016/j.actaastro.2010.06.032
Katzmarzyk, P. T., Church, T. S., Craig, C. L., and Bouchard, C. (2009). Sitting time and mortality from all causes, cardiovascular disease, and cancer. Med. Sci. Sports Exercise 41, 998–1005. doi: 10.1249/mss.0b013e3181930355
Duda, K. R., Jarchow, T., and Young, L. R. (2012). Squat exercise biomechanics during short-radius centrifugation. Aviat. Space Environ. Med. 83, 102–110.
Laing, C., Green, D. A., Mulder, E., Hinghofer−Szalkay, H., Blaber, A. P., Rittweger, J., et al. (2020). Effect of novel short−arm human centrifugation−induced gravitational gradients upon cardiovascular responses, cerebral perfusion and g−tolerance. J. Physiol. 598, 4237–4249. doi: 10.1113/jp273615
Lambert, E., and Lambert, G. W. (2014). Sympathetic dysfunction in vasovagal syncope and the postural orthostatic tachycardia syndrome. Front. Physiol. 5:280. doi: 10.3389/fphys.2014.00280
LeBlanc, A., Gogia, P., Schneider, V., Krebs, J., Schonfeld, E., and Evans, H. (1988). Calf muscle area and strength changes after five weeks of horizontal bed rest. Am. J. Sports Med. 16, 624–629. doi: 10.1177/036354658801600612
Lee, S. M. C., Feiveson, A. H., Stein, S., Stenger, M. B., and Platts, S. H. (2015). Orthostatic intolerance after ISS and spaceshuttle missions. Aerosp.Med.Hum.Perform. 86, 54–67. doi: 10.3357/AMHP.EC08.2015
Lehman, A., O’Rourke, N., Hatcher, L., and Stepanski, E. J. (2005). JMP for Basic Univariate and Multivariate Statistics: A Step-by-Step Guide. Cary, NC: SAS Institute. Inc.
Makivić, B., Nikić Djordjević, M., and Willis, M. S. (2013). Heart Rate Variability (HRV) as a tool for diagnostic and monitoring performance in sport and physical activities. J. Exercise Physiol. 16, 103–131.
Masatli, Z., Nordine, M., Maggioni, M. A., Mendt, S., Hilmer, B., Brauns, K., et al. (2018). Gender-specific cardiovascular reactions to + Gz interval training on a short arm human centrifuge. Front. Physiol. 9:1028. doi: 10.3389/fphys.2018.01028
McDonnell, A. C., Eiken, O., Frings−Meuthen, P., Rittweger, J., and Mekjavic, I. B. (2019). The LunHab project: muscle and bone alterations in male participants following a 10-day lunar habitat simulation. Exp. Physiol. 104, 1250–1261. doi: 10.1113/EP087482
Migeotte, P. F., Pattyn, N., Vanspauwen, R., Neyt, X., Acheroy, M., Van de Heyning, P., et al. (2009). Respiratory sinus arrhythmia on the ESA-short-arm human centrifuge. IEEE Eng. Med. Biol. Mag. 28, 86–91. doi: 10.1109/MEMB.2009.934618
Norsk, P. (2014). Blood pressure regulation IV: adaptive responses to weightlessness. Eur. J. Appl. Physiol. 114, 481–497. doi: 10.1007/s00421013-2797-2
Peterson, K., Ozawa, E. T., Pantalos, G. M., and Sharp, M. K. (2002). Numerical simulation of the influence of gravity and posture on cardiac performance. Ann. Biomed. Eng. 30, 247–259. doi: 10.1114/1.1451075
Powers, S. K., and Howley, E. T. (1997). Exercise Physiology: Theory and Application to Fitness and Performance. Edmonton: WCB.
Rittweger, J., and Frings-Meuthen, P. (2013). Bone loss in microgravity. Physiol. News 92, 38–41. doi: 10.36866/pn.92.38
Robertson, D. (1999). The epidemic of orthostatic tachycardia and orthostatic intolerance. Am. J. Med. Sci. 317, 75–77. doi: 10.1097/00000441-199902000-00001
Ruf, T. (1999). The lomb-scargle periodogram in biological rhythm research: analysis of incomplete and unequally spaced time-series. Biol. Rhythm Res. 30, 178–201. doi: 10.1076/brhm.30.2.178.1422
Saltin, B., Blomqvist, G., Mitchell, J. H., Johnson, R. L. Jr., Wildenthal, K., and Chapman, C. B. (1968). Response to exercise after bed rest and after training. Circulation 38, 71–78.
Sandler, H., and Vernikos, J. (1986). Inactivity: Physiological Effects. Orlando, FL: Academic Press Inc, doi: 10.1016/b978-0-12-618510-2.x5001-5
Sasaki, T., Iwasaki, K., Oka, T., Hisanaga, N., Ueda, T., Takada, Y., et al. (1999). Effect of working hours on cardiovascular-autonomic nervous functions in engineers in an electronics manufacturing company. Ind. Health 37, 55–61. doi: 10.2486/indhealth.37.55
Seps, B., Beckers, F., and Aubert, A. E. (2002). “Heart rate variability during gravity transitions,” in Computers in Cardiology (Piscataway, NJ: IEEE), 433–436.
Shoemaker, J. K., Hogeman, C. S., Khan, M., Kimmerly, D. S., and Sinoway, L. I. (2001). Gender affects sympathetic and hemodynamic response to postural stress. Am. J. Physiol. Heart Circ. Physiol. 281, H2028–H2035. doi: 10.1152/ajpheart.2001.281.5.H2028
Stenger, M. B., Evans, J. M., Knapp, C. F., Lee, S. M. C., Phillips, T. R., Perez, S. A., et al. (2012). Artificial gravity training reduces bed rest-induced cardiovascular deconditioning. Eur. J. Appl. Physiol. 112, 605–616. doi: 10.1007/s00421-0112005-1
Stenger, M. B., Evans, J. M., Patwardhan, A. R., Moore, F. B., Hinghofer-Szalkay, H., Rössler, A., et al. (2007). Artificial gravity training improves orthostatic tolerance in ambulatory men and women. Acta Astronaut. 60, 267–272. doi: 10.1016/j.actaastro.2006.08.008
Smith, J. J., Porth, C. M., and Erickson, M. (1994). Hemodynamic response to the upright posture. J. Clin. Pharmacol. 34, 375–386. doi: 10.1002/j.1552-4604.1994.tb04977.x
Tanaka, K., Nishimura, N., and Kawai, Y. (2016). Adaptation to microgravity, deconditioning, and countermeasures. J. Physiol. Sci. 67, 271–281. doi: 10.1007/s12576-016-0514-8
Verma, A. K., Xu, D., Bruner, M., Garg, A., Goswami, N., Blaber, A. P., et al. (2018). Comparison of autonomic control of blood pressure during standing and artificial gravity induced via short-arm human centrifuge. Front. Physiol. 9:712. doi: 10.3389/fphys.2018.00712
Vernikos, J. (1996). Human physiology in space. Bioessays 18, 1029–1037. doi: 10.1002/bies.950181215
Vernikos, J. (1997). Artificial gravity intermittent centrifugation as a space flight countermeasure. J. Gravit. Physiol. 4, 13–16.
Vernikos, J. (2017). “Gravity, sitting and health,” in Sedentary Behavior and Health, Concepts, Assessments and Interventions, eds W. Zhu and N. Owen (Urbana, IL: Human Kinetics), 13–29.
Vernikos, J., Ludwig, D. A., Ertl, A., Eade, C. E., Keil, L., and Ohara, D. B. (1996). Effect of standing or walking on physiological changes induced by head-down bedrest: implications for space flight. Aviat Space Environ. Med. 67, 1069–1079.
Vernikos, J., Walter, N., Worms, J. C., and Blanc, S. (2016). THESEUS: the European research priorities for human exploration of space. NPJ Microgravity 2:16034. doi: 10.1038/npjmgrav.2016.34
Wehrwein, E. A., and Joyner, M. J. (2013). Regulation of blood pressure by the arterial baroreflex and autonomic nervous system. Handb. Clin. Neurol. 117, 89–102. doi: 10.1016/b978-0-444-53491-0.00008-0
Wenner, M. M., Haddadin, A., Taylor, H. S., and Stachenfeld, N. S. (2013). Mechanisms contributing to low orthostatic tolerance in females: the influence of oestradiol. J. Physiol. 591, 2345–2355. doi: 10.1113/jphysiol.2012.247882
White, O., Barbiero, M., and Goswami, N. (2019). The effects of varying gravito-inertial stressors on grip strength and hemodynamic responses in men and women. Eur. J. Appl. Physiol. 119, 951–960. doi: 10.1007/s00421-019-04084-y
Wilmore, J. H., and Costill, D. L. (eds) (1994). “Cardiorespiratory adaptations to training,” in Physiology of Sport and Exercise, (Leeds: Human Kinetics), 214–238.
Winter, J., Laing, C., Johannes, B., Mulder, E., Brix, B., Roessler, A., et al. (2018). Galanin and adrenomedullin plasma responses during artificial gravity on a human short-arm centrifuge. Front. Physiol. 9:1956. doi: 10.3389/fphys.2018.01956
Yang, C. B., Zhang, S., Zhang, Y., Wang, B., Yao, Y. J., Wang, Y. C., et al. (2010). Combined short-arm centrifuge and aerobic exercise training improves cardiovascular function and physical working capacity in humans. Med. Sci. Monit. 16, 575–583.
Keywords: artificial gravity, blood gas analysis, cardiac output, deconditioning, heart rate variability, mean arterial pressure, short-arm human centrifuge
Citation: Kourtidou-Papadeli C, Frantzidis CA, Gilou S, Plomariti CE, Nday CM, Karnaras D, Bakas L, Bamidis PD and Vernikos J (2021) Gravity Threshold and Dose Response Relationships: Health Benefits Using a Short Arm Human Centrifuge. Front. Physiol. 12:644661. doi: 10.3389/fphys.2021.644661
Received: 21 December 2020; Accepted: 12 April 2021;
Published: 11 May 2021.
Edited by:
Marc-Antoine Custaud, Université d’Angers, FranceReviewed by:
Nandu Goswami, Medical University of Graz, AustriaAlan R. Hargens, University of California, San Diego, United States
Copyright © 2021 Kourtidou-Papadeli, Frantzidis, Gilou, Plomariti, Nday, Karnaras, Bakas, Bamidis and Vernikos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chrysoula Kourtidou-Papadeli, cGFwYWRjQGF1dGguZ3I=
†These authors have contributed equally to this work and share senior authorship