Quantitative analysis of the impact of respiratory state on the heartbeat-induced movements of the heart and its substructures

Abstract

Purpose

This study seeks to examine the influence of the heartbeat on the position, volume, and shape of the heart and its substructures during various breathing states. The findings of this study will serve as a valuable reference for dose-volume evaluation of the heart and its substructures in radiotherapy for treating thoracic tumors.

Methods

Twenty-three healthy volunteers were enrolled in this study, and cine four-dimensional magnetic resonance images were acquired during periods of end-inspiration breath holding (EIBH), end-expiration breath holding (EEBH), and deep end-inspiration breath holding (DIBH). The MR images were used to delineate the heart and its substructures, including the heart, pericardium, left ventricle (LV), left ventricular myocardium, right ventricle (RV), right ventricular myocardium (RVM), ventricular septum (VS), atrial septum (AS), proximal and middle portions of the left anterior descending branch (pmLAD), and proximal portion of the left circumflex coronary branch (pLCX). The changes in each structure with heartbeat were compared among different respiratory states.

Results

Compared with EIBH, EEBH increased the volume of the heart and its substructures by 0.25–3.66%, while the average Dice similarity coefficient (DSC) increased by − 0.25 to 8.7%; however, the differences were not statistically significant. Conversely, the VS decreased by 0.89 mm in the left–right (LR) direction, and the displacement of the RV in the anterior–posterior (AP) direction significantly decreased by 0.76 mm (p < 0.05). Compared with EIBH and EEBH, the average volume of the heart and its substructures decreased by 3.08–17.57% and 4.09–20.43%, respectively, during DIBH. Accordingly, statistically significant differences (p < 0.05) were observed in the volume of the heart, pericardium, LV, RV, RVM, and AS. The average DSC increased by 0–37.04% and − 2.6 to 32.14%, respectively, with statistically significant differences (p < 0.05) found in the right ventricular myocardium and interatrial septum. Furthermore, the displacements under DIBH decreased in the three directions (i.e.,− 1.73 to 3.47 mm and − 0.36 to 2.51 mm). In this regard, the AP displacement of the heart, LV, RV, RVM, LR direction, LV, RV, and AS showed statistically significant differences (p < 0.05). The Hausdorff distance (HD) of the heart and its substructures under the three breathing states are all greater than 11 mm.

Conclusion

The variations in the displacement and shape alterations of the heart and its substructures during cardiac motion under various respiratory states are significant. When assessing the dose-volume index of the heart and its substructures during radiotherapy for thoracic tumors, it is essential to account for the combined impacts of cardiac motion and respiration.

Introduction

The thoracic region is prone to developing malignancies such as lung, esophageal, and breast cancer, which can be treated with radiotherapy [1,2,3]. However, this poses a significant risk to patient survival time and quality of life, with radiation-induced heart damage being one of the most common complications [4,5,6,7]. Numerous studies have suggested that radiation-induced heart damage might offset the clinical benefits of radiotherapy, underlining the importance of early detection [8]. Current clinical methodology predicts such damage using dose-volume indices of the heart obtained from static simulation positioning 3DCT, which ignores the effects of breathing and cardiac motion on the heart's position, volume, and shape. Studies have indicated that the dose-volume indices differ significantly from the actual dose [9]. In past research, when conducting electrocardiogram-gated 4DCT scans under end-inspiration breath holding (EIBH), we observed significant cardiac shape and volume changes throughout the cycle, strongly affecting dosage evaluation accuracy [10]. When a patient's breathing state shifts, it changes the thoracic pressure, affecting the heart's range of motion and shape modulation. Unfortunately, few trials have sought to analyze the changes in the position, volume, and shape of the heart and its substructures under various respiratory states. Our study aims to examine these changes in three breathing conditions: EIBH, tidal end-expiratory breath holding (EEBH) and deep end-inspiration breath holding (DIBH).

Materials and methods

General study object information

The study included 23 participants, aged between 20 and 28 years, with a median age of 23 years, comprising 12 males and 11 females, all of whom demonstrated a breath-holding capacity of over 20 s during a single breath-hold in EIBH, EEBH and DIBH. Prior to participating in the study, all volunteers gave their informed consent, and the study was approved by the Ethics Committee of Shandong Cancer Hospital.

