Study of impaired cardiovascular and respiratory coupling during orthostatic stress based on joint symbolic dynamics.

The present study investigates the instantaneous coupling among the cardiac, vascular, and respiratory systems, using the heart rate, respiration, and systolic and diastolic blood pressure variability in 12 healthy and 16 vasovagal syncope female subjects during a head-up tilt (HUT) testing protocol at 70° This study contributes to the coupling analysis by using a nonlinear joint symbolic dynamics (JSD) in a high-temporal resolution scheme, based on 5 min segments of the time series that are shifted every minute. For each segment, a bivariate JSD matrix was constructed to obtain global and local coupling indices in accordance to Shannon's entropy and the probability of occurrence of various bivariate words, respectively. The novel approach revealed important findings in the coupling dynamics of the systems, thus allowing the detection of group differences during the early orthostatic phase, and during the HUT test, before the occurrence of any pre-syncopal symptoms. In patients, the global indices indicated a significant decrease of cardiovascular coupling, starting at 10 min after the tilt-up, manifested by reduced baroreflex sensitivity and cardiorespiratory coupling that was initiated 8 min after the onset of the orthostatic phase (OP). A decreased autonomic control on cardiovascular-respiratory couplings was further evidenced by increased alterations of the JSD indices during the OP compared to the supine position in patients compared to controls. Furthermore, findings based on local indices demonstrated that female patients showed reductions and disengagements in cardiovascular (p < 0.001) and cardiorespiratory (p < 0.01) couplings, as early as the first 5 min and during the complete OP.

Sistema de Apoyo para la Cuantificación Semiautomática de Fibrosis Pulmonar en Imágenes de TAC / Support System for Semiautomatic Quantification of Pulmonary Fibrosis in CT Images

Resumen: En este trabajo se presenta un método para calcular los niveles de fibrosis pulmonar en imágenes de tomografía axial computarizada. Se utilizó un algoritmo de segmentación semiautomática basado en el método de Chan-Vese. El método mostró similitudes de forma cualitativa en la región de la fibrosis con respecto al experto clínico. Sin embargo es necesario validar los resultados con una base de datos mayor. El método propuesto aproxima un porcentaje de fibrosis de forma fácil para apoyar su implementación en la práctica clínica minimizando la subjetividad del experto médico y generando una estimación cuantitativa de la región de fibrosis.

Abstract: A method to estimate the pulmonary fibrosis in computed tomography (CT) imaging is presented. A semi-automatic segmentation algorithm based on the Chan-Vese method was used. The proposed method shows a similar fibrosis región with respect to clinical expert. However, the results need to be validated in a bigger data base. The proposed method approximates a fibrosis percentage that allows to achieve this procedure easily in order to support its implementation in the clinical practice minimizing the clinical expert subjectivity and generating a quantitative estimation of fibrosis region.

Men and women should be separately investigated in studies of orthostatic challenge due to different gender-related dynamics of autonomic response.

In studies of autonomic regulation during orthostatic challenges only a few nonlinear methods have been considered without investigating the effect of gender in young controls. Especially, the temporal development of the autonomic regulation has not yet been explicitly analyzed using short-term segments in supine position, transition and orthostatic phase (OP). In this study, nonlinear analysis of cardiovascular and respiratory time series was performed to investigate how nonlinear indices are dynamically changing with respect to gender during orthostatic challenges. The analysis was carried out using shifted short-term segments throughout a head-up tilt test in 24 healthy subjects, 12 men (26 ± 4 years) and 12 age-matched women (26 ± 5 years), at supine position and during OP at 70°. The nonlinear methods demonstrated statistical differences in the autonomic regulation between males and females. Orthostatic stress caused significantly decreased heart rate variability due to increased sympathetic activity mainly in men, already at the beginning and during the complete OP, revealed by (a) increased occurrence of specific word types with constant fluctuations as pW111 from symbolic dynamics, (b) augmented fractal correlation properties by the short-term index alpha1 from detrended fluctuation analysis, (c) increased slope indices (21ati and 31ati) from auto-transinformation and (d) augmented time irreversibility indices demonstrating more temporal asymmetries and nonlinear dynamics in men than in women. After tilt-up, both men and women increased their sympathetic activity but in a different way. Time-dependent gender differences during orthostatic challenge were shown directly between men and women or indirectly comparing baseline and different temporal stages of OP. The proposed dynamical study of autonomic regulation has the advantage of screening the fluctuations of the sympathetic and vagal activities that can be quantified by the temporal behavior of nonlinear indices. The findings in this paper strongly suggest the need for gender separation in studies of the dynamics of autonomic regulation during orthostatic challenge.

