Low WSS and High OSI Measured by 3D Cine PC MRI Reflect High Pulmonary Artery Pressures in Suspected Secondary Pulmonary Arterial Hypertension

Reduced dynamic changes in pulmonary artery compliance during isometric handgrip exercise in patients with heart failure

Exercise intolerance is a debilitating symptom in heart failure (HF), adversely affecting both quality of life and long-term prognosis. Emerging evidence suggests that pulmonary artery (PA) compliance may be a contributing factor. This study aims to non-invasively assess PA compliance and its dynamic properties during isometric handgrip (HG) exercise in HF patients and healthy controls, using cardiovascular magnetic resonance (CMR). We prospectively enrolled 36 subjects, comprising 17 HF patients (NYHA class II and III) and 19 healthy controls. Participants performed an HG test, and we assessed changes in PA compliance and hemodynamic flow parameters using advanced CMR techniques. We also explored the relationship between CMR-derived PA compliance metrics and established clinical indicators, ensuring the validity of our findings through intra- and interobserver agreements. HF patients had significantly lower resting PA compliance compared to controls (28.9% vs. 50.1%, p < 0.01). During HG exercise, HF patients exhibited a dampened adaptability in PA compliance. Hemodynamic responses, including heart rate and blood pressure, were not significantly different between the groups. Further analyses revealed a significant correlation between changes in PA compliance and functional capacity, and an inverse relationship with NYHA class. Our study demonstrates a marked difference in PA vascular responses during HG exercise between HF patients and healthy controls. The compromised adaptability in PA compliance in HF patients is correlated with diminished functional capacity. These findings have significant clinical implications and may guide future interventional strategies in HF management.

Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstructs flow and increases pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the small pulmonary arteries due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we applied a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) approach to investigate the impact of CTEPH on pulmonary artery stiffening, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) in extralobar (main, right, and left) pulmonary arteries and intralobar (distal to the extralobar) arteries. Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar and intralobar arteries. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI, quantified by the novel parameter φ , which is related to thrombogenicity. Our analysis reveals a strong positive correlation between increases in mean pulmonary artery pressure (mPAP) and φ from baseline to CTEPH in individual subjects, which supports the suggestion that increased φ drives disease severity. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.

Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstruct flow and increase pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the pulmonary vasculature due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we developed a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) models to investigate the impact of CTEPH on pulmonary artery stiffening, time averaged wall shear stress (TAWSS), and oscillatory shear index (OSI). Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar (main, right and left pulmonary arteries) and intralobar vessels. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression, including the combination of low mean wall shear stress with high oscillation. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.

Clinical Application of 4D Flow MR Imaging to Pulmonary Hypertension

Pulmonary hypertension (PH) is characterized by elevated pulmonary arterial pressure (PAP). Although right-heart catheterization is the gold standard method for the diagnosis of PH by definition, various less-invasive imaging tests have been used for screening, detection of underlying diseases-causing PH, and monitoring of diseases. Among them, 4D flow MRI is an emerging and unique imaging test that allows for comprehensive visualization of blood flow in the right heart and proximal pulmonary arteries. The characteristic blood flow pattern observed in patients with PH is vortical flow formation in the main pulmonary artery. Recent studies have proposed the use of these findings to determine not only the presence of PH but also estimate the mean PAP. Other applications of 4D flow MRI for PH include measurement of wall shear stress, helicity, and 3D flow balance in the pulmonary arteries. It is worth noting that 4D flow has also the potential for longitudinal follow-ups. In this review, the clinical definition of PH, summary of conventional imaging tests, characteristics of pulmonary arterial flow as shown by 4D flow MRI, and clinical application of 4D flow MRI in the management of patients with PH will be discussed.

Technical Background for 4D Flow MR Imaging

