Real-time assessment of right and left ventricular volumes and function in children using high spatiotemporal resolution spiral bSSFP with compressed sensing

Background Real-time cardiovascular magnetic resonance (CMR) assessment of ventricular volumes and function enables data acquisition during free-breathing. The requirement for high spatiotemporal resolution in children necessitates the use of highly accelerated imaging techniques. Methods A novel real-time balanced steady state free precession (bSSFP) spiral sequence reconstructed using Compressed Sensing (CS) was prospectively validated against the breath-hold clinical standard for assessment of ventricular volumes in 60 children with congenital heart disease. Qualitative image scoring, quantitative image quality, as well as evaluation of biventricular volumes was performed. Standard BH and real-time measures were compared using the paired t-test and agreement for volumetric measures were evaluated using Bland Altman analysis. Results Acquisition time for the entire short axis stack (~ 13 slices) using the spiral real-time technique was ~ 20 s, compared to ~ 348 s for the standard breath hold technique. Qualitative scores reflected more residual aliasing artefact (p < 0.001) and lower edge definition (p < 0.001) in spiral real-time images than standard breath hold images, with lower quantitative edge sharpness and estimates of image contrast (p < 0.001). There was a small but statistically significant (p < 0.05) overestimation of left ventricular (LV) end-systolic volume (1.0 ± 3.5 mL), and underestimation of LV end-diastolic volume (− 1.7 ± 4.6 mL), LV stroke volume (− 2.6 ± 4.8 mL) and LV ejection fraction (− 1.5 ± 3.0%) using the real-time technique. We also observed a small underestimation of right ventricular stroke volume (− 1.8 ± 4.9 mL) and ejection fraction (− 1.4 ± 3.7%) using the real-time imaging technique. No difference in inter-observer or intra-observer variability were observed between the BH and real-time sequences. Conclusions Real-time bSSFP imaging using spiral trajectories combined with a compressed sensing reconstruction showed good agreement for quantification of biventricular metrics in children with heart disease, despite slightly lower image quality. This technique holds the potential for free breathing data acquisition, with significantly shorter scan times in children. Electronic supplementary material The online version of this article (10.1186/s12968-018-0500-9) contains supplementary material, which is available to authorized users.


Background
Evaluation of biventricular volumes and function is vital in the investigation of children with congenital heart disease. The clinical standard method of evaluating ventricular volumes is cardiovascular magnetic resonance (CMR)specifically, multi-slice cardiac gated balanced steady state free precession (bSSFP) cine imaging [1]. This technique is now routinely used in children and has an important role to play in disease management [2]. Unfortunately, it does require multiple breath-holds (BH) and this can be difficult in young children with heart disease. An alternative approach is real-time CMR in which each k-space is acquired in a single shot manner, rather than over several heart beats. The benefit of real-time CMR is that it can be acquired during free breathing and without cardiac gating. However, this comes at the cost of lower spatial and temporal resolution, which can affect accuracy [3].
One solution is to leverage accelerated imaging techniques (data undersampling or more efficient k-space filling) to improve resolution. Recently, a combination of undersampled radial k-space filling and compressed sensing (CS) has been shown to enable relatively high resolution real-time imaging [4,5]. Studies in adults have demonstrated good agreement between these radial real-time techniques and the breath-hold clinical standard for assessment of ventricular volumes [6][7][8][9]. However, imaging children requires much higher spatial and temporal resolution than available using conventional radial CS methods. Hence greater acceleration is needed to translate these techniques into the paediatric population. One possibility is to combine CS with more efficient spiral k-space filling. Spiral trajectories have previously been used to accelerate real-time phase contrast CMR with some success [10]. They can also be combined with golden angle spacing [11,12] to enforce incoherent aliasing necessary for CS reconstruction.
In this study, we implemented a novel real-time bSSFP spiral sequence reconstructed using CS. The purpose of this study was to validate this real-time sequence against the BH clinical standard for assessment of biventricular volumes in children with heart disease.

Study population
Between December 2017 and February 2018, 60 consecutive children referred for CMR and consented for additional research scans were included into this study. The local committee of the UK National Research Ethics Service approved the study (06/Q0508/124), and written consent was obtained from all subjects/guardians. Patients were excluded from this study if: a) there was a requirement for general anesthetic (15 children excluded during the recruitment period), b) the subject had an irregular heart rate or difficulty performing breath-holds (11 children excluded), or c) the subject had single ventricular anatomy which was excluded to ensure that biventricular volumes could be evaluated in all patients (4 children excluded).

