Claudin-Low Breast Cancer; Clinical & Pathological Characteristics

Abstract

Claudin-low breast cancer is a molecular type of breast cancer originally identified by gene expression profiling and reportedly associated with poor survival. Claudin-low tumors have been recognised to preferentially display a triple-negative phenotype, however only a minority of triple-negative breast cancers are claudin-low. We sought to identify an immunohistochemical profile for claudin-low tumors that could facilitate their identification in formalin fixed paraffin embedded tumor material. First, an in silico collection of ~1600 human breast cancer expression profiles was assembled and all claudin-low tumors identified. Second, genes differentially expressed between claudin-low tumors and all other molecular subtypes of breast cancer were identified. Third, a number of these top differentially expressed genes were tested using immunohistochemistry for expression in a diverse panel of breast cancer cell lines to determine their specificity for claudin-low tumors. Finally, the immunohistochemical panel found to be most characteristic of claudin-low tumors was examined in a cohort of 942 formalin fixed paraffin embedded human breast cancers with >10 years clinical follow-up to evaluate the clinico-pathologic and survival characteristics of this tumor subtype. Using this approach we determined that claudin-low breast cancer is typically negative for ER, PR, HER2, claudin 3, claudin 4, claudin 7 and E-cadherin. Claudin-low tumors identified with this immunohistochemical panel, were associated with young age of onset, higher tumor grade, larger tumor size, extensive lymphocytic infiltrate and a circumscribed tumor margin. Patients with claudin-low tumors had a worse overall survival when compared to patients with luminal A type breast cancer. Interestingly, claudin-low tumors were associated with a low local recurrence rate following breast conserving therapy. In conclusion, a limited panel of antibodies can facilitate the identification of claudin-low tumors. Furthermore, claudin-low tumors identified in this manner display similar clinical, pathologic and survival characteristics to claudin-low tumors identified from fresh frozen tumor material using gene expression profiling.

About the Authors:

Kay Dias

Affiliations Department of Pathology and Molecular Medicine, McMaster University, Hamilton, Ontario, Canada, Department of Oncology, McMaster University, Hamilton, Ontario, Canada

Anna Dvorkin-Gheva

Affiliation: Department of Biochemistry and Biomedical Sciences, McMaster University, Hamilton, Ontario, Canada

Robin M. Hallett

Affiliation: Department of Biochemistry and Biomedical Sciences, McMaster University, Hamilton, Ontario, Canada

Ying Wu

Affiliation: Department of Oncology, McMaster University, Hamilton, Ontario, Canada

John Hassell

Affiliations Department of Pathology and Molecular Medicine, McMaster University, Hamilton, Ontario, Canada, Department of Biochemistry and Biomedical Sciences, McMaster University, Hamilton, Ontario, Canada

Gregory R. Pond

Affiliation: Department of Oncology, McMaster University, Hamilton, Ontario, Canada

Mark Levine

Affiliation: Department of Oncology, McMaster University, Hamilton, Ontario, Canada

Tim Whelan

Affiliation: Department of Oncology, McMaster University, Hamilton, Ontario, Canada

Anita L. Bane

* E-mail: [email protected]

Affiliations Department of Pathology and Molecular Medicine, McMaster University, Hamilton, Ontario, Canada, Department of Oncology, McMaster University, Hamilton, Ontario, Canada

Introduction

In 2007, while conducting comparative gene expression analysis between transgenic mouse models of breast cancer and human breast cancer data sets, Herschkowitz et al discovered a novel molecular subtype of breast cancer which they named ‘claudin-low’ (CL) [1]. This subtype was characterized by the low expression of genes involved in tight junctions and epithelial cell-cell adhesion, including claudins 3, 4 and 7, occludin and E-cadherin. In addition, the human CL tumors showed low expression of luminal epithelial genes and high expression of lymphocyte and endothelial cell markers [1].

Subsequent to this report, a number of groups have further characterized this new tumor subtype and shown that CL tumors account for 7–14% of all invasive breast cancers, are enriched for genes associated with epithelial to mesenchymal transition (EMT), immune cell infiltration, IFNγ activation, mammary stem cells/breast tumor initiating cells and typically demonstrate high levels of genomic instability [2–5]. Pathologic examination of a limited number of tumors fitting this category, have shown a higher than expected prevalence of medullary-like and metaplastic special-type tumors and tumors with a triple negative (TN) phenotype. Clinically, they have been associated with a poor prognosis with some evidence that they may be relatively resistant to conventional chemotherapeutic agents [2, 3].

These studies have all identified the CL subtype by means of gene expression profiling, however this technique requires fresh frozen tumor material, which is not available for the majority of breast cancer patients. We sought to determine an immunohistochemical (IHC) profile that could identify CL tumors in formalin fixed paraffin embedded (FFPE) tumor specimens. Such a profile would enable us and others, to examine the pathologic and clinical significance of the CL tumors in larger cohorts of archived FFPE tumor specimens providing a more comprehensive analysis of the tumor subtype.