Recently, the hemodynamic assessments with 3D cine phase-contrast (PC) MRI (4D flow MRI) have attracted considerable attention from clinicians. Unlike 2D cine PC MRI, the technique allows for cardiac phase-resolved data acquisitions of flow velocity vectors within the entire FOV during a clinically viable period. Thus, the method has enabled retrospective flowmetry in the spatial and temporal axes, which are essential to derive hemodynamic parameters related to vascular homeostasis and those to the progression of the pathologies. Accelerations in imaging are critical for this technology to be clinically viable; however, a high SNR or velocity-to-noise ratio (VNR) is also vital for accurate flow measurements. In this chapter, the technologies enabling this difficult balance are discussed.

Patient-Specific Computational Analysis of Hemodynamics in Adult Pulmonary Hypertension

Pulmonary hypertension (PH) is a progressive disease characterized by elevated pressure and vascular resistance in the pulmonary arteries. Nearly 250,000 hospitalizations occur annually in the US with PH as the primary or secondary condition. A definitive diagnosis of PH requires right heart catheterization (RHC) in addition to a chest computed tomography, a walking test, and others. While RHC is the gold standard for diagnosing PH, it is invasive and posseses inherent risks and contraindications. In this work, we characterized the patient-specific pulmonary hemodynamics in silico for diverse PH WHO groups. We grouped patients on the basis of mean pulmonary arterial pressure (mPAP) into three disease severity groups: at-risk ([Formula: see text], denoted with A), mild ([Formula: see text], denoted with M), and severe ([Formula: see text], denoted with S). The pulsatile flow hemodynamics was simulated by evaluating the three-dimensional Navier-Stokes system of equations using a flow solver developed by customizing OpenFOAM libraries (v5.0, The OpenFOAM Foundation). Quasi patient-specific boundary conditions were implemented using a Womersley inlet velocity profile and transient resistance outflow conditions. Hemodynamic indices such as spatially averaged wall shear stress ([Formula: see text]), wall shear stress gradient ([Formula: see text]), time-averaged wall shear stress ([Formula: see text]), oscillatory shear index ([Formula: see text]), and relative residence time ([Formula: see text]), were evaluated along with the clinical metrics pulmonary vascular resistance ([Formula: see text]), stroke volume ([Formula: see text]) and compliance ([Formula: see text]), to assess possible spatiotemporal correlations. We observed statistically significant decreases in [Formula: see text], [Formula: see text], and [Formula: see text], and increases in [Formula: see text] and [Formula: see text] with disease severity. [Formula: see text] was moderately correlated with [Formula: see text] and [Formula: see text] at the mid-notch stage of the cardiac cycle when these indices were computed using the global pulmonary arterial geometry. These results are promising in the context of a long-term goal of identifying computational biomarkers that can serve as surrogates for invasive diagnostic protocols of PH.

Computational simulation-derived hemodynamic and biomechanical properties of the pulmonary arterial tree early in the course of ventricular septal defects

Untreated ventricular septal defects (VSDs) can lead to pulmonary arterial hypertension (PAH) characterized by elevated pulmonary artery (PA) pressure and vascular remodeling, known as PAH associated with congenital heart disease (PAH-CHD). Though previous studies have investigated hemodynamic effects on vascular mechanobiology in late-stage PAH, hemodynamics leading to PAH-CHD initiation have not been fully quantified. We hypothesize that abnormal hemodynamics from left-to-right shunting in early stage VSDs affects PA biomechanical properties leading to PAH initiation. To model PA hemodynamics in healthy, small, moderate, and large VSD conditions prior to the onset of vascular remodeling, computational fluid dynamics simulations were performed using a 3D finite element model of a healthy 1-year-old's proximal PAs and a body-surface-area-scaled 0D distal PA tree. VSD conditions were modeled with increased pulmonary blood flow to represent degrees of left-to-right shunting. In the proximal PAs, pressure, flow, strain, and wall shear stress (WSS) increased with increasing VSD size; oscillatory shear index decreased with increasing VSD size in the larger PA vessels. WSS was higher in smaller diameter vessels and increased with VSD size, with the large VSD condition exhibiting WSS >100 dyn/cm[Formula: see text], well above values typically used to study dysfunctional mechanotransduction pathways in PAH. This study is the first to estimate hemodynamic and biomechanical metrics in the entire pediatric PA tree with VSD severity at the stage leading to PAH initiation and has implications for future studies assessing effects of abnormal mechanical stimuli on endothelial cells and vascular wall mechanics that occur during PAH-CHD initiation and progression.