Imaging protocol
All imaging was performed on a 1.5 Tesla CMR scanner (Avanto, Siemens Healthineers, Erlangen, Germany) with maximum gradient amplitude of 40 mT/m and a maximum slew rate of 180 T/m/s. Signal detection was accomplished using two spine coils incorporated into the scan table and arranged in head-foot manner, and one anteriorly placed body-matrix coil (giving a total of 12 coil elements). A vector electrocardiographic system was used for cardiac gating.

Standard volumetric assessment
Standard cardiac-gated ventricular volume assessment was performed using a multi-slice retrospectively cardiac gated bSSFP Cartesian sequence. The temporal resolution was 29.5 ms and the in-plane spatial resolution was 1.5 × 1.5 mm (full sequence parameters are shown in Table 1). The imaging plane was in the ventricular short axis (SAX), planned using right ventricular (RV) and left ventricular (LV) long axes and 4-chamber images. Sufficient contiguous slices were acquired in the short axis to ensure coverage of the whole ventricle (1 3 ± 2 slices, range: 10 to 17). Each slice was acquired in a separate breath-hold, each lasting~5.5 ± 1.1 s (range: 3.6 to 8.8 s).

Real-time volumetric assessment
Real-time bSSFP imaging was performed using a 2D multi-slice uniform density spiral sequence. The spiral trajectories were designed using the method described by Hargreaves [13], assuming a field-of-view of 350 mm, a spatial resolution of~1.7 × 1.7 mm and 72 regularly spaced spiral interleaves for complete filling of k-space. A zeroth moment rewinder gradient was added at the end of the spiral readout to maintain bSSFP coherence (see Fig. 1 for pulse sequence diagram). In this study we did not use bipolar first moment rewinder gradients, due to the associated increase in TR [14]. Consecutive interleaves were rotated by the tiny golden angle (tGA,~47.26 o ) throughout data acquisition. This enforced incoherent aliasing required for the CS reconstruction, whilst being associated with lower eddy currents than traditional golden angle sampling patterns [15]. Each real-time frame was formed by combining 9 consecutive spiral interleaves (acceleration factor = 8), resulting in a temporal resolution of~30 ms. The full sequence parameters are shown in Table 1.
Each slice was acquired over two R-R intervals. The first was used to reach the steady state and the second was used for acquisition of data. Synchronisation of acquisition of data to the start of the R-wave for each slice also ensured the real-time images were temporally aligned for measurement purposes. After acquisition of each slice, the imaging plane was moved to the next contiguous slice. The imaging plane and the number of slices was the same as the clinical standard BH sequence. All real-time image acquisition was performed during free breathing.
The CS reconstruction was performed on an external GPU equipped computer (Tesla K40c, NVIDIA Corportation, Santa Clara, California, USA) with on-line communication to the native reconstructor [16,17]. Our implementation of CS solved a set of non-linear equations (representing the imaging process) through minimization of the following cost function: The first term enforces data consistency -where ρ is image data, E is the encoding matrix (the multi-coil non-uniform Fourier transform (FT) operator) and y is the acquired k-space data. The second term enforces sparse results through L 1 norm regularization. A temporal finite difference operator (TV t ), was used transform the image data into a sparser representation. This is similar to the previously described iterative Golden-angle Radial Sparse Parallel (GRASP) reconstruction [18]. The optimization was performed using a non-linear Conjugate Gradient algorithm, with 35 iterations. The regularization parameter (λ) were selected empirically on the first 10 patients, to optimize artefact removal without temporal blurring and was set at 0.05 for all patients. As the data for each slice is independent, reconstruction started immediately after the last phase for each slice. Coil sensitivity information was calculated from the time-averaged data (from each slice) [19].

Image processing and scoring
Evaluation of ventricular volumes and qualitative image scoring was performed by three clinical CMR specialists  (OT -12 years experience, MH -11 years experience, and KM -5 years experience). For both image quality scoring and assessment of LV and RV volumes, each clinician was the primary observer for 20 unique cases, of which 10 were re-evaluated to assess intra-observer variability. In addition, each observer assessed 10 cases from a different primary observer (5 cases from each of the other observers), to evaluate inter-observer variability. Thus, each observer scored and processed 40 cases. Overall 30 cases were used to evaluate intra-observer variability and the other 30 cases used to evaluate inter-observer variability.

