- Research article
- Open Access
Gene-expression profiling of microdissected breast cancer microvasculature identifies distinct tumor vascular subtypes
Breast Cancer Research volume 14, Article number: R120 (2012)
Angiogenesis represents a potential therapeutic target in breast cancer. However, responses to targeted antiangiogenic therapies have been reported to vary among patients. This suggests that the tumor vasculature may be heterogeneous and that an appropriate choice of treatment would require an understanding of these differences.
To investigate whether and how the breast tumor vasculature varies between individuals, we isolated tumor-associated and matched normal vasculature from 17 breast carcinomas by laser-capture microdissection, and generated gene-expression profiles. Because microvessel density has previously been associated with disease course, tumors with low (n = 9) or high (n = 8) microvessel density were selected for analysis to maximize heterogeneity for this feature.
We identified differences between tumor and normal vasculature, and we describe two subtypes present within tumor vasculature. These subtypes exhibit distinct gene-expression signatures that reflect features including hallmarks of vessel maturity. Potential therapeutic targets (MET, ITGAV, and PDGFRβ) are differentially expressed between subtypes. Taking these subtypes into account has allowed us to derive a vascular signature associated with disease outcome.
Our results further support a role for tumor microvasculature in determining disease progression. Overall, this study provides a deeper molecular understanding of the heterogeneity existing within the breast tumor vasculature and opens new avenues toward the improved design and targeting of antiangiogenic therapies.
The growth of tumors beyond a certain size requires the recruitment of an adequate blood supply, which is supplied by abnormal angiogenesis. This involves the triggering of an "angiogenic switch" , whereby the tumor microenvironment enters a proangiogenic mode in response to hypoxia. This process is accompanied by increased levels of multiple proangiogenic factors, including vascular endothelial growth factor A (VEGFA) and platelet-derived growth factor B (PDGFB), as well as decreases in antiangiogenic factors such as endostatin. This leads to an increase in the number of proliferating endothelial cells, along with enhanced endothelial cell recruitment and migration toward the tumor bed .
Because endothelial cells are considered to be genetically stable , modulation of angiogenic processes is considered to be a promising area for cancer therapy. However, clinical trials of antiangiogenic therapies in breast cancer have reported mixed results [3–5]. Heterogeneity among tumor vasculatures may partially explain this lack of consistency . The most common measure of vascular heterogeneity is microvessel density (MVD), which has weak prognostic value in breast cancer . MVD is significantly associated with histologic type, but does not appear to be linked to other features, including nodal status, tumor grade, estrogen receptor (ER) status, molecular subtype, or the presence of vascular invasion [7, 8]. However, other characteristics of the tumor vasculature, such as vascular proliferation and levels of circulating endothelial progenitor cells, also have clinical significance in cancer [9, 10]. Therefore, a more complete understanding of vascular heterogeneity is required.
Gene-expression profiling studies have provided a greater understanding of breast tumor heterogeneity [11, 12] and have highlighted the role of the stroma in influencing disease outcome [13–16]. However, previous studies of tumor vasculature have mostly been limited to small numbers of samples and have focused on direct comparisons between tumor-associated and normal vasculature [17–23]. Although these studies have identified alterations characteristic of the tumor vasculature per se, they have not addressed the question of vascular heterogeneity or investigated whether distinct intertumoral vascular subtypes exist. Here, we use laser-capture microdissection (LCM) to isolate tumor-associated and matched normal microvascular compartments separately from 17 human breast cancer samples. After gene-expression profiling, we identified and validated gene-expression signatures that define two distinct subtypes of tumor vasculature, an understanding of which would represent an important step toward the improved design and targeting of therapeutic modalities.
Materials and methods
Tissue samples from 21 patients undergoing surgery for primary invasive ductal carcinoma (IDC) with no prior neoadjuvant therapy, and from one patient undergoing reduction mammoplasty, were subjected to LCM (see Additional File 1, Table S1A, B). From this cohort, we obtained 17 samples of tumor-associated vasculature; in 15 cases, we also obtained normal vasculature. After determination of MVD (see Additional File 2, Figure S1A), PECAM1-positive cells (microvasculature) were collected by using LCM from within tumor beds (n = 17), as well as from morphologically normal tissue adjacent to breast tumors (n = 14) and from one reduction mammoplasty (sample V) (see Additional File 2, Figure S1B). Only PECAM1-positive cells present within capillaries and small vascular structures were collected, avoiding isolated PECAM1-positive cells and vessels with large lumens. Total RNA was extracted, amplified, labeled, and hybridized on Agilent 44K Whole Genome arrays (Agilent Technologies, Santa Clara, CA). The endothelial and epithelial cell content of isolated samples was assessed by quantitative real-time PCR (qRT-PCR) (see Additional File 2, Figure S1C).
Description of samples
Clinical information was collected on a prospective basis (median follow-up, 5.6 years; see Additional File 1, Table S1B). Occurrence of distant recurrence was determined by examination of medical records. Matched normal epithelial samples were obtained as previously described . This study was approved by the McGill University Health Centre Research Ethics Board (Protocols SDR-99-780 and SDR-00-966). All patients provided written, informed consent.
Laser-capture microdissection, RNA isolation, and microarray hybridization
Immunohistochemistry directed against PECAM1 for MVD quantitation was carried out as per the antibody manufacturer's instructions. Quantitation of PECAM1 staining density was performed by averaging the stained pixel intensity of three fields captured at 10 × magnification by using imageJ . Details of the anti-PECAM1 immunohistochemistry protocol used for guided LCM are presented in Supplementary Information (see Additional File 3, Supplementary Materials and Methods). LCM, RNA isolation and sample preparation, as well as microarray hybridization, were carried out as previously described . Microarray results have been submitted to the Gene Expression Omnibus (GSE15363).
Differential expression, class discovery, and class prediction
All analyses were performed using R/Bioconductor version 2.6 . Microarray data were extracted as described in Supplementary Information (see Additional File 3, Supplementary Materials and Methods). Initially, the set of most-variable probes (n = 1,168) was defined as those with an interquartile range above 2. Class discovery was performed by hierarchical clustering based on this set, using Ward's minimum-variance method with a correlation distance metric. The significance of the resulting clusters was calculated by using the pvclust package with 10,000 bootstrap iterations. Heatmaps were generated by scaling each row (gene) into Z-scores, subtracting the mean and dividing by its standard deviation. Differential expression between sample sets was determined by using Linear Models for Microarray Analysis (LIMMA) . For genes represented by multiple probes, only the probe with the largest interquartile range was used. Unless stated otherwise, a gene was considered differentially expressed if the FDR-adjusted P value was ≤0.05 .
Predictors of recurrence were generated as follows: gene expression for subtype B member profiles was adjusted by removing the average differences between subtype B and subtype A tumor vascular expression for each gene. Then, the six most significantly differentially expressed genes between recurrent and nonrecurrent tumor vasculature were identified by using a Wilcoxon rank-sum test. These genes were used to train a linear discriminant analysis classifier and tested for accuracy through leave-one-out cross-validation (LOOCV). All steps (including the adjustment of tumor vascular subtypes) were done under cross-validation. Gene sets of size two to 100 were tested for prediction, and six-gene predictors were found to exhibit the highest accuracy under LOOCV.
Quantitative real-time PCR and immunohistochemistry
Validation for selected genes and cognate proteins was performed by qRT-PCR and IHC, respectively, on subsets of samples. Full details and probe sequences are presented in Supplementary Information (see Additional File 3, Supplementary Materials and Methods).
Results and discussion
Because our goal was to investigate the heterogeneity of the tumor microvasculature in invasive ductal carcinomas, samples with high and low MVD were chosen to maximize the spectrum of heterogeneity present within the sample cohort. Within this set, we identified tumor samples with high (n = 8) or low (n = 9) MVD, as determined by immunohistochemistry directed against the endothelial marker PECAM1 (see Additional File 1; Table S1B; Additional File 2, Figure S1A).
Isolation of vasculature-enriched samples and gene-expression profiling
A significant enrichment of PECAM1 was observed in microdissected vasculature when compared with the epithelial cell marker cytokeratin 8 (KRT8); conversely, the latter marker was enriched in matched microdissected epithelial RNA samples (see Additional File 2, Figure S1C). These results indicate that the vascular samples are highly enriched in endothelial cell content.
Tumor and normal vasculatures possess distinct expression profiles
To investigate the differences among the expression profiles generated, we performed hierarchical clustering by using the 1,168 most variable genes (interquartile range >2) over the entire dataset. The normal and tumor vascular samples formed two distinct clusters (Figure 1A; see Additional File 2, Figure S2A). Further comparison of the tumor and normal vasculature samples identified 494 differentially expressed genes (FDR < 0.05, LIMMA; see Additional File 4, Sheet 1).
A significant overlap exists between the gene sets discriminating between tumor and normal vasculature identified by previous studies [17–23] and here (P < 0.05 for each, Fisher Exact test; Figure S3A). This overlap is significant even for the studies performed on different tumor types [19–23] or by using other isolation techniques [18–20, 23]. Common processes include extracellular matrix remodeling and cytokine signaling  (see Additional File 1, Table S2). Several genes previously reported to be elevated in tumor vasculature  are also elevated in our tumor vascular cohort, including osteonectin (SPARC), collagens 1A1, 3A1, 4A2, and 18A1 (COL1A1, COL3A1, COL4A2 and COL18A1), and Thy cell-surface antigen (THY1) (see Additional File 4, Sheet 1). This is consistent with the observation that tumor vasculature is constantly reorganizing [30, 31] and confirms that our normal and tumor vasculature are comparable to those isolated in previous studies.
The tumor vasculature also exhibits differential expression of genes related to oxidative phosphorylation/stress and MAPK signaling (Figure 1B; see Additional File 1, Table S2; Additional File 4, Sheet 1). It is known that reoxygenation of vasculature, consistent with the variable flow of blood in tumor vessels, leads to the creation of reactive oxygen species .
Tumor vascular expression profiles segregate into two distinct subtypes reflecting different biological processes
Interestingly, the hierarchical clustering revealed that the tumor samples form two distinct clusters, here termed A and B (Figure 1A; see Additional File 1, Figure S2A); this separation was also reflected in the first component when the set of tumor vascular samples were subjected to principal-component analysis (see Additional File 1, Figure S2B). Analysis of differential expression between these two clusters identifies 1,659 genes (FDR < 0.05, LIMMA; see Additional File 4, Sheet 2). This is more than 3 times larger than the number of genes differentially expressed between normal and tumor samples, emphasizing the significance of the clusters and indicating that homogeneity within each cluster is higher than that within the tumor vasculature as a whole. Thus, separating tumor samples by cluster may prove more informative than grouping them together.
The tumor vascular clusters are independent of major clinical variables, including ER and HER2 status, grade, and lymph node involvement (P ≥ 0.6; see Additional File 1, Table S1A) but are enriched for larger and smaller tumors as well as high and low MVD (P <0.05, t test). Interestingly, no statistically significant differential gene expression was found between high- and low-MVD samples (FDR < 0.05, LIMMA). This demonstrates that although MVD status is associated with the tumor vascular clusters identified here, additional properties of the samples are primarily responsible for the molecular profiles and observed clustering.
To identify biological functions specifically associated with each of the two subtypes of tumor vasculature, we identified the most highly differentially expressed genes (FDR < 0.0001, LIMMA, Figure 1C; FDR < 0.05, LIMMA; see Additional File 4, Sheet 2) and pathways (FDR < 0.05; Figure 1B; see Additional File 1, Table S2) between normal vasculature and each of the tumor vascular subtypes A and B.
Genes within Group 1 exhibit the highest expression in normal samples and the lowest expression in tumor-associated vasculature subtype A members (Figure 1C). Interestingly, these include AOC3 and FOXO1, both associated with vascular shear stress [33, 34]. Reduced shear stress suggests a decreased flow rate, consistent with the hypothesis that tumor vessels, and in particular those belonging to subtype A, lack appropriate perfusion . Many genes and processes related to metabolism and biosynthesis are also overrepresented in subtype A (see Additional File 1, Table S2; Additional File 4, Sheet 2), suggesting active proliferation. Consistent with this, Group 2, consisting of genes whose expression is elevated in subtype A (Figure 1C), includes cathepsin B (CTSB), associated with active vascular remodeling [35, 36].
Genes within Group 3 are more highly expressed in subtype B, and are enriched for elements involved in interactions with the extracellular matrix, collagen production, focal adhesion, hypoxia, glycolysis, immune response, and protein export, suggesting that this group may represent a more stable vasculature (Figure 1C; Additional File 1, Table S2). These include genes involved in antiangiogenic processes, such as thrombospondin 2 (THBS2) and the collagens COL18A1 and COL4A2, which can be cleaved to form endostatin and canstatin, respectively . Antiangiogenic signaling from these and other genes (see Additional File 1, Figure S3B) may contribute to the decreased overall MVD seen in these samples. Interestingly, NOTCH3 is also elevated in subtype B (see Additional File 4, Sheet 2); Notch signaling is implicated in vascular development and maturation [38, 39].
Tumor vasculature is often considered immature with respect to normal vasculature , the known functions of many of the individual genes identified here suggest that the two tumor vascular clusters may represent differences in degree of vascular maturity. This hypothesis was investigated by using several approaches. Pericytes line mature capillary vessels and would be co-isolated with PECAM1-positive endothelial cells during LCM, whereas the association of pericytes with immature vessels is relatively loose, and this cell type would therefore be co-isolated at a lower frequency in such vessels . Therefore, we explored whether clusters A and B differed with respect to their content of pericyte-specific signals, specifically the pericyte markers ACTA2, PDGFRβ, and RGS5 , as well as LAMB1, a surrogate of laminin-8 (α4β1γ1) implicated in vessel maturation  (Figure 1D). Interestingly, these markers are elevated in cluster B (P <0.05 for each; t test), and their expression pattern successfully recapitulates the clusters (Figure 1D). Differential expression of PDGFRβ was confirmed by qRT-PCR on amplified RNA prepared from tumor vascular samples (see Additional File 2, Figure S3C), whereas immunohistochemistry (IHC)-based validation was also performed for ACTA2 and LAMB1 (Figure 2A, B; see Additional File 2, Figure S4A). Serial staining directed against PECAM1 showed that most ACTA2 and LAMB1 expression is observed in proximity to the PECAM1-expressing vasculature (see Additional File 2, Figure S4B, C), while pixel counts over whole sections (normalized for PECAM1 staining, as assessed on a neighboring section) showed a significantly increased expression of ACTA2 and LAMB1 in tumors belonging to cluster B (P <0.05 for each, Wilcoxon test; Figure 2C; Additional File 1, Table S1A). Therefore, the differential gene expression observed between clusters is recapitulated at the protein level for these markers. In addition, assessment of LAMB1 mRNA content in selected samples by qRT-PCR significantly correlated with the microarray-derived values; this also held true for several other selected genes (see Additional File 2, Figure S5). Supporting this observation, a negative correlation between MVD and pericyte coverage has also been identified in endometrial cancer .
To assess whether the subtype B-specific gene-expression profile is driven by pericyte-specific genes, the sample orderings induced by pericyte marker expression and subtype-specific gene expression were compared (see Additional File 2, Figure S6). First, we confirmed that ordering of all samples by relative expression of known pericyte markers recapitulated the separation of subtype A and B members (Additional File 2, Figure S6A, B); this ordering was not correlated with expression of endothelial marker genes (Additional File 2, Figure S6C). Subtype B members were then ordered by relative expression of known pericyte markers (Additional File 2, Figure S6D). Then, retaining this ordering, the relative expression of the top 200 genes differentially expressed between the A and B subtypes (ranked by P value) was examined (see Additional File 2, Figure S6E). Importantly, no correlation was seen at the sample level between relative pericyte-marker content and relative expression levels of subtype-specific genes within the set of subtype B members, supporting that processes other than those linked to pericyte-marker content distinguish these two subtypes.
Pathway analysis highlights an increased expression of genes related to glycolysis and oxidative phosphorylation in subtype B (Figure 1B). This is intriguing, as elevation of either glycolysis or oxidative phosphorylation often leads to downregulation of the other . Coupled with an increased expression of ribosomal genes, the increase in oxidative phosphorylation could be evidence of increased metabolism in subtype B vasculature. Interestingly, the different tumor vascular subtypes also vary at the level of NK and T-cell marker expression, suggesting that vascular subtypes can affect the tumor microenvironment. Additionally, CXCL12/SDF-1, considered to play a role in pericyte recruitment , is decreased in subtype A tumor vasculature relative to subtype B, consistent with the presence of increased levels of pericyte markers in subtype B, whereas angiopoietin-2 (ANGPT2), thought to promote blood-vessel destabilization by acting as an ANGPT1 antagonist , is elevated in subtype B tumor vasculature.
Together, these results strongly support the hypothesis that the vascular clusters observed here reflect differences in vessel maturity, and indicate that related processes, including anti- and proangiogenic regulators and matrix remodeling, also contribute to the segregation of the vascular subtypes. We suggest that cluster B, associated with low MVD, may represent samples with a more mature tumor vasculature with respect to members of cluster A. Other potential factors driving cluster segregation may include unidirectional shear stress (USS), expected to be highest in vessels of subtype B. USS elicits an antiinflammatory response  and has been shown to induce miR-21 expression in endothelial cells, leading to decreased apoptosis and activation of the nitric oxide pathway, implicated in Notch signaling . All of these are consistent with our proposed subtypes (Figure 1B). Interestingly, genes previously identified as being commonly overexpressed in tumor endothelium and therefore proposed as specific markers for endothelial cell targeting, including SPARC and COL4A1 , are elevated in subtype B but not subtype A tumor vasculature, accentuating the importance of gaining a better understanding of the heterogeneity present in this compartment.
Most previous signatures of tumor vascular cells identify only tumor vascular subtype B
Previous studies (see Additional File 2, Figure S3A) have generated signatures that segregate tumor versus normal vascular samples, but have not identified distinct subtypes within the tumor vasculature. To test whether these signatures can differentiate between our tumor vascular subtypes, we examined the clustering induced by previously reported tumor vascular signatures [17, 18, 21–23] in our dataset. All previously published tumor vascular gene sets investigated could segregate members of subtype B from normal samples in our data (Figure 3; Additional File 2, Figure S3D-G); however, only the gene set described by Bhati et al.  generated from breast tumors using LCM, successfully segregates the subtype A and normal samples (Figure 3A). The inability of most existing datasets to generate signatures capable of distinguishing between tumor vascular subtypes may be partially explained by differences in the isolation technique used. Unlike LCM, immunomagnetic bead-based approaches using dispersed cells fail to isolate cells adjacent to the vessels that can provide additional information about the vascular microenvironment. Additionally, differences in the tissue of origin may contribute to these observations (see Additional File 2, Figure S3A), suggesting that subtype A may be enriched in breast tumors, and supporting the concept of tissue-specific antiangiogenic approaches.
The ability of the tumor vascular signature from Bhati et al.  to partially recapitulate our tumor vascular subtypes (Figure 3A) suggests that their samples contained features of both subtypes, and indicates that this vascular heterogeneity is likely to be common in breast cancer. Bhati et al.  used laser-capture microdissection to isolate factor VIII-expressing vascular cells from human breast cancers and normal breast tissue; however, they failed to identify subtypes within the tumor sample set. This may partly be due to the restricted sample size used in their study (n = 5), whereas the increased size of our patient cohort (n = 17), in combination with our selection based on MVD, likely enabled us to capture a broader range of interpatient variation in tumor vasculature. The small number of samples investigated in most previous studies (n = 1 to 5) [17–20, 22, 23] also raises the possibility that only members of one subtype were present in their study cohorts. Hence our identification and initial characterization of distinct breast cancer vasculature subtypes is a novel finding that emphasizes the concept that breast cancer heterogeneity extends to the different subcomponents of the tumor microenvironment.
Gene-expression differences induced by subtype membership mask differential gene expression associated with disease outcome
To determine whether differences in tumor vasculature gene expression were correlated with other aspects of tumor biology and disease course, we performed differential expression analysis between the five recurrent and 12 nonrecurrent tumor samples in our dataset. This identified a total of 128 genes (P < 0.01, no adjustment: see Additional File 1, Table S3), including TLR2 (involved in innate immunity)  and GAB1 (a scaffold protein linked to VEGF signaling and endothelial migration) . Subtype membership is not correlated with recurrence in our relatively small dataset. This agrees with previous reports that differences in vessel maturity are associated with lymph node status but not other clinical variables in breast cancer . However, we speculated that the differences in gene expression induced by subtype membership might be sufficiently large to mask differential gene expression associated with disease outcome. Similar conclusions have been drawn regarding gene expression-based predictors of breast cancer outcome based on whole-tumor data, which have generally been found to be more reliable when limited to specific tumor subtypes [51, 52]. Generation of subtype-specific predictors was not possible, given the small size of each group and the limited number of recurrences in each (see Additional File 1, Table S1A,B). Therefore, we investigated this hypothesis by removing the gene-expression changes that could be ascribed to subtype membership. This was achieved by calculating the difference in average expression between the members of the two subtypes for each probe, and subtracting this value from the relevant probe data for all members of one subtype. After this adjustment, we could successfully build a six-gene predictor that identified recurrent patient samples with 94% accuracy under cross-validation in our dataset (P = 0.002; Fisher Exact test; Figure 4A). Without applying this correction, no successful predictor could be constructed. The fact that the combination of these genes can predict recurrence under cross-validation within our dataset confirms their differential expression, despite their individual lack of significance after FDR adjustment.
This indicates that vascular markers of recurrence can be subtype-specific and that the subtypes contain different markers linked to recurrence. Additionally, the observation that successful construction of a predictor cannot be achieved without taking the inter-subtype differences into account suggests that these differences are sufficiently large to mask signals associated with other clinical variables, and supports that tumor vascular subtypes must be taken into account in future analyses of this compartment.
Our earlier work established that breast tumor stroma as a whole contains sufficient information to predict recurrence independently . This study highlights that features of patient-specific vasculatures can reflect differences in the tumor microenvironment predictive of recurrence. As such, the tumor vasculature, including its subtypes, should be considered when studying the determinants of recurrence. This adds another dimension to previous reports that different characteristics of vasculature are associated with tumor stage and metastasis [6, 9, 10, 53], and indicates that the effects of tumor vasculature on disease progression are complex. Further research will be necessary to identify successfully all relevant factors contributed by the vascular environment in tumor progression.
Drug targets are differentially expressed between subtypes
The existence of vascular subtypes has important implications for the choice of antiangiogenic treatment strategies. A substantial proportion of patients in such clinical trials fail to show clinical benefits [3, 4]; short-term treatment with one antiangiogenic drug led to an increased incidence of metastasis . Importantly, observed differences between our subtypes (for example, pericyte coverage) may reflect mechanisms of resistance to antiangiogenic therapy . Pericyte coverage can reduce tumor vessel vulnerability to VEGF inhibition, and in such cases, targeting treatments against both pericytes and endothelial cells may lead to improved efficacy .
The targets of several antiangiogenic drugs [56, 57] are also differentially expressed between subtypes (Figure 4B). Both PDGFRβ and MET (a receptor tyrosine kinase, whose natural ligand HGF is a potent angiogenic factor ), are elevated in subtype B; antagonists to αv integrin (ITGAV) , also elevated in subtype B, have been shown to block neovascularization. Vessel maturity has been linked to resistance to antiangiogenic therapy in melanoma , suggesting that this mechanism may be applicable across tumor types. The differential expression of drug targets indicates that subtype-specific approaches to targeting tumor angiogenesis may lead to improved response rates by permitting prestratification of patient populations. Modifications of the tumor vasculature, as suggested by a recent study identifying vessel normalization through ablation of the endothelial oxygen sensor PHD2 as a potential means to improve the delivery and efficacy of chemotherapeutic agents , may also have subtype-specific effects. The identification of biomarkers for prediction of response to antiangiogenic therapy is a current critical clinical need , and the vascular-subtype gene-expression profiles presented here provide new perspectives for the validation of such markers.
Overall, our data illustrate the heterogeneity of the breast tumor vasculature, define two tumor vascular subtypes that are associated with, but not defined by, microvessel density, and represent a first step toward establishing the landscape of vascular subtypes in breast cancer. Our previous work has shown that heterogeneity within the stromal compartment of breast tumors is linked to disease outcome . We now suggest that breast tumor vasculature heterogeneity also plays a role in determining disease course, and that further study thereof may be critical for improved patient stratification before selection of specific targeted therapeutic regimens.
false discovery rate
human epidermal growth factor receptor 2
Kyoto Encyclopedia of Genes and Genomes
quantitative real-time polymerase chain reaction.
Hanahan D, Folkman J: Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis. Cell. 1996, 86: 353-364. 10.1016/S0092-8674(00)80108-7.
Folkman J: Angiogenesis. Annu Rev Med. 2006, 57: 1-18. 10.1146/annurev.med.57.121304.131306.
Longo R, Gasparini G: Challenges for patient selection with VEGF inhibitors. Cancer Chemother Pharmacol. 2007, 60: 151-170. 10.1007/s00280-006-0403-6.
Bergers G, Hanahan D: Modes of resistance to anti-angiogenic therapy. Nat Rev Cancer. 2008, 8: 592-603. 10.1038/nrc2442.
Carpini JD, Karam AK, Montgomery L: Vascular endothelial growth factor and its relationship to the prognosis and treatment of breast, ovarian, and cervical cancer. Angiogenesis. 2010, 13: 43-58. 10.1007/s10456-010-9163-3.
Uzzan B, Nicolas P, Cucherat M, Perret G-Y: Microvessel density as a prognostic factor in women with breast cancer: a systematic review of the literature and meta-analysis. Cancer Res. 2004, 64: 2941-2955. 10.1158/0008-5472.CAN-03-1957.
Lopes N, Sousa B, Vieira D, Milanezi F, Schmitt F: Vessel density assessed by endoglin expression in breast carcinomas with different expression profiles. Histopathology. 2009, 55: 594-599. 10.1111/j.1365-2559.2009.03417.x.
Marinho A, Soares R, Ferro J, Lacerda M, Schmitt FC: Angiogenesis in breast cancer is related to age but not to other prognostic parameters. Pathol Res Pract. 1997, 193: 267-273. 10.1016/S0344-0338(97)80003-9.
Stefansson IM, Salvesen HB, Akslen LA: Vascular proliferation is important for clinical progress of endometrial cancer. Cancer Res. 2006, 66: 3303-3309. 10.1158/0008-5472.CAN-05-1163.
Naik RP, Jin D, Chuang E, Gold EG, Tousimis EA, Moore AL, Christos PJ, de Dalmas T, Donovan D, Rafii S, Vahdat LT: Circulating endothelial progenitor cells correlate to stage in patients with invasive breast cancer. Breast Cancer Res Treat. 2008, 107: 133-138.
Sørlie T, Perou C, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen M, van de Rijn M, Jeffrey S, Thorsen T, Quist H, Matese J, Brown P, Botstein D, Lønning PE, Børresen-Dale A: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.
Chin K, DeVries S, Fridlyand J, Spellman PT, Roydasgupta R, Kuo W-L, Lapuk A, Neve RM, Qian Z, Ryder T, Chen F, Feiler H, Tokuyasu T, Kingsley C, Dairkee S, Meng Z, Chew K, Pinkel D, Jain A, Ljung BM, Esserman L, Albertson DG, Waldman FM, Gray JW: Genomic and transcriptional aberrations linked to breast cancer pathophysiologies. Cancer Cell. 2006, 10: 529-541. 10.1016/j.ccr.2006.10.009.
Finak G, Bertos N, Pepin F, Sadekova S, Souleimanova M, Zhao H, Chen H, Omeroglu G, Meterissian S, Omeroglu A, Hallett M, Park M: Stromal gene expression predicts clinical outcome in breast cancer. Nat Med. 2008, 14: 518-527. 10.1038/nm1764.
Chang HY, Nuyten DS, Sneddon JB, Hastie T, Tibshirani R, Sorlie T, Dai H, He YD, van't Veer LJ, Bartelink H, van de Rijn M, Brown PO, van de Vijver MJ: Robustness, scalability, and integration of a wound-response gene expression signature in predicting breast cancer survival. Proc Natl Acad Sci USA. 2005, 102: 3738-3743. 10.1073/pnas.0409462102.
Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, Montgomery K, Chi JT, van de Rijn M, Botstein D, Brown PO: Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol. 2004, 2: E7-10.1371/journal.pbio.0020007.
Teschendorff AE, Miremadi A, Pinder SE, Ellis IO, Caldas C: An immune response gene expression module identifies a good prognosis subtype in estrogen receptor negative breast cancer. Genome Biol. 2007, 8: R157-10.1186/gb-2007-8-8-r157.
Bhati R, Patterson C, Livasy CA, Fan C, Ketelsen D, Hu Z, Reynolds E, Tanner C, Moore DT, Gabrielli F, Perou CM, Klauber-DeMore N: Molecular characterization of human breast tumor vascular cells. Am J Pathol. 2008, 172: 1381-1390. 10.2353/ajpath.2008.070988.
Parker BS, Argani P, Cook BP, Liangfeng H, Chartrand SD, Zhang M, Saha S, Bardelli A, Jiang Y, Martin TBS, Nacht M, Teicher BA, Klinger KW, Sukumar S, Madden SL: Alterations in vascular gene expression in invasive breast carcinoma. Cancer Res. 2004, 64: 7857-7866. 10.1158/0008-5472.CAN-04-1976.
Croix BS, Rago C, Velculescu V, Traverso G, Romans KE, Montgomery E, Lal A, Riggins GJ, Lengauer C, Vogelstein B, Kinzler KW: Genes expressed in human tumor endothelium. Science. 2000, 289: 1197-1202. 10.1126/science.289.5482.1197.
van Beijnum JR, Dings RP, van der Linden E, Zwaans BMM, Ramaekers FCS, Mayo KH, Griffioen AW: Gene expression of tumor angiogenesis dissected: specific targeting of colon cancer angiogenic vasculature. Blood. 2006, 108: 2339-2348. 10.1182/blood-2006-02-004291.
Buckanovich RJ, Sasaroli D, O'Brien-Jenkins A, Botbyl J, Hammond R, Katsaros D, Sandaltzopoulos R, Liotta LA, Gimotty PA, Coukos G: Tumor vascular proteins as biomarkers in ovarian cancer. J Clin Oncol. 2007, 25: 852-861. 10.1200/JCO.2006.08.8583.
Pen A, Moreno MJ, Martin J, Stanimirovic DB: Molecular markers of extracellular matrix remodeling in glioblastoma vessels: microarray study of laser-captured glioblastoma vessels. Glia. 2007, 55: 559-572. 10.1002/glia.20481.
Madden SL, Cook BP, Nacht M, Weber WD, Callahan MR, Jiang Y, Dufault MR, Zhang X, Zhang W, Walter-Yohrling J, Rouleau C, Akmaev VR, Wang CJ, Cao X, Martin TBS, Roberts BL, Teicher BA, Klinger KW, Stan R-V, Lucey B, Carson-Walter EB, Laterra J, Walter KA: Vascular gene expression in nonneoplastic and malignant brain. Am J Pathol. 2004, 165: 601-608. 10.1016/S0002-9440(10)63324-X.
Finak G, Sadekova S, Pepin F, Hallett M, Meterissian S, Halwani F, Khetani K, Souleimanova M, Zabolotny B, Omeroglu A, Park M: Gene expression signatures of morphologically normal breast tissue identify basal-like tumors. Breast Cancer Res. 2006, 8: R58-10.1186/bcr1608.
NIH ImageJ. [http://rsb.info.nih.gov/ij]
Suzuki R, Shimodaira H: Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006, 22: 1540-1542. 10.1093/bioinformatics/btl117.
Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B. 1995, 57: 289-300.
van Beijnum JR, Petersen K, Griffioen AW: Tumor endothelium is characterized by a matrix remodeling signature. Front Biosci (Schol Ed). 2009, 1: 216-225.
Jain RK: Molecular regulation of vessel maturation. Nat Med. 2003, 9: 685-693. 10.1038/nm0603-685.
Baluk P, Hashizume H, McDonald DM: Cellular abnormalities of blood vessels as targets in cancer. Curr Opin Genet Dev. 2005, 15: 102-111. 10.1016/j.gde.2004.12.005.
Houle F, Huot J: Dysregulation of the endothelial cellular response to oxidative stress in cancer. Mol Carcinog. 2006, 45: 362-367. 10.1002/mc.20218.
Chlench S, Disassa NM, Hohberg M, Hoffmann C, Pohlkamp T, Beyer G, Bongrazio M, Silva-Azevedo LD, Baum O, Pries AR, Zakrzewicz A: Regulation of Foxo-1 and the angiopoietin-2/Tie2 system by shear stress. FEBS Lett. 2007, 581: 673-680. 10.1016/j.febslet.2007.01.028.
Lalor PF, Edwards S, McNab G, Salmi M, Jalkanen S, Adams DH: Vascular adhesion protein-1 mediates adhesion and transmigration of lymphocytes on human hepatic endothelial cells. J Immunol. 2002, 169: 983-992.
Fukumura D, Jain RK: Tumor microenvironment abnormalities: causes, consequences, and strategies to normalize. J Cell Biochem. 2007, 101: 937-949. 10.1002/jcb.21187.
Premzl A, Turk V, Kos J: Intracellular proteolytic activity of cathepsin B is associated with capillary-like tube formation by endothelial cells in vitro. J Cell Biochem. 2006, 97: 1230-1240. 10.1002/jcb.20720.
Nyberg P, Xie L, Kalluri R: Endogenous inhibitors of angiogenesis. Cancer Res. 2005, 65: 3967-3979. 10.1158/0008-5472.CAN-04-2427.
Domenga V, Fardoux P, Lacombe P, Monet M, Maciazek J, Krebs LT, Klonjkowski B, Berrou E, Mericskay M, Li Z, Tournier-Lasserve E, Gridley T, Joutel A: Notch3 is required for arterial identity and maturation of vascular smooth muscle cells. Genes Dev. 2004, 18: 2730-2735. 10.1101/gad.308904.
Rehman AO, Wang C-Y: Notch signaling in the regulation of tumor angiogenesis. Trends Cell Biol. 2006, 16: 293-300. 10.1016/j.tcb.2006.04.003.
Mitchell TS, Bradley J, Robinson GS, Shima DT, Ng Y-S: RGS5 expression is a quantitative measure of pericyte coverage of blood vessels. Angiogenesis. 2008, 11: 141-151. 10.1007/s10456-007-9085-x.
Fujita M, Khazenzon NM, Bose S, Sekiguchi K, Sasaki T, Carter WG, Ljubimov AV, Black KL, Ljubimova JY: Overexpression of beta1-chain-containing laminins in capillary basement membranes of human breast cancer and its metastases. Breast Cancer Res. 2005, 7: R411-R421. 10.1186/bcr1011.
Williams KJ, Telfer BA, Airley RE, Peters HP, Sheridan MR, van der Kogel AJ, Harris AL, Stratford IJ: A protective role for HIF-1 in response to redox manipulation and glucose deprivation: implications for tumorigenesis. Oncogene. 2002, 21: 282-290. 10.1038/sj.onc.1205047.
Song N, Huang Y, Shi H, Yuan S, Ding Y, Song X, Fu Y, Luo Y: Overexpression of platelet-derived growth factor-BB increases tumor pericyte content via stromal-derived factor-1alpha/CXCR4 axis. Cancer Res. 2009, 69: 6057-6064.
Eklund L, Olsen BR: Tie receptors and their angiopoietin ligands are context-dependent regulators of vascular remodeling. Exp Cell Res. 2006, 312: 630-641. 10.1016/j.yexcr.2005.09.002.
Chien S: Mechanotransduction and endothelial cell homeostasis: the wisdom of the cell. Am J Physiol Heart Circ Physiol. 2007, 292: H1209-H1224.
Ishimura N, Bronk SF, Gores GJ: Inducible nitric oxide synthase up-regulates Notch-1 in mouse cholangiocytes: implications for carcinogenesis. Gastroenterology. 2005, 128: 1354-1368. 10.1053/j.gastro.2005.01.055.
van Beijnum JR, Petersen K, Griffioen AW: Tumor endothelium is characterized by a matrix remodeling signature. Front Biosci (Schol Ed). 2009, 1: 216-225.
Kim S, Takahashi H, Lin W-W, Descargues P, Grivennikov S, Kim Y, Luo J-L, Karin M: Carcinoma-produced factors activate myeloid cells through TLR2 to stimulate metastasis. Nature. 2009, 457: 102-106. 10.1038/nature07623.
Laramée M, Chabot C, Cloutier M, Stenne R, Holgado-Madruga M, Wong AJ, Royal I: The scaffolding adapter Gab1 mediates vascular endothelial growth factor signaling and is required for endothelial cell migration and capillary formation. J Biol Chem. 2007, 282: 7758-7769. 10.1074/jbc.M611327200.
Kakolyris S, Fox SB, Koukourakis M, Giatromanolaki A, Brown N, Leek RD, Taylor M, Leigh IM, Gatter KC, Harris AL: Relationship of vascular maturation in breast cancer blood vessels to vascular density and metastasis, assessed by expression of a novel basement membrane component, LH39. Br J Cancer. 2000, 82: 844-851. 10.1054/bjoc.1999.1010.
Sparano JA, Paik S: Development of the 21-gene assay and its application in clinical practice and clinical trials. J Clin Oncol. 2008, 26: 721-728. 10.1200/JCO.2007.15.1068.
Cardoso F, Veer LV, Rutgers E, Loi S, Mook S, Piccart-Gebhart MJ: Clinical application of the 70-gene profile: the MINDACT trial. J Clin Oncol. 2008, 26: 729-735. 10.1200/JCO.2007.14.3222.
Igreja C, Courinha M, Cachaço AS, Pereira T, Cabeçadas J, da Silva MG, Dias S: Characterization and clinical relevance of circulating and biopsy-derived endothelial progenitor cells in lymphoma patients. Haematologica. 2007, 92: 469-477. 10.3324/haematol.10723.
Ebos JM, Lee CR, Cruz-Munoz W, Bjarnason GA, Christensen JG, Kerbel RS: Accelerated metastasis after short-term treatment with a potent inhibitor of tumor angiogenesis. Cancer Cell. 2009, 15: 232-239. 10.1016/j.ccr.2009.01.021.
Sennino B, Kuhnert F, Tabruyn SP, Mancuso MR, Hu-Lowe DD, Kuo CJ, McDonald DM: Cellular source and amount of vascular endothelial growth factor and platelet-derived growth factor in tumors determine response to angiogenesis inhibitors. Cancer Res. 2009, 69: 4527-4536. 10.1158/0008-5472.CAN-08-3779.
Joyce JA: Therapeutic targeting of the tumor microenvironment. Cancer Cell. 2005, 7: 513-520. 10.1016/j.ccr.2005.05.024.
Gordon MS, Mendelson DS, Kato G: Tumor angiogenesis and novel antiangiogenic strategies. Int J Cancer. 2010, 126: 1777-1787.
You WK, McDonald DM: The hepatocyte growth factor/c-Met signaling pathway as a therapeutic target to inhibit angiogenesis. BMB Rep. 2008, 41: 833-839. 10.5483/BMBRep.2008.41.12.833.
Santulli RJ, Kinney WA, Ghosh S, Decorte BL, Liu L, Tuman RW, Zhou Z, Huebert N, Bursell SE, Clermont AC, Grant MB, Shaw LC, Mousa SA, Galemmo RA Jr, Johnson DL, Maryanoff BE, Damiano BP: Studies with an orally bioavailable alpha V integrin antagonist in animal models of ocular vasculopathy: retinal neovascularization in mice and retinal vascular permeability in diabetic rats. J Pharmacol Exp Ther. 2008, 324: 894-901.
Helfrich I, Scheffrahn I, Bartling S, Weis J, von Felbert V, Middleton M, Kato M, Ergun S, Schadendorf D: Resistance to antiangiogenic therapy is directed by vascular phenotype, vessel stabilization, and maturation in malignant melanoma. J Exp Med. 2010, 207: 491-503. 10.1084/jem.20091846.
Mazzone M, Dettori D, Leite de Oliveira R, Loges S, Schmidt T, Jonckx B, Tian YM, Lanahan AA, Pollard P, Ruiz de Almodovar C, De Smet F, Vinckier S, Aragones J, Debackere K, Luttun A, Wyns S, Jordan B, Pisacane A, Gallez B, Lampugnani MG, Dejana E, Simons M, Ratcliffe P, Maxwell P, Carmeliet P: Heterozygous deficiency of PHD2 restores tumor oxygenation and inhibits metastasis via endothelial normalization. Cell. 2009, 136: 839-851. 10.1016/j.cell.2009.01.020.
Jain RK, Duda DG, Willett CG, Sahani DV, Zhu AX, Loeffler JS, Batchelor TT, Sorensen AG: Biomarkers of response and resistance to antiangiogenic therapy. Nat Rev Clin Oncol. 2009, 6: 327-338. 10.1038/nrclinonc.2009.63.
The authors thank D. Fleiszer, A. Loutfi, C. Milne, F. Tremblay, M. Wexler (surgeons), F. Halwani (pathologist), D. Hori, L. Pasyuk, T. Vilhena (pathology staff), and the MUHC Anaesthesia Department for their assistance with obtaining samples, as well as S. Dumont for technical support with LCM, S. Cory for bioinformatics support, and R. Annan for editing. The microarrays were hybridized and scanned by the Genentech Microarray Facility and the Genentech Bioinformatics department. H.-P. Gerber coordinated the project at Genentech. This work was supported by grants from Meloche Monnex, the Québec Breast Cancer Foundation, and the Database and Tissue Bank Axis of the Réseau de recherche en cancer of the Fonds de la recherche en santé du Québec to MP; a Natural Science and Engineering Research Council of Canada Discovery Grants Program grant to MH; a Terry Fox New Frontiers Program grant to MP and MH; a U.S. Department of Defense Breast Cancer Predoctoral Traineeship Award to FP; Cedars Cancer Institute Fellowships to NB, JL, and SS; a L'Oréal Canada For Women in Science Research Excellence Fellowship, with the support of the Canadian Commission for UNESCO Program, to JL; and McGill University Health Centre Research Institute and Department of Medicine Fellowships to NB. We acknowledge infrastructure support and technical assistance from the Breast Cancer Functional Genomics Group, which is partially supported by funds from the Terry Fox New Frontiers Program. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. MP holds the Diane and Sal Guerrera Chair in Cancer Genetics at McGill University.
The authors declare that they have no competing interests.
FP participated in the design and coordination of the study, conducted bioinformatics analyses, and drafted the manuscript. NB supervised the later stages of the experimental work and helped to draft the manuscript. JL contributed to the later stages of the experimental work and helped to draft the manuscript. SS co-conceived of the study and supervised the initial phases of the experimental work. MS conducted immunohistochemistry and laser-capture microdissection experiments. HZ conducted RNA isolation, amplification, microarray, and qRT-PCR experiments. GF assisted with bioinformatics analyses and helped to draft the manuscript. SM provided clinical insights and assisted with obtaining samples. MTH co-supervised the design and coordination of the study, supervised bioinformatics analyses, and helped to draft the manuscript. MP conceived of the study, supervised its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
Nicholas Bertos, Julie Laferrière contributed equally to this work.
Electronic supplementary material
Additional file 1: Supplementary Tables S1 through S3. Supplementary Table S1A, B: Patient information. Supplementary Table S2. KEGG pathways differentially expressed between different categories of vasculature. Supplementary Table S3. List of genes differentially expressed between samples from recurrent and nonrecurrent cancer patients. (DOC 572 KB)
Additional file 2: Supplementary Figures S1 through S6. Supplementary Figure S1 Specific enrichment of the vasculature by laser capture microdissection. Supplementary Figure S2 Heatmap depicting differential gene expression after class discovery within the set of all microdissected endothelial samples and principal component analysis (PCA) of tumor endothelial samples. Supplementary Figure S3 Characteristics of the studies in which other tumor vascular signatures were generated, heatmaps representing expression of antiangiogenic genes and the pericyte marker PDGFRβ in tumor vasculature, and heatmaps of clustering induced by other tumor vascular signatures in the current dataset. Supplementary Figure S4 Images of anti-ACTA2 IHC, as well as of anti-ACTA2/LAMB1 and anti-PECAM1 IHC carried out on consecutive sections. Supplementary Figure S5 qRT-PCR validation for selected genes identified from analyses of microarray data. Supplementary Figure S6 Analyses of correlation between cell-type marker content and expression of genes differentially expressed between the A and B subtypes. (PDF 827 KB)
Additional file 3: Supplementary Materials and Methods. Supplementary Materials and Methods with corresponding references. (DOC 62 KB)
Additional file 4: Supplementary Datasets Sheets 1 and 2. Sheet 1: List of genes differentially expressed between normal and tumor vasculature. Sheet 2A-C: Lists of genes differentially expressed between normal vasculature and A and B tumor vascular-subtype members. (XLS 518 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Pepin, F., Bertos, N., Laferrière, J. et al. Gene-expression profiling of microdissected breast cancer microvasculature identifies distinct tumor vascular subtypes. Breast Cancer Res 14, R120 (2012). https://doi.org/10.1186/bcr3246
- Invasive Ductal Carcinoma
- Tumor Vasculature
- Antiangiogenic Therapy
- Pericyte Coverage
- Vessel Maturity