Acquisition of 4D MR images

All study participants underwent thoracic imaging with a GE 3.0 T superconducting MR scanner (Discovery 750 W, General Electric, USA) while maintaining EIBH, EEBH and DIBH, in which electrocardiogram gating was utilized to obtain 4D MR images. The scanning sequence used was a multislice 2D fast imaging employing steady-state acquisition (FIESTA). Scans had a slice thickness of 5 mm and a field of view measuring 50 cm and were performed without a gap. 4D MR images were reconstructed over 20 phases, ranging from 0 to 95% in 5% intervals, synchronized with the cardiac cycle.

Outlining of the heart and its substructures

Using MIM Maestro 7.1.5 (MIM, USA), the 4D MR images were manually segmented for different respiratory phases. In each image phase, the heart, pericardium, ventricular septum (VS), left and right ventricular myocardium (LVM and RVM, respectively), left ventricle (LV), right ventricle (RV), atrial septum (AS), proximal and middle portions of the left anterior descending coronary branch (pmLAD), and the proximal portion of the left circumflex coronary branch (pLCX) were segmented.

The upper boundary of the heart and pericardium is defined as the roof of the left atrium, and the lower boundary is the apex of the heart. The upper boundary of the left and right ventricles is the top of their respective chambers, and the lower boundary is the bottom of their chambers. The upper and lower boundaries of the left and right ventricular myocardium are the same as those of their corresponding chambers. The upper and lower boundaries of the interventricular septum are the same as those of the left ventricular myocardium. The upper boundary of the interatrial septum is the roof of the right atrium, and the lower boundary is the junction of the left atrium and left ventricle. The left anterior descending branch of the coronary artery starts at the bifurcation of the left main stem and ends at the bifurcation of the second diagonal branch. The left circumflex artery begins at the bifurcation of the left main stem and ends at the midpoint between the left margin of the heart and the posterior interventricular groove. All delineations were performed by two physicians jointly, and in case of any doubts, a third physician was involved in the discussion to determine the delineations.

Comparative analysis

We compared the displacement in the center of mass of the heart and its substructures in the left–right (LR), anterior–posterior (AP), and cranio-caudal (CC) directions during the cardiac cycle across the different respiratory states. Additionally, we compared the morphological and volumetric changes in the heart and its substructures. We calculated the volumetric variations by dividing the difference between the maximum and minimum volumes by the minimum volume. We evaluated the change in shape using the dice similarity coefficient (DSC).

$${\text{DSC}} = \frac{{2\left| {A \cap B} \right|}}{\left| A \right| + \left| B \right|}$$

here A refers to the volume of the region of interest (ROI) at the reference phase (0%), while B represents the volume of the ROI at other phases. The overlapping volume between the two is represented as |A ∩ B|.

The Hausdorff distance (HD) is a measure of the similarity between two sets of points. It is a definition of the distance between two point sets: Suppose there are two sets A = {a1,…, ap}, B = {b1,…,bq}; then, the HD between these two point sets is defined as:

$$HA,B = {\text{max}}\left( {\underbrace {{{\text{max}}}}_{a \in A}\left\{ {\underbrace {\min }_{b \in B}\left\| {a - b} \right\|} \right\},\underbrace {{{\text{max}}}}_{b \in B}\left\{ {\underbrace {{{\text{min}}}}_{a \in A}\left\| {b - a} \right\|} \right\}} \right)$$

HD is measured in millimeters (mm), and the higher the HD is, the greater the difference.

Statistical methods

SPSS 25 (IBM, Chicago, USA) software was used to analyze the data, and the results are presented as the mean ± standard deviation (${\overline{\text{x}}}$  ± s). We applied paired t tests to compare normally distributed data and Mann‒Whitney U tests to compare nonnormally distributed data. Significance was indicated if p < 0.05.

Results

Analysis of the movement of and morphological changes in the heart and its substructures during EIBH at end-expiration