Morphological and hemodynamic analysis of the patient-specific renal cell carcinoma models

Although the morbidity of renal cell carcinoma (RCC) has been increasing as the seventh most common tumours, to our knowledge, there is few studies foucsing on the hemodynamics of the renal artery (RA) with RCC. The objective of this study is to perform morphological and hemodynamic analysis of the RA and abdominal aorta artery (AAA) in the control healthy and RCC patient groups. Three-dimensional (3D) geometries are reconstructed from 18 control healthy subjects and 15 RCC patients based on Computed Tomography Angiography (CTA) images. There is higer in the lumen diameter of the RA (6.21 ± 0.89 mm) and curvature of the RA (1.2 ± 0.07) in the RCC patient group compared with the control healthy group (4.29 ± 1.08 mm, 1.1 ± 0.1), respectively. In the hemodynamic analysis, the surface area ratio (%) of low time-averaged wall shear stress (SAR-TAWSS) at the RA (10.65 ± 11.65) and AAA (48.49 ± 12.79) in the RCC patient group is significantly higher than that in the control healthy group (0.23 ± 0.22, 21.57 ± 20.5), respectively. It is found that RCC altered the morphology of the RA in the RCC patient group, which could deteriorate the hemodynamic environment of the RA and AAA. The finding in this study could enhance us to understand the progression of vascular disease caused by RCC.

Abnormal Flow Dynamics Result in Low Wall Shear Stress and High Oscillatory Shear Index in Abdominal Aortic Dilatation: Initial in vivo Assessment with 4D-flow MRI

PURPOSE

To characterize the non-laminar flow dynamics and resultant decreased wall shear stress (WSS) and high oscillatory shear index (OSI) of the infrarenal abdominal aortic dilatation, cardiac phase-resolved 3D phase-contrast MRI (4D-flow MRI) was performed.

METHODS

The prospective single-arm study was approved by the Institutional Review Board and included 18 subjects (median 67.5 years) with the dilated infrarenal aorta (median diameter 35 mm). 4D-flow MRI was conducted on a 1.5T MRI system. On 3D streamline images, laminar and non-laminar (i.e., vortex or helical) flow patterns were visually assessed both for the dilated aorta and for the undilated upstream aorta. Cardiac phase-resolved flow velocities, WSS and OSI, were also measured for the dilated aorta and the upstream undilated aorta.

RESULTS

Non-laminar flow represented by vortex or helical flow was more frequent and overt in the dilated aorta than in the undilated upstream aorta (P < 0.0156) with a very good interobserver agreement (weighted kappa: 0.82-1.0). The WSS was lower, and the OSI was higher on the dilated aortic wall compared with the proximal undilated segments. In mid-systole, mean spatially-averaged WSS was 0.20 ± 0.016 Pa for the dilated aorta vs. 0.68 ± 0.071 Pa for undilated upstream aorta (P < 0.0001), and OSI on the dilated aortic wall was 0.093 ± 0.010 vs. 0.041 ± 0.0089 (P = 0.013). The maximum values and the amplitudes of the WSS at the dilated aorta were inversely proportional to the ratio of dilated/undilated aortic diameter (r = -0.694, P = 0.0014).

CONCLUSION

4D-flow can characterize abnormal non-laminar flow dynamics within the dilated aorta in vivo. The wall of the infrarenal aortic dilatation is continuously and increasingly affected by atherogenic stimuli due to the flow disturbances represented by vortex or helical flow, which is reflected by lower WSS and higher OSI.

4D Flow when and how?