Orthostatic stress causes immediately increased blood pressure variability in women with vasovagal syncope.

The cardiovascular and respiratory autonomic nervous regulation has been studied mainly by hemodynamic responses during different physical stressors. In this study, dynamics of autonomic response to an orthostatic challenge was investigated by hemodynamic variables and by diverse linear and nonlinear indices calculated from time series of beat-to-beat intervals (BBI), respiratory cycle duration (RESP), systolic (SYS) and diastolic (DIA) blood pressure. This study included 16 young female patients (SYN) with vasovagal syncope and 12 age-matched female controls (CON). The subjects were enrolled in a head-up tilt (HUT) test, breathing normally, including 5min of baseline (BL, supine position) and 18min of 70° orthostatic phase (OP). To increase the time resolution of the analysis the time series were segmented in five-minute overlapping windows with a shift of 1min. Hemodynamic parameters did not show any statistical differences between SYN and CON. Time domain linear analysis revealed increased respiratory frequency and increased blood pressure variability (BPV) in patients during OP meaning increased sympathetic activity and vagal withdrawal. Frequency domain analysis confirmed a predominance of sympathetic tone by steadily increased values of low over high frequency power in BBI and of low frequency power in SYS and DIA in patients during OP. The nonlinear analysis by symbolic dynamics seemed to be highly suitable for differentiation of SYN and CON in the early beginning of OP, i.e., 5min after tilt-up. In particular the index SYS_plvar3 showed less patterns of low variability in patients reflecting a steadily increase in both BPV and sympathetic activity. The proposed dynamical analysis could lead to a better understanding of the temporal underlying mechanisms in healthy subjects and patients under orthostatic stress.

Assessment of time-frequency representation techniques for thoracic sounds analysis.

A step forward in the knowledge about the underlying physiological phenomena of thoracic sounds requires a reliable estimate of their time-frequency behavior that overcomes the disadvantages of the conventional spectrogram. A more detailed time-frequency representation could lead to a better feature extraction for diseases classification and stratification purposes, among others. In this respect, the aim of this study was to look for an omnibus technique to obtain the time-frequency representation (TFR) of thoracic sounds by comparing generic goodness-of-fit criteria in different simulated thoracic sounds scenarios. The performance of ten TFRs for heart, normal tracheal and adventitious lung sounds was assessed using time-frequency patterns obtained by mathematical functions of the thoracic sounds. To find the best TFR performance measures, such as the 2D local (ρ(mean)) and global (ρ) central correlation, the normalized root-mean-square error (NRMSE), the cross-correlation coefficient (ρ(IF)) and the time-frequency resolution (res(TF)) were used. Simulation results pointed out that the Hilbert-Huang spectrum (HHS) had a superior performance as compared with other techniques and then, it can be considered as a reliable TFR for thoracic sounds. Furthermore, the goodness of HHS was assessed using noisy simulated signals. Additionally, HHS was applied to first and second heart sounds taken from a young healthy male subject, to tracheal sound from a middle-age healthy male subject, and to abnormal lung sounds acquired from a male patient with diffuse interstitial pneumonia. It is expected that the results of this research could be used to obtain a better signature of thoracic sounds for pattern recognition purpose, among other tasks.

Linear and nonlinear analysis of base lung sound in extrinsic allergic alveolitis patients in comparison to healthy subjects.