During EIBH at end-expiration, the heart and its substructures were displaced by an average of 1.4–10.73 mm, 0.68–12.57 mm, and 2.31–13.18 mm in the LR, AP, and CC directions, respectively (Table 1). The volume variations ranged from 13.45 to 136.53% (Table 2), the DSC variations ranged from 6.75 to 2318.06% (Table 3), and the HD_MAX ranged from 11.8 to 27.28 mm (Table 4). The RVM exhibited the greatest displacement in the left–right direction (10.73 mm), while the pmLAD showed the greatest displacement in the AP and CC directions (12.57 mm and 13.18 mm, respectively), and the pericardium demonstrated the smallest displacement in all directions (1.4 mm, 0.68 mm, and 2.31 mm, respectively). The LV displayed the largest volume variations (136.53%), and the pericardium displayed the smallest (13.45%). The AS DSC exhibited the greatest change (2318.06%), and the pericardium DSC exhibited the least (6.75%). Statistically significant differences (p < 0.05) were observed between the maximum and minimum values of each structure's volume and DSCs (Tables 2 and 3). The LVM displayed the largest HDMAX value (27.28 mm), and the pericardium displayed the smallest HD_MAX value (11.8 mm) (Fig. 1).

Table 1 Comparison of the displacement of the centroid of the heart and its substructures under three respiratory states

Table 2 Comparison of volume changes in the heart and its substructures during the three respiratory states (X ± S, mm)

Table 3 Comparison of Dice similarity coefficient (DSC) changes in the heart and its substructures during the three respiratory states. (${\overline{\text{X}}}$ ± S, mm)

Table 4 Comparison of Hausdorff distance (HD) of the heart and its substructures under three respiratory states. (${\overline{\text{X}}}$ ± S, mm)

Fig. 1

Changes in the position and shape of the heart and its substructures at the end of a calm inhalation at the same level. The positions and shapes of the heart and its substructures are shown in the following time phases: A 0% phase, B 25% phase, C 50% phase, D 75% phase, and E the overall display of the heart and its substructures over the 20 time phases

Analysis of the movement and morphological changes of the heart and its substructures during EEBH at end-expiration

During EEBH at end-expiration, the heart and its substructures showed a mean displacement range of 1.19–11.69 mm, 0.64–11.53 mm, and 2.21–14.38 mm in the LR, AP, and CC directions, respectively (Table 1). The volume variations ranged from 12.74 to 123.75% (Table 2), and the DSC variations ranged from 6.9 to 1889.85% (Table 3). Except for pLCX, which exhibited the largest CC displacement at 14.38 mm, the RV volume variations were the greatest at 123.75%, which was different from the results of EIBH at end-expiration. The other structures remained unchanged regarding maximum and minimum displacement, volume variations, and DSC variations relative to the results of EIBH at end-expiration. Statistically significant differences (p < 0.05) were observed between the maximum and minimum values of each structure's volume and DSC (Tables 2 and 3). The LVM displayed the largest HDMAX value (24.32 mm), and the pericardium displayed the smallest HDMAX value (12.08 mm).

Compared to EIBH at end-expiration, the volume of the heart and its substructures increased by 0.25–3.66% under EEBH at end-expiration, but the differences were not statistically significant. The displacement of the heart and its substructures decreased in the LR direction by − 0.28 to 0.99 mm, in the AP direction by − 1.04 to 0.91 mm, and in the CC direction by − 1.37 to 2.52 mm. Only the displacement of the VS in the LR direction and the RV in the AP direction exhibited statistically significant differences (p < 0.05). The average DSC of the heart and its substructures increased by − 2.5 to 8.7%, but the differences were not statistically significant. HD_MAX has been reduced to a range of − 2.28 to 2.96 mm, with the maximum reduction being in the RVM (Fig. 2).

Fig. 2

Comparison of the position and shape of the heart and its substructures at the same level during the three respiratory states at the 0% phase. A Position and shape of the heart and its substructures at the end of calm exhalation, B position and shape of the heart and its substructures at the end of calm inhalation, C position and shape of the heart and its substructures at the end of deep inhalation, and D display of the position and shape of the heart and its substructures under the three respiratory states on an image obtained at the end of calm exhalation

Analysis of the movement and morphological changes of the heart and its substructures during DIBH at end-inhalation

