Evaluation of the ability of adjuvant tamoxifen-benefit gene signatures to predict outcome of hormone-naive estrogen receptor-positive breast cancer patients treated with tamoxifen in the advanced setting

Differential expression of NAT1 pharmacogene in hormone receptor positive vs. negative female breast tumors may affect drug treatment

Studies have reported overexpression of NAT1 gene for xenobiotic metabolizing arylamine N -acetyltransferase type 1 in estrogen receptor positive breast tumors, and this association has been linked to patient chemoresistance and response to tamoxifen. We probed the expression of NAT1 , using quantitative reverse transcription PCR to screen clinically characterized breast cancer tissue cDNA arrays. Primers detecting all NAT1 alternative transcripts were used, and the protocol and results are reported according to consensus guidelines. The clinical information about 166 tumor samples screened is provided, including tumor stage, estrogen and progesterone receptor status and HER2 expression. NAT1 was found to be significantly ( P  < 0.001) upregulated in hormone receptor positive vs. negative tumors. No correlation was apparent between NAT1 and tumor stage or HER2 expression. Our findings demonstrate a strong correlation between the expression of NAT1 and steroid hormone receptors in breast tumors, supporting its possible utility as a pharmacogenetic biomarker or drug target. Of the two polymorphic NAT genes, NAT1 is the one primarily expressed in breast tissue, and is subjected to regulation by two differential promoters and more than one polyadenylation signal. Hormonal factors may enhance NAT1 gene expression at the transcriptional or epigenetic level, and tamoxifen has additionally been shown to inhibit NAT1 enzymatic activity. The outcome of tamoxifen treatment is also more favorable in patients with NAT1 overexpressing tumors. The study adds to the growing body of evidence implicating NAT1 in breast cancer and its pharmacological treatment.

Comprehensive immunohistochemical analysis of RET, BCAR1, and BCAR3 expression in patients with Luminal A and B breast cancer subtypes

PURPOSE

Breast cancer (BC) is considered a heterogeneous disease composed of distinct subtypes with diverse clinical outcomes. Luminal subtype tumors have the best prognosis, and patients benefit from endocrine therapy. However, resistance to endocrine therapies in BC is an obstacle to successful treatment, and novel biomarkers are needed to understand and overcome this mechanism. The RET, BCAR1, and BCAR3 genes may be associated with BC progression and endocrine resistance.

METHODS

Aiming to evaluate the expression profile and prognostic value of RET, BCAR1, and BCAR3, we performed immunohistochemistry on tissue microarrays (TMAs) containing a cohort of 361 Luminal subtype BC.

RESULTS

Low expression levels of these three proteins were predominantly observed. BCAR1 expression was correlated with nuclear grade (p = 0.057), and BCAR3 expression was correlated with lymph node status (p = 0.011) and response to hormonal therapy (p = 0.021). Further, low expression of either BCAR1 or BCAR3 was significantly associated with poor prognosis (p = 0.005; p = 0.042). Pairwise analysis showed that patients with tumors with low BCAR1/low BCAR3 expression had a poorer overall survival (p = 0.013), and the low BCAR3 expression had the worst prognosis with RET high expression stratifying these patients into two different groups. Regarding the response to hormonal therapy, non-responder patients presented lower expression of RET in comparison to the responder group (p = 0.035). Additionally, the low BCAR1 expression patients had poorer outcomes than BCAR1 high (p = 0.015).

CONCLUSION

Our findings suggest RET, BCAR1, and BCAR3 as potential candidate markers for endocrine therapy resistance in Luminal BC.

BCAR3 promotes head and neck cancer growth and is associated with poor prognosis

