Omics-based profiling of carcinoma of the breast and matched regional lymph node metastasis

Proteomics of REPLICANT perfusate detects changes in the metastatic lymph node microenvironment

In breast cancer (BC), detecting low volumes of axillary lymph node (ALN) metastasis pre-operatively is difficult and novel biomarkers are needed. We recently showed that patient-derived ALNs can be sustained ex-vivo using normothermic perfusion. We now compare reactive (tumour-free; n = 5) and macrometastatic (containing tumour deposits >2 mm; n = 4) ALNs by combining whole section multiplex immunofluorescence with TMT-labelled LC-MS/MS of the circulating perfusate. Macrometastases contained significantly fewer B cells and T cells (CD4⁺/CD8⁺/regulatory) than reactive nodes (p = 0.02). Similarly, pathway analysis of the perfusate proteome (119/1453 proteins significantly differentially expressed) showed that immune function was diminished in macrometastases in favour of ‘extracellular matrix degradation’; only ‘neutrophil degranulation’ was preserved. Qualitative comparison of the perfusate proteome to that of node-positive pancreatic and prostatic adenocarcinoma also highlighted ‘neutrophil degranulation’ as a contributing factor to nodal metastasis. Thus, metastasis-induced changes in the REPLICANT perfusate proteome are detectable, and could facilitate biomarker discovery.

Molecular Concordance Between Primary Breast Cancer and Matched Metastases

Clinical management of breast cancer is increasingly personalized and based on molecular profiling. Often, primary tumors are used as proxies for systemic disease at the time of recurrence. However, recent studies have revealed substantial discordances between primary tumors and metastases, both with respect to traditional clinical treatment targets and on the genomic and transcriptomic level. With the increasing use of molecularly targeted therapy, discordance of actionable molecular targets between primary tumors and recurrences can result in nonoptimal treatment or unnecessary side effects. The purpose of this review is to illuminate the extent of cancer genome evolution through disease progression and the degree of molecular concordance between primary breast cancers and matched metastases. We present an overview of the most prominent studies investigating the expression of endocrine receptors, transcriptomics, and genome aberrations in primary tumors and metastases. In conclusion, biopsy of metastatic lesions at recurrence of breast cancer is encouraged to provide optimal treatment of the disease. Furthermore, molecular profiling of metastatic tissue provides invaluable mechanistic insight into the biology underlying metastatic progression and has the potential to identify novel, potentially druggable, drivers of progression.

Clonal expansion and linear genome evolution through breast cancer progression from pre-invasive stages to asynchronous metastasis

Evolution of the breast cancer genome from pre-invasive stages to asynchronous metastasis is complex and mostly unexplored, but highly demanded as it may provide novel markers for and mechanistic insights in cancer progression. The increasing use of personalized therapy of breast cancer necessitates knowledge of the degree of genomic concordance between different steps of malignant progression as primary tumors often are used as surrogates of systemic disease. Based on exome sequencing we performed copy number profiling and point mutation detection on successive steps of breast cancer progression from one breast cancer patient, including two different regions of Ductal Carcinoma In Situ (DCIS), primary tumor and an asynchronous metastasis. We identify a remarkable landscape of somatic mutations, retained throughout breast cancer progression and with new mutational events emerging at each step. Our data, contrary to the proposed model of early dissemination of metastatic cells and parallel progression of primary tumors and metastases, provide evidence of linear progression of breast cancer with relatively late dissemination from the primary tumor. The genomic discordance between the different stages of tumor evolution in this patient emphasizes the importance of molecular profiling of metastatic tissue directing molecularly targeted therapy at recurrence.

Changes in glycoprotein expression between primary breast tumour and synchronous lymph node metastases or asynchronous distant metastases

Background

Breast cancer is a very heterogeneous disease and some patients are cured by the surgical removal of the primary tumour whilst other patients suffer from metastasis and spreading of the disease, despite adjuvant therapy. A number of prognostic and treatment predictive factors have been identified such as tumour size, oestrogen (ER) and progesterone (PgR) receptor status, human epidermal growth factor receptor type 2 (HER2) status, histological grade, Ki67 and age. Lymph node involvement is also assessed during surgery to determine if the tumour has spread which requires dissection of the axilla and adjuvant treatment. The prognostic and treatment predictive factors assessing the nature of the tumour are all routinely based on the status of the primary tumour.

Results

We have analysed a unique tumour set of fourteen primary breast cancer tumours with matched synchronous axillary lymph node metastases and a set of nine primary tumours with, later developed, matched distant metastases from different sites in the body. We used a pairwise tumour analysis (from the same individual) since the difference between the same tumour-type in different patients was greater. Glycopeptide capture was used in this study to selectively isolate and quantify N-linked glycopeptides from tumours mixtures and the captured glycopeptides were subjected to label-free quantitative tandem mass spectrometry analysis. Differentially expressed proteins between primary tumours and matched lymph node metastasis and distant metastasis were identified. Two of the top hits, ATPIF1 and tubulin β-chain were validated by immunohistochemistry to be differentially regulated.

