Original Research ARTICLE
Long-Term Microgravity Exposure Increases ECG Repolarization Instability Manifested by Low-Frequency Oscillations of T-Wave Vector
- 1BSICoS Group, Aragón Institute of Engineering Research, IIS Aragón, Universidad de Zaragoza, Zaragoza, Spain
- 2Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, Milan, Italy
- 3CIBER en Bioingeniería, Biomateriales y Nanomedicina, Madrid, Spain
Ventricular arrhythmias and sudden cardiac death during long-term space missions are a major concern for space agencies. Long-duration spaceflight and its ground-based analog head-down bed rest (HDBR) have been reported to markedly alter autonomic and cardiac functioning, particularly affecting ventricular repolarization of the electrocardiogram (ECG). In this study, novel methods are developed, departing from previously published methodologies, to quantify the index of Periodic Repolarization Dynamics (PRD), an arrhythmic risk marker that characterizes sympathetically-mediated low-frequency oscillations in the T-wave vector. PRD is evaluated in ECGs from 42 volunteers at rest and during an orthostatic tilt table test recorded before and after 60-day –6° HDBR. Our results indicate that tilt test, on top of enhancing sympathetic regulation of heart rate, notably increases PRD, both before and after HDBR, thus supporting previous evidence on PRD being an indicator of sympathetic modulation of ventricular repolarization. Importantly, long-term microgravity exposure is shown to lead to significant increases in PRD, both when evaluated at rest and, even more notably, in response to tilt test. The extent of microgravity-induced changes in PRD has been associated with arrhythmic risk in prior studies. An exercise-based, but not a nutrition-based, countermeasure is able to partially reverse microgravity-induced effects on PRD. In conclusion, long-term exposure to microgravity conditions leads to elevated low-frequency oscillations of ventricular repolarization, which are potentiated following sympathetic stimulation and are related to increased risk for repolarization instabilities and arrhythmias. Tested countermeasures are only partially effective in counteracting microgravity effects.
After almost 60 years of human spaceflight, there is good evidence on detrimental effects on the human body associated with long-term space missions (Williams et al., 2009; Garrett-Bakelman et al., 2019). Two of the main causes underlying those effects are ionizing radiation and changes in gravity conditions. Specifically, microgravity-induced cardiac arrhythmias are a major concern for national space agencies, as very prolonged periods of time in the International Space Station or in a mission to Mars or the Moon might set the stage for the development of ventricular tachycardia or ventricular fibrillation that could end up in sudden cardiac death (Anzai et al., 2014; Caiani et al., 2016). Although the probability of undergoing serious cardiac arrhythmias in the course of a space mission is low, with the estimated probability of suffering a life-threatening event being of 1% per year in short to mid-duration spaceflights (Russomano et al., 2013), currently available data are limited and more sophisticated techniques should be employed to identify potential in-flight abnormalities in the electrical activity of the heart (Convertino, 2009).
Several factors may enhance predisposition to ventricular arrhythmias during spaceflight. Commonly reported bradycardia (Meck et al., 2001), changes in electrolyte composition of blood plasma (Smith and Zwart, 2008), psychological stress (Kanas et al., 2001) and, very relevantly, adaptation of cardiac autonomic modulation (Fritsch-Yelle et al., 1996) may all concur to adversely affect ventricular electrophysiology. In particular, reported alterations in the sympathetic nervous system might contribute to the documented increase in spatio-temporal inhomogeneity of ventricular repolarization, thus potentially providing an electrophysiological substrate for arrhythmias (Caiani et al., 2016). Nevertheless, further evidence on elevated arrhythmic risk during long-term space missions and its underlying mechanisms is yet to be established.
Studies assessing microgravity effects on ventricular repolarization during or immediately after spaceflight are limited. Major findings indicate that long-duration spaceflight prolongs cardiac repolarization, as measured by the QT corrected interval of the electrocardiogram (ECG) (D'Aunno et al., 2003). Due to the limited opportunities to obtain data from humans in space missions, mainly related to the hazards and high costs of spaceflight investigations, several ground-based models have been used to simulate space conditions, explore potential adverse effects associated with weightlessness and assess the effectiveness of proposed countermeasures. Long-term head-down (–6°) bed rest (HDBR) is a ground-based analog widely utilized to simulate microgravity effects on the human body (Pavy-Le Traon et al., 2007; Hargens and Vico, 2016). Relevant alterations in ventricular repolarization have been reported in HDBR studies. In a 90-day HDBR investigation, several subjects were reported to develop QRS-T angles above 100° (Sakowski et al., 2011), with these elevated values having been associated with 3- to 5-fold increased risk for cardiovascular mortality and sudden death in previous works (Kardys et al., 2003; Yamazaki et al., 2005). In another study of only 9- to 16-day HDBR, simulated microgravity was shown to lead to an increase in microvolt T-wave alternans (Grenon et al., 2005), a well-known marker of ventricular arrhythmias and sudden cardiac death. Of note, HDBR-induced increases in T-wave alternans correlated with changes in sympathetic function. In another short-term HDBR study (Mart́ın-Yebra et al., 2015) TWA was, however, shown not to increase during stress-test and tilt-table test after 5- and 21-day HDBR experiments. Interestingly, subjects suffering a more marked orthostatic intolerance after HDBR were found to be those presenting greater values of TWA, even before exposure to simulated microgravity.
A myriad of indices have been reported in the literature to assess ECG repolarization, including prolongation of the QTc interval (Mitchell and Meck, 2004), QT rate adaptation (Pueyo et al., 2004, 2008), QT interval variability (Piccirillo et al., 2009), T-wave alternans (Rosenbaum et al., 1994; Martinez and Olmos, 2005), or T-wave morphological variability (Adam et al., 1984; Badilini et al., 1997; Acar et al., 1999; Baumert et al., 2011; Ramírez et al., 2017a,b), among others. An index of Periodic Repolarization Dynamics (PRD) has been recently proposed to assess sympathetic modulation of ventricular repolarization by measuring low-frequency (below 0.1 Hz) oscillations in the T-wave vector (Rizas et al., 2014). PRD accounts for variability not limited to a specific time interval, as the QT interval or the T-peak-to-T-end interval, but more generally integrating all the spatio-temporal information in the T-wave vector, which can allow for a more robust characterization of beat-to-beat repolarization variations and can provide a better marker to anticipate electrical instabilities (Rizas et al., 2014, 2016).
In this study, ECG signals from healthy volunteers undergoing 60-day HDBR are analyzed. PRD is hypothesized to be able to characterize the effects of sustained simulated microgravity on ventricular repolarization, particularly in response to an orthostatic Tilt-Table Test (TTT), a common procedure used to assess autonomic nervous system function (Zygmunt and Stanczyk, 2010). To overcome identified issues related to angle quantification as part of the PRD technique, a number of updates on the originally reported methods (Rizas et al., 2014, 2016) are also proposed. Additionally, the effectiveness of two different countermeasures, based on exercise and nutrition, to mitigate or reduce microgravity-induced effects on ventricular repolarization during HDBR are assessed.
2. Materials and Methods
2.1. Study Population
Data from two 60-day –6° HDBR campaigns organized by the European Space Agency (ESA) as part of ESA bed rest studies were analyzed in this work. These studies were conducted between 2015 and 2017 in the :envihab facility of the Institute of Aerospace Medicine at the German Aerospace Center-DLR (Cologne, Germany) and at the Institute of Space Medicine and Physiology-MEDES (Toulouse, France).
For the experiment in Cologne, 22 male volunteers (29 ± 6 years, 181 ± 5 cm, 77 ± 7 kg) were enrolled and randomly distributed into either the countermeasure group (JUMP), who performed 48 training sessions of a varying number of countermovement jumps on a sledge jump system during the HDBR time period (Kramer et al., 2017), or the control group (CTRL), who did not perform any exercise. For the experiment in Toulouse, 20 male volunteers were enrolled (34 ± 7 years, 176 ± 4 cm, 73 ± 7 kg). They were randomly distributed into either the countermeasure group (NUTR), daily receiving a nutritional countermeasure consisting of a cocktail of anti-oxidants and vitamins (daily, 530 mg of polyphenol, 168 mg of vitamin E, 80 μg of Selenium-Solgar®, and 2.1 g of Omega-3—Omacor®), or the control group (CTRL), who did not receive this nutritional integration.
All subjects underwent prior comprehensive medical examination during the selection process and provided written informed consent to participate in the study, which was approved in advance by the respective Ethical Committees for Human Research at the host institutions.
2.2. Experimental Protocol
Both campaigns were divided into three phases: 15 days of PRE-HDBR baseline (BDC-15 to BDC-1), when the subjects became acclimated physiologically and psychologically to the facilities; 60 days of bed rest (HDT1 to HDT60), when subjects were in strict –6° HDBR (24 h/day); and 15 days of POST-HDBR recovery (R+0 to R+14). Figure 1 illustrates these three phases. From HDT1 to HDT60, subjects carried out all activities at –6° HDBR: eating, hygienic procedures (teeth brushing, bowel movement, showering) and free time activities (reading, watching, or using computer). Also, all subjects had the same scheduled wake-up (at 6:30 a.m. and 7:00 a.m. in DLR and MEDES campaign, respectively) and light off (at 11:00 p.m.). More information about the study protocol is available in Kramer et al. (2017).
Figure 1. Phases of the head-down bed rest (HDBR) campaign, with indication of the days when volunteers underwent tilt table tests: 2 days before the start of the HDBR period (BCD-2) and just after completing it (R+0).
Two TTTs were performed, one of them 2 days before the start of the HDBR period (BCD-2) and the other one just after completing it (R+0). In each TTT the subject was tilted head-up to an angle of 80° for up to 15 min. If the subject did not experience any presyncopal episode during that time, he was exposed to Lower Body Negative Pressures (LBNP) following a protocol of 3-min –10 mmHg steps for a maximum duration of 15 min. Thirty out of the 84 analyzed recordings did not present presyncopal episodes, of which 18 corresponded to PRE-HDBR and 12 to POST-HDBR.
High-resolution (1,000 Hz) 24-h Holter 12-lead ECG signals (Mortara Instrument) recorded at days BCD-2 and R+0 (both including a TTT) were available for this study. For each TTT, a 5-min interval prior to the start of the tilt phase, the first 5 min immediately following its start and the last 5 min of the tilt phase (possibly including LBNP) were analyzed (Figure 2). If the tilt phase lasted for less than 5 min, its whole duration was analyzed.
2.3. ECG Pre-processing
Raw ECG signals were pre-processed by a 50 Hz notch filter to remove powerline interference. Taking these pre-processed ECG signals as inputs, QRS detection and ECG wave delineation were performed by using a wavelet-based single-lead automatic system (Martinez et al., 2004). The outputs of the detection and delineation system were combined by using rules to obtain multi-lead ECG delineation marks (Martinez et al., 2004). Since subsequent analysis focused on the T-wave, a 40-Hz low-pass filter was applied to remove noise without altering the T-wave shape. Finally, cubic splines interpolation was applied to estimate and remove baseline wander. An example of an ECG recording as originally acquired and after application of different pre-processing steps is shown in Figure 3.
Figure 3. Example of an ECG recording as originally acquired and after application of different pre-processing steps.
2.4. Calculation of Angles Between Consecutive T Waves
An updated method based on the original method proposed in Rizas et al. (2014) was applied onto the pre-processed ECG signals to compute the angles between consecutive T-waves:
1. Orthogonal leads X, Y, Z were obtained from the 12-lead ECG by using the inverse Dower matrix (Edenbrandt and Pahlm, 1988).
2. The onset and end of each T-wave, denoted by Ton and Toff, were identified by the delineation system described above. When the delineation failed to identify a T-wave onset (T-wave offset, respectively) for a given beat, its location was set based on Ton (Toff, respectively) locations for adjacent beats with respect to their corresponding QRS positions. Specifically, Ton (or Toff) was located at a distance from the corresponding QRS fiducial point equal to the median interval between the QRS and Ton (or Toff) positions of 30 beats around that beat.
As T-wave boundaries change on a beat-to-beat basis, and may be influenced by delineation errors, the angle between each two consecutive T-waves was computed by defining a unique temporal window for both waves being analyzed. Specifically, for each angle calculation, the window onset was set at the latest Ton of both analyzed beats computed with respect to their QRS fiducial points, while the window end was set at the earliest Toff of both beats computed from their QRS fiducial points.
3. A constant value was subtracted from each T wave in each of the analyzed leads so that the amplitude at Toff was set to 0 mV. Subsequently, an average T-wave vector was calculated for each T-wave. The angle dT° between two consecutive T-waves, which is associated with the instantaneous degree of repolarization instability, was calculated by using the dot product of each pair of consecutive average T-wave vectors.
4. The dT° time series was filtered by using a 10th-order median filter to attenuate outliers and avoid very abrupt changes in the time series.
2.5. PRD Computation
Two different methods, based on Continuous Wavelet Transform (CWT) and Phase-Rectified Signal Averaging (PRSA), respectively, were developed based on the initial methodology proposed in Rizas et al. (2014, 2016). These methods were tested for quantification of the low-frequency components of the beat-to-beat dT° series. The steps followed in each of the two methods are depicted in Figure 4.
2.5.1. PRD Computation Using Continuous Wavelet Transform
CWT is one of the most widely-used tools for time-frequency analysis (Addison, 2005). Based on the dT° series calculated as described in section 2.4, the next steps were followed to compute PRD (Rizas et al., 2014):
5. The dT° series was linearly interpolated at 2 Hz and a 10-sample moving average filter was used to remove artifacts.
6. CWT was computed at all scales from 1 to 40 by using a 4th-order Gaussian wavelet to quantify low-frequency oscillations of dT°. Wavelet coefficients were obtained for each scale at each time point and an average wavelet coefficient was computed for each scale.
7. Scales (a) were converted to pseudo-frequencies (Fa, expressed in Hz) according to the following equation (Abry, 1997):
where Fc is the center frequency of the mother wavelet, in Hz, and Δ denotes the sampling period, in seconds.
PRDCWT was defined as the average wavelet coefficient in the frequency range between 0.025 and 0.1 Hz.
2.5.2. PRD Computation Using Phase-Rectified Signal Averaging
An alternative method to compute oscillatory fluctuations, with less computational requirements, has been proposed based on PRSA (Bauer et al., 2006). The following steps were followed to compute PRD from the dT° series (Rizas et al., 2016):
1. Anchor points were defined by comparing averages of M = 9 values of the dT° series previous and posterior to the anchor point candidate (xi). A beat i is considered an anchor point if:
2. Windows of 2 L values were defined around each anchor point. If an anchor point was so close to the beginning or to the end of the dT° series that there were not enough samples before or after it, it was disregarded. In this study, L = 20 was chosen because it was the minimum value to detect frequencies in the range of interest (0.025–0.1) Hz.
3. PRSA series were obtained by averaging the dT° series over all defined windows.
PRDPRSA was defined as the difference between maximum and minimum values of the PRSA series.
2.6. Heart Rate Variability Analysis
RR interval series were computed from the QRS detection marks obtained in section 2.3 for all analyzed 5-min segments at baseline as well as at the beginning and end of TTT. Instantaneous heart rate (HR) variability (HRV) series were calculated following the method described in Bailón et al. (2011). For each segment, the power spectral density (PSD) of HRV was computed by using the periodogram method. A high-frequency band (HF, [0.15, 0.4] Hz) and a low-frequency band (LF, [0.04, 0.15] Hz) were defined for HRV analysis in the frequency domain and the LF and HF powers were calculated by integrating the power spectrum in each of those two bands, respectively. The normalized LF power (LFn), the normalized HF power (HFn) and the ratio between the power in the LF and HF bands (LF/HF) were computed (Malik et al., 1996). Also the median HR (HRmedian) was computed.
2.7. Statistical Analysis
The Mann-Whitney U-test (or Wilcoxon rank-sum test) was used to compare independent samples, as when comparing each countermeasure (JUMP or NUTR) subgroup vs the corresponding CTRL subgroup. Wilcoxon signed-rank test was used for comparison of paired samples, as when comparing changes induced by HDBR or by TTT in a given group of subjects. Spearman's correlation coefficient ρ and Kendall's τ were used to quantify rank correlation between CWT and PRSA. All statistical analyses were carried out using MATLAB R2017a (9.2).
3.1. Comparison of PRD Computed by CWT- and PRSA-Based Methods
Figure 5 shows the two analyzed recordings, at PRE-HDBR and POST-HDBR, from a volunteer presenting small and large magnitudes of low-frequency oscillations in ventricular repolarization, respectively. The three plots represent the dT° series, the frequency pseudospectra (in terms of squared wavelet coefficients) and the PRSA series. The blue line corresponds to the ECG segment at PRE-HDBR and the red one to the ECG segment at POST-HDBR. Note that the case shown in blue presents low-frequency oscillations in dT° of small magnitude, which translates into low values of PRDCWT and PRDPRSA. The red case, in contrast, presents low-frequency oscillations in dT° of larger magnitude and is associated with considerably higher PRD values, both when measured by using CWT- and PRSA-based methods.
Figure 5. Examples of dT° series (upper), frequency pseudospectra (bottom left), and PRSA series (bottom right) for two ECG segments from a volunteer of the study, at PRE-HDBR and at POST-HDBR, presenting remarkably different magnitudes of low-frequency oscillations in ventricular repolarization. Associated PRD values are indicated in the bottom panels, as computed using CWT- and PRSA-based methods.
Figure 6 shows the correlation of PRD values computed by using the CWT-based method (X-axis) and the PRSA-based method (Y-axis) for all analyzed segments (baseline, beginning, and end of the tilt phase) in the CTRL group of DLR and MEDES campaigns, for both PRE-HDBR and POST-HDBR. The scatterplot shows a strong correlation between both methods. Rank correlation coefficients were: Spearman's ρ = 0.93 (p < 10−50), Kendall's τ = 0.79 (p < 10−35). In the following, all presented results use the PRSA-based method.
3.2. Tilt Test-Induced Effects on PRD
Figure 7, top panel, shows the results of the analysis of three 5-min segments of the dT° series, corresponding to baseline (prior to the tilt) as well as to beginning and end of the tilt phase, for all volunteers in the CTRL group of both campaigns (DLR and MEDES). Results are separately presented for PRE-HDBR (before HDBR) and POST-HDBR (after HDBR). As can be observed from the figure, PRD increased following tilt as compared to baseline, being the results statistically significant when the segment at the beginning of the tilt phase was analyzed. This was true for both PRE-HDBR and POST-HDBR. Results on the effects of tilt on the HRV indices LFn and LF/HF are presented in the bottom panels of Figure 7. Both indices showed significantly larger values in response to tilt, indicating increased sympathetic drive during orthostatic stress, both at PRE-HDBR and POST-HDBR. The effect of tilt on other HR and HRV indices is presented in Figure S1, which shows a significant increase of HRmedian and a significant decrease of HFn in response to tilt.
Figure 7. Boxplots of PRD, LFn and LF/HF at PRE-HDBR (in blue) and POST-HDBR (in red) evaluated at baseline and at the beginning and end of the tilt phase. ∙∙p < 0.01, ∙p < 0.05 (Wilcoxon signed-rank test).
3.3. Microgravity-Induced Effects on PRD
The effects of microgravity exposure on PRD obtained by comparing PRE-HDBR and POST-HDBR for CTRL group of DLR and MEDES campaigns can be observed from Figure 7. At baseline (before TTT), PRD was significantly increased at POST-HDBR with respect to PRE-HDBR, changing from 0.85 [0.31] deg at PRE-HDBR to 1.18 [0.51] deg at POST-HDBR (p < 0.01), as presented in the left columns of Figure 7.
Considering the analysis at the onset of the tilt phase, PRD was also increased at POST-HDBR with respect to PRE-HDBR, changing from 1.24 [0.72] deg at PRE-HDBR to 1.89 [0.91] deg at POST-HDBR (p < 0.05), as shown in the middle columns of Figure 7.
No statistically significant differences in PRD at POST-HDBR vs. PRE-HDBR were observed when analyzed at the end of the tilt phase (right columns in Figure 7), although there was a trend to increased PRD at POST-HDBR as compared to PRE-HDBR. Specifically, PRD increased from 1.18 [0.88] deg at PRE-HDBR to 1.27 [1.21] deg at POST-HDBR (n.s.).
In the case of the HRV indices LFn and LF/HF, statistically significant microgravity-induced increases were observed when the baseline period was analyzed (see Figure 7). Specifically, LFn changed from 0.67 [0.27] n.u. at PRE-HDBR to 0.78 [0.17] n.u. at POST-HDBR (p < 0.01). LF/HF changed from 2.12 [2.51] at PRE-HDBR to 3.74 [3.78] at POST-HDBR (p < 0.01). The index HFn significantly decreased from PRE-HDBR to POST-HDBR when evaluated at baseline, whereas HRmedian was significantly augmented due to microgravity either when evaluated at baseline or at the beginning and end of the tilt test (see Figure S1).
3.4. PRD and HRV Relation
Figure S2 shows the relationship between tilt-induced changes in PRD and in HR or HRV indexes (LFn, LF/HF and HRmedian) at PRE-HDBR and POST-HDBR. No significant correlation could be found between PRD and HR or HRV, with Spearman's correlation coefficient ρ being below 0.15 in all evaluated cases (n.s.).
3.5. Effectiveness of Exercise-Based Countermeasure
The ability of a jump-based countermeasure to reverse the effects of microgravity was evaluated by comparing PRD values at PRE-HDBR and POST-HDBR in each of the CTRL and JUMP subgroups of the DLR campaign. The values for PRD measured for each phase of the TTT are presented in Table 1. Although there were increases in PRD values from PRE-HDBR to POST-HDBR in both CTRL and JUMP subgroups, the increase was much more attenuated in the JUMP subgroup. Significant differences were found at the beginning of tilt for both CTRL and JUMP subgroups. Illustration of the effects of the JUMP countermeasure are presented in Figure 8 (left panel), which shows values of ΔPRD, calculated as the PRD value at POST-HDBR minus the PRD value at PRE-HDBR for each subject. From the figure it is clear that whereas values of ΔPRD were clearly positive in the CTRL subgroup, particularly during the tilt phase, values were remarkably closer to 0 in the JUMP subgroup.
Table 1. PRD values (median [IQR]) at all phases of TTT for PRE-HDBR and POST-HDBR in CTRL and JUMP subgroup.
Figure 8. ΔPRD values (median ± IQR) measured respect to PRE-HDBR for the CTRL subgroups (in orange) and the countermeasure (JUMP or NUTR) subgroups (in green) at baseline, beginning and end of the tilt phase. (n.s., Wilcoxon rank-sum test).
3.6. Effectiveness of Nutrition-Based Countermeasure
Results on the effectiveness of a nutrition-based countermeasure are presented in Table 2. Figure 8 (right panel) illustrates these results in terms of ΔPRD (differences between POST-HDBR and PRE-HDBR calculated for each subject in the analyzed subgroups). As can be observed from Table 2, baseline PRD increased significantly from PRE-HDBR to POST-HDBR in both the CTRL and NUTR subgroups of the MEDES campaign. When evaluation was performed at the beginning of the tilt test, PRD increased at POST-HDBR with respect to PRE-HDBR, although differences were not statistically significant. At the end of the tilt phase, PRD showed a trend of increase in the NUTR group but not in the CTRL group. Results shown in Figure 8 (right panel) confirm the lack of effectiveness of the evaluated nutrition-based countermeasure.
Table 2. PRD values (median [IQR]) at all phases of TTT for PRE-HDBR and POST-HDBR in CTRL and NUTR subgroup.
This study aimed at investigating alterations in ventricular repolarization associated with long-term exposure to simulated microgravity conditions elicited by 60-day HDBR. Two methods have been developed for quantification of low-frequency oscillations in the T-wave of the ECG, departing from the original methodology proposed in Rizas et al. (2014, 2016). These methods, one using CWT and the other one using PRSA, have been shown to render concordant results in terms of the index of Periodic Repolarization Dynamics, PRD, a marker of low-frequency repolarization oscillations whose increase has been shown to be predictor of ventricular arrhythmias and sudden cardiac death (Rizas et al., 2014, 2016). This study has proved that microgravity remarkably enhances PRD, particularly when evaluated in response to sympathetic stimulation induced by tilt test. An exercise-based countermeasure has been shown to partially reverse microgravity-induced effects on PRD, whereas a nutrition-based countermeasure has been shown not to be effective at all.
The methods developed in this study for PRD quantification departed from the CWT- and PRSA-based methods proposed in Rizas et al. (2014, 2016), respectively. Whereas, the CWT-based method in Rizas et al. (2014) used spherical coordinates, our method used Cartesian coordinates, which rendered improved results for cases where T-wave vectors were close to the axes. Also, our method introduced a refinement on the temporal window used for T-wave definition so as to guarantee that the two consecutive T-waves involved in each angle computation had comparable T-wave window beginnings and ends with respect to their corresponding QRS fiducial points. An additional difference regards the number of samples used for the moving average filter, which was 30 in Rizas et al. (2014) and 10 in our study to minimize distortion of relevant information in the frequency band of interest. For our updated CWT- and PRSA-based methods, correlation analysis has confirmed a strong agreement between them. Although our PRD values are notably different from those obtained in Rizas et al. (2014, 2016), the agreement between CWT- and PRSA-based methods is in concordance with the findings reported in Rizas et al. (2016), where an approach based on PRSA was presented as an alternative to the approach using the CWT technique. The advantage of the PRSA approach is that it highly reduces the computational cost associated with PRD computation.
The analysis conducted in this work has shown that the autonomic changes induced by TTT are manifested as an increase in PRD, both when measured at PRE-HDBR and at POST-HDBR. Such PRD changes could be attributable to an increased sympathetic drive, as indicated by increases in the HRV indices LFn and LF/HF, in line with many other HRV studies, including the ones pioneering spectral HRV analysis during TTT (Pagani et al., 1986, 1988). It is well-known that sympathetic stimulation influences ventricular repolarization and modifies the characteristics of the T-wave in the ECG (Ramirez et al., 2011). Our results showing an increase in PRD in response to TTT are in line with the changes in PRD reported in response to variations in sympathetic activity or β-adrenergic modulation (Rizas et al., 2014, 2016). In our study, those changes are shown not to be explained by HRV changes but to reflect direct autonomic modulation of the ventricular myocardium, in accordance with the findings reported in Rizas et al. (2014, 2016). In vivo studies in patients have demonstrated that the same low-frequency oscillatory behavior of ventricular repolarization occurs locally, as measured from activation recovery intervals (ARIs) obtained from unipolar epicardial electrograms during ventricular pacing (Hanson et al., 2014; Porter et al., 2018). In those studies, heightened arousal of the sympathetic nervous system was elicited and maintained by mental stress or by Valsalva maneuver, which allowed characterization of low-frequency oscillations in ARI, a surrogate of action potential duration (APD), showing that those oscillations are coupled to oscillations in systolic and diastolic blood pressure (Hanson et al., 2014; Porter et al., 2018). Computational studies have provided insight into the mechanisms underlying sympathetically-mediated low-frequency oscillations of APD and the observed inter-individual differences (Pueyo et al., 2016; Sampedro-Puente et al., 2019). Specifically, phasic changes in both β-adrenergic stimulation and hemodynamic loading, a known accompaniment of enhanced sympathetic activity, have been demonstrated to contribute to low-frequency oscillations in APD, with these two actions being synergistic (Pueyo et al., 2016). Ionic differences in the densities of the L-type calcium (ICaL), rapid delayed rectifier potassium (IKr), and inwardly rectifier potassium (IK1) currents have been identified as the main drivers of inter-individual differences in the magnitude of low-frequency APD oscillations (Sampedro-Puente et al., 2019).
Importantly, our results have provided evidence on significant effects of long-duration microgravity simulation on cardiac electrical activity. In line with previously published studies, this work has confirmed that microgravity markedly alters ventricular repolarization (D'Aunno et al., 2003; Grenon et al., 2005; Sakowski et al., 2011; Bolea et al., 2012, 2013; Caiani et al., 2013), with those alterations being more manifested when evaluated in response to sympathetic stimulation. This study adds one more T-wave characteristic to the list of ECG repolarization properties proved to be modulated by microgravity. The quantified PRD index represents a form of temporal variability in ventricular repolarization, specifically focused on oscillations of frequencies below 0.1 Hz. Although other measures of ECG temporal variability have been investigated during or immediately after simulated microgravity exposure (Sakowski et al., 2011; Bolea et al., 2013), PRD can provide a more robust characterization of repolarization instability by encompassing global T-wave vector information. Also, the PRD index, by accounting for frequencies below 0.1 Hz, has been proven to be related to sympathetic modulation of ventricular repolarization (Rizas et al., 2016). On the basis that augmented sympathetic activity is associated with adverse outcomes in different patient populations (Verrier and Antzelevitch, 2004), the evaluated PRD index is of great interest for risk prediction. The enhancement of spatial and/or temporal ventricular heterogeneities observed in this and other studies in relation to long-term exposure to microgravity conditions suggest that microgravity could accentuate repolarization instability and thus increase ventricular arrhythmic risk, especially immediately upon gravity restoration. In particular, this study has shown that PRD quantified following 60-day HDBR is highly elevated, up to 50% at rest and up to 100% in response to TTT, with respect to PRE-HDBR values. The extent of change in PRD values measured immediately after 60-day HDBR could be associated with high arrhythmic risk taking as a reference previous studies on risk assessment in post-myocardial infarction patients, where those extents of change were found in patients who died vs. those who survived during follow-up (Rizas et al., 2014, 2017). This is in line with other studies that have reported on subjects presenting long-term microgravity-induced changes in ECG repolarization of an extent similar to those associated with more than 3-fold increased hazard ratio for sudden cardiac death in general populations (Sakowski et al., 2011).
Additionally, this study has assessed two countermeasures in their ability to counteract microgravity-induced effects on ventricular repolarization. The first applied countermeasure, based on an exercise training protocol, although markedly attenuated microgravity effects as measured by changes in the PRD index, it was not able to completely reverse them. These results on partial effectiveness of exercise-based countermeasures are in line with the findings reported in Kramer et al. (2017), Maggioni et al. (2018), and Caiani et al. (2018), which investigated the same jump-based countermeasure to reverse musculoskeletal and cardiovascular deconditioning. In other studies, exercise-based countermeasures have shown to be very effective in preserving bone and muscular conditions (McRae et al., 2012; Kramer et al., 2017; Maggioni et al., 2018).
The second tested countermeasure, a nutritional supplementation composed of an anti-oxidant and anti-inflammatory dietary mix, has been shown to be far from being effective in reducing microgravity-induced effects on ventricular repolarization. This is in agreement with other studies pointing out to lack of effectiveness of this countermeasure in counteracting microgravity exposure effects on bone turnover (Austermann et al., 2019). Importantly, the intake of omega-3 fatty acids, which are components of the dietary mix, and their possible protection of cardiovascular health should additionally be viewed in relation to the potentially increased risk for ventricular arrhythmias. Such a relation is, nevertheless, controversial, with some studies suggesting that they have detrimental arrhythmogenic effects, whereas other postulate minimal effects or highly anti-arrhythmic potential (Albert, 2012; von Schacky, 2012; Coronel, 2017; Tribulova et al., 2017). Although one reason to include this type of acid in a dietary support was its protective effects on bones (Zwart et al., 2010), the findings of the present study point out that the dietary mix could not reduce adverse cardiac effects of microgravity simulation. Further studies including larger number of subjects are needed to confirm or refute these findings. Also, it is relevant to note that, when evaluating the effects of the tested countermeasures, the CTRL subgroups of the JUMP and NUTR studies did not share the same ventricular repolarization characteristics as evaluated by PRD, despite the subjects of both studies having similar physical conditions. Specifically, subjects in the CTRL subgroup of the JUMP study presented higher values of PRD, both at PRE-HDBR and POST-HDBR. Because of that reason, our results on countermeasure effects were assessed in relative terms. Nevertheless, the inclusion of a larger number of subjects would definitely allow more robust analysis of absolute and relative microgravity-induced changes. Additionally, future studies could test other types of nutritional supplements to improve the ability to counteract deleterious effects associated with long-term microgravity exposure (Cena et al., 2003). Based on the results of this study and the concordance with the outcomes of other studies, a modified jump training or a combination of exercise- and nutrition-based countermeasure (Schneider et al., 2009; Konda et al., 2019; Kramer et al., 2019), possibly including other components like pharmacological agents or artificial gravity (Evans et al., 2018), would be suggested to compensate for adverse microgravity-induced effects on ventricular repolarization.
The effects of long-duration microgravity on ventricular repolarization have been assessed by evaluation of the PRD index, a marker of low-frequency repolarization oscillations whose increase is related to high risk for ventricular arrhythmias and sudden cardiac death. Two methods have been developed for robust quantification of PRD, which have shown to present very good agreement. Long-term microgravity exposure has been proven to markedly elevate PRD, particularly when evaluated in response to enhanced sympathetic activity induced by a tilt table test. A countermeasure based on exercise training has been shown to partially counteract microgravity-induced changes in ventricular repolarization as assessed immediately upon gravity restoration.
Data Availability Statement
The datasets analyzed in this article are not publicly available. Requests to access the datasets should be directed to European Space Agency.
The studies involving human participants were reviewed and approved by Institute of Aerospace Medicine—German Aerospace Center-DLR and by Institute of Space Medicine and Physiology-MEDES. The patients/participants provided their written informed consent to participate in this study.
EP and JM devised the project, the main conceptual ideas and proof outline, and were responsible for overseeing the research and providing critical insight and recommendations regarding the focus, structure, and content of the paper. SP performed computational simulations and analyzed the data results. EC was responsible for the definition of the bed rest data acquisition protocols and contributed with technical details and analysis support. FL contributed by managing on-site data acquisition in both campaigns. All authors participated in writing and proofreading throughout the publication process.
This work was supported by projects ERC-2014-StG 638284 (ERC), DPI2016-75458-R (MINECO) and Aragón Government (Reference Group BSICoS T39-17R and project LMP124-18) cofunded by FEDER 2014-2020 Building Europe from Aragon. We also acknowledge the support of the Italian Space Agency (contract 2018-7-U.0, PI and recipient EC).
Conflict of Interest
The 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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2019.01510/full#supplementary-material
Adam, D. R., Smith, J. M., Akselrod, S., Nyberg, S., Powell, A. O., and Cohen, R. J. (1984). Fluctuations in t-wave morphology and susceptibility to ventricular fibrillation. J. Electrocardiol. 17, 209–218.
Austermann, K., Baecker, N., Zwart, S. R., Smith, S. M., and Heer, M. (2019). Effects of antioxidants on bone turnover markers in 6° head-down tilt bed rest. Front. Physiol. 9:48. doi: 10.3389/conf.fphys.2018.26.00048
Badilini, F., Fayn, J., Maison-Blanche, P., Leenhardt, A., Forlini, M. C., Denjoy, I., et al. (1997). Quantitative aspects of ventricular repolarization: relationship between three-dimensional T wave loop morphology and scalar qt dispersion. Ann. Noninvasive Electrocardiol. 2, 146–157.
Bailón, R., Laouini, G., Grao, C., Orini, M., Laguna, P., and Meste, O. (2011). The integral pulse frequency modulation model with time-varying threshold: application to heart rate variability analysis during exercise stress testing. IEEE Trans. Biomed. Eng. 58, 642–652. doi: 10.1109/TBME.2010.2095011
Bauer, A., Kantelhardt, J. W., Bunde, A., Barthel, P., Schneider, R., Malik, M., et al. (2006). Phase-rectified signal averaging detects quasi-periodicities in non-stationary data. Phys. A Stat. Mech. Appl. 364, 423–434. doi: 10.1016/j.physa.2005.08.080
Baumert, M., Lambert, E., Vaddadi, G., Sari, C. I., Esler, M., Lambert, G., et al. (2011). Cardiac repolarization variability in patients with postural tachycardia syndrome during graded head-up tilt. Clin. Neurophysiol. 122, 405–409. doi: 10.1016/j.clinph.2010.06.017
Bolea, J., Caiani, E. G., Pueyo, E., Laguna, P., and Almeida, R. (2012). “Microgravity effects on ventricular response to heart rate changes,” in Engineering in Medicine and Biology Society (EMBC), 2012 Annual International Conference of the IEEE (San Diego, CA: IEEE), 3424–3427.
Bolea, J., Laguna, P., Caiani, E. G., and Almeida, R. (2013). “Heart rate and ventricular repolarization variabilities interactions modification by microgravity simulation during head-down bed rest test,” in Computer-Based Medical Systems (CBMS), 2013 IEEE 26th International Symposium on (Porto: IEEE), 552–553.
Caiani, E., Landreani, F., Costantini, L., Mulder, E., Gerlach, D., Vaïda, P., et al. (2018). “Effectiveness of high-intensity jump training countermeasure on mitral and aortic flow after 58-days head-down bed-rest assessed by phase-contrast MRI,” in International Astronautical Congress, Vol. 2018 (Bremen), 135–141.
Caiani, E. G., Martin-Yebra, A., Landreani, F., Bolea, J., Laguna, P., and Vaïda, P. (2016). Weightlessness and cardiac rhythm disorders: current knowledge from space flight and bed-rest studies. Front. Astron. Space Sci. 3:27. doi: 10.3389/fspas.2016.00027
Caiani, E. G., Pellegrini, A., Bolea, J., Sotaquira, M., Almeida, R., and Vaïda, P. (2013). Impaired T-wave amplitude adaptation to heart-rate induced by cardiac deconditioning after 5-days of head-down bed-rest. Acta Astron. 91, 166–172. doi: 10.1016/j.actaastro.2013.05.016
Convertino, V. A. (2009). Status of cardiovascular issues related to space flight: implications for future research directions. Respir. Physiol. Neurobiol. 169, S34–S37. doi: 10.1016/j.resp.2009.04.010
D'Aunno, D. S., Dougherty, A. H., DeBlock, H. F., and Meck, J. V. (2003). Effect of short-and long-duration spaceflight on QTc intervals in healthy astronauts. Am. J. Cardiol. 91, 494–497. doi: 10.1016/s0002-9149(02)03259-9
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
Garrett-Bakelman, F. E., Darshi, M., Green, S. J., Gur, R. C., Lin, L., Macias, B. R., et al. (2019). The NASA Twins Study: a multidimensional analysis of a year-long human spaceflight. Science 364:eaau8650. doi: 10.1126/science.aau8650
Grenon, S. M., Xiao, X., Hurwitz, S., Ramsdell, C. D., Sheynberg, N., Kim, C., et al. (2005). Simulated microgravity induces microvolt T wave alternans. Ann. Noninvasive Electrocardiol. 10, 363–370. doi: 10.1111/j.1542-474X.2005.00654.x
Hanson, B., Child, N., Van Duijvenboden, S., Orini, M., Chen, Z., Coronel, R., et al. (2014). Oscillatory behavior of ventricular action potential duration in heart failure patients at respiratory rate and low frequency. Front. Physiol. 5:414. doi: 10.3389/fphys.2014.00414
Kardys, I., Kors, J. A., van der Meer, I. M., Hofman, A., van der Kuip, D. A., and Witteman, J. C. (2003). Spatial QRS-T angle predicts cardiac death in a general population. Eur. Heart J. 24, 1357–1364. doi: 10.1016/S0195-668X(03)00203-3
Konda, N. N., Karri, R. S., Winnard, A., Nasser, M., Evetts, S., Boudreau, E., et al. (2019). A comparison of exercise interventions from bed rest studies for the prevention of musculoskeletal loss. npj Microgravity 5:12. doi: 10.1038/s41526-019-0073-4
Kramer, A., Kümmel, J., Mulder, E., Gollhofer, A., Frings-Meuthen, P., and Gruber, M. (2017). High-intensity jump training is tolerated during 60 days of bed rest and is very effective in preserving leg power and lean body mass: an overview of the Cologne RSL Study. PLOS ONE 12:e0169793. doi: 10.1371/journal.pone.0169793
Kramer, A., Poppendieker, T., and Gruber, M. (2019). Suitability of jumps as a form of high-intensity interval training: effect of rest duration on oxygen uptake, heart rate and blood lactate. Eur. J. Appl. Physiol. 119, 1149–1156. doi: 10.1007/s00421-019-04105-w
Maggioni, M. A., Castiglioni, P., Merati, G., Brauns, K., Gunga, H.-C., Mendt, S., et al. (2018). High-intensity exercise mitigates cardiovascular deconditioning during long-duration bed rest. Front. Physiol. 9:1553. doi: 10.3389/fphys.2018.01553
Malik, M., Bigger, J. T., Camm, A. J., Kleiger, R. E., Malliani, A., Moss, A. J., et al. (1996). Heart rate variability: standards of measurement, physiological interpretation, and clinical use. Eur. Heart J. 17, 354–381.
Martinez, J. P., Almeida, R., Olmos, S., Rocha, A. P., and Laguna, P. (2004). A wavelet-based ECG delineator: evaluation on standard databases. IEEE Trans. Biomed. Eng. 51, 570–581. doi: 10.1109/TBME.2003.821031
Martín-Yebra, A., Caiani, E. G., Monasterio, V., Pellegrini, A., Laguna, P., and Martínez, J. P. (2015). Evaluation of T-wave alternans activity under stress conditions after 5 d and 21 d of sedentary head-down bed rest. Physiol. Meas. 36, 2041–2055. doi: 10.1088/0967-3334/36/10/2041
McRae, G., Payne, A., Zelt, J. G., Scribbans, T. D., Jung, M. E., Little, J. P., et al. (2012). Extremely low volume, whole-body aerobic–resistance training improves aerobic fitness and muscular endurance in females. Appl. Physiol. Nutr. Metab. 37, 1124–1131. doi: 10.1139/h2012-093
Meck, J. V., Reyes, C. J., Perez, S. A., Goldberger, A. L., and Ziegler, M. G. (2001). Marked exacerbation of orthostatic intolerance after long- vs. short-duration spaceflight in veteran astronauts. Psychosom. Med. 63, 865–873. doi: 10.1097/00006842-200111000-00003
Pagani, M., Lombardi, F., Guzzetti, S., Rimoldi, O., Furlan, R., Pizzinelli, P., et al. (1986). Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ. Res. 59, 178–193.
Pagani, M., Malfatto, G., Pierini, S., Casati, R., Masu, A. M., Poli, M., et al. (1988). Spectral analysis of heart rate variability in the assessment of autonomic diabetic neuropathy. J. Auton. Nervous Syst. 23, 143–153.
Pavy-Le Traon, A., Heer, M., Narici, M. V., Rittweger, J., and Vernikos, J. (2007). From space to Earth: advances in human physiology from 20 years of bed rest studies (1986–2006). Eur. J. Appl. Physiol. 101, 143–194. doi: 10.1007/s00421-007-0474-z
Piccirillo, G., Magrì, D., Ogawa, M., Song, J., Chong, V. J., Han, S., et al. (2009). Autonomic nervous system activity measured directly and QT interval variability in normal and pacing-induced tachycardia heart failure dogs. J. Am. Coll. Cardiol. 54, 840–850. doi: 10.1016/j.jacc.2009.06.008
Porter, B., van Duijvenboden, S., Bishop, M. J., Orini, M., Claridge, S., Gould, J., et al. (2018). Beat-to-beat variability of ventricular action potential duration oscillates at low frequency during sympathetic provocation in humans. Front. Physiol. 9:147. doi: 10.3389/fphys.2018.00147
Pueyo, E., Malik, M., and Laguna, P. (2008). A dynamic model to characterize beat-to-beat adaptation of repolarization to heart rate changes. Biomed. Signal Process. Control 3, 29–43. doi: 10.1016/j.bspc.2007.09.005
Pueyo, E., Orini, M., Rodríguez, J. F., and Taggart, P. (2016). Interactive effect of beta-adrenergic stimulation and mechanical stretch on low-frequency oscillations of ventricular action potential duration in humans. J. Mol. Cell. Cardiol. 97, 93–105. doi: 10.1016/j.yjmcc.2016.05.003
Pueyo, E., Smetana, P., Caminal, P., de Luna, A. B., Malik, M., and Laguna, P. (2004). Characterization of QT interval adaptation to RR interval changes and its use as a risk-stratifier of arrhythmic mortality in amiodarone-treated survivors of acute myocardial infarction. IEEE Trans. Biomed. Eng. 51, 1511–1520. doi: 10.1109/TBME.2004.828050
Ramírez, J., Orini, M., Mincholé, A., Monasterio, V., Cygankiewicz, I., Bayés de Luna, A., et al. (2017a). Sudden cardiac death and pump failure death prediction in chronic heart failure by combining ECG and clinical markers in an integrated risk model. PLOS ONE 12:e0186152. doi: 10.1371/journal.pone.0186152
Ramírez, J., Orini, M., Tucker, J. D., Pueyo, E., and Laguna, P. (2017b). Variability of ventricular repolarization dispersion quantified by time-warping the morphology of the T-waves. IEEE Trans. Biomed. Eng. 64, 1619–1630. doi: 10.1109/TBME.2016.2614899
Ramirez, R. J., Ajijola, O. A., Zhou, W., Holmström, B., Lüning, H., Laks, M. M., et al. (2011). A new electrocardiographic marker for sympathetic nerve stimulation: modulation of repolarization by stimulation of stellate ganglia. J. Electrocardiol. 44, 694–699. doi: 10.1016/j.jelectrocard.2011.07.030
Rizas, K. D., Hamm, W., Kääb, S., Schmidt, G., and Bauer, A. (2016). Periodic repolarisation dynamics: a natural probe of the ventricular response to sympathetic activation. Arrhythm. Electrophysiol. Rev. 5, 31–36. doi: 10.15420/aer.2015:30:2
Rizas, K. D., McNitt, S., Hamm, W., Massberg, S., Kääb, S., Zareba, W., et al. (2017). Prediction of sudden and non-sudden cardiac death in post-infarction patients with reduced left ventricular ejection fraction by periodic repolarization dynamics: MADIT-II substudy. Eur. Heart J. 38, 2110–2118. doi: 10.1093/eurheartj/ehx161
Rizas, K. D., Nieminen, T., Barthel, P., Zürn, C. S., Kähönen, M., Viik, J., et al. (2014). Sympathetic activity-associated periodic repolarization dynamics predict mortality following myocardial infarction. J. Clin. Invest. 124, 1770–1780. doi: 10.1172/JCI70085
Russomano, T., Baers, J. H., Velho, R., Cardoso, R. B., Ashcroft, A., Rehnberg, L., et al. (2013). A comparison between the 2010 and 2005 basic life support guidelines during simulated hypogravity and microgravity. Extreme Physiol. Med. 2:11. doi: 10.1186/2046-7648-2-11
Sakowski, C., Starc, V., Smith, S. M., and Schlegel, T. T. (2011). Sedentary long-duration head-down bed rest and ECG repolarization heterogeneity. Aviat. Space Environ. Med. 82, 416–423. doi: 10.3357/asem.2945.2011
Sampedro-Puente, D. A., Fernandez-Bes, J., Porter, B., van Duijvenboden, S., Taggart, P., and Pueyo, E. (2019). Mechanisms underlying interactions between low-frequency oscillations and beat-to-beat variability of celullar ventricular repolarization in response to sympathetic stimulation: implications for arrhythmogenesis. Front. Physiol. 10:916. doi: 10.3389/fphys.2019.00916
Schneider, S. M., Lee, S. M. C., Macias, B. R., Watenpaugh, D. E., and Hargens, A. R. (2009). WISE-2005: exercise and nutrition countermeasures for upright VO2pk during bed rest. Med. Sci. Sports Exerc. 41, 2165–2176. doi: 10.1249/mss.0b013e3181aa04e5
Tribulova, N., Szeiffova Bacova, B., Egan Benova, T., Knezl, V., Barancik, M., and Slezak, J. (2017). Omega-3 index and anti-arrhythmic potential of omega-3 PUFAs. Nutrients 9:1191. doi: 10.3390/nu9111191
Zwart, S. R., Pierson, D., Mehta, S., Gonda, S., and Smith, S. M. (2010). Capacity of omega-3 fatty acids or eicosapentaenoic acid to counteract weightlessness-induced bone loss by inhibiting nf-κb activation: from cells to bed rest to astronauts. J. Bone Miner. Res. 25, 1049–1057. doi: 10.1359/jbmr.091041
Keywords: microgravity, periodic repolarization dynamics (PRD), ventricular repolarization, autonomous nervous system, electrocardiogram (ECG) processing, tilt table test
Citation: Palacios S, Caiani EG, Landreani F, Martínez JP and Pueyo E (2019) Long-Term Microgravity Exposure Increases ECG Repolarization Instability Manifested by Low-Frequency Oscillations of T-Wave Vector. Front. Physiol. 10:1510. doi: 10.3389/fphys.2019.01510
Received: 26 August 2019; Accepted: 29 November 2019;
Published: 17 December 2019.
Edited by:George E. Billman, The Ohio State University, United States
Reviewed by:Ruben Coronel, University of Amsterdam, Netherlands
Emilio Vanoli, University of Pavia, Italy
Copyright © 2019 Palacios, Caiani, Landreani, Martínez and Pueyo. 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: Saúl Palacios, firstname.lastname@example.org
†These authors have contributed equally to this work as last authors