OBJECTIVE: Pulmonary disorders are frequently characterized by the presence of adventitious sounds added to the breathing or base lung sound (BLS). The aim of this work was to assess the features of BLS in extrinsic allergic alveolitis (EAA) patients in comparison to healthy subjects, applying linear and nonlinear analysis techniques. METHODS: We investigated the multichannel lung sounds on the posterior chest of 16 females, 8 healthy and 8 EAA patients, when breathing at 1.5 L/s. BLS linear features were obtained from the power spectral density (PSD) while nonlinear features were extracted by the concepts of irregularity and complexity, i.e., spectral, sample and multiscale entropy. RESULTS: The results demonstrated that spectral percentiles of BLS were lower in EAA patients than in healthy subjects but statistical significance (p<0.05) was obtained only for expiration at the left apical and both basal regions. Also, the maximum amplitude of the PSD in patients reached statistical significance ( p < 0.05) for the expiratory phase at basal regions. In the case of nonlinear techniques, significant lower values ( p < 0.05) were obtained for EAA patients during both respiratory phases at left apical and both basal regions. CONCLUSION: In conclusion, we found that BLS in chronic EAA patients is characterized by lower spectral percentiles, lower irregularity and lower complexity than in healthy subjects suggesting the feasibility of its clinical usefulness by screening its temporal alteration.

Effects of fetal respiratory movements on the short-term fractal properties of heart rate variability.

We evaluated the effect of fetal respiratory movements (RM) on the heart rate (HR) fractal dynamics.Abdominal ECG recordings were collected from low-middle-risk pregnant woman at rest. Mean gestational age was 34.8 ± 3.7 weeks. Ultrasound images were simultaneously acquired determining if RM were exhibited by fetuses. 13 pairs of HR series were compared. Each pair included 5 min of data from the same fetus either during the manifestation of RM or when there was no persistent indication of them. Detrended fluctuation analysis was applied to these series for obtaining the scaling exponent α1. HR series were also assessed using the conventional parameters RMSSD and HF power.The main findings of this contribution were the lack of significant changes in the scaling exponent α1 of fetal HR fluctuations as a result of RM. By contrast, HF power and RMSSD did show significant changes associated with the manifestation of fetal RM (p < 0.001 and p < 0.05, respectively). Yet the scaling exponent was the only parameter showing a significant relationship with the particular frequency of fetal RM (r s  = 0.6, p < 0.03). Given the invariability of α1 regarding the manifestation of fetal RM, we consider that the HR short-term fractal properties are convenient for assessing the cardiovascular prenatal regulation.

Assessment of ICA algorithms for the analysis of crackles sounds.

Blind source separation by independent component analysis has been applied extensively in the biomedical field for extracting different contributing sources in a signal. Regarding lung sounds analysis to isolate the adventitious sounds from normal breathing sound is relevant. In this work the performance of FastICA, Infomax, JADE and TDSEP algorithms was assessed using different scenarios including simulated fine and coarse crackles embedded in recorded normal breathing sounds. Our results pointed out that Infomax obtained the minimum Amari index (0.10037) and the maximum signal to interference ratio (1.4578e+009). Afterwards, Infomax was applied to 25 channels of recorded normal breathing sound where simulated fine and coarse crackles were added including acoustic propagation effects. A robust blind crackle separation could improve previous results in generating an adventitious acoustic thoracic imaging.

Fractal and nonlinear changes in the long-term baseline fluctuations of fetal heart rate.

The interpretation of heart rate patterns obtained by fetal monitoring relies on the definition of a baseline, which is considered as the running average heart rate in the absence of external stimuli during periods of fetal rest. We present a study along gestation of the baseline's fluctuations, in relation to fractal and nonlinear properties, to assess these fluctuations according with time-varying attracting levels introduced by maturing regulatory mechanisms. A low-risk pregnancy was studied weekly from the 17th to 38th week of gestation during long-term recording sessions at night (>6 h). Fetal averaged pulse rate samples and corresponding baseline series were obtained from raw abdominal ECG ambulatory data. The fractal properties of these series were evaluated by applying detrended fluctuation analysis. The baseline series were also explored to evaluate nonlinear properties and time ordering by applying the scaling magnitude and sign analyses. Our main findings are that the baseline shows fractal and even nonlinear anticorrelated fluctuations. This condition was specially the case before mid-gestation, as revealed by α values near to unit, yet becoming significantly more complex after 30 weeks of gestation as indicated by α(mag) values >0.5. The structured (i.e. not random) fluctuations and particular nonlinear changes that we found thus suggest that the baseline provides on itself information concerning the functional integration of cardiac regulatory mechanisms.

