A Framework for Assessing the Effect of Cardiac and Respiratory Motion for Stereotactic Arrhythmia Radioablation Using a Digital Phantom With a 17-Segment Model: A STOPSTORM.eu Consortium Study.

PURPOSE: The optimal motion management strategy for patients receiving stereotactic arrhythmia radioablation (STAR) for the treatment of ventricular tachycardia (VT) is not fully known. We developed a framework using a digital phantom to simulate cardiorespiratory motion in combination with different motion management strategies to gain insight into the effect of cardiorespiratory motion on STAR. METHODS AND MATERIALS: The 4-dimensional (4D) extended cardiac-torso (XCAT) phantom was expanded with the 17-segment left ventricular (LV) model, which allowed placement of STAR targets in standardized ventricular regions. Cardiac- and respiratory-binned 4D computed tomography (CT) scans were simulated for free-breathing, reduced free-breathing, respiratory-gating, and breath-hold scenarios. Respiratory motion of the heart was set to population-averaged values of patients with VT: 6, 2, and 1 mm in the superior-inferior, posterior-anterior, and left-right direction, respectively. Cardiac contraction was adjusted by reducing LV ejection fraction to 35%. Target displacement was evaluated for all segments using envelopes encompassing the cardiorespiratory motion. Envelopes incorporating only the diastole plus respiratory motion were created to simulate the scenario where cardiac motion is not fully captured on 4D respiratory CT scans used for radiation therapy planning. RESULTS: The average volume of the 17 segments was 6 cm3 (1-9 cm3). Cardiac contraction-relaxation resulted in maximum segment (centroid) motion of 4, 6, and 3.5 mm in the superior-inferior, posterior-anterior, and left-right direction, respectively. Cardiac contraction-relaxation resulted in a motion envelope increase of 49% (24%-79%) compared with individual segment volumes, whereas envelopes increased by 126% (79%-167%) if respiratory motion also was considered. Envelopes incorporating only the diastole and respiration motion covered on average 68% to 75% of the motion envelope. CONCLUSIONS: The developed LV-segmental XCAT framework showed that free-wall regions display the most cardiorespiratory displacement. Our framework supports the optimization of STAR by evaluating the effect of (cardio)respiratory motion and motion management strategies for patients with VT.

Stereotactic Arrhythmia Radioablation (STAR): Assessment of cardiac and respiratory heart motion in ventricular tachycardia patients - A STOPSTORM.eu consortium review.

AIM: To identify the optimal STereotactic Arrhythmia Radioablation (STAR) strategy for individual patients, cardiorespiratory motion of the target volume in combination with different treatment methodologies needs to be evaluated. However, an authoritative overview of the amount of cardiorespiratory motion in ventricular tachycardia (VT) patients is missing. METHODS: In this STOPSTORM consortium study, we performed a literature review to gain insight into cardiorespiratory motion of target volumes for STAR. Motion data and target volumes were extracted and summarized. RESULTS: Out of the 232 studies screened, 56 provided data on cardiorespiratory motion, of which 8 provided motion amplitudes in VT patients (n = 94) and 10 described (cardiac/cardiorespiratory) internal target volumes (ITVs) obtained in VT patients (n = 59). Average cardiac motion of target volumes was < 5 mm in all directions, with maximum values of 8.0, 5.2 and 6.5 mm in Superior-Inferior (SI), Left-Right (LR), Anterior-Posterior (AP) direction, respectively. Cardiorespiratory motion of cardiac (sub)structures showed average motion between 5-8 mm in the SI direction, whereas, LR and AP motions were comparable to the cardiac motion of the target volumes. Cardiorespiratory ITVs were on average 120-284% of the gross target volume. Healthy subjects showed average cardiorespiratory motion of 10-17 mm in SI and 2.4-7 mm in the AP direction. CONCLUSION: This review suggests that despite growing numbers of patients being treated, detailed data on cardiorespiratory motion for STAR is still limited. Moreover, data comparison between studies is difficult due to inconsistency in parameters reported. Cardiorespiratory motion is highly patient-specific even under motion-compensation techniques. Therefore, individual motion management strategies during imaging, planning, and treatment for STAR are highly recommended.

