
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
2 as a Potential Novel Prognostic Marker in Breast Cancer
Authors' Affiliations: 1 Institute of Pathology, University Hospital of the RWTH Aachen, Aachen, Germany; 2 Institute of Pathology, Charité University Hospital; 3 metaGen Pharmaceuticals i.L., Berlin, Germany; 4 Institute of Pathology, University of Regensburg; 5 Central Tumor Registry, Regensburg, Germany; 6 Signature Diagnostics AG, Potsdam, Germany; 7 University Hospital Carl Gustav Carus, Department of Surgery, Dresden, Germany; 8 Department of Gynecology, Friedrich Schiller University, Jena, Germany; 9 Department of Pathology, University of Vermont College of Medicine, Burlington, Vermont; and 10 Institute of Pathology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
Requests for reprints: Edgar Dahl, Institute of Pathology, University Hospital of the RWTH Aachen, Pauwelsstrasse 30, 52047 Aachen, Germany. Phone: 49-241-8088431; Fax: 49-241-8082439; E-mail: edahl{at}ukaachen.de.
| Abstract |
|---|
|
|
|---|
Experimental Design: Twenty-four matched pairs of invasive ductal breast cancer and corresponding benign breast tissue were investigated by a combination of laser microdissection and gene expression profiling. Differential expression of candidate genes was validated by dot blot analysis of cDNA in 50 pairs of matching benign and malignant breast tissue. Cellular expression of candidate genes was further validated by RNA in situ hybridization, quantitative reverse transcription-PCR, and immunohistochemistry using tissue microarray analysis of 272 nonselected breast cancers. Multivariate analysis of factors on overall survival and recurrence-free survival was done.
Results: Fifty-four genes were found to be up-regulated and 78 genes were found to be down-regulated. Dot blot analysis reduced the number of up-regulated genes to 15 candidate genes that showed at least a 2-fold overexpression in >15 of 50 (30%) tumor/normal pairs. We selected phosphatidic acid phosphatase type 2 domain containing 1A (PPAPDC1A) and karyopherin
2 (KPNA2) for further validation. PPAPDC1A and KPNA2 RNA was up-regulated (fold change >2) in 84% and 32% of analyzed tumor/normal pairs, respectively. Nuclear protein expression of KPNA2 was significantly associated with shorter overall survival and recurrence-free survival. Testing various multivariate Cox regression models, KPNA2 expression remained a highly significant, independent and adverse risk factor for overall survival.
Conclusions: Gene expression profiling of laser-microdissected breast cancer tissue revealed novel genes that may represent potential molecular targets for breast cancer therapy and prediction of outcome.
Expression profiling by DNA microarray technology has extensively been applied to identify breast cancerassociated genes (37). Many of such studies, however, are limited by exclusive expression analysis of bulk tissue containing a mixture of tumor, stromal, endothelial, and inflammatory cells.
The aim of this study was to identify potential diagnostic marker genes and molecular targets for breast cancer therapy using a combination of oligonucleotide arrays, multi-tissue Northern blots, and tissue microarray technologies in laser-microdissected pure tissue samples.
| Materials and Methods |
|---|
|
|
|---|
|
Poly(A)+ RNA was isolated from lysed tissue cells by magnetic separation (PolyATract System 1000; Promega, Heidelberg, Germany) according to the specifications of the manufacturer. cDNA synthesis was followed by three cycles of in vitro transcription as previously described (10). In brief, RNA was primed with the Affymetrix T7-oligo-dT promoter-primer combination (5'-GGCCAGTGAATTGTAATACGACTCACTATAGGGAGGCGGT24-3'). TaqMan PCR tested the cDNA of each round of amplification for its integrity and cDNAs of low quality were excluded from further processing. Biotinylated nucleotides were incorporated into the aRNA during the last in vitro transcription. Hybridization and detection of labeled aRNA used the metg001A Affymetrix GeneChip array according to the instructions of the manufacturer.
Chip design and bioinformatics analysis. The custom-designed Affymetrix oligonucleotide array (Gene Chip metg001A) consists of 6,117 probe sets representing 3,950 cDNAs based on the annotation of the probe sets with the GoldenPath assembly11 and has previously been described (11). In short, we selected genes from signal transduction pathways known to be involved in the progression of human tumors (e.g., transforming growth factor ß, RAS, and WNT pathways) and cDNA fragments derived from systematic screening of EST libraries for genes that are differentially expressed in normal and tumor tissues (12).
Preprocessing and chip-level corrections of the Affymetrix metg001A gene chip data has previously been described (13). In brief, Gene Chips were scanned using an Agilent GeneArray Scanner (Agilent Technologies, Palo Alto, CA). Raw intensity values were extracted from the CEL-files. Using the 75% percentile of the perfect match intensities for each probe set (PMQ value) generated a representative expression value. For each probe set, a nonparametric Wilcoxon test was calculated by comparing the intensities of the perfect match and mismatch probes to test the probe sets for presence or absence of an expression signal. To minimize data perturbation caused by experimental variation, a model-fitting algorithm was applied to the PMQ data creating a theoretical reference chip. Expression data of each individual chip were compared with the reference chip. For further analysis, data were transformed into log space (ln).
Probe sets were ranked according to three different variables as described by Kristiansen et al. (11): first, according to the ratio of the median of expression values in tumor samples divided by the median of expression values in normal samples (change fold median); second, according to the percentage of patients who had a tumor/normal fold change exceeding 2; and third, according to the Golub criteria (14). The cross section of the top 5% of genes for each method was selected for further analysis.
Breast cancer profiling array. The breast cancer profiling array (BD Clontech, Heidelberg, Germany) contains 50 pairs of cDNAs generated from matching tumor and normal tissue samples of individual patients, spotted on a nylon membrane (http://www.clontech.com/techinfo/manuals/pdf/pt3424-1.pdf). Hybridizations using 25 ng of a gene-specific 32P-labeled cDNA probe digested from Unigene cDNA clones were done according to the recommendations of the manufacturer. The tumor/normal intensity ratio was calculated using a STORM-860 phosphoimager (Molecular Dynamics, Eugene, OR) and normalized against the background. We defined a candidate gene as up-regulated in breast cancer by this technique if a >2-fold up-regulation was detectable in at least 30% of analyzed tumor/normal pairs (Table 2 ).
|
RNA in situ hybridization. IMAGE consortium plasmid clones were linearized with HindIII and EcoRI for sense and antisense probes, respectively. A 500-bp cDNA fragment (GenBank accession no. N25267) was used for PPAPDC1A. Riboprobes were generated and digoxigenin labeled using the Dig RNA labeling kit (Roche Applied Science, Mannheim, Germany). Sections of paraffin-embedded breast tissues were deparaffinized, rehydrated, and processed according to the instructions of the manufacturer (Roche Applied Science). Hybridized probes were detected using alkaline phosphataseconjugated anti-DIG antibodies and BM Purple as substrate (Roche Applied Science).
Immunohistochemistry. Immunohistochemical studies for the expression of TP53 and HER2 used an avidin-biotin peroxidase method with a 3,3'-diaminobenzidine chromatogen. For karyopherin
2 (KPNA2), the ChemMate detection kit (DAKO, Glostrup, Denmark) was used. After antigen retrieval [microwave oven for 30 minutes (TP53, HER2) or 10 minutes (KPNA2) at 250 W], immunohistochemistry was carried out in a NEXES immunostainer (Ventana, Tucson, AZ) following the instructions of the manufacturer. The following primary antibodies were used: anti-TP53 (mouse monoclonal Bp53-12, Santa Cruz Biotechnology, Inc., Santa Cruz, CA; dilution 1:1,000), anti-KPNA2 (goat polyclonal SC6917, Santa Cruz Biotechnology; dilution 1:200), and anti-HER2 (rabbit polyclonal A0485, DAKO; dilution 1:400). Normal testicular parenchyma was chosen as internal positive control for KPNA2 immunohistochemistry. One surgical pathologist (A.H.) did a blinded evaluation of the slides without knowledge of clinical data. Causes of noninterpretable results included lack of tumor tissue and presence of necrosis or crush artifact. TP53 and KNPA2 positivity was defined as strong nuclear staining in at least 10% cells. HER2 expression was scored according to the DAKO HercepTest (16).
Breast cancer tissue microarray. A tissue microarray was constructed as previously described (17) and contained 289 nonselected formalin-fixed, paraffin-embedded primary breast cancers (stage I -IIIC) together with matched normal breast tissue of each patient. An experienced surgical pathologist (A.H.) evaluated H&E-stained slides of all specimens before construction of the tissue microarray to identify representative tumor areas. Clinical follow-up, provided by the Central Tumor Registry Regensburg, Germany was available for 272 breast cancer patients with a median follow-up period of 79 months (0-148 months). Clinicopathologic variables of breast cancer cases included in the tissue microarray are summarized in Table 3 .
|
| Results |
|---|
|
|
|---|
Up-regulated and down-regulated genes were ranked according to our recently described stringent bioinformatics method (11), which requires that a candidate gene is selected among the top 5% of genes in any of three methods that are regularly used for gene expression analysis. Specifically, the Golub approach, the fold change 2 approach, and the fold change median approach identified 54 up-regulated and 78 down-regulated genes in our samples (Supplementary Figures 2 and 3). Differential expression of all up-regulated candidate genes and two of the down-regulated candidate genes was validated by dot blot analysis on a nylon array containing spotted cDNAs derived from 50 matching pairs of normal and tumor breast tissue (Fig. 1 ). The specificity of each gene probe was controlled on a multi-tissue Northern blot (BD Clontech) containing poly(A)+ RNA samples from 16 human tissues before hybridization to the breast cancer profiling array (Fig. 2A-C ). Candidate genes needed to show an at least 2-fold up-regulation in >30% of analyzed normal/tumor pairs to be at least as differentially expressed as HER2 (Fig. 1). Table 2 summarizes the 15 up-regulated genes fulfilling these criteria and compares their grade of differential expression in DNA array analysis (24 matched pairs), dot blot analysis (50 matched pairs), and reverse transcription-PCR analysis (6 matched pairs). The molecular function of these 15 highly overexpressed genes was annotated according to the Panther classification system (Applied Biosystems). We found that these genes are involved in important biological processes, such as cell cycle control, intracellular signaling, and cell adhesion (Table 2). Our 15 genes reached overexpression values ranging from 32% to 84% whereas HER2 was overexpressed in only 26% of analyzed tissue pairs (Fig. 1). Seven of the 15 genes have never been associated with breast tumorigenesis before.
|
|
|
KPNA2 is a nuclear and cytoplasmic protein also known as importin
1 or Rch1. KPNA2 is thought to connect karyophilic proteins to the nuclear import machinery (24). We found uniform and moderate expression of KPNA2 mRNA in most human normal tissues (Fig. 2C) and prominent expression in testis. The human KPNA2 transcript was determined to be 2.1 kb in size. KPNA overexpression in matching tissues was 67% in the DNA array experiment and 32% on the breast cancer profiling array (Table 2). Because human normal breast tissue was not present on the multi-tissue Northern blot from BD Clontech, we further analyzed KPNA2 expression by real-time PCR in eight microdissected breast normal tissues in comparison with other human normal tissues also present on the multi-tissue Northern blot (Fig. 4
). This analysis showed variable KPNA2 expression in human normal tissue and a mean expression level that was comparable to that found in thymus.
|
|
10% (P = 0.010), and positive HER2 status (P = 0.002) were significantly associated with KPNA2 expression.
Overall survival and recurrence-free survival were compared between KPNA2-negative and KPNA2-positive cases by univariate log-rank statistics (Table 3). Patients with KPNA2-positive tumors (
10%) had an estimated mean overall survival time of 101 months (95% confidence interval, 90-112) compared with 120 months (95% confidence interval, 110-129) in patients with negative KPNA2 staining (P = 0.0047; Fig. 6A
). KPNA2 expression was also associated with shorter recurrence-free survival (P = 0.0013).
|
50 years (P = 0.041), positive lymph node status (P = 0.015), and negative estrogen receptor status (P = 0.021) were significant (Table 4). After reverse selection, the same three variables and KPNA2 immunohistochemistry remained in the model. Nuclear KPNA2 expression was an independent negative risk factor for overall survival (hazard ratio, 2.420; 95% confidence interval 1.200-4.882; P = 0.014). Because of the assumption of proportional hazards, the probability of death was consistently valid during the entire observation period. Finally, multiplicative terms of interaction (Inter 1-7) were considered, representing interactions between nuclear KPNA2 expression and dichotomous covariables. In the global model, only the interaction between lymph node status and KPNA2 immunohistochemistry was significant (P = 0.031). After reverse selection, a model containing age
50 years (P = 0.023), negative estrogen receptor status (P = 0.003), and the interaction between lymph node status and KPNA2 immunohistochemistry (Inter 3; hazard ratio, 4.558; 95% confidence interval, 2.397-8.664; P < 0.001) was found.
|
| Discussion |
|---|
|
|
|---|
3,950 cDNA fragments) found 15 genes (Table 2) that are more up-regulated than HER2. Differential expression was evaluated by hybridization of gene-specific probes to a dot blot array containing 50 cDNA pairs derived from breast tumor and matching normal breast tissue (Breast Cancer Profiling Array, BD Clontech). Whereas HER2 RNA was up-regulated at least 2-fold in 26% of matching tumor/normal pairs, the 15 genes reached scores between 32% (KPNA2) and 84% (PPAPDC1A). Several recently discovered breast cancerassociated genes were included in the 15 gene list [e.g., aurora kinase A (AURKA; ref. 25) and calgranulin B (S100A9; ref. 26)]. However, seven genes have not been implied in mammary carcinogenesis thus far and may increase the pool of potentially useful molecules for breast cancer diagnosis and therapy.
In 1999, Blobel and colleagues showed that proteins of the karyopherin
and ß families play a central role in nucleocytoplasmic transport (reviewed in ref. 27). KPNA2 is an adaptor protein within the classic nuclear protein import machinery which mediates the import of signaling factors in the nucleus and the export of response molecules to the cytoplasm (24). The classic nuclear protein import requires a nuclear localization signal within the cargo protein to be imported. This signal is recognized by the adaptor protein karyopherin/importin
. In blood lymphocytes, activation of cellular signaling was associated with strongly increased KPNA2 expression and redistribution of the KPNA2 protein from the cytoplasm to both the nucleus and the plasma membrane (28). Functional studies of BRCA1 by Thakur et al. have provided conclusive evidence that BRCA1 is a nuclear protein transported into the nucleus by the karyopherin pathway (29), possibly linking KPNA2 expression with breast carcinogenesis. The present study observed strong nuclear staining of KPNA2 in breast tumor cells compared with a weak or absent staining in normal breast tissue, consistent with the hypothesis that increased nuclear signaling is present in a high percentage of human breast tumors. Lack of KPNA2 expression in normal breast tissue and overexpression in the majority of carcinomas exposes KPNA2 immunohistochemistry as a potential diagnostic marker. Investigation of KPNA2 expression in preneoplastic lesions of the breast is mandatory to establish the time point of KPNA2 up-regulation in the multistep process of mammary carcinogenesis.
Many studies have focused on the potential of gene expression profiles to predict the clinical outcome of breast cancer (46, 3032). However, results of former studies are limited by exclusive analysis of bulk tissue containing a mixture of tumor, stromal, endothelial, and inflammatory cells. Some of the genes in the signatures presented to date seem to be derived from nonepithelial components of the tumor (31), suggesting that stromal elements represent an important contributing factor to the metastatic phenotype.
Survival differences were also noted between the different subtypes of breast tumors as defined by expression patterns (32, 33). Five distinct gene expression patterns were distinguished (3234), including one basal-like, one HER2-overexpressing, two luminal-like, and one normal breast tissuelike subgroup. Cluster analysis of two published independent data sets, representing different patient cohorts from different laboratories, uncovered the same breast cancer subtypes (32). The patients with basal-like and HER2-subtypes were associated with the shortest survival. This strongly supported the idea that many of these breast tumor subtypes represent biologically distinct disease entities with different clinical outcome.
Supervised clustering with patient survival as the supervising variable by Sorlie et al. (33) resulted in a final list of 264 cDNA clones. This 264-clone set was then used for a hierarchical clustering analysis. The resulting diagram showed that almost all of the 264 cDNA clones selected fell into the three main gene expression clusters: the luminal/estrogen receptor positive cluster, the basal epithelial cluster, and the proliferation cluster. Interestingly, KPNA2 was strongly expressed in tumors of the proliferation cluster (i.e., in the group of tumors with the shortest overall and recurrence-free survival). Notably, Dressman et al. (35) defined molecular signatures that correlate with response to neoadjuvant chemotherapy. Karyopherin
6 (KPNA6) was among the genes that predicted clinical response of stage IIB/III breast cancer patients (35). In our study, KPNA2 expression (
10%) was significantly associated with shorter recurrence free survival (P = 0.0013). Due to missing retrospective data on chemotherapy and radiation therapy, no conclusions about therapy stratification of breast cancer patients based on KPNA2 expression can be drawn. A prospective clinical trial is currently planned together with the Department of Gynecology, University of Regensburg, Regensburg, Germany, to address this issue.
A remarkable feature of the expression signatures identified in previous studies is that they usually involve fewer than 100 genes (5, 30), in one instance even only 17 genes (31). However, the incomplete overlap between the different sets of defined genes is confusing (6, 32). Comparison of our 54 up-regulated and 78 down-regulated genes with the gene signatures of van de Vijver et al. (4) and Wang et al. (7) is difficult because of differences in patients and techniques. No overlap with the study of van de Vijver et al. (ref. 4; 70 genes) and only a single match (KPNA2) to the study by Wang et al. (ref. 7; 76 genes) were found. However, these studies compared breast tumors from long-term and short-term breast cancer survivors whereas our study compared normal breast tissue and invasive breast cancer. Therefore, a large overlap, which was even not detectable between the two closely related studies mentioned above, was not expected. Of note is that KPNA2, the gene we have focused on in our study, was also found by Wang et al. (7).
Sortiriou et al. (36) have examined whether histologic grade was associated with gene expression profiles of breast cancer and whether such profiles could be used to improve histologic grading; expression profiles between histologic grade 1 and histologic grade 3 tumors were compared. Interestingly, KPNA2 was among the top overexpressed genes that were associated with histologic grade (36). Sortiriou et al. concluded that grade based on gene expression may reclassify grade 2 tumors into two groups with high and low risk of recurrence, improving the accuracy of histologic tumor grading and thus its prognostic value. In our tissue microarray study, KPNA2 expression and high tumor grade coassociated significantly (P < 0.001). However, in the subgroup of grade 2 tumors, no difference in tumor-related survival (P = 0.1807) and tumor recurrence (P = 0.3156) on KPNA2 expression could be observed. Previous studies have shown that mutations of the TP53 gene predict poor prognosis (reviewed in ref. 37) and are associated with poor response to systemic chemotherapy (38). Overexpression of the HER2 protein is a well-known prognostic factor associated with poor survival in breast cancer, which also was found for the HER2-positive group defined by Sorlie et al. (33). In our study, high rates of KNPA2 expression were significantly associated with positive TP53 and HER2 immunohistochemistry and a high proliferation index (Table 3). KPNA2 seems to be characteristic of the basal-like subtype of breast cancers, possibly representing a different clinical entity of breast tumors, which is associated with shorter survival times and a high frequency of TP53 mutations.
In summary, novel genes with clinical utility to select patients more likely to develop aggressive disease have been identified using a combination of laser microdissection of matched tumor/normal pairs and array expression analysis. We illustrated an independent negative correlation between KPNA2 expression in the primary tumor and overall survival in node-positive breast cancer patients. KPNA2 expression may represent a potential diagnostic marker to predict differential clinical responses to treatment and disease outcome. Prospective studies are currently conducted to validate our findings.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
P.J. Wild and A. Rosenthal share the senior authorship for this work.
Received 11/ 1/05; revised 3/23/06; accepted 4/27/06.
| References |
|---|
|
|
|---|
nuclear localization signal (NLS)-binding groove: functional requirements for NLS binding. J Biol Chem 2003;278:4194753.
(Rch1) at the plasma membrane and subcellular redistribution during lymphocyte activation. Chromosoma 2003;112:8795.[Medline]This article has been cited by other articles:
![]() |
P. Kalab and R. Heald The RanGTP gradient - a GPS for the mitotic spindle J. Cell Sci., May 15, 2008; 121(10): 1577 - 1586. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Douard, X.-L. Cui, P. Soteropoulos, and R. P. Ferraris Dexamethasone Sensitizes the Neonatal Intestine to Fructose Induction of Intestinal Fructose Transporter (Slc2A5) Function Endocrinology, January 1, 2008; 149(1): 409 - 423. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. J. Terry, E. B. Shows, and S. R. Wente Crossing the Nuclear Envelope: Hierarchical Regulation of Nucleocytoplasmic Transport Science, November 30, 2007; 318(5855): 1412 - 1416. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |