A quantitative high resolution voxel-wise assessment of myocardial blood flow from contrast-enhanced first-pass magnetic resonance perfusion imaging: microsphere validation in a magnetic resonance compatible free beating explanted pig heart model

High-resolution free-breathing automated quantitative myocardial perfusion by cardiovascular magnetic resonance for the detection of functionally significant coronary artery disease

AIMS

Current assessment of myocardial ischaemia from stress perfusion cardiovascular magnetic resonance (SP-CMR) largely relies on visual interpretation. This study investigated the use of high-resolution free-breathing SP-CMR with automated quantitative mapping in the diagnosis of coronary artery disease (CAD). Diagnostic performance was evaluated against invasive coronary angiography (ICA) with fractional flow reserve (FFR) measurement.

METHODS AND RESULTS

Seven hundred and three patients were recruited for SP-CMR using the research sequence at 3 Tesla. Of those receiving ICA within 6 months, 80 patients had either FFR measurement or identification of a chronic total occlusion (CTO) with inducible perfusion defects seen on SP-CMR. Myocardial blood flow (MBF) maps were automatically generated in-line on the scanner following image acquisition at hyperaemic stress and rest, allowing myocardial perfusion reserve (MPR) calculation. Seventy-five coronary vessels assessed by FFR and 28 vessels with CTO were evaluated at both segmental and coronary territory level. Coronary territory stress MBF and MPR were reduced in FFR-positive (≤0.80) regions [median stress MBF: 1.74 (0.90-2.17) mL/min/g; MPR: 1.67 (1.10-1.89)] compared with FFR-negative regions [stress MBF: 2.50 (2.15-2.95) mL/min/g; MPR 2.35 (2.06-2.54) P < 0.001 for both]. Stress MBF ≤ 1.94 mL/min/g and MPR ≤ 1.97 accurately detected FFR-positive CAD on a per-vessel basis (area under the curve: 0.85 and 0.96, respectively; P < 0.001 for both).

CONCLUSION

A novel scanner-integrated high-resolution free-breathing SP-CMR sequence with automated in-line perfusion mapping is presented which accurately detects functionally significant CAD.

Adenosine Triphosphate Stress Myocardial Strain in Ischemic Heart Disease: An Animal Study with Histological Validation

RATIONALE AND OBJECTIVES

It is still challenging for cardiac magnetic resonance (CMR) to detect ischemic heart disease (IHD) without the use of gadolinium contrast. We aimed to evaluate the potential value of adenosine triphosphate (ATP) stress myocardial strain derived from feature tracking (FT) as a novel method for detecting IHD in a swine model.

MATERIALS AND METHODS

CMR cines, myocardial perfusion imaging at rest and during ATP stress, and late gadolinium enhancement were obtained in both control and IHD swine. Normal, remote, ischemic, and infarcted myocardium were analyzed. The diagnostic accuracy of myocardial strain for infarction and ischemia was assessed using coronary angiography and pathology as reference.

RESULTS

Eleven IHD swine and five healthy control swine were enrolled in this study. Strain parameters, even at rest, were associated with myocardial ischemia and infarction(all p < 0.05). The area under receiver operating characteristic curve (AUC) values of all strain parameters for detecting infarcted myocardium exceeded 0.900 (all p < 0.05). The AUC values for detecting ischemic myocardium were as follows: 0.906 and 0.847 for stress and rest radial strain, 0.763 and 0.716 for stress and rest circumferential strain, 0.758 and 0.663 for stress and rest longitudinal strain (all p < 0.001). Heat maps demonstrated that all strain parameters showed mild to moderate correlations with the stress myocardial blood flow and myocardial perfusion reserve (all p < 0.05).

CONCLUSION

CMR-FT-derived ATP stress myocardial strain shows promise as a noninvasive method for detecting myocardial ischemia and infarction in an IHD swine model, with rest strain parameters offering potential as a needle-free diagnostic option.

Development of a cardiovascular magnetic resonance-compatible large animal isolated heart model for direct comparison of beating and arrested hearts