Breast cancer anti-estrogen resistance protein 3 (BCAR3) is involved in anti-estrogen resistance and other important aspects of breast cancer. However, the role of BCAR3 in other solid tumors remains unclear. The relationship between the clinicopathologic characteristics of head and neck squamous cell carcinoma (HNSCC) patients and BCAR3 was analyzed using the Wilcoxon’s signed-rank test and logistic regression. The association between BCAR3 expression and clinicopathologic features and survival was analyzed using Cox regression and the Kaplan–Meier method. In vivo and in vitro assays were performed to validate the effect of BCAR3 on HNSCC growth. BCAR3-related mRNAs were determined by calculating the Pearson’s correlation coefficient based on The Cancer Genome Atlas (TCGA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, and gene set enrichment analysis (GSEA) were used to predict the potential functions of BCAR3. BCAR3 expression is overexpressed in HNSCC and was shown to be associated with perineural invasion (PNI) and poor survival. BCAR3 silencing significantly attenuated the proliferation of HNSCC cells, whereas BCAR3 depletion inhibited tumor growth in vitro. GO and KEGG functional enrichment analyses, and GSEA showed that BCAR3 expression in HNSCC was associated with biological processes, such as cell adhesion, actin binding, cadherin binding, and angiogenesis. BCAR3, which promotes HNSCC growth, is associated with perineural invasion and may be a potential molecular prognostic marker of poor survival in HNSCC.

Circular RNA in Chemonaive Lymph Node Negative Colon Cancer Patients

Circular RNAs (circRNAs) appear important in tumor progression of colon cancer (CC). We identified an extensive catalog of circRNAs in 181 chemonaive stage I/II colon tumors, who underwent curative surgery between 2007 and 2014. We identified circRNAs from RNAseq data, investigated common biology related to circRNA expression, and studied the association between circRNAs and relapse status, tumor stage, consensus molecular subtypes (CMS), tumor localization and microsatellite instability (MSI). We identified 2606 unique circRNAs. 277 circRNAs (derived from 260 genes) were repeatedly occurring in at least 20 patients of which 153 showed a poor or even negative (R < 0.3) correlation with the expression level of their linear gene. The circular junctions for circSATB2, circFGD6, circKMT2C and circPLEKHM3 were validated by Sanger sequencing. Multiple correspondence analysis showed that circRNAs were often co-expressed and that high diversity in circRNAs was associated with favorable disease-free survival (DFS), which was confirmed by Cox regression analysis (Hazard Ratio (HR) 0.60, 95% CI 0.38-0.97, p = 0.036). Considering individual circRNAs, absence of circMGA was significantly associated with relapse, whereas circSATB2, circNAB1, and circCEP192 were associated with both MSI and CMS. This study represents a showcase of the potential clinical utility of circRNAs for prognostic stratification in patients with stage I-II colon cancer and demonstrated that high diversity in circRNAs is associated with favorable DFS.

Promyelocytic Leukemia (PML) gene regulation: implication towards curbing oncogenesis

Dysregulation of PML, a significant tumor suppressor is linked with cancers of different histological origins, with a decreased expression observed with a higher tumor grade. This necessitates studying the mechanisms to maintain a stable expression of PML. However much less is known about the transcriptional regulation of PML, more so in the context of breast carcinoma. ERβ has emerged as a critical factor in understanding breast cancer, especially since a huge proportion of breast cancers are ERα⁻ and thus insensitive to tamoxifen therapy. This study aims to uncover an unidentified mechanism of PML gene regulation and its stabilization in breast cancer via ERβ signalling and the impact on cellular apoptosis. We found that clinical expression of PML positively correlates with that of ERβ both in normal and breast carcinoma samples and inversely correlates with markers of cellular proliferation, hinting towards a possible mechanistic interdependence. Both mRNA and protein expression of PML were increased in response to ERβ overexpression on multiple human breast cancer cell lines. Mechanistically, luciferase reporter assays and chromatin-immunoprecipitation assays demonstrated that ERβ can interact with the PML promoter via ERE and AP1 sites to enhance its transcription. ERβ induced stable PML expression causes a decline of its target protein Survivin and simultaneously provides a stable docking platform leading to stabilisation of its target Foxo3a, further causing transcriptional upregulation of pro-apoptotic factors p21 and p27. Immunohistochemical analyses of cancer and normal breast tissues and functional assays conducted corroborated the findings. Collectively, our study identifies ERβ signalling as a novel mechanism for PML gene regulation in ERα⁻ breast cancer. It also reveals bi-directional downstream effect in which ‘ERβ-PML-(Foxo3a/Survivin)’ network acts as a therapeutic axis by suppressing cellular survival and promoting cellular apoptosis in breast carcinoma.

The DNA deaminase APOBEC3B interacts with the cell-cycle protein CDK4 and disrupts CDK4-mediated nuclear import of Cyclin D1

Apolipoprotein B mRNA editing enzyme catalytic subunit-like protein 3B (APOBEC3B or A3B), as other APOBEC3 members, is a single-stranded (ss)DNA cytosine deaminase with antiviral activity. A3B is also overexpressed in multiple tumor types, such as carcinomas of the bladder, cervix, lung, head/neck, and breast. A3B generates both dispersed and clustered C-to-T and C-to-G mutations in intrinsically preferred trinucleotide motifs (TCA/TCG/TCT). A3B-catalyzed mutations are likely to promote tumor evolution and cancer progression and, as such, are associated with poor clinical outcomes. However, little is known about cellular processes that regulate A3B. Here, we used a proteomics approach involving affinity purification coupled to MS with human 293T cells to identify cellular proteins that interact with A3B. This approach revealed a specific interaction with cyclin-dependent kinase 4 (CDK4). We validated and mapped this interaction by co-immunoprecipitation experiments. Functional studies and immunofluorescence microscopy experiments in multiple cell lines revealed that A3B is not a substrate for CDK4-Cyclin D1 phosphorylation nor is its deaminase activity modulated. Instead, we found that A3B is capable of disrupting the CDK4-dependent nuclear import of Cyclin D1. We propose that this interaction may favor a more potent antiviral response and simultaneously facilitate cancer mutagenesis.

Population variability of rhesus macaque (Macaca mulatta) NAT1 gene for arylamine N-acetyltransferase 1: Functional effects and comparison with human

Human NAT1 gene for N-acetyltransferase 1 modulates xenobiotic metabolism of arylamine drugs and mutagens. Beyond pharmacogenetics, NAT1 is also relevant to breast cancer. The population history of human NAT1 suggests evolution through purifying selection, but it is unclear whether this pattern is evident in other primate lineages where population studies are scarce. We report NAT1 polymorphism in 25 rhesus macaques (Macaca mulatta) and describe the haplotypic and functional characteristics of 12 variants. Seven non-synonymous single nucleotide variations (SNVs) were identified and experimentally demonstrated to compromise enzyme function, mainly through destabilization of NAT1 protein and consequent activity loss. One non-synonymous SNV (c.560G > A, p.Arg187Gln) has also been characterized for human NAT1 with similar effects. Population haplotypic and functional variability of rhesus NAT1 was considerably higher than previously reported for its human orthologue, suggesting different environmental pressures in the two lineages. Known functional elements downstream of human NAT1 were also differentiated in rhesus macaque and other primates. Xenobiotic metabolizing enzymes play roles beyond mere protection from exogenous chemicals. Therefore, any link to disease, particularly carcinogenesis, may be via modulation of xenobiotic mutagenicity or more subtle interference with cell physiology. Comparative analyses add the evolutionary dimension to such investigations, assessing functional conservation/diversification among primates.

The circular RNome of primary breast cancer

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

Association of microRNA-7 and its binding partner CDR1-AS with the prognosis and prediction of 1st-line tamoxifen therapy in breast cancer

The large number of non-coding RNAs (ncRNAs) and their breadth of functionalities has fuelled many studies on their roles in cancer. We previously linked four microRNAs to breast cancer prognosis. One of these microRNAs, hsa-miR-7, was found to be regulated by another type of ncRNA, the circular non-coding RNA (circRNA) CDR1-AS, which contains multiple hsa-miR-7 binding sites. Based on this finding, we studied the potential clinical value of this circRNA on breast cancer prognosis in a cohort based on a cohort that was previously analysed for hsa-miR-7 and in an adjuvant hormone-naïve cohort for 1^st-line tamoxifen treatment outcomes, in which we also analysed hsa-miR-7. A negative correlation was observed between hsa-miR-7 and CDR1-AS in both cohorts. Despite associations with various clinical metrics (e.g., tumour grade, tumour size, and relapse location), CDR1-AS was neither prognostic nor predictive of relevant outcomes in our cohorts. However, we did observe stromal CDR1-AS expression, suggesting a possible cell-type specific interaction. Next to the known association of hsa-miR-7 expression with poor prognosis in primary breast cancer, we found that high hsa-miR-7 expression was predictive of an adverse response to tamoxifen therapy and poor progression-free and post-relapse overall survival in patients with recurrent disease.

Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons

BACKGROUND

Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data.

RESULTS

We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality.

CONCLUSIONS

We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

High mRNA expression of splice variant SYK short correlates with hepatic disease progression in chemonaive lymph node negative colon cancer patients

Objective

Overall and splice specific expression of Spleen Tyrosine Kinase (SYK) has been posed as a marker predicting both poor and favorable outcome in various epithelial malignancies. However, its role in colorectal cancer is largely unknown. The aim of this study was to explore the prognostic role of SYK in three cohorts of colon cancer patients.

Methods

Total messenger RNA (mRNA) expression of SYK, SYK(T), and mRNA expression of its two splice variants SYK short (S) and SYK long (L) were measured using quantitative reverse transcriptase (RT-qPCR) in 240 primary colon cancer patients (n = 160 patients with chemonaive lymph node negative [LNN] and n = 80 patients with adjuvant treated lymph node positive [LNP] colon cancer) and related to microsatellite instability (MSI), known colorectal cancer mutations, and disease-free (DFS), hepatic metastasis-free (HFS) and overall survival (OS). Two independent cohorts of patients with respectively 48 and 118 chemonaive LNN colon cancer were used for validation.

Results

Expression of SYK and its splice variants was significantly lower in tumors with MSI, and in KRAS wild type, BRAF mutant and PTEN mutant tumors. In a multivariate Cox regression analysis, as a continuous variable, increasing SYK(S) mRNA expression was associated with worse HFS (Hazard Ratio[HR] = 1.83; 95% Confidence Interval[CI] = 1.08–3.12; p = 0.026) in the LNN group, indicating a prognostic role for SYK(S) mRNA in patients with chemonaive LNN colon cancer. However, only a non-significant trend between SYK(S) and HFS in one of the two validation cohorts was observed (HR = 4.68; 95%CI = 0.75–29.15; p = 0.098).

Conclusion

In our cohort, we discovered SYK(S) as a significant prognostic marker for HFS for patients with untreated LNN colon cancer. This association could however not be confirmed in two independent smaller cohorts, suggesting that further extensive validation is needed to confirm the prognostic value of SYK(S) expression in chemonaive LNN colon cancer.

The DNA cytosine deaminase APOBEC3B promotes tamoxifen resistance in ER-positive breast cancer

Breast tumors often display extreme genetic heterogeneity characterized by hundreds of gross chromosomal aberrations and tens of thousands of somatic mutations. Tumor evolution is thought to be ongoing and driven by multiple mutagenic processes. A major outstanding question is whether primary tumors have preexisting mutations for therapy resistance or whether additional DNA damage and mutagenesis are necessary. Drug resistance is a key measure of tumor evolvability. If a resistance mutation preexists at the time of primary tumor presentation, then the intended therapy is likely to fail. However, if resistance does not preexist, then ongoing mutational processes still have the potential to undermine therapeutic efficacy. The antiviral enzyme APOBEC3B (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3B) preferentially deaminates DNA C-to-U, which results in signature C-to-T and C-to-G mutations commonly observed in breast tumors. We use clinical data and xenograft experiments to ask whether APOBEC3B contributes to ongoing breast tumor evolution and resistance to the selective estrogen receptor modulator, tamoxifen. First, APOBEC3B levels in primary estrogen receptor-positive (ER+) breast tumors inversely correlate with the clinical benefit of tamoxifen in the treatment of metastatic ER+ disease. Second, APOBEC3B depletion in an ER+ breast cancer cell line results in prolonged tamoxifen responses in murine xenograft experiments. Third, APOBEC3B overexpression accelerates the development of tamoxifen resistance in murine xenograft experiments by a mechanism that requires the enzyme's catalytic activity. These studies combine to indicate that APOBEC3B promotes drug resistance in breast cancer and that inhibiting APOBEC3B-dependent tumor evolvability may be an effective strategy to improve efficacies of targeted cancer therapies.

Clinical Relevance of Targeting the Gastrin-Releasing Peptide Receptor, Somatostatin Receptor 2, or Chemokine C-X-C Motif Receptor 4 in Breast Cancer for Imaging and Therapy

Imaging and therapy using radioligands targeting receptors overexpressed on tumor cells is successfully applied in neuroendocrine tumor patients. Because expression of the gastrin-releasing peptide receptor (GRPR), somatostatin receptor 2 (SSTR2), and chemokine C-X-C motif receptor 4 (CXCR4) has been demonstrated in breast cancer, targeting these receptors using radioligands might offer new imaging and therapeutic opportunities for breast cancer patients. The aim of this study was to correlate messenger RNA (mRNA) expression of GRPR, SSTR2, and CXCR4 with clinicopathologic and biologic factors, and with prognosis and prediction to therapy response, in order to identify specific breast cancer patient groups suited for the application of radioligands targeting these receptors.

METHODS

First, we studied GRPR and SSTR2 expression in 13 clinical breast cancer specimens by in vitro autoradiography and correlated this with corresponding mRNA levels to investigate whether mRNA levels reliably represent cell surface expression. Next, GRPR, SSTR2, and CXCR4 mRNA levels were measured by quantitative reverse transcriptase polymerase chain reaction in 915 primary breast cancer tissues and correlated with known clinicopathologic and biologic factors, disease-free survival, distant metastasis-free survival, and overall survival (DFS, MFS, and OS, respectively). In 224 adjuvant hormonal treatment-naïve estrogen receptor (ER, ESR1)-positive patients who received tamoxifen as first-line therapy for recurrent or metastatic disease, the expression levels of the receptors were correlated with progression-free survival.

RESULTS

Our results showed a significant positive correlation between GRPR and SSTR2 expression analyzed by in vitro autoradiography and by quantitative reverse transcriptase polymerase chain reaction (Spearman's rank correlation coefficient [Rs]=0.94, P<0.001, and Rs=0.73, P=0.0042, respectively). Furthermore, high GRPR and SSTR2 mRNA levels were observed more frequently in ESR1-positive specimens, whereas high CXCR4 expression was associated with ESR1-negative specimens. Also, high mRNA expression of CXCR4 was associated with a prolonged DFS, MFS, and OS (multivariate hazard ratio MFS=0.76 [95% confidence interval, 0.64-0.90], P=0.001), whereas high mRNA levels of GRPR were associated with a prolonged progression-free survival after the start of first-line tamoxifen treatment (multivariate hazard ratio=0.68 [95% confidence interval, 0.48-0.97], P=0.031).

CONCLUSION

Our data indicate that imaging and therapy using GRPR or SSTR2 radioligands might especially be beneficial for ESR1-positive breast cancer and CXCR4 radioligands for ESR1-negative breast cancer.

Tamoxifen induces a pluripotency signature in breast cancer cells and human tumors

Tamoxifen is the treatment of choice in estrogen receptor alpha breast cancer patients that are eligible for adjuvant endocrine therapy. However, ∼50% of ERα-positive tumors exhibit intrinsic or rapidly acquire resistance to endocrine treatment. Unfortunately, prediction of de novo resistance to endocrine therapy and/or assessment of relapse likelihood remain difficult. While several mechanisms regulating the acquisition and the maintenance of endocrine resistance have been reported, there are several aspects of this phenomenon that need to be further elucidated. Altered metabolic fate of tamoxifen within patients and emergence of tamoxifen-resistant clones, driven by evolution of the disease phenotype during treatment, appear as the most compelling hypotheses so far. In addition, tamoxifen was reported to induce pluripotency in breast cancer cell lines, in vitro. In this context, we have performed a whole transcriptome analysis of an ERα-positive (T47D) and a triple-negative breast cancer cell line (MDA-MB-231), exposed to tamoxifen for a short time frame (hours), in order to identify how early pluripotency-related effects of tamoxifen may occur. Our ultimate goal was to identify whether the transcriptional actions of tamoxifen related to induction of pluripotency are mediated through specific ER-dependent or independent mechanisms. We report that even as early as 3 hours after the exposure of breast cancer cells to tamoxifen, a subset of ERα-dependent genes associated with developmental processes and pluripotency are induced and this is accompanied by specific phenotypic changes (expression of pluripotency-related proteins). Furthermore we report an association between the increased expression of pluripotency-related genes in ERα-positive breast cancer tissues samples and disease relapse after tamoxifen therapy. Finally we describe that in a small group of ERα-positive breast cancer patients, with disease relapse after surgery and tamoxifen treatment, ALDH1A1 (a marker of pluripotency in epithelial cancers which is absent in normal breast tissue) is increased in relapsing tumors, with a concurrent modification of its intra-cellular localization. Our data could be of value in the discrimination of patients susceptible to develop tamoxifen resistance and in the selection of optimized patient-tailored therapies.

Evaluation of the ability of adjuvant tamoxifen-benefit gene signatures to predict outcome of hormone-naive estrogen receptor-positive breast cancer patients treated with tamoxifen in the advanced setting

To identify molecular markers indicative of response to tamoxifen and easily implemented in the routine setting, we recently reported three gene signatures that could stratify post-menopausal tamoxifen-treated, estrogen receptor-positive (ER+) patients according to outcome in the adjuvant setting. Here, we evaluated the predictive potential of the total of 14 genes included in the 3 gene signatures using 2 hormone-naïve Dutch ER+ cohorts of a total of 285 recurrent breast cancer patients treated with first-line tamoxifen. mRNA levels were measured by reverse transcriptase quantitative PCR (RT-qPCR) and the length of progression-free survival (PFS) was used as the primary endpoint. A Mann-Whitney U test was used to select for differentially expressed genes between tumors of patients who showed or did not show progressive disease within 6 months after start of tamoxifen treatment. Cox univariate and multivariate regression analysis for PFS were used to further assess their (independent) predictive potential. Five (BCAR3, BCL2, ESR1, IGF1R, and NCOA1) of the 14 genes analyzed showed significantly higher mRNA levels in tumors of patients who showed no disease progression within 6 months. Only BCAR3, BCL2 and NAT1 were significantly associated with a favorable PFS in multivariate analysis that included the traditional predictive factors: age, dominant relapse site, disease-free interval, ER and progesterone receptor (PGR), and adjuvant chemotherapy. This study shows that BCAR3, BCL2 and NAT1 in particular exhibit predictive promise regarding the efficacy of tamoxifen treatment in recurrent disease, in addition to the previously shown favorable outcome in the adjuvant setting.