Skip to main content

The reliability and feasibility of non-contrast adenosine stress cardiovascular magnetic resonance T1 mapping in patients on haemodialysis



Identifying coronary artery disease (CAD) in patients with end-stage renal disease (ESRD) is challenging. Adenosine stress native T1 mapping with cardiovascular magnetic resonance (CMR) may accurately detect obstructive CAD and microvascular dysfunction in the general population. This study assessed the feasibility and reliability of adenosine stress native T1 mapping in patients on haemodialysis.


The feasibility of undertaking rest and adenosine stress native T1 mapping using the single-shot Modified Look-Locker inversion recovery (MOLLI) sequence was assessed in 58 patients on maintenance haemodialysis using 3 T CMR. Ten patients underwent repeat stress CMR within 2 weeks for assessment of test-retest reliability of native T1, stress T1 and delta T1 (ΔT1). Interrater and intrarater agreement were assessed in 10 patients. Exploratory analyses were undertaken to assess associations between clinical variables and native T1 values in 51 patients on haemodialysis.


Mean age of participants was 55 ± 15 years, 46 (79%) were male, and median dialysis vintage was 21 (8; 48) months. All patients completed the scan without complications. Mean native T1 rest, stress and ΔT1 were 1261 ± 57 ms, 1297 ± 50 ms and 2.9 ± 2.5%, respectively. Interrater and intrarater agreement of rest T1, stress T1 and ΔT1 were excellent, with intraclass correlation coefficients (ICC) > 0.9 for all. Test-retest reliability of rest and stress native T1 were excellent or good (CoV 1.2 and 1.5%; ICC, 0.79 and 0.69, respectively). Test-retest reliability of ΔT1 was moderate to poor (CoV 27.4%, ICC 0.55). On multivariate analysis, CAD, diabetes mellitus and resting native T1 time were independent determinants of ΔT1 (β = − 0.275, p = 0.019; β = − 0.297, p = 0.013; β = − 0.455; p < 0.001, respectively).


Rest and adenosine stress native T1 mapping is feasible and well-tolerated amongst patients with ESRD on haemodialysis. Although rater agreement of the technique is excellent, test-retest reliability of ΔT1 is moderate to poor. Prospective studies should evaluate the relationship between this technique and established methods of CAD assessment and association with outcomes.


Patients with chronic kidney disease and end stage renal disease (ESRD) experience excess mortality from cardiovascular disease [1].. Developing safe and robust investigations that can reliably identify coronary artery disease (CAD) to enable optimal risk stratification, intervention and management is a priority for this population [2]. Contrast-enhanced cardiovascular magnetic resonance (CMR) is a validated and established technique for the assessment of CAD, allowing evaluation of both myocardial perfusion defects and infarct burden [3, 4]. However, this requires administration of gadolinium-based contrast, which is contraindicated in patients with advanced kidney disease, due to the rare but serious risk of nephrogenic systemic fibrosis [5].

Native T1 mapping is a non-contrast parametric mapping technique on CMR that allows quantitative tissue characterisation on a pixel-by-pixel basis [6]. T1 relaxation time increases with higher free water content [7]. The total myocardial T1 value is derived from the contribution of the intracellular, extracellular and blood compartments [8]. Under normal circumstances, adenosine vasodilator stress causes an increase in myocardial blood volume and, consequently, an increase in myocardial native T1 [9]. Myocardial stress T1 reactivity (ΔT1) is, therefore, thought to represent the percentage increase in myocardial blood volume during vasodilator stress [10]. As disturbances in absolute quantification of myocardial blood volume changes can detect significant coronary artery stenoses [11], a blunted ΔT1 could be a non-invasive non-contrast measure to differentiate ischaemic, infarcted and normal myocardium [10]. Adenosine stress T1 mapping has recently been demonstrated to accurately detect obstructive CAD and microvascular dysfunction in patients with stable angina and healthy subjects [12]. However, this technique has never been tested in a population with ESRD. These patients have a complex cardiovascular disease phenotype, marked by both coronary and non-coronary disease-related processes, high levels of interstitial myocardial fibrosis [13, 14] and significantly higher myocardial native T1 times compared to control subjects [15, 16]. Whether stress T1 mapping could be an appropriate test to assess CAD in this patient group is therefore not known.

In this study we assessed the feasibility of non-contrast rest and stress T1 mapping in patients with ESRD on haemodialysis. We aimed to assess the ΔT1, its reliability with test-retest measurements, interrater and intrarater agreement in patients on haemodialysis and we conducted an exploratory analysis to investigate associations between clinical variables and ΔT1.


Study population