The autonomic condition of children with congenital hypothyroidism as indicated by the analysis of heart rate variability.

Early treatment by thyroid hormone in congenital hypothyroidism (CH) prevents for mental retardation, but the autonomic condition could well remain impaired. The aim of this study was to get insights about the cardiovascular regulation of CH children by analysing their heart rate variability (HRV) data. RR-interval data of 34 early treated CH children, 24 to 48 months of age, were collected at rest in a clinical environment. In all these cases, CH resulted from an abnormal development of thyroid gland. Data from 29 age matched-controls were also collected in a scholar environment. Short time series were then analysed in time and frequency domains. Detrended fluctuation analysis (DFA) was also applied to obtain the scaling exponent α(1). We found that the HRV obtained from children with CH is diminished as suggested by statistical differences (p<0.05) in the following parameters: RMSDD, pNN20, total power, LF and HF. However, no significant differences were found in the LF/HF ratio, the central frequencies of the LF and HF bands, as well as in α(1). Given that the significant differences were basically reflected in static parameters rather than in the frequency distribution or scale-invariant behaviour of the studied data, we consider that the analysis points toward the effect of treatment on rising heart rate by the thyroid hormone, and not necessarily toward any residual autonomic dysfunction, as the dominant underlying cause of the diminished HRV exhibited by CH children.

Assessment of multichannel lung sounds parameterization for two-class classification in interstitial lung disease patients.

This work deals with the assessment of different parameterization techniques for lung sounds (LS) acquired on the whole posterior thoracic surface for normal versus abnormal LS classification. Besides the conventional technique of power spectral density (PSD), the eigenvalues of the covariance matrix and both the univariate autoregressive (UAR) and the multivariate autoregressive models (MAR) were applied for constructing feature vectors as input to a supervised neural network (SNN). The results showed the effectiveness of the UAR modeling for multichannel LS parameterization, using new data, with classification accuracy of 75% and 93% for healthy subjects and patients, respectively.

Short-term heart rate dynamics of pregnant women.

Aiming to detect the stage of gestation where dynamical changes of the RR fluctuations may occur, we assessed short-term fluctuations of low risk pregnant women. Ninety six, 10min ECG recordings were collected along gestation (7 to 39 weeks). Corresponding RR fluctuations series were analysed to obtain the RMSSD, α(1), α(1(mag)) and α(1(sign)) parameters. Four groups covering first, second and last trimesters of gestation were conformed. No significant changes in α(1), which was close to unit, and α(1(sign)) among gestational groups were identified. But, in accordance with previous findings, we did find a significant reduction of RMSSD along gestation, and significant short-term changes that indicate a higher degree of nonlinearity after about 26 weeks of gestation (α(1(mag))>0.5)). These results suggest that the short-term heart rate dynamics of low risk pregnant women do not become compromised during gestation, despite the increased haemodynamic demands and other ongoing adaptations. Yet the complexity of the mechanisms involved in the cardiac regulation of pregnant women does seem to increase from mid-pregnancy, possibly owing to new short-term control influences or to modifications regardless the strength of the regulatory interactions.

Short-term heart rate dynamics of women during labor.

We studied 10 minutes segments of heartbeat interval fluctuations from 18 young women in labor with normal outcome of pregnancy. Data of each studied case were classified into two distinct groups. One group involving segments where the uterine activity was observable (three or more contractions in ten minutes), and the other group of reference having segments with fewer uterine activity or not presenting contractions at all. For comparison, we also included segments collected during the last trimester of gestation prior to labor from a third group of women. Corresponding RR interval series were analyzed to estimate RR(mean), RMSSD, α(1), α(1(MAG)) and α(1(SIGN)) parameters. No significant differences among groups were identified in RMSSD, α(1) and α(1(MAG)) Nevertheless, α(1(SIGN)) did present significant differences in comparison with the last trimester results (p<0.007), revealing a subtle change in the temporal organization of maternal RR series during labor. Results of these parameters then suggest that during labor, despite preserving a concomitant non-linear influence, the maternal short-term autonomic cardiac regulation behaves with less antagonism.

Imaging the thoracic distribution of normal breath sounds.