Image quality
The mid-ventricular short-axis cine loops from each technique were scored on a 5-point Likert scale in three categories: sharpness of endocardial border (1 = nondiagnostic, 2 = poor, 3 = adequate, 4 = good, 5 = excellent), temporal fidelity (or blurring) of wall motion (1 = nondiagnostic, 2 = poor, 3 = adequate, 4 = good, 5 = excellent) and residual aliasing artefacts, which appear as spiral streaks across image (1 = non-diagnostic, 2 = severe, 3 = moderate, 4 = mild, 5 = minimal). Each observer scored the patient data that they had previously segmented (N = 40 cases each), at a separate session. All loops were presented in a random manner using a custom-built python application, with the observers blinded to diagnosis, patient number and type of sequence.

Ventricular function
Quantification of LV and RV volumes was performed in a similar manner for each technique using the OsiriX open source DICOM viewing platform (OsiriX v.9.0, OsiriX Foundation, Geneva, Switzerland) [20]. Firstly, the end-diastolic and end-systolic phases were identified for each ventricle through visual inspection of the midventricular cine. The endocardial borders of all slices at end systole and diastole were then traced manually (including papillary muscles and trabeculation in the myocardial mass). This allowed calculation of end-diastolic volume (EDV) and end-systolic volume (ESV). Stroke volume (SV) was obtained by subtracting ESV from EDV and ejection fraction (EF) = SV/EDV × 100. In addition, LV epicardial borders were traced in end-diastole and combined with endocardial borders to obtain LV mass. Observers were presented with each anonymized volume (including repeated volumes for intra-observer variability) in a random order, blinded to diagnosis, patient number and type of sequence.

Quantitative image assessment
Calculation of signal-to noise ratio (SNR) and contrastto-noise ratio (CNR) in images reconstructed using CS is nontrivial due to the uneven distribution of noise.
Therefore, we calculated blood pool-to-myocardial signal intensity ratio [3] as a quantitative measure of image contrast, which is important when segmenting data. In all patients, the average blood pool and myocardial signal intensity values were measured in the end diastolic volume (from the LV endocardial and LV epicardial borders drawn for volumetric calculation, across the entire volume). The image contrast equalled blood pool signal intensity divided by myocardial signal intensity. Edge sharpness of the mid-ventricular short-axis slice was evaluated by measuring the maximum gradient of the normalized pixel intensities across the border of the septum, as previously described [21]. To reduce noise, which results in artificially high gradients (representing sharp edges), the pixel intensities were fit to a tenth order polynomial, before differentiation. Edge sharpness was calculated in six positions across the septum, for all cardiac phases, and the average value was used for comparison.

Statistics
Statistical analyses were performed by using STATA software (STATA SE, v.14.2, Stata Corporation, College Station, Texas, USA). All results are expressed as the mean ± standard deviation. Mean ventricular volumes, function and mass measured using the standard BH and spiral real-time techniques were compared using a paired t-test. For assessment of agreement of ventricular volumes and function, the standard BH data was used as the clinical standard for Bland-Altman analysis. Association of the percentage differences between the real-time and BH data against both heart rate and body surface area were tested using Spearman's rank correlation coefficient. Inter and intra-observer variability was assessed using the coefficient of variation (CoV) calculated using the logarithmic method described by Bland and Altman [22]. In addition, inter-observer and intra-observer variability were also assessed using one-way intraclass correlation coefficient (ICC). Both CoV and ICC were displayed with their 95% confidence intervals. Qualitative and quantitative image scores were also compared using the paired t-test. This was done as previous work has shown that the paired t-test has a lower Type II error rate compared to non-parametric tests for Likert scale data [23]. It is therefore more likely to detect differences between the two techniques. A p-value of less than 0.05 indicated a significant difference.

Demographics and feasibility
The mean age of population was 13.6 ± 2.7 years (median: 13.5, range: 7.0-18.3) and 33 (55%) were female. The mean heart rate was 82.3 ± 15.9 bpm (median: 82, range: 50-114). The full demographic information and patient diagnoses are shown in Table 2. Real-time and standard breath hold imaging was successfully performed in all subjects. The total time acquisition time for the full breath hold stack was 348 ± 79 s (median: 335, range: 239-573) compared to 19.8 ± 5.8 s for the acquisition of the real-time spiral stack (median: 19 s, range: 11-36). The online CS reconstruction for each real-time slice (all phases) was~1.5 s, with an add-itional~3.5 s per slice to send the data to the external computer and receive the reconstructed series back to the scanner database. In our buffered implementation, this resulted in all phases from all slices being available for view-ing~28 s after the end of the acquisition.