To this end we collated an in silico data base comprising the expression profiles of approximately 1600 individual breast cancers. Tumors comprising this database were classified into the known molecular subtypes; luminal A, luminal B, HER-2 enriched, basal-like, normal-like, molecular apocrine and CL using previously published classifiers [2, 6–8]. Genes differentially expressed between the CL and all other molecular subtypes of breast cancer were identified. The protein products of some of the top differentially expressed genes were examined using IHC for their ability to identify the CL subtype in a diverse panel of breast cancer cell lines. Finally, the panel of IHC markers that in combination optimally discriminated between CL tumor cell lines and cell lines of other molecular subtypes was examined in a large tissue microarray (TMA) of primary invasive human breast cancers. The tumors identified in this TMA cohort as CL using the surrogate IHC profile described were compared with tumors of other molecular subtypes for clinical-pathologic characteristics of known prognostic importance, disease free survival (DFS), overall survival (OS) and local recurrence rates (LRR) (S1 Fig).

Materials and Methods

In silico data collection

In the course of our study we analyzed the gene expression profiles in silico of 7 external datasets, obtained using Affymetrix HG-U133A GeneChip arrays. These profiles were deposited in the Gene Expression Omibus (GEO) (accession numbers of the datasets are: GSE3494, GSE1456, GSE7390, GSE2034, GSE6532, GSE17705 and GSE25066) and comprise a total of 2,027 samples (S1 Table). Redundant samples were removed as previously described reducing the number of unique samples to 1,695 [8]. All samples used for our study were normalized with frozen Robust Multi-array Analysis (fRMA) [9], technical variation was removed using the ComBat and DWD (Distance-Weighted Discrimination) methods [10, 11]. After combining all datasets Spearman correlation coefficients for pair-wise comparisons of samples using 68 house-keeping probe sets were computed, and only samples exhibiting a correlation higher than 0.95 with at least half of the dataset were selected for further classification. The latter filtering method yielded a dataset comprising 1,593 human breast tumor sample transcript profiles.

In silico molecular subtype assignment

The 1,593 tumors were assigned to one of 7 molecular subtypes (luminal A, luminal B, HER-2-enriched, basal-like, normal-like, molecular apocrine or CL) using 710 genes obtained from previously published gene classifiers [6–8, 12]. In brief the standardized centroid was computed for each subtype by taking the average expression of each gene across the subtype and dividing it by the standard deviation of expression of that gene across that subtype. Spearman rank correlation coefficient was computed for each sample relative to each of the 7 reference centroids, and the subtype was assigned based on the highest correlation coefficient. For the assignment we used a coefficient cut-off of 0.3; therefore 1,196 samples were classified into the 7 established molecular subtypes as previously described [8]. This classification yielded 80 (6.69%) samples defined as CL.

Identification of genes differentially expressed between CL tumors and all other molecular subtypes

To identify gene patterns unique to CL tumors we used the “limma” package (Bioconductor; [13]) to compare the expression profiles of samples assigned to the CL subtype to those of the all other subtypes (6 pair-wise comparisons in total). To this end the moderated F-statistic was used, followed by Benjamini-Yekutieli adjustment for multiple testing [14]. The 710 genes belonging to the molecular subtype classifier were used for this analysis and only genes differentially expressed with at least a 2 fold change were examined further.

Breast tumor cell lines representative of the molecular subtypes of breast cancer

A total of 9 breast cancer cell lines known to replicate the luminal, basal and CL subtypes of primary human breast tumor samples were grown using the recommended culture conditions [15], these included 5 luminal (MCF7, ZR751, SKBR3, BT474, MDA-MB-361), 2 basal-like (BT20, HCC 1954) and 2 CL cell lines (BT549, MDA-MB-231) [15]. Cells were fixed and paraffin embedded as detailed in the supplementary information (S1 File). The paraffin embedded cell lines were stained immunohistochemically for ER, PR, HER2, CK5, EGFR, E-cadherin, claudin 3, claudin 4, claudin 7 and CD24 using methods as listed in S2 Table.

Human FFPE breast tumor material

942 T1 or T2, node negative breast cancers treated with breast conserving therapy, which had been accrued as part of the Accelerated Hypofractionated Whole Breast Irradiation (AHWBI) trial (see S1 File) were available for analysis [16, 17]. This cohort had 10 years of clinical follow-up available including, LR, DFS and OS.

A single hematoxylin and eosin (H&E) stained section, representative of each invasive carcinoma was reviewed by the study pathologist (ALB). Tumors were assessed for tumor type, grade, lympho-vascular space invasion (LVI), extensive lymphocytic infiltrate and margin circumscription. Tumors were classified according to the World Health Organization (WHO) histologic classification of breast tumors [18] and graded using the Nottingham grading system [19]. The lymphocytic tumor infiltrate was graded on a four point scale; none, minimal, moderate and extensive (see S1 File). The study pathologist was blinded to the patient outcome during the review process.

The invasive tumor component of each H&E stained section was encircled with permanent ink for TMA construction. Three 0.6mm cores of tissue were taken from the paraffin tumor block and used for TMA construction (Pathology Device, Sun Praire, WI) as previously described [16].

Four μm sections were cut from all TMAs and immunohistochemical staining for ER, PR, HER2, CK5, EGFR, Ki67, claudin 3, claudin 4, claudin 7, E-cadherin, CD24, CD44, ALDH1 and CD8 was performed using methods as listed in S2 Table. Microwave antigen retrieval was carried out in a Micromed T/T Mega Microwave Processing Lab Station (ESBE Scientific, Markham, Ontario, Canada). Sections were developed with diaminobenzidinetetrahydrochloride (DAB) and counterstained in Mayer’s hematoxylin.

