Front. Oncol.Frontiers in OncologyFront. Oncol.2234-943XFrontiers Media S.A.10.3389/fonc.2025.1482112OncologyOriginal ResearchUnderstanding cellular proliferation activity in breast cancer using multi-compartment model of transverse relaxation time mapping on 3T MRINkondeKangwa Alex12†CheungSai Man1†SennNicholas3HeJiabao1*1Translational and Clinical Research Institute, Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne, United Kingdom2Department of Physics, School of Natural and Applied Sciences, Mulungushi University, Kabwe, Zambia3Institute of Medical Sciences, School of Medicine, Medical Sciences and Nutrition, University of Aberdeen, Aberdeen, United Kingdom
Edited by: Michael Albert Thomas, University of California, Los Angeles, United States
Reviewed by: Glaucia Maria Machado-Santelli, University of São Paulo, Brazil
Yen-Peng Liao, University of Texas Southwestern Medical Center, United States
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction
Precise understanding of proliferative activity in breast cancer holds significant value in the monitoring of neoadjuvant treatment, while current immunostaining of Ki-67 from biopsy or resected tumour suffers from partial sampling error. Multi-compartment model of transverse relaxation time has been proposed to differentiate intra- and extra-cellular space and biochemical environment but susceptible to noise, with recent development of Bayesian algorithm suggested to improve robustness. We hence hypothesise that intra- and extra-cellular transverse relaxation times using Bayesian algorithm might be sensitive to proliferative activity.
Materials and methods
Twenty whole tumour specimens freshly excised from patients with invasive ductal carcinoma were scanned on a 3 T clinical scanner. The overall transverse relaxation time was computed using a single-compartment model with the non-linear least squares algorithm, while intra- and extra-cellular transverse relaxation times were computed using a multi-compartment model with the Bayesian algorithm. Immunostaining of Ki-67 was conducted, yielding 9 and 11 cases with high and low proliferating activities respectively.
Results
For single-compartment model, there was a significant higher overall transverse relaxation time (p = 0.031) in high (83.55 ± 7.38 ms) against low (73.30 ± 11.30 ms) proliferating tumours. For multi-compartment model, there was a significant higher intra-cellular transverse relaxation time (p = 0.047) in high (73.52 ± 10.92 ms) against low (61.30 ± 14.01 ms) proliferating tumours. There was no significant difference in extra-cellular transverse relaxation time (p = 0.203) between high and low proliferating tumours.
Conclusions
Overall and Bayesian intra-cellular transverse relaxation times are associated with proliferative activities in breast tumours, potentially serving as a non-invasive imaging marker for neoadjuvant treatment monitoring.
neoadjuvant treatmentBayesianKi-67intra-cellularextra-cellularbiochemical environmentNewcastle University10.13039/501100000774section-in-acceptanceCancer Imaging and Image-directed InterventionsIntroduction
New pharmaceutical options and the improvement in screening services have improved the 5-year survival rate of breast cancer (1, 2), the most prevalent cancer worldwide (3), from 79.9% between 2000 – 2004 to 86.3% between 2014 – 2018 in the UK. Neoadjuvant therapies, often a combination of chemotherapy and hormonal therapy, are increasingly used to improve surgical outcomes by downstaging large tumours to facilitate breast conservation surgery (4). Neoadjuvant therapy is typically lengthy and costly, exposing non-responding patients to unnecessary potential adverse effects (5, 6). Rapid tumour growth, a central prognostic indicator of breast cancer, can be estimated using the proliferative activity marker Ki-67, highlighting nuclear protein expressed during cell division (7). However, Ki-67 using biopsy is invasive with narrow spatial coverage, suboptimal for accurately estimating intrinsically heterogeneous proliferative activity across the whole tumour during treatment (8, 9). Hence, an imaging marker of proliferative activity is central to effective neoadjuvant treatment monitoring in breast cancer.
PET tracer of 3’-deoxy-3’-[18F] fluorothymidine (FLT) has been shown to correlate with Ki-67 in breast cancer (10), but causes myelosuppression, peripheral neuropathy and nausea at high doses (7). Peri-tumoural texture features of kurtosis, skewness and entropy, based on maximum enhancement from dynamic contrast-enhanced (DCE) MRI are associated with proliferative activity in large cohorts (11, 12), but the primary sensitivity of DCE MRI towards angiogenesis compromises specificity (13). Apparent diffusion coefficient from diffusion weighted imaging (DWI) can distinguish between low and high proliferating breast tumours (11, 14), but is susceptible to biological noise (15). Both peri-tumoural lipid composition from chemical shift-encoded imaging (CSEI) (16) and perfusion fraction from intravoxel incoherent motion (IVIM) MRI (17) have shown correlations against proliferative activity, but the primary sensitivity to tissue composition and perfusion limits specificity (18). Hence, an imaging approach specific to proliferative activity with low susceptibility to biological noise is highly desirable.
Transverse relaxation time mapping has been shown to highlight tumour pathology (19), but only offers a crude picture of tissue microenvironment (20). Rapidly proliferating tumours are associated with elevated angiogenesis to support accelerated metabolic activities (21) exacerbating vascular permeability and fluid efflux (22), and the consequent increased free water pool and a diluted biochemical microenvironment (21, 23) lead to elevated transverse relaxation time in high proliferating breast tumours as observed in vivo (24). Intra-cellular transverse relaxation time, although more susceptible to noise compared to overall transverse relaxation time, is more specific to alterations in intra-cellular processes (25), potentially revealing tumour response to chemotherapy targeting at rapid cell division at an early stage (26). Recent application of multi-compartment model differentiates lumen, intra- and extra-cellular spaces, as demonstrated in luminal water imaging in prostate cancer (27), classification of adipocytic tumour in soft tissue (28), differentiation of layers in articular cartilage (29) and estimation of myelin volume in the brain (30). Although multi-compartment model has been suggested as a marker of proliferative activity in breast cancer (31), the susceptibility to erroneous attribution of noise (20, 29) significantly curtailed the clinical application. Bayesian algorithm, using probabilistic constraints between neighbouring image voxels, was subsequently developed to reduce susceptibility to noise (32–34), and we have recently demonstrated the value in improving imaging methods for neoadjuvant treatment monitoring (17). We hence hypothesise that multi-compartment model of transverse relaxation time mapping using Bayesian algorithm might be a sensitive marker of proliferative activity.
Materials and methods
To probe this hypothesis, a cross-sectional study was conducted on 20 breast tumour specimens freshly excised from patients. Multi-compartment model of transverse relaxation time maps using Bayesian algorithm was performed to derive the intra- and extra-cellular transverse relaxation times and the volume ratio of each breast tumour specimen (Figure 1). The study was approved by the North-West – Greater Manchester East Research Ethics Committee (REC Reference: 16/NW/0221), and signed written informed consent was obtained from all participants prior to the study.
Study design. Twenty specimens freshly excised from patients with invasive ductal carcinoma grades II and III were imaged on a clinical 3T MRI scanner using quantitative transverse relaxation time mapping, with subsequent histopathological confirmation of low and high cellular proliferation activity. The overall transverse relaxation time was computed voxel-wise using a single-compartment model, while the intra- and extra-cellular transverse relaxation times and volume ratio were calculated voxel-wise using a multi-compartment model.
Clinical procedure
Fresh tumour specimens were excised from twenty female patients (age mean 57 years, range 35 to 78 years) with invasive ductal carcinoma (IDC) (10 grade II and 10 grade III). The study included patients with tumour size larger than 1.5 cm in diameter on mammography. Patients with previous malignancies and undergoing chemotherapy or radiotherapy prior to surgery were excluded from the study. Upon excision, the tumour specimens were immobilised in a 10% buffered formalin solution in a sealed container before imaging. Routine histopathological examination was performed to determine the histological tumour diameter, grade and Nottingham Prognostic Index (NPI), with oestrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) (35). Tumour cellular proliferation activity marker Ki-67 (36) was assessed semi-quantitatively after single-batch immunostaining, with high Ki-67 indicating that more than 14% (37) of tumour cell nuclei staining positive above the background (38, 39), and there were 9 and 11 cases with high and low proliferative activities respectively.
Transverse relaxation time mapping
Quantitative transverse relaxation time images were acquired from each specimen on a clinical 3T MRI scanner (Achieva TX, Philips Healthcare, Netherlands) using a 32-channel receiver coil for signal detection and the body coil for uniform transmission. The images were acquired using a multishot gradient and spin echo (GRASE) sequence (40), with 24 echoes, initial echo time of 13 ms, echo spacing of 13 ms, repetition time of 9943 ms, field of view of 141 × 141 mm2, slice thickness of 2.2 mm, and image resolution of 2.2 × 2.2 mm2. The overall transverse relaxation time was computed voxel-wise using single-compartment model with the non-linear least squares method in MATLAB (R2021a, MathWorks Inc., Natick, MA, USA). The intra- and extra-cellular transverse relaxation times and volume ratio were calculated voxel-wise using the multi-compartment model with the Bayesian algorithm (33). The whole tumour was delineated for each specimen using MRIcron (v1.0.20190902, University of South Carolina, Colombia, USA) on conventional DWI images acquired at b value at 800 smm-2 using pulsed gradient spin echo (PGSE) sequence, with the necrotic regions excluded from the analysis. The overall transverse relaxation time, intra- and extra-cellular transverse relaxation times, and intra- and extra-cellular volumes were computed as the mean within the whole tumour from corresponding quantitative maps for each specimen, and the volume ratio subsequently calculated as the ratio between intra-cellular volume against combined intra- and extra-cellular volumes (27).
Statistical analysis
All statistical analysis was performed in the SPSS software 27.0 (IBM Corp, Armonk, NY, USA). Shapiro-Wilk test for normality was performed on all the outcome measures. Independent sample t-tests or Mann-Whitney U tests were conducted to compare the overall, intra- and extra-cellular transverse relaxation times and volume ratio between tumours with high and low proliferative activities. Pearson’s correlation tests were performed between the overall transverse relaxation time, intra-cellular transverse relaxation time and volume ratio against the tumour diameter, with Spearman’s rank correlation test performed between the extra-cellular transverse relaxation time against the tumour diameter. Spearman’s rank correlation tests were performed between the transverse relaxation times and volume ratio against NPI. A p-value < 0.05 was considered statistically significant.
Results
The patient demography and statistical findings are shown in Tables 1, 2 respectively. There were no significant differences in patient characteristics between the high and low Ki-67 groups.
Patient demography.
Demographic
All (n = 20)
Ki-67 expression
p-value
Low (n = 11)
High (n = 9)
Age (years)^
58 ± 14
63 ± 11
51 ± 12
0.052
Tumour Size (mm)^
24.9 ± 5.7
24.1 ± 6.6
25.9 ± 4.7
0.489
Nottingham Prognostic Index (NPI) †
4.44 (3.50 – 4.59)
3.70 (3.50 – 4.57)
4.44 (4.42 – 4.60)
0.199
Histological grade*
II
10
8
2
0.070
III
10
3
7
Lymphovascular Invasion (LVI)*
LVI Positive
12
7
5
1
LVI Negative
8
4
4
Necrosis*
Present
9
6
3
0.406
Absent
11
5
6
Tumour Infiltrating Lymphocytes (TILs)*
Low
13
7
6
1
Intermediate
4
3
1
0.591
High
3
1
2
0.566
^Independent samples t-test for continuous variables; †Mann-Whitney U-test for continuous non-parametric variables; *Fisher's exact test for categorical variables.
Breast cancer patients with high and low proliferative activity marker Ki-67 expression are shown for each group and the entire cohort. Values are expressed as mean and standard deviation for normally distributed data, and median and interquartile range for non-normally distributed data. Histopathological entries are expressed as the number of positive observations.
The overall transverse relaxation time (T2, MONO), intra- and extra-cellular transverse relaxation times (T2S, T2L) and volume ratio (f) between molecular subtypes, and low and high Ki-67 expression are shown. The correlations of T2, MONO, T2S, T2L and f against tumour diameter and Nottingham Prognostic Index (NPI) are also shown. Statistically significant differences (p < 0.05) are in bold and marked by '*'.
There was a significant higher overall transverse relaxation time (p = 0.031, Figure 2A) in high Ki-67 tumours (83.55 ± 7.38 ms) against low Ki-67 tumours (73.30 ± 11.30 ms). There was a significant higher intra-cellular transverse relaxation time (p = 0.047, Figure 2B) in high Ki-67 tumours (73.52 ± 10.92 ms) against low Ki-67 tumours (61.30 ± 14.01 ms). There was no significant difference (p = 0.203, Figure 2C) in the extra-cellular transverse relaxation time between high Ki-67 tumours (147.38 ± 8.84 ms) and low Ki-67 tumours (156.56 ± 19.16 ms). There was no significant difference (p = 0.073, Figure 2D) in volume ratio between high Ki-67 tumours (33.64 ± 8.33%) and low Ki-67 tumours (41.65 ± 10.08%).
Group comparison of relaxometry properties between low and high Ki-67 expression group in breast cancer (n = 11, 9). (A) Overall transverse relaxation times (T2, MONO) are significantly higher in high proliferating tumours than in low proliferating tumours, potentially indicating a diluted biochemical environment. (B) Intra-cellular transverse relaxation times (T2S) are significantly higher in high proliferating tumours than in low proliferating tumours, potentially indicating a dilution of the intra-cellular biochemical environment. (C) There was no significant difference in extra-cellular transverse relaxation times (T2L) between high and low proliferating tumours, suggesting biochemical homeostasis. (D) The near-significantly higher volume ratio (f), in low proliferating tumours than in high proliferating tumours, suggests possible physical homeostasis between intra- and extra-cellular fluid volumes. Each dot represents the mean ROI of the parameter, arranged by low (Ki-67 ≤ 14%) and high (Ki-67 > 14%) proliferative activity. Error bars represent the mean ± standard deviation. Two-tailed independent t-tests were conducted, with p-values provided; statistically significant results (p < 0.05) are marked by ‘*’.
Against tumour diameter, there was a significant negative correlation in extra-cellular transverse relaxation time (ρ = -0.50, p = 0.025, Figure 3A), but not in intra-cellular transverse relaxation time (p = 0.979), volume ratio (p = 0.922) and overall transverse relaxation time (p = 0.360) (Figure 3). There were no significant correlations between extra-cellular transverse relaxation time (p = 0.947), intra-cellular transverse relaxation time (p = 0.828), volume ratio (p = 0.473) and overall transverse relaxation time (p = 0.806) against NPI (Figure 4). Typical breast tumour specimens of high and low Ki-67 expression and the corresponding maps of overall and intra-cellular transverse relaxation times are shown in Figure 5, and maps of extra-cellular transverse relaxation time and volume ratio are shown in Supplementary Figure S1 (Supplementary Material 1).
Correlation analysis between relaxometry properties and breast tumour diameter (n = 20). (A) Extra-cellular transverse relaxation time (T2L). (B) Intra-cellular transverse relaxation time (T2S). (C) Volume ratio (f). (D) Overall transverse relaxation time (T2, MONO). Statistically significant correlations (p < 0.05) are marked by ‘*’.
Correlation analysis between relaxometry properties and Nottingham Prognostic Index (NPI) in breast cancer (n = 20). (A) Extra-cellular transverse relaxation time (T2L). (B) Intra-cellular transverse relaxation time (T2S). (C) Volume ratio (f). (D) Overall transverse relaxation time (T2, MONO).
High (Ki-67 > 14%) and low (Ki-67 ≤ 14%) proliferating breast tumour specimens in haematoxylin and eosin (H & E) staining and the corresponding transverse relaxation time maps. Sections are taken from the greatest dimension of the tumour diameter. Magnification, x20. (A) A high Ki-67 expression of 49.18%. (B) A low Ki-67 expression of 4.96%. (C) Overall transverse relaxation time (T2, MONO) map of the high Ki-67 specimen. (D) T2, MONO map of the low Ki-67 specimen. (E) Intra-cellular transverse relaxation time (T2S) map of the high Ki-67 specimen. (F) T2S map of the low Ki-67 specimen.
Discussion
In this study, we found significant differences between high and low proliferating tumours in overall and intra-cellular transverse relaxation times, but not in extra-cellular transverse relaxation time and volume ratio. We found a significant negative correlation of the extra-cellular transverse relaxation time against tumour diameter, but not in the remaining three relaxometry properties. No significant correlation of the four relaxometry properties against NPI was found.
The elevated overall transverse relaxation time in high proliferating tumours indicates a reduced rate of signal dissipation potentially due to a diluted biochemical environment (23). The rapid metabolic activity in high proliferating tumours results in increased free water pools in the tumour (21), while the enhanced angiogenesis exacerbates vascular permeability and fluid leakage into the tumour (22). Hence, the central characteristics of rapidly proliferating tumours (41), of increased metabolic activity and enhanced angiogenesis, contribute to the diluted biochemical environment in high proliferating tumours (42). The elevated intra-cellular transverse relaxation time in high proliferating tumours indicates a reduced rate of signal dissipation potentially due to a diluted intra-cellular biochemical environment (43). The rapid cell division in high proliferating tumours, as highlighted by Ki-67 (21), demands upregulated transportation of amino acids from the cytoplasmic space to sustain biosynthesis in the nucleus (44), leading to a diluted intra-cellular biochemical environment. The stronger statistical significance in overall transverse relaxation time compared to intra-cellular transverse relaxation time might be the result of lower intra-cellular signal as a fraction of the overall signal, but might also indicate the presence of feedforward and feedback pathways. Angiogenesis exacerbates vascular permeability to support accelerated metabolic activity, compounding fluid efflux in both the intra- and extra-cellular compartments (22), potentially leading to a more significant elevation in the overall transverse relaxation time in high proliferating tumours (23). In contrast, the elevated transportation of amino acids into the nucleus from the cytoplasmic space of the intra-cellular compartment to support accelerated biosynthesis in rapidly proliferating tumours partially cancels out the fluid efflux, potentially leading to a more subdued increase in the intra-cellular transverse relaxation time (44). Chemotherapies, targeting at rapid cell division, significantly alter the intra-cellular biochemical environment (26), hence sensitive imaging markers might contribute to early response monitoring (25). The lack of significant difference in the extra-cellular transverse relaxation time indicates an undisturbed rate of signal dissipation potentially reflecting biochemical homeostasis, although unnaturalised free radicals as by-products from aerobic glycolysis (45) are actively expelled to extra-cellular volume (46), promoting fibrosis (47), angiogenesis (22) and reduction of free water (48). The lack of significant difference in the volume ratio reflects the physical homeostasis of intra- and extra-cellular fluid volumes (49) critical for minimal disruption of osmotic pressure and transportation across cellular compartments in tumours (21). Vesicles, carrying oxidative stress signal transmitters, are released in the extra-cellular volume, but are actively transported back into the intra-cellular volume to modulate oxidative stress and fluid imbalance (50). The near-significant difference in the volume ratio compared to the significant difference in overall and intra-cellular transverse relaxation times might be the result of intrinsic lower effect size compared to measurement error, but might potentially indicate the secondary effects of structural changes (51) following biochemical alteration (52). Protective homeostatic mechanism (49), modulating fluid exchange through vesicles between the compartments (50), may dampen the effects of proliferative activities on the volume ratio with the same trend but weaker significance against transverse relaxation times (53).
The negative correlation between extra-cellular transverse relaxation time and tumour diameter might be due to the increased cellularity and reduced extra-cellular free water pools associated with larger tumours (48). The significant negative correlation between extra-cellular transverse relaxation time and tumour size indicates a more concentrated biochemical environment or less free water at larger tumour size. Larger tumours suffer from hypoxia due to inadequate vascular supply (54) and extra-cellular matrix remodelling (55), and the increase in the macromolecules to water ratio in the extra-cellular compartment (56) lead to a reduction in extra-cellular transverse relaxation time (57). There was no significant correlation between overall transverse relaxation time, intra-cellular transverse relaxation time and volume ratio with tumour diameter, reflecting the insensitivity of the three relaxometry properties to morphological tumour size. The lack of significant correlation between overall transverse relaxation time, intra-cellular and extra-cellular transverse relaxation times and volume ratio against NPI indicated that the four relaxometry properties, although potential treatment monitoring markers, might not be sensitive to recurrence and metastatic risks revealed by NPI (58).
Although an understanding of the heterogeneity in the overall and intra-cellular transverse relaxation times required thorough quantitative analysis across the whole tumour using appropriate texture features, however conjectural visual exploration might serve to develop future hypothesis. Intra-cellular transverse relaxation time showed sparse focal elevation at the edges of high proliferating tumours but more homogeneity in low proliferating tumours, while overall transverse relaxation time showed more homogeneous elevation in high proliferating tumours but exhibited bands of elevation at the edges of low proliferating tumours (Figure 5). The focal increase in intra-cellular free water content, manifested as a lengthening of intra-cellular transverse relaxation time, in high proliferating tumours might reflect the elevated transportation of amino acids into the nucleus from the cytoplasmic space to sustain biosynthesis in the nucleus at the advancing edge of the tumour (44). The bands of increased free water content, manifested as the elevation of overall transverse relaxation time, at the edge of low proliferating tumours might reflect the presence of angiogenesis and oedema concentrated at the tumour margin (22, 42), in contrast to their widespread presence across the whole tumour in high proliferating tumours (24). Precise spatial correlation between imaging and pathology, demanding co-localisation at a fine cellular resolution, was not feasible during this study, and should be performed in future investigations.
To our knowledge, this is the first study investigating the relaxation properties as a marker of breast tumour proliferative activity using Bayesian algorithm in conjunction with multi-compartment model to mitigate the impact from noise. This study only recruited patients with large breast tumours to take into account the reliance of Bayesian algorithm on neighbouring voxels, and necrotic tumour cores were excluded from analysis to avoid confounding factors from non-viable tissue in large tumours. The tumours were freshly excised to eliminate biological noise and further improve analysis accuracy, and imaging was performed immediately to avoid the alteration to relaxation properties from formalin submersion (59). The cohort size was small due to inclusion criteria on large tumour to accommodate novel imaging methods, and future large cohort studies on patients are critical to fully understand the value of the imaging marker. Future simulation study on the variability and high resolution imaging from ultra-high field MRI can support the establishment of a confidence range in sensitivity and an estimation of the highest effective image resolution, allowing studies to be carried out at a higher image resolution and in smaller tumours.
Conclusions
Overall and Bayesian intra-cellular transverse relaxation times are associated with proliferative activities in breast tumours, potentially serving as a non-invasive imaging marker for neoadjuvant treatment monitoring.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by North-West - Greater Manchester East Research Ethics Committee (REC Reference: 16/NW/0221). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This project was funded by NHS Grampian Endowment Research Fund (15/1/052). KAN’s PhD study is supported by the Commonwealth Scholarship Commission UK (ZMCS-2020-870). SMC was funded by Cancer Research UK (C68628/A28312). The funding sources were not involved in the study design, in the collection, analysis and interpretation of data, in the writing of the report nor in the decision to submit the article for publication.
Acknowledgments
The authors would like to thank Dr. Matthew Clemence (Philips Healthcare Clinical Science, UK) for clinical scientist support, Ms. Bolanle Brikinns for patient recruitment support and Ms. Dawn Younie for logistic support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2025.1482112/full#supplementary-material
ReferencesJohnSBroggioJ. Cancer survival in England adult, stage at diagnosis and childhood - patients followed up to 2018 (2019). Office for National Statistics. Available online at: https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/bulletins/cancersurvivalinengland/stageatdiagnosisandchildhoodpatientsfollowedupto2018 (Accessed Feb 23, 2024).Healthcare Quality Indicators: Cancer Care. OECD Data Explorer- Archive (2024). Available online at: https://stats.oecd.org/Index.aspx?QueryId=51882 (Accessed Apr 30, 2024).FerlayJColombetMSoerjomataramIParkinDMPiñerosMZnaorA. Cancer statistics for the year 2020: An overview. Int J Cancer. (2021) 149:778–89. doi: 10.1002/ijc.33588AsselainBBarlowWBartlettJBerghJBergsten-NordströmEBlissJ. Long-term outcomes for neoadjuvant versus adjuvant chemotherapy in early breast cancer: meta-analysis of individual patient data from ten randomised trials. Lancet Oncol. (2018) 19:27–39. doi: 10.1016/S1470-2045(17)30777-5BaltussenJCde GlasNAvan HolsteinYvan der ElstMTrompetSUit den BoogaardA. Chemotherapy-related toxic effects and quality of life and physical functioning in older patients. JAMA Netw Open. (2023) 6:e2339116. doi: 10.1001/jamanetworkopen.2023.39116National Cancer Institute. Hormone Therapy for Breast Cancer Fact Sheet (2022). Available online at: https://www.cancer.gov/types/breast/breast-hormone-therapy-fact-sheetwhat-are-the-side-effects-of-hormone-therapy (Accessed Feb 26, 2024).BeresfordMJWilsonGDMakrisA. Measuring proliferation in breast cancer: Practicalities and applications. Breast Cancer Res. (2006) 8(6):216. doi: 10.1186/bcr1618BilousM. Breast core needle biopsy: issues and controversies. Mod Pathol. (2010) 23 Suppl 2:S36–45. doi: 10.1038/modpathol.2010.34ZhangYJWeiLLiJZhengYQLiXR. Status quo and development trend of breast biopsy technology. Gland Surg. (2013) 2:15–24. doi: 10.3978/j.issn.2227-684X.2013.02.01KennyLMVigushinDMAl-NahhasAOsmanSLuthraSKShoushaS. Quantification of cellular proliferation in tumor and normal tissues of patients with breast cancer by [18F]fluorothymidine-positron emission tomography imaging: evaluation of analytical methods. Cancer Res. (2005) 65:10104–12. doi: 10.1158/0008-5472.CAN-04-4297EL-MetwallyDMonierDHassanAHelalAM. Preoperative prediction of Ki-67 status in invasive breast carcinoma using dynamic contrast-enhanced MRI, diffusion-weighted imaging and diffusion tensor imaging. Egyptian J Radiol Nucl Med. (2023) 54:62. doi: 10.1186/s43055-023-01007-yHeXZhouJYeSChengQMiaoHXuN. Differences in tumour heterogeneity based on dynamic contrast-enhanced MRI between tumour and peritumoural stroma for predicting Ki-67 status of invasive ductal carcinoma. Clin Radiol. (2021) 76:470.e13–470.e22. doi: 10.1016/j.crad.2020.12.008MannRMKuhlCKKinkelKBoetesC. Breast MRI: guidelines from the european society of breast imaging. Eur Radiol. (2008) 18:1307–18. doi: 10.1007/s00330-008-0863-7ZhangLHaoJGuoJZhaoXYinX. Predicting of ki-67 expression level using diffusion-weighted and synthetic magnetic resonance imaging in invasive ductal breast cancer. Breast J. (2023) 2023:1–8. doi: 10.1155/2023/6746326FowlerAMMankoffDAJoeBN. Imaging neoadjuvant therapy response in breast cancer. Radiology. (2017) 285:358–75. doi: 10.1148/radiol.2017170180ChanKSCheungSMSennNHusainEMasannatYHeysS. Peri-tumoural spatial distribution of lipid composition and tubule formation in breast cancer. BMC Cancer. (2022) 22:285. doi: 10.1186/s12885-022-09362-1CheungSMWuWSSennNSharmaRMcGoldrickTGagliardiT. Towards detection of early response in neoadjuvant chemotherapy of breast cancer using Bayesian intravoxel incoherent motion. Front Oncol. (2023) 13:1277556. doi: 10.3389/fonc.2023.1277556Szubert-FranczakAENaduk-OstrowskaMPasiczKPodgórskaJSkrzyńskiWCieszanowskiA. Intravoxel incoherent motion magnetic resonance imaging: basic principles and clinical applications. Pol J Radiol. (2020) 85:624–35. doi: 10.5114/pjr.2020.101476MerchantTEThelissenGRPde GraafPWNieuwenhuizenCWEAKievitHCEden OtterW. Application of a mixed imaging sequence for MR imaging characterization of human breast disease. Acta Radiol. (1993) 34:356–61. doi: 10.1080/02841859309173257AbabnehZBeloeilHBerdeCBGambarotaGMaierSEMulkernRV. Biexponential parameterization of diffusion and T2 relaxation decay curves in a rat muscle edema model: Decay curve components and water compartments. Magn Reson Med. (2005) 54:524–31. doi: 10.1002/mrm.20610ZhuJThompsonCB. Metabolic regulation of cell growth and proliferation. Nat Rev Mol Cell Biol. (2019) 20:436–50. doi: 10.1038/s41580-019-0123-5FoxSBGeneraliDGHarrisAL. Breast tumour angiogenesis. Breast Cancer Res. (2007) 9:216. doi: 10.1186/bcr1796MarquesMPMde CarvalhoALMBMartinsCBSilvaJDSarterMGarcía SakaiV. Cellular dynamics as a marker of normal-to-cancer transition in human cells. Sci Rep. (2023) 13:21079. doi: 10.1038/s41598-023-47649-wMatsudaMKidoTTsudaTOkadaKShiraishiYSuekuniH. Utility of synthetic MRI in predicting the Ki-67 status of oestrogen receptor-positive breast cancer: a feasibility study. Clin Radiol. (2020) 75:398.e1–8. doi: 10.1016/j.crad.2019.12.021TőkésTTőkésAMSzentmártoniGKisznerGMadarasLKulkaJ. Expression of cell cycle markers is predictive of the response to primary systemic therapy of locally advanced breast cancer. Virchows Archiv. (2016) 468:675–86. doi: 10.1007/s00428-016-1925-xQuFGuYXueMHeMZhouFWangG. Impact of therapy on cancer metabolism in high-risk localized prostate cancer treated with neoadjuvant docetaxel and androgen deprivation therapy. Prostate. (2021) 81:560–71. doi: 10.1002/pros.24134SabouriSChangSDSavdieRZhangJJonesECGoldenbergSL. Luminal water imaging: a new MRI T2 mapping technique for prostate cancer diagnosis. Radiology. (2017) 284:451–9. doi: 10.1148/radiol.2017161687NikiforakiKIoannidisGSLagoudakiEManikisGHde BreeEKarantanasA. Multiexponential T2 relaxometry of benign and Malignant adipocytic tumours. Eur Radiol Exp. (2020) 4:45. doi: 10.1186/s41747-020-00175-0ReiterDAIrrechukwuOLinPCMoghadamSvon ThaerSPleshkoN. Improved MR-based characterization of engineered cartilage using multiexponential T2 relaxation and multivariate analysis. NMR Biomed. (2012) 25:476–88. doi: 10.1002/nbm.1804BjörkMZachariahDKullbergJStoicaP. A multicomponent T2 relaxometry algorithm for myelin water imaging of the brain. Magn Reson Med. (2016) 75:390–402. doi: 10.1002/mrm.25583RomanovRLadutenkoKZubkovMAndreychenkoA. (2022). Method of moments for relaxation-based signal separation in magnetic resonance breast cancer detection, in: 2022 IEEE International Multi-Conference on Engineering, Computer and Information Sciences (SIBIRCON), Yekaterinburg, Russia, 2022 Nov 11 – 13. pp. 630–3. New Jersey, USA: IEEE.BretthorstGLHuttonWCGarbowJRAckermanJJH. Exponential parameter estimation (in NMR) using bayesian probability theory. Concepts Magn Reson Part A. (2005) 27A:55–63. doi: 10.1002/cmr.a.20043BouhraraMReiterDASpencerRG. Bayesian analysis of transverse signal decay with application to human brain. Magn Reson Med. (2015) 74:785–802. doi: 10.1002/mrm.25457GustafssonOMonteliusMStarckGLjungbergM. Impact of prior distributions and central tendency measures on Bayesian intravoxel incoherent motion model fitting. Magn Reson Med. (2018) 79:1674–83. doi: 10.1002/mrm.26783ElstonCWEllisIO. pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. Histopathology. (1991) 19:403–10. doi: 10.1111/j.1365-2559.1991.tb00229.xTuominenVJRuotoistenmäkiSViitanenAJumppanenMIsolaJ. ImmunoRatio: a publicly available web application for quantitative image analysis of estrogen receptor (ER), progesterone receptor (PR), and Ki-67. Breast Cancer Res. (2010) 12:R56. doi: 10.1186/bcr2615GoldhirschAWinerEPCoatesASGelberRDPiccart-GebhartMThürlimannB. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Ann Oncol. (2013) 24:2206–23. doi: 10.1093/annonc/mdt303YerushalmiRWoodsRRavdinPMHayesMMGelmonKA. Ki67 in breast cancer: prognostic and predictive potential. Lancet Oncol. (2010) 11:174–83. doi: 10.1016/S1470-2045(09)70262-1DenkertCLoiblSMüllerBMEidtmannHSchmittWDEiermannW. Ki67 levels as predictive and prognostic parameters in pretherapeutic breast cancer core biopsies: a translational investigation in the neoadjuvant GeparTrio trial. Ann Oncol. (2013) 24:2786–2793. doi: 10.1093/annonc/mdt350ChuMLChienCPWuWCChungHW. Gradient- and spin-echo (GRASE) MR imaging: a long-existing technology that may find wide applications in modern era. Quant Imaging Med Surg. (2019) 9:1477–84. doi: 10.21037/qims.2019.09.13AmanoMFujitaSTakeiNSanoKWadaASatoK. Feasibility of quantitative MRI using 3D-QALAS for discriminating immunohistochemical status in invasive ductal carcinoma of the breast. J Magnetic Resonance Imaging. (2023) 58:1752–9. doi: 10.1002/jmri.28683ParkerDOuld IsmailAAWolfRBremSAlexanderSHodgesW. Freewater estimatoR using iNtErpolated iniTialization (FERNET): Characterizing peritumoral edema using clinically feasible diffusion MRI data. PloS One. (2020) 15:e0233645. doi: 10.1371/journal.pone.0233645PilatusUShimHArtemovDDavisDVan ZijlPCMGlicksonJD. Intracellular volume and apparent diffusion constants of perfused cancer cell cultures, as measured by NMR. Magn Reson Med. (1997) 37:825–32. doi: 10.1002/mrm.1910370605PalmWThompsonCB. Nutrient acquisition strategies of mammalian cells. Nature. (2017) 546:234–42. doi: 10.1038/nature22379YangNLiHCaoCZhaoLSongXWangW. Tumor microenvironment-activated theranostic nanoreactor for NIR-II Photoacoustic imaging-guided tumor-specific photothermal therapy. Fundam Res. (2024) 4:178–87. doi: 10.1016/j.fmre.2022.04.021UchiyamaHZhaoQLHassanMAAndocsGNojimaNTakedaK. EPR-spin trapping and flow cytometric studies of free radicals generated using cold atmospheric argon plasma and X-ray irradiation in aqueous solutions and intracellular milieu. PloS One. (2015) 10:e0136956. doi: 10.1371/journal.pone.0136956MorryJNgamcherdtrakulWYantaseeW. Oxidative stress in cancer and fibrosis: Opportunity for therapeutic intervention with antioxidant compounds, enzymes, and nanoparticles. Redox Biol. (2017) 11:240–253. doi: 10.1016/j.redox.2016.12.011DewhirstMWCaoYMoellerB. Cycling hypoxia and free radicals regulate angiogenesis and radiotherapy response. Nat Rev Cancer. (2008) 8:425–37. doi: 10.1038/nrc2397MorishitaKWatanabeKIchijoH. Cell volume regulation in cancer cell migration driven by osmotic water flow. Cancer Sci. (2019) 110:2337–47. doi: 10.1111/cas.14079DoradoEDoriaMLNagelkerkeAMcKenzieJSManeta-StavrakakiSWhittakerTE. Extracellular vesicles as a promising source of lipid biomarkers for breast cancer detection in blood plasma. J Extracell Vesicles. (2024) 13:e12419. doi: 10.1002/jev2.12419YoukJHSonEJChungJKimJAKimEk. Triple-negative invasive breast cancer on dynamic contrast-enhanced and diffusion-weighted MR imaging: comparison with other breast cancer subtypes. Eur Radiol. (2012) 22:1724–34. doi: 10.1007/s00330-012-2425-2AkgülSPatchAMD’souzaRCJMukhopadhyayPNonesKKempeS. Intratumoural heterogeneity underlies distinct therapy responses and treatment resistance in glioblastoma. Cancers. (2019) 11:190. doi: 10.3390/CANCERS11020190GuoMPegoraroAFMaoAZhouEHAranyPRHanY. Cell volume change through water efflux impacts cell stiffness and stem cell fate. PNAS. (2017) 114(41): E8618–E8627. doi: 10.1073/pnas.1705179114SeršaIBajdFSavarinMJesenkoTČemažarMSeršaG. Multiparametric high-resolution MRI as a tool for mapping of hypoxic level in tumors. Technol Cancer Res Treat. (2018) 17:1–8. doi: 10.1177/1533033818797066YangZXueFLiMZhuXLuXWangC. Extracellular matrix characterization in gastric cancer helps to predict prognosis and chemotherapy response. Front Oncol. (2021) 11:753330. doi: 10.3389/fonc.2021.753330DrainAPZahirNNortheyJJZhangHHuangPJMallerO. Matrix compliance permits NF-κB activation to drive therapy resistance in breast cancer. J Exp Med. (2021) 218(5):e20191360. doi: 10.1084/jem.20191360SouzaFFD’AmatoGJonczakEECostaPTrentJCRosenbergAE. MRI T2 mapping assessment of T2 relaxation time in desmoid tumors as a quantitative imaging biomarker of tumor response: preliminary results. Front Oncol. (2023) 13:1286807. doi: 10.3389/fonc.2023.1286807Diagnostics Advisory committee members and NICE project team. Tumour profiling tests to guide adjuvant chemotherapy decisions in early breast cancer Diagnostics guidance (2024). NICE guidance. Available online at: www.nice.org.uk/guidance/dg58 (Accessed Aug 15, 2024).EbataKNorikiSInaiKKimuraH. Changes in magnetic resonance imaging relaxation time on postmortem magnetic resonance imaging of formalin-fixed human normal heart tissue. BMC Med Imaging. (2021) 21:134. doi: 10.1186/s12880-021-00666-5