1. Introduction
Magnetic resonance cholangiopancreatography (MRCP) is a non-invasive imaging method used in everyday clinical practice to assess anatomical features and abnormalities of the intrahepatic and extrahepatic bile ducts and of the pancreatic ducts [
1,
2,
3,
4,
5]. MRCP does not involve the risks associated with endoscopic retrograde cholangiopancreatography (ERCP) such as acute pancreatitis, bowel perforation, infections, and bleeding [
6]. Recent years have witnessed the development of three-dimensional (3D) imaging, which provides better image quality and greater diagnostic confidence compared to two-dimensional (2D) imaging [
7].
Two 3D MRCP breathing management methods are currently available. The free-breathing (FB) method relies on respiratory gating has well-established diagnostic performance and spatial resolution characteristics when used to evaluate diseases of the bile ducts and pancreas [
2]. However, this method requires quite long acquisition time. The other method, in which the images are acquired during a single breath-hold (BH), is being actively developed.
Research into means of improving MRCP sequences has several goals, of which the most important is reduction in the acquisition time in order to limit motion artifacts, notably those due to breathing, which are particularly challenging when imaging the abdomen. A shorter acquisition time also improves the comfort of the patient, who must remain immobile throughout the acquisition. Another advantage is the imaging of a greater number of patients during a given time, which decreases costs and wait-list times.
In practice, the acquisition time and diagnostic quality of FB MRCP with respiratory triggering (RT) are variable and difficult to predict. Only part of the k-space is acquired during each breathing cycle. Diagnostic performance may be adversely affected by irregular breathing due, for instance, to abdominal pain or failure of the machine to detect breaths. However, irregular breathing can also occur in patients with a World Health Organization (WHO) performance status of 0 for unknown reasons [
8]. FB with RT acquisition times of 7 min [
9], and 6 min [
10], have been reported. Longer acquisition times are often associated with poorer image quality [
9,
11].
3D MRCP is obtained using a T2-weighted fast spin-echo sequence with the variable-flip-angle technique, such as SPACE (Sampling Perfection with Application-optimized Contrast using different flip-angle Evolutions), CUBE, and VISTA (Volume Isotropic Turbo spin echo Acquisition). The high flip angles of radiofrequency pulses for the lines near the center of the Fourier space ensure a good contrast-to-noise ratio (CNR), while the lower angles for the lines at the periphery of the Fourier space provide good spatial resolution. The clinical relevance of the SPACE technique for MRCP has been demonstrated [
12].
Methods that have been evaluated to decrease the acquisition time include parallel imaging (PI) and compressed sensing (CS), which combines data undersampling with iterative reconstruction. PI has entered the mainstream of clinical practice. The decrease in acquisition time is achieved by undersampling the k-space; for instance, by skipping every other line. The main drawbacks are the moderate acceleration factors in the 2–4 range and a decrease in the CNR [
13,
14]. CS has undergone considerable development in recent years as a method for imaging moving targets, such as the heart and abdominal organs [
15,
16,
17,
18,
19]. CS relies on the sparsity of redundant data within standard magnetic resonance images. Random undersampling of the compressible representation of an image is performed, and the image is then restored by iterative reconstruction [
20].
Several studies have evaluated the clinical feasibility of CS with FB or BH acquisition, with interesting results that favored BH sequences [
9,
21,
22,
23]. A prototype CS-BH-SPACE MRCP sequence produced by Siemens Healthcare (Erlangen, Germany) seeks to eliminate motion artifacts during a BH.
The objective of this study was to prospectively compare two MRCP sequences at 3T, namely, the conventional 3D RT-SPACE sequence (designated RT-MRCP hereafter) and the prototype 3D CS-BH-SPACE sequence (designated CS-BH-MRCP hereafter), in terms of qualitative and quantitative image quality and radiologist’s diagnostic confidence for detecting common bile duct (CBD) lithiasis, biliary anastomosis stenosis in liver-transplant recipients, and communication of pancreatic cyst with the main pancreatic duct (MPD), the latter allowing the non-invasive diagnosis of branch-duct intraductal papillary mucinous neoplasm of the pancreas (BD-IPMN).
4. Discussion
In our study, based on quantitative and qualitative evaluations, the CS-BH-MRCP sequence was preferred over the conventional RT-MRCP sequence in terms of image quality at 3T with the advantage of a much shorter acquisition time. The CBD to PBT CR was significantly better. The CS-BH-MRCP sequence demonstrated significantly better overall image quality, fewer artefacts, better background noise suppression, a better visualization of the distal CBD, of the right and left primary IHBD and the distal and proximal parts of the MPD. No significant difference was found regarding the following qualitative criteria: the visualization of the proximal CBD, cystic duct confluence, secondary IHBDs, the central part of the MPD, the entire biliary system and entire MPD.
Overall, our results are consistent with previous reports. Three studies found a significantly better overall image quality with the CS-BH-MRCP sequence at 3T [
9,
24,
31]. Studies also showed a significantly better visualization of the CBD [
9], the primary IHBD [
9,
24], and the cystic duct [
24], with the CS-BH-MRCP sequence at 3T. Nevertheless, the literature seems to show a discrepancy concerning the visualization of the MDP. With 200 patients scanned at 3T, the study of Blaise et al. showed a significantly better visualization of the MPD [
24], although no segmental analysis was performed. This is in contrast with Zhu et al.’s prospective study including 80 patients, which found a worse visualization of the MPD with a CS-BH-MRCP sequence in comparison with a conventional NT-MRCP sequence at 3T [
10]. CS-BH-MRCP had thus lower diagnostic sensitivity [
10]. This finding prompted the same investigators to conduct another study evaluating a modified CS-BH-MRCP sequence with a smaller field of view (FOV) and higher spatial resolution that achieved better visualization of the MPD and secondary IHBDs than the “original” CS-BH-MRCP. This modified protocol also showed higher sensitivity for detecting pancreatic duct abnormalities [
11]. Another optimized CS-BH-MRCP sequence at 3T with decreased accelerator factor and a reduced FOV and matrix without changes in spatial resolution, was proposed by Song et al., and demonstrated comparable or even better image quality than conventional MRCP [
32]. Overall, the MPD was also better visualized with the optimized sequence [
32].
Although our study’s results and previous reports suggested the superiority of the CS-BH-MRCP sequence at 3T, the same was not observed at 1.5T in some studies [
24,
31]. Taron et al.’s study suggested a better overall image quality with the conventional NT-MRCP at 1.5T, although the results were not significant [
31]. At 1.5T, in Blaise et al.’s study, the conventional RT-MRCP acquisition showed a significant superior overall image quality with better visualization of the biliopancreatic ducts, whereas only sharpness was improved with BH-CS-MRCP [
24]. In a recent study, a short single BH CS-MRCP sequence, that allowed a reduced acquisition time of 8 s, demonstrated higher scores for image quality, duct sharpness and duct visualization than the conventional NT-MRCP, a CS-NT-MRCP, and a long single BH CS-MRCP (acquisition time of 17 s) sequences, the results being not always significant for all criteria and sequence to sequence comparison [
33]. This highlights the potential superiority of the CS-BH-MRCP sequence, even at 1.5T and with an even shorter acquisition time.
Unlike the proximal and distal MDP that were more clearly visualized with the CS-BH-MRCP, no difference was found concerning the central MPD. The same was observed for the secondary IHBDs. It might be partially explained by the lower number of coronal slices acquired with the CS-BH-MRCP than with the RT-MRCP (64 vs. 120 slices), resulting in a smaller acquisition volume. The central MPD and secondary IHBDs were indeed less frequently imaged with the CS-BH-MRCP, resulting in a bad duct visualization score of 1. Most of the patients were referred for suspected choledocholithiasis, the field of view was thus most likely centered on the CBD with less care being taken to cover the other pancreato-biliary ducts. Great care is, therefore, required when choosing the 3D imaging volume position most appropriate for the suspected diagnosis. The radiologic technologist must also ensure that the acquisition covers as many ducts as technically possible.
A significant difference in the CBD to periductal tissues CR between the two sequences that favored the CS-CH-MRCP sequence was found in our study. Although the CBD to periductal CR values were similar to values previously reported by Seo et al., i.e., 0.92 ± 0.03 for MRCP with PI and 0.91 ± 0.03 for MRCP with PI and CS, in their study, the significant difference in CR favored the MRCP sequence without CS [
23]. Note that, unlike our study, both sequences were acquired using free-breathing navigator-triggered method and the two sequences acquisition parameters differed only for the acceleration factor and repetition time [
23]. In Song et al.’s study, a significantly better CBD to PTB tissues CR was achieved using an optimized BH-CS-MRCP sequence compared to conventional MRCP (0.99 ± 0.01 versus 0.94 ± 0.04,
p < 0.001), with slightly higher CR values compared to our study for both sequence [
32].
In our study, 3D MRCP with CS was successfully acquired during a 17 s BH in 68 patients, i.e., 97% of the original cohort of 70 patients. This result was obtained although the patients had required hospital admission and exhibited multiple comorbidities likely to cause greater difficulty with maintaining a BH compared to the general population. Our study cohort was thus representative of everyday practice. Similar success rates have been reported [
10,
11]. The 10 min period needed for image reconstruction precluded immediate evaluation of image quality with repeated acquisition or the performance of additional sequences if needed. This point is a limitation to the use of the BH sequence instead of the free-breathing sequence. We used the reconstruction time to prepare the next patient. However, with the recently marketed latest version of the MRI machine software, the reconstruction time is only 20 s, allowing for the fast evaluation of the image quality of the acquisition and its repetition if necessary.
In the present study, radiologist diagnostic confidence was significantly better with the CS-BH-MRCP sequence that is certainly linked to the image quality. These results tend to support the use of the CS-BH-MRCP sequence for the diagnosis of biliary and ductal pancreatic diseases. Image quality indeed needs to be good enough in order to make a diagnosis with a very high degree of certainty. However, according to the senior reader, image quality was optimal to ensure complete confidence in the diagnosis (diagnostic confidence score of 3) in only 11.8% of patients with the RT-MRCP sequence and 26.5% with the CS-BH-MRCP sequence. Moderate confidence (diagnostic confidence score of 2) was achieved in 69.1% with the RT-MRCP and in 57.4% with the CS-BH-MRCP.
ERCP was rarely performed in our study and was often delayed. Furthermore, the sometimes lengthy times between MRCP and ERCP might explain the discrepancy between their findings, since spontaneous migration of the stone to the digestive tract might occur before ERCP. In addition, for patients who underwent ERCP in other centers, the results were not always available. We were therefore unable to assess and compare the diagnostic performance of the MRCP sequences. However, our analysis did not show a significant difference for the detection of bile duct lithiasis between the two sequences. Of the 11 patients with detected lithiasis, 6 had stones visible on the CS-BH-MRCP sequence. Given the short acquisition time of this sequence, it would be of interest to determine the sensitivity of a second acquisition in the event of bad image quality on the first acquisition. Tokoro et al.’s study suggested that the addition of the CS-BH-MRCP to the conventional MRCP protocol at 3T added value to the MRCP examination, since the CS-BH-MRCP could compensate for the image deterioration of the RT-MRCP caused by motion artefacts, although the image quality of the CS-BH-MRCP was not better than the RT-MRCP [
34].
We were unable to further analyze the data on pancreatic cystic lesions, due to the weak sample size of this subgroup and the absence of available ERCP results. Nonetheless, the BH sequence visualized the proximal and distal parts of the MPD more clearly compared to the free-breathing sequence. In addition, a communication between the pancreatic cyst and the MDP was slightly more often visualized with the CS-BH-sequence, allowing the diagnosis of BD-IPMN. Therefore, the BH sequence may be relevant for evaluating pancreatic duct disorders, provided the acquisition volume is well centered on the MPD. The optimized CS-BH-MRCP proposed by Song et al. showed very interesting results and significantly better demonstrated the communication between the pancreatic cyst and the MPD as compared to the conventional MRCP [
32]. With 41 patients included for the evaluation BD-IPMN using MRCP at 1.5T, the short single BH CS-MRCP sequence at 1.5T proposed by Henninger et al. demonstrated significantly higher scores in all the diagnostic approach criteria (lesion conspicuity, confidence, communication) compared to the conventional NT-MRCP, a CS-NT-MRCP, and a long single BH CS-MRCP sequences [
33]. CS-BH-MRCP sequences that are specifically optimized for pancreatic ducts diseases assessment, could therefore improve the diagnostic performance in this indication.
The strengths of this study included the prospective design, reading of the images by two observers, exhaustive analysis of image parameters, and large number of patients compared to the published reports. However, our study presented several limitations. First, it was a single-center study. Second, the image quality analysis was mainly subjective. Nevertheless, the analysis was performed by two readers and interobserver agreement for the assessment of image quality qualitative evaluation variables ranged from moderate to substantial. Third, the excessively small subgroup sizes of patients with pancreatic cyst and liver-transplant recipients did not provide us enough data to allow a meaningful statistical analysis. Fourth, ERCP was rarely performed to confirm the diagnosis or the results were not available. We were therefore unable to assess the diagnostic performance of the MRCP sequences. However, this study mainly focused on image quality assessment. Fifth, blinded reading of the image parameters was biased by the recognizable appearance of CS-BH-MRCP images. Sixth, many acquisition parameters differed between the two sequences, that may make comparison challenging. However, this study compared a RT sequence with a BH-CS sequence. These are fundamentally different scan procedures that cannot be performed with identical protocol parameters. As such, these differences are not a true limitation of the study design but an inherent consequence of the applied techniques.