Each of the immunohistochemical TMA and tumor cell line stained sections was scored using Allred’s scoring method [20], which adds scores for the intensity of staining (absent: 0, weak: 1, moderate: 2, and strong: 3) to the percentage of cells stained (none: 0, <1%: 1, 1–10%: 2, 11–33%: 3, 34–66%: 4 and 67–100%: 5) to yield a ‘raw’ score of 0 or 2–8. Previously validated cut-offs for ER and PR were used (0, 2 = negative, 3–8 = positive) [21, 22]. Strong complete membranous staining was assessed for HER2 and the cut-off of ≥ 6 was used to indicate positivity [23]. For CK5, EGFR, CD24, CD44 and ALDH1 a score of ≥ 4, was considered positive. For claudin 3, 4, 7 and E-cadherin a score of ≤4 was considered ‘low’ expression. For Ki67 a minimum of 100 tumor nuclei were counted per core and the tumor was considered Ki67 ‘low’ if the percentage of positively stained nuclei was <14% and Ki67 ‘high’ is the percentage of positively stained nuclei was ≥14% [12, 24]. The raw score data from the TMAs were reformatted using a TMA deconvoluter software program into a format suitable for statistical analysis [25]. The highest score from each TMA tumor triplicate was entered into the statistical analysis.

Tumors that had an Allred score of 4 or 5 for HER2 were considered equivocal or indeterminate for HER2 overexpression and fluorescent in situ hybridization (FISH) was performed on representative tumor sections using the HER2 DNA probe kit (Path-Vysion, Vysis) as previously described [16]. A HER2 to centromere 17 ratio of ≥2 was considered to indicate amplification in accordance with guidelines [23].

Tumors were classified as luminal A if they expressed ER or PR and were negative for HER2 and were Ki67 ‘low’; luminal B if they expressed ER or PR and were either HER2 positive or were Ki67 ‘high’; HER2 enriched if they did not express ER or PR but were positive for HER2; basal-like if they did not express ER, PR or HER2 (triple negative) but expressed CK5 and/or EGFR and CL if they did not express ER or PR or HER2 (TN) and had low expression of at least two of the following markers E-cadherin, claudin 3, claudin 4 and claudin 7 [16, 26]. Tumors were considered unclassified for molecular subtype when results of one or more IHC markers were unavailable due to loss of invasive tumor on sequential TMA slides.

Ethics statement

The AHWBI FFPE samples were obtained with research ethics board (REB) approval. We did not pursue individual patient consent for tumor samples and this was not required by our REB process for a number of reasons, including. 1. The trial was performed many years ago from 1993 to 1996. As such we recognized that many patients had died and many others were likely to have changed residence. 2. We believed that contacting patients' families would be difficult and likely upsetting given the limited study we were conducting. 3. The analysis we planned is limited to IHC testing and while this is linked to the patients' original data base from the trial all the data was grouped and made anonymous.

Statistical analysis

Summary statistics were used to describe the patient cohort and outcomes. The Kaplan-Meier method was used to estimate time-to-event outcomes. Comparison between different subtypes was performed using the log-rank tests, χ2 test, Cochran-Armitage test for trend or Kruskal-Wallis test as appropriate. All statistical tests were two-sided and statistical significance was defined as a p-value of 0.05 or less. Statistical analyses were performed in SAS version 9.0 (SAS Institute, Cary, NC) and figures were plotted using R version 3.2.2 (www.r-project.org).

Results

Identification of CL subtype using gene expression profiling

We compiled gene expression profiles from 7 independent datasets for which clinical follow-up was available. Together these data sets represent 1,593 non-redundant tumors. Using established subtype classifiers [2, 6–8] 1,196 (75.1%) out of the 1,593 samples were classified into one of the seven molecular subtypes; luminal A (n = 340, 28.43%), luminal B (n = 216, 18.06%), HER2-enriched (n = 81, 6.77%), basal (n = 280, 23.41%), normal-like (n = 179, 14.97%), molecular apocrine (n = 20, 1.68%) and CL (n = 79, 8.4%). The OS and DFS of the molecular subtypes are presented in S2 Fig.

Identification of genes differentially expressed between the CL subtype and all other molecular subtypes of breast cancer

Using a pair wise comparison 60 genes were identified as being differentially expressed by CL tumors relative to all other molecular subtypes of breast cancer (Table 1, Fig 1). Some of the genes we found to be expressed at significantly lower levels in the CL subtype relative to all other subtypes included E-cadherin, claudin 3, claudin 4 and genes associated with luminal epithelial differentiation including CD24, CK8 and CK18 as previously described [1, 2, 5]. As other authors have identified low expression of claudin 7 as being a characteristic feature of CL tumors, we examined the expression of claudin 7 across the breast cancer subtypes in our in silico database. We found that claudin 7 was expressed at a lower level in CL tumors relative to all other subtypes with the exception of basal-like tumors (S3 Fig)

[Figure omitted. See PDF.]

Fig 1. Hierarchical Clustering of 1,593 breast tumor samples using the 60 genes identified as differentially expressed between CL tumors and tumors of all other molecular subtypes