Patients with ESRD on haemodialysis were recruited as part of the CYCLE-HD study (ISRCTN11299707). The study was approved by the NHS Research Ethics Committee East Midlands (Northampton; REC ref.: 14/EM/1190). All patients provided written informed consent to participate in the study. Adult patients undergoing in-centre maintenance haemodialysis for at least 3 months were eligible for inclusion in the study, provided they were able to undertake exercise and undergo CMR scanning. Inclusion and exclusion criteria are as previously described [17].

CMR image protocol

CMR imaging was performed on a 3 T CMR platform (Skyra, Siemens Healthineers, Erlangen, Germany) with an 18-channel phased-array anterior coil. To ensure consistency of fluid status, all patients underwent CMR assessment on a non-dialysis day, not after the two-day break, between 18 and 24 h of the last haemodialysis session. The same standardised CMR protocol was used for all study participants as previously described [16], in line with internationally recognised standards [18]. Briefly, mid-ventricular, short-axis native T1 maps were acquired using the single-shot Modified Look-Locker inversion recovery (MOLLI) sequence [19] with the 3(3)3(3)5 sampling pattern. Images were acquired at end-expiration, using free-breathing with motion correction and electrocardiogram (ECG)-gating. Typical parameters were: slice thickness 8.0 mm, field of view 300 × 400 mm, flip angle 500, minimum T1 120ms, inversion-time increment 80 ms. The MOLLI sequence was chosen due to local expertise [20] and the excellent reliability demonstrated in this patient population at 3 T [21]. All participants were asked to avoid caffeine on the day of their scan. At the study visit, they had the choice to undergo adenosine stress scanning or not. If they agreed to it, after the acquisition of rest images, adenosine was infused intravenously at an initial dose of 140 μg/kg/min for 3 minutes. Adequate haemodynamic response was defined as the occurrence of two of the following: ≥10% increase in heart rate; ≥10 mmHg decrease in systolic blood pressure; reporting of typical symptoms (e.g. flushing, shortness of breath, chest discomfort). In case of inadequate haemodynamic response, adenosine dose was increased incrementally to 170 μg/kg/min and to a maximum of 210 μg/kg/min [22]. Stress T1 maps were acquired in a mid-ventricular short-axis slice in the same slice position as resting T1 maps.

CMR image analysis

All scans were anonymised and analysed offline by blinded observers. The software package CMR42 (version 5.10.1, Circle Cardiovascular Imaging, Calgary, Alberta, Canada) was used for both left ventricular (LV) and native T1 analyses. LV quantification was performed by a single reader (MGB) as previously described [16]. A separate reader completed T1 map analysis (FP) as previously described [21]. Briefly, endocardial and epicardial contours were manually drawn for native T1 maps. Care was taken to allow sufficient margins of separation between the myocardium and other tissue interfaces to avoid partial blood pool volume or epicardial fat contaminations. After the definition of the anterior right ventricular insertion point, the mid-ventricular short axis slice was automatically divided into 6 segments, according to the American Heart Association 16-segment model (Fig. 1). Infarcted segments, which were too thin to contour (with the inability to differentiate myocardium from blood pool) or with clear akinesis on cine imaging were excluded, as were myocardial segments with artefact. If certain segments required exclusion, the average myocardial T1 time was calculated as the mean of the remaining segments. If no segments required exclusion, the analysis software generated the average myocardial T1 value using a three-parameter least-squares fitting technique with heart rate correction. ΔT1 was calculated, as previously described using the following formula [8, 10, 12, 23,24,25]:

$$ \Delta T1=\frac{stress\ T1- rest\ T1}{rest\ T1}\times 100 $$
Fig. 1
figure 1

T1 map analysis and segmentation of mid-ventricular slice (a) at rest and (b) during adenosine stress, based on definition of right ventricular insertion point (blue dot)

Image quality was graded as excellent, moderate or non-analysable.

An additional, descriptive analysis comparing regional wall motion (RWM) scores to ΔT1 scores was completed. The mid-ventricular, short-axis LV cine image corresponding to the mid-ventricular native T1 maps for each subject were scored for RWM abnormalities. Each segment was scored 1 (normal), 2 (hypokinetic), 3 (akinetic) or 4 (dyskinetic) (AA). Segments were then compared to corresponding ΔT1 segment for analysis.

Test-retest reliability, interrater and intrarater agreement

Ten patients underwent an identical CMR scan within 2 weeks of the baseline scan as part of a reliability sub-study. Test-retest reliability was assessed by a single, blinded observer (FP). Interrater agreement was assessed by analysis of the same 10 scans by two independent, blinded readers (FP; MGB). Intrarater agreement was assessed by re-analysis of 10 scans by a single blinded reader (FP), 1 month apart.

This sample size of 10 was supported by a power calculation previously done to detect a 2.5% difference in native T1 between test-retest scans with 90% power [21].

Exploratory analyses