Conclusions

We show that the expression of a large number of glycosylated proteins change between primary tumours and matched lymph node metastases and distant metastases, confirming that cancer cells undergo a molecular transformation during the spread to a secondary site. The proteins are part of important pathways such as cell adhesion, migration pathways and immune response giving insight into molecular changes needed for the tumour to spread. The large difference between primary tumours and lymph node and distant metastases also suggest that treatment should be based on the phenotype of the lymph node and distant metastases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12014-015-9084-7) contains supplementary material, which is available to authorized users.

Pathway-centric integrative analysis identifies RRM2 as a prognostic marker in breast cancer associated with poor survival and tamoxifen resistance

Breast cancer (BCa) molecular subtypes include luminal A, luminal B, normal-like, HER-2-enriched, and basal-like tumors, among which luminal B and basal-like cancers are highly aggressive. Biochemical pathways associated with patient survival or treatment response in these more aggressive subtypes are not well understood. With the limited availability of pathologically verified clinical specimens, cell line models are routinely used for pathway-centric studies. We measured the metabolome of luminal and basal-like BCa cell lines using mass spectrometry, linked metabolites to biochemical pathways using Gene Set Analysis, and developed a novel rank-based method to select pathways on the basis of their enrichment in patient-derived omics data sets and prognostic relevance. Key mediators of the pathway were then characterized for their role in disease progression. Pyrimidine metabolism was altered in luminal versus basal BCa, whereas the combined expression of its associated genes or expression of one key gene, ribonucleotide reductase subunit M2 (RRM2) alone, associated significantly with decreased survival across all BCa subtypes, as well as in luminal patients resistant to tamoxifen. Increased RRM2 expression in tamoxifen-resistant patients was verified using tissue microarrays, whereas the metabolic products of RRM2 were higher in tamoxifen-resistant cells and in xenograft tumors. Both genetic and pharmacological inhibition of this key enzyme in tamoxifen-resistant cells significantly decreased proliferation, reduced expression of cell cycle genes, and sensitized the cells to tamoxifen treatment. Our study suggests for evaluating RRM2-associated metabolites as noninvasive markers for tamoxifen resistance and its pharmacological inhibition as a novel approach to overcome tamoxifen resistance in BCa.

Increased copy number of the DLX4 homeobox gene in breast axillary lymph node metastasis

DLX4 is a homeobox gene strongly implicated in breast tumor progression and invasion. Our main objective was to determine the DLX4 copy number status in sentinel lymph node (SLN) metastasis to assess its involvement in the initial stages of the axillary metastatic process. A total of 37 paired samples of SLN metastasis and primary breast tumors (PBT) were evaluated by fluorescence in situ hybridization, quantitative polymerase chain reaction and array comparative genomic hybridization assays. DLX4 increased copy number was observed in 21.6% of the PBT and 24.3% of the SLN metastasis; regression analysis demonstrated that the DLX4 alterations observed in the SLN metastasis were dependent on the ones in the PBT, indicating that they occur in the primary tumor cell populations and are maintained in the early axillary metastatic site. In addition, regression analysis demonstrated that DLX4 alterations (and other DLX and HOXB family members) occurred independently of the ones in the HER2/NEU gene, the main amplification driver on the 17q region. Additional studies evaluating DLX4 copy number in non-SLN axillary lymph nodes and/or distant breast cancer metastasis are necessary to determine if these alterations are carried on and maintained during more advanced stages of tumor progression and if could be used as a predictive marker for axillary involvement.

Biological resonance for cancer metastasis, a new hypothesis based on comparisons between primary cancers and metastases

Many hypotheses have been proposed to try to explain cancer metastasis. However, they seem to be contradictory and have some limitations. Comparisons of primary tumors and matched metastases provide new insight into metastasis. The results show high concordances and minor differences at multiple scales from organic level to molecular level. The concordances reflect the commonality between primary cancer and metastasis, and also mean that metastatic cancer cells derived from primary cancer are quite conservative in distant sites. The differences reflect variation that cancer cells must acquire new traits to adapt to foreign milieu during the course of evolving into a new tumor in second organs. These comparisons also provided new information on understanding mechanism of vascular metastasis, organ-specific metastasis, and tumor dormancy. The collective results suggest a new hypothesis, biological resonance (bio-resonance) model. The hypothesis has two aspects. One is that primary cancer and matched metastasis have a common progenitor. The other is that both ancestors of primary cancer cells and metastatic cancer cells are under similar microenvironments and receive similar or same signals. When their interactions reach a status similar to primary cancer, metastasis will occur. Compared with previous hypotheses, the bio-resonance hypothesis seems to be more applicable for cancer metastasis to explain how, when and where metastasis occurs. Thus, it has important implications for individual prediction, prevention and treatment of cancer metastasis.