https://doi.org/10.1371/journal.pone.0168669.g001

[Figure omitted. See PDF.]

Table 1. Genes differentially expressed between CL tumors and all other molecular subtypes of breast cancer using gene-expression profiling.

https://doi.org/10.1371/journal.pone.0168669.t001

Genes expressed at significantly higher levels in CL tumors relative to other molecular subtypes included many genes involved in immune response, host defence and apoptosis including ADAMDEC1, BTN3A3, CD3D, COLEC12, CXCL9, LTB, PSMB10 and MAF (Table 1). The preponderance of immune related genes in CL tumors is thought to reflect considerable immune cells infiltrate in these tumors [5, 27].

Genes associated with tumor invasiveness and epithelial to mesenchymal transition (EMT); CTSK and PLAC8 respectively (Table 1), were also identified as being upregulated in CL tumors [28]. An enrichment for gene signatures associated with EMT has previously been demonstrated for CL tumors [4, 29–31].

Evaluation of immunohistochemical markers of CL tumors in breast cancer cell lines

We next sought to capture the CL phenotype identified in silico with a simple IHC based assay. We first compiled a panel of breast cancer cell lines of known molecular subtype including 2 CL cell lines; these were fixed in formalin and embedded in paraffin blocks to ensure that they were processed in a manner analagous to human tumor samples and profiled for ER, PR, HER2, CK5/6 and EGFR as previously described [16]. Taking into consideration the availability of high quality antibodies and our in silico results (Table 1) five antibodies were selected for testing; E-cadherin, claudin 3, claudin 4, claudin 7 and CD24. All of these 5 markers had been shown in silico to be differentially expressed in CL tumors relative to other molecular subtypes of breast cancer.

As illustrated in Table 2, those cell lines known to be CL had absent or low expression of the luminal epithelial markers (ER, PR, HER-2) and the epithelial cell-cell adhesion markers (claudin 3, claudin 4, claudin 7 and E-cadherin) at the cut points described. In contrast, the luminal cell lines were positive for at least one of the luminal epithelial cell markers (ER, PR or HER2) together with at least two of the epithelial cell-cell adhesion markers (claudin 3, claudin 4, claudin 7 and E-cadherin). The basal cell lines were negative for all luminal epithelial cell markers (ER, PR & HER2), positive for myoepithelial cell markers (CK5/6 & EGFR) and positive for at least 3 of the epithelial cell-cell adhesion proteins (claudin 3, claudin 4, claudin7 and E-cadherin). None of the profiled cell lines expressed CD24, which may possibly be an artifact of cell culture.

[Figure omitted. See PDF.]

Table 2. Expression of Immunohistochemical Markers in Breast Cancer Cell Lines.

https://doi.org/10.1371/journal.pone.0168669.t002

From these experiments we concluded that a surrogate IHC panel for the identification of tumors belonging to the CL molecular subtype would be TN (ER-, PR- and HER2-), together with low or absent expression of at least 2 of the 4 epithelial cell-cell adhesion proteins, claudin 3, claudin 4, claudin 7 and E-cadherin.

Clinical-pathologic characteristics of CL tumors using TMAs of human breast cancers

We next sought to examine the clinical-pathologic tumor features and survival characteristics of CL tumors identified using the surrogate IHC panel in a large cohort of invasive breast cancer with long-term clinical outcome data.

942 primary invasive breast tumors arrayed in triplicate in TMAs were examined for the expression of a panel of IHC markers; ER, PR, HER-2, Ki67, CK5/6, EGFR, claudin 3, claudin 4, claudin 7 and E-cadherin to approximate the known molecular subtypes of breast cancer [16, 26]. 776 (82.4%) of the 942 tumors could be classified into one of the five molecular subtypes; luminal A (n = 389, 41.3%), luminal B (n = 234, 24.8%), HER2 enriched (n = 21, 2.2%), basal-like (n = 53, 5.6%) and CL (n = 79, 8.4%) as described in the material and methods (Fig 2). 166 (17.6%) cases could not be classified into one of the molecular subtypes due to unavailable IHC data and were placed into an ‘unclassified’ category for purposes of analysis. A suitable surrogate IHC profile for the molecular apocrine group or the normal-like group is not available and these subtypes were not considered further in this data set.

[Figure omitted. See PDF.]

Fig 2.

Representative H&E stained image of luminal A [A], luminal B [B], HER2 [C], basal [D] and Claudin-low [E] tumors. ER staining in luminal A [F], luminal B [G], HER2 [H], basal [I] and Claudin-low [J] tumors. PR staining in luminal A [K], luminal B [L], HER2 [M], basal [N] and Claudin-low [O] tumors. HER2 staining in luminal A [P], luminal B [Q], HER2 [R], basal [S] and Claudin-low [T] tumors. Ki67 staining in luminal A [U], luminal B [V], HER2 [W], basal [X] and Claudin-low [Y] tumors. EGFR staining in luminal A [Z], luminal B [AA], HER2 [BB], basal [CC] and Claudin-low [DD] tumors. CK5 staining in luminal A [EE], luminal B [FF], HER2 [GG], basal [HH] and Claudin-low [II] tumors. Claudin-3 staining in luminal A [JJ], luminal B [KK], HER2 [LL], basal [MM] and Claudin-low [NN] tumors. Claudin-4 staining in luminal A [OO], luminal B [PP], HER2 [QQ], basal [RR] and Claudin-low [SS] tumors. Claudin-7 staining in luminal A [TT], luminal B [UU], HER2 [VV], basal [WW] and Claudin-low [XX] tumors. E-cadherin staining in luminal A [YY], luminal B [ZZ], HER2 [AAA], basal [BBB] and Claudin-low [CCC] tumors.