Cardiac motion results in image artefacts and quantification errors in many cardiovascular magnetic resonance (CMR) techniques, including microstructural assessment using diffusion tensor cardiovascular magnetic resonance (DT-CMR). Here, we develop a CMR-compatible isolated perfused porcine heart model that allows comparison of data obtained in beating and arrested states. Ten porcine hearts (8/10 for protocol optimisation) were harvested using a donor heart retrieval protocol and transported to the remote CMR facility. Langendorff perfusion in a 3D-printed chamber and perfusion circuit re-established contraction. Hearts were imaged using cine, parametric mapping and STEAM DT-CMR at cardiac phases with the minimum and maximum wall thickness. High potassium and lithium perfusates were then used to arrest the heart in a slack and contracted state, respectively. Imaging was repeated in both arrested states. After imaging, tissue was removed for subsequent histology in a location matched to the DT-CMR data using fiducial markers. Regular sustained contraction was successfully established in six out of 10 hearts, including the final five hearts. Imaging was performed in four hearts and one underwent the full protocol, including colocalised histology. The image quality was good and there was good agreement between DT-CMR data in equivalent beating and arrested states. Despite the use of autologous blood and dextran within the perfusate, T2 mapping results, DT-CMR measures and an increase in mass were consistent with development of myocardial oedema, resulting in failure to achieve a true diastolic-like state. A contiguous stack of 313 5-μm histological sections at and a 100-μm thick section showing cell morphology on 3D fluorescent confocal microscopy colocalised to DT-CMR data were obtained. A CMR-compatible isolated perfused beating heart setup for large animal hearts allows direct comparisons of beating and arrested heart data with subsequent colocalised histology, without the need for onsite preclinical facilities.

Effects of High Fat Versus Normal Diet on Extracellular Vesicle-Induced Angiogenesis in a Swine Model of Chronic Myocardial Ischemia

Background Mesenchymal stem cell-derived extracellular vesicles (EVs) promote angiogenesis in the ischemic myocardium. This study examines the difference in vascular density, myocardial perfusion, molecular signaling, and gene expression between normal diet (ND) and high fat diet (HFD) groups at baseline and following intramyocardial injection of EVs. Methods and Results Intact male Yorkshire swine fed either an ND (n=17) or HFD (n=14) underwent placement of an ameroid constrictor on the left circumflex coronary artery. Subsequently, animals received either intramyocardial injection of vehicle-saline as controls; (ND-controls n=7, HFD-controls, n=6) or EVs; (ND-EVs n=10, HFD-EVs n=8) into the ischemic territory. Five weeks later, myocardial function, perfusion, vascular density, cell signaling, and gene expression were examined. EVs improved indices of myocardial contractile function, myocardial perfusion, and arteriogenesis in both dietary cohorts. Interestingly, quantification of alpha smooth muscle actin demonstrated higher basal arteriolar density in HFD swine compared with their ND counterparts; whereas EVs were associated with increased CD31-labeled endothelial cell density only in the ND tissue, which approached significance. Levels of total endothelial nitric oxide synthase, FOXO1 (forkhead box protein O1) , transforming growth factor-β, phosphorylated VEGFR2 (vascular endothelial growth factor receptor 2), and phosphorylated MAPK ERK1/ERK2 (mitogen-activated protein kinase) were higher in ischemic myocardial lysates from ND-controls compared with HFD-controls. Conversely, HFD-control tissue showed increased expression of phosphorylated endothelial nitric oxide synthase, phosphorylated FOXO1, VEGFR2, and MAPK ERK1/ERK2 with respect to ND-controls. Preliminary gene expression studies indicate differential modulation of transcriptional activity by EVs between the 2 dietary cohorts. Conclusions HFD produces a profound metabolic disorder that dysregulates the molecular mechanisms of collateral vessel formation in the ischemic myocardium, which may hinder the therapeutic angiogenic effects of EVs.

3D-Printed Coronary Implants Are Effective for Percutaneous Creation of Swine Models with Focal Coronary Stenosis

Reliable, closed-chest methods for creating large animal models of acute myocardial hypoperfusion are limited. We demonstrated the feasibility and efficacy of using magnetic resonance (MR)-compatible 3D-printed coronary implants for establishing swine models of myocardial hypoperfusion. We designed, manufactured, and percutaneously deployed implants in 13 swine to selectively create focal coronary stenosis. To test the efficacy of the implants to cause hypoperfusion or ischemia in the perfused territory, we evaluated regional wall motion, myocardial perfusion, and infarction using MR imaging. The overall swine survival rate was 85% (11 of 13). The implant retrieval rate was 92% (12 of 13). Fluoroscopic angiography confirmed focal stenosis. Cine and perfusion MRI showed regional wall motion abnormalities and inducible ischemia, respectively. Late gadolinium enhancement and histopathology showed no myocardial infarction. Our minimally invasive technique has promising applications for validation of new diagnostic methods in cardiac MR. Graphical abstract Our new minimally invasive, percutaneous method for creating swine models of acute focal coronary stenosis can be used for magnetic resonance imaging studies of myocardial ischemia. Comparable to existing methods in its efficacy and reliability, this rapid prototyping technique will allow researchers to more easily conduct translational cardiac imaging studies of coronary artery disease in large animal models.