Exploratory analyses were undertaken to establish potential associations between ΔT1 and clinical indicators of CAD, including age, sex, heart rate (HR), blood pressure (BP), smoking status, previous history of CAD, diabetes mellitus, previous renal transplant, hypertension and dyslipidaemia. All CMR images used for the exploratory analysis were acquired as part of the baseline assessment, with patient demographic details, comorbidities, medications, haematological and biochemical data were recorded prospectively. Total time on renal replacement therapy (RRT), time with a functioning kidney transplant and haemodialysis vintage were prospectively collected and calculated from electronic hospital records. The presence of CAD was defined as any of: reported diagnosis of CAD on baseline assessment; fixed or reversible ischaemia on myocardial perfusion scan; significant stenosis on invasive coronary angiography; infarcted myocardium on CMR; previous coronary revascularisation (coronary artery bypass graft or percutaneous coronary intervention). Individuals who had a non-analysable stress scan were excluded from this analysis.

Statistical analysis

Statistical analyses were performed using SPSS (v 25.0, Statistical Package for Social Sciences, International Business Machines, Inc.,, Armonk, New York, USA) and GraphPad Prism (v 8.1.1, GraphPad Software, Inc., La Jolla, California, USA). Categorical data are reported as frequency (%) of observation. Normality was assessed using histograms and Shapiro-Wilk tests for continuous data. Normally distributed data are expressed as mean ± standard deviation (SD). Non-normally distributed data are expressed as median (P25, 75). Test-retest reliability, interrater and intrarater agreement were assessed using coefficients of variability (CoV), intraclass correlation coefficients (ICC) and Bland-Altman analyses [26]. Reliability was considered excellent for CoV < 10%, ICC > 0.90; good for CoV < 15%, ICC 0.75–0.90; moderate for CoV < 20%, ICC 0.50–0.74 and poor for CoV > 20%, ICC < 0.50. Comparisons between groups were made using independent sample t-tests data and Mann-Whitney U tests for normally and non-normally distributed data, respectively. Differences between nominal variables were assessed using Chi-squared and Fisher’s exact tests. Correlations between variables were assessed with Pearson’s and Spearman’s rank analysis for normally and non-normally distributed data, respectively. Multivariable linear regression analyses were conducted to assess independent determinants of ΔT1. Variables inputted into regression models were limited to avoid overfitting. Statistical significance level was defined as a two-tailed p-value < 0.05. The exploratory analysis comparing RWM to regional ΔT1 are presented as descriptive statistics due to the imbalance in numbers between groups.


Feasibility and analysability

The baseline scans of 58 of the 130 participants of the CYCLE-HD study agreed to have both rest and adenosine stress native T1 mapping. Baseline characteristics are shown in Table 1. Additional haematological, biochemical and medication data can be found in Appendix 1. A consort diagram showing patients included in the study is shown in Fig. 2. Only two of the 124 (1.6%) resting T1 maps were not analysable due to poor image quality. Image quality was either excellent (n = 93) or moderate (n = 29) for the remaining 122 resting T1 maps. Seven of 58 (12%) stress T1 maps were not analysable. Image quality for the remaining 51 stress T1 maps was either excellent (n = 32) or moderate (n = 18). All patients with analysable stress T1 maps also had analysable matched resting T1 maps. Of the 58 participants who underwent stress testing, six required an increase in adenosine dose. Only one participant had an inadequate haemodynamic response despite maximal dose of adenosine. All patients tolerated and successfully completed the scan with no complications necessitating termination of the scan. After the administration of adenosine, there was a significant increase in HR (71 ± 12 bpm versus 85 ± 14 bpm; p < 0.001) and a significant drop in systolic BP (152 ± 34 mmHg versus 136 ± 35 mmHg; p < 0.001). 86% of the patients experienced mild symptoms (Table 2). There was no significant association between heart rate and native T1 values at rest (r = − 0.096; p = 0.476) or during stress (r = − 0.118; p = 0.387).

Table 1 Baseline characteristics
Fig. 2
figure 2

Consort diagram

Table 2 Summary of stress response


The mean interval between test and retest reliability scans was 7 ± 4 days. All ten patients had test-retest resting T1 maps as part of the reliability sub-study, but only nine underwent stress testing on the retest scan (one had caffeine). For test-retest scans, all T1 images were analysable. Image quality was either excellent (rest, n = 18; stress, n = 11) or moderate (rest, n = 2; stress, n = 4). The interrater and intrarater agreement of all parameters were excellent (ICC > 0.9). Test-retest reliability was good-excellent for rest and stress native T1, while it was moderate-poor for ΔT1 (Table 3 for global measurements and Table 4 for segmental measurements). Segmental values for ΔT1 test-retest reliability were not analysed given the already moderate-poor reliability of global measurements. Bland-Altman plots did not show evidence of systematic bias, although some data points were outside the 95% confidence interval (Fig. 3).