https://doi.org/10.1371/journal.pone.0168669.g002

The clinical-pathologic features are listed in Table 3. There were significant differences in median patient age, tumor size, grade, extensive lymphocytic infiltrate and margin circumscription between subtypes. Patients also varied according to adjuvant systemic therapy received.

[Figure omitted. See PDF.]

Table 3. Clinical-pathologic Characteristics of Invasive Breast Cancers According to Molecular Subtypes in the TMA cohort.

https://doi.org/10.1371/journal.pone.0168669.t003

When compared to all other subtypes combined, CL tumors were more likely to occur at a younger age (mean = 52.8 versus 59.0, p<0.001) be high grade (68.1% versus 15.2% were grade III, p<0.001), larger size (32.9% versus 18.5% were ≥2cm, p = 0.005) and to be characterized by an extensive lymphocytic infiltrate (41.9% versus 7.9%, p<0.001) and to have pushing or circumscribed tumor margins (44.3% versus 20.8%, p<0.001) (Table 3, Fig 3).

[Figure omitted. See PDF.]

Fig 3. Disease-Free Survival (DFS) of Patients enrolled in the AHWBI trial by Tumor Subtype

https://doi.org/10.1371/journal.pone.0168669.g003

CL tumors and association with markers of breast cancer stem cells/tumor initiating cells

We examined the association between CL tumors in our TMA cohort and known markers of breast cancer stem cells/tumor initiating cells including ALDH1 and CD44hi/CD24-/low [32–35]. When compared to all other subtypes combined, CL tumors were more likely to express ALDH1 (17.1% versus 7.3%, p = 0.010) and they showed a trend towards an association with the CD44hi/CD24-/low phenotype (31.3% versus 21.7%, p = 0.090) (Table 3, S4 Fig).

Prognosis

In the final cohort of 942 patients, 348 had a disease-related event, 61 had a local recurrence and 204 deaths were observed. With respect to DFS, luminal A and CL cancers had the best prognosis at 10 years (72.2% and 72.5% respectively) whereas basal-like and HER2 enriched had the worst (61.2 and 56.7% respectively), (p = 0.002; Table 3, Fig 3). The majority of the recurrences associated with CL tumors occurred within the first 5 years of diagnosis. Although not statistically significant (p = 0.40) the groups with the worst OS were patients with basal-like and CL tumors (Table 3, Fig 4). With regards to local recurrence (LR), only 1 of 79 (1.3%) patients with a CL tumor had a LR as compared to 21 (5.4%) of luminal A, 15 (6.8%) of luminal B, 6 (28.6%) of the HER2 enriched patients, 4 (7.5%) of the basal-like (p<0.001, Table 3, Fig 5).

[Figure omitted. See PDF.]

Fig 4. Overall Survival (OS) of Patients enrolled in the AHWBI trial by Tumor Subtype

https://doi.org/10.1371/journal.pone.0168669.g004

[Figure omitted. See PDF.]

Fig 5. Local Recurrence (LR) of Patients enrolled in the AHWBI trial by Tumor Subtype

https://doi.org/10.1371/journal.pone.0168669.g005

Discussion

An appropriate IHC surrogate approach for the classification of CL breast tumors has not been heretofore rigorously identified limiting our ability to recognize and study this subtype in further detail. To this end, using an in silico dataset of 1,196 breast tumors each of which could be assigned to one of 7 molecular subtypes; luminal A, luminal B, HER2 enriched, basal-like, normal-like, molecular apocrine and CL we identified that 6.69% (n = 80) of the tumors were CL. These CL tumors differed from other molecular subtypes by their lower expression of epithelial cell-cell adhesion factors, markers of luminal epithelial cell differentiation and by their elevated expression of genes involved in immunity and host defence and tumor cell invasiveness and EMT. These results reaffirm the dominant biological pathways functioning in CL tumors. Breast cancer cell lines reflective of the CL subtype were also found to express low or absent levels of the epithelial cell-cell adhesion proteins examined (E-cadherin, claudin 3, claudin 4 and claudin 7) and markers of breast luminal epithelial cell differentiation (ER, PR and HER-2). This unique IHC profile (TN and cell adhesion protein low/absent) could distinguish the CL subtype in cell line preparations from the other breast cancer molecular subtypes studied.

Furthermore, when we examined a large retrospective collection of 942 primary human breast cancers with clinical, pathologic and outcome data for the surrogate CL IHC profile we identified a tumor group with distinguishing morphologic and outcome characteristics. We observed that the incidence of the CL subtype was 8.4% similar to the 7–14% incidence reported previously [2, 5]. Phenotypic characteristics of this subtype when compared to other subtypes included an association with high tumor grade, large tumor size, an extensive lymphocytic infiltrate and circumscribed/pushing tumor margins. Many of these features are component features of medullary or atypical medullary breast cancer, a subtype of TN breast cancer. An association between these medullary-like features and CL tumors has been previously identified [2]. The extensive lymphocytic infiltrate associated with CL tumors in our primary breast tumor cohort correlates well with the preponderance of immune associated genes upregulated in this subtype as identified by our in silico experiments.