High-Resolution Cardiac Magnetic Resonance Imaging Techniques for the Identification of Coronary Microvascular Dysfunction

OBJECTIVES

This study assessed the ability to identify coronary microvascular dysfunction (CMD) in patients with angina and nonobstructive coronary artery disease (NOCAD) using high-resolution cardiac magnetic resonance (CMR) and hypothesized that quantitative perfusion techniques would have greater accuracy than visual analysis.

BACKGROUND

Half of all patients with angina are found to have NOCAD, while the presence of CMD portends greater morbidity and mortality, it now represents a modifiable therapeutic target. Diagnosis currently requires invasive assessment of coronary blood flow during angiography. With greater reliance on computed tomography coronary angiography as a first-line tool to investigate angina, noninvasive tests for diagnosing CMD warrant validation.

METHODS

Consecutive patients with angina and NOCAD were enrolled. Intracoronary pressure and flow measurements were acquired during rest and vasodilator-mediated hyperemia. CMR (3-T) was performed and analyzed by visual and quantitative techniques, including calculation of myocardial blood flow (MBF) during hyperemia (stress MBF), transmural myocardial perfusion reserve (MPR: MBF_HYPEREMIA / MBF_REST), and subendocardial MPR (MPR_ENDO). CMD was defined dichotomously as an invasive coronary flow reserve <2.5, with CMR readers blinded to this classification.

RESULTS

A total of 75 patients were enrolled (57 ± 10 years of age, 81% women). Among the quantitative perfusion indices, MPR_ENDO and MPR had the highest accuracy (area under the curve [AUC]: 0.90 and 0.88) with high sensitivity and specificity, respectively, both superior to visual assessment (both p < 0.001). Visual assessment identified CMD with 58% accuracy (41% sensitivity and 83% specificity). Quantitative stress MBF performed similarly to visual analysis (AUC: 0.64 vs. 0.60; p = 0.69).

CONCLUSIONS

High-resolution CMR has good accuracy at detecting CMD but only when analyzed quantitatively. Although omission of rest imaging and stress-only protocols make for quicker scans, this is at the cost of accuracy compared with integrating rest and stress perfusion. Quantitative perfusion CMR has an increasingly important role in the management of patients frequently encountered with angina and NOCAD.

Optimal Use of Vasodilators for Diagnosis of Microvascular Angina in the Cardiac Catheterization Laboratory

Among patients with angina and nonobstructive coronary artery disease, those with coronary microvascular dysfunction have a poor outcome. Coronary microvascular dysfunction is usually diagnosed by assessing flow reserve with an endothelium-independent vasodilator like adenosine, but the optimal diagnostic threshold is unclear. Furthermore, the incremental value of testing endothelial function has never been assessed before. We sought to determine what pharmacological thresholds correspond to exercise pathophysiology and myocardial ischemia in patients with coronary microvascular dysfunction.

Deep-Learning-Based Preprocessing for Quantitative Myocardial Perfusion MRI

BACKGROUND

Quantitative myocardial perfusion cardiac MRI can provide a fast and robust assessment of myocardial perfusion status for the noninvasive diagnosis of myocardial ischemia while being more objective than visual assessment. However, it currently has limited use in clinical practice due to the challenging postprocessing required, particularly the segmentation.

PURPOSE

To evaluate the efficacy of an automated deep learning (DL) pipeline for image processing prior to quantitative analysis.

STUDY TYPE

Retrospective.

POPULATION

In all, 175 (350 MRI scans; 1050 image series) clinical patients under both rest and stress conditions (135/10/30 training/validation/test).

FIELD STRENGTH/SEQUENCE

3.0T/2D multislice saturation recovery T1 -weighted gradient echo sequence.

ASSESSMENT

Accuracy was assessed, as compared to the manual operator, through the mean square error of the distance between landmarks and the Dice similarity coefficient of the segmentation and bounding box detection. Quantitative perfusion maps obtained using the automated DL-based processing were compared to the results obtained with the manually processed images.

STATISTICAL TESTS

Bland-Altman plots and intraclass correlation coefficient (ICC) were used to assess the myocardial blood flow (MBF) obtained using the automated DL pipeline, as compared to values obtained by a manual operator.

RESULTS