Image quality
Representative images are shown in Figs. 2 and 3, and the corresponding movies can be seen in Additional files 1, 2, 3, and 4. These demonstrate that the spiral real-time images are of good quality. However, Fig. 4 shows examples of cases with residual aliasing artefacts and reduced edge definition (corresponding movies can be seen in Additional file 5 and Additional file 6), which is reflected in the qualitative image scores. Specifically, spiral real-time images had more residual aliasing artefacts than standard breath hold images (artefact score: 3.8 ± 0.1 vs 4.8 ± 0.1, p < 0.001) and lower edge definition (edge score: 4.4 ± 0.7 vs 4.8 ± 0.5, p < 0.001). There was a trend towards poorer motion fidelity (motion score: 4.6 ± 0.1 vs 4.8 ± 0.1, p = 0.06). There were no significant differences between qualitative image scores for the different observers (p > 0.1).
Quantitative edge sharpness was also significantly lower (p < 0.001) for the spiral real-time images compared to the standard BH images (0.68 ± 0.19 mm − 1 vs 0.53 ± 0.16 mm − 1 , respectively). Additionally, quantitative estimates of image contrast (blood pool-to-myocardial signal intensity ratio) were significantly lower (p < 0.001) for the spiral real-time images compared to the standard BH images (2.7 ± 0.4 vs. 3.2 ± 0.4, respectively).

Ventricular volume quantification
Ventricular metrics measured using spiral real-time and clinical standard breath hold imaging are shown in Table  3. Bland-Altman plots for LV metrics are shown in Figs. 5 and 6. There was a small (1.7 mL) but significant   5%). There was also a small but significant overestimation of LV mass (3.1 g). Nevertheless, the limits of agreement between real-time and breath hold data for all LV metrics were narrow. Body surface area correlated with the percentage difference between real-time and breath hold derived LV mass (rho = − 0.27, p = 0.041). Differences for all other LV metrics were not associated (p > 0.05) with either heart rate or body surface area. Bland-Altman plots for RV metrics are shown in Figs. 6 and 7. There were no significant differences in RV EDV or ESV, however there was a small underestimation of RV SV (1.8 mL) and RV EF (1.4%). However, the limits of agreement between real-time and BH data for all RV metrics were narrow. Differences in RV metrics were not significantly (p > 0.05) associated with either heart rate or body surface area.
The inter-observer CoV's and ICC's are shown in Table 4 and intra-observer CoV's and ICC's are shown in Table 5. The largely overlapping confidence intervals demonstrated that there were no significant differences in inter-observer and intra-observer variability between spiral real-time and standard breath hold derived RV and LV metrics.

Discussion
The main findings of this study were that: i) real-time bSSFP imaging using tiny golden-angle spiral trajectories combined with a CS reconstruction was feasible, ii) the image quality of the real-time technique was slightly lower than the standard breath hold technique, iii) there was good agreement for quantification of both LV and RV metrics between the spiral real-time and breath hold clinical standard techniques.

Spiral bSSFP sequence with CS reconstruction
Real-time CMR is particularly useful when imaging children, as it can be performed quickly and without breath holds. In this study, we implemented a novel real-time sequence that leveraged both efficient spiral k-space filling and CS reconstruction. Spiral k-space filling has previously  been used to successfully accelerate spoiled gradient echo sequences (i.e. phase contrast CMR [24] and CMR angiography [25]). However, there are some problems associated with combining spiral trajectories with bSSFP readouts [14,26,27]. These mainly relate to off-resonance effects and longer repetitions times. Previous studies have shown that these effects can be partly mitigated using both zeroth and first order moment rewinders to balance the imaging gradients [14,28]. The problem with this approach is that it significantly increases repetition time and lowers temporal resolution. Therefore, we chose to use a shorter conventional zeroth order moment rewinder, combined with a relatively short readout (1.86 ms) to reduce off-resonance effects. However, this increases the number of interleaves required to fully sample k-space. Consequently, high acceleration factors were required to ensure adequate temporal resolution and we utilized CS to reconstruct artefact-free images. Our CS implementation shares many characteristics with the previously validated GRASP reconstruction for radial trajectories [18]. Specifically, we utilized tiny golden angle spacing to produce temporal incoherency and penalized temporal finite differences in the reconstruction. Combining spiral bSSFP with CS resulted in a real-time sequence with only slightly lower spatial resolution, and similar temporal resolution, to clinical standard breath hold imaging.