CL tumors have previously been shown to be enriched for genes associated with mammary stem cells/breast cancer tumor initiating cells, from which it has been inferred that this tumor subtype may be enriched for these primitive cells types and even potentially derived from the malignant transformation of a mammary epithelial stem cell [29, 36–40]. To test if our IHC profile accurately captures this aspect of CL tumor biology we examined the expression of known cancer stem cell/breast tumor initiating cell markers; ALDH1 and CD44hi/CD24low/- in our TMA cohort of human tumors [33–35]. The CL subtype was significantly more likely to express the mammary stem cell/breast tumor initiating cell marker ALDH1 (p = 0.01) than non-CL tumors and there was a trend for association between CL tumors and the CD44hi/CD24low/- phenotype (p = 0.09). These results suggest that the surrogate IHC profile used is accurately identifying the CL subtype of tumors.

CL tumors as identified by gene expression profiling have been shown to have an outcome intermediate between that of luminal A and poor prognostic subtypes such as luminal B, basal-like and HER2 enriched subtypes [2, 5]. Using the IHC definition described we show that while not statistically significant (p = 0.04) patients with CL had a worse OS at 10 years (81.6%) compared to patients with luminal A disease (85.8%) (Table 3). The relatively good outcome described for patients with CL tumors in our cohort may be attributable to the presence of an extensive lymphocytic infiltrate in many of these tumors. Increased quantities of tumor infiltrating lymphocytes have been shown by a number of investigators to be associated with good outcome in breast cancer in general and TN tumors in particular [41–44]. The relative good outcome for CL tumors in our cohort could also be ascribed to the low tumor stage of all patients eligible for entry onto this trial (T1 /T2 and N0).

This is the first study to report on the LR rates for CL tumors. The CL subtype had the lowest rate of LR (1.3% at 5 and 10 years) of any molecular subtype studied. Given that all patients in this cohort were treated with breast conserving surgery and whole breast irradiation this finding may suggest that CL tumors are particularly sensitive to radiation. However, given the low number of CL tumors and the overall low LR rate in the study population this result should be considered hypothesis generating and would need to be validated in other data sets.

In summary, we have taken a two step approach to identify a surrogate IHC profile to identify CL tumors in FFPE tumor samples. We tested this profile (TN and low expression of at least two of four epithelial cell-cell adhesion markers; claudin 3, claudin 4, claudin 7 and E-cadherin) in a large cohort of breast tumors with long-term follow up. We have demonstrated that approximately 8% of all invasive breast cancer fall into this unique molecular subtype and the tumor type is characterized by distinguishing morphologic features including high tumor grade, large size and some of the characteristic features of medullary-type cancers. Uniquely we demonstrate that CL tumors have a low incidence of LR following breast conserving therapy (BCT). In addition, CL tumors show an association with known cancer stem cell markers when compared with all other molecular subtypes of breast cancer. While our results are encouraging and suggest that using the IHC panel described, CL tumors can be identified it would be valuable to validate these findings in an independent breast cancer cohort with long-term follow-up,

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. Flow diagram of work.

https://doi.org/10.1371/journal.pone.0168669.s001

(TIF)

S1 Table. Files used for in silico analysis.

https://doi.org/10.1371/journal.pone.0168669.s002

(DOCX)

S1 File. Supplementary information for materials and methods.

https://doi.org/10.1371/journal.pone.0168669.s003

(DOC)

S2 Table. Antibodies and conditions of use.

https://doi.org/10.1371/journal.pone.0168669.s004

(DOCX)

S2 Fig. Overall survival (OS) and disease-free survival (DFS) by molecular subtype of the patients belonging to the in silico cohort of 1,593 breast cancers.

https://doi.org/10.1371/journal.pone.0168669.s005

(DOCX)

S3 Fig. Box plot illustrating the expression of claudin 7 across the different molecular subtypes of breast cancer.

https://doi.org/10.1371/journal.pone.0168669.s006

(DOCX)

S4 Fig. Composite image of luminal A, luminal B, HER2, Basal and Claudin-low tumors stained with H&E and ALDH1, CD44, CD24 and CD8.

https://doi.org/10.1371/journal.pone.0168669.s007

(TIF)

Acknowledgments

We would like to acknowledge the technical contributions of Ying Wu, who performed all the immunohistochemical staining described.

Author Contributions

1. Conceptualization: KD ALB.

2. Data curation: ADG GRP RMH.

3. Formal analysis: GRP ML TW JH ADG.

4. Funding acquisition: ALB JH.

5. Investigation: ALB KD.

6. Methodology: KD GRP YW ALB.

7. Project administration: ALB.

8. Resources: ALB.

9. Software: ADG RMH.

10. Supervision: ALB.

11. Validation: ALB YW.

12. Visualization: ALB.

13. Writing – original draft: KD ALB.

14. Writing – review & editing: KD ALB TW ML JH.

Citation: Dias K, Dvorkin-Gheva A, Hallett RM, Wu Y, Hassell J, Pond GR, et al. (2017) Claudin-Low Breast Cancer; Clinical & Pathological Characteristics. PLoS ONE 12(1): e0168669. https://doi.org/10.1371/journal.pone.0168669