The mean (SD) error in the detection of the time of peak signal enhancement in the left ventricle was 1.49 (1.4) timeframes. The mean (SD) Dice similarity coefficients for the bounding box and myocardial segmentation were 0.93 (0.03) and 0.80 (0.06), respectively. The mean (SD) error in the RV insertion point was 2.8 (1.8) mm. The Bland-Altman plots showed a bias of 2.6% of the mean MBF between the automated and manually processed MBF values on a per-myocardial segment basis. The ICC was 0.89, 95% confidence interval = [0.87, 0.90].

DATA CONCLUSION

We showed high accuracy, compared to manual processing, for the DL-based processing of myocardial perfusion data leading to quantitative values that are similar to those achieved with manual processing.

LEVEL OF EVIDENCE

3 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:1689-1696.

Pixel-wise assessment of cardiovascular magnetic resonance first-pass perfusion using a cardiac phantom mimicking transmural myocardial perfusion gradients

PURPOSE

Cardiovascular magnetic resonance first-pass perfusion for the pixel-wise detection of coronary artery disease is rapidly becoming the clinical standard, yet no widely available method exists for its assessment and validation. This study introduces a novel phantom capable of generating spatially dependent flow values to enable assessment of new perfusion imaging methods at the pixel level.

METHODS

A synthetic multicapillary myocardial phantom mimicking transmural myocardial perfusion gradients was designed and manufactured with high-precision 3D printing. The phantom was used in a stationary flow setup providing reference myocardial perfusion rates and was scanned on a 3T system. Repeated first-pass perfusion MRI for physiological perfusion rates between 1 and 4 mL/g/min was performed using a clinical dual-sequence technique. Fermi function-constrained deconvolution was used to estimate pixel-wise perfusion rate maps. Phase contrast (PC)-MRI was used to obtain velocity measurements that were converted to perfusion rates for validation of reference values and cross-method comparison. The accuracy of pixel-wise maps was assessed against simulated reference maps.

RESULTS

PC-MRI indicated excellent reproducibility in perfusion rate (coefficient of variation [CoV] 2.4-3.5%) and correlation with reference values (R² = 0.985) across the full physiological range. Similar results were found for first-pass perfusion MRI (CoV 3.7-6.2%, R² = 0.987). Pixel-wise maps indicated a transmural perfusion difference of 28.8-33.7% for PC-MRI and 23.8-37.7% for first-pass perfusion, matching the reference values (30.2-31.4%).

CONCLUSION

The unique transmural perfusion pattern in the phantom allows effective pixel-wise assessment of first-pass perfusion acquisition protocols and quantification algorithms before their introduction into routine clinical use.

Coronary Microvascular Dysfunction Is Associated With Myocardial Ischemia and Abnormal Coronary Perfusion During Exercise

BACKGROUND

Coronary microvascular dysfunction (MVD) is defined by impaired flow augmentation in response to a pharmacological vasodilator in the presence of nonobstructive coronary artery disease. It is unknown whether diminished coronary vasodilator response correlates with abnormal exercise physiology or inducible myocardial ischemia.

METHODS

Patients with angina and nonobstructive coronary artery disease had simultaneous coronary pressure and flow velocity measured using a dual sensor-tipped guidewire during rest, supine bicycle exercise, and adenosine-mediated hyperemia. Microvascular resistance (MR) was calculated as coronary pressure divided by flow velocity. Wave intensity analysis quantified the proportion of accelerating wave energy (perfusion efficiency). Global myocardial blood flow and subendocardial:subepicardial perfusion ratio were quantified using 3-Tesla cardiac magnetic resonance imaging during hyperemia and rest; inducible ischemia was defined as hyperemic subendocardial:subepicardial perfusion ratio <1.0. Patients were classified as having MVD if coronary flow reserve <2.5 and controls if coronary flow reserve ≥2.5, with researchers blinded to the classification.

RESULTS

Eighty-five patients were enrolled (78% female, 57±10 years), 45 (53%) were classified as having MVD. Of the MVD group, 82% had inducible ischemia compared with 22% of controls (P<0.001); global myocardial perfusion reserve was 2.01±0.41 and 2.68±0.49 (P<0.001). In controls, coronary perfusion efficiency improved from rest to exercise and was unchanged during hyperemia (59±11% vs 65±14% vs 57±18%; P=0.02 and P=0.14). In contrast, perfusion efficiency decreased during both forms of stress in MVD (61±12 vs 44±10 vs 42±11%; both P<0.001). Among patients with a coronary flow reserve <2.5, 62% had functional MVD, with normal minimal MR (hyperemic MR<2.5 mmHg/cm/s), and 38% had structural MVD with elevated hyperemic MR. Resting MR was lower in those with functional MVD (4.2±1.0 mmHg/cm/s) than in those with structural MVD (6.9±1.7 mmHg/cm/s) or controls (7.3±2.2 mmHg/cm/s; both P<0.001). During exercise, the structural group had a higher systolic blood pressure (188±25 mmHg) than did those with functional MVD (161±27 mmHg; P=0.004) and controls (156±30 mmHg; P<0.001). Functional and structural MVD had similar stress myocardial perfusion and exercise perfusion efficiency values.