Clinical utility
Clinical utility depends on the ability of the reader to process the images to accurately measure ventricular volumes. In this study, images were processed by three independent CMR specialists (not involved in the development of the sequence) This is a more 'real-world' test of the sequence, in contrast to many studies in which a single observer (often involved in the development of the technology) is used.
We demonstrated that spiral real-time images were of diagnostic quality. However, they did have slightly lower edge definition and greater amount of residual aliasing artefact. This is not surprising as the real-time images had lower spatial resolution and were acquired with significant acceleration. Interestingly, although CS with temporal finite difference sparsity can result in temporal blurring, there was only a trend towards reduced motion fidelity.
Irrespective of differences in image quality, from a clinical point it is the ability to accurately measure ventricular volumes that is paramount. In this study, we found that LV EDV was smaller and LV ESV was bigger with real-time imaging compared to standard BH imaging. We also found similar differences in RV ESV, although this did not reach statistical significance. These biases were probably due to inaccurate segmentation of the real-time images, secondary to lower edge sharpness and myocardial blood pool contrast. In diastole, we believe these inaccuracies resulted in blood pool being included in the myocardial mass and EDV being underestimated. We believe the converse occurred in systole, with myocardial mass being included in the blood pool and ESV being overestimated.  An important question in children, is whether the differences between real-time and standard breath hold imaging were associated with heart rate or body size. In this study, the only significant association we found was between percentage overestimation of LV mass and decreasing BSA. This was probably due to errors in epicardial segmentation in smaller children, possibly due to factors such as less pericardial fat and geometric considerations. Nevertheless, the biases were all below 5% and would be expected to have minimal effect on clinical decision making. More importantly, the limits of agreement were narrow compared to previously described real-time methods (including those that utilize CS) [6,29]. We believe that this is due to the high spatial and temporal resolution of our spiral real-time sequence, which was only slightly lower than the standard breath hold sequence. In this study, we also examined inter-observer and intra-observer variability (half the populations tested for each). For the standard breath hold sequence the coefficients of variation and ICC's were similar to previously published data [30], with RV data having greater variability than LV data. This is to be expected due to the greater difficulty in segmenting the more trabeculated and complex shaped RV. Importantly, the inter-observer and intra-observer variability of real-time sequences were similar to standard breath hold imaging. This is an important finding, as demonstrating reliability is vital for clinical translation. We believe that our findings show that spiral real-time imaging can successfully be used for assessment of paediatric heart disease. The main benefit over standard breath hold imaging is that it can be acquired during free breathing. Another important secondary benefit is that a whole ventricular stack can be acquired in approximately 20 s, as opposed to almost 6 min for standard breath hold imaging. Thus, this technique holds the potential to significantly shorten CMR scan times in children.

Limitations
The main limitation of this study was the minimum age in the population was 7 years, which limits the translatability of our findings to younger children and infants. The reason for this age limit was that younger children and infants would have found it difficult to perform the standard breath hold sequence necessary for validation. Nevertheless, future studies concentrating on optimization and validation of this technique for infants would be useful. Another limitation was that we did not assess RV mass in this study. Several studies have suggested that RV mass may be an important predictive marker in certain diseases [31]. However, RV mass is difficult to accurately quantify in patients without RV hypertrophy, even when using standard breath hold imaging. Therefore, we did not include it this study, although it could be included in future studies concentrating on patient with RV disease. It should also be noted that the real-time imaging was  performed during free breathing and physiological variation may account for some discrepancies with standard breath hold imaging. However, the differences probably remain small because the slices in the real-time short axis stack are all acquired at in different points in the respiratory cycle. Thus, any discrepancies may 'average out' when summed across the whole ventricle. In terms of intra-observer variability, there was no prescribed time delay between processing of data, as volumes (including repeated volumes) were presented data in a randomized manner. However, both the standard breath hold and real-time data was treated in the same manner and the variability measures were in keeping with previously published data. Thus, we believe the similar intra-observer variability of the real-time and standard breath hold data does reflect the clinical reality. Finally, one disadvantage of compressed sensing is long reconstruction times. In this study, we used an online GPU based reconstruction that fed images back to the scanner in a clinically acceptable time frame (< 30 s for the entire short axis stack). Nevertheless, more work is required so that reconstruction times can be reduced to that of standard Cartesian imaging.

Conclusion
In conclusion, we implemented a novel spiral real-time bSSFP sequence with CS reconstruction. We showed good agreement between biventricular volumes measured using spiral real-time and standard breath hold techniques in children with heart disease. Thus, we believe that this technique could be used successfully in the paediatric population.