Table 3 Test-retest reliability, interrater and intrarater agreement of global rest native T1, stress T1 and ΔT1
Table 4 Test-retest reliability, interrater and intrarater agreement for individual mid-ventricular myocardial segments
Fig. 3
figure 3

Bland-Altman plots for test-retest reliability of (a) resting native T1, (b) stress native T1, (c) ΔT1; interrater agreement of (d) resting native T1, (e) stress native T1, (f) ΔT1; intrarater agreement of (g) resting native T1, (h) stress native T1, (i) ΔT1

Associations between clinical variables and myocardial stress T1 reactivity

For the 51 participants with analysable stress T1 maps, an exploratory analysis was conducted to establish associations between ΔT1 and clinical characteristics. ΔT1 was not significantly different between males and females (p = 0.51). ΔT1 decreased with increasing age and increasing haemoglobin A1c (HbA1c) (Fig. 4). ΔT1 was blunted in patients with evidence of CAD (1.2% ± 1.8 versus 3.4% ± 2.5, p = 0.004) or diabetes (1.9% ± 2.1 versus 3.6% ± 2.6; p = 0.019) (Fig. 5). Patients who had a previous kidney transplant showed a higher ΔT1 compared to those who solely had dialysis as RRT (4.3% ± 2.9 versus 2.2% ± 2.0; p = 0.004) (Fig. 5). However, the proportion of total RRT time spent with a functioning kidney transplant as opposed to dialysis was not significantly associated with ΔT1. There was no clear relationship between ΔT1 in all-time smokers compared to non-smokers (2.4% ± 2.4 versus 3.3% ± 2.8; p = 0.23). ΔT1 was not significantly different in patients with or without hypertension or dyslipidaemia. ΔT1 was not significantly associated with HR and BP. ΔT1 was inversely related to resting native T1 (r = − 0.533; p < 0.01) but it was not significantly correlated with stress native T1.

Fig. 4
figure 4

Scatter plots displaying correlations between (left) ΔT1 and age; (right) ΔT1 and glycosylated haemoglobin (HbA1c)

Fig. 5
figure 5

Differences in myocardial stress T1 reactivity (ΔT1) between groups. ΔT1 was significantly lower in patients with coronary artery disease; diabetes mellitus compared to patients without these pathologies. Patient who previously had a kidney transplant displayed significantly higher ΔT1 compared to patients who had only been on dialysis

On multivariate analysis, CAD, diabetes and resting native T1 remained independent determinants of ΔT1, after adjustment for age and proportion of total RRT time with a functioning graft (Table 5).

Table 5 Multivariate linear regression model for the independent predictors of myocardial T1 reactivity (ΔT1)

Regional wall motion scores and ΔT1

For the 51 participants with analysable stress T1 maps 306 segments were analysed on corresponding mid-ventricular LV cine images. 279 segments were scored ‘1’, 19 segments were scored ‘2’ and 8 segments were scored ‘3’. All but one of the segments scored ‘3’ corresponded to an area on T1 maps identified as infarct which was too thin to contour. Mean ΔT1 in segments scored ‘1’ was 3.01% (3.5), mean ΔT1 in segments scored ‘2’ was 1.6% (3.4) and the only segment scored ‘3’ with an analysable T1 maps had a ΔT1 of 1.5%.


This is the first study to assess the feasibility of stress T1 mapping in patients on haemodialysis. We have shown that adenosine stress T1 mapping is a safe, feasible and well-tolerated technique for the potential assessment of myocardial ischaemia in a cohort of 58 patients with ESRD on haemodialysis and that there were no adverse events associated with its use. We have shown native T1 mapping at rest and during stress to be a highly reproducible technique and have described the reliability and agreement of ΔT1 for the first time in this population, which was moderate. ΔT1 was independently associated with CAD, diabetes and resting native T1.