CONCLUSION

In patients with angina and nonobstructive coronary artery disease, diminished coronary flow reserve characterizes a cohort with inducible ischemia and a maladaptive physiological response to exercise. We have identified 2 endotypes of MVD with distinctive systemic vascular responses to exercise; whether endotypes have a different prognosis or require different treatments merits further investigation.

Danish study of Non-Invasive testing in Coronary Artery Disease 2 (Dan-NICAD 2): Study design for a controlled study of diagnostic accuracy

BACKGROUND

Coronary computed tomography angiography (CTA) is the preferred primary diagnostic modality when examining patients with low to intermediate pre-test probability of coronary artery disease (CAD). Only 20-30% of these have potentially obstructive CAD. Because of the relatively poor positive predictive value of coronary CTA, unnecessary invasive coronary angiographies (ICAs) are conducted with the costs and risks associated with the procedure. Hence, an optimized diagnostic CAD algorithm may reduce the numbers of ICAs not followed by revascularization. The Dan-NICAD 2 study has 3 equivalent main aims: (1) To examine the diagnostic precision of a sound-based diagnostic algorithm, The CADScor®System (Acarix A/S, Denmark), in patients with a low to intermediate pre-test risk of CAD referred to a primary examination by coronary CTA. We hypothesize that the CADScor®System provides better stratification prior to coronary CTA than clinical risk stratification scores alone. (2) To compare the diagnostic accuracy of 3T cardiac magnetic resonance imaging (3T CMRI), ⁸²rubidium positron emission tomography (⁸²Rb-PET), and CT-derived fractional flow reserve (FFR_CT) in patients where obstructive CAD cannot be ruled out by coronary CTA using ICA fractional flow reserve (FFR) as reference standard. (3) To compare the diagnostic performance of quantitative flow ratio (QFR) and ICA-FFR in patients with low to intermediate pre-test probability of CAD using ⁸²Rb-PET as reference standard.

METHODS

Dan-NICAD 2 is a prospective, multicenter, cross-sectional study including approximately 2,000 patients with low to intermediate pre-test probability of CAD and without previous history of CAD. Patients are referred to coronary CTA because of symptoms suggestive of CAD, as evaluated by a cardiologist. Patient interviews, sound recordings, and blood samples are obtained in connection with the coronary CTA. If coronary CTA does not rule out obstructive CAD, patients will be examined by 3T CMRI ⁸²Rb-PET, FFR_CT, ICA, and FFR. Reference standard is ICA-FFR. Obstructive CAD is defined as an FFR ≤0.80 or as high-grade stenosis (>90% diameter stenosis) by visual assessment. Diagnostic performance will be evaluated as sensitivity, specificity, predictive values, likelihood ratios, calibration, and discrimination. Enrolment started January 2018 and is expected to be completed by June 2020. Patients are followed for 10 years after inclusion.

DISCUSSION

The results of the Dan-NICAD 2 study are expected to contribute to the improvement of diagnostic strategies for patients suspected of CAD in 3 different steps: risk stratification prior to coronary CTA, diagnostic strategy after coronary CTA, and invasive wireless QFR analysis as an alternative to ICA-FFR.

Dynamic Changes in Microvascular Flow Conductivity and Perfusion After Myocardial Infarction Shown by Image-Based Modeling

