Contrast-Enhanced Magnetic Resonance Imaging Based T1-Mapping and Extracellular Volume Fractions are Associated with Peripheral Artery Disease

Introduction

A natural outgrowth of the past decade’s advancements in magnetic resonance imaging (MRI) and computed tomography (CT) has been the enhanced diagnostic accuracy of myriad ailments owing to the more robust processing power of state-of-the-art imaging modalities¹. Notwithstanding the reproducibly more accurate detection of pathologies previously encumbered by limitations in imaging quality and efficiency, paradigm-changing advances in the disease monitoring and diagnosis realms—ones that detect diseases during their clinically silent stages rather than the conventional presentation when curative interventions are at best a long shot— continue to elude clinical fruition.

With advances in resolution and processing power occurring amid a deep learning and artificial intelligence (AI) revolution², the opportunity has emerged to, at still-more granular radiographic levels, define pathophysiologic underpinnings that not only more accurately delineate whether an anomaly is present but stratify the abnormality based on the most probable culprit central process or disease^3,4. Such advancements in precision medicine and radiography have, for example, already been realized in the breast cancer diagnosis realm with significantly improved early detection of what would otherwise, in the setting of conventional surveillance modalities, prove clandestine tumors—and with superior efficacy capable of surmounting historically challenging clinical diagnostic conundrums^5,6.

To this end, radiomics-driven solutions are continually demonstrating game-changing potential in the effort to deliver early diagnostic and disease monitoring solutions to the cardiovascular disease arena down to previously unseen pathophysiologic levels. Notable advancements to this end include Bayesian networks that leverage data from stress echocardiography and CT angiography (CTA) to accurately predict future cardiovascular events such as myocardial infarction (MI) or death^7,8 a machine learning algorithm that outperforms conventional diagnostic strategies for accurate detection of heart failure with preserved ejection fraction⁹ among myriad other novel algorithms that show promise in the diagnosis of obstructive coronary stenoses¹⁰, myocardial infarction¹¹, and peripheral vascular disease (PVD)¹². Owing to the grim 10-year mortality rate facing patients with limb claudication that eclipses 50% when considering age and comorbid factors¹³ and the increased rate of myocardial infarction, stroke, and cardiovascular death in patients with both symptomatic and asymptomatic peripheral artery disease (PAD)¹⁴, it is imperative that one such paradigm change occurs in the surveillance for the earliest warning signs of limb ischemia.

One avenue of promise has emerged and stems from advances in quantifying the extracellular volume fraction (ECV) in the myocardium measured with contrast-enhanced magnetic resonance imaging (CE-MRI). The T1 relaxation time, a measure of how fast the nuclear spin magnetization equilibrates after a radiofrequency pulse in MRI, is an important tool for highlighting soft-tissue contrast. T1 times are specific to and span myriad pathologies, from fibrosis to fatty infiltration. Owing to the pathological underpinnings of fibrosis, arising from derailments in and expansion of extracellular matrix (ECM) remodeling¹⁵, contrast-enhanced quantification of ECV has proven a useful means for detecting microscopic areas of ischemia, myocardial scarring, and perfusion^13,16.

Advances in cardiac MRI leverage pre- and post-contrast injection in T1 sequences to disclose precise, focal areas of pathology by measuring the rate of change of T1 relaxation relative to the degree of contrast pooling, with T1 times corresponding to pathologic states. For example, acute membrane rupture emblematic of acute MI and myocarditis allows aberrant infiltration of gadolinium contrast into myocytes, resulting in increased gadolinium concentration, prolonged T1 relaxation times, and the corresponding hyper-enhancement hallmark of acute MI. By the same token, in chronic MI, scarring prevails, with a corresponding expansion of the interstitium that leads to increased gadolinium concentration, hyper-enhancement, and correspondingly prolonged T1 times.

Previous works^17,18 have demonstrated the potential reproducibility of T1 mapping beyond central cardiovascular to peripheral vascular pathophysiologic areas¹⁹. In this study, we investigated T1 times in PAD patients and matched controls. We hypothesized that differences in T1 relaxation times will be observed in PAD patients compared with appropriately matched controls.

Methods