Our group has previously demonstrated the excellent reliability of rest native T1 mapping in patients with ESRD, independent of fluid status [21]. In our study, ΔT1 appears to be a less reliable measure. This may be attributable to the methodology for ΔT1 calculation, as it is based on two independent measurements, amplifying the possibility of variability. Liu et al. assessed reliability with Bland-Altman analysis [12], therefore our study is the first to describe ΔT1 reliability in terms of CoV and ICC, precluding direct comparison of these values with previously reported data. The test-retest bias on Bland-Altman analysis for ΔT1 we report (0.6% ± 1.8) is poorer compared to that described by Liu et al. (0.18% ± 0.36) [12]. This discrepancy might be attributable to differences in the selection of patients, as our study focusses on prevalent ESRD patients, some with comorbidities such as CAD, hypertension and diabetes, whereas Liu et al. assessed ΔT1 reliability in healthy subjects [12]. The difference in acquisition technique (MOLLI versus ShMOLLI) might also contribute to the increased variability, however a head-to-head comparison of the two techniques is needed to confirm this. Cardiac and respiratory motion remain a concern for MOLLI-based acquisition during stress as they lead to loss of spatial resolution, in the presence of high HR, RR variability and residual respiratory movement, which are not always amenable to inline motion correction [24, 27]. The lack of a significant correlation between HR and native T1 values at rest or during stress in our study suggests HR did not impact T1 value itself, despite the HR sensitivity of the classic MOLLI 3(3)3(3)5 technique [19]. This is probably due to the fact that only rarely the HR during stress was above 100 bpm, where the MOLLI sequence is known to be less reliable. When compared to the test-retest reliability reported for the assessment of myocardial perfusion reserve on CMR by semi-quantitative (myocardial perfusion reserve index) or quantitative (myocardial perfusion reserve) analysis, ΔT1 as a method for assessment of myocardial blood volume change performs comparably. Test-retest reliability of myocardial perfusion reserve index has been reported to have a CoV of 19 to 27% in healthy subjects and patients with CAD [28, 29]. For myocardial perfusion reserve, studies have reported test-retest CoVs varying from of 13.3 to 35% [29, 30] in healthy subjects and ICCs varying from 0.26 to 0.88 [30, 31].

We have shown the mean ΔT1 for this population is 2.9 ± 2.5%, substantially below the ΔT1 previously reported in healthy population (6.3 ± 1.1%) [10] and equivalent to that for myocardium with microvascular dysfunction (3.0 ± 0.9%) [12]. This is as expected as patients with ESRD on haemodialysis are known to have a high preponderance of micro and macrovascular dysfunction. Moreover, ΔT1 appears to be able to differentiate between patients with and without CAD and diabetes, further demonstrating the potential clinical application of the technique. Our finding that a history of CAD and diabetes were both independent determinants of ΔT1 in multivariate analysis is in keeping with the blunted ΔT1 reported in the presence of macrovascular epicardial stenosis [10, 12] or microvascular disease [12, 23]. However, although our study does suggest that the presence of CAD or diabetes are clinical determinants a blunted ΔT1, the wide standard deviations we report in this study question the applicability of the ΔT1 thresholds proposed by Liu et al. to detect obstructive CAD and microvascular dysfunction (1.5 and 4.0%, respectively) [12] to patients with ESRD. The description of putative differences between segmental RWM scores and ΔT1 should be viewed as exploratory and hypothesis generating but do not contradict the hypothesis that ΔT1 is blunted in areas of reduced myocardial blood supply.

We acknowledge that we have analysed ΔT1 of the whole mid-ventricular short axis slice, as opposed to individual coronary territories meaning that the ΔT1 might average across normal and non-normal segments, possibly leading to a higher degree of variation. We also found that increased resting myocardial native T1 was an independent determinant of a blunted ΔT1 amongst patients with ESRD. Two possible explanations may account for this finding. Firstly, a higher resting native T1 might signify increased levels of myocardial fibrosis or other pathological processes related to excess water, such as oedema or inflammation [6]. Perivascular myocardial fibrosis is known to predispose to diminished coronary reserve [32, 33] and it could, therefore, limit the increase in myocardial blood volume and result in a blunted ΔT1. Secondly, an elevated resting native T1 could indicate compensatory microcirculatory vasodilation downstream of obstructive epicardial coronary stenosis to preserve myocardial oxygen supply [34]. This would cause an expansion of the intravascular space and subsequent myocardial blood volume increase [35, 36], reflected in the elevated resting T1, in keeping with what was observed by Liu et al. [10, 12]. Our study did find a trend towards higher resting native T1 amongst patients with CAD compared to patients without (1279 ms ± 40 versus 1255 ms ± 60), however, this was not statistically significant, possibly due to inadequate sample size and a type II error. The lack of a trend towards higher resting native T1 in individuals with DM possibly suggests that microvascular dysfunction does not cause a similar expansion of baseline intravascular volume and a subsequent elevation of resting native T1, in keeping with the results of previous studies [10, 12, 37].

With regards to concerns of the influence of body water content on native T1 measurement, our group have previously demonstrated that myocardial native T1 was unaffected by changes in fluid status [21]. Furthermore, a subsequent study assessed changes in myocardial native T1 pre and post haemodialysis. Whilst they found a statistically significant difference, the absolute mean change in native T1 was only 13 ms, representing approximately a 1% change [38].

Study limitations