Background Microcirculation is a decisive factor in tissue reperfusion inadequacy following myocardial infarction ( MI ). Nonetheless, experimental assessment of blood flow in microcirculation remains a bottleneck. We sought to model blood flow properties in coronary microcirculation at different time points after MI and to compare them with healthy conditions to obtain insights into alterations in cardiac tissue perfusion. Methods and Results We developed an image-based modeling framework that permitted feeding a continuum flow model with anatomical data previously obtained from the pig coronary microvasculature to calculate physiologically meaningful permeability tensors. The tensors encompassed the microvascular conductivity and were also used to estimate the arteriole-venule drop in pressure and myocardial blood flow. Our results indicate that the tensors increased in a bimodal pattern at infarcted areas on days 1 and 7 after MI while a nonphysiological decrease in arteriole-venule drop in pressure was observed; contrary, the tensors and the arteriole-venule drop in pressure on day 3 after MI , and in remote areas, were closer to values for healthy tissue. Myocardial blood flow calculated using the condition-dependent arteriole-venule drop in pressure decreased in infarcted areas. Last, we simulated specific modes of vascular remodeling, such as vasodilation, vasoconstriction, or pruning, and quantified their distinct impact on microvascular conductivity. Conclusions Our study unravels time- and region-dependent alterations of tissue perfusion related to the structural changes occurring in the coronary microvasculature due to MI . It also paves the way for conducting simulations in new therapeutic interventions in MI and for image-based microvascular modeling by applying continuum flow models in other biomedical scenarios.

Simulation of the Perfusion of Contrast Agent Used in Cardiac Magnetic Resonance: A Step Toward Non-invasive Cardiac Perfusion Quantification

This work presents a new mathematical model to describe cardiac perfusion in the myocardium as acquired by cardiac magnetic resonance (CMR) perfusion exams. The combination of first pass (or contrast-enhanced CMR) and late enhancement CMR is a widely used non-invasive exam that can identify abnormal perfused regions of the heart via the use of a contrast agent (CA). The exam provides important information to the diagnosis, management, and prognosis of ischemia and infarct: perfusion on different regions, the status of microvascular structures, the presence of fibrosis, and the relative volume of extracellular space. This information is obtained by inferring the spatiotemporal dynamics of the contrast in the myocardial tissue from the acquired images. The evaluation of these physiological parameters plays an important role in the assessment of myocardial viability. However, the nature of cardiac physiology poses great challenges in the estimation of these parameters. Briefly, these are currently estimated qualitatively via visual inspection of images and comparison of relative brightness between different regions of the heart. Therefore, there is a great urge for techniques that can help to quantify cardiac perfusion. In this work, we propose a new mathematical model based on multidomain flow in porous media. The model is based on a system of partial differential equations. Darcy's law is used to obtain the pressure and velocity distribution. CA dynamics is described by reaction-diffusion-advection equations in the intravascular space and in the interstitial space. The interaction of fibrosis and the CA is also considered. The new model treats the domains as anisotropic media and imposes a closed loop of intravascular flow, which is necessary to reproduce the recirculation of the CA. The model parameters were adjusted to reproduce clinical data. In addition, the model was used to simulate different scenarios: normal perfusion; endocardial ischemia due to stenosis in a coronary artery in the epicardium; and myocardial infarct. Therefore, the computational model was able to correlate anatomical features, stenosis and the presence of fibrosis, with functional ones, cardiac perfusion. Altogether, the results suggest that the model can support the process of non-invasive cardiac perfusion quantification.

Pixel-wise quantification of myocardial perfusion using spatial Tikhonov regularization

Quantification of myocardial perfusion by contrast-enhanced cardiovascular magnetic resonance imaging (CMR) aims for an observer independent and reproducible risk assessment of cardiovascular disease. Currently, the data used for the pixel-wise analysis of cardiac perfusion are either filtered prior to a fitting procedure, which inherently reduces the spatial resolution of data; or all pixels are considered without any regularization or prior filtering, which yields an unstable fit in the presence of low signal-to-noise ratio. Here, we propose a new pixel-wise analysis based on spatial Tikhonov regularization which exploits the spatial smoothness of the data and ensures accurate quantification even for images with low signal-to-noise ratio. The regularization parameter is determined automatically by an L-curve criterion. We study the performance of our method on a numerical phantom and demonstrate that the method reduces significantly the root-mean square error in the perfusion estimate compared to a non-regularized fit. In patient data our method allows us to recover the myocardial perfusion and to distinguish between healthy and ischemic regions.

Feasibility of simultaneous PET-MR perfusion using a novel cardiac perfusion phantom

Background

PET-MR scanners are beginning to be employed for quantitative myocardial perfusion imaging. In order to examine simultaneous perfusion calculations, this work describes a feasibility study of simultaneous PET-MR of gadolinium-based contrast agent (GBCA) and PET radiotracer in a novel cardiac perfusion phantom.

Results