Optic nerve sheath diameter assessment by neurosonology: A review of methodologic discrepancies.

BACKGROUND AND PURPOSE: Reported cutoff values of the optic nerve sheath diameter (ONSD) for the diagnosis of elevated intracranial pressure (ICP) are inconsistent. This hampers ONSD as a possible noninvasive bedside monitoring tool for ICP. Because the influence of methodological differences on variations in cutoff values is unknown, we performed a narrative review to identify discrepancies in ONSD assessment methodologies and to investigate their effect on reported ONSD values. METHODS: We used a structured and quantitative approach in which each ONSD methodology found in the reviewed articles was categorized based on the characteristic appearance of the ultrasound images and ultrasound marker placement. Subsequently, we investigated the influence of the different methodologies on ONSD values by organizing the ONSDs with respect to these categories. RESULTS: In a total of 63 eligible articles, we could determine the applied ONSD assessment methodology. Reported ultrasound images either showed the optic nerve and its sheath as a dark region with hyperechoic striped band at its edges or as a single dark region surrounded by lighter retrobulbar fat. Four different ultrasound marker positions were used to delineate the optic nerve sheath, which resulted in different ONSD values and more importantly, different sensitivities to changes in ICP. CONCLUSIONS: Based on our observations, we recommend to place ultrasound markers at the outer edges of the hyperechoic striped bands or at the transitions from the single dark region to the hyperechoic retrobulbar fat because these locations yielded the highest sensitivity of ONSD measurements for increased ICP.

An Automated Algorithm for Optic Nerve Sheath Diameter Assessment from B-mode Ultrasound Images.

BACKGROUND AND PURPOSE: The optic nerve sheath diameter (ONSD) is a promising surrogate marker for the detection of raised intracranial pressure (ICP). However, inconsistencies in manual ONSD assessment are thought to affect ONSD and the corresponding ONSD cutoff values for the diagnosis of elevated ICP, hereby hampering the full potential of ONSD. In this study, we developed an image intensity-invariant algorithm to automatically estimate ONSD from B-mode ultrasound images at multiple depths. METHODS: The outcomes of the algorithm were validated against manual ONSD measurements by two human experts. Each expert analyzed the images twice (M1 and M2) in unknown order. RESULTS: The algorithm proved capable of segmenting the ONSD in 39 of 42 images, hereby showing mean differences of -.08 ± .45 and -.05 ± .41 mm compared to averaged ONSD values (M1 + M2/2) of Operator 1 and Operator 2, respectively, whereas the mean difference between the two experts was .03 ± .26 mm. Moreover, differences between algorithm-derived and expert-derived ONSD values were found to be much smaller than the 1 mm difference that is expected between patients with normal and elevated ICP, making it likely that our algorithm can distinguish between these patient groups. CONCLUSIONS: Our algorithm has the potential to improve the accuracy of ONSD as a surrogate marker for elevated ICP because it has no intrinsic variability. However, future research should be performed to validate if the algorithm does indeed result in more accurate noninvasive ICP predictions.

The Role of One-Dimensional Model-Generated Inter-Subject Variations in Systemic Properties on Wall Shear Stress Indices of Intracranial Aneurysms.

Variations in systemic properties of the arterial tree, such as aging-induced vessel stiffness, can alter the shape of pressure and flow waveforms. As a consequence, the hemodynamics around a cerebral aneurysm change, and therefore, also the corresponding in- and outlet boundary conditions (BCs) used for three-dimensional (3D) calculations of hemodynamic indices. In this study, we investigate the effects of variations in systemic properties on wall shear stress (WSS) indices of a cerebral aneurysm. We created a virtual patient database by varying systemic properties within physiological ranges. BCs for 3D-CFD simulations were derived using a pulse wave propagation model for each realization of the virtual database. WSS indices were derived from the 3D simulations and their variabilities quantified. Variations in BCs, caused by changes in systemic properties, yielded variabilities in the WSS indices that were of the same order of magnitude as differences in these WSS indices between ruptured and unruptured aneurysms. Sensitivity analysis showed that the systemic properties impacted both in- and outlet BCs simultaneously and altered the WSS indices. We conclude that the influence of variations in patient-specific systemic properties on WSS indices should be evaluated when using WSS indices in multidisciplinary rupture prediction models.

