Abstract
The aim of this study was to evaluate the dynamic variations in the quantitative parameters of diffusion-weighted imaging (DWI) at different b-value combinations in a prostate cancer (PCa) mouse model for noninvasive monitoring of histopathological changes. Twenty-five male C57BL/6J mice were randomly allocated into a control group (n=5) or an experimental group (n=20). The experimental groups were used to establish the PCa model. On days 9, 12, 15, and 18 post-modeling, 5 mice were randomly selected for MRI, including T1WI, T2WI, T2WI SPIR, and DWI. The b-values were set at 0, 500, 1000, 1500, and 2000 s/mm2. Apparent diffusion coefficient (ADC) and exponential apparent diffusion coefficient (EADC) values from different b-value combinations were measured. Post-MRI, tumors were excised for histopathological analysis. DWI quantitative parameters, tumor nuclear fraction, and Ki-67 area fraction were compared on different days, along with correlation analysis. ADC values gradually decreased as tumor progressed, whereas EADC values gradually increased. Tumor nuclear fraction increased over time. Ki-67 increased first and then decreased. Tumor nuclear fraction was negatively correlated with the ADC value and positively correlated with the EADC value. The Ki-67 was positively correlated with the ADC value and negatively correlated with the EADC value. ADC values at b=1000, 1500 s/mm2 and the EADC values at b=0, 500 s/mm2 demonstrated the strongest correlations with the tumor nuclear fraction; the ADC and EADC values at b=500, 1000 s/mm2 were more strongly correlated with Ki-67, being potential noninvasive imaging biomarkers for monitoring changes in tumor histopathology.
Prostate cancer; Diffusion-weighted imaging; Apparent diffusion coefficient; Exponential apparent diffusion coefficient; Tumor nuclear fraction
Introduction
Prostate cancer (PCa) represents one of the most prevalent malignant tumors in males globally, ranking second to lung cancer in terms of mortality (1,2). PCa is characterized by rapid progression and marked heterogeneity (3,4). Accurately estimating tumor proliferation conditions is crucial for grading tissue, tailoring treatments to individual patients, and evaluating treatment efficacy (5,6). However, the procedure for obtaining human pathological biopsies is invasive (7,8), complicating the dynamic observation of histopathological alterations in tumors. Consequently, there is a critical demand for noninvasive methodologies capable of monitoring histopathological alterations, facilitating continuous observation and supporting clinical interventions.
Magnetic resonance imaging (MRI) is a noninvasive approach with superior soft tissue resolution. Multiparametric magnetic resonance imaging (mp-MRI) is extensively applied to diagnose, stage, and monitor prostate cancer (9- 11). Diffusion-weighted imaging (DWI), a functional imaging sequence, allows for the routine calculation of apparent diffusion coefficient (ADC) values to assess the diffusion of water molecules in tissue (12,13). These quantitative parameters will effectively evaluate histopathological and biological changes (14- 16). However, previous studies used only a fixed combination of b-values in DWI for quantitative parameter determination (17,18) and did not explore the significance of the optimal diffusion sensitivity factor recommended by the latest Prostate Imaging Reporting and Data System (PI-RADS v2.1) or the combination of multi b-values (19). In addition, no study has investigated the relationship between the exponential apparent diffusion coefficient (EADC) and the pathologic progression of prostate cancer.
In this study, we established a C57BL/6J mouse model of prostate cancer and conducted multi b-value DWI at various stages of tumor progression. The objective of this study was to identify which combination of b-values is more effective for noninvasive and dynamic monitoring of tumor growth and pathological changes, thereby providing valuable insights for preclinical research, active surveillance, and assessment of therapeutic interventions in prostate cancer.
Material and Methods
Experimental animals
The Animal Ethical Care Committee of Qiqihar Medical University approved the animal experimental protocol. In this study, we selected 25 C57BL/6J male mice, aged 6-8 weeks and weighing 18-22 g, sourced from Liaoning Changsheng Biotechnology Co., Ltd. [License No. SCXK (Liao) 2020-0001; China]. The animal facility maintained regular ventilation with constant temperature and humidity (temperature: 20-24°C, humidity: 45-55%). The mice had ad libitum access to food and water. The mice were randomly allocated to the control group (n=5) or the experimental groups (n=20) via a random number table. Within the experimental group, five mice were selected randomly for further experiments on the 9th, 12th, 15th, and 18th days post-modeling.
Preparation of the orthotopic PCa animal model
The mouse prostate cancer cell line RM-1 was cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum and 1% penicillin-streptomycin solution (100 U/mL penicillin and 0.1 mg/mL streptomycin) and incubated in a 5% CO2 atmosphere at 37°C. The cells in the logarithmic growth phase were harvested, digested, and resuspended at a concentration of 2×107 cells/mL. The mice were anesthetized with 4% pentobarbital sodium and immobilized in a supine position on cardboard. The abdominal skin was disinfected before surgery. Under a stereomicroscope, a transverse incision approximately 1.0 cm in length was made in the midline position at the base of the thigh, adjacent to both sides of the lower abdomen, to access the abdominal cavity. The bladder was located and repositioned cranially via tweezers, displacing surrounding adipose tissue to reveal the prostate. A microinjector introduced 10 μL of the RM-1 cell suspension into the prostate capsule, with a visible bulge indicating successful injection. The organs were repositioned, and the abdominal muscle and skin were sutured separately via 6-0 monofilament sutures. Post-operation, the wounds were disinfected with iodophor. The mice in the control group received standard care.
MRI techniques
All the scans were performed via a 3.0 T MR scanner (Achieva; Philips Healthcare, The Netherlands) with an 8-channel knee coil. Within the experimental group, five mice were randomly chosen for data collection and underwent conventional MRI and multi b-value DWI sequences on the 9th, 12th, 15th, and 18th days post-modeling. Five control mice were subjected to repeated scanning. Post-anesthesia, mice were positioned prone on a custom water model for scanning. The sequence and parameter details were as follows: coronal fast spin echo (FSE)-T1-weighted imaging (T1WI): repetition time (TR) 633 ms, echo time (TE) 10 ms, field of view (FOV) 100×100 mm, matrix 192×155, slice thickness 2 mm, slice gap 0.2 mm, number of excitations (NEX) 2; coronal FSE-T2-weighted imaging (T2WI): TR 3000 ms, TE 80 ms, FOV 100×100 mm, matrix 192×155, slice thickness 2 mm, slice gap 0.2 mm, NEX 2; axial T1WI, T2WI, T2WI spectral pre-saturation with inversion recovery (SPIR): matrix 212×161, with other parameters consistent with the aforementioned settings; and coronal echo planar imaging (EPI)-DWI: TR 2000 ms, TE 67 ms, FOV 120×120 mm, matrix 64×63, slice thickness 2 mm, slice gap 0.2 mm, NEX 5, and diffusion sensitivity factor (b-values) set at 0, 500, 1000, 1500, and 2000 s/mm2.
MRI post-processing and quantitative measurements
The Philips IntelliSpace Portal post-processing workstation facilitated image processing, generating ADC and EADC maps via MR Diffusion software (Achieva; Philips Healthcare). The delineation of the region of interest (ROI) and measurement of pertinent parameters were collaboratively conducted by two experienced radiologists, with consultations from senior radiologists sought in instances of differing opinions. Both were blinded to the grouping of the mice and the histopathological results. The methodology for calculating the volume of MRI-measured tumors involved measuring the maximum anterior-posterior diameter (a), maximum transverse diameter (b), and maximum longitudinal diameter (c) of the prostate tumor via coronal and axial T2WI. The tumor volume was computed via the following formula: V = a × b × c × 0.52 mm3. T2WI served as the anatomical reference for ROI delineation. On DWI with a b-value of 0 s/mm2, the largest layer of the tumor was chosen, and three circular ROIs of 3 mm2 each were positioned to circumvent areas of hemorrhage and necrosis. ADC and EADC reconstructions were conducted at b-values combinations of 500, 1000, 1500, 2000 s/mm2 and 0 s/mm2, respectively, with the average ADC and EADC values of the ROIs automatically calculated. By modifying the combinations of b-values (b=500, 1000/1500/2000 s/mm2, b=1000, 1500/2000 s/mm2, and b=1500, 2000 s/mm2), the ADC and EADC values of different combinations of b-values were obtained. The quantitative parameters were derived from the monoexponential equation (20):
where S denotes the signal intensity at different b-values, b denotes the b-value, and EADC = S2 / S1. The three ROIs delineated on DWI were transferred to the coronal T2WI, and three ROIs of equivalent size were delineated on the right thigh muscle within the same plane. The signal intensity ratio between them was designated by the signal intensity ratio (SIR) of prostate cancer, which signifies the T2WI signal of the tumor. The measurement formula is as follows: SIR = Stumor / Smuscle.
Histopathological analysis
Following MRI examination, tumor tissues from the mice in the experimental group were excised after cardiac perfusion with 4% paraformaldehyde. Normal prostate tissue from the mice in the control group was excised in the same way after an MRI examination on the 18th day. These tissues were subsequently rinsed with saline. Three-dimensional diameters were measured via calipers to determine the tumor volume. The tissues were then fixed in 4% paraformaldehyde for 48 h, dehydrated, and embedded in paraffin. Paraffin-embedded sections, approximately 4-µm thick, were subjected to hematoxylin and eosin (HE) staining to facilitate observation of viable tissue, necrosis, and tumor invasion in prostate samples under a microscope (Nikon, Japan). In the experimental group, three sections per mouse sample were chosen at random, and within each section, three fields were selected for analysis at 400× magnification to compute the mean tumor nuclear fraction. The methodology for assessing the tumor nuclear fraction (6,21) involved the use of ImageJ software (v.1.53; National Institutes of Health, USA) to convert HE-stained images to 8-bit grayscale images, setting the threshold from 0 to 110 to isolate the nuclear area. The tumor nuclear fraction was calculated as the ratio of the nuclear area to the total field area, reported as a percentage.
Ki-67 immunohistochemical staining
The paraffin-embedded sections were dewaxed, hydrated, heated at high temperature for 10 min, inactivated for 10 min, and the serum was blocked for 1 h. A 1:400 dilution of rabbit anti-Ki-67 polyclonal antibody (Bioss, China) was incubated at 4°C overnight, the goat anti-rabbit antibody was incubated at room temperature for 1 h, streptavidin was added at room temperature for 30 min, color was developed, restaining was performed, and, after differentiation, the water was washed back to blue, and the tablets were sealed. A phosphoric acid buffer salt solution was used as a negative control, a normal spleen was used as a positive control, and brown staining of the cell nucleus was used as a positive control. Three sections were selected for each sample in the experimental group, and three fields of view were randomly selected for each section under 400× magnification. Positive areas were extracted via color deconvolution via ImageJ software, and the area fraction of the three fields was measured to calculate the average value, which was used to characterize the proliferation ability of the tumor cells.
Statistical analysis
Data analysis was conducted via SPSS software (v.26.0; IBM, USA). All the data are reported as means±SD. The Shapiro-Wilk test was used to assess the normality of the data. Unpaired t-tests were used for two-group comparisons on normally distributed data, and comparisons across multiple groups were performed via one-way analysis of variance (ANOVA). Two-group comparisons were performed with the Mann-Whitney U test on non-normally distributed data, whereas the Kruskal-Wallis test was employed for multiple-group comparisons. Spearman's rank correlation coefficient was used for correlation analysis. A P-value of less than 0.05 was considered statistically significant.
Results
General situations
All the mice survived until the conclusion of their respective experiments. Initially, the experimental group exhibited satisfactory health conditions, characterized by smooth fur and normal eating behaviors, mirroring those observed in the control group. However, beginning on day 12 post-modeling, a noticeable decline in the health status of the experimental mice was observed. This decline was marked by the thinning and coarsening of their fur, a reduced frequency of eating, drinking, and defecation, and a significant decrease in activity levels. Palpation revealed a hard mass in the lower abdomen. The body weights in the experimental group initially increased post-modeling but significantly decreased in the later stages. The weights of the mice in the control group showed a steady upward trend, without a sharp upward or downward trend (Figure 1).
Changes in the body weight of mice in the experimental and control groups after modeling. Data are reported as means and SD.
MRI results
MRI of the mice in the experimental group at different time points revealed tumor formation in the lower abdomen (Figure 2), which was challenging to differentiate from the normal prostate. The prostate tumor displayed a low signal on T1WI, closely resembling the signal of the testes, and a medium to high signal on T2WI and T2WI SPIR, slightly lower than that of the testes, with a well-defined margin and homogeneous signal. In cross-sectional views, the tumor was confined within the pelvic cavity. The tumor size progressively increased, with scattered areas of long T1 and long T2 necrotic signals emerging within the tumor on days 15 and 18. By day 18, the tumor had expanded to fill the lower abdomen and exert pressure on adjacent tissues, occupying the pelvic cavity in cross-sectional views. In contrast, the prostates of the control group mice, which were situated below the bladder, presented medium to low signals on both T1WI and T2WI, with well-defined margins. No significant difference was found in the images on the different days.
Standard MRI sequences of the mice. A-E, Coronal T1WI images of mice from the control and experimental groups. F-J, Coronal T2WI images from each group. K-O, Axial T2WI SPIR from each group. Red arrows indicate testicles, yellow arrows denote normal prostate tissue, and white arrows point to prostate tumors. Prostate tumors in the experimental group were discernible across all MRI sequences, with clear delineation from adjacent tissues and a noted increase in size over the observation period. MRI: magnetic resonance imaging; T1WI: T1-weighted imaging; T2WI: T2-weighted imaging; T2WI SPIR: T2-weighted imaging with spectral pre-saturation with inversion recovery. Scale bar 500 μm.
The tumor exhibited a markedly high signal intensity on DWI (Figure 3), marginally lower than the testes. Over time, the tumor's signal intensity increased, creating a stark contrast with the adjacent tissue. As the b-values increased, the tumor signal intensity progressively decreased, yet it remained substantially greater than the surrounding tissue. In the control group, the mouse prostate displayed a slightly elevated signal on DWI.
Diffusion-weighted imaging (DWI) of mice at various b-values. A-E, Coronal DWI of mice from the control and experimental groups at a b-value of 500 s/mm2. F-J, Coronal DWI at a b-value of 1000 s/mm2 for each group. K-O, Coronal DWI at a b-value of 1500 s/mm2 for each group. P-T, Coronal DWI at a b-value of 2000 s/mm2 for each group. Red arrows point to testicles, yellow arrows indicate normal prostate tissue, and white arrows indicate prostate tumors. DWI revealed that prostate tumors in the experimental group presented high signal intensity, which became more pronounced over time. As the b-values increase, the signal intensity of prostate tumors decreases at the same time points. Conversely, the DWI signals of normal prostate tissue in the control group were significantly lower than the testicles and prostate tumors. Scale bar 500 μm.
MRI volume measurements revealed a progressive increase in tumor size over time (Table 1), with the volumetric differences between groups reaching statistical significance (P<0.05).
Comparison of tumor volume and signal intensity ratio (SIR) of prostate cancer at each time in the experimental group (n=5).
Changes in the ADC and EADC values
The ADC values at b-value combinations (b=0, 500/1000/1500/2000 s/mm2) for the experimental group at various time points were consistently lower than those of the control group (P<0.001, Table 2), indicating a decreasing trend over time (Figure 4).
Comparison of apparent diffusion coefficient (ADC) values from b-value combinations (b=0, 500/1000/1500/2000 s/mm2) in each group (n=5) (×10-3 mm2/s).
Temporal evolution of quantitative DWI parameters in mice at days 9, 12, 15, and 18. A, ADC value variation across different b-value combinations (b=0, 500/1000/1500/2000 s/mm2) in the experimental group. B, EADC value variation across different b-value combinations (b=0, 500/1000/1500/2000 s/mm2) in the experimental group. C, Variations in the ADC values across different b-value combinations [(b=500, 1000/1500/2000 s/mm2), (b=1000, 1500/2000 s/mm2), (b=1500, 2000 s/mm2)] in the experimental group. D, Variation in EADC values across different b-value combinations [(b=500, 1000/1500/2000 s/mm2), (b=1000, 1500/2000 s/mm2), (b=1500, 2000 s/mm2)] in the experimental group. The ADC values of each b-value combination gradually decreased with time, whereas the EADC values of each b-value combination gradually increased. Data are reported as means and SD. Different letters indicate statistically significant differences (P<0.05) (ANOVA and Kruskal-Wallis test). DWI: diffusion-weighted imaging; ADC: apparent diffusion coefficient; EADC: exponential apparent diffusion coefficient.
In the experimental group, the ADC values of the remaining b-value combinations (b=500, 1000/1500/2000 s/mm2; b=1000, 1500/2000 s/mm2; and b=1500, 2000 s/mm2) also exhibited a decreasing trend over time (P<0.001, Supplementary Table S1, Figure 4). Except the ADC values at b=1500 and 2000 s/mm2, there were no significant differences between the experimental group and the control group at the other b-value combinations on the 9th day (P>0.05). On the 12th day, there were no significant differences in the ADC values between the experimental and control groups when b=500, 1000/2000 s/mm2 (P>0.05). Moreover, statistically significant differences emerged between the experimental and control groups at other time points (P<0.001).
No significant differences were observed in the ADC values at different b-value combinations for the control group at various time points (P>0.05, Table 2 and Supplementary Table S1).
The EADC values at b-value combinations (b=0, 500/1000/1500/2000 s/mm2) for the experimental group at various time points tended to increase over time (Figure 4) and were consistently greater than those of the control group (P<0.001, Table 3).
Comparison of exponential apparent diffusion coefficient (EADC) values from b-value combinations (b=0, 500/1000/1500/2000 s/mm2) in each group (n=5).
The EADC values of the other b-value combinations (b=500, 1000/1500/2000 s/mm2; b=1000, 1500/2000 s/mm2; and b=1500, 2000 s/mm2) in the experimental group also exhibited a consistent increasing trend over time (P<0.001, Supplementary Table S2, Figure 4). Specifically, at b=1500 and 2000 s/mm2, there was no significant difference in the EADC values between the control group and the experimental group on the 9th day (P>0.05); however, statistically significant differences were observed at other time points (P<0.01).
There were no significant differences in the EADC values at different b-value combinations for the control group across various time points (P>0.05, Table 3, Supplementary Table S2).
Changes in the tumor T2WI signal
In the experimental group, there was no statistically significant difference in the signal intensity ratio of prostate cancer at any time point (P>0.05, Table 1).
Histopathological analysis and changes in tumor volume
Upon dissection of the mice in the experimental group, the tumors below the bladder were identified as irregular solid masses. These masses were connected posteriorly to the seminal vesicles and encapsulated, featuring a tough texture (Figure 5). Partial hemorrhage was observed within the capsule. In the central section of the tumor, the tissue appeared gray-white and interspersed with areas of hemorrhage and necrosis, and distinguishing the residual prostate from the tumor tissue was difficult. Anatomical measurements of the tumor volume revealed a progressive increase over time (Table 1), with statistically significant differences in volume at each time point (P<0.05).
Comparative anatomical and HE-stained images of prostate tumors in mice. A-E, General anatomical observation of prostate and prostate tumors from the control and experimental groups. The red arrow indicates the seminal vesicle, the green arrow indicates the bladder, the yellow arrow indicates the normal prostate, and the white arrow indicates the prostate tumor. F-J, HE staining of prostate and prostate tumors from the control and experimental groups (×200, scale bar 100 µm). K-O, HE staining of prostate and prostate tumors from the control and experimental groups (×400, scale bar 50 µm). *Residual normal prostate gland in the tumor; the general anatomy shows that the tumor volume gradually increased with time. HE staining revealed that the tumor cells continued to invade the normal prostate gland, and hemorrhage, necrosis, and structural disorders occurred in the center of the tumor on day 18. HE: hematoxylin and eosin.
HE staining revealed that the tumor cells were densely packed, exhibited irregular shapes, with visible nucleoli and pronounced nuclear staining and pleomorphism. This cellular arrangement was associated with neutrophil infiltration and abnormal mitotic figures (Figure 5). On the 9th day, HE staining revealed the presence of residual normal prostate glands within the tumor cells, indicating the destruction of basal cells and invasion into the surrounding normal prostate tissue. By the 12th day, a significant reduction in the number of normal prostate glands within the tumor area was observed, with tumor cells continuing to invade surrounding glands and expand their invasion area. On the 15th day, the glandular structure within the tumor area was completely obliterated, the central areas of the tumor displayed scattered hemorrhages, and vacuolation began at the tumor periphery. By the 18th day, lamellar hemorrhage and necrosis were evident in the central area of the tumor, disrupting the internal structure of the tumor and leading to vacuolation and hyalinization.
Changes in the tumor nuclear fraction
In the experimental group, the tumor nuclear fraction increased in a time-dependent manner (Figure 6, Table 4). This increase was not statistically significant between the 15th and 18th days (P>0.05). However, significant differences were observed at all other measured time points (P<0.001).
Ki-67 immunohistochemical staining and quantitative analysis of the experimental group. A-D, Ki-67 immunohistochemical staining of prostate tumors from the experimental group (×400, scale bar 50 µm). Comparative analysis of the tumor nuclear fraction (E) and the Ki-67 area fraction (F) at days 9, 12, 15, and 18. Ki-67-positive expression was obvious in the tumor cell nucleus and was most intense on day 12. Data are reported as means and SD. Different letters indicate a statistically significant difference (P<0.05; ANOVA).
Comparison of tumor nuclear fraction and Ki-67 area fraction at each time point in the experimental group (n=5).
Changes in the Ki-67 area fraction
Ki-67-positive immunohistochemical staining revealed the presence of Ki-67-positive prostate tumor cell nuclei in the experimental group; the tumor cell nuclei were obviously stained at all time points, and the positive nuclei were the densest on the 12th day (Figure 6). The results of the quantitative analysis revealed that the Ki-67 area fraction gradually increased from days 9 to 12 and then began to decrease beginning on day 15 (Figure 6). There was no statistically significant difference between days 15 and 18 (P>0.05), while there was a statistically significant difference between the other time points (P<0.001) (Table 4).
Correlation analysis
MRI volume measurements of the mice in the experimental group revealed a strong positive correlation with the volume measured via calipers (r=0.994, P<0.001). The tumor volume was positively and negatively correlated with the tumor nuclear fraction and Ki-67 area fraction, respectively (r=0.895, r=-0.729, P<0.001). The tumor nuclear fraction was negatively correlated with the Ki-67 area fraction (r=-0.650, P<0.001; Figure 7).
Scatter plot of the correlation analysis of pathological indicators. A, Correlation analysis of the experimental group MRI volume and anatomical volume. B, Correlation analysis of the experimental group tumor nuclear fraction and anatomical volume. C, Correlation analysis of the experimental group Ki-67 area fraction and tumor anatomical volume. D, Correlation analysis of the experimental group Ki-67 area fraction and the tumor nuclear fraction. MRI: magnetic resonance imaging.
The tumor nuclear fraction showed a negative correlation with the ADC values at different b-value combinations, while it had a positive correlation with the EADC values at each b-value combination (Table 5).
The Ki-67 area fraction was positively correlated with the ADC values at different b-value combinations, while it was negatively correlated with the EADC values at each b-value combination (Table 5).
Discussion
Investigating prostate cancer requires appropriate animal models. The orthotopic prostate cancer model has a high tumor formation rate and notable repeatability, circumventing the prolonged and inconsistent tumor development time associated with genetically engineered mice. Moreover, this approach maintains the integrity of the mouse immune system. Concurrently, modeling prostate cancer in its native location facilitates accurate replication of the tumor growth and development microenvironment (22), a method that has been extensively employed in related research (23,24). In the present study, C57BL/6J mice were injected in situ with a mouse-derived prostate cancer cell suspension (RM-1) to establish a model. The experimental group demonstrated a 100% tumor formation rate with stable tumor growth, providing a solid foundation for subsequent investigations. After tumor implantation, the experimental group initially gained weight due to tumor growth. However, in the later stages, their condition deteriorated progressively, leading to a significant decrease in weight, which was markedly different from the stable weight observed in the control group.
DWI is a noninvasive functional MRI technique that evaluates water molecule diffusion in living tissues and enables the calculation of quantitative parameters such as ADC and EADC values. A brighter DWI signal correlates with lower ADC values, suggesting that a greater degree of water molecule diffusion is limited (12). The diffusion process of water molecules in living tissues is influenced by various structural characteristics of the tissue, including cell density, microvascular perfusion, and the presence of proteins that transport water molecules through the cell membrane. Research indicates that quantitative parameters derived from DWI may offer insights into pathological and biological features, such as cell quantity and tissue microstructure (17,25). These parameters may be effective indicators for evaluating histopathology changes and biological behavior (14- 16). However, previous studies did not monitor DWI quantitative parameters under multi b-value combinations during the dynamic process of prostate cancer development (17,18). Therefore, this study was designed to investigate the variations in DWI quantitative parameters across different b-values combinations during tumor progression to identify imaging indices that can accurately reflect changes in tumor histopathology.
In this investigation, DWI of prostate tumors revealed high signals, with the tumor signal intensifying over time. As the b-values increased, the DWI signal progressively diminished, and the signal-to-noise ratio of the image decreased. Nonetheless, the tumor DWI signal remained significantly greater than the surrounding tissues. The signal strength corresponded to the tumor area delineated on T1WI and T2WI, suggesting restricted water molecule diffusion within the tumor region. The ADC values at b-value combinations (b=0, 500/1000/1500/2000 s/mm2) in the experimental group were lower than the control group, demonstrating a decreasing trend concurrent with tumor growth. This aligns with the progressive brightening of the DWI signal over time, indicating an escalating restriction of water molecule diffusion in the tumor tissue as the tumor expands. This observation aligns with the findings of Hill et al. (26). In contrast, the EADC values at b-value combinations (b=0, 500/1000/1500/2000 s/mm2) for the experimental group exceeded the control group, presenting an upward trajectory with tumor progression that was consistent with DWI signal changes. This trend was inversely related to the changes in the ADC values. The EADC values decreased with increasing b-values at the same scanning time point, mirroring the reduction in the tumor signal visible to the naked eye. These findings suggest that EADC values can be used to quantify the DWI signal and furnish more precise information than visual assessment. This accuracy may be attributed to the method of EADC calculation, which involves dividing the high b-value DWI signal by the low b-value DWI signal to determine the signal intensity ratio of images at two distinct b-values. Since low b-value DWI partially reflects T2WI characteristics, EADC values can mitigate the T2WI projection effect in DWI images, yielding a more accurate DWI signal (27- 29).
In addition, we innovated in the combination of b-values and included two nonzero b-values in the mono-exponential model for quantitative parameter calculation. With the exceptions of b=1500, 2000 s/mm2, there was no significant difference in the ADC values generated by other b-value combinations between the experimental and control groups on day 9. This finding indicated that these combinations of b-values may be less effective in distinguishing early malignant tissue from normal tissue than the traditional method using b=0 s/mm2 for calculating quantitative parameters. Interestingly, unlike the ADC values, the EADC values at b=1500, 2000 s/mm2 also faced challenges distinguishing between early tumor and normal prostate tissue, which we attributed to differences in the information represented by the ADC and EADC values. In addition to the fundamental diffusion information of water molecules, the ADC values derived from higher b-value combinations also encompass data on aquaporins that mediate the diffusion and transportation of water molecules across cell membranes. However, the calculation of EADC values is contingent upon tissue signal strength, and further data are needed to ascertain its potential in offering more detailed micro-level information for discriminating between benign and malignant tissues.
For the selection of b-values, this study integrated the range suggested by PI-RADS v2.1 with b-values commonly employed in clinical practice, setting them at 0, 500, 1000, 1500, and 2000 s/mm2. These values represent low, medium, and high b-values for evaluating their applicability. Irrespective of the b-value combinations, the ADC and EADC values in the experimental group decreased and increased over time, respectively. The trend of changes across the various b-value combinations was consistent (Figure 4), suggesting that the dynamic alterations in the ADC and EADC values in the prostate cancer mouse model were stable. Moreover, we observed that as the tumor advanced, there was no statistically significant difference in the signal intensity ratio of the tumor at different time points. This finding indicates that the alterations in DWI signal intensity and ADC and EADC values were independent of changes in the T2WI signal, thus demonstrating that the trend of DWI quantitative parameters genuinely and effectively reflected the alteration of water molecule diffusion and microstructure in tumor tissues. In the control group, the repeated scanning results revealed that the ADC and EADC values of the different b-value combinations did not change with time. This stability might serve as an effective metric for the longitudinal assessment of prostate cancer onset and progression. Nonetheless, there was no marked difference in the ADC and EADC values at some adjacent time points in the experimental group, mirroring findings by Hectors et al. (25), who reported insignificant ADC value variances among patients with differing Gleason scores and between the pT2 and pT3 stages. This could be attributed to the influence of the tumor microenvironment on water molecule diffusion, where factors such as the presence of collagen, matrix swelling, and stromal disintegration within the tumor might impact the quantitative DWI parameters. Alternatively, the lack of significant differences could result from the study's limited sample size and brief data collection interval. These aspects will be addressed and improved upon in future research endeavors.
In the experimental group, prostate tumor cells rapidly proliferated, compromising the integrity of basal cells and invading normal glandular structures. This behavior aligns with hallmark features of prostate cancer, such as the degeneration of normal gland architecture, heightened tumor cell density, and irregular distribution of these cells (4). Typically, tumors at more advanced stages are characterized by greater cell density (17). Consequently, this study employed the tumor nuclear fraction observed via HE staining as a proxy for tumor cell density, which indirectly reflects alterations in tumor cell quantity. The findings revealed that the tumor nuclear fraction increased as the tumor progressed, mirroring the tumor's aberrant proliferation characteristic. Between days 15 and 18, no significant variance was observed in the tumor nuclear fraction, potentially attributable to the swift expansion of tumors in early stages compared with the decelerated growth in later stages. The occurrence of hemorrhage, necrosis, vacuolation, and hyalinization within the tumors may also contribute to the decrease in cell quantity observed (30). In addition, the dynamic expression of Ki-67 was detected by immunohistochemistry to reflect tumor proliferation. The results revealed that the tumors had proliferation ability at all time points. Quantitative analysis showed that the Ki-67 area fraction decreased on day 15 indicating that the proliferation ability of tumor cells decreased. This finding also explains why the tumor nuclear fraction measured by HE staining did not increase significantly between days 15 and 18.
Correlation analysis revealed a strong positive relationship between tumor volume determined through anatomical measurements and MRI within the experimental group. This finding underscores the efficacy of MRI as a noninvasive, accurate tool for monitoring tumor growth and evaluating treatment efficacy, which aligns with the findings of Ni et al. (31). Additionally, a positive correlation was observed between tumor nuclear fraction and tumor volume assessed by anatomical measurements, suggesting that MRI-based volume measurements not only precisely reflect changes in tumor size but also may infer shifts in tumor cell quantity. A negative correlation was detected between the tumor nuclear fraction and the ADC values across all b-value combinations, implying that tumor cell proliferation constricts the interstitial and intraluminal spaces, thereby restricting water molecule diffusion (32,33). Conversely, a positive relationship was found between the tumor nuclear fraction and the EADC values at all b-value combinations, indicating the ability of the EADC to accurately quantify the DWI signal intensity and, in conjunction with the ADC values, delineate water molecule diffusion. The tumor nuclear fraction significantly increased from days 9 to 15, followed by a minor increase from days 15 to 18, mirroring the decreasing trend of the ADC over the same period. These results suggest a biophysical relationship between the ADC and EADC values and the pathological structure of prostate cancer. In this investigation, the ADC values at b=1000, 1500 s/mm2, as well as the EADC values at b=0, 500 s/mm2, demonstrated the strongest correlations with the tumor nuclear fraction (r=-0.859 and 0.884, P<0.001, respectively). The correlation between the ADC values and the tumor nuclear fraction consistently remained close to -0.800 across all b-value combinations, indicating the robustness of the ADC values in indirectly reflecting cell density. However, for combinations of b-values other than 0 s/mm2, the correlation between EADC values and the tumor nuclear fraction gradually decreased with increasing adopted b-values and reached its minimum at b=1500, 2000 s/mm2 (r=0.562). This finding diverges from previous studies (34- 36) that showed that high b-values possess greater application value, potentially due to variations in research subjects, equipment, and evaluative indices.
In addition, the results revealed that Ki-67 expression was negatively correlated with tumor volume and nuclear fraction, positively correlated with the ADC value at each b-value combination, and negatively correlated with the EADC value. This may be due to the gradual increase in the density of tumor cells and the increased restriction of the water molecule diffusion in the tumor with progression. As a result, the proliferation ability of tumor cells, which urgently need water, blood, and nutrients to maintain growth, is also limited; therefore, the proliferation ability of tumor cells is reduced. In this study, the ADC and EADC values at b=500, 1000 s/mm2 were more strongly correlated with Ki-67 than the other b-value combinations (r=0.748, -0.701, P<0.001), and the ability of DWI-based quantitative parameters to characterize cell proliferation was less intuitive than that of cell density. In general, there is a close correlation between tumor nuclear fraction, Ki-67 expression, and DWI quantitative parameters, which has not been discussed in other studies and needs further exploration.
In conclusion, the longitudinal assessment of ADC and EADC values facilitated tumor progression monitoring, presenting a potential alternative imaging biomarker for quantifying the quantity and proliferation ability of prostate cancer cells. The ADC values at b=1000, 1500 s/mm2 and the EADC values at b=0, 500 s/mm2 may be most valuable for evaluating cell quantity. The ADC and EADC values at b=500, 1000 s/mm2 may be most valuable for evaluating tumor proliferation ability. Future research should explore the relationships between the quantitative parameters of DWI sequences and the characteristics of luminal spaces and more pathological indicators to investigate the potential of MRI for the non-invasive estimation of pathological alterations in prostate cancer tissues.
Supplementary Material
Click to view [pdf].
-
Funding
This work was supported by the Health Commission of Heilongjiang Province (Grant No. 20230909010030) and the Qiqihar Science and Technology Bureau Joint Guidance Project (Grant No. LSFGG-2022067).
References
-
1 Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Piãeros M, et al. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer. 2019;144:1941-53, doi: 10.1002/ijc.31937.
» https://doi.org/10.1002/ijc.31937 -
2 Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209-49, doi: 10.3322/caac.21660.
» https://doi.org/10.3322/caac.21660 -
3 Rebello RJ, Oing C, Knudsen KE, Loeb S, Johnson DC, Reiter RE, et al. Prostate cancer. Nat Rev Dis Primers. 2021;7:9, doi: 10.1038/s41572-020-00243-0.
» https://doi.org/10.1038/s41572-020-00243-0 -
4 Compérat E. Evolution of prostate cancer histopathology. Curr Opin Urol. 2019;29:587-92, doi: 10.1097/MOU.0000000000000669.
» https://doi.org/10.1097/MOU.0000000000000669 -
5 Sun Y, Reynolds HM, Wraith D, Williams S, Finnegan ME, Mitchell C, et al. Voxel-wise prostate cell density prediction using multiparametric magnetic resonance imaging and machine learning. Acta Oncol. 2018;57:1540-6, doi: 10.1080/0284186X.2018.1468084.
» https://doi.org/10.1080/0284186X.2018.1468084 -
6 Chen B, Tao J, Wu T, Feng H, Xie J, Zhong L, et al. Correlation between 18F-FDG PET-derived parameters and quantitative pathological characteristics of soft tissue sarcoma. Quant Imaging Med Surg. 2023;13:7842-53, doi: 10.21037/qims-23-412.
» https://doi.org/10.21037/qims-23-412 -
7 Sekhoacha M, Riet K, Motloung P, Gumenku L, Adegoke A, Mashele S. Prostate cancer review: genetics, diagnosis, treatment options, and alternative approaches. Molecules. 2022;27:5730, doi: 10.3390/molecules27175730.
» https://doi.org/10.3390/molecules27175730 -
8 Matuszczak M, Schalken JA, Salagierski M. Prostate cancer liquid biopsy biomarkers' clinical utility in diagnosis and prognosis. Cancers (Basel). 2021;13:3373, doi: 10.3390/cancers13133373.
» https://doi.org/10.3390/cancers13133373 -
9 Fernandes MC, Yildirim O, Woo S, Vargas HA, Hricak H. The role of MRI in prostate cancer: current and future directions. MAGMA. 2022;35:503-21, doi: 10.1007/s10334-022-01006-6.
» https://doi.org/10.1007/s10334-022-01006-6 -
10 Lee CH, Tan TW, Tan CH. Multiparametric MRI in active surveillance of prostate cancer: an overview and a practical approach. Korean J Radiol. 2021;22:1087-99, doi: 10.3348/kjr.2020.1224.
» https://doi.org/10.3348/kjr.2020.1224 -
11 Würnschimmel C, Chandrasekar T, Hahn L, Esen T, Shariat SF, Tilki D. MRI as a screening tool for prostate cancer: current evidence and future challenges. World J Urol. 2023;41:921-8, doi: 10.1007/s00345-022-03947-y.
» https://doi.org/10.1007/s00345-022-03947-y -
12 Ueno Y, Tamada T, Sofue K, Murakami T. Diffusion and quantification of diffusion of prostate cancer. Br J Radiol. 2022;95:20210653, doi: 10.1259/bjr.20210653.
» https://doi.org/10.1259/bjr.20210653 -
13 Tamada T, Ueda Y, Ueno Y, Kojima Y, Kido A, Yamamoto A. Diffusion-weighted imaging in prostate cancer. MAGMA. 2022;35:533-47, doi: 10.1007/s10334-021-00957-6.
» https://doi.org/10.1007/s10334-021-00957-6 -
14 Fliedner FP, Engel TB, El-Ali HH, Hansen AE, Kjaer A. Diffusion weighted magnetic resonance imaging (DW-MRI) as a non-invasive, tissue cellularity marker to monitor cancer treatment response. BMC Cancer. 2020;20:134, doi: 10.1186/s12885-020-6617-x.
» https://doi.org/10.1186/s12885-020-6617-x -
15 Qin X, Mu R, Zheng W, Li X, Liu F, Zhuang Z, et al. Comparison and combination of amide proton transfer magnetic resonance imaging and the apparent diffusion coefficient in differentiating the grades of prostate cancer. Quant Imaging Med Surg. 2023;13:812-24, doi: 10.21037/qims-22-721.
» https://doi.org/10.21037/qims-22-721 -
16 Guo R, Yang SH, Lu F, Han ZH, Yan X, Fu CX, et al. Evaluation of intratumoral heterogeneity by using diffusion kurtosis imaging and stretched exponential diffusion-weighted imaging in an orthotopic hepatocellular carcinoma xenograft model. Quant Imaging Med Surg. 2019;9:1566-78, doi: 10.21037/qims.2019.08.18.
» https://doi.org/10.21037/qims.2019.08.18 -
17 Hill DK, Heindl A, Zormpas-Petridis K, Collins DJ, Euceda LR, Rodrigues DN, et al. Non-invasive prostate cancer characterization with diffusion-weighted MRI: insight from in silico studies of a transgenic mouse model. Front Oncol. 2017;7:290, doi: 10.3389/fonc.2017.00290.
» https://doi.org/10.3389/fonc.2017.00290 -
18 Jung DC, Lee HJ, Seo JW, Park SY, Lee SJ, Lee JH, et al. Diffusion-weighted imaging of a prostate cancer xenograft model seen on a 7 Tesla animal MR scanner: comparison of ADC values and pathologic findings. Korean J Radiol. 2012;13:82-9, doi: 10.3348/kjr.2012.13.1.82.
» https://doi.org/10.3348/kjr.2012.13.1.82 -
19 Turkbey B, Rosenkrantz AB, Haider MA, Padhani AR, Villeirs G, Macura KJ, et al. Prostate imaging reporting and data system version 2.1: 2019 Update of prostate imaging reporting and data system version 2. Eur Urol. 2019;76:340-51, doi: 10.1016/j.eururo.2019.02.033.
» https://doi.org/10.1016/j.eururo.2019.02.033 -
20 Zhang J, Suo S, Liu G, Zhang S, Zhao Z, Xu J, et al. Comparison of monoexponential, biexponential, stretched-exponential, and kurtosis models of diffusion-weighted imaging in differentiation of renal solid masses. Korean J Radiol. 2019;20:791-800, doi: 10.3348/kjr.2018.0474.
» https://doi.org/10.3348/kjr.2018.0474 -
21 Zheng X, Chen Y, Zheng D, Xiao Y, Chen J, Pan J, et al. Diffusion kurtosis imaging and tumour microstructure for monitoring response to radiotherapy in human nasopharyngeal carcinoma xenografts. Jpn J Clin Oncol. 2020;50:548-55, doi: 10.1093/jjco/hyaa002.
» https://doi.org/10.1093/jjco/hyaa002 -
22 Valkenburg KC, Pienta KJ. Drug discovery in prostate cancer mouse models. Expert Opin Drug Discov. 2015;10:1011-24, doi: 10.1517/17460441.2015.1052790.
» https://doi.org/10.1517/17460441.2015.1052790 -
23 Xia H, Luo X, Yin W. Inhibition of prostate cancer growth by immunization with a GM-CSF-modified mouse prostate cancer RM-1 cell vaccine in a novel murine model. Oncol Lett. 2018;15:538-44, doi: 10.3892/ol.2017.7332.
» https://doi.org/10.3892/ol.2017.7332 -
24 Zhang D, Wang Y, Liu L, Li Z, Yang S, Zhao W, et al. Establishment and evaluation of ectopic and orthotopic prostate cancer models using cell sheet technology. J Transl Med. 2022;20:381, doi: 10.1186/s12967-022-03575-5.
» https://doi.org/10.1186/s12967-022-03575-5 -
25 Hectors SJ, Semaan S, Song C, Lewis S, Haines GK, Tewari A, et al. Advanced diffusion-weighted imaging modeling for prostate cancer characterization: correlation with quantitative histopathologic tumor tissue composition-a hypothesis-generating study. Radiology. 2018;286:918-28, doi: 10.1148/radiol.2017170904.
» https://doi.org/10.1148/radiol.2017170904 -
26 Hill DK, Kim E, Teruel JR, Jamin Y, Widerøe M, Søgaard CD, et al. Diffusion-weighted MRI for early detection and characterization of prostate cancer in the transgenic adenocarcinoma of the mouse prostate model. J Magn Reson Imaging. 2016;43:1207-17, doi: 10.1002/jmri.25087.
» https://doi.org/10.1002/jmri.25087 -
27 Provenzale JM, Engelter ST, Petrella JR, Smith JS, MacFall JR. Use of MR exponential diffusion-weighted images to eradicate T2 “shine-through” effect. AJR Am J Roentgenol. 1999;172:537-9, doi: 10.2214/ajr.172.2.9930819.
» https://doi.org/10.2214/ajr.172.2.9930819 -
28 Park SY, Kim CK, Park JJ, Park BK. Exponential apparent diffusion coefficient in evaluating prostate cancer at 3 T: preliminary experience. Br J Radiol. 2016;89:20150470, doi: 10.1259/bjr.20150470.
» https://doi.org/10.1259/bjr.20150470 -
29 Sprinkart AM, Marx C, Träber F, Block W, Thomas D, Schild H, et al. Evaluation of exponential ADC (eADC) and computed DWI (cDWI) for the detection of prostate cancer. Rofo. 2018;190:758-66, doi: 10.1055/a-0637-9980.
» https://doi.org/10.1055/a-0637-9980 -
30 Kroemer G, Pouyssegur J. Tumor cell metabolism: cancer's Achilles' heel. Cancer Cell. 2008;13:472-82, doi: 10.1016/j.ccr.2008.05.005.
» https://doi.org/10.1016/j.ccr.2008.05.005 -
31 Ni J, Bongers A, Chamoli U, Bucci J, Graham P, Li Y. In vivo 3D MRI measurement of tumour volume in an orthotopic mouse model of prostate cancer. Cancer Control. 2019;26:1073274819846590, doi: 10.1177/1073274819846590.
» https://doi.org/10.1177/1073274819846590 -
32 Kobus T, van der Laak JA, Maas MC, Hambrock T, Bruggink CC, Hulsbergen-van de Kaa CA, et al. Contribution of histopathologic tissue composition to quantitative MR spectroscopy and diffusion-weighted imaging of the prostate. Radiology. 2016;278:801-11, doi: 10.1148/radiol.2015142889.
» https://doi.org/10.1148/radiol.2015142889 -
33 Wichtmann BD, Zöllner FG, Attenberger UI, Schönberg SO. Multiparametric MRI in the diagnosis of prostate cancer: physical foundations, limitations, and prospective advances of diffusion-weighted MRI. Rofo. 2021;193:399-409, doi: 10.1055/a-1276-1773.
» https://doi.org/10.1055/a-1276-1773 -
34 Woo S, Suh CH, Kim SY, Cho JY, Kim SH. Head-to-head comparison between high- and standard-b-value DWI for detecting prostate cancer: a systematic review and meta-analysis. AJR Am J Roentgenol. 2018;210:91-100, doi: 10.2214/AJR.17.18480.
» https://doi.org/10.2214/AJR.17.18480 - 35 Werner S, Zinsser D, Esser M, Nickel D, Nikolaou K, Othman AE. Enhanced image processing using complex averaging in diffusion-weighted imaging of the prostate: the impact on image quality and lesion detectability. Diagnostics (Basel). 2023;13:2325.
-
36 Zhang G, Xu Z, Zheng J, Wang M, Ren J, Wei X, et al. Ultra-high b-value DWI in predicting progression risk of locally advanced rectal cancer: a comparative study with routine DWI. Cancer Imaging. 2023;23:59, doi: 10.1186/s40644-023-00582-7.
» https://doi.org/10.1186/s40644-023-00582-7
Edited by
-
Editor
Roberto César Lima-Júniorhttps://orcid.org/0000-0002-7033-655X