During DIBH at end-inhalation, the heart and its substructures were displaced by an average of 0.94–8.22 mm, 0.86–11.72 mm, and 2.26–13.54 mm in the LR, AP, and CC directions, respectively (Table 1). The volume variations ranged from 14.07 to 131.31% (Table 2), the DSC variations ranged from 6.82 to 736.97% (Table 3), and the HD_MAX ranged from 11.53 to 21.91 mm (Table 4). The structures with the largest and smallest displacement, volume variations, and DSC variations remained unchanged from the results of EIBH at end-inhalation. Statistically significant differences (p < 0.05) were observed between the maximum and minimum values of each structure's volume and DSC (Tables 2 and 3). The RVM displayed the largest HDMAX value (21.91 mm), and the AS displayed the smallest HDMAX value (11.53 mm).

Compared to EIBH at end-inhalation, the volume of the heart and its substructures decreased by 3.08–17.57% during DIBH at end-inhalation. Statistically significant differences (p < 0.05) were observed in the heart, pericardium, LV, RV, RVM, and AS but not in the LVM and VS. Displacements in the LR direction decreased by − 0.16 to 3.47 mm, in the AP direction by − 0.22 to 2.56 mm, and in the CC direction by − 1.73 to 3.25 mm. Statistically significant differences (p < 0.05) were observed in the LR direction of the heart, pericardium, LV, RV, RVM, and VS and in the AP direction of the LV, RVM, and AS. The average DSC of the heart and its substructures increased by 0–37.04%, with statistically significant differences (p < 0.05) for the LVM, RVM, and AS, but not for the heart, pericardium, LV, RV, or VS. HD_MAX was reduced to a range of − 1.16 to 5.37 mm, with the maximum reduction being in the RVM (Fig. 2).

Relative to those during EEBH at end-expiration, the volume of the heart and its substructures decreased by 4.09–20.43% during DIBH at end-inhalation, with statistically significant differences (p < 0.05) in the heart, pericardium, LV, RV, RVM, VS, and AS, but not in LVM. Displacements in the LR direction decreased by 0.12–2.51 mm, in the AP direction by − 0.18 to 2.05 mm, and in the CC direction by − 0.36 to 0.73 mm. Statistically significant differences (p < 0.05) were observed in the LR direction of the heart, LV, RV, RVM, and VS, as well as in the AP direction of LV, RV, RVM, AS, and pLCX. In the CC direction, only pLCX exhibited statistically significant differences. The average DSC of the heart and its substructures increased by − 2.6 to 32.14%, with statistically significant differences (p < 0.05) observed only in the RVM and AS. HD_MAX has been reduced to a range of − 1.7 to 2.63 mm, with the maximum reduction being in the pLCX.

Figure 3 illustrates the volume changes of the heart and its substructures at different time points. There was good consistency in the volume changes in the heart, pericardium, LV, RV, and RVM. Additionally, there was good consistency in the volume changes in the LVM and the AS. Among the three breathing states, the volume of the heart and its substructures was at its lowest during DIBH.

Fig. 3

Volume changes in the heart and its substructures at different time phases during the three respiratory states

Discussion

This study is the only quantitative analysis of the changes in the position, volume, and shape of the heart and its substructures during three different breathing states, EIBH, EEBH, and DIBH, over the cardiac cycle.

Radiation-induced cardiac injury is a leading cause of nontumor-related death in patients with thoracic tumors receiving radiation therapy. According to Gagliardi et al. [11], adverse effects of radiation therapy on the heart are usually related to the radiation dose volume, but significant variations have been observed in the dose-volume indices proposed by various researchers, and there is no unified standard to date. This issue primarily stems from the fact that the dose-volume indices obtained from 3DCT simulation positioning under voluntary breathing disregard the influence of breathing and heartbeats, leading to substantial calculation errors for the dose-volume indices. Respiratory and cardiac motion can displace the heart, which can cause a decline in target coverage or an elevation of the dose on normal tissues, ultimately increasing the risk of treatment failure or complications [12,13,14]. Motion dose-volume index assessment research of the heart and its substructures is currently in its infancy.