Intima heterogeneity in stress assessment of atherosclerotic plaques.

Atherosclerotic plaque rupture is recognized as the primary cause of cardiac and cerebral ischaemic events. High structural plaque stresses have been shown to strongly correlate with plaque rupture. Plaque stresses can be computed with finite-element (FE) models. Current FE models employ homogeneous material properties for the heterogeneous atherosclerotic intima. This study aimed to evaluate the influence of intima heterogeneity on plaque stress computations. Two-dimensional FE models with homogeneous and heterogeneous intima were constructed from histological images of atherosclerotic human coronaries (n = 12). For homogeneous models, a single stiffness value was employed for the entire intima. For heterogeneous models, the intima was subdivided into four clusters based on the histological information and different stiffness values were assigned to the clusters. To cover the reported local intima stiffness range, 100 cluster stiffness combinations were simulated. Peak cap stresses (PCSs) from the homogeneous and heterogeneous models were analysed and compared. By using a global variance-based sensitivity analysis, the influence of the cluster stiffnesses on the PCS variation in the heterogeneous intima models was determined. Per plaque, the median PCS values of the heterogeneous models ranged from 27 to 160 kPa, and the PCS range varied between 43 and 218 kPa. On average, the homogeneous model PCS values differed from the median PCS values of heterogeneous models by 14%. A positive correlation (R2 = 0.72) was found between the homogeneous model PCS and the PCS range of the heterogeneous models. Sensitivity analysis showed that the highest main sensitivity index per plaque ranged from 0.26 to 0.83, and the average was 0.47. Intima heterogeneity resulted in substantial changes in PCS, warranting stress analyses with heterogeneous intima properties for plaque-specific, high accuracy stress assessment. Yet, computations with homogeneous intima assumption are still valuable to perform sensitivity analyses or parametric studies for testing the effect of plaque geometry on PCS. Moreover, homogeneous intima models can help identify low PCS, stable type plaques with thick caps. Yet, for thin cap plaques, accurate stiffness measurements of the clusters in the cap and stress analysis with heterogeneous cap properties are required to characterize the plaque stability.

Biomechanical changes during abdominal aortic aneurysm growth.

The biomechanics-based Abdominal Aortic Aneurysm (AAA) rupture risk assessment has gained considerable scientific and clinical momentum. However, such studies have mainly focused on information at a single time point, and little is known about how AAA properties change over time. Consequently, the present study explored how geometry, wall stress-related and blood flow-related biomechanical properties change during AAA expansion. Four patients with a total of 23 Computed Tomography-Angiography (CT-A) scans at different time points were analyzed. At each time point, patient-specific properties were extracted from (i) the reconstructed geometry, (ii) the computed wall stress at Mean Arterial Pressure (MAP), and (iii) the computed blood flow velocity at standardized inflow and outflow conditions. Testing correlations between these parameters identified several nonintuitive dependencies. Most interestingly, the Peak Wall Rupture Index (PWRI) and the maximum Wall Shear Stress (WSS) independently predicted AAA volume growth. Similarly, Intra-luminal Thrombus (ILT) volume growth depended on both the maximum WSS and the ILT volume itself. In addition, ILT volume, ILT volume growth, and maximum ILT layer thickness correlated with PWRI as well as AAA volume growth. Consequently, a large ILT volume as well as fast increase of ILT volume over time may be a risk factor for AAA rupture. However, tailored clinical studies would be required to test this hypothesis and to clarify whether monitoring ILT development has any clinical benefit.