OBJECTIVE: To assess the feasibility to generate a confident image of normal breath sounds (BS) based on the quantitative analysis of multichannel sensors and imaging them in three known clinical classes, i.e., tracheal, bronchial and vesicular, identifying their spatial distribution with high resolution on the posterior thoracic surface. METHODS: Three parametrization techniques, the percentile frequencies, the univariate AR modeling, and the eigenvalues of the covariance matrix were evaluated when applied to BS. These sounds were acquired in twelve healthy subjects by a 5x5 sensor array on the posterior thoracic surface plus the sound at the tracheal position, to obtain feature vectors that fed a supervised multilayer neural network. Based on BS classification rate, the spatial distribution of each BS class was obtained by constructing an image using deterministic interpolation. RESULTS: The univariate AR modeling was the best parametrization technique producing a classification performance of 96% during the validation phase and just 4% of not classified feature vectors. Corresponding values for the percentile frequencies were 92% and 7.7%, whereas for the eigenvalues were 91% and 9.0%. CONCLUSION: This work shows that it is possible to generate confident images associated with the distribution of normal BS classes. Therefore, a detailed image about the spatial distribution of BS in humans might be helpful for detecting lung diseases.

Repeatability of heart rate variability in congenital hypothyroidism as analysed by detrended fluctuation analysis.

The analysis of heart rate fluctuations, or heart rate variability (HRV), may be applied to explore children's neurodevelopment. However, previous studies have reported poor reliability (repeatability) of HRV measures in children at rest and during light exercise. Whether the reliability can be improved by controlling variables such as physical activity, breathing rate and tidal volume, or by selecting non-conventional techniques for analysing the data remains as an open question. We evaluated the short-term repeatability of RR-interval data from medicated children with congenital hypothyroidism (CH). The alpha(1) exponents, obtained by detrended fluctuation analysis (DFA), from the data of 21 children collected at two different sessions were compared. Elapsed days between sessions were 59 +/- 33, and data were obtained during 10 min, trying to restrict the children's activity while being seated. We found statistical agreement between the means of alpha(1) exponents for each session (p = 0.94) and no bias with a low-coefficient variation (9.1%); an intraclass correlation coefficient ri = 0.48 ([0.14 0.72], 95% confidence interval) was also estimated. These findings, which were compared with results obtained by conventional time and frequency techniques, indicate the existence of agreement between the alpha(1) exponents obtained at each session, thereby providing support concerning the repeatability of HRV data as analysed by DFA in children with congenital hypothyroidism. Of particular interest was also the agreement found by using the central frequency of the high-frequency band and the parameter pNN20, both showing better or similar ri than alpha(1) (0.77 [0.57 0.89] and 0.51 [0.17 0.74], respectively), yet considerably better repeatability than other conventional time and frequency parameters.

Applying fractal analysis to short sets of heart rate variability data.

The aim of this study was to explore the interchangeability of fractal scaling exponents derived from short- and long-term recordings of real and synthetic data. We compared the alpha(1) exponents as obtained by detrended fluctuation analysis from RR-interval series (9 am to 6 pm) of 54 adults in normal sinus rhythm, and the alpha(1) estimated from shorted segments of these series involving only 50, 100, 200 and 300 RR intervals. Three series of synthetic data were also analysed. The principal finding of this study is the lack of individual agreement between alpha(1) derived from long and short segments of HRV data as indicated by the existence of bias and low intraclass correlation coefficient (r(i) = 0.158). The extent of variation in the estimation of alpha(1) from real data does not only appear related to segments' length, but also to different dynamics among subjects or lack of uniform scaling behaviour. However, we did find statistical agreement between the means of alpha(1) exponents from long and short segments, even for segments involving just 50 RR intervals. According to results of synthetic series, the 95% confidence interval found for the variation of alpha(1) using segments with 300 samples is [-0.1783 + 0.1828]. Caution should be taken concerning the use of short segments to obtain representative exponents of fractal RR dynamics; a circumstance not fully considered in several studies.

Time-Frequency Representations for second heart sound analysis.