Differential proteomic analysis of a human breast tumor and its matched bone metastasis identifies cell membrane and extracellular proteins associated with bone metastasis

The classical fate of metastasizing breast cancer cells is to seed and form secondary colonies in bones. The molecules closely associated with these processes are predominantly present at the cell surface and in the extracellular space, establishing the first contacts with the target tissue. In this study, we had the rare opportunity to analyze a bone metastatic lesion and its corresponding breast primary tumor obtained simultaneously from the same patient. Using mass spectrometry, we undertook a proteomic study on cell surface and extracellular protein-enriched material. We provide a repertoire of significantly modulated proteins, some with yet unknown roles in the bone metastatic process as well as proteins notably involved in cancer cell invasiveness and in bone metabolism. The comparison of these clinical data with those previously obtained using a human osteotropic breast cancer cell line highlighted an overlapping group of proteins. Certain differentially expressed proteins are validated in the present study using immunohistochemistry on a retrospective collection of breast tumors and matched bone metastases. Our exclusive set of selected proteins supports the setup of further investigations on both clinical samples and experimental bone metastasis models that will help to reveal the finely coordinated expression of proteins that favor the development of metastases in the bone microenvironment.

Tumor heterogeneity in neoplasms of breast, colon, and skin

BACKGROUND

Different cell subpopulations in a single tumor may show diverse capacities for growth, differentiation, metastasis formation, and sensitivity to treatments. Thus, heterogeneity is an important feature of tumors. However, due to limitations in experimental and analytical techniques, tumor heterogeneity has rarely been studied in detail.

PRESENTATION OF THE HYPOTHESIS

Different tumor types have different heterogeneity patterns, thus heterogeneity could be a characteristic feature of a particular tumor type.

TESTING THE HYPOTHESIS

We applied our previously published mathematical heterogeneity model to decipher tumor heterogeneity through the analysis of genetic copy number aberrations revealed by array CGH data for tumors of three different tissues: breast, colon, and skin. The model estimates the number of subpopulations present in each tumor. The analysis confirms that different tumor types have different heterogeneity patterns. Computationally derived genomic copy number profiles from each subpopulation have also been analyzed and discussed with reference to the multiple hypothetical relationships between subpopulations in origin-related samples.

IMPLICATIONS OF THE HYPOTHESIS

Our observations imply that tumor heterogeneity could be seen as an independent parameter for determining the characteristics of tumors. In the context of more comprehensive usage of array CGH or genome sequencing in a clinical setting our study provides a new way to realize the full potential of tumor genetic analysis.

Variations in stromal signatures in breast and colorectal cancer metastases

The tumour microenvironment (TME) plays an important role in tumour survival and growth, but little is known about the degree of preservation between different stromal response patterns found in primary tumours and their metastases. We have previously identified gene expression profiles for two distinct stromal signatures in breast carcinoma of fibroblast (aka DTF) and macrophage (aka CSF1) response and found them to be correlated with clinicopathological features, including outcome. In this study, we compare the DTF fibroblast and CSF1 macrophage stromal response patterns in primary breast and colorectal cancers to their matched lymph node metastases. In both breast and colorectal cancer, there was a significant positive correlation between the CSF1 macrophage signature in the primary tumours and the matched lymph node metastases, as assessed by immunohistochemical markers. No such correlation was observed for the DTF fibroblast signature. A similar result was seen in independent analysis of two published gene expression microarray datasets. The variations of these stromal reaction patterns from the primary to the metastasis shed light on the relationship between the neoplastic cells and the non-neoplastic cells in the TME. The preservation of the CSF1 macrophage response pattern in metastases lends support to targeting the CSF1 pathway in cancer.

Frequent genetic differences between matched primary and metastatic breast cancer provide an approach to identification of biomarkers for disease progression

Breast cancer is a major cause of morbidity and mortality in women and its metastatic spread is the principal reason behind the fatal outcome. Metastasis-related research of breast cancer is however underdeveloped when compared with the abundant literature on primary tumors. We applied an unexplored approach comparing at high resolution the genomic profiles of primary tumors and synchronous axillary lymph node metastases from 13 patients with breast cancer. Overall, primary tumors displayed 20% higher number of aberrations than metastases. In all but two patients, we detected in total 157 statistically significant differences between primary lesions and matched metastases. We further observed differences that can be linked to metastatic disease and there was also an overlapping pattern of changes between different patients. Many of the differences described here have been previously linked to poor patient survival, suggesting that this is a viable approach toward finding biomarkers for disease progression and definition of new targets useful for development of anticancer drugs. Frequent genetic differences between primary tumors and metastases in breast cancer also question, at least to some extent, the role of primary tumors as a surrogate subject of study for the systemic disease.