References

1. Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol 2007;8(5):R76. pmid:17493263

2. Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, et al. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res 2010 Sep 2;12(5):R68. pmid:20813035

3. Prat A, Perou CM. Deconstructing the molecular portraits of breast cancer. Mol Oncol 2011 Feb;5(1):5–23. pmid:21147047

4. Hennessy BT, Gonzalez-Angulo AM, Stemke-Hale K, Gilcrease MZ, Krishnamurthy S, Lee JS, et al. Characterization of a naturally occurring breast cancer subset enriched in epithelial-to-mesenchymal transition and stem cell characteristics. Cancer Res 2009 May 15;69(10):4116–4124. pmid:19435916

5. Sabatier R, Finetti P, Guille A, Adelaide J, Chaffanet M, Viens P, et al. Claudin-low breast cancers: clinical, pathological, molecular and prognostic characterization. Mol Cancer 2014 Oct 2;13:228-4598-13-228.

6. Farmer P, Bonnefoi H, Becette V, Tubiana-Hulin M, Fumoleau P, Larsimont D, et al. Identification of molecular apocrine breast tumours by microarray analysis. Oncogene 2005 Jul 7;24(29):4660–4671. pmid:15897907

7. Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature 2000 Aug 17;406(6797):747–752. pmid:10963602

8. Dvorkin-Gheva A, Hassell JA. Identification of a novel luminal molecular subtype of breast cancer. PLoS One 2014 Jul 30;9(7):e103514. pmid:25076125

9. McCall MN, Bolstad BM, Irizarry RA. Frozen robust multiarray analysis (fRMA). Biostatistics 2010 Apr;11(2):242–253. pmid:20097884

10. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2007 Jan;8(1):118–127. pmid:16632515

11. Benito M, Parker J, Du Q, Wu J, Xiang D, Perou CM, et al. Adjustment of systematic microarray data biases. Bioinformatics 2004 Jan 1;20(1):105–114. pmid:14693816

12. Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol 2009 Mar 10;27(8):1160–1167. pmid:19204204

13. Smyth GK. Limma: linear models for microarray data. In: Gentleman R, Carey V, Duboit S, Irizarry R, Huber W, editors. Bioinformatics and Computational Biology Solutions using R and Bioconductor New York: Springer; 2005.

14. Benjamini YY,D. The control of the false discovery rate in multiple testing under dependency. Annals of Statistics 2001;29:1165–1188.

15. Heiser LM, Sadanandam A, Kuo WL, Benz SC, Goldstein TC, Ng S, et al. Subtype and pathway specific responses to anticancer compounds in breast cancer. Proc Natl Acad Sci U S A 2012 Feb 21;109(8):2724–2729. pmid:22003129

16. Bane AL, Whelan TJ, Pond GR, Parpia S, Gohla G, Fyles AW, et al. Tumor Factors Predictive of Response to Hypofractionated Radiotherapy in a Randomized Trial Following Breast Conserving Therapy. Ann Oncol 2014 Feb 22.

17. Whelan TJ, Pignol JP, Levine MN, Julian JA, MacKenzie R, Parpia S, et al. Long-term results of hypofractionated radiation therapy for breast cancer. N Engl J Med 2010 Feb 11;362(6):513–520. pmid:20147717

18. Lakhani SR, Ellis IO, Schnitt SJ, Tan PH, Van De Vijver MJ. WHO Classification of Tumours of the Breast.; 2012.

19. Elston CW, Ellis IO. 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 Nov;19(5):403–410. pmid:1757079

20. Allred DC, Harvey JM, Berardo M, Clark GM. Prognostic and predictive factors in breast cancer by immunohistochemical analysis. Mod Pathol 1998 Feb;11(2):155–168. pmid:9504686

21. Harvey JM, Clark GM, Osborne CK, Allred DC. Estrogen receptor status by immunohistochemistry is superior to the ligand-binding assay for predicting response to adjuvant endocrine therapy in breast cancer. J Clin Oncol 1999 May;17(5):1474–1481. pmid:10334533

22. Mohsin SK, Weiss H, Havighurst T, Clark GM, Berardo M, Roanh le D, et al. Progesterone receptor by immunohistochemistry and clinical outcome in breast cancer: a validation study. Mod Pathol 2004 Dec;17(12):1545–1554. pmid:15272277

23. Wolff AC, Hammond ME, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. J Clin Oncol 2013 Nov 1;31(31):3997–4013. pmid:24101045

24. Cheang MC, Chia SK, Voduc D, Gao D, Leung S, Snider J, et al. Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. J Natl Cancer Inst 2009 May 20;101(10):736–750. pmid:19436038

25. Liu CL, Montgomery KD, Natkunam Y, West RB, Nielsen TO, Cheang MC, et al. TMA-Combiner, a simple software tool to permit analysis of replicate cores on tissue microarrays. Mod Pathol 2005 Dec;18(12):1641–1648. pmid:16258508

26. Nielsen TO, Hsu FD, Jensen K, Cheang M, Karaca G, Hu Z, et al. Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin Cancer Res 2004 Sep 15;10(16):5367–5374. pmid:15328174

