- Open Access
Feasibility study of dual parametric 2D histogram analysis of breast lesions with dynamic contrast-enhanced and diffusion-weighted MRI
Journal of Translational Medicinevolume 16, Article number: 325 (2018)
This study aimed to investigate the diagnostic value of a dual-parametric 2D histogram classification method for breast lesions.
This study included 116 patients with 72 malignant and 44 benign breast lesions who underwent CAIPIRINHA-Dixon-TWIST-VIBE dynamic contrast-enhanced (CDT-VIBE DCE) and readout-segmented diffusion-weighted magnetic resonance examination. The volume of interest (VOI), which encompassed the entire lesion, was segmented from the last phase of DCE images. For each VOI, a 1D histogram analysis (mean, median, 10th percentile, 90th percentile, kurtosis and skewness) was performed on apparent diffusion coefficient (ADC) and volume transfer constant (Ktrans) maps; a 2D histogram image (Ktrans-ADC) was generated from the pixelwise aligned maps, and its kurtosis and skewness were calculated. Each parameter was correlated with pathological results using the Mann–Whitney test and receiver operating characteristic curve analysis.
For the Ktrans histogram, the area under the curve (AUC) of the mean, median, 90th percentile and kurtosis had statistically diagnostic values (mean: 0.760; median: 0.661; 90th percentile: 0.781; and kurtosis: 0.620). For the ADC histogram, the AUC of the mean, median, 10th percentile, skewness and kurtosis had statistically diagnostic values (mean: 0.661; median: 0.677; 10th percentile: 0.656; skewness: 0.664; and kurtosis: 0.620). For the 2D Ktrans-ADC histogram, the skewness and kurtosis had statistically higher diagnostic values (skewness: 0.831, kurtosis: 0.828) than those of the 1D histogram (all P < 0.05).
The dual-parametric 2D histogram analysis revealed better diagnostic accuracy for breast lesions than single parametric histogram analysis of either Ktrans or ADC maps.
Dynamic contrast-enhanced (DCE) MRI detects sequential changes in contrast enhancement that can be interpreted for quantitative analysis [1, 2]. Pharmacokinetic parameters based on the Tofts model  derived from DCE MRI have shown to have a good correlation with tumor angiogenesis [4, 5], and among these parameters, the volume transfer constant (Ktrans) has potential for differentiating of breast lesions. Malignant tumors tend to have higher Ktrans values than benign tumors [4,5,6,7]. However, the reported Ktrans values for breast lesions and the diagnostic accuracy of these values fluctuates with different MRI acquisition environment and postprocessing techniques. One particular issue with DCE MRI of previous studies is a lack of simultaneously high temporal and high spatial resolution DCE sequences due to technological limitations [4, 7,8,9]. According to the model, higher temporal resolution helps to modify the dynamic curve to be closer to the real physiological situation. Though not directly comparable, studies have verified that high temporal resolution sequences have better diagnostic performances than low temporal resolution sequences [5, 10]. Recently, high temporal resolution can be achieved while spatial resolution is preserved with controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA)-Dixon-time-resolved angiography with stochastic trajectories (TWIST)-volumetric interpolated breath-hold examination (VIBE) DCE and this sequence has demonstrated promising results [11, 12].
Diffusion-weighted imaging (DWI) is another widely studied MRI technique based on the assumption that the movement of water molecules is more restricted in malignant tumors [13,14,15,16]. The apparent diffusion coefficient (ADC), which is derived from DWI, has shown to be useful for differentiating breast lesions. The reported mean ADC values and their diagnostic value varies depending on different DWI scanning environments [13,14,15,16,17,18] and image interpretations. Nevertheless, there is substantial overlap of ADC values between benign and malignant lesions. With the introduction of readout-segmented echo-planar imaging (RS-EPI) for breast MR imaging, DWI image quality has significantly improved and lesion conspicuity of RS-EPI DWI is better than that of single-shot echo-planar imaging (SS-EPI) with 3 T MRI .
One possible approach is to combine DCE and DWI to obtain better performance in terms of differential diagnosis. Histogram analysis is a mathematical method [12, 20,21,22] that excels in describing tumor heterogeneity by providing quantitative metrics. With the two state of the art acquisition sequences, CAIPIRINHA-Dixon-TWIST (CDT)-VIBE DCE and readout-segmented diffusion-weighted (RS-DWI) MRI, we proposed a dual parametric 2D histogram to combine Ktrans and ADC in this study and hypothesized that this histogram would perform better than single parametric histogram for differentiating breast lesions.
The institutional review board approved this prospective study and waived the informed consent requirement. From April 2014 to August 2017, 130 female patients with breast masses palpated or indicated by mammography or ultrasound were scanned with CDT-VIBE DCE and RS-DWI MRI. Fourteen patients were excluded for the following reasons: (a) no obvious lesion was detected on MRI (n = 8), (b) pathology revealed a borderline phyllode tumor (n = 3) or lymphoma (n = 1), and (c) there was incomplete fat suppression on DCE images (n = 2). The enrolled patients underwent resection surgery after MRI examinations. For patients with multiple lesions, the largest one was evaluated. Therefore, the final study cohort included 116 lesions from 116 patients. Morphological features and enhancing patterns were described according to the Breast Imaging Reporting and Data System (BI-RADS).
MR imaging protocol
All patients underwent breast MRI in the prone position with a 3 T MR scanner (Skyra; Siemens Healthcare, Erlangen, Germany) using a dedicated 16-channel bilateral breast coil. Dynamic imaging was performed first with the prototype CDT-VIBE DCE sequence, which consisted of two sequences: 39 s for T1 mapping (using two flip angles: 2 and 14) and 4.44 min for 40 phases of DCE imaging. The two sequences shared the same spatial resolution and field of view (FOV). Gadolinium–diethylene triamine pentaacetic acid (Gd-DTPA) (Magnevist, Bayer Healthcare, Berlin, Germany) was administered intravenously with a power injector at the beginning of the 4th phase at a dose of 0.1 mmol/kg of body weight and a rate of 2 mL/sr. Axial (5 min after contrast injection) and sagittal (7 min after contrast injection) contrast-enhanced images were also obtained for routine clinical diagnosis. Then, the RS-DWI sequence was performed (17 min after contrast injection). Detailed parameters of the sequences are listed in Table 1.
Pharmacokinetic analysis was performed using the DCE software package (Tissue 4D, version: syngo MR D13, Siemens Healthcare, Germany). After motion correction of the 40 phases of DCE images, pharmacokinetic parameters were calculated based on the two-compartment Tofts model  with a population average arterial input function (slow type) provided by Tissue 4D.
ADC maps were generated after DWI scans and exported to a prototype registration package based on the Insight Toolkit (ITK). Each ADC map was rigidly registered and resliced to be aligned with the last phase of the DCE images. The volume of interest (VOI) was segmented using the active contour segmentation mode of the ITK-SNAP software (version 3.2, University of Pennsylvania, Philadelphia, USA)  to cover the entire breast lesion on each slice. Cystic areas, vessels, calcifications and artifacts were removed. All the VOIs were drawn by one junior radiologist and confirmed by another senior radiologist.
For each ADC and Ktrans 1D histogram, the following parameters were calculated: mean, median, 10th percentile, 90th percentile, kurtosis and skewness. A normalized 2D histogram (Ktrans-ADC) map was calculated from the registered and resliced Ktrans and ADC maps for each VOI in a prototype dual-parameter mapping package, and its kurtosis and skewness were calculated. Each parameter was correlated with pathological results.
Histogram parameters and other continuous data are expressed as the mean (± standard deviation) and were compared with the Mann–Whitney test. Categorical data was compared with the Chi-squared test. Receiver operating characteristic (ROC) curve analysis was used to evaluate the effectiveness of histogram parameters for differentiating breast lesions. The optimal thresholds were chosen at the point closest to the top-left part of the plot with perfect sensitivity or specificity. Two ROC curves were compared using the DeLong method. Statistical analyses were performed with MedCalc (Version 16.8, MedCalc Software, Mariakerke, Belgium). P values of less than 0.05 were deemed to indicate statistical significance.
Pathological results revealed that 72 lesions were malignant tumors, comprising 62 invasive ductal carcinomas (IDCs), 7 ductal carcinomas in situ (DCIS), 1 invasive lobular carcinoma (ILC), 1 solid papillary carcinoma (SPC), and 1 myoepithelial carcinoma (MEC), while 44 lesions were benign breast diseases, comprising 21 fibroadenomas, 10 intraductal papillomas (IPs), 6 cases of benign adenosis, 5 benign phyllodes tumors and 2 cases of mastitis. The lesions varied in long diameter from 9 to 74 mm and there was no significant difference in lesion size between the two groups. The detailed morphological features and enhancing patterns are listed in Table 2.
Histogram analysis and comparisons between benign and malignant lesions are shown in Table 3. Malignant tumors had higher Ktrans values and lower ADC values than benign tumors. For Ktrans, the mean (P < 0.0001), median (P = 0.003), 90th percentile (P < 0.0001) were significantly different between benign and malignant tumors, and for ADC, the mean (P = 0.013), 10th percentile (P = 0.024) and median (P = 0.006) were significantly different between benign and malignant tumors. When comparing distribution parameters, the kurtosis of Ktrans (P = 0.028) and the skewness (P = 0.011) of ADC were significantly different between benign and malignant tumors. In the 2D histogram analysis, both skewness (P < 0.0001) and kurtosis (P < 0.0001) were significantly different between benign and malignant tumors. Figure 1 shows the two situations in which Ktrans and ADC resolved each other’s ambiguity. Figure 2 shows the intuitive differences in skewness and kurtosis in Ktrans and ADC histograms between malignant tumor and benign lesions.
ROC curve analysis results
The ROC curve analysis results are shown in Table 4 and Fig. 3. For Ktrans histograms, the 10th percentile revealed the highest AUC (0.781), while the median of the ADC histogram revealed the highest AUC (0.677). The AUC of the 2D Ktrans-ADC histogram-derived skewness and 2D kurtosis were were 0.831 and 0.828, respectively, and were significantly higher than all of the Ktrans and ADC histogram parameters (all P < 0.05), reaching a sensitivity of 68.18% and 75% for skewness and kurtosis, respectively, and a specificity of 84.72% and 75.0% for skewness and kurtosis, respectively.
In this study, we combined the Ktrans derived from DCE and the ADC from DWI to form a 2D histogram and tested its feasibility for differentiating breast lesions. The results showed that the skewness and kurtosis of the 2D Ktrans-ADC histogram were significantly different between malignant and benign breast lesions and had better diagnostic performance than 1D histogram parameters.
Ktrans refers to the coefficient of transendothelial transport of contrast media from the vascular compartment to the tissue interstitium and is correlated to blood flow and permeability. In a malignant tumor microenvironment, neoangiogenesis is more active and intense, and myriads of tiny leaky vessels develop, resulting in increased blood flow and permeability . Thus, the Ktrans values of malignant tumors are higher than those of benign tumors. In this study, with a fairly high temporal resolution of 5.6 s, the 90th percentile of Ktrans had the highest AUC, which was also statistically higher than the 10th percentile, probably because voxels with higher Ktrans values represent the most aggressive regions of tumors with abundant blood perfusion. This is consistent with the results of a similar Ktrans histogram analysis from another study . Of the Tofts model parameters, Ktrans and the extravascular extracellular space (EES) fractional volume (Ve) are directly related to fundamental physiology and require access to the absolute values of tracer concentration, whereas Kep is the ratio of Ktrans to Ve (Kep = Ktrans/Ve), which can be derived from data regarding the shape of the tracer concentration vs time. Ktrans is more widely studied than Kep for differentiation of lesions, as the clinical value of Kep is discrepant between studies [5, 7, 25, 26]. Considering the above mentioned reasons, we chose to use Ktrans to conduct this study.
The ADC histogram analysis showed that voxels with lower ADC values had better diagnostic value than those with higher values, which is consistent with the results reported by Suo et al.  that the minimal ADC value had the highest AUC. This trend may be reflected by the histological heterogeneity of tumors; for instance, IDC usually consists of a variety of components such as cancer nests, stroma, intratumoral fibrosis or necrosis and intraductal components, and thus presents radiological heterogeneity on ADC maps. The 10th percentile of the ADC histogram may indicate the more invasive part of malignant breast tumors. The ADC values of breast lesions in this study were relatively lower than those reported in other study that used RS-DWI . One possible reason is that we used a whole lesion VOI, which involved components such as stroma and intratumoral fibrosis that have higher ADC values than tumor nests. Another issue is that the ADC map was linearly resliced according to the Ktrans map to generate the pixelwise Ktrans-ADC map, which might have resulted in a slight decrease in ADC values.
Multiparametric imaging analysis has raised a lot of interests in recent years. In early pilot studies, the ADC was taken as a supplementary parameter to the DCE imaging to obtain a more definitive diagnosis [17, 27, 28]. However, in these studies, data were not quantified, and diagnostic accuracy was still restricted by the experience of radiologists. To develop quantitative and objective methods for interpreting data, several statistical models, such as the Local Hyperplane-based relief feature weighting scheme , the linear discriminant analysis  and the step-wise multivariate logistic regression model,  were used to incorporate DCE and DWI information.
Histograms depict the distribution of datasets; kurtosis measures how extreme observations are, and skewness measures the asymmetry of distribution. From our results, the kurtosis of the 1D and 2D histograms was higher (although the kurtosis of the ADC histogram revealed no significant differences) for benign lesions than for malignant lesions. This is probably because cells are more uniform and regularly aligned with balanced blood flow in benign lesions, where most of the voxels tend to crowd together, and form a sharper peak. Positive skewness is called right-skewed, and indicates that the tail on the right side is longer than that on the left side, and the mass of the distribution is concentrated on the left, making the histogram curve lean to the left. In our results, the ADC skewness of malignant tumors was positive, and that of benign lesions was negative. Thus, malignant tumor histogram curves lean to the left while benign lesion histogram curves lean to the right. This is because most of the voxels in malignant tumors have low ADC values and lie to the left side of the histogram, while voxels in benign lesion tend to lie to the right side. Similar results were also obtained in another ADC histogram analysis study . If the skewness of two parameters has the same sign, the parameter with the higher absolute value is more skewed. As the results revealed, the Ktrans skewness of benign lesions was higher than that of malignant tumors, probably because the Ktrans values of benign lesions were lower than those of malignant lesions, and most voxels lie to the left, making the curve lean more to the left. Such differences might not be statistically significant in 1D histograms, as the kurtosis of ADC and the skewness of Ktrans in this study and other studies [12, 20, 32], but these differences are amplified in the 2D histogram.
This study had several limitations. First, there was misregistration between ADC map and DCE images, which could have influenced VOIs, especially for small lesions and non-mass lesions. Another limitation is that DWI was performed approximately 17 min after the administration of contrast media. Although a prior study has shown that there are no significant effects on ADC values after the administration of a gadolinium-based contrast agent , it may be preferable to acquire the DWI sequence before contrast injection to avoid any confounding effects. Finally, the sample size of benign lesions was relatively small. Thus, this method needs to be tested on a larger scale.
In summary, the combination of Ktrans and ADC performed better than single parameters alone using histogram analysis for differentiating breast lesions. The 2D kurtosis and 2D skewness are both candidates to help characterize breast lesions.
magnetic resonance imaging
controlled aliasing in parallel imaging results in higher acceleration
time-resolved angiography with stochastic trajectories
volumetric interpolated breath-hold examination
readout-segmented diffusion-weighted imaging
volume of interest
apparent diffusion coefficient
volume transfer constant
area under the curve
Breast Imaging Reporting and Data System
receiver operating characteristic
invasive ductal carcinomas
ductal carcinomas in situ
invasive lobular carcinoma
solid papillary carcinoma
Kuhl CK, Schild HH, Morakkabati N. Dynamic bilateral contrast-enhanced MR imaging of the breast: trade-off between spatial and temporal resolution. Radiology. 2005;236(3):789–800.
Flickinger FW, Allison JD, Sherry RM, Wright JC. Differentiation of benign from malignant breast masses by time-intensity evaluation of contrast enhanced MRI. Magn Reson Imaging. 1993;11(5):617–20.
Tofts PS, Brix G, Buckley DL, Evelhoch JL, Henderson E, Knopp MV, et al. Estimating kinetic parameters from dynamic contrast-enhanced T(1)-weighted MRI of a diffusable tracer: standardized quantities and symbols. J Magn Reson Imaging. 1999;10(3):223–32.
Furman-Haran E, Schechtman E, Kelcz F, Kirshenbaum K, Degani H. Magnetic resonance imaging reveals functional diversity of the vasculature in benign and malignant breast lesions. Cancer. 2005;104(4):708–18.
El Khouli RH, Macura KJ, Kamel IR, Jacobs MA, Bluemke DA. 3-T dynamic contrast-enhanced MRI of the breast: pharmacokinetic parameters versus conventional kinetic curve analysis. AJR Am J Roentgenol. 2011;197(6):1498–505.
Hauth EAM, Jaeger HJ, Maderwald S, Muehler A, Kimmig R, Forsting M. Quantitative 2- and 3-dimensional analysis of pharmacokinetic model-derived variables for breast lesions in dynamic, contrast-enhanced MR mammography. Eur J Radiol. 2008;66(2):300–8.
Veltman J, Stoutjesdijk M, Mann R, Huisman HJ, Barentsz JO, Blickman JG, et al. Contrast-enhanced magnetic resonance imaging of the breast: the value of pharmacokinetic parameters derived from fast dynamic imaging during initial enhancement in classifying lesions. Eur Radiol. 2008;18(6):1123–33.
Li L, Wang K, Sun X, Wang K, Sun Y, Zhang G, et al. Parameters of dynamic contrast-enhanced MRI as imaging markers for angiogenesis and proliferation in human breast cancer. Med Sci Monit. 2015;21:376–82.
Yim H, Kang DK, Jung YS, Jeon GS, Kim TH. Analysis of kinetic curve and model-based perfusion parameters on dynamic contrast enhanced MRI in breast cancer patients: correlations with dominant stroma type. Magn Reson Imaging. 2016;34(1):60–5.
Amarnath J, Sangeeta T, Mehta SB. Role of quantitative pharmacokinetic parameter (transfer constant: K(trans)) in the characterization of breast lesions on MRI. Indian J Radiol Imaging. 2013;23(1):19–25.
Michaely HJ, Morelli JN, Budjan J, Riffel P, Nickel D, Kroeker R, et al. CAIPIRINHA-Dixon-TWIST (CDT)-volume-interpolated breath-hold examination (VIBE): a new technique for fast time-resolved dynamic 3-dimensional imaging of the abdomen with high spatial resolution. Invest Radiol. 2013;48(8):590–7.
Li Z, Ai T, Hu Y, Yan X, Nickel MD, Xu X, et al. Application of whole-lesion histogram analysis of pharmacokinetic parameters in dynamic contrast-enhanced MRI of breast lesions with the CAIPIRINHA-Dixon-TWIST-VIBE technique. J Magn Reson Imaging. 2018;47(1):91–6.
Marini C, Iacconi C, Giannelli M, Cilotti A, Moretti M, Bartolozzi C. Quantitative diffusion-weighted MR imaging in the differential diagnosis of breast lesion. Eur Radiol. 2007;17(10):2646–55.
Rubesova E, Grell A-S, De Maertelaer V, Metens T, Chao S-L, Lemort M. Quantitative diffusion imaging in breast cancer: a clinical prospective study. J Magn Reson Imaging. 2006;24(2):319–24.
Guo Y, Cai Y-Q, Cai Z-L, Gao Y-G, An N-Y, Ma L, et al. Differentiation of clinically benign and malignant breast lesions using diffusion-weighted imaging. J Magn Reson Imaging. 2002;16(2):172–8.
Sinha S, Lucas-Quesada FA, Sinha U, DeBruhl N, Bassett LW. In vivo diffusion-weighted MRI of the breast: potential for lesion characterization. J Magn Reson Imaging. 2002;15(6):693–704.
Partridge SC, DeMartini WB, Kurland BF, Eby PR, White SW, Lehman CD. Quantitative diffusion-weighted imaging as an adjunct to conventional breast MRI for improved positive predictive value. AJR Am J Roentgenol. 2009;193(6):1716–22.
Bogner W, Gruber S, Pinker K, Grabner G, Stadlbauer A, Weber M, et al. Diffusion-weighted MR for differentiation of breast lesions at 3.0 T: how does selection of diffusion protocols affect diagnosis? Radiology. 2009;253(2):341–51.
Bogner W, Weber M, Helbich TH, Trattnig S. Readout-segmented echo-planar imaging improves the diagnostic performance of diffusion-weighted MR breast examinations at 30 T. Radiology. 2012;263(1):64–76.
Suo S, Zhang K, Cao M, Suo X, Hua J, Geng X, et al. Characterization of breast masses as benign or malignant at 3.0T MRI with whole-lesion histogram analysis of the apparent diffusion coefficient. J Magn Reson Imaging. 2016;43(4):894–902.
Zhang Y-D, Wu C-J, Wang Q, Zhang J, Wang X-N, Liu X-S, et al. Comparison of utility of histogram apparent diffusion coefficient and R2* for differentiation of low-grade from high-grade clear cell renal cell carcinoma. AJR Am J Roentgenol. 2015;205(2):W193–201.
Kim EJ, Kim SH, Park GE, Kang BJ, Song BJ, Kim YJ, et al. Histogram analysis of apparent diffusion coefficient at 3.0 T: correlation with prognostic factors and subtypes of invasive ductal carcinoma. J Magn Reson Imaging. 2015.
Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. Neuroimage. 2006;31(3):1116–28.
Kuhl CK, Mielcareck P, Klaschik S, Leutner C, Wardelmann E, Gieseke J, et al. Dynamic breast MR imaging: are signal intensity time course data useful for differential diagnosis of enhancing lesions? Radiology. 1999;211(1):101–10.
Schabel MC, Morrell GR, Oh KY, Walczak CA, Barlow RB, Neumayer LA. Pharmacokinetic mapping for lesion classification in dynamic breast MRI. J Magn Reson Imaging. 2010;31(6):1371–8.
Jena A, Mehta SB, Taneja S. Optimizing MRI scan time in the computation of pharmacokinetic parameters (K(trans)) in breast cancer diagnosis. J Magn Reson Imaging. 2013;38(3):573–9.
Yabuuchi H, Matsuo Y, Okafuji T, Kamitani T, Soeda H, Setoguchi T, et al. Enhanced mass on contrast-enhanced breast MR imaging: lesion characterization using combination of dynamic contrast-enhanced and diffusion-weighted MR images. J Magn Reson Imaging. 2008;28(5):1157–65.
Kul S, Cansu A, Alhan E, Dinc H, Gunes G, Reis A. Contribution of diffusion-weighted imaging to dynamic contrast-enhanced MRI in the characterization of breast tumors. AJR Am J Roentgenol. 2011;196(1):210–7.
Cai H, Liu L, Peng Y, Wu Y, Li L. Diagnostic assessment by dynamic contrast-enhanced and diffusion-weighted magnetic resonance in differentiation of breast lesions under different imaging protocols. BMC Cancer. 2014;14:366.
Fusco R, Sansone M, Filice S, Granata V, Catalano O, Amato DM, et al. Integration of DCE-MRI and DW-MRI quantitative parameters for breast lesion classification. Biomed Res Int. 2015;2015:237863.
Rahbar H, Partridge SC, Demartini WB, Gutierrez RL, Allison KH, Peacock S, et al. In vivo assessment of ductal carcinoma in situ grade: a model incorporating dynamic contrast-enhanced and diffusion-weighted breast MR imaging parameters. Radiology. 2012;263(2):374–82.
Shin HJ, Kim SH, Lee HJ, Gong G, Baek S, Chae EY, et al. Tumor apparent diffusion coefficient as an imaging biomarker to predict tumor aggressiveness in patients with estrogen-receptor-positive breast cancer. NMR Biomed. 2016;29(8):1070–8.
Arponen O, Sudah M, Sutela A, Taina M, Masarwah A, Liimatainen T, et al. Gadoterate meglumine decreases ADC values of breast lesions depending on the b value combination. Sci Rep. 2018;8(1):87.
BW, YG and HL formulated the research goals and aims. HL conceived the prototypes for data analysis. DW did the MRI scan and preserved the data. YC conceived data collection, interpreted the data, and drafted and revised the manuscript. All authors read and approved the final manuscript.
We thank all staffs involved in the acquisition of data. We are grateful to all involved cancer patients for their participation to the study. Also, we thank American Journal Experts for their support in linguistic edit of the manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated during the current study are not publicly available due to the patient confidentiality but are available from the corresponding author on reasonable request.
Consent for publication
There were MRI images presented in Fig. 1 of this article from patients. Consent for publication was obtained.
Ethics approval and consent to participate
The institutional review board approved this prospective study and waived the informed consent requirement.
This study was funded by Shanghai Shenkang Hospital Development Center, Project Number: SKJS1502; and China National Center for Biotechnology Development, Project Number: 2016YFC1303004.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.