[¹⁸F]F⁻ and GBCA were injected simultaneously into a cardiac phantom using a range of ground-truth myocardial perfusion rates of 1 to 5 ml/g/min. PET quantification of K₁ (ml/g/min) was performed using a single tissue compartment model. MR perfusion was calculated using a model-independent signal deconvolution technique. PET and MR signal traces from the phantom aorta and myocardial sections show true simultaneous PET and MR arterial input functions (AIF) and myocardial uptake respectively at each perfusion rate. Calculation of perfusion parameters showed both K₁ and h(t = 0) (PET and MR perfusion parameters respectively) to be linearly related with the ground truth perfusion rate (P_T), and also linearly related to each other (R² = 0.99). The highest difference in perfusion values between K₁ and P_T was 16% at 1 ml/g/min, and the mean difference for all other perfusion rates was <3%.

Conclusions

The perfusion phantom allows accurate and reproducible simulation of the myocardial kinetics for simultaneous PET-MR imaging, and may find use in protocol design and development of PET-MR based quantification techniques and direct comparison of quantification of the two modalities.

Estimating extraction fraction and blood flow by combining first-pass myocardial perfusion and T1 mapping results

Background

Quantifying myocardial perfusion is complicated by the complexity of pharmacokinetic model being used and the reliability of perfusion parameter estimates. More complex modeling provides more information about the underlying physiology, but too many parameters in complex models introduce a new problem of reliable estimation. To overcome the problem of multiple parameters, we have developed a technique that combines knowledge from two different cardiac magnetic resonance (MR) imaging techniques: dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and T1 mapping. Using extracellular volume (ECV) estimates from T1 mapping may allow more robust model parameter estimates.

Methods

Simulations and human scans were performed. The myocardial perfusion scans used an ungated saturation recovery prepared TurboFLASH pulse sequence. Four short-axis (SA) slices were acquired after a single saturation pulse with a saturation recovery time of ~25 ms before the first slice. Gadoteridol was injected and ~240 frames were acquired over a minute with shallow breathing and no electrocardiograph (ECG) gating. This was followed 20±5 minutes later by an injection of regadenoson to induce hyperemia. The data were acquired using an under-sampled golden angle radial acquisition. Modified look-locker inversion recovery (MOLLI) T1 mapping was performed in 3 slices pre- and post-contrast. The pre- and post-contrast T1 maps were used for ECV estimation. Quantification of perfusion was done using a 4-parameter model with additional information about ECV supplied during model fitting. Phase contrast scans of the coronary sinus (CS) were acquired at rest and immediately after the stress perfusion acquisition to estimate global flow.

Results

Without ECV information, the 5-parameter model fails to converge to a unique solution and often gives incorrect estimates for the perfusion parameters. The myocardial blood flow (MBF) estimates during rest and stress were 0.9±0.1 and 2.3±0.6 mL/min/g, respectively. The extraction fraction estimates were 0.49±0.04 and 0.34±0.05 during rest and stress, respectively.

Conclusions

These results show that it is possible to successfully fit a dynamic perfusion model with an extraction fraction parameter by using information from T1 mapping scans. This hybrid approach is especially important when the 5-parameter model alone fails to converge on a unique solution. This work is a good example of exploiting information overlaps between various cardiac MR imaging techniques.

Validation of myocardial perfusion quantification by dynamic CT in an ex-vivo porcine heart model

To test the accuracy of quantification of myocardial perfusion imaging (MPI) using computed tomography (CT) in ex-vivo porcine models. Five isolated porcine hearts were perfused according to Langendorff. Hearts were perfused using retrograde flow through the aorta and blood flow, blood pressure and heart rate were monitored throughout the experiment. An inflatable cuff was placed around the circumflex (Cx) artery to create stenosis grades which were monitored using a pressure wire, analysing perfusion at several fractional flow reserve values of 1.0, 0.7, 0.5, 0.3, and total occlusion. Second-generation dual-source CT was used to acquire dynamic MPI in shuttle mode with 350 mAs/rot at 100 kVp. CT MPI was performed using VPCT myocardium software, calculating myocardial blood flow (MBF, ml/100 ml/min) for segments perfused by Cx artery and non-Cx myocardial segments. Microspheres were successfully infused at three stenosis grades in three of the five hearts. Heart rate ranged from 75 to 134 beats per minute. Arterial blood flow ranged from 0.5 to 1.4 l min and blood pressure ranged from 54 to 107 mmHg. MBF was determined in 400 myocardial segments of which 115 were classified as ‘Cx-territory’. MBF was significantly different between non-Cx and Cx segments at stenosis grades with an FFR ≤0.70 (Mann–Whitney U test, p < 0.05). MBF showed a moderate correlation with microsphere MBF for the three individual hearts (Pearson correlation 0.62–0.76, p < 0.01). CT MPI can be used to determine regional differences in myocardial perfusion parameters, based on severity of coronary stenosis. Significant differences in MBF could be measured between non-ischemic and ischemic segments.