27. Rody A, Karn T, Liedtke C, Pusztai L, Ruckhaeberle E, Hanker L, et al. A clinically relevant gene signature in triple negative and basal-like breast cancer. Breast Cancer Res 2011 Oct 6;13(5):R97. pmid:21978456

28. Li C, Ma H, Wang Y, Cao Z, Graves-Deal R, Powell AE, et al. Excess PLAC8 promotes an unconventional ERK2-dependent EMT in colon cancer. J Clin Invest 2014 May;124(5):2172–2187. pmid:24691442

29. Taube JH, Herschkowitz JI, Komurov K, Zhou AY, Gupta S, Yang J, et al. Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes. Proc Natl Acad Sci U S A 2010 Aug 31;107(35):15449–15454. pmid:20713713

30. Asiedu MK, Ingle JN, Behrens MD, Radisky DC, Knutson KL. TGFbeta/TNF(alpha)-mediated epithelial-mesenchymal transition generates breast cancer stem cells with a claudin-low phenotype. Cancer Res 2011 Jul 1;71(13):4707–4719. pmid:21555371

31. Morel AP, Hinkal GW, Thomas C, Fauvet F, Courtois-Cox S, Wierinckx A, et al. EMT inducers catalyze malignant transformation of mammary epithelial cells and drive tumorigenesis towards claudin-low tumors in transgenic mice. PLoS Genet 2012;8(5):e1002723. pmid:22654675

32. Charafe-Jauffret E, Ginestier C, Iovino F, Tarpin C, Diebel M, Esterni B, et al. Aldehyde dehydrogenase 1-positive cancer stem cells mediate metastasis and poor clinical outcome in inflammatory breast cancer. Clin Cancer Res 2010 Jan 1;16(1):45–55. pmid:20028757

33. Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF. Prospective identification of tumorigenic breast cancer cells. Proc Natl Acad Sci U S A 2003 Apr 1;100(7):3983–3988. pmid:12629218

34. Ginestier C, Hur MH, Charafe-Jauffret E, Monville F, Dutcher J, Brown M, et al. ALDH1 Is a Marker of Normal and Malignant Human Mammary Stem Cells and a Predictor of Poor Clinical Outcome. Cell Stem Cell 2007 Nov 15;1(5):555–567. pmid:18371393

35. Bane A, Viloria-Petit A, Pinnaduwage D, Mulligan AM, O'Malley FP, Andrulis IL. Clinical-pathologic significance of cancer stem cell marker expression in familial breast cancers. Breast Cancer Res Treat 2013 Jul;140(1):195–205. pmid:23813303

36. Creighton CJ, Chang JC, Rosen JM. Epithelial-mesenchymal transition (EMT) in tumor-initiating cells and its clinical implications in breast cancer. J Mammary Gland Biol Neoplasia 2010 Jun;15(2):253–260. pmid:20354771

37. Herschkowitz JI, Zhao W, Zhang M, Usary J, Murrow G, Edwards D, et al. Comparative oncogenomics identifies breast tumors enriched in functional tumor-initiating cells. Proc Natl Acad Sci U S A 2012 Feb 21;109(8):2778–2783. pmid:21633010

38. Bruna A, Greenwood W, Le Quesne J, Teschendorff A, Miranda-Saavedra D, Rueda OM, et al. TGFbeta induces the formation of tumour-initiating cells in claudinlow breast cancer. Nat Commun 2012;3:1055. pmid:22968701

39. Hollier BG, Tinnirello AA, Werden SJ, Evans KW, Taube JH, Sarkar TR, et al. FOXC2 expression links epithelial-mesenchymal transition and stem cell properties in breast cancer. Cancer Res 2013 Mar 15;73(6):1981–1992. pmid:23378344

40. Knezevic J, Pfefferle AD, Petrovic I, Greene SB, Perou CM, Rosen JM. Expression of miR-200c in claudin-low breast cancer alters stem cell functionality, enhances chemosensitivity and reduces metastatic potential. Oncogene 2015 Mar 9.

41. Denkert C, von Minckwitz G, Brase JC, Sinn BV, Gade S, Kronenwett R, et al. Tumor-Infiltrating Lymphocytes and Response to Neoadjuvant Chemotherapy With or Without Carboplatin in Human Epidermal Growth Factor Receptor 2-Positive and Triple-Negative Primary Breast Cancers. J Clin Oncol 2014 Dec 22.

42. Loi S, Michiels S, Salgado R, Sirtaine N, Jose V, Fumagalli D, et al. Tumor infiltrating lymphocytes are prognostic in triple negative breast cancer and predictive for trastuzumab benefit in early breast cancer: results from the FinHER trial. Ann Oncol 2014 Apr 29.

43. Denkert C, Loibl S, Noske A, Roller M, Muller BM, Komor M, et al. Tumor-associated lymphocytes as an independent predictor of response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol 2010 Jan 1;28(1):105–113. pmid:19917869

44. Adams S, Gray RJ, Demaria S, Goldstein L, Perez EA, Shulman LN, et al. Prognostic Value of Tumor-Infiltrating Lymphocytes in Triple-Negative Breast Cancers From Two Phase III Randomized Adjuvant Breast Cancer Trials: ECOG 2197 and ECOG 1199. J Clin Oncol 2014 Jul 28.