The study was approved by the institutional review board and conducted in keeping with the Helsinki Declaration of 1973. Study participants provided informed consent, and the study obtained approval from the institutional review boards at the Houston Methodist Hospital and the Michael E. DeBakey Veterans Affairs Medical Center, Houston, TX. Suitable patients from the Houston Methodist Hospital and the Michael E. DeBakey Veterans Affairs Medical Center in Houston, TX were identified using inclusion and exclusion criteria in the medical record. PAD patients and controls without PAD were recruited and matched on the basis of demographic and clinical characteristics. Standard of care continued uninterrupted throughout study subject participation in this observational imaging work. Patients with contraindications to MRI and those with an estimated glomerular filtration rate (eGFR) ≤40 mL/min/1.73.m² were excluded from this study.

Participants were positioned on the MRI table feet first in the supine position, and imaging was performed at the mid-calf level. Efforts to minimize motion artifact were employed using a bilateral leg coil employed during the imaging runtime. A gadolinium-based contrast agent was administered intravenously (gadopentetate dimeglumine [Magnevist, Bayer Healthcare, Whippany, New Jersey, USA at 0.2 mmol/kg, or gadobutrol [Gadavist, Bayer Healthcare, at 0.1 mmol/kg) with flow rates of 2-4 ml/sec proceeded by a 20 ml saline flush. MRI scans were saved in the DICOM format prior to rendering and further analysis. The arterial lumens of the anterior tibialis (AT), posterior tibial (PT), and peroneal artery (PE) were segmented and visualized as permitted by contrast enhancement using Sante DICOM Editor Version 3.0 (Santesoft, Nicosia, Cyprus).

The reproducibility and quality of the lumen frame analysis and interpretation were assessed through inter- and intraobserver reproducibility analyses. Cross-sectional leg muscle area (CSLMA) was measured using our previously described methodology.¹⁸

T1-mapping was done before and after gadolinium contrast administration with a motion-corrected Modified Look-Locker Inversion Recovery (MOLLI) pulse sequence. T1 values are reported in millisecond (ms) and were calculated with a 3-parameter Levenberg-Marquardt curve fitting algorithm. We reported the peak T1 values across the five calf muscle compartments including the anterior (AM), lateral (LM), and deep posterior (DM) muscle groups; soleus (SM), and gastrocnemius (GM) muscles. Averaged peak blood pool T1 values were obtained from the posterior (PT) and anterior tibialis (AT) and peroneal artery (PE). PAD ECV was calculated for each muscle group as (1-hematocrit in %) x [ (1/T1m post) - (1/T1m pre)] / [ (1/T1b post) - (1/T1b pre)], where T1m denotes the maximum T1 value of the respective muscle compartments and T1b the averaged peak arterial blood T1 value.

Variable normality was determined with the Shapiro-Wilk test. Group differences for categorical variables were analyzed with the Chi-square or Fisher exact test. Continuous variables were analyzed with an independent samples Student’s t-test, and the Mann-Whitney-Wilcoxon test was used for non-normal variables. Associations between T1 times and markers of PAD were determined using linear regression. Sante DICOM Editor 3D was used for muscle compartment and lower extremity tracings, T1 extrapolation, and T1 mapping. Intra- and interclass correlations, performed on Stata with the ICC command using DICOM outputs from Sante, were conducted to assess variability in tracing and MRI data analysis between and within operators.

All statistical tests were two-sided, and a p-value <0.05 was used as the cutoff for statistically significant phenomena. The statistical analyses were performed with Stata Statistical Software: Release 13 (College Station, Texas, StataCorp LP).

Results

Discussion

The major finding of this work was that native maximum T1 times averaged over five skeletal muscle groups were higher in PAD patients compared to matched controls. This creates an opportunity for future prospective studies of T1 times and PAD presence and severity.

Prior works have demonstrated the value of T1-mapping in coronary artery disease and left ventricular dysfunction. The findings of this study add to the growing evidence that T1-mapping is of interest in assessing cardiovascular disease. Our study has limitations. The heterogeneity of findings, namely with regard to T1 relaxation times. Additionally, the small sample size of this work precludes generalizability of its findings but does provide early groundwork that future works with larger patient cohorts can expound upon to further assess the role of native peak T1 and ECV in the study of PAD.

In conclusion, native peak T1 values and ECV fractions of several calf muscle compartments are higher in PAD patients compared with matched controls. These data suggest a potential utility of ECV to noninvasively determine skeletal muscle fibrosis in PAD.

Conflicts of Interest

The authors declare no conflicts of interest.