Several researches have tried to provide a means to analyze the second heart sound (S2) in an attempt to understand the functional mechanisms in its genesis and for diagnosis purposes. In this work we tested Time-Frequency Representation (TFR) for simulated S2 selecting and applying classical and modern TFRs such as the Spectrogram, the Wigner-Ville Distribution, the Time Varying Autoregressive (TVAR) model, the Scalogram, and the Hilbert-Huang Spectrum (HHS) by Empirical Mode Decomposition. Two performance measures are proposed, the first one based on local 2D correlations (rho) between the ideal and the estimated TFRs images, while the second one based on time moments of the TFR images to provide the normalized root-mean-square error (NRMSE). Under no noise conditions, the TFRs by HHS and the TVAR modeling, by the Burg algorithm, resulted in a rho(average) of 0.788 and 0.812, and NRMSE of 0.172 and 0.195, respectively. Therefore, based on the lowest NRMSE, HHS was considered the TFR with the best performance. Afterward, HHS was applied to real S2 acquired at the aortic and pulmonary focal points.

Analysis of discontinuous adventitious lung sounds by Hilbert-Huang spectrum.

It is now widely accepted that crackles are associated with different pulmonary pathologies and different efforts have been done to detect and to extract them. Consequently, due to the difficulty for their characterization, the selection of an adequate time-frequency representation (TFR) for the analysis of their time-frequency dynamics is important. Traditionally, normal and abnormal lung sounds have been analyzed by the Spectrogram (SP). However, this analysis tool has certain disadvantages when one deals with nonstationary signals. As an effort to point out the appropriate analysis tool for crackles, this paper shows the performance of the Hilbert-Huang spectrum (HHS) for the analysis of fine and coarse crackles, simulated and real ones. The HHS allowed to analyze the evolving time-frequency of crackle sounds straightforward with good resolution compared with SP. Beside this enhanced time-frequency course, HHS could be useful to establish a signature to detect and separate fine from coarse crackles in order to screen pathologies and their progress during medication.

Respiratory acoustic thoracic imaging (RATHI): assessing intrasubject variability.

Respiratory acoustic thoracic imaging (RATHI) permits analysing lung sounds (LS) temporal and spatial distribution, however, a deep understanding of RATHI repeatability associated with the pulmonary function is necessary. As a consequence, in the current work intrasubject variability of RATHI is evaluated at different airflows. For generating RATHIs, LS were acquired at the posterior thoracic surface. The associated image was computed at the inspiratory phases by interpolation through a Hermite function. The acoustic information of eleven subjects was considered at airflows of 1.0, 1.5 and 2.0 L/s. Several RATHIs were generated for each subject according to the number of acquired inspiratory phases. Quadratic mutual information based on Cauchy-Schwartz inequality (I(CS)) was used to evaluate the degree of similitude between intrasubject RATHIs. The results indicated that, for the same subject, I(CS) averaged 0.893, 0.897, and 0.902, for airflows of 1.0, 1.5, and 2 L/s, respectively. In addition, when the airflow was increased, increments in intensity values and in the dispersion of the spatial distribution reflected in RATHI were observed. In conclusion, since the intrasubject variability of RATHI was low for airflows between 1.0 and 2.0 L/s, the pattern of sound distribution during airflow variations is repeatable but differences in sound intensity should be considered.

Asymmetry in lung sound intensities detected by respiratory acoustic thoracic imaging (RATHI) and clinical pulmonary auscultation.

RATHI was introduced as an attempt to further improve the association between anatomical zones and specific breathing activity, both spatially and temporally. This work compares RATHI with clinical pulmonary auscultation (PA) to assess the concordance between both procedures to detect asymmetries in lung sound (LS) intensities. Twelve healthy young males participated in the study and were auscultated by two experts. RATHI consisted in the acquisition of acoustical signals with an array of 5x5 sensors, while experts auscultated and described the intensity of LS heard using the same stethoscope on each sensor's position within the array. Comparisons were established looking for intensity asymmetries between apical vs. basal pulmonary regions and right vs. left hemithorax. By RATHI, most of the subjects showed asymmetries between apical and basal regions higher than 20%, whereas between left and right hemithorax asymmetries higher than 20% occurred only half of the time. RATHI and PA agreed 83 to 100% when apical to base acoustical information was compared, but when left to right asymmetries were considered these figures were about 40 to 50%. We concluded that RATHI has advantages as it gave more detailed and measurable information on LS than clinicians, who could not detect intensity asymmetries mainly below 20%.