Motion-corrected multiparametric renal arterial spin labelling at 3 T: reproducibility and effect of vasodilator challenge

Objectives We investigated the feasibility and reproducibility of free-breathing motion-corrected multiple inversion time (multi-TI) pulsed renal arterial spin labelling (PASL), with general kinetic model parametric mapping, to simultaneously quantify renal perfusion (RBF), bolus arrival time (BAT) and tissue T1. Methods In a study approved by the Health Research Authority, 12 healthy volunteers (mean age, 27.6 ± 18.5 years; 5 male) gave informed consent for renal imaging at 3 T using multi-TI ASL and conventional single-TI ASL. Glyceryl trinitrate (GTN) was used as a vasodilator challenge in six subjects. Flow-sensitive alternating inversion recovery (FAIR) preparation was used with background suppression and 3D-GRASE (gradient and spin echo) read-out, and images were motion-corrected. Parametric maps of RBF, BAT and T1 were derived for both kidneys. Agreement was assessed using Pearson correlation and Bland-Altman plots. Results Inter-study correlation of whole-kidney RBF was good for both single-TI (r2 = 0.90), and multi-TI ASL (r2 = 0.92). Single-TI ASL gave a higher estimate of whole-kidney RBF compared to multi-TI ASL (mean bias, 29.3 ml/min/100 g; p <0.001). Using multi-TI ASL, the median T1 of renal cortex was shorter than that of medulla (799.6 ms vs 807.1 ms, p = 0.01), and mean whole-kidney BAT was 269.7 ± 56.5 ms. GTN had an effect on systolic blood pressure (p < 0.05) but the change in RBF was not significant. Conclusions Free-breathing multi-TI renal ASL is feasible and reproducible at 3 T, providing simultaneous measurement of renal perfusion, haemodynamic parameters and tissue characteristics at baseline and during pharmacological challenge. Key points • Multiple inversion time arterial spin labelling (ASL) of the kidneys is feasible and reproducible at 3 T. • This approach allows simultaneous mapping of renal perfusion, bolus arrival time and tissue T 1 during free breathing. • This technique enables repeated measures of renal haemodynamic characteristics during pharmacological challenge. Electronic supplementary material The online version of this article (10.1007/s00330-018-5628-3) contains supplementary material, which is available to authorized users.


Introduction
The auto-regulatory feedback mechanisms that maintain renal blood flow (RBF) and glomerular filtration rate (GFR) over a wide range of haemodynamic pathophysiological conditions can be disturbed in metabolic and inflammatory disease states promoting the progression of nephropathy [1][2][3]. Measuring and interpreting altered renal perfusion at baseline and in response to physiological stressors would facilitate assessment of individuals at risk of kidney disease or its complications, and assist in the development of pharmacological interventions [4]. Measurement of renal perfusion currently relies on exogenous administration of contrast medium or clearance of radioisotopes which have limited availability, may be contraindicated when renal function is impaired, and do not allow repeated realtime assessment during pharmacological challenge [5].
Magnetic resonance imaging (MRI) is emerging as a noninvasive approach for assessing renal structure and function including physiological parameters such as tissue perfusion, oxygenation, and water diffusion [6]. Arterial spin labelling (ASL) is a method that uses flowing blood as an endogenous contrast agent that has been widely used in neuro-imaging applications [7], and has also been shown to be feasible in healthy, transplanted and diseased kidneys [8]. ASL uses a radiofrequency (RF) pulse to magnetically label water protons in blood, so that they act as a diffusible tracer. Subtraction of the labelled images from control images allows perfusion maps to be quantified by fitting a kinetic model to the data [9]. In standard single-inversion time (TI) ASL techniques, a single inversion labelling pulse is applied, and a perfusion value is calculated using a simplified model which neglects variations in renal tissue T 1 and makes assumptions concerning the bolus arrival characteristics and blood T 1 [10,11]. As the signal is dependent on the chosen TI, which specifies the interval between labelling and image acquisition, dynamically acquiring data at multiple-TIs would allow analysis using a general kinetic model that accounts for the effects of variable transit delays as well as tissue relation properties [12].
In this study, we investigated the feasibility and reproducibility of a free-breathing motion-corrected multi-TI pulsed renal ASL (PASL) sequence, with general kinetic model parametric mapping, enabling perfusion to be quantified with correction for bolus arrival time (BAT) and renal T 1 . We also assessed its potential to evaluate renal perfusion during pharmacological challenge. For this we used the vasodilator glyceryl trinitrate (GTN) which in animal models leads to a transient fall in systemic blood pressure but maintained kidney perfusion due to renal auto-regulatory mechanisms [13].