DNA copy number aberrations in breast cancer by array comparative genomic hybridization

Array comparative genomic hybridization (CGH) has been popularly used for analyzing DNA copy number variations in diseases like cancer. In this study, we investigated 82 sporadic samples from 49 breast cancer patients using 1-Mb resolution bacterial artificial chromosome CGH arrays. A number of highly frequent genomic aberrations were discovered, which may act as “drivers” of tumor progression. Meanwhile, the genomic profiles of four “normal” breast tissue samples taken at least 2 cm away from the primary tumor sites were also found to have some genomic aberrations that recurred with high frequency in the primary tumors, which may have important implications for clinical therapy. Additionally, we performed class comparison and class prediction for various clinicopathological parameters, and a list of characteristic genomic aberrations associated with different clinicopathological phenotypes was compiled. Our study provides clues for further investigations of the underlying mechanisms of breast carcinogenesis.

Molecular markers of breast axillary lymph node metastasis

In breast cancer, axillary lymph node status is one of the most important prognostic variables and a crucial component to the staging system. Several clinico-histopathological parameters are considered to be strong predictors of metastasis; however, they fail to accurately classify breast tumors according to their clinical behavior and to predict which patients will have disease recurrence. Methods based on genome-wide microarray analyses have been used to identify molecular markers with respect to the development of axillary lymph node metastasis. Most of these markers can be detected in the primary tumors, which can potentially lead to the ability to identify patients at the time of diagnosis who are at high risk for lymph node metastasis, allowing for early intervention and more suitable adjuvant treatments.

Estimation of tumor heterogeneity using CGH array data

Background

Array-based comparative genomic hybridization (CGH) is a commonly-used approach to detect DNA copy number variation in whole genome-wide screens. Several statistical methods have been proposed to define genomic segments with different copy numbers in cancer tumors. However, most tumors are heterogeneous and show variation in DNA copy numbers across tumor cells. The challenge is to reveal the copy number profiles of the subpopulations in a tumor and to estimate the percentage of each subpopulation.

Results

We describe a relation between experimental data and exact DNA copy number and develop a statistical method to reveal the heterogeneity of tumors containing a mixture of different-stage cells. Furthermore, we validate the method on simulated data and apply the method to 29 pairs of breast primary tumors and their matched lymph node metastases.

Conclusion

We demonstrate a new method for CGH array analysis that allows a tumor sample to be classified according to its heterogeneity. The method gives an interpretable series of copy number profiles, one for each major subpopulation in a tumor. The profiles facilitate identification of copy number alterations in cancer development.

A single lysis solution for the analysis of tissue samples by different proteomic technologies

Cancer, being a major healthcare concern worldwide, is one of the main targets for the application of emerging proteomic technologies and these tools promise to revolutionize the way cancer will be diagnosed and treated in the near future. Today, as a result of the unprecedented advances that have taken place in molecular biology, cell biology and genomics there is a pressing need to accelerate the translation of basic discoveries into clinical applications. This need, compounded by mounting evidence that cellular model systems are unable to fully recapitulate all biological aspects of human dissease, is driving scientists to increasingly use clinically relevant samples for biomarker and target discovery. Tissues are heterogeneous and as a result optimization of sample preparation is critical for generating accurate, representative, and highly reproducible quantitative data. Although a large number of protocols for preparation of tissue lysates has been published, so far no single recipe is able to provide a "one-size fits all" solubilization procedure that can be used to analyse the same lysate using different proteomics technologies. Here we present evidence showing that cell lysis buffer 1 (CLB1), a lysis solution commercialized by Zeptosens [a division of Bayer (Schweiz) AG], provides excellent sample solubilization and very high 2D PAGE protein resolution both when using carrier ampholytes and immobilized pH gradient strips. Moreover, this buffer can also be used for array-based proteomics (reverse-phase lysate arrays or direct antibody arrays), allowing the direct comparison of qualitative and quantitative data yielded by these technologies when applied to the same samples. The usefulness of the CLB1 solution for gel-based proteomics was further established by 2D PAGE analysis of a number of technically demanding specimens such as breast carcinoma core needle biopsies and problematic tissues such as brain cortex, cerebellum, skeletal muscle, kidney cortex and tongue. This solution when combined with a specific sample preparation technique - cryostat sectioning of frozen specimens - simplifies tissue sample preparation and solves most of the difficulties associated with the integration of data generated by different proteomic technologies.