4D Flow is an emerging MR technique enabling three-dimensional and cardiac phase-resolved flowmetry with ECG-gated phase-contrast MRI that increased the speed of data acquisitions, accuracy and robustness. The method is promoting researches in areas that have not been fully addressed before in the cardiovascular system, such as flowmetry of the bloodstream across the valves, within the heart chambers, complexed flow dynamics such as vortex, helical or retrograde. Wall shear stress and other potential biomarkers derived from 4D Flow are known to be related to vascular wall diseases such as atherosclerosis. In this review, fundamental concepts of 4D Flow technique and post-processing, benefits and limitations as well as its clinical applications are discussed, and the importance of quality control and validation of the method is emphasized. New ideas inspired by 4D Flow can help clinicians and MR scientists further understand the role of flow dynamics in health sciences, diseases and various aspects of cardiovascular physiology.

A canonical correlation analysis of the relationship between clinical attributes and patient-specific hemodynamic indices in adult pulmonary hypertension

Pulmonary hypertension (PH) is a progressive disease affecting approximately 10-52 cases per million, with a higher incidence in women, and with a high mortality associated with right ventricle (RV) failure. In this work, we explore the relationship between hemodynamic indices, calculated from in silico models of the pulmonary circulation, and clinical attributes of RV workload and pathological traits. Thirty-four patient-specific pulmonary arterial tree geometries were reconstructed from computed tomography angiography images and used for volume meshing for subsequent computational fluid dynamics (CFD) simulations. Data obtained from the CFD simulations were post-processed resulting in hemodynamic indices representative of the blood flow dynamics. A retrospective review of medical records was performed to collect the clinical variables measured or calculated from standard hospital examinations. Statistical analyses and canonical correlation analysis (CCA) were performed for the clinical variables and hemodynamic indices. Systolic pulmonary artery pressure (sPAP), diastolic pulmonary artery pressure (dPAP), cardiac output (CO), and stroke volume (SV) were moderately correlated with spatially averaged wall shear stress (0.60 ≤ R² ≤ 0.66; p < 0.05). Similarly, the CCA revealed a linear and strong relationship (ρ = 0.87; p << 0.001) between 5 clinical variables and 2 hemodynamic indices. To this end, in silico models of PH blood flow dynamics have a high potential for predicting the relevant clinical attributes of PH if analyzed in a group-wise manner using CCA.

Assessment of the Hemodynamics of Autogenous Arteriovenous Fistulas With 4D Phase Contrast-Based Flow Quantification MRI in Dialysis Patients

BACKGROUND

Regular monitoring of autogenous arteriovenous fistulas (AVFs) for hemodialysis patients has importance. Hence, 4D flow MRI may be an alternative for assessing the hemodynamics of AVFs.

PURPOSE

To compare the hemodynamics of AVFs using Doppler ultrasound (DUS) and 4D-MRI in renal dialysis patients.

STUDY TYPE

Case-control study from October 2017 to April 2018.

POPULATION

Fifty patients (age [range] = 59.52 [39-71] years) with AVFs were included.

FIELD STRENGTH/SEQUENCE

Black-blood MRI and 4D flow MRI at 3.0T and AVF ultrasonography were also performed.

ASSESSMENT

The hemodynamics acquired from 4D flow MRI and ultrasonography by two radiologists were compared. The AVF anatomy was described through an examination of the black-blood MRI.

STATISTICAL TESTS

The consistency of AVF anatomy and hemodynamics and the consistency of the hemodynamics of AVFs from 4D flow MRI and ultrasound were analyzed by paired t-tests. The morphological parameters of AVFs acquired from black-blood MRI were used for a Pearson correlation analysis with the hemodynamic parameters obtained from 4D flow MRI data.

RESULTS

The consistency of the morphological and hemodynamic parameters measured from MRI by the two radiologists was good (all P < 0.01). The velocities and flow volumes from the 4D flow MRI and vascular ultrasound of AVFs were in moderate agreement (all P < 0.05, r = 0.292-0.569), except for the peak flow velocity at the anastomosis (P = 0.366, r = -0.078). The flow volume and WSS near the anastomotic site were closely related to the morphology of the AVFs (all P < 0.05). The hemodynamics of the complications group were significantly different from those of patients without any complications (normal patients group) (all P < 0.01).

DATA CONCLUSION