Subjects
This single-centre prospective observational and interventional study was approved by the Health Research Authority (17/ EE/0068). Twelve healthy participants were recruited between March and July 2017. Standard published safety contraindications to MRI were applied [14]. No participants had known or suspected cardiovascular or renal disease, and none were pregnant. All participants gave written informed consent.

Study protocol
A single-TI and multi-TI sequence were acquired consecutively in each participant. Reproducibility was assessed by repeating this protocol twice on the same day. In six participants, multi-TI ASL was used to assess the effect of administering a vasodilator. Heart-rate and blood pressure were monitored using an MRI-compatible system (TeslaM3; MIPM, Mammendorf, Germany). Following two consecutive baseline ASL measurements, 800 μg of GTN was administered sublingually. The first post-GTN image was acquired at 5 min, followed by three further measurements at 7-min intervals.

Imaging protocol
Supine imaging was performed on a 3T MAGNETOM Prisma (Siemens Healthcare, Erlangen, Germany) with XR gradients (80 mT/m at 200 T/m/s) using an 18-channel anterior coil and a six-channel posterior coil. Oblique coronal imaging of both kidneys was performed with the volume positioned to exclude the abdominal aorta. A proton density image (M 0 ) was acquired as a reference for quantitative modelling, followed by control and label images. Each ASL acquisition employed a prototype flow sensitive alternating inversion recovery (FAIR) [15] preparation with a three-dimensional gradient and spin echo (3D-GRASE) [16] read-out (Fig. 1). The field of view (FOV) was 64 × 150 × 300 mm with a voxel size of 4 × 4.7 × 4.7 mm. All the k-space data for one TI were acquired at each excitation. The repetition (TR) and echo times (TE) for the single-TI protocol were TR/TE = 3,500/ 20.82 ms, with a receiver bandwidth (BW) of 2,298 Hz/pixel and 40 signal averages, giving a total acquisition time of 4 m 42 s; for the multi-TI protocol, the corresponding parameters were TR/TE = 3,500/16.88, BW = 3,720 Hz/pixel, 4 signal averages and 6 m 33 s. No partial-Fourier encoding, segmented acquisition or parallel imaging acceleration was used in any direction.
Background suppression (BS) was applied, consisting of an initial pre-saturation pulse, immediately after the labelling pulse, followed by a pair of separated adiabatic inversion pulses, each of 15.36 ms duration; this nulled signal from tissues with T 1 of values 230 ms and 460 ms, working as a double-inversion preparation module. The timing of the inversion pulses was chosen such that the sum of the remaining tissue signals was significantly suppressed [17]. The BS pulses were followed by a delay of 30 ms, to avoid negative longitudinal magnetisation. The bolus length (BL) was limited to 800 ms using a train of modulated saturation pulses. A fatsaturation module was employed immediately before image acquisition, and a series of modulated saturation pulses was applied, anterior and posterior to the imaging volume, in order to limit the length of the labelled blood bolus and to reduce the effect of intra-vascular blood flowing into the volume of interest. Control and label images were corrected for motion by performing retrospective 2D elastic registration using proprietary vendor software. Images were acquired during free breathing. For the single-TI protocol, a T 1 value for blood (T 1b ) of 1,250 ms was assumed for protons in water; the single TI was also 1,250 ms, whereas for the multi-TI protocol, there were 14 TIs, starting at 200 ms and increasing by uniform steps of 175 ms up to 2,475 ms. The range of TIs was chosen to sample both the onset of the renal perfusion time-course and the final signal decay.

Parametric mapping
The perfusion weighted image (PWI) was constructed from control and label images, and the parameter maps (motioncorrected M 0 and f for the single-TI protocol; also, BAT for the multi-TI acquisition) were calculated on-line from this, using the prior M 0 image for signal calibration and image co-registration.
Voxel-wise single-TI perfusion values were calculated according to this formula (derived by means of simplifying assumptions from the full Buxton model) [18]: Perfusion values inferred from the multi-TI scans were calculated in accordance with the Buxton general kinetic model [19]: For PASL, Buxton gives a delivery function c(t) which is non-zero within a fixed window, determined by the BAT Δt and the labelling pulse duration τ, modified by an inversion efficiency α and governed by T 1b : The residue (or venous clearance) function r(t) is governed, in turn, by the exchange of protons between water and tissue:  Fig. 1 A Pulse sequence diagram for the multi-TI ASL sequence: radiofrequency (RF) preparation consists of an ASL labelling pulse (red) followed by four slice-selective (SS) pulses (green) on the labelling volume. A pair of background suppression (BS) pulses (light blue) occurs during a series of SS amplitude-modulated saturation pulses (orange) prior to a fat saturation prepulse (pink). The multi-inversion time readout begins with an SS excitation (red), followed by SS refocusing pulses (dark blue) with echo-planar imaging (EPI; grey). B Slice positioning of the imaging volume in renal ASL measurements. Aorta is not incorporated in the imaging volume to ensure proper labelling of blood flow into the region of interest Finally, the relaxation function is controlled by the T 1 of tissue: The new terms introduced here are:

Kidney segmentation
Whole kidneys were manually segmented slice by slice in MATLAB 2016b (MathWorks, Natick, MA, USA). For the avoidance of partial-volume effects, only the central eight sections of each 16-section 3D image were used. In order to allow a consistent comparison between single-TI and multi-TI studies, the cortex was considered as the outer 3 voxels bounded by the renal segmentation. We also demonstrated segmentation of the cortex from the medulla on multi-TI datasets by propagating a mask, obtained by applying a threshold to the T 1 map histogram, to the perfusion map [8]. We adopted a semi-automated approach where an observer interactively set a T 1 threshold while viewing the resulting segmentation on the 3D volume of images.

Statistical analysis
Model fitting and parametric mapping were performed in MATLAB 2016b (MathWorks) and statistical analysis in IBM SPSS version 22 (IBM, Armonk, NY, USA). Categorical variables were expressed as percentages and continuous variables as mean ± standard deviation (SD). The correlation between repeated measurements was assessed using Pearson's coefficient with major axis regression. The agreement between methods was assessed using Bland-Altman plots with 95% confidence intervals (CI) for the mean difference and 95% CI for the limits of agreement. Continuous variables were compared with the Wilcoxon signed-rank test.
The difference in means of repeated observations in the interventional study was assessed with a repeated-measures analysis of variance (ANOVA). In all tests, a p value <0.05 was taken as significant.