The history of the microsphere method for measuring blood flows with special reference to myocardial blood flow: a personal memoir

We use many types of equipment and technologies to make our measurements but give little thought to how they developed. Evolution was once described as a series of recoils from blind alleys, and this is exemplified by the gradual development of the microsphere method of measuring blood flows. The microsphere method is one of the most frequently used methods for measuring blood flow to organs and portions of organs. The method can measure myocardial blood flow with reasonable accuracy (within 10%) down to samples weighing >50 mg but probably will not do so for samples weighing 1-10 mg. Microspheres with diameters from 10 to 15 μm provide the best compromise between accurate flow measurement and retention in tissue. Radioactive labels have been almst entirely replaced by fluorescent labels, but colored microspheres and neutron-activated labels are also used.NEW & NOTEWORTHY The contributions of the various individuals who developed the microsphere method of measuring regional blood flows and how these advances took place are brought to light in this paper.

Microvascular ischemia in hypertrophic cardiomyopathy: new insights from high-resolution combined quantification of perfusion and late gadolinium enhancement

BACKGROUND

Microvascular ischemia is one of the hallmarks of hypertrophic cardiomyopathy (HCM) and has been associated with poor outcome. However, myocardial fibrosis, seen on cardiovascular magnetic resonance (CMR) as late gadolinium enhancement (LGE), can be responsible for rest perfusion defects in up to 30% of patients with HCM, potentially leading to an overestimation of the ischemic burden. We investigated the effect of left ventricle (LV) scar on the total LV ischemic burden using novel high-resolution perfusion analysis techniques in conjunction with LGE quantification.

METHODS

30 patients with HCM and unobstructed epicardial coronary arteries underwent CMR with Fermi constrained quantitative perfusion analysis on segmental and high-resolution data. The latter were corrected for the presence of fibrosis on a pixel-by-pixel basis.

RESULTS

High-resolution quantification proved more sensitive for the detection of microvascular ischemia in comparison to segmental analysis. Areas of LGE were associated with significant reduction of myocardial perfusion reserve (MPR) leading to an overestimation of the total ischemic burden on non-corrected perfusion maps. Using a threshold MPR of 1.5, the presence of LGE caused an overestimation of the ischemic burden of 28%. The ischemic burden was more severe in patients with fibrosis, also after correction of the perfusion maps, in keeping with more severe disease in this subgroup.

CONCLUSIONS

LGE is an important confounder in the assessment of the ischemic burden in patients with HCM. High-resolution quantitative analysis with LGE correction enables the independent evaluation of microvascular ischemia and fibrosis and should be used when evaluating patients with HCM.

Perfusion dyssynchrony analysis

AIMS

We sought to describe perfusion dyssynchrony analysis specifically to exploit the high temporal resolution of stress perfusion CMR. This novel approach detects differences in the temporal distribution of the wash-in of contrast agent across the left ventricular wall.

METHODS AND RESULTS

Ninety-eight patients with suspected coronary artery disease (CAD) were retrospectively identified. All patients had undergone perfusion CMR at 3T and invasive angiography with fractional flow reserve (FFR) of lesions visually judged >50% stenosis. Stress images were analysed using four different perfusion dyssynchrony indices: the variance and coefficient of variation of the time to maximum signal upslope (V-TTMU and C-TTMU) and the variance and coefficient of variation of the time to peak myocardial signal enhancement (V-TTP and C-TTP). Patients were classified according to the number of vessels with haemodynamically significant CAD indicated by FFR <0.8. All indices of perfusion dyssynchrony were capable of identifying the presence of significant CAD. C-TTP >10% identified CAD with sensitivity 0.889, specificity 0.857 (P < 0.0001). All indices correlated with the number of diseased vessels. C-TTP >12% identified multi-vessel disease with sensitivity 0.806, specificity 0.657 (P < 0.0001). C-TTP was also the dyssynchrony index with the best inter- and intra-observer reproducibility. Perfusion dyssynchrony indices showed weak correlation with other invasive and non-invasive measurements of the severity of ischaemia, including FFR, visual ischaemic burden, and MPR.

CONCLUSION

These findings suggest that perfusion dyssynchrony analysis is a robust novel approach to the analysis of first-pass perfusion and has the potential to add complementary information to aid assessment of CAD.