The cross-sectional design of this feasibility study precludes from any causal relationship to be inferred from the findings we present and does not allow the assessment of clinical endpoints including cardiac events and mortality. Our study was conducted at 3 T using the classic MOLLI sequence, therefore our findings in terms of reliability are not necessarily generalisable to a 1.5 T platform, other MOLLI variants or saturation-recovery single-shot acquisition sequences. Despite not finding a significant correlation between HR and T1 values, we recognise that the HR dependency of the MOLLI sequence may have led to underestimation of stress T1 values and, consequently ΔT1, although maximum HR was rarely elevated above 100 bpm. While this study offers important information with regards to the feasibility, reliability and clinical determinants of ΔT1 in patients with ESRD, the lack of basal and apical slices meant that we could not build the complete American Heart Association 16-segment model for analysis of individual coronary territories and global myocardium. We were unable to provide a direct comparison of stress T1 mapping findings against invasive coronary angiography or non-invasive investigations for CAD, as there was no simultaneous screening for CAD, meaning that we cannot be sure we captured all patients with CAD. This also meant that CAD was assessed as a binary variable, without taking into account extent or severity of disease. Additionally, although obviously infarcted segments of myocardium were excluded from analysis, we acknowledge small areas of replacement fibrosis might not have been detected and, therefore, excluded. As high levels of replacement fibrosis would account for a blunted ΔT1, this is another reason why direct comparison with established investigations for CAD is required.


Rest-vasodilator stress native T1 mapping is feasible in ESRD although test-retest reliability is moderate. These data support the development of an adequately powered prospective study to evaluate the diagnostic accuracy of stress T1 mapping in an ESRD population through direct comparison to invasive coronary angiography and non-invasive measures of CAD.

Availability of data and materials

The datasets generated and/or analysed during the current study are not publicly available but are available from the corresponding author on reasonable request.



Blood pressure


Coronary artery disease


Cardiovascular magnetic resonance


Coefficient of variability




End stage renal disease


Heart rate


Intraclass correlation coefficient


Left ventricle/left ventricular


Modified Look-Locker inversion recovery


Renal replacement therapy


Regional wall motion