Study participant demographics
All 12 subjects successfully completed the protocol. All data sets were analysed and included in the final analysis. There were five men (age, 26 Fig. 6, with the corresponding histogram in Supplementary Fig. 1, and the resulting derived data given in Supplementary Table 1. Using multi-TI ASL, averaged over both kidneys, mean RBF was 188.54 ± 31.91 ml/min/100 g in the whole kidney and 190.72 ± 32.35 ml/min/100 g in the cortex. Inter-study correlation (r 2 ) was 0.91 for whole kidney perfusion and 0.90 for cortical perfusion. The T 1 of the renal cortex was shorter than that of the medulla (799.61 vs 807.11 ms, p = 0.01). Mean BAT in the whole kidney was 269.73 ± 58.43 ms.

GTN challenge
GTN administration was well tolerated by all participants. Mean systolic (p = 0.02) and diastolic blood pressure (p = 0.01) differed significantly between time points, but the effect on RBF (p = 0.33) and BAT (p = 0.13) was not significant (Supplementary Fig. 2). the tissue-blood partition coefficient of water, given as 0.9 ml/g t in-flow or inversion time (between the end of the labelling pulse and the beginning of read-out) t′ a dummy variable, such that 0 ≤ t′ ≤ t

Discussion
Free-breathing multi-TI ASL acquired at 3 T allows simultaneous measurement of RBF, BAT and tissue T 1 using a general kinetic model for parameter fitting. Such an approach obviates the need for the parameter assumptions of single-TI ASL, is well-tolerated in healthy subjects and is highly reproducible. This technique enables repeated measures of renal haemodynamic characteristics during pharmacological challenge and offers a robust method for assessing renal physiological mechanisms over time.
Respiratory motion has posed a significant challenge in renal ASL due to the potential for mis-registration of pixels between temporally-adjacent labelled and control imageswhich is compounded by the longer readouts required for multi-TI ASL. Previous work has relied on acquiring data at multiple inversion times during separate breath-holds at 3 T [20,21], or using prospective respiratory gating at 1.5 T [22,23]. The disadvantage of prospective gating is that it prolongs acquisition time unnecessarily if separation of respiratory phases is not required. Alternatively, sequential breath-hold acquisitions may be up to 24 s [20], leading to a technical failure rate of approximately 8% [24,25]. Whilst affine registration can adjust for bulk motion [26], our non-linear elastic model may be more suited to the non-rigid motion typical of abdominal imaging and enabled a time-efficient free-breathing acquisition for multi-TI ASL.
Our ASL values for cortical renal perfusion are in agreement with those obtained at 3 T using breath-hold multi-TI sequences: 151 ± 37 to 215 ± 65 ml/min/100 g [20,21], but  Fig. 3 Parametric mapping using fitting of multi-TI ASL data to a general kinetic model. Coronal images of both kidneys are shown with M 0 , T 1 , BAT and perfusion maps even in healthy volunteers at this field strength single-TI values range between 229 ± 41 to 327 ± 63 ml/min/100 g [10,11,26]. This variation can be attributed to different acquisition schemes, kinetic models and study populations. In this regard, multi-TI ASL could be advantageous since it takes delivery, clearance and relaxation effects fully into account to better represent the dynamics of renal blood flow. The longer T 1 of tissues at high field strength also allows for longer postlabelling delays and therefore higher signal in PWI [27]. Previous studies using multiple post-labelling delays have demonstrated good correlation between ASL-derived parameters and dynamic contrast-enhancement (DCE) [20], although more modest agreement was observed with 99m Tc-MAG3 scintigraphy-suspected to be a consequence of respiratory motion [25]. Variations in renal ATT have been observed between young and old age groups, and multi-TI ASL enables patient-specific modelling of RBF [25]. Beyond correction of absolute RBF estimates, the value of BAT as a parameter of renal haemodynamic state is unknown and worthy of further investigation. Kidney interstitial fibrosis is a histological finding that represents a common final pathway in many renal diseases and T 1 values correlate well with fibrosis and inflammation in both animal models and humans [28]. Multi-TI ASL enables renal T 1 values to be modelled as part of a multi-parametric assessment of renal function and avoids the need to acquire and co-  Fig. 4 Scatterplot with major axis (MA) regression of inter-study perfusion measurement of cortex and whole kidney using single-TI and multi-TI ASL register a separate renal T 1 mapping sequence [11]. T 1 maps that are spatially aligned with other parametric maps enable an objective and consistent segmentation of cortex from medulla, which is important as altered cortico-medullary physiology is a characteristic of renal insufficiency [29]. Our T 1 values for renal cortex and medulla were lower than previously-reported data at 3 T, using an inversion recovery method, with less differentiation between these tissues. [30] A plausible explanation would be variation in effective T 1 * due to the influence of different read-out pulses on recovery of longitudinal magnetisation. Recently, a T 1 mapping phantom has been proposed for cardiac MRI to ensure measurement stability over time at individual sites across vendors, software versions and imaging sequences [31], which could be adapted to the clinically relevant ranges overserved in the kidney. Otherwise, if quantitative values are reported reference ranges obtained from healthy volunteers at each imaging centre may be necessary [32]. Previous studies have shown the sensitivity of ASL to map dynamic changes in renal perfusion during administration of adenosine [33], protein load [34] and furosemide [35]. GTN is known to be a potent venodilator with the acute effect mainly being a reduction of systolic blood pressure due to an increase in arterial compliance and a decrease in venous return resulting in reflex tachycardia [36]. In contrast to other smooth muscle vasodilators, renal autoregulation maintains stable arterial perfusion during GTN administration in animal models [13]. We observed a typical haemodynamic response to sublingual GTN but the effect on RBF was not significant consistent with stable blood flow despite a fluctuation in perfusion pressure. This experiment shows as a proof-of-principle how simultaneous dynamic assessment of multiple renal perfusion and tissue characteristics (i.e. BAT, RBF, T 1 ) can be readily performed during multi-TI renal ASL.
Our study has limitations. We used a relatively small sample of young healthy adults and this may not be representative of broader populations and those with renal impairment. However, the acquisition technique is not dependent on breath-hold capability and would be expected to be well-

[ms]
A B Fig. 6 Parametric mapping of multi-TI ASL data enables cortex to be differentiated from medulla by the difference in T 1 (a), and for this to be propagated to the perfusion map for anatomic segmentation (b). The corresponding histogram and threshold is shown in Supplementary Fig. 2 tolerated in patient groups. We did not have a Bgold standardô f renal perfusion to compare with as administration of either gadolinium-based contrast agents or radio-tracers to healthy subjects is avoided if possible. Future work on quantitative evaluation of ASL could be validated against phantoms that simulate capillary blood flow (Qasper; Gold Standard Phantoms, London, UK).
In conclusion, free-breathing multi-TI renal ASL is a reproducible and well-tolerated approach for simultaneous parametric mapping of kidney perfusion and tissue characteristics at 3 T.
Funding The study was funded by the Medical Research Council; and the National Institute for Health Research (NIHR) Biomedical Research Centre based at Imperial College Healthcare NHS Trust and Imperial College London.

Compliance with ethical standards
Guarantor The scientific guarantor of this publication is Declan O'Regan.

Conflict of interest
The authors had a research agreement with Siemens for use of this prototype sequence.
Statistics and biometry No complex statistical methods were necessary for this paper.
Informed consent Written informed consent was obtained from all subjects in this study.
Ethical approval Healthy Research Authority approval was obtained.

Methodology
• prospective • observational and experimental • performed at one institution Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.