Compared with ultrasonography, 4D flow MRI is a promising technique to noninvasively estimate the AVF hemodynamics of renal dialysis patients.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:1272-1280.

Elevated diastolic wall shear stress in regurgitant semilunar valvular lesions

BACKGROUND

Alterations in wall shear stress (WSS) assessed using 4D flow MRI have been shown to play a role in various vascular pathologies, such as bicuspid aortic valve aortopathy. Most studies have focused on systolic WSS, whereas altered diastolic hemodynamics in regurgitant semilunar valvular lesions have not so far been well characterized.

PURPOSE

To investigate diastolic WSS in aortic and pulmonary regurgitation.

STUDY TYPE

Retrospective data analysis.

POPULATION

Thirty tetralogy of Fallot patients, 19 bicuspid aortic valve patients, 11 healthy volunteers.

FIELD STRENGTH/SEQUENCE

5 T, 3D time-resolved phase-contrast MRI with 3D velocity encoding.

ASSESSMENT

Estimation of WSS and its axial and circumferential vector components along cardiac cycle timeframes in the proximal main pulmonary artery in pulmonary regurgitation (PR) and in the proximal ascending aorta in aortic regurgitation (AR) as well as in healthy volunteers.

STATISTICAL TESTS

Wilcoxon matched pairs test was used for intra-group comparisons and Mann-Whitney test for intergroup comparisons. Correlations were assessed using Spearman correlation.

RESULTS

WSS along the entire cardiac cycle was higher in PR and AR in comparison with controls (mean WSS 0.381 ± 0.070 vs. 0.220 ± 0.018, P < 0.0001; 0.361 ± 0.099 vs. 0.212 ± 0.030, P < 0.0001; respectively). Peak diastolic WSS was significantly higher than the mean WSS in AR and PR (P < 0.0001-0.005). The severity of PR correlated with the peak diastolic axial WSS (Spearman's r _s  = 0.454, P = 0.018), whereas the severity of AR correlated with both peak systolic and diastolic tangential WSS (Spearman's r _s  = 0.458, P = 0.049; r _s  = 0.539, P = 0.017, respectively).

DATA CONCLUSION

Elevated diastolic WSS is a component of the altered flow hemodynamics in AR and PR. This may give more insight into the pathophysiologic role of WSS in vascular remodeling in AR and PR.

LEVEL OF EVIDENCE

4 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2019;50:763-770.

4D Flow MR Imaging of Ophthalmic Artery Flow in Patients with Internal Carotid Artery Stenosis

Background and Purpose:

To assess the clinical feasibility of time-resolved 3D phase contrast (4D Flow) MRI assessment of the ophthalmic artery (OphA) flow in patients with internal carotid artery stenosis (ICS).

Materials and Methods:

Twenty-one consecutive patients with unilateral ICS were recruited. 4D Flow MRI and acetazolamide-stress brain perfusion single photon emission computed tomography (SPECT) were performed. The flow direction on the affected-side OphA was categorized into native flow (anterograde or unclear) and non-native flow (retrograde flow) based on 4D Flow MRI. In the affected-side middle cerebral artery (MCA) territory, the ratio of rest cerebral blood flow to normal control (RCBF_MCA) and cerebral vascular reserve (CVR_MCA) were calculated from SPECT dataset. High-risk patients were defined based on the previous large cohort study (RCBF_MCA < 80% and CVR_MCA < 10%).

Results:

Eleven patients had native OphA flow (4 anterograde, 7 unclear) and the remaining 10 had non-native OphA flow. RCBF_MCA and CVR_MCA each were significantly lower in non-native flow group (84.9 ± 18.9% vs. 69.8 ± 7.3%, P < 0.05; 36.4 ± 20.6% vs. 17.0 ± 15.0%, P < 0.05). Four patients in the non-native flow group and none in the native flow group were confirmed as high-risk (Sensitivity/Specificity, 1.00/0.65).

Conclusion:

The 6 min standard 4D Flow MRI assessment of OphA in patients with ICS can predict intracranial hemodynamic impairment.