Myocardial stress T1 reactivity


  1. Saran R, Robinson B, Abbott KC, Agodoa LY, Bragg-Gresham J, Balkrishnan R, et al. US Renal Data System 2017 Annual Data Report: epidemiology of kidney disease in the United States 2018.

    Google Scholar 

  2. Poli FE, Gulsin GS, McCann GP, Burton JO, Graham-Brown MP. The assessment of coronary artery disease in patients with end-stage renal disease. Clin Kidney J. 2019;12(5):721–34.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Nagel E, Greenwood JP, McCann GP, Bettencourt N, Shah AM, Hussain ST, et al. Magnetic resonance perfusion or fractional flow Reserve in Coronary Disease. N Engl J Med. 2019;380(25):2418–28.

    Article  PubMed  Google Scholar 

  4. von Knobelsdorff-Brenkenhoff F, Schulz-Menger J. Cardiovascular magnetic resonance imaging in ischemic heart disease. J Magn Reson Imaging. 2012;36(1):20–38.

    Article  Google Scholar 

  5. Kay J. Nephrogenic systemic fibrosis: a gadolinium-associated fibrosing disorder in patients with renal dysfunction. Ann Rheum Dis. 2008;67(Suppl 3):iii66–9.

    Article  PubMed  Google Scholar 

  6. Moon JC, Messroghli DR, Kellman P, Piechnik SK, Robson MD, Ugander M, et al. Myocardial T1 mapping and extracellular volume quantification: a Society for Cardiovascular Magnetic Resonance (SCMR) and CMR working Group of the European Society of cardiology consensus statement. J Cardiovasc Magn Reson. 2013;15(1):92.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Piechnik SK, Neubauer S, Ferreira VM. State-of-the-art review: stress T1 mapping—technical considerations, pitfalls and emerging clinical applications. MAGMA. 2018;31(1):131–41.

    Article  PubMed  Google Scholar 

  8. Mahmod M, Piechnik SK, Levelt E, Ferreira VM, Francis JM, Lewis A, et al. Adenosine stress native T1 mapping in severe aortic stenosis: evidence for a role of the intravascular compartment on myocardial T1 values. J Cardiovasc Magn Reson. 2014;16(1):92.

    Article  PubMed  PubMed Central  Google Scholar 

  9. McCommis KS, Zhang H, Goldstein TA, Misselwitz B, Abendschein DR, Gropler RJ, et al. Myocardial blood volume is associated with myocardial oxygen consumption: an experimental study with cardiac magnetic resonance in a canine model. JACC Cardiovasc Imaging. 2009;2(11):1313–20.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Liu A, Wijesurendra RS, Francis JM, Robson MD, Neubauer S, Piechnik SK, et al. Adenosine stress and rest T1 mapping can differentiate between ischemic, infarcted, remote, and normal myocardium without the need for gadolinium contrast agents. JACC Cardiovasc Imaging. 2016;9(1):27–36.

    Article  PubMed  PubMed Central  Google Scholar 

  11. McCommis KS, Goldstein TA, Abendschein DR, Misselwitz B, Pilgram T, Gropler RJ, et al. Roles of myocardial blood volume and flow in coronary artery disease: an experimental MRI study at rest and during hyperemia. Eur Radiol. 2010;20(8):2005–12.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Liu A, Wijesurendra RS, Liu JM, Greiser A, Jerosch-Herold M, Forfar JC, et al. Gadolinium-free cardiac MR stress T1-mapping to distinguish epicardial from microvascular coronary disease. J Am Coll Cardiol. 2018;71(9):957–68.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Mall G, Huther W, Schneider J, Lundin P, Ritz E. Diffuse intermyocardiocytic fibrosis in uraemic patients. Nephrol Dial Transplant. 1990;5(1):39–44.

    Article  CAS  PubMed  Google Scholar 

  14. Aoki J, Ikari Y, Nakajima H, Mori M, Sugimoto T, Hatori M, et al. Clinical and pathologic characteristics of dilated cardiomyopathy in hemodialysis patients. Kidney Int. 2005;67(1):333–40.

    Article  PubMed  Google Scholar 

  15. Rutherford E, Talle MA, Mangion K, Bell E, Rauhalammi SM, Roditi G, et al. Defining myocardial tissue abnormalities in end-stage renal failure with cardiac magnetic resonance imaging using native T1 mapping. Kidney Int. 2016;90(4):845–52.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Graham-Brown MP, March DS, Churchward DR, Stensel DJ, Singh A, Arnold R, et al. Novel cardiac nuclear magnetic resonance method for noninvasive assessment of myocardial fibrosis in hemodialysis patients. Kidney Int. 2016;90(4):835–44.

    Article  PubMed  Google Scholar 

  17. Graham-Brown M, March D, Churchward D, Young H, Dungey M, Lloyd S, et al. Design and methods of CYCLE-HD: improving cardiovascular health in patients with end stage renal disease using a structured programme of exercise: a randomised control trial. BMC Nephrol. 2016;17(1):69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Kramer CM, Barkhausen J, Flamm SD, Kim RJ, Nagel E. Standardized cardiovascular magnetic resonance (CMR) protocols 2013 update. J Cardiovasc Magn Reson. 2013;15(1):91.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Messroghli DR, Radjenovic A, Kozerke S, Higgins DM, Sivananthan MU, Ridgway JP. Modified look-locker inversion recovery (MOLLI) for high-resolution T1 mapping of the heart. Magn Reson Med. 2004;52(1):141–6.

    Article  PubMed  Google Scholar 

  20. Singh A, Horsfield MA, Bekele S, Khan JN, Greiser A, McCann GP. Myocardial T1 and extracellular volume fraction measurement in asymptomatic patients with aortic stenosis: reproducibility and comparison with age-matched controls. Eur Heart J Cardiovasc Imaging. 2015;16(7):763–70.

    Article  PubMed  Google Scholar 

  21. Graham-Brown MP, Rutherford E, Levelt E, March DS, Churchward DR, Stensel DJ, et al. Native T1 mapping: inter-study, inter-observer and inter-center reproducibility in hemodialysis patients. J Cardiovasc Magn Reson. 2017;19(1):21.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Karamitsos TD, Ntusi NA, Francis JM, Holloway CJ, Myerson SG, Neubauer S. Feasibility and safety of high-dose adenosine perfusion cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2010;12(1):66.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Levelt E, Piechnik SK, Liu A, Wijesurendra RS, Mahmod M, Ariga R, et al. Adenosine stress CMR T1-mapping detects early microvascular dysfunction in patients with type 2 diabetes mellitus without obstructive coronary artery disease. J Cardiovasc Magn Reson. 2017;19(1):81.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Kuijpers D, Prakken NH, Vliegenthart R, van Dijkman PR, van der Harst P, Oudkerk M. Caffeine intake inverts the effect of adenosine on myocardial perfusion during stress as measured by T1 mapping. Int J Card Imaging. 2016;32(10):1545–53.

    Article  Google Scholar 

  25. Van Dijk R, Kuijpers D, Kaandorp T, van Dijkman P, Vliegenthart R, van der Harst P, et al. Effects of caffeine intake prior to stress cardiac magnetic resonance perfusion imaging on regadenoson-versus adenosine-induced hyperemia as measured by T1 mapping. Int J Card Imaging. 2017;33(11):1753–9.

    Article  Google Scholar 

  26. Bland JM, Altman D. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;327(8476):307–10.

    Article  Google Scholar 

  27. Kellman P, Hansen MS. T1-mapping in the heart: accuracy and precision. J Cardiovasc Magn Reson. 2014;16(1):2.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Chih S, Macdonald PS, Feneley MP, Law M, Graham RM, McCrohon JA. Reproducibility of adenosine stress cardiovascular magnetic resonance in multi-vessel symptomatic coronary artery disease. J Cardiovasc Magn Reson. 2010;12(1):42.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Larghat AM, Maredia N, Biglands J, Greenwood JP, Ball SG, Jerosch-Herold M, et al. Reproducibility of first-pass cardiovascular magnetic resonance myocardial perfusion. J Magn Reson Imaging. 2013;37(4):865–74.

    Article  PubMed  Google Scholar 

  30. Brown LA, Onciul SC, Broadbent DA, Johnson K, Fent GJ, Foley JR, et al. Fully automated, inline quantification of myocardial blood flow with cardiovascular magnetic resonance: repeatability of measurements in healthy subjects. J Cardiovasc Magn Reson. 2018;20(1):48.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Likhite D, Suksaranjit P, Adluru G, Hu N, Weng C, Kholmovski E, et al. Interstudy repeatability of self-gated quantitative myocardial perfusion MRI. J Magn Reson Imaging. 2016;43(6):1369–78.

    Article  PubMed  Google Scholar 

  32. Schwartzkopff B, Motz W, Frenzel H, Vogt M, Knauer S, Strauer BE. Structural and functional alterations of the intramyocardial coronary arterioles in patients with arterial hypertension. Circulation. 1993;88(3):993–1003.

    Article  CAS  PubMed  Google Scholar 

  33. Díez J, López B, González A, Querejeta R. Clinical aspects of hypertensive myocardial fibrosis. Curr Opin Cardiol. 2001;16(6):328–35.

    Article  PubMed  Google Scholar 

  34. Salerno M, Beller GA. Noninvasive assessment of myocardial perfusion. Circ Cardiovasc Imaging. 2009;2(5):412–24.

    Article  PubMed  Google Scholar 

  35. Le DE, Jayaweera AR, Wei K, Coggins MP, Lindner JR, Kaul S. Changes in myocardial blood volume over a wide range of coronary driving pressures: role of capillaries beyond the autoregulatory range. Heart. 2004;90(10):1199–205.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Lindner JR, Skyba DM, Goodman NC, Jayaweera AR, Kaul S. Changes in myocardial blood volume with graded coronary stenosis. Am J Phys. 1997 Jan;272(1 Pt 2):H567–75.

    CAS  Google Scholar 

  37. Zorach B, Shaw PW, Bourque J, Kuruvilla S, Balfour PC, Yang Y, et al. Quantitative cardiovascular magnetic resonance perfusion imaging identifies reduced flow reserve in microvascular coronary artery disease. J Cardiovasc Magn Reson. 2018;20(1):14.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Kotecha T, Martinez-Naharro A, Yoowannakul S, Lambe T, Rezk T, Knight DS, et al. Acute changes in cardiac structural and tissue characterisation parameters following haemodialysis measured using cardiovascular magnetic resonance. Sci Rep. 2019;9(1):1–8.

    Article  CAS  Google Scholar 

