Quantitative stain-free imaging and digital profiling of collagen structure reveal diverse survival of triple negative breast cancer patients
Breast Cancer Research volume 22, Article number: 42 (2020)
Stromal and collagen biology has a significant impact on tumorigenesis and metastasis. Collagen is a major structural extracellular matrix component in breast cancer, but its role in cancer progression is the subject of historical debate. Collagen may represent a protective layer that prevents cancer cell migration, while increased stromal collagen has been demonstrated to facilitate breast cancer metastasis.
Stromal remodeling is characterized by collagen fiber restructuring and realignment in stromal and tumoral areas. The patients in our study were diagnosed with triple-negative breast cancer in Singapore General Hospital from 2003 to 2015. We designed novel image processing and quantification pipelines to profile collagen structures using numerical imaging parameters. Our solution differentiated the collagen into two distinct modes: aggregated thick collagen (ATC) and dispersed thin collagen (DTC).
Extracted parameters were significantly associated with bigger tumor size and DCIS association. Of numerical parameters, ATC collagen fiber density (CFD) and DTC collagen fiber length (CFL) were of significant prognostic value for disease-free survival and overall survival for the TNBC patient cohort. Using these two parameters, we built a predictive model to stratify the patients into four groups.
Our study provides a novel insight for the quantitation of collagen in the tumor microenvironment and will help predict clinical outcomes for TNBC patients. The identified collagen parameters, ATC CFD and DTC CFL, represent a new direction for clinical prognosis and precision medicine. We also compared our result with benign samples and DICS samples to get novel insight about the TNBC heterogeneity. The improved understanding of collagen compartment of TNBC may provide insights into novel targets for better patient stratification and treatment.
Breast cancer is the second leading cause of cancer-associated mortality in women  and represents a serious global health concern. WHO has estimated that globally, > 1.7 million new breast cancer cases and 522,000 breast cancer-associated deaths occurred in 2012 . Breast cancer can be classified as either non-invasive ductal carcinoma in situ (DCIS) or invasive ductal carcinoma (IDC). Invasive tumors can be further categorized into different molecular subtypes, including luminal A, luminal B, Her2+, and basal-like/triple-negative breast cancers (TNBCs), based on expression profiling or protein surrogates [2,3,4]. TNBC is defined as ER(−), PR(−), and HER2(−) and is typically more heterogeneous and aggressive than other subtypes, resulting in a relatively poor prognosis. Scant data and lack of biomarkers limit us to build robust models that can stratify patients and predict clinical outcomes [5,6,7,8,9]. Differences between ethnic groups and uneven reporting serve to further complicate matters. For example, while the metabolic profile of African-American TBNC patients has been reported , Asian women are subjected to different risk factors due to lifestyle and genetic differences, which can modulate treatment responses and impact disease outcomes [11,12,13,14,15,16]. Despite the association between breast stromal collagen density and invasive breast cancer, which is second only to deleterious germline BRCA1 and BRCA2 mutations [17, 18], stromal biology remains poorly delineated compared with other cancer cell compartments and immune cell populations. Although some previous breast cancer studies incorporate Asian women, large-scale, in-depth, and specific collagen structural data from non-Caucasian women, especially those with TNBC, are lacking. In this work, our group sought to investigate the effect of collagen structure on TNBC outcomes in a cohort of Asian women.
The PD-L1/PD-1 treatment is promising, while the overall response rate is not ideal including breast cancer. Many factors may impact the patient outcome. Mesenchymal cells, immune cells, extracellular matrix (ECM) components, lymphatics, and vasculature are all present in breast cancer stroma. Immune cell infiltration, disruption of the basal membrane and myoepithelial cell layer, and remodeling of stromal collagen are among the earliest key events in the development of invasive breast cancer [19, 20]. Recent studies have also demonstrated that stromal reaction influences the efficacy of particular treatments, including immunotherapy [21, 22]. Furthermore, the importance of stromal and collagen biology in breast tumor progression is highlighted by differential immune, angiogenic, and fibroblastic responses, and an array of stromal genes may be used to predict the clinical outcome [23, 24]. Collagen fibers represent the major structural ECM component in breast tumors, and increases in stromal collagen fibers have been demonstrated to facilitate breast tumor formation, invasion, and metastasis. Collagen is secreted by cancer-associated fibroblasts (CAFs), which are involved in tumor stromal activation, and may lead to tumor progression through multiple mechanisms, including neoangiogenesis, tumor cell proliferation, and invasion [25,26,27,28,29]. CAFs also affect tumor progression by reprogramming the tumor microenvironment at both the metabolic and immune levels and by promoting adaptive resistance to chemotherapy . In this work, we only focus on the quantification of the collagen remodeling and explore its impact on the TNBC patient survival.
Stromal collagen remodeling is characterized by collagen realignment in the stromal compartment, and multiple recent studies have demonstrated that collagen structure, profiles, and patterns within the tumor stromal microenvironment have diagnostic and prognostic value [31,32,33,34,35,36,37,38,39,40]. However, the role of collagen remodeling in tumor formation and metastasis is complex, and while correlations between breast tissue density and tumor formation have been reported , the relationships between collagen density, structure, and cancer cell migration and invasion remain to be fully understood. Collagen in breast tissue was historically thought to form a physical barrier that prevents tumor cell migration, with collagen degradation and deposition being prerequisites for this process [31, 33, 34, 38, 39]. Furthermore, three tumor-associated collagen signatures (TACs) have been reported , and radial alignment of collagen fibers around the boundaries of tumors is associated with tumor cell invasion. As a result, this topic remains the subject of debate within the field, and no conclusive statement has been reached.
This issue may, at least partially, be due to historical methodological limitations. The remodeling and quantification of collagen have traditionally been studied using biochemical staining techniques, such as Picrosirius Red (PSR) staining or Masson’s trichrome (MT) staining. However, this analysis is highly dependent on the staining protocol and color deconvolution pipeline in digital image analysis. Due to the limitations presented by these imaging techniques, only the total amount of collagen present in tissue samples is typically assessed, with the quantitative structure of collagen fibers remaining beyond their scope. However, the development of quantitative stain-free collagen imaging techniques, such as second harmonic generation (SHG), permits the quantification of collagen structure at a finer level of detail. In this study, we used two-photon excitation (TPE) and SHG to scan breast cancer tissue microarrays (TMAs). SHG is a multiphoton, laser-based, quantitative non-linear optical imaging technique used to identify fibrillary collagen in fixed tissues. Due to its physical principles, it is highly sensitive to changes in collagen fibril and fiber structure and also to the remodeling of connective tissue.
We developed a fully automatic digital collagen profiling platform based on TPE/SHG imaging techniques to quantify collagen using four parameter categories: intensity/area, textural, structural, and fiber distribution features. Numerical image features were extracted from each image. The extracted features were selected using feature selection algorithms, and associations between these features and the existing clinicopathological parameters were investigated across the whole patient cohort. Bioinformatics models were designed for the purpose of classification (diagnosis) and prediction (prognosis). Our results provide a computational solution to classify collagen into two distinct modes based on fiber SHG signal intensity, texture, and morphology: aggregated thick collagen (ATC) and dispersed thin collagen (DTC). Several imaging features were strongly correlated with clinicopathological characteristics, and ATC collagen fiber density (CFD) and DTC collagen fiber length (CFL) were revealed to be of prognostic value based on our patient cohort and their clinical outcomes. Separation of ATC and DTC provides a novel understanding of collagen remodeling during cancer progression, and our results may help to resolve the debate over whether collagen has a role in inhibiting or promoting patient survival. All extracted parameters are listed in Supp. Table 01, and they are quantified in ATC and DTC region separately.
Materials and methods
Formalin-fixed paraffin-embedded (FFPE) human breast cancer tissues were obtained from the Department of Anatomical Pathology, Division of Pathology, Singapore General Hospital (SGH). The present study utilized archival specimens from a cohort of 388 TNBC breast cancer patients diagnosed at the Division of Pathology of SGH from 2003 to 2015. The Centralized Institutional Review Board (IRB) of SingHealth provided ethical approval for the use of patient material in this study (CIRB ref: 2018/2910-2013/664/F). Clinicopathological parameters including age, ethnicity, tumor size, histologic grade, subtype, associated ductal carcinoma in situ, lymphovascular invasion, axillary lymph node status, tumor borders, and growth patterns were reviewed.
The tissues were fixed and processed in Surgipath 10% neutral buffered formalin (Leica Biosystems, Wetzlar, Germany) for 6–48 h in accordance with the American Society of Clinical Oncology/College of American Pathologists (ASCO/CAP) Clinical Practice Guidelines [41, 42].
Hematoxylin and eosin (H&E)-stained tumor sample slides were reviewed. Three tumor cores (1 mm diameter) were punched from corresponding donor FFPE blocks. The three cores are reviewed by clinical pathologists, and the selected representative core is used to construct the TMA assays using an MTA-1 Manual Tissue Arrayer (Beecher Instruments, Inc., Sun Prairie, WI, USA), as previously described .
Image acquisition with SHG
The constructed TMA assays were deparaffinized and then scanned using the Genesis 200™ system (HistoIndex Pte. Ltd., Singapore). The laser power of the system was set at “High,” and binning was set at two frames to reduce the noise level. The brightness of each pixel represents the collagen density in the fiber.
A total of eight TMAs was scanned, and corresponding images were acquired, containing data from all 388 TNBC patients. These images were then processed according to the computational pipeline jointly developed by the pathologists and computer scientists. The details of this pipeline are described later. Each image of a given TMA sample was named according to its pathological identification, linking the given image to clinicopathological parameters.
A representative TNBC TMA image is presented in Fig. 1a, including the TPE/SHG image and the H&E-stained image in Fig. 1b. After a TPE/SHG image has been taken from the unstained slides, the same slide can be stained using H&E. The collagen compartments can then be merged with the H&E image, as shown in Fig. 1c. The merging of collagen with H&E provides information concerning collagen remodeling at different regions in the tissue. Although PSR or MT staining is often used to quantify collagen in tissue samples, these are not truly quantitative imaging solutions, as the red color in PSR and the blue color in MT staining do not necessarily indicate collagen fiber structure. SHG, however, is a laser-based quantitative imaging technique. Once the system is well calibrated, the physical principles result in a numerical reading of collagen density, as shown in Fig. 1e–h. Fiber density is a critical indication of the collagen packing structure. A representative stronger and weaker fiber signal is shown in Fig. 1e, g; their respective SHG signal intensity is shown in Fig. 1f, h. Measurements of SHG intensity are taken starting from the blue dot and ending at the red dot along the path of the yellow dotted line. A brighter collagen fiber has much higher SHG intensity in Fig. 1f, and a weaker collagen fiber has lower SHG intensity. Besides the collagen fiber density, named as CFD in this work, it is possible to extract other collagen structural information from the images.
The processing of the TMA image is as follows: A sample mask was first defined to cover the area occupied by the TMA tissue. Within this given region, a two-mode Gaussian mixture model (GMM) was applied to determine whether a given pixel signaled positive for collagen. Based on the GMM model, an optimized thresholding value was determined to binarize the image, creating a collagen binary image. A similar scheme was applied to the red channel to create a tissue binary image. The collagen area ratio (CAR) was then measured based on the area of the collagen binary image and sample mask. The collagen fiber density (CFD) is the ratio between the summation of the total collagen density signal within the collagen binary image and the area of the collagen binary image. Skeletonization was applied to the collagen binary image in order to identify the central line of the collagen network, as shown by the red solid line in Fig. 1i. The junction points were then identified, as shown by the yellow dots in Fig. 1i. The number of junction points is an indication of the complexity of the collagen network. The ratio between the number of junction points and the total length of the skeleton provides the collagen reticular index (CRI). Other structural features, such as collagen fiber length (CFL) and collagen fiber thickness (CFT), were also measured accordingly. A list of key parameters and their definitions are provided in Supp. Table 02. All parameters in this work are measured at the ATC, DTC, and ALL tissue regions, which is the combination of both ATC and DTC.
In order to separate ATC and DTC in the tissue, an intensity and morphology filtering pipeline was developed. Based on the intensity and morphology (the degree of collagen aggregation), it is possible to differentiate the sample mask into two distinct modes, as presented in Fig. 1d. We applied an intensity morphological opening operation with a diameter bigger than the 3 times of collagen fiber thickness to identify the aggregated collagen area, i.e., ATC area which is red in Fig. 1d. The rest of the sample mask area, i.e., blue area in Fig. 1d, will be defined as the DTC area. We also tested the diameter for the opening operation, and it is robust to identify isolated fibers. Intensity/area, textural, structural, and fiber distribution features were measured separately for both ATC and DTC. Once the features were extracted from the images, their relationships with clinicopathological parameters were assessed. Associations with low p values are presented in corresponding figures. The selected features were then used to build a prediction model for patient survival.
Statistical analysis pipeline
Follow-up data were obtained from medical records for the entire cohort. Disease-free survival (DFS) and overall survival (OS) were defined as time from the point of diagnosis to recurrence or to death/the date of the last follow-up, respectively. Statistical analysis was performed using the R statistical language (version 3.3.1). The relationships between clinicopathological parameters and collagen profiles were tested using χ2 and Fisher’s exact tests. Survival analyses were conducted using univariate Cox regression models based on binary thresholded imaging parameters. The final threshold was selected from a range of thresholds determined using the percentiles of the imaging parameters in 10% increments that gave the best separation between the two groups.
Numerical collagen structural parameters are associated with clinicopathological parameters
Numerical collagen features were extracted from digital images as described in the “Materials and methods” section. Detailed formulations and definitions of the key parameters are provided in Supp. Table 02. The association between numerical imaging features and clinicopathological parameters was studied. All defined features are measured separately in the ATC area (red area in Fig. 1d), DTC area (blue area in Fig. 1d), and ALL area (the combined red and blue area in Fig. 1d). For example, ATC CFD, DTC CFD, and ALL CFD were measured as different regions for CFD.
Tumor size is a key clinicopathological parameter used to evaluate the likelihood of cancer progression. We screened the possible associations between our numerical collagen features and tumor size. As shown in Fig. 2a, there are significant differences in both ATC CAR and DTC CAR between patients with tumors of ≤ 20 mm and > 20 mm. (Note that the significant difference in ALL CAR in Fig. 2a is a weighted measurement between ATC and DTC.) The ATC area ratio, which represents the ratio between the red area and the whole TMA area as shown in Fig. 1d, was also higher of tumor > 20 mm, but no significance (data not presented). Notably, there was a significant difference in ATC CFD, DTC CFD, and ALL CFD between patients with tumors of ≤ 20 mm and > 20 mm in Fig. 2b. This is a strong indication that larger tumors tended to have a stronger impact in microenvironment remodeling around tumor nests. Furthermore, there was a significant difference in both ATC CFL and DTC CFL in Fig. 2c. The latter points towards an interplay between tumor size and the surrounding stromal environment remodeling.
The correlation between the quantitative collagen features and the degree of DCIS association for those TNBC patients were of particular interest. The DCIS association is a pathological feature to describe if the TNBC patient is associated with DICS subtype. First, patients with associated DCIS had a greater ATC area ratio compared with those with no associated DCIS. ATC CAR and DTC CAR were also significantly higher in patients with a positive DCIS association (Fig. 3a). The data indicated that ALL CAR was significantly different in these two sub groups. ATC CFD was higher in patient samples with associated DCIS compared with those with no associations (Fig. 3b). This indicated that ATC collagen may have formed a stronger protective barrier surrounding the DCIS. There was a strong significant difference in DTC CFD (Fig. 3b) and DTC CFL (Fig. 3c), as well as in ATC CFT and DTC CFT (Fig. 3d) between patients with associated DCIS and those with no associated DCIS. This may imply that TNBCs which harbored associated DCIS potentially have different fiber packing structures, which may cause differences in fiber brightness, length, and thickness. The aggregation of the collagen is relevant with DCIS association.
Certainly, not all clinicopathological parameters are relevant with the extracted quantitative features, which describe the properties of collagen. For example, tumor grade is a critical pathological parameter associated with tumor behavior and is evaluated based on nuclear pleomorphism, mitotic activity, and tubule formation. However, we did not find a strong association between our numerical collagen parameters and tumor grades in our patient cohort. One potential rational is that although tumor grading is a key feature used in prognostication, it is not directly linked with collagen remodeling. The CFL and CFD results are presented in Supp. Figure 1A and Supp. Fig. 1B, respectively.
In order to evaluate the prognostic value of the collagen parameters, we performed univariate Cox regression univariate analyses. Most of the parameters listed in Supp. Table 02 are able to stratify the patients in their own way. For example, we found out that “higher” ATC CFD was associated with improved DFS and OS in TNBC patients (DFS, p = 0.0163; OS, p = 0.0085). All cutoff values of the parameters were determined as described in the “Materials and methods” section. The Kaplan-Meier survival plots for ATC CFD are shown in Supp. Figure 2A-B. We further explore if any combinations of two quantitative collagen parameters are able to stratify the patients in more detail. We found ATC CFD combined with DTC CFL was of significant prognostic value for patient DFS and OS. (The survival curves of DTC CFL are shown in Supp Fig. 2C-D (DFS, p = 0.3; OS, p = 0.24)). Predictive prognostic models are essential for clinical practice to stratify the patients in the healthcare system. Figure 4a, b indicates that when the 177 patients were divided into four distinct groups based on ATC CFD and DTC CFL, they experienced significant differences in survival (DFS, p = 0.0046; OS, p = 0.0062). This further demonstrates the heterogeneity of the TNBC subtype. The number of patients in each group is also moderately well balanced (Fig. 4). “Brighter” ATC CFD and “longer” DTC CFL indicate the best prognosis, as shown by the blue curve, while “weaker” ATC CFD and “longer” DTC CFL indicate the worst, as shown by the red curve. The other two conditions, represented by the black and green curves, fall between the aforementioned extremes and result in intermediate patient survival. The impacts of different significant parameters on these four groups of patients are summarized in Supp. Table 3.
Our results may help resolve the historical dilemma concerning the role of collagen in breast cancer. Proposed mechanisms underlying the effect of ATC CFD and DTC CFL on survival are presented in Fig. 5a–d. The differences may be due to the collagen bonds aggregating in different ways: for example, ATC CFD(+) represents a protective structure formed by the strongly aggregated collagen. Depending on how strong this first protective structure is, we hypothesize that DTC CFL(+) may represent a potential “highway” that facilitates either immune cell migration or cancer cell migration. As presented in Fig. 5a, ATC CFD(+) may manifest as a strong layer of collagen fibers around the tumor nest. It would potentially be difficult for malignant cells to penetrate this layer and escape from their original site. Even when cells escape, in the DTC CFL(+) condition, long collagen fibers facilitate migration of immune cells to respond and deal with malignant cells. This combination, therefore, provides the best survival outcome for the patients as illustrated in Fig. 5a, with cancer cell migration potentially inhibited. Figure 4d, however, shows the opposite extreme. In the ATC CFD(−) condition, the collagen may be unable to form a strong protective barrier, meaning that cancer cells are able to escape easily. Once this occurs, longer collagen fibers in the DTC area (DTC CFL(+)) may facilitate their migration more efficiently and therefore overcome the capacity of immune cells to respond to the invasion. Figure 4b, c shows the two intermediate conditions. In this proposed model, patient survival is primarily impacted by the following two factors: (i) how easily the malignant cells are able to escape from their original site (This is relevant with the ATC CFD and ATC area ratio.) and (ii) how easily they can migrate once this occurs, which is associated with the DTC CFL. We show here how the stromal collagen can have a bipartite role in cancer progression promoting either malignant cell migration or immune cell response depending on the state of the ATC CFD and DTC CFL.
The 2D distribution of patient ATC CFD and DTC CFL is presented in Fig. 6a. The color-coded spots indicate the patient groups in Fig. 4. However, the mechanism underlying collagen aggregation remains unclear. The ATC area ratio is defined as the ratio between the red area in Fig. 1d and the total tissue area (the combination of red and blue areas in Fig. 1d). This means the percentage of the given tissue area that is occupied by ATC. The ATC area ratio of patients with the best survival, i.e., ATC CFD(+) and DTC CFL(+) patients presented by the blue survival curve in Fig. 3a, b, was significantly higher than in the other three conditions as shown in Fig. 6b. Even if we only compared with ATC CFD(+) and DTC CFL(−) patients, represented by the green bar in Fig. 6b, the ATC area ratio of patients with the best survival was almost threefold higher. Based on our patient grouping, the ATC CFD of the patients with the best survival is also significantly higher than the group with the worse survival. The following two conditions were therefore associated with optimal survival: (i) an increased amount of aggregated collagen in the tissue and (ii) brighter collagen fibers, which provide indirect evidence of different collagen fiber microstructures. These results provide evidence suggesting that collagen aggregation itself plays a key role in patient survival. Figure 6d shows the DTC CFL for the four patient groups. Differences in DTC may cause different degrees of cell migration, including both immune cells and malignant cells, in the tissue, as well as subsequent invasion, with longer collagen fibers facilitating cancer cell migration to remote areas. Thus, collagen fiber length may also have impacted patient survival.
Conclusion and discussion
Increased computational power and the integration of computational solutions with artificial intelligence (AI) and machine learning (ML) have resulted in substantial improvements to clinical diagnosis and prognosis. This has had significant implications for the field of histopathology too, where traditional assessments typically involve qualitative characterization on a numerical basis. The resultant discrete numbers are often not sufficient to build robust mathematical models for clinical diagnosis or prognosis, particularly when generated from manual assessments with certain subjective inconsistency. Thus, scanning and digitization of glass slides represent the future for AI-assisted pathology. Imaging techniques are also changing, shifting from basic H&E staining to other techniques, including immunohistochemistry, immunofluorescence, and stain-free imaging techniques, such as SHG. In this work, we applied TPE/SHG to investigate the effect of collagen structures on the prognosis of TNBC patients, which is an aggressive but heterogeneous subtype. SHG imaging quantified the collagen component, including collagen types I and III. Our team developed an automatic image analysis solution to extract different collagen features from the SHG images, including collagen structural information.
The literature has reported that while collagen can potentially form a protective barrier and prevent cancer cells from escaping their original site, collagen fibers have also been found to serve as a “highway” that facilitates cancer cell migration to remote locations, impacting patient survival [44,45,46,47]. Previously, the mechanisms underlying these contradictory reports were unclear . This historical lack of understanding was due to the paucity of methods to differentiate between collagen structures. As demonstrated by a number of previous studies, it is an oversimplification to characterize collagen remodeling as a simple increase or decrease in total ECM collagen content . It is more important to understand the structural differences in various regions under diverse conditions. For example, the orientation of collagen fibers, such as TACs, is a key parameter . Fibrils are packed to form collagen fibers. The organization and distribution of these fibers within the tissue are also important, for example, whether the collagen forms any aggregations. One unique parameter generated by SHG in the present study was CFD, which provides a strong indication of collagen fiber strength and fibril packing structure, as the laser power and system calibration generate a correlation between pixel brightness and collagen strength. In our study, we showed some extracted parameters were strongly associated with certain clinicopathological features, including tumor size and DCIS associations. However, tumor grade was not relevant with the investigated collagen features as we expected. This is because grade assessment is mainly based on nuclei morphological and packing features. It is not linked to collagen characteristics. Regardless, the results of our study suggest the importance of assessing collagen structure, as it appears to represent a potential key prognostic parameter. According to our data, tumor grade and tumor size are also poor predictors of prognosis, as shown in Supp Fig. 3.
Furthermore, our results suggest that the collagen compartment has two distinct modes, named ATC and DTC, and differences between these modes underpin the different roles of collagen. Collagen remodeling can either prevent or promote cancer cell migration, and while the mechanism remains to be fully understood, we demonstrated that ATC CFD and DTC CFL together have strong prognostic value. Higher ATC CDF indicates the presence of stronger collagen fibers surrounding the tumor nest, through which it is difficult for cancer cells to escape (Fig. 5a). Therefore, this condition results in a more favorable prognosis. DTC CFL also affects patient survival in a complex way. If a strong protective layer is present to prevent the invasion of malignant cells, the longer collagen fibers facilitate immune cells to promptly respond to combat those malignant cells. On the other hand, if the first layer of protection is weak leading to accumulation of a large number of malignant cells in the stromal area, longer collagen fibers may provide a more effective “highway” which assists cancer cell migration. Once the cancer cell attaches to the long collagen fibers, it is easier for them to “grasp” it and subsequently migrate long distances. Figure 5 schematizes the four possible conditions, and survival curves for these groups are presented in Fig. 4.
Furthermore, due to limitations of previous techniques, the quantification of collagen structure and investigation of the relationships between numerical collagen structure parameters and patient survival were not previously feasible. We demonstrated that patients with the best survival (ATC CFD(+) and DTC CFL(+)) had increased aggregated collagen area (Fig. 6b) and increased collagen brightness, representing density (Fig. 6c), compared with other groups. However, increases in collagen brightness may have been caused by collagen aggregation. The presence of more collagen in a given small area may also have resulted in a stronger signal. However, in Supp. Fig. 4A, high DTC CFD was observed in patients with the best survival. This suggests that in the DTC area, the collagen fiber was also brighter. Through combining the above data, we conclude that the collagen fiber packing structure differed in such a way that it caused high CFD in both the ATC and DTC regions. Potentially, the collagen fiber packing structure is also associated with the collagen aggregation patterns. In both the DTC and ATC regions, the CFT average value was the same. This further supports the observation that the CFD is not due to a simple clustering of nearby fibers but a real increase of intensity of the fibers (the denser the fibrils within the fibers are, the higher the SHG intensity is, directly linking the intensity of the fiber with the density of fibrils that compose it). As summarized in Supp. Table 3, the ATC is the main actor in the difference between the best and worse survival groups whereas the DTC is the key to highlight differences among all 4 groups and notably the intermediate groups. ATC area ratio is a key parameter to differentiate the best survival patients with the rest of condition, while other parameters contribute to further classify patients on other aspects. Although the results of this study demonstrate that the aggregation of collagen and collagen brightness (represented by CFD and indicating different collagen packing structures) and the DTC CFL are two key collagen-associated parameters that may impact patient survival, the generic mechanism that controls the differences described in our data remains unclear and requires further investigation.
TNBC patients represent a heterogeneous group . Lymph node status is the only prognostic feature for TNBC patients, also supported by our existing data, as shown in Supp. Fig. 5. The use of this prognostic factor is dependent on the detection of cancer cells in the lymph nodes. The percentage and number of patients with positive lymph node status in the four groups is presented in Supp. Fig. 6A and B. In total, 55% of ATC CFD (−) and DTC CFL(+) patients were lymph node positive, last row in Supp. Fig. 6A; the prognostic model presented in the study shows some correlation with the lymph node status, i.e., 55% of worse survival patients were lymph node positive whereas only 67% of the best survival patients were lymph node negative. Theoretically, it would be possible to build a model based on the collagen parameters investigated in this study to predict patient survival even before metastasis.
The collagen structure parameters extracted in the present study were closely associated with certain pathological parameters, and the mechanism presented in Fig. 5 warrants further investigation. We demonstrated that collagen structure differs in three aspects: (i) the degree of aggregation, (ii) collagen fiber packing structure, and (iii) the length of the collagen in the DTC regions. Although this evidence suggests a potential answer to the historical debate on the role of collagen in breast cancer survival, there may or may not have been causality. Other cellular factors may affect collagen structure, for example, migrating cell behavior, including both immune cell and the migratory capability of cancer cells followed by epithelial-mesenchymal transition (EMT). For example, NK cells have antifibrotic properties which decrease with progression of fibrosis . Tumor necrosis factor α (TNF-α) , transforming growth factor β (TGF-β) [51,52,53], IL-11 [54, 55], and OSM [56, 57] are some of the known cytokines involved in the fibrogenesis pathways. However, if any of these factors will impact the abovementioned three aspects and then affect the TNBC patient survival still need further investigation.
Finally, we compared our results with the benign samples and DCIS samples. The data of ATC area ratio, ATC CFD, and DTC CFL are presented in Supp. Fig. 7, 8, 9, respectively. The original data points are presented in panel A, and the mean differences, taking benign as a reference, are presented in panel B. As we mentioned in the “Results” section, ATC CFD(+) and DTC CFL(+) have the best survival. According to the ATC area ratio, ATC CFD, and DTC CFL of these three features, the original data of this best survival population has very similar pattern as patients who were diagnosed as DICS and the mean differences are also very limited. Although other conditions may have one or two parameters closer to DICS patient, they have worse survival as we discussed in Fig. 5.
Although the finding in this work is promising for better stratification of TBNC patients in healthcare system, we must emphasize that it is still preliminary and not yet replicated using whole slide imaging data from an independent cohort. Well-designed rigid validation must be conducted to further confirm its performance before the application in clinical work. The model presented in this work is based on ATC CFD and DTC CFL as highlighted by the blue and green colors in Supp. Table 3; however, other parameters may also contribute to the patient stratification indicating the underline governing mechanism of collagen role in the TNBC patients might be more complicated. It requires further elucidation and is worthy of additional large-scale study. On the other hand, ATC area ratio, highlighted by red in Supp. Table 3, is able to nicely differentiate the best survival patients with other three groups. In practical, this more simplified model can possibly provide some immediate insight in clinic. Tumor progression is an outcome of mutations in multiple different genes, and the relationships between collagen characteristics and other critical tumor compartments, such as the immune response, EMT, adipose, and angiogenesis, should be investigated in a more systematic study.
Availability of data and materials
The images and numerical data of this work will be available upon request.
Ferlay J, Soerjomataram I, Dikshit R, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136:E359–86.
Thike AA, Cheok PY, Jara-Lazaro AR, et al. Triple-negative breast cancer: clinicopathological characteristics and relationship with basal-like breast cancer. Mod Pathol. 2009;23:123–33.
Iqbal J, Ginsburg O, Rochon PA, et al. DIfferences in breast cancer stage at diagnosis and cancer-specific survival by race and ethnicity in the United States. JAMA. 2015;313:165–73.
Foulkes WD, Smith IE, Reis-Filho JS. Triple-negative breast cancer. N Engl J Med. 2010;363:1938–48.
Dent R, Trudeau M, Pritchard KI, et al. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res. 2007;13:4429–34.
Kassam F, Enright K, Dent R, et al. Survival outcomes for patients with metastatic triple-negative breast cancer: implications for clinical practice and trial design. Clin Breast Cancer. 2009;9:29–33.
Thike AA, Yong-Zheng Chong L, Cheok PY, et al. Loss of androgen receptor expression predicts early recurrence in triple-negative and basal-like breast cancer. Mod Pathol. 2014;27:352–60.
Cheng CL, Thike AA, Tan SY, et al. Expression of FGFR1 is an independent prognostic factor in triple-negative breast cancer. Breast Cancer Res Treat. 2015;151:99–111.
Matsumoto H, Koo SL, Dent R, et al. Role of inflammatory infiltrates in triple negative breast cancer. J Clin Pathol. 2015.
Kanaan YM, Sampey BP, Beyene D, Esnakula AK, Naab TJ, Ricks-Santi LJ, Dasi S, Day A, Blackman KW, Frederick W, Copeland RL Sr, Gabrielson E, Dewitty RL Jr. Metabolic profile of triple-negative breast cancer in African-American women reveals potential biomarkers of aggressive disease. Cancer Genomics Proteomics. 2014;11(6):279–94.
Sung H, Rosenberg PS, Chen WQ, et al. Female breast cancer incidence among Asian and Western populations: more similar than expected. J Natl Cancer Inst. 2015;107(7). https://doi.org/10.1093/jnci/djv107. Print 2015 Jul.
Leong SPL, Shen Z-Z, Liu T-J, et al. Is breast cancer the same disease in Asian and Western countries? World J Surg. 2010;34:2308–24.
Bhoo-Pathy N, Yip C-H, Hartman M, et al. Breast cancer research in Asia: adopt or adapt Western knowledge? Eur J Cancer. 2013;49:703–9.
de Bruin MA, Kwong A, Goldstein BA, et al. Breast cancer risk factors differ between Asian and white women with BRCA1/2 mutations. Familial Cancer. 2012;11:429–39.
Kim H, Choi DH. Distribution of BRCA1 and BRCA2 mutations in Asian patients with breast cancer. J Breast Cancer. 2013;16:357–65.
Wong ESY, Shekar S, Met-Domestici M, et al. Inherited breast cancer predisposition in Asians: multigene panel testing outcomes from Singapore. Npj Genomic Medicine. 2016;1:15003.
Hwang ES, Miglioretti DL, Ballard-Barbash R, et al. Association between breast density and subsequent breast cancer following treatment for ductal carcinoma in situ. Cancer Epidemiol Biomark Prev. 2007;16:2587–93.
Boyd NF, Dite GS, Stone J, et al. Heritability of mammographic density, a risk factor for breast cancer. N Engl J Med. 2002;347:886–94.
Khamis ZI, Sahab ZJ, Sang QX. Active roles of tumour stroma in breast cancer metastasis. Int J Breast Cancer. 2012;2012:574025.
Kalluri R, Zeisberg M. Fibroblasts in cancer. Nat Rev Cancer. 2006;6:392–401.
Turley SJ, Cremasco V, Astarita JL. Immunological hallmarks of stromal cells in the tumour microenvironment. Nat Rev Immunol. 2015;15:669–82.
Jiang H, Hegde S, Knolhoff BL, et al. Targeting focal adhesion kinase renders pancreatic cancers responsive to checkpoint immunotherapy. Nat Med. 2016;22:851–60.
Finak G, Bertos N, Pepin F, et al. Stromal gene expression predicts clinical outcome in breast cancer. Nat Med. 2008;14:518–27.
Lisanti MP, Sotgia F, Pestell RG, et al. Stromal glycolysis and MCT4 are hallmarks of DCIS progression to invasive breast cancer. Cell Cycle. 2013;12:2935–6.
Bhowmick NA, Neilson EG, Moses HL. Stromal fibroblasts in cancer initiation and progression. Nature. 2004;432:332–7.
Orimo A, Weinberg RA. Stromal fibroblasts in cancer: a novel tumour-promoting cell type. Cell Cycle. 2006;5:1597–601.
Mercier I, Casimiro MC, Wang C, et al. Human breast cancer-associated fibroblasts (CAFs) show caveolin-1 downregulation and RB tumour suppressor functional inactivation: implications for the response to hormonal therapy. Cancer Biol Ther. 2008;7:1212–25.
Giannoni E, Bianchini F, Masieri L, et al. Reciprocal activation of prostate cancer cells and cancer-associated fibroblasts stimulates epithelial-mesenchymal transition and cancer stemness. Cancer Res. 2010;70:6945–56.
Erez N, Truitt M, Olson P, et al. Cancer-associated fibroblasts are activated in incipient neoplasia to orchestrate tumour-promoting inflammation in an NF-kappaB-dependent manner. Cancer Cell. 2010;17:135–47.
Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16:582–98.
Acerbi I, Cassereau L, Dean I, et al. Human breast cancer invasion and aggression correlates with ECM stiffening and immune cell infiltration. Integr Biol. 2015;7:1120–34.
Tan WJ, Yan J, Xu S, et al. Second harmonic generation microscopy is a novel technique for differential diagnosis of breast fibroepithelial lesions. J Clin Pathol. 2015;68:1033–5.
Burke K, Brown E. The use of second harmonic generation to image the extracellular matrix during tumour progression. IntraVital. 2014;3:e984509.
Tilbury K, Campagnola PJ. Applications of second-harmonic generation imaging microscopy in ovarian and breast cancer. Perspectives Med Chem. 2015;7:21–32.
Robinson BK, Cortes E, Rice AJ, et al. Quantitative analysis of 3D extracellular matrix remodelling by pancreatic stellate cells. Biology Open. 2016;5:875–82.
Tilbury K, Campagnola PJ. Applications of second-harmonic generation imaging microscopy in ovarian and breast cancer. Perspect Med Chem. 2015;7:21–32.
Keikhosravi A, Bredfeldt JS, Sagar AK, et al. Second-harmonic generation imaging of cancer. Methods Cell Biol. 2014;123:531–46.
Birk JW, Tadros M, Moezardalan K, et al. Second harmonic generation imaging distinguishes both high-grade dysplasia and cancer from normal colonic mucosa. Dig Dis Sci. 2014;59:1529–34.
Nadiarnykh O, LaComb RB, Brewer MA, et al. Alterations of the extracellular matrix in ovarian cancer studied by Second Harmonic Generation imaging microscopy. BMC Cancer. 2010;10:94.
Bredfeldt JS, Liu Y, Conklin MW, et al. Automated quantification of aligned collagen for human breast carcinoma prognosis. J Pathol Inform. 2014;5:28.
Allison KH, Hammond MEH, Dowsett M, McKernin SE, Carey LA, Fitzgibbons PL, Hayes DF, Lakhani SR, Chavez-MacGregor M, Perlmutter J, Perou CM, Regan MM, Rimm DL, Symmans WF, Torlakovic EE, Varella L, Viale G, Weisberg TF, McShane LM, Wolff AC. Estrogen and Progesterone Receptor Testing in Breast Cancer: American Society of Clinical Oncology/College of American Pathologists Guideline Update. Arch Pathol Lab Med. 2020;144(5):545-563. https://doi.org/10.5858/arpa.2019-0904-SA. Epub 2020 Jan 13.
Wolff AC, Hammond MEH, Allison KH, Harvey BE, Mangu PB, Bartlett, JMS, Bilous M, Ellis IO, Fitzgibbons P, Hanna W, Jenkins RB, Press MF, Spears PA, Vance GH, Viale G, McShane LM, Dowsett M. Human Epidermal Growth Factor Receptor 2 Testing in Breast Cancer: American Society of Clinical Oncology/College of American Pathologists Clinical Practice Guideline Focused Update. J Clin Oncol. 2018;36(20):2105-2122. https://doi.org/10.1200/JCO.2018.77.8738. Epub, 2018 May 30.
Thike AA, Yong-Zheng Chong L, Cheok PY, Li HH, Wai-Cheong Yip G, Huat Bay B, Tse GM, Iqbal J, Tan PH. Loss of androgen receptor expression predicts early recurrence in triple-negative and basal-like breast cancer. Mod Pathol. 2014;27(3):352-60. https://doi.org/10.1038/modpathol.2013.145. Epub 2013 Aug 9.
Conklin MW, Keely PJ. Why the stroma matters in breast cancer: insights into breast cancer patient outcomes through the examination of stromal biomarkers. Cell Adhes Migr. 2012;6(3):249–60. https://doi.org/10.4161/cam.20567 Epub 2012 May 1.
Provenzano PP, Eliceiri KW, Campbell JM, Inman DR, White JG, Keely PJ. Collagen reorganization at the tumour-stromal interface facilitates local invasion. BMC Med. 2006;4(1):38.
Conklin MW, Eickhoff JC, Riching KM, Pehlke CA, Eliceiri KW, Provenzano PP, Friedl A, Keely PJ. Aligned collagen is a prognostic signature for survival in human breast carcinoma. Am J Pathol. 2011;178(3):1221–32. https://doi.org/10.1016/j.ajpath.2010.11.076.
Takai K, Le A, Weaver VM, Werb Z. Targeting the cancer-associated fibroblasts as a treatment in triple-negative breast cancer. Oncotarget. 2016;7(50):82889–901. https://doi.org/10.18632/oncotarget.12658.
Fang M, Yuan J, Peng C, Li Y. Collagen as a double-edged sword in tumour progression. Tumour Biol. 2014;35(4):2871–82. https://doi.org/10.1007/s13277-013-1511-7 Epub 2013 Dec 15.
Jeong WI, Park O, Gao B. Abrogation of the antifibrotic effects of natural killer cells/interferon-gamma contributes to alcohol acceleration of liver fibrosis. Gastroenterology. 2008;134:248–58. https://doi.org/10.1053/j.gastro.2007.09.034.
Yang YM, Seki E. TNFα in liver fibrosis. Curr Pathobiol Rep. 2015;3(4):253–61 Epub 2015 Sep 30.
Isabel F and Daniel CD, Transforming growth factor-β-induced cell plasticity in liver fibrosis and hepatocarcinogenesis. Front. Oncol., 2018 https://doi.org/10.3389/fonc.2018.00357.
Fabregat I, Moreno-Càceres J, Sánchez A, Dooley S, Dewidar B, Giannelli G, Ten Dijke IT-LIVER Consortium. TGF-β signalling and liver disease. FEBS J 2016;283(12):2219–2232. doi: https://doi.org/10.1111/febs.13665. Epub 2016 Feb 22.
Neuschwander-Tetri BA, Loomba R, Sanyal AJ, Lavine JE, Van Natta ML, Abdelmalek MF, Chalasani N, Dasarathy S, Diehl AM, Hameed B, Kowdley KV, McCullough A, Terrault N, Clark JM, Tonascia J, Brunt EM, Kleiner DE, Doo E. NASH clinical research network. Farnesoid X nuclear receptor ligand obeticholic acid for non-cirrhotic, non-alcoholic steatohepatitis (FLINT): a multicentre, randomised, placebo-controlled trial. Lancet. 2015;385(9972):956–65.
Widjaja AA, Singh BK, Adami E, Viswanathan S, Dong J, D’Agostino GA, Ng B, Lim WW, Tan J, Paleja BS, Tripathi M, Lim SY, Shekeran SG, Chothani SP, Rabes A, Sombetzki M, Bruinstroop E, Min LP, Sinha RA, Albani S, Yen PM, Schafer S, Cook SA. Inhibiting interleukin 11 signaling reduces hepatocyte death and liver fibrosis, inflammation, and steatosis in mouse models of nonalcoholic steatohepatitis. Gastroenterology. 2019. doi: https://doi.org/10.1053/j.gastro.2019.05.002. [Epub ahead of print].
Anissa A. Widjaja, Brijesh K. Singh, Eleonora Adami, Sivakumar Viswanathan, Giuseppe A. D’Agostino, Jinrui Dong, Benjamin Ng, Jessie Tan, Bhairav S. Paleja, Madhulika Tripathi, Sze Yun Lim, Sonia P. Chothani, Wei Wen Lim, Anne Rabes, Martina Sombetzki, Eveline Bruinstroop, Rohit A. Sinha, Salvatore Albani, Paul M. Yen, Sebastian Schafer, Stuart A. Cook. IL-11 neutralising therapies target hepatic stellate cell-induced liver inflammation and fibrosis in NASH. Gastroenterology doi: https://doi.org/10.1053/j.gastro.2019.05.002.
Ehab A. Ayaub, Anisha Dubey, Jewel Imani, Fernando Botelho, Martin R. J. Kolb, Carl D. Richards & Kjetil Ask. Overexpression of OSM and IL-6 impacts the polarization of pro-fibrotic macrophages and the development of bleomycin-induced lung fibrosis, Scientific Reportsvolume 7, Article number: 13281 (2017).
Michitaka M, Shinya T,Naoko M, Eiko S, Hitoshi O, Atsushi M, Minoru T. Oncostatin M causes liver fibrosis by regulating cooperation between hepatic stellate cells and macrophages in mice. Hepatology. 2018;67(1):296-312. https://doi.org/10.1002/hep.29421. Epub 2017 Nov 15.
This research was jointly funded by the A*STAR Biomedical Research Council, National Medical Research Council Stratified Medicine Programme Office (SMPO201302) awarded to Dr. Puay Hoon Tan, and A*STAR Biomedical Research Council (BMRC IAF 111170) awarded to Dr. Yu Weimiao. Dr. Jabed Iqbal is a recipient of the Transition Award from the Singapore National Medical Research Council (NMRC/TA/0041/2015). The authors would like to thank HistoIndex Pte. Ltd. for complimentary scanning the TNBC slides without obligation.
Ethics approval and consent of participate
The ethics approvals of this study are granted by SingHealth Centralized Institutional Review Board (CIRB) with ref: 2018/2910-2013/664/F.
Consent for publication
All authors of this work agreed to publish with Breast Cancer Research once accepted.
The authors declare no potential competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Collagen length and density in comparison with tumour size and grades. Supplementary Figure 2. Kaplan Meier survival curves for Two individual parameters of prognostic value. Supplementary Figure 3. Kaplan Meier survival curves for tumour size and grade. Supplementary Figure 4. Some other key aspects on the collagen structural differences in the four given groups. Supplementary Figure 5. Kaplan Meier Survival curves for the Lymph node status. Supplementary Figure 6. Lymph node status(+/-) distributions in the four groups of patients. Supplementary Figure 7. Comparison of the ATC Area Ratio between the four groups of patients, Benign samples and the DCIS patient cohort. Supplementary Figure 8. Comparison of the ATC CFD between the four groups of patients, Benign samples and the DCIS patient cohort. Supplementary Figure 9. Comparison of the DTC CFL between the four groups of patients, Benign samples and the DCIS patient cohort. Supplementary Table 1. Comprehensive Parameters List measured in SHG/TPE images for TNBC patient cohort. Supplementary Table 2. Key Parameters List and their definition for TNBC patient cohort. Supplementary Table 3. Pairwise P-values for the set of features measured in 4 groups.
About this article
Cite this article
Gole, L., Yeong, J., Lim, J.C.T. et al. Quantitative stain-free imaging and digital profiling of collagen structure reveal diverse survival of triple negative breast cancer patients. Breast Cancer Res 22, 42 (2020). https://doi.org/10.1186/s13058-020-01282-x
- Triple-negative breast cancers
- Collagen profile
- Quantitative imaging
- Second harmonic generation microscopy