Tong et al. [15] utilized calm inspiration electrocardiogram-gated 4D-CT to investigate volume variations in the heart, pericardium, and LVM during relaxation. The findings revealed volume variations of 16.49 ± 3.85%, 12.62 ± 3.94%, and 24.23 ± 11.35% in response to heartbeats. These findings are similar to the results of this study, in which the volume variations of the heart, pericardium, and LVM during calm inspiration at end-exhalation were 17.79 ± 3.79%, 13.45 ± 3.33%, and 28.14 ± 8.39%, respectively. Johnson et al. [16] examined 15 patients with coronary heart disease using digital subtraction X-ray angiography. The study found that the left anterior descending and left circumflex coronary arteries exhibited significant motion in the head-foot direction relative to the LR and front-back directions, which is comparable with the results of the present study. Li et al. [17] identified larger motion amplitudes in the AP and head-foot directions of the AS than those of the VS in their AS study. This finding aligns with the observed motion patterns of the VS and AS in response to heartbeats during the three breathing states evaluated in this research.

This study revealed a fundamental similarity in the volume and shape variations of the heart, pericardium, LV, RV, and RVM across all breathing states. Similarly, the volume and shape variations of the LVM and AS also displayed a basic consistency. However, the AS exhibited a distinct pattern of changes relative to the heart and pericardium. These findings suggest that relying on conventional dose-volume indices for the heart or pericardium to represent cardiac substructures can result in substantial errors.

The comparison of heart and substructure volumes revealed basic consistency in the volume and shape change patterns between EIBH and EEBH. The primary dissimilarity lies in the LR motion of the VS and the front-back motion of the RV. However, for patients with esophageal cancer or central lung cancer undergoing radiation dose evaluation, the effects of LR and front-back direction movements must be considered. Conversely, for breast cancer patients undergoing radiation therapy, increased front-back direction movements can cause an increase in the volume entering or approaching the irradiation range, resulting in a potentially significant dose escalation.

Compared with EIBH, DIBH led to a significant reduction in the volume of the heart and its substructures—except for the LVM—as well as a notable decrease in LR displacement of the heart, LV, RV, RVM, and VS. This effect primarily results from heightened thoracic pressure during DIBH, which causes compression of the heart and its substructures, resulting in decreased volume and LR displacement. Simultaneously, there is a significant reduction in the magnitude of the front-back motion of the LV, RVM, and AS, which can limit the entry of the heart and its substructures into the radiation range for therapies for thoracic wall tumors, breast cancer, and other cancers.

Observations of displacement changes in the heart and substructures during the three breathing states attest to the highly irregular motion present throughout the cardiac cycle. The most noteworthy motions of the heart, pericardium, and pLCX occur in the head-foot direction, whereas the LR direction displacement of the RV and RVM is greatest, and the front-back motion is the most pronounced for the LV and LVM. The motion patterns of the pmLAD and AS are similar in both the front-back and head-foot directions. The displacement of the VS is similar across all three directions.

The HD_MAX of the heart and its substructures are both greater than 11 mm. HD is very sensitive to abnormal changes. Breath holding at the end of deep inspiration can reduce the HD_MAX of cardiac substructures except for the pericardium and left ventricular cavity. The displacement and deformation of cardiac substructures vary in each direction. Overall, the deformation in the apical direction is greater than that in the base direction. The changes in the ventricles and myocardium are relatively large and are mainly related to heart function.

The participants in this study comprised a healthy population. We employed magnetic resonance imaging, a radiation-free scanning method that requires multiple repeated breath holdings over a lengthy duration, which demands a high level of physical fitness. Young and healthy participants are able to better control respiratory repeatability in different states, resulting in higher-quality MR images. By obtaining healthy cardiac images, a cardiac motion model was developed, which provides a reliable basis for analyzing the assessment errors in cardiac doses. In future studies, we plan to apply the established motion model of the heart and its substructures in different breathing states to patients with thoracic tumors, such as breast cancer, esophageal cancer, and lung cancer. Our goal is to examine the differences in dose-volume indices caused by different breathing states and to facilitate practical clinical applications.

Conclusion

Our quantitative analysis of the effects of cardiac motion during three breathing states on the motion state of the heart and its substructures indicated significant differences in the displacement and shape changes of these structures. When evaluating radiation doses for thoracic tumor radiotherapy, it is essential to consider the combined effects of cardiac motion and breathing on the heart and its substructures when evaluating dose-volume indices.