- Open Access
Collagen fiber features and COL1A1: are they associated with elastic parameters in breast lesions, and can COL1A1 predict axillary lymph node metastasis?
BMC Cancer volume 22, Article number: 1004 (2022)
This study aimed to explore whether collagen fiber features and collagen type I alpha 1 (COL1A1) are related to the stiffness of breast lesions and whether COL1A1 can predict axillary lymph node metastasis (LNM).
Ninety-four patients with breast lesions were consecutively enrolled in the study. Amongst the 94 lesions, 30 were benign, and 64 were malignant (25 were accompanied by axillary lymph node metastasis). Ultrasound (US) and shear wave elastography (SWE) were performed for each breast lesion before surgery. Sirius red and immunohistochemical staining were used to examine the shape and arrangement of collagen fibers and COL1A1 expression in the included tissue samples. We analyzed the correlation between the staining results and SWE parameters and investigated the effectiveness of COL1A1 expression levels in predicting axillary LNM.
The optimal cut-off values for Emax, Emean, and Eratio for diagnosing the benign and malignant groups, were 58.70 kPa, 52.50 kPa, and 3.05, respectively. The optimal cutoff for predicting axillary LNM were 107.5 kPa, 85.15 kPa, and 3.90, respectively. Herein, the collagen fiber shape and arrangement features in breast lesions were classified into three categories. One-way analysis of variance (ANOVA) showed that Emax, Emean, and Eratio differed between categories 0, 1, and 2 (P < 0.05). Meanwhile, elasticity parameters were positively correlated with collagen categories and COL1A1 expression. The COL1A1 expression level > 0.145 was considered the cut-off value, and its efficacy in benign and malignant breast lesions was 0.808, with a sensitivity of 66% and a specificity of 90%. Furthermore, when the COL1A1 expression level > 0.150 was considered the cut-off, its efficacy in predicting axillary LNM was 0.796, with sensitivity and specificity of 96% and 59%, respectively.
The collagen fiber features and expression levels of COL1A1 positively correlated with the elastic parameters of breast lesions. The expression of COL1A1 may help diagnose benign and malignant breast lesions and predict axillary LNM.
The incidence of female breast cancer in China is increasing annually, seriously affecting the quality of life of women . The axillary lymph nodes are early metastatic sites of breast cancer, and proper preoperative evaluation is key to the formulation of rational treatment regimens. Furthermore, stiffness is an important characteristic of tissues and organs and is the basis for clinicians to diagnose breast abnormalities. The advent of shear wave elastography (SWE) offers the possibility of achieving quantitative judgment of tissue stiffness and can facilitate appropriate degradation of masses with low suspicion of breast imaging reporting and data system (BI-RADS) category 4A at soft borders . Recently, several studies have reported that SWE showed good performance in assessing benign and malignant lesions [3,4,5]. The stromal reaction that occurs because of invasive behavior results in abnormal tumor-associated collagen deposition, potentially leading to increased elastic modulus values in breast cancer . Previous studies have shown that the elastic parameters of SWE in breast tumors can be used to predict the lymph node status of invasive breast cancer (IBC); however, the pathological basis is unclear .
Previous studies on SWE and breast tissue stroma have focused on the relationship between the quantity of collagen or fibrosis and elastic parameters . However, the influence of the shape and arrangement of collagen fibers on the stiffness of the extracellular matrix (ECM) has rarely been considered. Malignant lesions are stiffer than benign lesions because of the presence of desmoplastic ECM, with type I collagen being a major structural component . The primary function of collagen type I is to support organs and tissues and cause tissue to form tension . Additionally, collagen type I is closely related to cell growth, proliferation, differentiation, and inflammation . Moreover, collagen type I consists of a triple-helix structure consisting of α1(I) chains (COL1A1) and α2(I) chains (COL1A2). These two polypeptide chains are encoded by different genes and are generally synthesized at a ratio of 2:1 . Nevertheless, little is known about the functional role of COL1A1 in elastography feature differences and axillary lymph node metastasis (LNM) of breast lesions. This study aimed to explore the relationship between the shape and arrangement of collagen fibers and the elastic parameters of breast lesions, which were evaluated using SWE. Additionally, we evaluated the role of COL1A1 expression in the differential diagnosis of benign and malignant breast lesions and analyzed the feasibility of COL1A1 in predicting axillary LNM.
From May 2021 to December 2021, 94 patients with breast lesions were consecutively enrolled in this prospective study, and the 94 breast lesions were examined using conventional ultrasound (US) and SWE before surgical resection. The inclusion criterion was hospitalized women with clear pathological diagnosis. In contrast, the exclusion criteria included patients who underwent biopsy or neoadjuvant chemotherapy before US and SWE examinations and patients with pathological multifocal IBC. The pathology results obtained by surgery were regarded as gold criteria. Furthermore, this trial was conducted in accordance with the Declaration of Helsinki (revised in 2013). The study was approved by the Ethics Committee of the Chinese PLA General Hospital (No. S2021-683–01), and informed consent was obtained from all participants. Based on postoperative pathology, the patients were divided into benign and malignant groups. The malignant group was divided into the negative axillary lymph node metastasis (LNM-) group and the positive axillary lymph node metastasis (LNM +) group.
Preoperative US and SWE examination
SWE was performed after conventional US examination using the AixPlorer US system (SuperSonic Imagine, Aix-en-Provence, France) with an L14-5 linear array probe operating at 4–15 MHz. Furthermore, the same sonographer, with more than 10 years of clinical experience, performed all US and SWE examinations for the breast masses. The probe was applied as lightly as possible to avoid excessive pressure and was maintained as steadily as possible for at least 10–20 s during elastic image acquisition. Additionally, the participants were asked to hold their breath to prevent motion artifacts. Furthermore, to enable a comparison with histopathological data, the maximum imaging section of the lesion was obtained to measure the stiffness of the breast lesion. Conventional US was used to obtain the maximum imaging section of the breast lesion and then start the SWE check, freezing, and measurement. Conversely, attention was paid while adjusting the region of interest (ROI) to include any halo or inelastic abnormal edge lesions while minimizing any normal tissue. The Eratio was obtained by placing the ROI in the normal breast tissue at the same depth as the lesion. This process was repeated three times, and the final values of Emax, Emean, Emin, and Eratio were the average values of these three repetitions.
Collagen staining and categorization
Sirius red staining was used to assess the shape and arrangement of collagen fibers, and immunohistochemical staining was used to examine the expression of COL1A1. All staining was performed on paraffin sections made from the largest plane of the surgical specimen of the breast tumor. An Olympus BX53P microscope equipped with a polarized filter was used to examine the collagen fibers at 400 × magnification. In this study, the collagen fiber shape and arrangement feature in breast lesions were classified into three categories. The definitions were described according to the following descriptions: category 0, wavy collagen fibers similar to collagen fibers in normal breast tissue; category 1, taut parallel collagen fibers around tumor nests; category 2, straightened and aligned collagen fibers tending to be perpendicular to the tumor boundary [13,14,15,16]. Following the classification principle we adopted, five fields of view were randomly selected. Collagen category was determined when the characteristic collagen arrangement exceeded 80% within the five fields. Analysis and classification of collagen fibers in each section were carried out by a physician with five years of experience in histopathology, who was blinded to the SWE and surgical pathology results.
COL1A1 expression was detected by immunohistochemistry (IHC). IHC staining results were analyzed using Image-Pro Plus 6.0 software (Media Cybernetics, Rockville, MD, USA). Five areas of interest were taken from each slice under a 200 × field of view and photographed. After measuring the integrated optical density (IOD) and area, the average optical density (AOD = IOD/Area) was used to evaluate the expression of COL1A1.
All statistical analyses were performed using SPSS (version 25.0; IBM Corp., Armonk, NY, USA). The measurement data are expressed as mean ± standard deviation (SD) and count data as percentages (%). One-way analysis of variance (ANOVA) was used to compare the data between the three groups. The chi-square test was used to compare the rates between the groups. Furthermore, receiver operating characteristic (ROC) curves were constructed for SWE values and COL1A1 expression levels to differentiate benign from malignant breast lesions and predict axillary LNM. The corresponding sensitivity and specificity were recorded (95% confidence interval).
Additionally, the correlation between COL1A1 expression level and elastic parameters in breast lesions was evaluated using Pearson correlation analysis, and the coefficients were defined as follows: < 0.40 indicated poor reliability; 0.40–0.75 indicated good reliability; and > 0.75 indicated excellent reliability . In addition, Kendall's tau-b correlation coefficient was used to determine the correlation between different collagen fiber features and SWE parameters in the 94 cases of breast lesions. A value of P < 0.05 was considered statistically significant.
Overall, 94 breast tumors were surgically removed and histopathologically approved, including 30 benign and 64 malignant lesions (Fig. 1). Benign lesions included fibroadenomas (n = 17), sclerosing adenoses (n = 5), inflammatory lesions (n = 3, including 2 acute suppurative mastitis and 1 granulomatous lobular mastitis), intraductal papillomas (n = 4), and benign phyllodes tumor (n = 1). The malignant lesions included invasive carcinomas (n = 59, including 45 invasive ductal carcinomas, 6 invasive lobular carcinomas, 4 invasive ductal and lobular carcinomas and 4 mucinous carcinomas), and ductal carcinomas in situ (n = 5). Furthermore, the presence of axillary lymph nodes was determined based on postoperative pathology, and malignant breast tumors were divided into LNM + and LNM- groups. Amongst the 64 malignant tumors, 25 (39%) were positive for axillary lymphatic metastases, and 39 (61%) were negative (Table 1).
Elastic parameters, collagen fiber features, and COL1A1 expression
The Ultrasonic elastography, Sirius red staining, and COL1A1 immunohistochemical staining of the included 94 samples are shown in Fig. 2. The ROC curves of Emax, Emean, and Eratio elasticity for the diagnosis of the malignant and benign groups are shown in Fig. 3. Furthermore, the optimal cut-off values yielding the maximum sensitivity and specificity were greater than 58.70 kPa, 52.50 kPa, and 3.05. The areas under the ROC curves were 0.853 (0.779–0.928), 0.823 (0.737–0.910), and 0.882 (0.813–0.951), respectively.
The best cut-off values of SWE for diagnosing lymphatic metastasis were as follows: 107.5 kPa for Emax, 85.15 kPa for Emean, and 3.90 for Eratio. Eratio showed the best diagnostic performance according to the area under the curve (AUC), with a 95% confidence interval (CI) of 0.845 (0.754–0.937). A pairwise comparison of SWE parameters of the benign, LNM-, and LNM + groups showed that the difference was statistically significant (Fig. 4).
In this study, the collagen fiber shape and arrangement features in the breast lesions were classified into three categories. For all 94 lesions, one-way ANOVA showed that Emax, Emean, and Eratio differed between categories 0, 1, and 2 (all P < 0.05; Table 2). This shows that Emax, Emean, and Eratio increase when the collagen fiber class increases from 0 to 2. Furthermore, Kendall's tau-b correlation coefficient calculation showed that the correlation coefficient between Emax and collagen category was 0.318, that between Emean and collagen category was 0.261, and that between Eratio and collagen category was 0.349 (all P ≤ 0.001; Table 3).
COL1A1 expression levels were used for differential diagnosis of benign and malignant breast lesions. The mean expression level of COL1A1 in benign breast lesions was 0.105 ± 0.323; in malignant breast lesions, the mean expression level was 0.159 ± 0.047 (P = 0.012). Furthermore, an average optical density > 0.145 was the cut-off value, the effective rate was 0.808 (0.723–0.894), and the sensitivity and specificity were 66% and 90%, respectively. When the expression level of COL1A1 was used to diagnose LNM- and LNM + malignant breast lesions, the average optical density > 0.150 was taken as the cut-off value, the effectiveness was 0.80 (0.69–0.91), and the sensitivity and specificity were 96% and 59%, respectively.
A pairwise comparison of COL1A1 expression levels in the benign, LNM-, and LNM + groups showed that the difference was statistically significant (Fig. 4). Pearson correlation analysis showed a positive correlation between Emax, Emean, Eratio of breast lesions, and the expression level of COL1A1 (r = 0.406, 0.362, and 0.425, respectively). (P < 0.001; Fig. 5).
SWE is a new US diagnostic technology. It uses shear waves induced by acoustic radiation pulses to visualize and quantify tissue stiffness in a real-time, reliable, and reproducible manner . Tissue stiffness is an important parameter for diagnosing potential malignancies or other diseases . Our study showed that compared with benign breast lesions, malignant lesions had higher Emax, Emean, and Eratio, which was consistent with the results of several previous studies [20,21,22].
Furthermore, considering that preoperative knowledge of axillary lymph node status is crucial for the use of neoadjuvant chemotherapy and appropriate surgical treatment, it is necessary to develop imaging tools for LNM diagnosis [23, 24]. As expected from current trials, the main role of US on the breast in the near future may be to exclude the presence of axillary LNM [25, 26]. Our study showed that the best cut-off values of the SWE parameters for diagnosing lymphatic metastasis were calculated as follows: 107.5 kPa for Emax, 85.15 kPa for Emean, and 3.9 for Eratio.
According to ROC analysis, Eratio exhibited the best diagnostic performance with an AUC of 0.85 (0.75–0.94). Xin Wen et al.  reported that lesion stiffness was a predictor of axillary lymph node metastasis in breast cancer. Additionally, the optimal cut-off values of SWE parameters for predicting LNM were calculated as follows: 111.05 kPa for Emax, 79.80 kPa for Emean, and 6.89 for EmeanR. According to ROC analysis, Emax exhibited the best diagnostic performance with an AUC of 0.82 (0.76–0.87). This difference may have been caused by the following three factors: fewer cases were included in our study, the setting of the instrument and how Eratio was obtained (there is no relevant information in Wen's report), and operator-dependent changes.
Similarly, the present study demonstrated that the shape and arrangement of ECM collagen fibers are divided into three categories, which positively correlate with the Emax, Emean, and Eratio of breast lesion stiffness assessed by SWE. According to the results of all 94 lesions, Emax, Emean, and Eratio showed an increasing tendency from category 0 to category 2 (all P < 0.05).
In normal breast tissue, collagen fibers showed wavy or sinusoidal waves [13, 14]. Thus, this type of collagen fiber was assigned to category 0. Brabrand et al. reported that collagen in the peritumoral region was arranged in a more parallel alignment and aligned or linearized collagen correlated with increased stiffness of the breast tissue or lesions [10, 28]. We classified this type as class 1. Type 2 is characterized by bundles of straightened collagen fibers perpendicular to the tumor boundary, which can change into a honeycomb shape and show high hardness . These results indicate that breast lesion stiffness increases when ECM collagen changes into taut, straightened and arranged parallel fibers. Additionally, a more complicated structure, collagen fibers connected in a honeycomb arrangement, indicated the highest stiffness. This indicates that ECM collagen fiber shape and arrangement may also contribute to the stiffness variance in malignant breast lesions.
The survival of patients with breast cancer is correlated with stromal biology, including the reorganization of the ECM to promote cancer invasion and migration . Uddin et al. suggested that the expression of the stromal gene COL1A1 is associated with the progression of breast cancer and recurrence-free survival in breast cancer patients . In our cohort, 30 tissue samples in the benign group and 64 tissue samples in the malignant group were consecutively included for COL1A1 immunohistochemical staining. Pearson correlation analysis showed a positive correlation between Emax, Emean, Eratio of breast lesions, and the expression level of COL1A1 (r = 0.406, 0.362, and 0.425, respectively; P < 0.001).
Furthermore, the average expression level of COL1A1 in benign breast lesions was significantly lower than that in malignant breast lesions. COL1A1 expression can be used to predict lymph node metastasis. An average optical density > 0.150 was taken as the cut-off value, the effectiveness was 0.80 (0.69–0.91), and the sensitivity and specificity were 96% and 59%, respectively. Moreover, many studies [31, 32] have also proved the correlation between collagen fiber content and SWE in breast lesions and proved that changes in collagen fiber content lead to changes in tumor tissue hardness. Li et al.  studied 148 triple-negative breast cancers and found that high COL1A1 expression in triple-negative breast cancer was an independent prognostic factor.
In addition, protein–protein interaction network analysis confirmed that COL1A1 is a prognostic matrix gene in breast cancer, and its expression is related to the progression of breast cancer . These results suggest that COL1A1 may regulate tumor metastasis through some signaling pathways, and we have begun to investigate the possible signaling mechanisms using in vivo and in vitro methods.
This study had four limitations that need to be addressed. First, this study did not include the correlation between stiffness and other components of the ECM, such as elastic fibers, laminin, and fibronectin; however, collagen is the most important structural protein in the ECM, and its correlation with the elastic modulus may be the most important. Second, the number of included cases was small. Third, limited to the type of US application, examinations were performed by only one operator (no evaluation of inter-observer variability).
The collagen fiber features and expression levels of COL1A1 were positively correlated with the elastic parameters of breast lesions. The expression of COL1A1 may be helpful in diagnosing benign and malignant breast lesions and predicting axillary lymph node metastasis.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Collagen type I alpha 1
Collagen type I alpha 2
Shear wave elastography
Invasive breast cancers
Lymph node metastasis
- LNM -:
Negative axillary lymph node metastasis
- LNM + :
Positive axillary lymph node metastasis
Integrated optical density
Receiver operating characteristic
Area under the curve
Li T, Mello-Thoms C, Brennan PC. Descriptive epidemiology of breast cancer in China: incidence, mortality, survival and prevalence. Breast Cancer Res Treat. 2016;159(3):395–406.
Berg WA. Reducing unnecessary biopsy and follow-up of benign cystic breast lesions. Radiology. 2020;295(1):52–3.
Suvannarerg V, et al. Diagnostic performance of qualitative and quantitative shear wave elastography in differentiating malignant from benign breast masses, and association with the histological prognostic factors. Quant Imaging Med Surg. 2019;9(3):386–98.
Wang ZL, et al. Study of quantitative elastography with supersonic shear imaging in the diagnosis of breast tumours. Radiol Med. 2013;118(4):583–90.
Zhang MK, et al. TGF-β1: is it related to the stiffness of breast lesions and can it predict axillary lymph node metastasis? Ann Transl Med. 2021;9(10):870.
Seewaldt V. ECM stiffness paves the way for tumor cells. Nat Med. 2014;20(4):332–3.
Evans A, et al. Does shear wave ultrasound independently predict axillary lymph node metastasis in women with invasive breast cancer? Breast Cancer Res Treat. 2014;143(1):153–7.
Chamming’s F, et al. Shear wave elastography of tumour growth in a human breast cancer model with pathological correlation. Eur Radiol. 2013;23(8):2079–86.
Lee SH, et al. Shear-wave elastographic features of breast cancers: comparison with mechanical elasticity and histopathologic characteristics. Invest Radiol. 2014;49(3):147–55.
McConnell JC, et al. Increased peri-ductal collagen micro-organization may contribute to raised mammographic density. Breast Cancer Res. 2016;18(1):5.
van ’t Veer LJ, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415(6871):530–6.
Prockop DJ. Mutations that alter the primary structure of type I collagen. The perils of a system for generating large structures by the principle of nucleated growth. J Biol Chem. 1990;265(26):15349–52.
Burke K, Tang P, Brown E. Second harmonic generation reveals matrix alterations during breast tumor progression. J Biomed Opt. 2013;18(3):31106.
Falzon G, Pearson S, Murison R. Analysis of collagen fibre shape changes in breast cancer. Phys Med Biol. 2008;53(23):6641–52.
Provenzano PP, et al. Collagen density promotes mammary tumor initiation and progression. BMC Med. 2008;6:11.
Shi XQ, et al. Correlation between elastic parameters and collagen fibre features in breast lesions. Clin Radiol. 2018;73(6):595.e1-595.e7.
Anvari A, Halpern EF, Samir AE. Essentials of statistical methods for assessing reliability and agreement in quantitative imaging. Acad Radiol. 2018;25(3):391–6.
Bamber J, et al. EFSUMB guidelines and recommendations on the clinical use of ultrasound elastography. Part 1: Basic principles and technology. Ultraschall Med. 2013;34(2):169–84.
Riching KM, et al. 3D collagen alignment limits protrusions to enhance breast cancer cell persistence. Biophys J. 2014;107(11):2546–58.
Qu S, et al. Osterix promotes the migration and angiogenesis of breast cancer by upregulation of S100A4 expression. J Cell Mol Med. 2019;23(2):1116–27.
Link T, et al. Clinical relevance of circulating MACC1 and S100A4 transcripts for ovarian cancer. Mol Oncol. 2019;13(5):1268–79.
Jia W, et al. Shear wave elastography and pulsed doppler for breast lesions: Similar diagnostic performance and positively correlated stiffness and blood flow resistance. Eur J Radiol. 2022;147: 110149.
Meng L, et al. Development of a prediction model based on LASSO regression to evaluate the risk of non-sentinel lymph node metastasis in Chinese breast cancer patients with 1–2 positive sentinel lymph nodes. Sci Rep. 2021;11(1):19972.
He Z, et al. Identification of risk factors associated with axillary lymph node metastasis for sentinel lymph node-positive breast cancer patients. J Oncol. 2020;2020:8884337.
van Roozendaal LM, et al. Clinically node negative breast cancer patients undergoing breast conserving therapy, sentinel lymph node procedure versus follow-up: a Dutch randomized controlled multicentre trial (BOOG 2013–08). BMC Cancer. 2017;17(1):459.
Gentilini O, Veronesi U. Abandoning sentinel lymph node biopsy in early breast cancer? A new trial in progress at the European Institute of oncology of Milan (SOUND: sentinel node vs observation after axillary UltraSouND). Breast. 2012;21(5):678–81.
Wen X, et al. Quantitative shear wave elastography in primary invasive breast cancers, based on collagen-S100A4 pathology, indicates axillary lymph node metastasis. Quant Imaging Med Surg. 2020;10(3):624–33.
Brabrand A, et al. Alterations in collagen fibre patterns in breast cancer. A premise for tumour invasiveness? Apmis. 2015;123(1):1–8.
Conklin MW, Keely PJ. Why the stroma matters in breast cancer: insights into breast cancer patient outcomes through the examination of stromal biomarkers. Cell Adh Migr. 2012;6(3):249–60.
Uddin MN, Wang X. Identification of key tumor stroma-associated transcriptional signatures correlated with survival prognosis and tumor progression in breast cancer. Breast Cancer. 2022;29(3):541–61.
Liu G, et al. Shear wave elasticity of breast lesions: would it be correlated with the extracellular matrix components? Gland Surg. 2019;8(4):399–406.
Wang ZL, et al. Relationship between elasticity and collagen fiber content in breast disease: a preliminary report. Ultrasonics. 2015;57:44–9.
Li FF, Fan JJ, Ma BL. [Effect of high expression of type I collagen α1 chain protein on the prognosis of triple negative breast cancer]. Zhonghua Zhong Liu Za Zhi. 2020;42(2):122–6.
This study was supported by grants from the National Natural Science Foundation of China (No.82071925).
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Chinese PLA General Hospital (No. S2021-683–01), and informed consent was taken from all individual participants.
Consent for publication
All authors disclosed no relevant relationships.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jiang, Y., Wang, B., Li, J.K. et al. Collagen fiber features and COL1A1: are they associated with elastic parameters in breast lesions, and can COL1A1 predict axillary lymph node metastasis?. BMC Cancer 22, 1004 (2022). https://doi.org/10.1186/s12885-022-10092-7
- Breast lesions
- Lymph node metastasis