Download references


This study is part of the research portfolio supported by the NIHR Leicester Biomedical Research Centre and the Leicester Clinical Research Facility based at the University Hospitals of Leicester and the University of Leicester. FEP received a Wolfson Intercalated Award administered by the Royal College of Physicians and a Kidney Research UK Intercalated Degree Award.


This study is independent research arising from a Clinician Scientist Award (Dr James Burton, CS-2013-13-014) supported by the NIHR, and a grant from Kidney Research UK (Research Innovation Grant IN02/2013). This work is also supported by the NIHR Leicester Biomedical Research Centre, and the NIHR Clinical Research Facility. The views expressed in this publication are those of the authors and not necessarily those of the National Health Service, the National Institute for Health Research or the Department of Health.

Author information

Authors and Affiliations



FEP – data collection, analysis, figure preparation, manuscript preparation and revision. GSG – manuscript revision and figure preparation. DSM – data acquisition and manuscript revision. AA – completion of regional wall motion scoring and manuscript revision. KSP – scanning and image acquisition. JW – scanning and image acquisition. GPM – study conception and manuscript revision. JOB – study conception and design, manuscript revision. MGB – study conception, oversight of data collection and acquisition, revision and final approval of manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Matthew PM Graham-Brown.

Ethics declarations

Ethics approval and consent to participate

Patients were recruited as a part of the CYCYLE-HD study (ISRCTN 11299707) [17]. The study was given ethical approval by the NHS Research Ethics Committee East Midlands (Northampton; REC ref.: 14/EM/1190). All participants gave written and informed consent.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix 1

Appendix 1

Table 6

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Poli, F., Gulsin, G., March, D. et al. The reliability and feasibility of non-contrast adenosine stress cardiovascular magnetic resonance T1 mapping in patients on haemodialysis. J Cardiovasc Magn Reson 22, 43 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: